matlab 问题

nkkopd

Error using lmiterm (line 281)
lhs of LMI #1, block (1,1): term dimensions incompatible with other
terms in row #1 or column #1
9 q$ T9 U; e6 O3 I0 t
; N0 z& I$ W1 {/ m4 ]Error in SCP_03 (line 103)
7 I, l! q6 s: E( w- ^3 Slmiterm([1 1 1 P1],-1,1)
. ]1 z; Z7 C: L/ U" Z1 R/ X( @+ V
) J" _! H% [7 s" u9 s: P& Vclear
clc

& B; c0 O: U+ ~% B: ]0 Zgamma = 0.95;
; o: ^+ s! ?7 |' c0 \
sigma1 = 2;
sigma2 = 1.98;
sigma3 = 1.86;
/ `; k5 t* u+ d! j! X, x6 `. [sigma = sigma1+sigma2+sigma3;

) C: C* |1 \2 W) A) m9 ?" Jgamma1 = 0.89;
6 O; Z& l& l' Y# p. a& ^5 Jgamma2 = 0.88;
gamma3 = 0.85;
+ I6 ^  p. z* KGamma = [gamma1, gamma1, gamma1; gamma2, gamma2, gamma2; gamma3, gamma3, gamma3];%3*3

Lambda = diag([2, 2]);%2*2
  X4 B; ]5 t/ _
+ M  m, y% q0 N9 n( A; g2 EW = [kron(Gamma, Lambda), zeros(6,3); zeros(3,6), zeros(3)];%9*9
: o* L5 t, {  b, N
beta1 = 0.88;
2 X. f. X: ]* w0 _, g2 Z( Sbeta2 = 0.85;
beta3 = 0.87;
7 n. a* c+ p0 B  K/ ^% IXi = 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
L = diag([0.34, 0.34, 0.34, 0.34, 0.34, 0.34, 0, 0, 0]);%花L1+L2  9*9

# l6 r/ v. O" cN = [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
( e5 ~% W' [# W- L) ~/ Y, cK3 =[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
3 z, @" \- E6 d- V: `# M
7 M, t: v: y1 O) I' e7 r/ w# OH = [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

- c) }: o/ e$ M" E! o5 p+ x2 E" jE = [0.3;0.3;0.3];%3*1

5 U5 k' {$ N) z% \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

theta(1) = 2;
theta(2) = 3;
theta(3) = 1;
, T8 h) m7 O; o. s# {: Ktheta(4) = 3;



for k = 1:3
    i = theta(k);
; ~- Q  O3 ~2 S/ D+ @4 C# x    j = theta(k+1);
% A/ P9 M0 R4 B   
0 @% I' S" ^# {. ]+ Y- j  Phi = zeros(3);
  
for t = 1:3
/ V. y/ S1 j7 P9 E  f) N     
# n) f6 Y& b% k- q& S    if (i-t) ~= 0;
' r, O5 y! G$ V! s; I1 A6 C          delta(t) = 0;
! ?0 L% f/ I! {- Q    elseif (i-t) == 0;
1 \( f* v+ j: D4 }          delta(t) = 1;
         Phi(t,t) = delta(t);
1 e: O2 t9 N% S) F% |* Q" ~    end   
  
! J- A) ]  l$ j- R$ F6 Cend

  barA = [A, zeros(6,3); (kron(Phi,eye(1)))*C, kron((eye(3)-Phi),eye(1))];%9*9
& I4 S. H" F5 E% |9 O& r: e  barB = [zeros(6,1), B; (kron(Phi,eye(1)))*E, (kron(Phi,eye(1)))*D];%9*4
7 w. `2 n& @7 l" u: l  R  barC = [(kron(Phi,eye(1)))*C, kron((eye(3)-Phi),eye(1))];%3*9
4 c! l, e. n7 Q) _  barD = [(kron(Phi,eye(1)))*E, kron(Phi,eye(1))*D];%3*4
5 r. ~- B4 ?: H5 Q) K1 [8 i+ rend
  % 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]);
  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
$ p" i( b5 n: B) @) B3 i2 w4 l  P12 = lmivar(1,[2 1]);
  P22 = lmivar(2,[2 1]);
  P32 = lmivar(1,[1 1]);
  S0 J% g$ I$ |! @( y  P2 = lmivar(3,[kron(eye(3),P12),kron(eye(3),P22);kron(eye(3),P22'),kron(eye(3),P32)]);%9*9
5 r# Q% B* u; W2 N) T2 ~( b, g) z  P13 = lmivar(1,[2 1]);
4 U" u* T: y2 O3 r/ U/ m% N( Q  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
/ p* \) T  M9 {- ~. i  Z1 = lmivar(2,[9 3]);
; Q  T, }& c' Y1 `; x  Z2 = lmivar(2,[9 3]);
  Z3 = lmivar(2,[9 3]);
0 W: Y; d$ {4 `/ v& a" s- h  Q = lmivar(1,[9 1]);
  mu1 = lmivar(1,[1,1]);
  mu2 = lmivar(1,[1,1]);
  mu3 = lmivar(1,[1,1]);
  s = lmivar(1,[1,1]);
+ v* w4 ~) X; I9 n" d" e; p' _  barGamma = lmivar(1,[9 1]);
  barLambda = lmivar(1,[9 1]);
  
lmiterm([1 1 1 barLambda],2,1)
# E* g4 j# ^. @lmiterm([1 1 1 P1],-1,1)
8 j2 D7 }! Q2 n. E: V6 flmiterm([1 1 1 mu1],-1,U+U')
! _2 Z, n/ E2 W0 n% _8 r! ~& Llmiterm([1 1 1 mu2],-1,L+L')
lmiterm([1 1 1 s],1,N'*N)
3 w3 q! z+ o+ `8 r# _6 T# y- mlmiterm([1 1 1 0],H'*H)
) b7 L! _1 J. v; vlmiterm([1 1 5 mu1],1,U')
' \. G" r/ g: p! L- Q  clmiterm([1 1 7 mu2],1,L')
. z9 T* I7 J! B  D/ s6 p; c9 blmiterm([1 1 12 P2],barA'+W',1)
0 P/ U) m4 l9 L: ^2 M1 \# F, g4 Ylmiterm([1 2 2 P2],barC'*K1',K1*barC)
lmiterm([1 2 2 P2],barC'*K2',K2*barC)
3 y. }) l" K2 i/ [lmiterm([1 2 2 P2],barC'*K3',K3*barC)
lmiterm([1 2 2 P1],-1,1)
lmiterm([1 2 2 mu1],-1,U+U')
6 r' @1 E* n# j& B% b$ b" \. jlmiterm([1 2 2 mu2],-1,L+L')
lmiterm([1 2 6 mu1],1,U')
lmiterm([1 2 8 mu2],1,L')
- q6 p8 z4 Y7 }- v4 n& Klmiterm([1 2 11 P2],barC'*K1',K1*barD)
lmiterm([1 2 11 P2],barC'*K2',K2*barD)
1 J2 Z: |( \7 P0 Y2 B2 z8 P/ rlmiterm([1 2 11 P2],barC'*K3',K3*barD)
$ y0 U8 |5 Z( L) O' @4 Vlmiterm([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')
lmiterm([1 4 4 mu3],-1,L+L')
lmiterm([1 4 10 mu3],1,L')
lmiterm([1 5 5 mu1],-2,eye(9))
3 h9 F% F2 z% |6 Z* `lmiterm([1 5 12 P2],1,1)
lmiterm([1 6 6 mu1],-2,eye(9))
! h1 {. c/ ^2 S; Q4 s5 g/ klmiterm([1 6 13 P2],1,1)
6 ]5 M1 S) w: A4 Z2 Z. clmiterm([1 7 7 Q],1,1)
1 o" v0 p; g3 T0 flmiterm([1 7 7 mu2],-2,eye(9))
5 ~  {6 d6 z2 Q5 Tlmiterm([1 8 8 Q],1,1)
lmiterm([1 8 8 mu2],-2,eye(9))
- e* a* S( P0 w* Jlmiterm([1 9 9 barGamma],2,1)
3 `6 A; U" I% G; M% P  l9 ylmiterm([1 9 9 Q],-1,1)
lmiterm([1 9 9 mu3],-2,eye(9))
7 z0 q$ b6 F+ llmiterm([1 9 12 P2],barXi',1)
5 }/ m. Y9 f# qlmiterm([1 10 10 Q],-1,1)
/ j6 Q2 _2 C; i' M1 Mlmiterm([1 10 10 mu3],-2,eye(9))
+ a, F" w; _4 V3 g2 llmiterm([1 10 13 P2],barXi',1)
5 x! m! l8 V1 D0 Q) llmiterm([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)
lmiterm([1 11 11 0],-(gamma^2)*eye(4))
lmiterm([1 11 12 P2],barB',1)
lmiterm([1 11 13 P2],barB',1)
! \+ ?9 x- `$ A+ G) Blmiterm([1 11 13 -Z2],-barD',1)
; ^4 z0 v+ [: }) V0 K/ h, nlmiterm([1 12 14 P2],1,M)
lmiterm([1 13 14 P2],1,M)

lmiterm([-1 1 1 P1],1,1)
# A, d$ b# l$ v9 _lmiterm([-2 1 1 P2],1,1)
lmiterm([-3 1 1 P3],1,1)
1 N' v- W' |  d6 I3 X
lmiterm([2 1 1 barLambda],2,1)
4 S2 P) l, v- f) |2 jlmiterm([2 1 1 P2],-1,1)
' h3 ]# q  ~# n. n& qlmiterm([2 1 1 mu1],-1,U+U')
5 p) l! ?3 }! y! `9 E5 x& @lmiterm([2 1 1 mu2],-1,L+L')
5 o3 z* I/ {# T- J; D: b8 elmiterm([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)
lmiterm([2 2 2 P3],barC'*K1',K1*barC)
lmiterm([2 2 2 P3],barC'*K2',K2*barC)
lmiterm([2 2 2 P3],barC'*K3',K3*barC)
lmiterm([2 2 2 P2],-1,1)
; c7 I' s8 H, w( I; p4 Mlmiterm([2 2 2 mu1],-1,U+U')
lmiterm([2 2 2 mu2],-1,L+L')
9 b1 u& b$ A  x- g8 P1 almiterm([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)
7 i% e; k; _# L( X4 c( @lmiterm([2 2 11 P3],barC'*K3',K3*barD)
6 Y  ?- v( H( y9 D7 Dlmiterm([2 2 13 P3],barA'+W',1)
lmiterm([2 2 13 -Z3],-barC',1)
4 ^  q) R+ ?: X. Z: Jlmiterm([2 3 3 mu3],-1,L')
: ^4 c; i- R( X! M" Klmiterm([2 3 9 mu3],1,L+L')
lmiterm([2 4 4 mu3],-1,L+L')
lmiterm([2 4 10 mu3],1,L')
lmiterm([2 5 5 mu1],-2,eye(9))
lmiterm([2 5 12 P3],1,1)
9 M- O. m0 ^7 B: x# }lmiterm([2 6 6 mu1],-2,eye(9))
lmiterm([2 6 13 P3],1,1)
lmiterm([2 7 7 Q],1,1)
* R5 x9 B6 i- f5 L3 wlmiterm([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)
% @+ L" ~9 m: N# K/ }  n8 ]lmiterm([2 9 9 Q],-1,1)
/ j; J! t% |9 ^# \$ J1 n4 m6 r4 ~lmiterm([2 9 9 mu3],-2,eye(9))
8 b0 [' B$ H' C& V4 ylmiterm([2 9 12 P3],barXi',1)
lmiterm([2 10 10 Q],-1,1)
% N' P6 q' _! q# {$ d, G: i5 o" u  ~lmiterm([2 10 10 mu3],-2,eye(9))
lmiterm([2 10 13 P3],barXi',1)
6 {7 n1 a# k5 t6 q0 d: V+ t6 _0 M, nlmiterm([2 11 11 P3],barD'*K1',K1*barD)
  p; v" F4 j, f) F5 Zlmiterm([2 11 11 P3],barD'*K2',K2*barD)
lmiterm([2 11 11 P3],barD'*K3',K3*barD)
* [0 w% P: m( jlmiterm([2 11 11 0],-(gamma^2)*eye(4))
lmiterm([2 11 12 P3],barB',1)
lmiterm([2 11 13 P3],barB',1)
( H, y4 q& q( J+ y% Z$ U! Blmiterm([2 11 13 -Z3],-barD',1)
lmiterm([2 12 14 P3],1,M)
lmiterm([2 13 14 P3],1,M)
  A: x" A4 x; \+ d
G=getlmis;
[tmin b]=feasp(G);

P11 = dec2mat(G,b,P11);
P21 = dec2mat(G,b,P21);
" T/ k/ C* L! Q5 E2 F7 i" [% DP31 = dec2mat(G,b,P31);
! X/ E2 b. r5 q/ SP1 = [kron(eye(3), P11), kron(eye(3), P21); kron(eye(3), P21'), kron(eye(3), P31)]
$ N' m$ w$ ^5 Z8 ]4 a! i) tP12 = dec2mat(G,b,P12);
6 V: E& Q7 ?2 j4 iP22 = dec2mat(G,b,P22);
% f/ i6 K3 g4 H+ y7 EP32 = dec2mat(G,b,P32);
- h3 l, R$ c+ M9 UP2 = [kron(eye(3), P12), kron(eye(3), P22); kron(eye(3), P22'), kron(eye(3), P32)]
" X* r5 c7 K0 n. k1 d5 HP13 = dec2mat(G,b,P13);
P23 = dec2mat(G,b,P23);
5 t: `" t3 g" n" _$ B8 |1 e0 mP33 = dec2mat(G,b,P33);
( x9 M2 G4 C- \2 v  iP3 = [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)];
) h( @2 [3 o6 [; M# ?( L3 Q9 QbarLambda = [kron(sigma*eye(3), Lambda'*(P(1i).P(i))*Lambda), zeros(6,3); zeros(3,6), zeros(3)];
" E) _! n! \3 V2 X; Z( |) A! MQ = dec2mat(G,b,Q)
) n2 }) S) |/ ?  Ymu1 = dec2mat(G,b,mu1)
) s( Z. D6 w! }. kmu2 = dec2mat(G,b,mu2)
mu3 = dec2mat(G,b,mu3)
2 V& L2 f' v8 Y2 G2 Rs = dec2mat(G,b,s)
: N7 r8 P6 [% s

zaiyiaaaa

zzz.dan

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

IBB-EUT

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