matlab 问题

nkkopd

Error using lmiterm (line 281)
8 K0 t3 s: O( i9 Y2 f# U& s0 x7 p+ Glhs of LMI #1, block (1,1): term dimensions incompatible with other
! x  F- j" H/ F  q6 Nterms in row #1 or column #1
$ ^2 s  |4 l: u4 b  b
& ]+ f. [! a. FError in SCP_03 (line 103)
9 _, z# I- \# M1 [+ b! [. |/ }+ elmiterm([1 1 1 P1],-1,1)

clear
6 V6 K" f- l: d3 q2 v% w5 @, O, Y4 qclc

gamma = 0.95;

sigma1 = 2;
  h( q1 _9 c) H( u) ~5 ~/ Ksigma2 = 1.98;
% Y2 x7 z& X: n9 v4 ~  B& u; Hsigma3 = 1.86;
0 _; |: C+ |& D2 Q) d% q% ?+ csigma = sigma1+sigma2+sigma3;
. q0 S; O" P: Q* l6 p  i) C
9 q$ ?% U: d: N, M" R, h9 mgamma1 = 0.89;
+ m$ F5 ~. \! @8 Tgamma2 = 0.88;
% N: R6 ~3 L! e. _+ ?: z7 {gamma3 = 0.85;
Gamma = [gamma1, gamma1, gamma1; gamma2, gamma2, gamma2; gamma3, gamma3, gamma3];%3*3
) U) v8 \' w( _4 D- Y1 ~  t2 q
2 {9 `- G" J; mLambda = diag([2, 2]);%2*2
  T. k" j! T  ]7 U) u; M( L5 R% F
' y" r* g% T: j. r. F, v8 c. X$ KW = [kron(Gamma, Lambda), zeros(6,3); zeros(3,6), zeros(3)];%9*9
+ A& {! L4 f( c5 p
- u: d2 H' V* d8 {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

U = diag([1.68, 1.5, 1.68, 1.5, 1.68, 1.5, 0, 0, 0]);%花U1+U2   9*9
0 D) e& x, \1 [: }5 n! iL = diag([0.34, 0.34, 0.34, 0.34, 0.34, 0.34, 0, 0, 0]);%花L1+L2  9*9

- Q8 z  h9 J- m+ n3 ]6 ]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
0 O; [0 B* S3 fM = [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
  `! h: H2 O' V! L, J
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

( C" W  p. g! `3 AK1 = [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
6 u7 I& @0 E! C9 I& Y; EK2 = [-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
% N  i' X! C& ]# z4 t* p2 a$ k
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
2 t; D6 B5 R* g' w4 K
6 n  k9 D6 m9 a) ZE = [0.3;0.3;0.3];%3*1

" D2 O7 }8 E; y2 b- m) [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

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

/ o$ _( ^8 Y! t3 ]theta(1) = 2;
$ @/ g- A6 q4 I5 K& w; p3 \; s: ftheta(2) = 3;
! T% A7 d6 {& gtheta(3) = 1;
* h0 K0 N9 L+ g( A" N( _; I2 stheta(4) = 3;

2 O/ F7 V2 I  L8 ~1 |3 w2 ?" X
1 d+ v* ~+ N* r3 K
for k = 1:3
) H0 v' p  m1 ]    i = theta(k);
7 z2 S& `2 r, W) f. D4 \$ X- k( I    j = theta(k+1);
   
( ~$ p8 i* F3 O  O  Phi = zeros(3);
, }1 h5 [( ], B4 R: U  
for t = 1:3
     
    if (i-t) ~= 0;
7 F9 }  G# b+ k( V2 {3 M% w9 _          delta(t) = 0;
    elseif (i-t) == 0;
          delta(t) = 1;
2 f% @8 f0 y2 S- Q# ~7 l* Q# x% O         Phi(t,t) = delta(t);
, Z( [' k# Y6 ?0 ^0 W2 \    end   
  
end

' v" W2 ^8 _7 T+ s. [3 I+ H  barA = [A, zeros(6,3); (kron(Phi,eye(1)))*C, kron((eye(3)-Phi),eye(1))];%9*9
: E0 c: }' ]  p" d2 ^! S/ q. Z2 c  barB = [zeros(6,1), B; (kron(Phi,eye(1)))*E, (kron(Phi,eye(1)))*D];%9*4
0 f, t: J, c7 B8 c4 M  barC = [(kron(Phi,eye(1)))*C, kron((eye(3)-Phi),eye(1))];%3*9
+ J3 M) u4 T  s5 g8 X2 J" E- k  barD = [(kron(Phi,eye(1)))*E, kron(Phi,eye(1))*D];%3*4
end
% F6 F1 S* W8 S  % N=3  n=2  m=1 p=1 r=1 q=1  M=3
; ^9 o) \' b" C0 d  p2 H  setlmis([])
  P11 = lmivar(1,[2 1]);
  [& Y$ b$ X! W" x5 G  P21 = lmivar(2,[2 1]);
0 k$ a% |+ |  n  P31 = lmivar(1,[1 1]);
9 z" a4 H  ]! M: i# Y  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]);
! }7 |& s( [/ E' e  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]);
7 _; ~+ B9 g' T  V+ K  P23 = lmivar(2,[2 1]);
/ Y  ]1 B7 s0 A3 F  P33 = lmivar(1,[1 1]);
, T( c/ |: g5 G3 b  Z  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]);
  Q = lmivar(1,[9 1]);
  mu1 = lmivar(1,[1,1]);
3 D' L& o3 [. t% j. T  mu2 = lmivar(1,[1,1]);
  mu3 = lmivar(1,[1,1]);
# t( U; F3 s/ M$ ~/ W% h0 R; e  s = lmivar(1,[1,1]);
! ?, k, J, I7 O( Q+ S  barGamma = lmivar(1,[9 1]);
  barLambda = lmivar(1,[9 1]);
  
lmiterm([1 1 1 barLambda],2,1)
lmiterm([1 1 1 P1],-1,1)
lmiterm([1 1 1 mu1],-1,U+U')
8 T9 S2 u% X& A) M% ?lmiterm([1 1 1 mu2],-1,L+L')
* p. }, ^0 B% g* I/ S% Nlmiterm([1 1 1 s],1,N'*N)
lmiterm([1 1 1 0],H'*H)
4 P. |. V0 v5 Rlmiterm([1 1 5 mu1],1,U')
( Y& E2 i: O1 }7 u  o. ^lmiterm([1 1 7 mu2],1,L')
lmiterm([1 1 12 P2],barA'+W',1)
lmiterm([1 2 2 P2],barC'*K1',K1*barC)
* O: B# x6 Y8 Q# V6 t% almiterm([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')
lmiterm([1 2 2 mu2],-1,L+L')
( e* o. f/ l% P, _8 {) hlmiterm([1 2 6 mu1],1,U')
# w4 l3 H& g& D- n4 ilmiterm([1 2 8 mu2],1,L')
lmiterm([1 2 11 P2],barC'*K1',K1*barD)
lmiterm([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)
lmiterm([1 2 13 -Z2],-barC',1)
* r+ J$ k2 h& M7 i$ K! rlmiterm([1 3 3 mu3],-1,L')
lmiterm([1 3 9 mu3],1,L+L')
lmiterm([1 4 4 mu3],-1,L+L')
lmiterm([1 4 10 mu3],1,L')
lmiterm([1 5 5 mu1],-2,eye(9))
lmiterm([1 5 12 P2],1,1)
" j5 O8 d! F- o* M0 J% Clmiterm([1 6 6 mu1],-2,eye(9))
7 v! f+ t# C/ z  Q/ c) \. wlmiterm([1 6 13 P2],1,1)
lmiterm([1 7 7 Q],1,1)
lmiterm([1 7 7 mu2],-2,eye(9))
4 @# N1 x+ ~( h6 n- v+ _1 `: Qlmiterm([1 8 8 Q],1,1)
/ ^1 N+ W) C3 O1 n5 R# {. clmiterm([1 8 8 mu2],-2,eye(9))
! X" w' A8 [" ?+ y" v2 M- glmiterm([1 9 9 barGamma],2,1)
lmiterm([1 9 9 Q],-1,1)
! ~: `1 w3 ^, `2 u2 N# Y6 G+ B4 _7 |lmiterm([1 9 9 mu3],-2,eye(9))
/ v9 h5 Y" E2 D5 O. V8 Hlmiterm([1 9 12 P2],barXi',1)
. q" [- ^$ u5 I9 s' w' e2 l* ~lmiterm([1 10 10 Q],-1,1)
lmiterm([1 10 10 mu3],-2,eye(9))
; F% b7 ?/ `, S1 y8 Tlmiterm([1 10 13 P2],barXi',1)
lmiterm([1 11 11 P2],barD'*K1',K1*barD)
9 U, m$ V3 O; c6 s6 H6 Qlmiterm([1 11 11 P2],barD'*K2',K2*barD)
lmiterm([1 11 11 P2],barD'*K3',K3*barD)
$ K7 i+ o; m5 L0 }' Plmiterm([1 11 11 0],-(gamma^2)*eye(4))
lmiterm([1 11 12 P2],barB',1)
6 b  v: N% R0 y( w# w) c: R7 w  hlmiterm([1 11 13 P2],barB',1)
* L; f  e5 q% ]/ O0 Glmiterm([1 11 13 -Z2],-barD',1)
lmiterm([1 12 14 P2],1,M)
+ X& C2 G! o# J6 {) L7 O$ Olmiterm([1 13 14 P2],1,M)

% Z; z: G5 U2 Y  ~# \lmiterm([-1 1 1 P1],1,1)
$ Z2 S6 |) l1 U( v+ a  o2 A: {7 blmiterm([-2 1 1 P2],1,1)
/ G6 `. F1 t' {+ _3 M" y! Blmiterm([-3 1 1 P3],1,1)

& Y4 f1 ?% Q/ d  f: Q) Clmiterm([2 1 1 barLambda],2,1)
lmiterm([2 1 1 P2],-1,1)
lmiterm([2 1 1 mu1],-1,U+U')
& o7 S" _; m2 g1 ?: h, V6 _" blmiterm([2 1 1 mu2],-1,L+L')
lmiterm([2 1 1 s],1,N'*N)
lmiterm([2 1 1 0],H'*H)
3 |* S2 Z7 P- `$ Y5 Q1 ~lmiterm([2 1 5 mu1],1,U')
& f) K8 G: `; e3 Plmiterm([2 1 7 mu2],1,L')
lmiterm([2 1 12 P3],barA'+W',1)
lmiterm([2 2 2 P3],barC'*K1',K1*barC)
lmiterm([2 2 2 P3],barC'*K2',K2*barC)
3 v/ s( w9 s' Y( f7 Q$ Flmiterm([2 2 2 P3],barC'*K3',K3*barC)
lmiterm([2 2 2 P2],-1,1)
+ W) {6 I( K4 r8 tlmiterm([2 2 2 mu1],-1,U+U')
& T' U4 q4 d* g/ y3 c4 V. Glmiterm([2 2 2 mu2],-1,L+L')
3 u& a9 Y* m7 i- Y6 R4 u. slmiterm([2 2 6 mu1],1,U')
lmiterm([2 2 8 mu2],1,L')
lmiterm([2 2 11 P3],barC'*K1',K1*barD)
: ?* [; l- V1 P6 Y% v9 y4 Clmiterm([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')
lmiterm([2 4 4 mu3],-1,L+L')
9 ], [! B' C' h% n3 Mlmiterm([2 4 10 mu3],1,L')
lmiterm([2 5 5 mu1],-2,eye(9))
lmiterm([2 5 12 P3],1,1)
2 I; T2 F% K6 d  P3 Xlmiterm([2 6 6 mu1],-2,eye(9))
lmiterm([2 6 13 P3],1,1)
+ f' e9 V# t% b, slmiterm([2 7 7 Q],1,1)
3 A/ F1 ?/ b& s' P( jlmiterm([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)
: U, E* t4 n' Llmiterm([2 9 9 Q],-1,1)
( M2 \! k/ s" v! U! f7 zlmiterm([2 9 9 mu3],-2,eye(9))
' z3 H4 \' m/ m7 E4 f; [lmiterm([2 9 12 P3],barXi',1)
lmiterm([2 10 10 Q],-1,1)
lmiterm([2 10 10 mu3],-2,eye(9))
+ V  y/ g: D& X5 O0 F& f  B. tlmiterm([2 10 13 P3],barXi',1)
lmiterm([2 11 11 P3],barD'*K1',K1*barD)
7 f0 p$ p4 \. X9 ^  n' wlmiterm([2 11 11 P3],barD'*K2',K2*barD)
lmiterm([2 11 11 P3],barD'*K3',K3*barD)
3 z7 t6 `, [: h9 Glmiterm([2 11 11 0],-(gamma^2)*eye(4))
0 w' V6 q+ Y/ q1 `- `5 Xlmiterm([2 11 12 P3],barB',1)
" ]( T+ u4 f2 s& w( v2 klmiterm([2 11 13 P3],barB',1)
& f& u% X+ A6 x+ h- p" r" Z+ j" Nlmiterm([2 11 13 -Z3],-barD',1)
lmiterm([2 12 14 P3],1,M)
lmiterm([2 13 14 P3],1,M)
2 ~1 Q% U1 P' Z. B9 [9 w+ y
G=getlmis;
[tmin b]=feasp(G);

P11 = dec2mat(G,b,P11);
% t4 V( c6 f/ dP21 = dec2mat(G,b,P21);
4 Z2 V, h0 {. F* |$ x5 J4 y- M1 oP31 = dec2mat(G,b,P31);
P1 = [kron(eye(3), P11), kron(eye(3), P21); kron(eye(3), P21'), kron(eye(3), P31)]
8 I5 a8 J- Z, U. u* }  PP12 = dec2mat(G,b,P12);
1 @7 E0 W. ~1 n& H2 ?P22 = dec2mat(G,b,P22);
8 i# Q% y) b" N; x: g% e& R. s+ zP32 = dec2mat(G,b,P32);
3 i7 H" E  V8 x$ F# EP2 = [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);
P3 = [kron(eye(3), P13), kron(eye(3), P23); kron(eye(3), P23'), kron(eye(3), P33)]
Z(1) = dec2mat(G,b,Z(1));
- ?. [- c) l( W7 f3 S! }  y) QZ(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)];
* [) o7 A2 f) rQ = dec2mat(G,b,Q)
9 ~/ n  ?1 Y/ fmu1 = dec2mat(G,b,mu1)
) l7 _2 @: g; p0 W; ?  I- ^mu2 = dec2mat(G,b,mu2)
# |) F6 l! r* i; z/ g  Smu3 = dec2mat(G,b,mu3)
) N0 |  A' W( {* c7 gs = dec2mat(G,b,s)
5 {" W2 ~" P- z: [

zaiyiaaaa

zzz.dan

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

IBB-EUT

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