matlab 问题

nkkopd

Error using lmiterm (line 281)
lhs of LMI #1, block (1,1): term dimensions incompatible with other
4 J  I% V# p+ v1 N: C. L: lterms in row #1 or column #1

Error in SCP_03 (line 103)
lmiterm([1 1 1 P1],-1,1)
; T- }: }( u) @) n0 ~8 p' N
clear
clc
" d" c& I* }( a; A1 B
% t; b( h8 ?( S  Ggamma = 0.95;

% M0 Q% F7 p( U& C  qsigma1 = 2;
9 B8 A# U7 z/ V. o$ u- |# ^  {5 Gsigma2 = 1.98;
sigma3 = 1.86;
' Q; l+ M- q* b) c$ `* Msigma = sigma1+sigma2+sigma3;
$ [# W$ a6 S( d- f) P" E* q
gamma1 = 0.89;
- U' [/ ^" X7 K. g- @' \% H- egamma2 = 0.88;
( `' D% n( E! [gamma3 = 0.85;
Gamma = [gamma1, gamma1, gamma1; gamma2, gamma2, gamma2; gamma3, gamma3, gamma3];%3*3

Lambda = diag([2, 2]);%2*2
) v1 `6 k7 k# P6 \" E' h
* a8 }0 v$ s) s' Z" @" ~8 w, uW = [kron(Gamma, Lambda), zeros(6,3); zeros(3,6), zeros(3)];%9*9

3 t# ]/ L& }3 ?  u" n2 A5 I& U; |beta1 = 0.88;
beta2 = 0.85;
beta3 = 0.87;
Xi = diag([beta1, beta2, beta3]);%3*3
barXi = [kron(Xi, eye(2)), zeros(6,3); zeros(3,6), zeros(3)];%9*9
0 h* j7 h3 o/ J5 ]3 M
U = diag([1.68, 1.5, 1.68, 1.5, 1.68, 1.5, 0, 0, 0]);%花U1+U2   9*9
/ A, v* t, r2 }7 PL = diag([0.34, 0.34, 0.34, 0.34, 0.34, 0.34, 0, 0, 0]);%花L1+L2  9*9
: R& i$ K7 l$ b: b$ ?3 |% }* S
N = [0.1, 0.1, 0, 0, 0, 0, 0, 0, 0; 0, 0, 0.1, 0.1, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0.1, 0.1, 0, 0, 0];%\barN   3*9
M = [0.1, 0, 0; 0.2, 0, 0; 0, 0.1, 0; 0, 0.1, 0; 0, 0, 0.2; 0, 0, 0.1;zeros(3)];%\bARM   3*9

C =[0.5, -0.4, 0, 0, 0, 0; 0, 0, 0.3, 0.2, 0, 0; 0, 0, 0, 0, 0.4, -0.5];%3*6
, X; A# M; T9 T( w! a: F
K1 = [0.1, -0.2, 0.1; 0.3, -0.2, 0.1; 0.25, 0.1, 0.1; -0.2, 0.1, -0.3; -0.1, 0.4 0.25; 0.1, 0.2, 0.3; 0.3, -0.2, 0.2; 0.1, 0.25, 0.3; 0.15, 0.4, 0.25];%9*3
8 Q# T! y) b. ]6 [; kK2 = [-0.1, 0.25, 0.1; 0.23, -0.2, 0.15; 0.25, -0.1, 0.1; 0.2, 0.15, 0.3; -0.1, 0.3 0.25; 0.1, 0.2, 0.3; 0.3, -0.2, 0.2; 0.1, 0.25, -0.3; 0.1, 0.4, 0.25];%9*3
9 Z# o: S$ g% \7 qK3 =[0.1, -0.2, 0.15; 0.3, -0.2, 0.1; 0.25, 0.1, 0.1; -0.2, 0.1, 0.3; -0.1, 0.4 0.2; -0.15, 0.2, 0.3; 0.3, -0.2, 0.25; -0.1, 0.25, 0.3; 0.15, 0.35, 0.25];%9*3

2 g+ k) x6 \$ O( K5 BH = [0.2, 0.8, 0, 0, 0, 0, 0, 0, 0; 0, 0, 0.2, 0.8, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0.2, 0.8, 0, 0, 0];%\barH   3*9

E = [0.3;0.3;0.3];%3*1
7 E& I) v8 S$ c4 f3 b
D = diag([0.5, 0.7, 0.4]);%3*3

A = [0.2, 0, 0, 0, 0, 0; 0.1, -0.3, 0, 0, 0, 0; 0, 0, 0.2, 0.1, 0, 0; 0, 0, 0.1, -0.2, 0, 0; 0, 0, 0, 0, -0.3, 0.1; 0, 0, 0, 0, 0, -0.2];%6*6

0 k/ i7 ], P5 `" P% v- S1 M  F& ]B = [0.6, 0, 0; 0.3, 0, 0; 0, 0.4, 0; 0, 0.5, 0; 0, 0, 0.8; 0, 0, 0.5];%6*3

theta(1) = 2;
! F& S/ j- m; s, z9 q4 q) Q5 htheta(2) = 3;
& _5 f0 g8 G; F0 |0 r! x6 htheta(3) = 1;
theta(4) = 3;

5 O2 h8 A+ V5 V* J
+ I+ _1 D8 y0 M8 Q- z7 B
& T) |9 _! c, M7 n; W+ cfor k = 1:3
) H& r- A8 ?7 [* B) Z4 N) _9 I    i = theta(k);
$ K* q: E0 D, }& T" L: K6 c    j = theta(k+1);
9 o& ?! j7 U7 b( V   
3 S$ G( y: z; x5 \  Phi = zeros(3);
  
& J2 x6 n& i& u( T9 j; `for t = 1:3
     
8 D! ?, A8 m0 T    if (i-t) ~= 0;
( W- L( P4 @- H8 p1 Z3 Q: ^$ _8 `          delta(t) = 0;
/ V0 G- i% K, K" }) f! Y; ?    elseif (i-t) == 0;
, Q" H% G7 |3 L7 z8 M/ a6 ~9 U          delta(t) = 1;
         Phi(t,t) = delta(t);
    end   
$ U4 A* F/ c; ~( `  
* P3 J- q- N2 @$ k1 wend
( u3 \$ f4 j# b" z3 h
  barA = [A, zeros(6,3); (kron(Phi,eye(1)))*C, kron((eye(3)-Phi),eye(1))];%9*9
$ P  l# u8 I$ d1 b  barB = [zeros(6,1), B; (kron(Phi,eye(1)))*E, (kron(Phi,eye(1)))*D];%9*4
. t. S$ X" |- P: K4 t  barC = [(kron(Phi,eye(1)))*C, kron((eye(3)-Phi),eye(1))];%3*9
  barD = [(kron(Phi,eye(1)))*E, kron(Phi,eye(1))*D];%3*4
. J7 [7 i' F+ j, zend
  % N=3  n=2  m=1 p=1 r=1 q=1  M=3
  setlmis([])
2 N( U1 b* ~, o) O  P11 = lmivar(1,[2 1]);
3 r1 K2 j+ s1 M& G2 _  P21 = lmivar(2,[2 1]);
  P31 = lmivar(1,[1 1]);
  P1 = lmivar(3,[kron(eye(3),P11),kron(eye(3),P21);kron(eye(3),P21'),kron(eye(3),P31)]);%9*9
( H. h2 `# ?" D0 @. r4 e  P12 = lmivar(1,[2 1]);
: y6 d" \/ m. c  P22 = lmivar(2,[2 1]);
5 c  S& t3 a$ l% r  P32 = lmivar(1,[1 1]);
  P2 = lmivar(3,[kron(eye(3),P12),kron(eye(3),P22);kron(eye(3),P22'),kron(eye(3),P32)]);%9*9
  P13 = lmivar(1,[2 1]);
  P23 = lmivar(2,[2 1]);
  P33 = lmivar(1,[1 1]);
& O/ }% ~% e. x( ^3 ~2 t2 ^  P3 = lmivar(3,[kron(eye(3),P13),kron(eye(3),P23);kron(eye(3),P23'),kron(eye(3),P33)]);%9*9
. b. v" o7 }8 ?9 i( p% o" S* {  Z1 = lmivar(2,[9 3]);
4 B# g$ U. `9 a' _, t/ s( U- G  Z2 = lmivar(2,[9 3]);
+ @; q6 B. d. {* l1 l  Z3 = lmivar(2,[9 3]);
  Q = lmivar(1,[9 1]);
  mu1 = lmivar(1,[1,1]);
; ^: z" d! E$ h1 K5 u7 R4 U1 Z  mu2 = lmivar(1,[1,1]);
  mu3 = lmivar(1,[1,1]);
  s = lmivar(1,[1,1]);
  z' B6 V! M+ P/ K  barGamma = lmivar(1,[9 1]);
  barLambda = lmivar(1,[9 1]);
  
lmiterm([1 1 1 barLambda],2,1)
1 H/ t- `5 B- F; Z' zlmiterm([1 1 1 P1],-1,1)
lmiterm([1 1 1 mu1],-1,U+U')
lmiterm([1 1 1 mu2],-1,L+L')
lmiterm([1 1 1 s],1,N'*N)
2 }9 x/ X$ y4 g+ zlmiterm([1 1 1 0],H'*H)
lmiterm([1 1 5 mu1],1,U')
. ^/ y, G3 O/ W, clmiterm([1 1 7 mu2],1,L')
1 e1 X# ?5 h$ N, w  Slmiterm([1 1 12 P2],barA'+W',1)
/ ?2 {! X* [5 Y: Zlmiterm([1 2 2 P2],barC'*K1',K1*barC)
lmiterm([1 2 2 P2],barC'*K2',K2*barC)
lmiterm([1 2 2 P2],barC'*K3',K3*barC)
lmiterm([1 2 2 P1],-1,1)
lmiterm([1 2 2 mu1],-1,U+U')
$ U5 f6 R: l: }lmiterm([1 2 2 mu2],-1,L+L')
8 R: S7 ~$ d" Y" ?lmiterm([1 2 6 mu1],1,U')
lmiterm([1 2 8 mu2],1,L')
lmiterm([1 2 11 P2],barC'*K1',K1*barD)
lmiterm([1 2 11 P2],barC'*K2',K2*barD)
, u* |7 D1 ^7 q9 ~8 ulmiterm([1 2 11 P2],barC'*K3',K3*barD)
lmiterm([1 2 13 P2],barA'+W',1)
lmiterm([1 2 13 -Z2],-barC',1)
' j0 b& a1 Y" e% O. {3 |lmiterm([1 3 3 mu3],-1,L')
lmiterm([1 3 9 mu3],1,L+L')
8 D' w6 g, }3 f! g8 T/ dlmiterm([1 4 4 mu3],-1,L+L')
( C/ s+ l' F* |5 K/ I: ulmiterm([1 4 10 mu3],1,L')
$ o' r& x1 G. v# klmiterm([1 5 5 mu1],-2,eye(9))
! Z5 o2 e+ i7 F; L/ Ilmiterm([1 5 12 P2],1,1)
$ w7 b1 H( b* ?0 C- Flmiterm([1 6 6 mu1],-2,eye(9))
lmiterm([1 6 13 P2],1,1)
9 R* \+ A- M$ Z  M. t) glmiterm([1 7 7 Q],1,1)
lmiterm([1 7 7 mu2],-2,eye(9))
% L; S  p) d0 T4 Qlmiterm([1 8 8 Q],1,1)
! b1 K5 q$ x, m) Mlmiterm([1 8 8 mu2],-2,eye(9))
$ A! b+ _0 k; x+ D% k) ~lmiterm([1 9 9 barGamma],2,1)
lmiterm([1 9 9 Q],-1,1)
( a$ t  t% |  o* ~lmiterm([1 9 9 mu3],-2,eye(9))
lmiterm([1 9 12 P2],barXi',1)
% y5 `" s2 J0 @' ?lmiterm([1 10 10 Q],-1,1)
2 u1 d( X6 `! \8 ^lmiterm([1 10 10 mu3],-2,eye(9))
lmiterm([1 10 13 P2],barXi',1)
' f4 m( @+ o7 T2 N) Ylmiterm([1 11 11 P2],barD'*K1',K1*barD)
lmiterm([1 11 11 P2],barD'*K2',K2*barD)
lmiterm([1 11 11 P2],barD'*K3',K3*barD)
- k+ n* e( x! D5 m/ V5 B" Jlmiterm([1 11 11 0],-(gamma^2)*eye(4))
lmiterm([1 11 12 P2],barB',1)
lmiterm([1 11 13 P2],barB',1)
lmiterm([1 11 13 -Z2],-barD',1)
4 l( K- t) G7 u. D- y" qlmiterm([1 12 14 P2],1,M)
lmiterm([1 13 14 P2],1,M)
' |3 w% k, v  y) d9 O! `, T
lmiterm([-1 1 1 P1],1,1)
: d; p: Q( N2 e+ ~lmiterm([-2 1 1 P2],1,1)
& B# n6 Y& C3 E3 f  {; O& N2 plmiterm([-3 1 1 P3],1,1)
9 ^; @  a) D" ?5 F7 q3 a2 r
# y3 q+ |1 `) H- b/ N% V* zlmiterm([2 1 1 barLambda],2,1)
+ Y$ G& x- y2 R  |( Hlmiterm([2 1 1 P2],-1,1)
lmiterm([2 1 1 mu1],-1,U+U')
1 k  r  A2 u! d; Q9 T# d1 ulmiterm([2 1 1 mu2],-1,L+L')
2 Q; d3 |1 R1 {; Y" `" H& J( ilmiterm([2 1 1 s],1,N'*N)
) s" i8 Y2 M  S# e: flmiterm([2 1 1 0],H'*H)
, X. b# l6 q' |# T1 V. Clmiterm([2 1 5 mu1],1,U')
+ w7 N' R# ^8 L- X1 [lmiterm([2 1 7 mu2],1,L')
' t# p3 g# @" d) H( E, \6 @# Almiterm([2 1 12 P3],barA'+W',1)
lmiterm([2 2 2 P3],barC'*K1',K1*barC)
8 u" z* \: V, f! Jlmiterm([2 2 2 P3],barC'*K2',K2*barC)
lmiterm([2 2 2 P3],barC'*K3',K3*barC)
lmiterm([2 2 2 P2],-1,1)
lmiterm([2 2 2 mu1],-1,U+U')
( P* k* A! e1 Wlmiterm([2 2 2 mu2],-1,L+L')
lmiterm([2 2 6 mu1],1,U')
lmiterm([2 2 8 mu2],1,L')
8 S/ t& r/ N2 j% w8 Y; r5 Olmiterm([2 2 11 P3],barC'*K1',K1*barD)
lmiterm([2 2 11 P3],barC'*K2',K2*barD)
lmiterm([2 2 11 P3],barC'*K3',K3*barD)
lmiterm([2 2 13 P3],barA'+W',1)
  `  b3 Y. s' P& Jlmiterm([2 2 13 -Z3],-barC',1)
lmiterm([2 3 3 mu3],-1,L')
7 o* X9 _$ T% |. g; Xlmiterm([2 3 9 mu3],1,L+L')
lmiterm([2 4 4 mu3],-1,L+L')
) T/ Z8 T+ R6 R; Jlmiterm([2 4 10 mu3],1,L')
lmiterm([2 5 5 mu1],-2,eye(9))
: ~! x1 ~5 r5 m/ p- Y# N- jlmiterm([2 5 12 P3],1,1)
lmiterm([2 6 6 mu1],-2,eye(9))
lmiterm([2 6 13 P3],1,1)
1 Z: F( B! ]8 L( B" _8 ?lmiterm([2 7 7 Q],1,1)
lmiterm([2 7 7 mu2],-2,eye(9))
lmiterm([2 8 8 Q],1,1)
2 Y1 z* h. _  a6 j. `  Jlmiterm([2 8 8 mu2],-2,eye(9))
lmiterm([2 9 9 barGamma],2,1)
& P- m8 Q4 u3 v2 U- [1 ?: Ilmiterm([2 9 9 Q],-1,1)
lmiterm([2 9 9 mu3],-2,eye(9))
lmiterm([2 9 12 P3],barXi',1)
! w  \9 l# m/ I1 M6 Llmiterm([2 10 10 Q],-1,1)
lmiterm([2 10 10 mu3],-2,eye(9))
" m, G! r. }# |# \5 O4 llmiterm([2 10 13 P3],barXi',1)
lmiterm([2 11 11 P3],barD'*K1',K1*barD)
/ U( q7 n+ \/ H1 O: S' f- ]* G; s8 blmiterm([2 11 11 P3],barD'*K2',K2*barD)
lmiterm([2 11 11 P3],barD'*K3',K3*barD)
lmiterm([2 11 11 0],-(gamma^2)*eye(4))
' j7 \1 o: V+ J/ i2 qlmiterm([2 11 12 P3],barB',1)
lmiterm([2 11 13 P3],barB',1)
) }( R$ Y) y& ]/ D; ilmiterm([2 11 13 -Z3],-barD',1)
. a$ M' S* Q% E$ hlmiterm([2 12 14 P3],1,M)
& e5 X0 |3 f$ X7 I0 mlmiterm([2 13 14 P3],1,M)

G=getlmis;
- x1 J* B5 b2 V. Q( {- b[tmin b]=feasp(G);
* W/ R# O. H+ S  v* F
P11 = dec2mat(G,b,P11);
P21 = dec2mat(G,b,P21);
P31 = dec2mat(G,b,P31);
% n& j& {: B5 S6 A2 DP1 = [kron(eye(3), P11), kron(eye(3), P21); kron(eye(3), P21'), kron(eye(3), P31)]
P12 = dec2mat(G,b,P12);
; J1 d( @4 E- L: tP22 = dec2mat(G,b,P22);
9 f7 O) a$ P7 ?5 YP32 = dec2mat(G,b,P32);
. Q- K( T' o% O7 f( r/ gP2 = [kron(eye(3), P12), kron(eye(3), P22); kron(eye(3), P22'), kron(eye(3), P32)]
P13 = dec2mat(G,b,P13);
+ B8 q2 n. _9 s, xP23 = dec2mat(G,b,P23);
P33 = dec2mat(G,b,P33);
: t" r2 h; e( [. I' x3 bP3 = [kron(eye(3), P13), kron(eye(3), P23); kron(eye(3), P23'), kron(eye(3), P33)]
Z(1) = dec2mat(G,b,Z(1));
Z(2) = dec2mat(G,b,Z(2));
Z(3) = dec2mat(G,b,Z(3));
barGamma = [kron(Xi*(eye(3)-Xi), P(1i).P(i)), zeros(6,3); zeros(3,6), zeros(3)];
barLambda = [kron(sigma*eye(3), Lambda'*(P(1i).P(i))*Lambda), zeros(6,3); zeros(3,6), zeros(3)];
( w0 j5 E9 B7 o. J$ PQ = dec2mat(G,b,Q)
1 M: ~3 r% G. ^: p# omu1 = dec2mat(G,b,mu1)
9 h8 ?- K& p' j$ n+ p5 |mu2 = dec2mat(G,b,mu2)
+ h: @: l: r" ]mu3 = dec2mat(G,b,mu3)
s = dec2mat(G,b,s)
9 m. j% }) X/ `) o8 k' R

zaiyiaaaa

zzz.dan

找到对应行，调试一下就知道了

IBB-EUT

定义变量时的P1维数不一致，P21和P31不一致