matlab 问题

nkkopd

Error using lmiterm (line 281)
; r+ R; L& w: `0 Hlhs of LMI #1, block (1,1): term dimensions incompatible with other
terms in row #1 or column #1

1 ^' K/ V% A& r! r3 _& nError in SCP_03 (line 103)
% t- C' _- z. R! L) B7 j$ Vlmiterm([1 1 1 P1],-1,1)
- f3 i# z4 @2 P# I9 g+ J- t6 n
clear
' m+ i1 Q  e1 D7 Iclc
" g# ^3 F3 A5 H7 E+ n5 S# Q% A; t1 s. ^
gamma = 0.95;

sigma1 = 2;
sigma2 = 1.98;
6 [* l3 q; I1 Tsigma3 = 1.86;
& I# ]! k7 a5 V  `sigma = sigma1+sigma2+sigma3;

gamma1 = 0.89;
gamma2 = 0.88;
gamma3 = 0.85;
! Y/ n( {7 A% f" a2 _. L- b7 x, y$ xGamma = [gamma1, gamma1, gamma1; gamma2, gamma2, gamma2; gamma3, gamma3, gamma3];%3*3

$ D$ b8 D3 b/ }* E) GLambda = diag([2, 2]);%2*2
5 N6 y. f6 _3 F* P8 ]
W = [kron(Gamma, Lambda), zeros(6,3); zeros(3,6), zeros(3)];%9*9

beta1 = 0.88;
/ G% e5 J. k! ~/ A$ I3 S4 \" ?( w! Wbeta2 = 0.85;
. s' Q5 X) Y) E1 {# I, O. ebeta3 = 0.87;
9 ?7 n; \: @) pXi = diag([beta1, beta2, beta3]);%3*3
barXi = [kron(Xi, eye(2)), zeros(6,3); zeros(3,6), zeros(3)];%9*9
) w& q* u% H$ @+ k' i
# |3 M' `( v* N- \" qU = diag([1.68, 1.5, 1.68, 1.5, 1.68, 1.5, 0, 0, 0]);%花U1+U2   9*9
L = diag([0.34, 0.34, 0.34, 0.34, 0.34, 0.34, 0, 0, 0]);%花L1+L2  9*9

4 i* T2 Y! `7 ?7 o* q. eN = [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
9 _* f1 w) i+ }( k! W1 o* J& v
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
; q$ c2 X" g3 v6 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
/ X6 U0 m3 b& G" C$ V1 Z: b
7 G9 W# x# O" n0 sH = [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

D = diag([0.5, 0.7, 0.4]);%3*3
/ F/ {+ a7 E; J: \* k  u* T
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
! I5 k& C# t4 a/ A& N1 j+ u
0 M0 f" B, E9 qB = [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
; U' n& y5 S4 |& G
# c# F* S: @# K8 r, l; K0 ^+ s0 rtheta(1) = 2;
, @  z- b; Y" N. f. q* C# ptheta(2) = 3;
theta(3) = 1;
theta(4) = 3;
4 I! L" A3 w; E; n+ ]( Y! Z


for k = 1:3
; P, j0 R5 X4 d8 l9 f4 b0 Y; d* E' X    i = theta(k);
3 l, m1 T# K& G( I+ W2 Q8 x7 Q    j = theta(k+1);
+ t% ~( N( _6 N6 q   
  Phi = zeros(3);
  
for t = 1:3
- g' ]1 u3 u9 G5 I/ [; O: F6 h6 r     
3 j, P8 M4 h. u: |+ g' b! S    if (i-t) ~= 0;
          delta(t) = 0;
' B. ?4 Q6 e& I% g4 F    elseif (i-t) == 0;
% L- S8 j7 i; a3 U          delta(t) = 1;
         Phi(t,t) = delta(t);
    end   
  
2 t' g: z' Q: L; M2 cend

  barA = [A, zeros(6,3); (kron(Phi,eye(1)))*C, kron((eye(3)-Phi),eye(1))];%9*9
  ~. V: q% |5 N+ ^4 l2 T& i- \  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
4 z3 w' I8 Z  }6 ?' p! n  barD = [(kron(Phi,eye(1)))*E, kron(Phi,eye(1))*D];%3*4
end
# r1 M* y8 C! T  % N=3  n=2  m=1 p=1 r=1 q=1  M=3
  a* L  t, o: v2 B  setlmis([])
  P11 = lmivar(1,[2 1]);
8 O( n0 L0 @! I$ t, r  _+ R" y' S/ u  P21 = lmivar(2,[2 1]);
0 Q) b, T9 j: B  U9 P$ v  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
3 W4 w4 `4 ^9 K+ g) y  P12 = lmivar(1,[2 1]);
  P22 = lmivar(2,[2 1]);
; ~2 w, i" Z* |' E$ i  `  P32 = lmivar(1,[1 1]);
* o3 C3 U* x+ d% O' S0 d  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]);
+ n9 L6 D$ S( Q, J  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]);
9 `+ l+ L; H# z2 V( }5 k6 @) R  Q = lmivar(1,[9 1]);
/ G- F3 L% Y, m. Q# K4 n  mu1 = lmivar(1,[1,1]);
9 |% g, d! T7 {  mu2 = lmivar(1,[1,1]);
2 z/ B5 L. X0 b- b$ T  `- u; s  mu3 = lmivar(1,[1,1]);
5 f. ?, I! I; d1 m- K  s = lmivar(1,[1,1]);
  barGamma = lmivar(1,[9 1]);
  barLambda = lmivar(1,[9 1]);
  
lmiterm([1 1 1 barLambda],2,1)
lmiterm([1 1 1 P1],-1,1)
; C) t9 v) @, {. p% [1 @lmiterm([1 1 1 mu1],-1,U+U')
0 \1 q% t- x* n9 Y6 p! Blmiterm([1 1 1 mu2],-1,L+L')
* r- s( c& j6 t. F4 f8 ~lmiterm([1 1 1 s],1,N'*N)
lmiterm([1 1 1 0],H'*H)
lmiterm([1 1 5 mu1],1,U')
lmiterm([1 1 7 mu2],1,L')
lmiterm([1 1 12 P2],barA'+W',1)
lmiterm([1 2 2 P2],barC'*K1',K1*barC)
lmiterm([1 2 2 P2],barC'*K2',K2*barC)
$ \( B! h2 i% u. i& @0 wlmiterm([1 2 2 P2],barC'*K3',K3*barC)
  l# S( y" ]( s  y+ D4 m2 wlmiterm([1 2 2 P1],-1,1)
. o2 E$ z/ J# g' tlmiterm([1 2 2 mu1],-1,U+U')
+ s* N1 p6 X6 T' n6 e( ?lmiterm([1 2 2 mu2],-1,L+L')
lmiterm([1 2 6 mu1],1,U')
lmiterm([1 2 8 mu2],1,L')
lmiterm([1 2 11 P2],barC'*K1',K1*barD)
" e. g2 Z2 L2 x; b, F, m# p1 M4 Almiterm([1 2 11 P2],barC'*K2',K2*barD)
lmiterm([1 2 11 P2],barC'*K3',K3*barD)
lmiterm([1 2 13 P2],barA'+W',1)
( t. }: w$ C7 |8 T( n4 Zlmiterm([1 2 13 -Z2],-barC',1)
lmiterm([1 3 3 mu3],-1,L')
lmiterm([1 3 9 mu3],1,L+L')
lmiterm([1 4 4 mu3],-1,L+L')
6 f8 n+ Z$ E0 |# ?% f7 Ulmiterm([1 4 10 mu3],1,L')
7 s* M+ o* L7 W, wlmiterm([1 5 5 mu1],-2,eye(9))
lmiterm([1 5 12 P2],1,1)
lmiterm([1 6 6 mu1],-2,eye(9))
6 l4 y& H8 v  Z; H; A1 d$ Q' almiterm([1 6 13 P2],1,1)
$ Z9 I$ G% L& _$ dlmiterm([1 7 7 Q],1,1)
lmiterm([1 7 7 mu2],-2,eye(9))
lmiterm([1 8 8 Q],1,1)
lmiterm([1 8 8 mu2],-2,eye(9))
8 W$ n5 Q$ H6 u$ {+ b3 Olmiterm([1 9 9 barGamma],2,1)
lmiterm([1 9 9 Q],-1,1)
" v( k) u; r/ T$ X$ ]0 m  S  U. t) ]lmiterm([1 9 9 mu3],-2,eye(9))
lmiterm([1 9 12 P2],barXi',1)
1 |/ j& O  ]& N6 K1 t  Slmiterm([1 10 10 Q],-1,1)
+ d3 w7 M2 V+ u3 l! K( N1 ?; elmiterm([1 10 10 mu3],-2,eye(9))
' `4 n/ s( ^" F" }& olmiterm([1 10 13 P2],barXi',1)
lmiterm([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)
  Q8 i, j4 {1 c; T* d. Mlmiterm([1 11 11 0],-(gamma^2)*eye(4))
lmiterm([1 11 12 P2],barB',1)
- P6 z  ?1 F3 b. Mlmiterm([1 11 13 P2],barB',1)
lmiterm([1 11 13 -Z2],-barD',1)
5 w# q, {: n8 U' R& W2 J1 Hlmiterm([1 12 14 P2],1,M)
7 N( A) U2 k/ t6 P1 p: Y5 `lmiterm([1 13 14 P2],1,M)

0 y+ E5 d) ?+ X0 R8 G4 T& l# xlmiterm([-1 1 1 P1],1,1)
; J, T  r( a/ _* c0 v; E! nlmiterm([-2 1 1 P2],1,1)
7 j- j0 _% p0 {3 Ylmiterm([-3 1 1 P3],1,1)

lmiterm([2 1 1 barLambda],2,1)
lmiterm([2 1 1 P2],-1,1)
& g% C6 l! i9 c9 ^9 hlmiterm([2 1 1 mu1],-1,U+U')
lmiterm([2 1 1 mu2],-1,L+L')
7 m7 D! T4 m5 s" ulmiterm([2 1 1 s],1,N'*N)
+ @2 Q' a- y5 M/ v+ glmiterm([2 1 1 0],H'*H)
5 p+ k! ~: k8 D$ C- |- ?lmiterm([2 1 5 mu1],1,U')
% P- _4 d% B1 b2 l8 rlmiterm([2 1 7 mu2],1,L')
lmiterm([2 1 12 P3],barA'+W',1)
& g2 _/ A" m5 Y' b+ V7 `: rlmiterm([2 2 2 P3],barC'*K1',K1*barC)
lmiterm([2 2 2 P3],barC'*K2',K2*barC)
3 C* @( w: S+ [1 B/ slmiterm([2 2 2 P3],barC'*K3',K3*barC)
lmiterm([2 2 2 P2],-1,1)
9 F  O+ U  f( Zlmiterm([2 2 2 mu1],-1,U+U')
lmiterm([2 2 2 mu2],-1,L+L')
% p# P% f) o, a7 Flmiterm([2 2 6 mu1],1,U')
& D" a- j4 @$ Dlmiterm([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')
lmiterm([2 3 9 mu3],1,L+L')
" I+ @0 h$ ?& a; w1 a0 p+ Clmiterm([2 4 4 mu3],-1,L+L')
( ~& ?# J. v4 i$ v; ~lmiterm([2 4 10 mu3],1,L')
5 l8 X# N6 A$ ]8 dlmiterm([2 5 5 mu1],-2,eye(9))
* j4 w( `, e8 l! E" Wlmiterm([2 5 12 P3],1,1)
lmiterm([2 6 6 mu1],-2,eye(9))
lmiterm([2 6 13 P3],1,1)
2 I9 q! q* U' e' Vlmiterm([2 7 7 Q],1,1)
' ~. p" ^7 c( }lmiterm([2 7 7 mu2],-2,eye(9))
1 a& L. N/ @" x2 Slmiterm([2 8 8 Q],1,1)
* v; w. l+ p# r- q1 m' ~- Xlmiterm([2 8 8 mu2],-2,eye(9))
" r( y" w% O" r9 {0 f3 o* {) P* s: t- olmiterm([2 9 9 barGamma],2,1)
: ]6 A6 N% {" ?. A5 j  `lmiterm([2 9 9 Q],-1,1)
lmiterm([2 9 9 mu3],-2,eye(9))
+ J5 h0 ]7 W  _. ]7 almiterm([2 9 12 P3],barXi',1)
. W7 k7 Z% \' {6 K8 olmiterm([2 10 10 Q],-1,1)
lmiterm([2 10 10 mu3],-2,eye(9))
lmiterm([2 10 13 P3],barXi',1)
lmiterm([2 11 11 P3],barD'*K1',K1*barD)
lmiterm([2 11 11 P3],barD'*K2',K2*barD)
. v) u0 J+ m$ I% u9 Blmiterm([2 11 11 P3],barD'*K3',K3*barD)
lmiterm([2 11 11 0],-(gamma^2)*eye(4))
9 E4 U0 T( {+ w4 j3 vlmiterm([2 11 12 P3],barB',1)
  f9 ]1 y! ]% p: r' e9 S( [9 hlmiterm([2 11 13 P3],barB',1)
) K8 S5 l* X+ U- x/ jlmiterm([2 11 13 -Z3],-barD',1)
lmiterm([2 12 14 P3],1,M)
4 G/ N$ y" F% Llmiterm([2 13 14 P3],1,M)

G=getlmis;
; S& g. w, b) E& W& s! r[tmin b]=feasp(G);
2 x4 A/ R& a" q% y% `0 i8 K
P11 = dec2mat(G,b,P11);
' b% w$ H5 E4 F1 oP21 = dec2mat(G,b,P21);
7 k$ G7 i% J' h. }" EP31 = dec2mat(G,b,P31);
! ]  P) W" x$ M$ G. w$ dP1 = [kron(eye(3), P11), kron(eye(3), P21); kron(eye(3), P21'), kron(eye(3), P31)]
2 h! V) G" _, r7 ^/ eP12 = dec2mat(G,b,P12);
' }; f: W1 H+ r# d* wP22 = dec2mat(G,b,P22);
P32 = dec2mat(G,b,P32);
/ P# P  R) z7 ^  {7 X. I+ g* L! D5 MP2 = [kron(eye(3), P12), kron(eye(3), P22); kron(eye(3), P22'), kron(eye(3), P32)]
# k5 ^$ G0 u7 B; AP13 = dec2mat(G,b,P13);
/ F" }' ]9 W4 r' G6 A8 Y! YP23 = dec2mat(G,b,P23);
- ~2 t+ j0 I) A3 a* K; T3 a+ JP33 = dec2mat(G,b,P33);
P3 = [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));
8 O: p4 N7 X( S3 n6 g+ [3 EZ(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)];
, }7 N  z' r+ g0 D" H* k% P- F* Z: }* C4 hQ = dec2mat(G,b,Q)
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不一致