matlab 问题

nkkopd

Error using lmiterm (line 281)
lhs of LMI #1, block (1,1): term dimensions incompatible with other
2 l/ I4 G+ |; B7 W" pterms in row #1 or column #1

, ?7 [' I4 U2 |1 F1 l8 A4 v) ^6 }Error in SCP_03 (line 103)
lmiterm([1 1 1 P1],-1,1)

7 B- `# [7 `6 U6 Lclear
: w( a+ V/ I) Y6 xclc

gamma = 0.95;

# c1 l& [5 f$ t5 T- B) ?7 X. Ysigma1 = 2;
sigma2 = 1.98;
sigma3 = 1.86;
+ l& D0 u$ ]5 k, n3 Y: ]) C( }5 J6 jsigma = sigma1+sigma2+sigma3;

. K2 N! o% p1 e3 Zgamma1 = 0.89;
3 `" `" f% v( S" i6 ?8 Vgamma2 = 0.88;
0 U+ e& W1 q  t2 V- sgamma3 = 0.85;
* K7 ?: a: x6 Z/ Q3 iGamma = [gamma1, gamma1, gamma1; gamma2, gamma2, gamma2; gamma3, gamma3, gamma3];%3*3

Lambda = diag([2, 2]);%2*2

9 _9 M( L! V/ u: UW = [kron(Gamma, Lambda), zeros(6,3); zeros(3,6), zeros(3)];%9*9
/ k& H9 a2 H+ [6 ]5 b
6 J3 L; _" P; F$ T0 F3 R8 Zbeta1 = 0.88;
0 {3 d# ?( g# D- cbeta2 = 0.85;
beta3 = 0.87;
Xi = diag([beta1, beta2, beta3]);%3*3
: A* S/ I2 A( r4 Z/ ubarXi = [kron(Xi, eye(2)), zeros(6,3); zeros(3,6), zeros(3)];%9*9

0 z5 s; @; o6 B- }* u  B, @4 V: ]U = diag([1.68, 1.5, 1.68, 1.5, 1.68, 1.5, 0, 0, 0]);%花U1+U2   9*9
; Z5 T! H, q3 S9 v& A3 YL = diag([0.34, 0.34, 0.34, 0.34, 0.34, 0.34, 0, 0, 0]);%花L1+L2  9*9
* P/ F( U* P7 v) g0 u( A/ z8 W4 d
* ~1 Y+ b5 L) z1 q/ nN = [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

8 {, n! {+ n' D5 W% r) _9 \8 RK1 = [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
4 S+ t1 _  p$ `- }; t1 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
K3 =[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

+ Q: V1 f* n5 U8 h' y% q6 [( v' iH = [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

6 N% E- m  K8 z/ N, H% ?9 aE = [0.3;0.3;0.3];%3*1

  Y/ E0 w0 f7 l7 G8 \$ {0 AD = diag([0.5, 0.7, 0.4]);%3*3
4 R; q6 j2 u: P/ V8 e' d9 T( f- w
" v& F% h$ p& Z' k# SA = [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

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;
' q8 w9 _' O3 l7 s. ctheta(2) = 3;
theta(3) = 1;
1 a, P% \6 |! g2 \: r. ]theta(4) = 3;
3 Q# U7 e9 P3 z5 D. q

8 ^% U; Z( [9 w& v7 A1 T
for k = 1:3
    i = theta(k);
    j = theta(k+1);
   
6 {7 S4 ^# T2 C; E  Phi = zeros(3);
  
for t = 1:3
     
& B8 X3 B+ @% e5 B, L8 h5 ?5 q3 X    if (i-t) ~= 0;
          delta(t) = 0;
4 y+ [( k+ s2 R" s    elseif (i-t) == 0;
          delta(t) = 1;
         Phi(t,t) = delta(t);
0 j4 X. p2 g2 v* p    end   
( V( {( K( i1 q9 e' l9 U$ _  
end
4 y7 ]! k2 t  Y! x- K6 Y7 u: j, z) p
  barA = [A, zeros(6,3); (kron(Phi,eye(1)))*C, kron((eye(3)-Phi),eye(1))];%9*9
# g4 r( V' o% h1 X  barB = [zeros(6,1), B; (kron(Phi,eye(1)))*E, (kron(Phi,eye(1)))*D];%9*4
  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
) h: x* G' f: x/ r' W9 Vend
  % N=3  n=2  m=1 p=1 r=1 q=1  M=3
  setlmis([])
% n5 B4 z" \; D. ~9 ~% x: u  P11 = lmivar(1,[2 1]);
! F; V2 u& X& s; r; Q  G* B  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
  P12 = lmivar(1,[2 1]);
  P22 = lmivar(2,[2 1]);
- {# I4 W  E3 r' h( M  c  P32 = lmivar(1,[1 1]);
8 M1 q1 E1 x/ h  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]);
6 g: g+ c7 [8 k8 L$ q  P23 = lmivar(2,[2 1]);
! ~) u7 S" M, I" i# L* h' q  P33 = lmivar(1,[1 1]);
  P3 = lmivar(3,[kron(eye(3),P13),kron(eye(3),P23);kron(eye(3),P23'),kron(eye(3),P33)]);%9*9
  Z1 = lmivar(2,[9 3]);
  Z2 = lmivar(2,[9 3]);
  Z3 = lmivar(2,[9 3]);
5 o4 {7 U+ |% }$ N+ |  Q = lmivar(1,[9 1]);
  mu1 = lmivar(1,[1,1]);
  mu2 = lmivar(1,[1,1]);
% H3 a& h# m& w2 [  mu3 = lmivar(1,[1,1]);
  s = lmivar(1,[1,1]);
6 u- J( ~2 R2 m! z  barGamma = lmivar(1,[9 1]);
  barLambda = lmivar(1,[9 1]);
7 e8 G' e. b* Y+ ~6 b+ f  A  
lmiterm([1 1 1 barLambda],2,1)
/ k; P9 `  h6 q, _lmiterm([1 1 1 P1],-1,1)
lmiterm([1 1 1 mu1],-1,U+U')
lmiterm([1 1 1 mu2],-1,L+L')
4 L; E4 ~# u3 X, L8 s) ], r) Clmiterm([1 1 1 s],1,N'*N)
lmiterm([1 1 1 0],H'*H)
9 J1 @1 L% I  U' olmiterm([1 1 5 mu1],1,U')
! S% y7 C( Y) j* W2 f0 T/ Clmiterm([1 1 7 mu2],1,L')
4 \3 h9 `. u4 J0 ^$ s! ]lmiterm([1 1 12 P2],barA'+W',1)
7 h; l/ ^, ?% n+ }8 [) k  a. Qlmiterm([1 2 2 P2],barC'*K1',K1*barC)
5 q8 M. r2 I) |" _! y4 M5 t3 d' xlmiterm([1 2 2 P2],barC'*K2',K2*barC)
0 ~: E+ W9 ?7 v6 `9 @8 X# Qlmiterm([1 2 2 P2],barC'*K3',K3*barC)
8 G1 b4 Z% o3 U; n# }lmiterm([1 2 2 P1],-1,1)
+ i: ]3 i0 j& Q$ P) c2 `) K0 c$ k6 \4 Slmiterm([1 2 2 mu1],-1,U+U')
# M. h3 h$ M  e2 k: V" blmiterm([1 2 2 mu2],-1,L+L')
9 u$ n8 U, Q5 s% ulmiterm([1 2 6 mu1],1,U')
0 v) \! B: `% `, S+ }1 ~( vlmiterm([1 2 8 mu2],1,L')
# M* m5 y: s9 W# D- z5 ulmiterm([1 2 11 P2],barC'*K1',K1*barD)
) Q# }% N0 ]8 s1 p4 ~7 J5 Z  D( _lmiterm([1 2 11 P2],barC'*K2',K2*barD)
: E% @$ g; J& J' Z$ rlmiterm([1 2 11 P2],barC'*K3',K3*barD)
# X8 p# I6 s: K0 @1 R0 l4 w3 Klmiterm([1 2 13 P2],barA'+W',1)
0 T- p# C" ^: R: N% e3 K# e% k3 u& j  Ulmiterm([1 2 13 -Z2],-barC',1)
lmiterm([1 3 3 mu3],-1,L')
2 I* Z2 D. R; Clmiterm([1 3 9 mu3],1,L+L')
lmiterm([1 4 4 mu3],-1,L+L')
( E, b" S8 I7 n7 w( klmiterm([1 4 10 mu3],1,L')
lmiterm([1 5 5 mu1],-2,eye(9))
lmiterm([1 5 12 P2],1,1)
lmiterm([1 6 6 mu1],-2,eye(9))
0 ~6 @$ r! E: I+ {; u/ {lmiterm([1 6 13 P2],1,1)
  R( W* ~4 {3 Ilmiterm([1 7 7 Q],1,1)
lmiterm([1 7 7 mu2],-2,eye(9))
) _  X" z, n$ y7 ylmiterm([1 8 8 Q],1,1)
lmiterm([1 8 8 mu2],-2,eye(9))
lmiterm([1 9 9 barGamma],2,1)
$ x$ N& ^8 I9 V5 b& B- Tlmiterm([1 9 9 Q],-1,1)
* _5 P. I. X) @& r2 d, ilmiterm([1 9 9 mu3],-2,eye(9))
lmiterm([1 9 12 P2],barXi',1)
lmiterm([1 10 10 Q],-1,1)
4 U) s" m) U, B+ Elmiterm([1 10 10 mu3],-2,eye(9))
" ~: \" z- {% D: m7 g1 P. plmiterm([1 10 13 P2],barXi',1)
/ X) }+ n) h, n# ]0 f5 vlmiterm([1 11 11 P2],barD'*K1',K1*barD)
lmiterm([1 11 11 P2],barD'*K2',K2*barD)
- ~& ^4 i! E2 S  M3 H! mlmiterm([1 11 11 P2],barD'*K3',K3*barD)
lmiterm([1 11 11 0],-(gamma^2)*eye(4))
lmiterm([1 11 12 P2],barB',1)
$ h! A7 X+ z) D) o- j$ r7 M& b# Plmiterm([1 11 13 P2],barB',1)
lmiterm([1 11 13 -Z2],-barD',1)
lmiterm([1 12 14 P2],1,M)
lmiterm([1 13 14 P2],1,M)
. b- J& c- v* B. p: A
lmiterm([-1 1 1 P1],1,1)
lmiterm([-2 1 1 P2],1,1)
lmiterm([-3 1 1 P3],1,1)
# m# g, `4 t! ?
0 c& B1 C4 j7 G# W6 l( {+ V5 G/ rlmiterm([2 1 1 barLambda],2,1)
lmiterm([2 1 1 P2],-1,1)
lmiterm([2 1 1 mu1],-1,U+U')
lmiterm([2 1 1 mu2],-1,L+L')
  _0 g, b4 ^2 C% E% [lmiterm([2 1 1 s],1,N'*N)
lmiterm([2 1 1 0],H'*H)
7 S. g4 |- c4 ?0 rlmiterm([2 1 5 mu1],1,U')
lmiterm([2 1 7 mu2],1,L')
5 s+ j* t- l+ Z/ |  }# mlmiterm([2 1 12 P3],barA'+W',1)
/ z/ h2 A) V8 |. Clmiterm([2 2 2 P3],barC'*K1',K1*barC)
/ N. j8 x" B/ [" Qlmiterm([2 2 2 P3],barC'*K2',K2*barC)
lmiterm([2 2 2 P3],barC'*K3',K3*barC)
3 m+ N6 T4 s6 W% Z3 Glmiterm([2 2 2 P2],-1,1)
" Q! G9 j6 s; L  Tlmiterm([2 2 2 mu1],-1,U+U')
lmiterm([2 2 2 mu2],-1,L+L')
lmiterm([2 2 6 mu1],1,U')
lmiterm([2 2 8 mu2],1,L')
lmiterm([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)
lmiterm([2 2 13 -Z3],-barC',1)
lmiterm([2 3 3 mu3],-1,L')
9 N  Y: o/ J; f2 ~9 Nlmiterm([2 3 9 mu3],1,L+L')
lmiterm([2 4 4 mu3],-1,L+L')
% `6 n6 L4 |1 Y! K  ^' Wlmiterm([2 4 10 mu3],1,L')
lmiterm([2 5 5 mu1],-2,eye(9))
+ v1 ^& F& `) y, Nlmiterm([2 5 12 P3],1,1)
lmiterm([2 6 6 mu1],-2,eye(9))
) t5 j1 `4 y# @0 H; [1 m2 C$ G% @+ ]lmiterm([2 6 13 P3],1,1)
lmiterm([2 7 7 Q],1,1)
lmiterm([2 7 7 mu2],-2,eye(9))
lmiterm([2 8 8 Q],1,1)
lmiterm([2 8 8 mu2],-2,eye(9))
/ @7 {6 G5 F0 ^' a! j8 H9 N: olmiterm([2 9 9 barGamma],2,1)
+ B/ y/ y: N. A) i' I4 b) Almiterm([2 9 9 Q],-1,1)
1 T1 }7 \* l( U! j9 @  ~lmiterm([2 9 9 mu3],-2,eye(9))
/ o: w: r. a) q. B/ Flmiterm([2 9 12 P3],barXi',1)
9 d5 ~& l2 J' d( |- n& E% Slmiterm([2 10 10 Q],-1,1)
3 d7 o0 d  h+ v4 ilmiterm([2 10 10 mu3],-2,eye(9))
lmiterm([2 10 13 P3],barXi',1)
lmiterm([2 11 11 P3],barD'*K1',K1*barD)
2 w; b0 }# C1 x" R( H5 plmiterm([2 11 11 P3],barD'*K2',K2*barD)
lmiterm([2 11 11 P3],barD'*K3',K3*barD)
" M/ [* s; t5 w6 l$ s1 {lmiterm([2 11 11 0],-(gamma^2)*eye(4))
7 r2 n! }2 H4 @4 O5 R7 I5 T# j/ qlmiterm([2 11 12 P3],barB',1)
lmiterm([2 11 13 P3],barB',1)
8 d+ Q% x6 D) L" v3 }lmiterm([2 11 13 -Z3],-barD',1)
lmiterm([2 12 14 P3],1,M)
2 B8 W; G5 o; ^: _& Ilmiterm([2 13 14 P3],1,M)
/ i: ]) F/ L% r3 ^, e
: t( \) M* c- K- }4 V( iG=getlmis;
[tmin b]=feasp(G);
4 p4 H5 ^3 z( p/ y6 }
% k" R  @6 o6 v4 DP11 = dec2mat(G,b,P11);
P21 = dec2mat(G,b,P21);
P31 = dec2mat(G,b,P31);
P1 = [kron(eye(3), P11), kron(eye(3), P21); kron(eye(3), P21'), kron(eye(3), P31)]
! H' o- u  O1 r4 aP12 = dec2mat(G,b,P12);
4 n: y8 j3 Y( k+ j! t& u' RP22 = dec2mat(G,b,P22);
8 L  a3 h* r2 o0 eP32 = dec2mat(G,b,P32);
9 Y" @. C# c4 m6 \8 ?; SP2 = [kron(eye(3), P12), kron(eye(3), P22); kron(eye(3), P22'), kron(eye(3), P32)]
& R# {" o+ U& T/ h7 v  \P13 = dec2mat(G,b,P13);
P23 = dec2mat(G,b,P23);
' W* k# T2 x; [P33 = dec2mat(G,b,P33);
P3 = [kron(eye(3), P13), kron(eye(3), P23); kron(eye(3), P23'), kron(eye(3), P33)]
+ R3 z* O( W+ JZ(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)];
: @  }0 n$ j& @3 cQ = dec2mat(G,b,Q)
1 t& R0 B) U1 j% J( j. M6 X  L$ {mu1 = dec2mat(G,b,mu1)
mu2 = dec2mat(G,b,mu2)
mu3 = dec2mat(G,b,mu3)
s = dec2mat(G,b,s)

zaiyiaaaa

zzz.dan

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

IBB-EUT

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