matlab 问题

nkkopd

Error using lmiterm (line 281)
3 ^* M" j6 Y3 w: \$ n7 o$ {lhs of LMI #1, block (1,1): term dimensions incompatible with other
terms in row #1 or column #1
% B: X5 {: A0 m1 ~  h
) u; E: Q. L- j4 R4 W4 b  q% IError in SCP_03 (line 103)
& K, i& C6 N* A+ J8 k' u7 klmiterm([1 1 1 P1],-1,1)
" ^0 N& ^/ F* T
) b9 ?$ I2 h, Xclear
% y+ m* l/ v0 d6 p! {9 mclc

gamma = 0.95;

9 D5 y4 C3 z) Q( ~0 Vsigma1 = 2;
4 t& D/ J3 H; P6 s0 asigma2 = 1.98;
1 B: @5 @7 g3 y$ P( jsigma3 = 1.86;
sigma = sigma1+sigma2+sigma3;

gamma1 = 0.89;
4 s7 B) O$ s/ qgamma2 = 0.88;
gamma3 = 0.85;
Gamma = [gamma1, gamma1, gamma1; gamma2, gamma2, gamma2; gamma3, gamma3, gamma3];%3*3
9 N# Q7 [. Q7 A; [9 V0 k: G
/ S% o, H1 {% zLambda = diag([2, 2]);%2*2

W = [kron(Gamma, Lambda), zeros(6,3); zeros(3,6), zeros(3)];%9*9

+ ~- A6 Q, L# Y% G5 }beta1 = 0.88;
7 {5 t2 e) s* w9 `  z# }: Bbeta2 = 0.85;
. m  s& T2 _$ U3 ^  H, f# c% Cbeta3 = 0.87;
Xi = diag([beta1, beta2, beta3]);%3*3
1 c! \( g6 y7 _* Q8 v( qbarXi = [kron(Xi, eye(2)), zeros(6,3); zeros(3,6), zeros(3)];%9*9
$ b6 h1 e! q$ C% v; A
- |8 u& }- s7 b7 c' o/ x1 rU = diag([1.68, 1.5, 1.68, 1.5, 1.68, 1.5, 0, 0, 0]);%花U1+U2   9*9
2 l$ l, m8 G: Y! F7 GL = diag([0.34, 0.34, 0.34, 0.34, 0.34, 0.34, 0, 0, 0]);%花L1+L2  9*9

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

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
K2 = [-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 n  t6 Y- p; K. sK3 =[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
+ d& T) q! B$ M/ N4 p( z! ~
H = [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 t) o3 T# q& nE = [0.3;0.3;0.3];%3*1
& n. d- m% f8 k4 O, i
" ]3 [: f# S3 q/ P2 l% W; oD = diag([0.5, 0.7, 0.4]);%3*3
' a/ K# f+ g1 P7 P/ f
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
' G# ~/ E! a. \% J9 s
/ _- L7 f  A8 \8 B: t  D' z) wB = [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
- ]7 r  O$ K0 T: L5 v
) c3 @3 W6 q& ?- i( u- Atheta(1) = 2;
theta(2) = 3;
theta(3) = 1;
; G, b) |; c7 Ntheta(4) = 3;

$ z6 N: u6 ?: D8 I3 c
2 O( @* R1 d; ~, I
for k = 1:3
$ O9 G& |# ~# v: E; ]' E& E) K: Z/ z    i = theta(k);
( S0 x9 T# `. w! ^    j = theta(k+1);
" Q; \1 q$ {0 l: F6 c  Q& ~   
( v$ E; D7 J. f8 P7 \8 O  Phi = zeros(3);
  
5 F7 \/ A! Q8 j8 W1 dfor t = 1:3
+ C+ \6 w+ r- D5 w, t     
# `  w9 M& P8 K! j0 W3 P    if (i-t) ~= 0;
7 \/ p& [6 S8 h- r/ O3 \/ [  |          delta(t) = 0;
% o8 v3 h/ |* b1 G' ~7 q; Y    elseif (i-t) == 0;
          delta(t) = 1;
! t- _7 h0 B; f/ |& ]         Phi(t,t) = delta(t);
$ x' e& D8 |! T3 M! `$ s5 R( h    end   
4 W5 }0 G5 P- S$ |0 |! ?; Z% w3 S  
9 V4 x0 n$ e  Q3 c* wend

  barA = [A, zeros(6,3); (kron(Phi,eye(1)))*C, kron((eye(3)-Phi),eye(1))];%9*9
1 o$ A- {& l. _& g  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
" p- W. z  p% M9 _/ g5 ?' v* Bend
7 R- x' x! K" K" p. o: R  % N=3  n=2  m=1 p=1 r=1 q=1  M=3
  setlmis([])
  P11 = lmivar(1,[2 1]);
  P21 = lmivar(2,[2 1]);
. c" k( o9 y! ~: A( C  P31 = lmivar(1,[1 1]);
0 }( j$ C5 L- c( n: X% k  P1 = lmivar(3,[kron(eye(3),P11),kron(eye(3),P21);kron(eye(3),P21'),kron(eye(3),P31)]);%9*9
( e5 D& {9 v7 E/ C7 @0 X3 {4 i  P12 = lmivar(1,[2 1]);
  P22 = lmivar(2,[2 1]);
  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
" h" y7 }4 h( T  k& d7 ?! t  P13 = lmivar(1,[2 1]);
# {( p( G  {* h6 H% W  P23 = lmivar(2,[2 1]);
  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
4 u, z2 k; |* r' c( {  Z1 = lmivar(2,[9 3]);
  Z2 = lmivar(2,[9 3]);
  Z3 = lmivar(2,[9 3]);
4 Z& M6 _0 e0 F9 W/ s% y  Q = lmivar(1,[9 1]);
  mu1 = lmivar(1,[1,1]);
  mu2 = lmivar(1,[1,1]);
  mu3 = lmivar(1,[1,1]);
4 E8 y$ o$ n( w  s = lmivar(1,[1,1]);
  barGamma = lmivar(1,[9 1]);
0 X& l2 Z: x$ O* H+ n8 {8 J7 v  barLambda = lmivar(1,[9 1]);
" D- K7 x/ l; A8 }6 R1 ^' P  
. i4 W+ j$ J+ ]lmiterm([1 1 1 barLambda],2,1)
lmiterm([1 1 1 P1],-1,1)
( ^, g8 J. U  `3 i6 Ylmiterm([1 1 1 mu1],-1,U+U')
, S; t$ y6 K& Rlmiterm([1 1 1 mu2],-1,L+L')
lmiterm([1 1 1 s],1,N'*N)
' W7 F! @- A4 |3 t* Ilmiterm([1 1 1 0],H'*H)
  @1 b- L5 q9 K) X# J7 ~  Olmiterm([1 1 5 mu1],1,U')
3 S+ b' `6 K/ w5 g7 G9 E9 n, ilmiterm([1 1 7 mu2],1,L')
lmiterm([1 1 12 P2],barA'+W',1)
8 U' P7 M& y3 e3 k/ l  j* Ilmiterm([1 2 2 P2],barC'*K1',K1*barC)
8 n* Q" {" l' Q' nlmiterm([1 2 2 P2],barC'*K2',K2*barC)
8 o" H/ s0 s" W& J. j% {9 `3 S6 l. rlmiterm([1 2 2 P2],barC'*K3',K3*barC)
: ~2 E6 l2 r1 M; {4 J  flmiterm([1 2 2 P1],-1,1)
lmiterm([1 2 2 mu1],-1,U+U')
* U9 ?7 @4 m, S; E9 {, a$ f# slmiterm([1 2 2 mu2],-1,L+L')
* K6 h( h$ M' V! M) Plmiterm([1 2 6 mu1],1,U')
$ z& J8 r+ h) W7 M( @+ flmiterm([1 2 8 mu2],1,L')
: e# u) s" |' |  H% z& l3 n. y4 olmiterm([1 2 11 P2],barC'*K1',K1*barD)
. k# n' m/ @5 w9 E- |2 N) jlmiterm([1 2 11 P2],barC'*K2',K2*barD)
: W4 a+ M* K3 c5 Y' Ilmiterm([1 2 11 P2],barC'*K3',K3*barD)
, B4 _$ |) ?) [1 almiterm([1 2 13 P2],barA'+W',1)
lmiterm([1 2 13 -Z2],-barC',1)
lmiterm([1 3 3 mu3],-1,L')
lmiterm([1 3 9 mu3],1,L+L')
! c& q3 D" Q4 w; ?5 c* m* K/ elmiterm([1 4 4 mu3],-1,L+L')
% |  Q1 P& o& m9 N8 Ilmiterm([1 4 10 mu3],1,L')
lmiterm([1 5 5 mu1],-2,eye(9))
lmiterm([1 5 12 P2],1,1)
9 m' t! J+ u' `0 u' J3 }lmiterm([1 6 6 mu1],-2,eye(9))
: }  Q8 t( D0 `% k0 Hlmiterm([1 6 13 P2],1,1)
+ o. L" Y# U  ?. w6 \: blmiterm([1 7 7 Q],1,1)
  l7 ]9 a+ r& I" s* o$ W2 Dlmiterm([1 7 7 mu2],-2,eye(9))
lmiterm([1 8 8 Q],1,1)
7 ?, f' D. j* \" b# U1 o$ a. M$ s/ glmiterm([1 8 8 mu2],-2,eye(9))
3 ^7 O% h2 o8 x" flmiterm([1 9 9 barGamma],2,1)
lmiterm([1 9 9 Q],-1,1)
6 d* H. k' l5 n0 e  b& o) G) Tlmiterm([1 9 9 mu3],-2,eye(9))
lmiterm([1 9 12 P2],barXi',1)
1 D( y* X% }+ |% b* glmiterm([1 10 10 Q],-1,1)
lmiterm([1 10 10 mu3],-2,eye(9))
' ~# Z0 r; j/ ^+ M7 r  F8 qlmiterm([1 10 13 P2],barXi',1)
lmiterm([1 11 11 P2],barD'*K1',K1*barD)
* _, A; |  Q( q$ {lmiterm([1 11 11 P2],barD'*K2',K2*barD)
3 X4 X1 T4 l* {+ vlmiterm([1 11 11 P2],barD'*K3',K3*barD)
lmiterm([1 11 11 0],-(gamma^2)*eye(4))
/ o3 W/ Q" u, ?2 ylmiterm([1 11 12 P2],barB',1)
lmiterm([1 11 13 P2],barB',1)
. _& k$ ?) Y# g. dlmiterm([1 11 13 -Z2],-barD',1)
lmiterm([1 12 14 P2],1,M)
lmiterm([1 13 14 P2],1,M)

lmiterm([-1 1 1 P1],1,1)
8 x7 j+ O8 u5 I4 Q" [0 R! b7 xlmiterm([-2 1 1 P2],1,1)
) J$ ~$ \" r+ y3 M+ }% V$ elmiterm([-3 1 1 P3],1,1)

) `) J  h6 D( m( e, Glmiterm([2 1 1 barLambda],2,1)
lmiterm([2 1 1 P2],-1,1)
4 r+ K0 X8 K* U; W) Ilmiterm([2 1 1 mu1],-1,U+U')
lmiterm([2 1 1 mu2],-1,L+L')
lmiterm([2 1 1 s],1,N'*N)
lmiterm([2 1 1 0],H'*H)
lmiterm([2 1 5 mu1],1,U')
lmiterm([2 1 7 mu2],1,L')
lmiterm([2 1 12 P3],barA'+W',1)
; e# g: W0 V  c0 y/ n5 rlmiterm([2 2 2 P3],barC'*K1',K1*barC)
% L3 s7 t" g7 L( D- z' flmiterm([2 2 2 P3],barC'*K2',K2*barC)
lmiterm([2 2 2 P3],barC'*K3',K3*barC)
% |- k4 k* G) j- E+ x5 }lmiterm([2 2 2 P2],-1,1)
, w- Z8 F* a* u9 S# Wlmiterm([2 2 2 mu1],-1,U+U')
lmiterm([2 2 2 mu2],-1,L+L')
+ c. Q4 L) L! U: r) ^. E- k( Elmiterm([2 2 6 mu1],1,U')
  H+ g! ^$ k" }lmiterm([2 2 8 mu2],1,L')
lmiterm([2 2 11 P3],barC'*K1',K1*barD)
+ ]* j+ `, W( X) d4 _lmiterm([2 2 11 P3],barC'*K2',K2*barD)
lmiterm([2 2 11 P3],barC'*K3',K3*barD)
; S5 f& g. U: q3 Ilmiterm([2 2 13 P3],barA'+W',1)
lmiterm([2 2 13 -Z3],-barC',1)
lmiterm([2 3 3 mu3],-1,L')
3 o  t# {6 S; T( |; D7 k9 k  Rlmiterm([2 3 9 mu3],1,L+L')
2 n0 V* ]( z" o: |  X' xlmiterm([2 4 4 mu3],-1,L+L')
! z% S( Y1 D+ o! I! `lmiterm([2 4 10 mu3],1,L')
lmiterm([2 5 5 mu1],-2,eye(9))
6 G1 b8 V) {) h' w# }lmiterm([2 5 12 P3],1,1)
& g6 g; \6 a% c! Ilmiterm([2 6 6 mu1],-2,eye(9))
3 v# @) n, G) o' A4 d. glmiterm([2 6 13 P3],1,1)
7 r5 w: Z" I$ f+ _* U, Llmiterm([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))
lmiterm([2 9 9 barGamma],2,1)
lmiterm([2 9 9 Q],-1,1)
# o9 V$ i1 g& m- D4 O, d+ O3 Blmiterm([2 9 9 mu3],-2,eye(9))
) f/ t5 O# P0 O5 F) Q4 Olmiterm([2 9 12 P3],barXi',1)
lmiterm([2 10 10 Q],-1,1)
' p) G# H8 U3 d; a* `, xlmiterm([2 10 10 mu3],-2,eye(9))
3 E# c0 T! s! U* @6 g7 p  Q' Qlmiterm([2 10 13 P3],barXi',1)
. k* u4 ?- [) c9 T9 Flmiterm([2 11 11 P3],barD'*K1',K1*barD)
9 r6 D* v3 R: \5 f5 nlmiterm([2 11 11 P3],barD'*K2',K2*barD)
lmiterm([2 11 11 P3],barD'*K3',K3*barD)
. s* P7 z" [: j2 z; x2 m: T/ clmiterm([2 11 11 0],-(gamma^2)*eye(4))
8 @2 o7 E3 P, s1 [lmiterm([2 11 12 P3],barB',1)
5 S! x  W5 @! F3 o* u, X7 O- dlmiterm([2 11 13 P3],barB',1)
! C# N: @0 n/ J3 vlmiterm([2 11 13 -Z3],-barD',1)
lmiterm([2 12 14 P3],1,M)
' Z: d( A) u/ Rlmiterm([2 13 14 P3],1,M)
: f& N: D. x! x" t
G=getlmis;
[tmin b]=feasp(G);

P11 = dec2mat(G,b,P11);
2 S3 H" e* z5 @/ W3 tP21 = dec2mat(G,b,P21);
P31 = dec2mat(G,b,P31);
  A8 t9 |5 }1 u0 X' cP1 = [kron(eye(3), P11), kron(eye(3), P21); kron(eye(3), P21'), kron(eye(3), P31)]
P12 = dec2mat(G,b,P12);
P22 = dec2mat(G,b,P22);
P32 = dec2mat(G,b,P32);
P2 = [kron(eye(3), P12), kron(eye(3), P22); kron(eye(3), P22'), kron(eye(3), P32)]
P13 = dec2mat(G,b,P13);
P23 = dec2mat(G,b,P23);
P33 = dec2mat(G,b,P33);
: X: s! v, p" I& E/ XP3 = [kron(eye(3), P13), kron(eye(3), P23); kron(eye(3), P23'), kron(eye(3), P33)]
0 L. M) A  c8 m% v: i& q8 X) 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)];
Q = dec2mat(G,b,Q)
. b- l& F( t: p% f. F# _  ymu1 = dec2mat(G,b,mu1)
/ Y% f  K3 J+ Y4 e9 ?+ p3 ~7 D. B- ~mu2 = dec2mat(G,b,mu2)
mu3 = dec2mat(G,b,mu3)
) P! I+ X4 g* ^( ?& ps = dec2mat(G,b,s)
* y1 `. ~7 ]+ m5 f( J  W  Q

zaiyiaaaa

zzz.dan

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

IBB-EUT

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