matlab 问题

nkkopd

Error using lmiterm (line 281)
lhs of LMI #1, block (1,1): term dimensions incompatible with other
( a( `7 `, J- L. P  Q( }6 c( nterms in row #1 or column #1
7 H! m4 |* q( o* `) y, A' B0 l
Error in SCP_03 (line 103)
2 ^) G* H: l' _+ q0 c0 y5 A3 L' o4 vlmiterm([1 1 1 P1],-1,1)
+ u9 m: f1 J6 w- M% ?6 Q
" }+ |3 I7 Q0 T% z) Y& Q" t* m8 bclear
0 k9 z& ?1 `8 f' W: C! ?clc
. @! ]* y$ P8 H: ~
( G; M4 `2 K$ v- _1 R2 z& z* {gamma = 0.95;
5 F4 q! h9 Q  {: I$ h9 [2 Z
' i* a+ }3 S, V) Z5 q, C5 n+ D7 _2 wsigma1 = 2;
sigma2 = 1.98;
" l5 l, ~" a7 S4 `- G7 N4 j% Ysigma3 = 1.86;
8 {, y0 M! S5 o2 l: Bsigma = sigma1+sigma2+sigma3;
# ^  x0 n" @/ g3 W/ O+ G/ G
! R7 G$ g* K! [+ j; Vgamma1 = 0.89;
% S8 L- |2 h1 T: @: ygamma2 = 0.88;
% K+ z1 r$ s0 u( l% L# o+ `8 Lgamma3 = 0.85;
Gamma = [gamma1, gamma1, gamma1; gamma2, gamma2, gamma2; gamma3, gamma3, gamma3];%3*3
) H) c  A! X0 o
Lambda = diag([2, 2]);%2*2
1 t. I4 }& X) W- O
( D$ V) k$ E$ I0 P+ {  u$ SW = [kron(Gamma, Lambda), zeros(6,3); zeros(3,6), zeros(3)];%9*9
9 j9 ^8 a# T; ~  ~% y, f3 b
beta1 = 0.88;
beta2 = 0.85;
3 \2 Z; v% ?" i/ i! `% mbeta3 = 0.87;
. U6 a; n, _8 v5 g2 R0 oXi = diag([beta1, beta2, beta3]);%3*3
, G6 C% M$ _3 m# p/ }( D7 HbarXi = [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
, h+ s9 H+ H9 V$ H! Q% r7 M# [L = diag([0.34, 0.34, 0.34, 0.34, 0.34, 0.34, 0, 0, 0]);%花L1+L2  9*9
7 {  H3 c1 c4 v, ?; U% _( Z
! d0 a' F7 |' G& B) {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
# D7 R4 P/ l/ z' b$ S0 i6 h9 UM = [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
+ j; r4 z# [+ w
) y2 e' \. \( q1 n% QC =[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
% L+ V3 s( ]& R/ m3 _
5 X7 O+ `8 R) n4 `( i# ^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
1 X6 G5 C$ }, q  F9 E; pK2 = [-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
* y2 ?/ B/ C6 {" X, {# wK3 =[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
6 I! X; m( f7 d; c1 o1 f2 Q; m
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
1 e, ]  `/ ]0 C8 I7 J0 p" [! g
" f9 W; i- _. R1 RE = [0.3;0.3;0.3];%3*1

, m0 u" b7 ?, U) M" k+ z0 g' nD = diag([0.5, 0.7, 0.4]);%3*3

/ o1 L: Q2 m0 ?% h( aA = [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
0 t2 D1 N# T! U9 e! y- W" Z3 J
' y; w4 @, c4 LB = [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
1 [! N! o% o: v) p  F
theta(1) = 2;
theta(2) = 3;
theta(3) = 1;
- L! f9 h5 V) v  }: @! rtheta(4) = 3;
2 q+ L* Q% v9 q1 I7 ]$ }( v: O

, y' e. o1 y4 e3 _3 c
3 s7 L* w$ m% g% {* Qfor k = 1:3
( n/ Z7 d9 b, a0 b. d! P; B0 m4 s    i = theta(k);
3 u: b9 Z* R) L' \' Q    j = theta(k+1);
   
  Phi = zeros(3);
  
for t = 1:3
6 l3 J) C# i& p. F+ n+ E5 Q1 \     
0 w) `5 {' ~2 {, V7 r+ q" C    if (i-t) ~= 0;
9 S6 E; C) e+ \( @$ ?( S* }0 |' F; G          delta(t) = 0;
9 C  m( u+ r6 v    elseif (i-t) == 0;
          delta(t) = 1;
         Phi(t,t) = delta(t);
; V- B: {' j7 l" m2 l: m% ?    end   
" }! T' |* ^1 [2 R! D  
8 B: K0 J+ ]) Q8 r- L- @* ^  d0 ^% Fend
# B! A! d6 E3 _9 q/ {" L! Q
. p' F' i( b" c! X" T6 A  barA = [A, zeros(6,3); (kron(Phi,eye(1)))*C, kron((eye(3)-Phi),eye(1))];%9*9
  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
6 s' L& ]5 A9 M) A( c) e& o  barD = [(kron(Phi,eye(1)))*E, kron(Phi,eye(1))*D];%3*4
end
1 t. z( D" F1 z  % N=3  n=2  m=1 p=1 r=1 q=1  M=3
  setlmis([])
. b4 o2 J' [8 U  j# h  P11 = lmivar(1,[2 1]);
  P21 = lmivar(2,[2 1]);
$ Q$ K' u: l' g, y% q5 v  P31 = lmivar(1,[1 1]);
( ^1 c2 N) v7 `9 y; O5 A3 ]  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]);
# A4 K( U% a: ^% d) Z& p/ n  P32 = lmivar(1,[1 1]);
! A2 ?! S0 {- g5 X7 \6 s5 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]);
$ \7 Y$ p1 H5 E/ G; e+ T  P33 = lmivar(1,[1 1]);
8 w* m1 z+ V3 O1 W  P3 = lmivar(3,[kron(eye(3),P13),kron(eye(3),P23);kron(eye(3),P23'),kron(eye(3),P33)]);%9*9
" f& d+ R* {9 J% M5 ^, F6 X  Z1 = lmivar(2,[9 3]);
/ f- Q7 \- @2 C5 O7 h; v' f  Z2 = lmivar(2,[9 3]);
+ t- f4 T* V$ [; |! `$ N$ V1 }  Z3 = lmivar(2,[9 3]);
  Q = lmivar(1,[9 1]);
" B) r( g, }# [& M) I. F2 B& F  mu1 = lmivar(1,[1,1]);
9 a8 L/ v$ g/ W4 |2 R; x  mu2 = lmivar(1,[1,1]);
: _4 P' L7 t, M8 x: S  mu3 = lmivar(1,[1,1]);
2 b/ M: Q% \' Z0 P% N. o  s = lmivar(1,[1,1]);
  barGamma = lmivar(1,[9 1]);
  barLambda = lmivar(1,[9 1]);
/ U8 |' U5 H: s$ N& [" q6 u' g& D  
lmiterm([1 1 1 barLambda],2,1)
lmiterm([1 1 1 P1],-1,1)
lmiterm([1 1 1 mu1],-1,U+U')
) C$ E* u' z$ y! W  @/ `lmiterm([1 1 1 mu2],-1,L+L')
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)
5 K& x4 O- D0 q$ glmiterm([1 2 2 P2],barC'*K2',K2*barC)
2 M+ ^) j! S2 w3 I! M: v! h1 almiterm([1 2 2 P2],barC'*K3',K3*barC)
lmiterm([1 2 2 P1],-1,1)
9 N7 K% g  Y1 t9 Y/ Vlmiterm([1 2 2 mu1],-1,U+U')
lmiterm([1 2 2 mu2],-1,L+L')
lmiterm([1 2 6 mu1],1,U')
1 S( n( S# K/ \( G2 A" Ilmiterm([1 2 8 mu2],1,L')
lmiterm([1 2 11 P2],barC'*K1',K1*barD)
4 z8 |# Z/ t: E2 |/ M; ~+ qlmiterm([1 2 11 P2],barC'*K2',K2*barD)
6 ?0 {9 F4 h/ V, K3 \) almiterm([1 2 11 P2],barC'*K3',K3*barD)
lmiterm([1 2 13 P2],barA'+W',1)
( \" b3 P1 H' }' t2 s1 \lmiterm([1 2 13 -Z2],-barC',1)
1 z" L1 b: n! A& B- J6 H  Z; plmiterm([1 3 3 mu3],-1,L')
' O9 s- d8 f0 Elmiterm([1 3 9 mu3],1,L+L')
lmiterm([1 4 4 mu3],-1,L+L')
lmiterm([1 4 10 mu3],1,L')
* P* F* L- g; D5 [* |' x. Zlmiterm([1 5 5 mu1],-2,eye(9))
lmiterm([1 5 12 P2],1,1)
* @. ?4 a. K# j5 |6 Almiterm([1 6 6 mu1],-2,eye(9))
lmiterm([1 6 13 P2],1,1)
1 ]5 \9 i( v& Qlmiterm([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))
- k* D. X# K" T1 @/ N9 x/ Qlmiterm([1 9 9 barGamma],2,1)
lmiterm([1 9 9 Q],-1,1)
lmiterm([1 9 9 mu3],-2,eye(9))
" P1 W6 z7 {  i5 @lmiterm([1 9 12 P2],barXi',1)
lmiterm([1 10 10 Q],-1,1)
1 }# E( k7 y. _4 ~lmiterm([1 10 10 mu3],-2,eye(9))
0 o0 f" N/ e/ K9 Elmiterm([1 10 13 P2],barXi',1)
7 W. p5 ^7 R" ]0 f6 D  a% K% ~' r- glmiterm([1 11 11 P2],barD'*K1',K1*barD)
, x9 ~  {' e% y* z! ylmiterm([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)
lmiterm([1 11 13 -Z2],-barD',1)
, ?, c+ h+ e' M* g+ \4 b  elmiterm([1 12 14 P2],1,M)
lmiterm([1 13 14 P2],1,M)
) J- z' A, n: x- K
lmiterm([-1 1 1 P1],1,1)
lmiterm([-2 1 1 P2],1,1)
lmiterm([-3 1 1 P3],1,1)

lmiterm([2 1 1 barLambda],2,1)
lmiterm([2 1 1 P2],-1,1)
8 l6 V4 ?# X8 t3 zlmiterm([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')
# v9 I  ~6 W* p  U3 |1 x# Tlmiterm([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)
+ W6 V% i6 E& Y- h1 T9 k# Y9 Dlmiterm([2 2 2 P3],barC'*K3',K3*barC)
5 d0 G- Z) b& w) ?lmiterm([2 2 2 P2],-1,1)
9 d: |. i2 Z( g- x3 r/ mlmiterm([2 2 2 mu1],-1,U+U')
lmiterm([2 2 2 mu2],-1,L+L')
8 J( M7 ?7 F& q; @7 Klmiterm([2 2 6 mu1],1,U')
lmiterm([2 2 8 mu2],1,L')
lmiterm([2 2 11 P3],barC'*K1',K1*barD)
. j. W7 m1 X# R0 v5 Q; f1 qlmiterm([2 2 11 P3],barC'*K2',K2*barD)
lmiterm([2 2 11 P3],barC'*K3',K3*barD)
& Q$ |. G; U* ^9 O/ _lmiterm([2 2 13 P3],barA'+W',1)
lmiterm([2 2 13 -Z3],-barC',1)
, K4 b8 u9 U! plmiterm([2 3 3 mu3],-1,L')
lmiterm([2 3 9 mu3],1,L+L')
lmiterm([2 4 4 mu3],-1,L+L')
lmiterm([2 4 10 mu3],1,L')
4 n2 Q( s! Z" M  X) y8 elmiterm([2 5 5 mu1],-2,eye(9))
lmiterm([2 5 12 P3],1,1)
4 W/ s0 [& L/ R# [lmiterm([2 6 6 mu1],-2,eye(9))
. E5 ^4 {, E# F9 d$ x0 U+ i# o# m3 Qlmiterm([2 6 13 P3],1,1)
lmiterm([2 7 7 Q],1,1)
lmiterm([2 7 7 mu2],-2,eye(9))
+ [0 V' M" d7 [6 xlmiterm([2 8 8 Q],1,1)
) g5 i9 M' N  q  y  Ylmiterm([2 8 8 mu2],-2,eye(9))
lmiterm([2 9 9 barGamma],2,1)
lmiterm([2 9 9 Q],-1,1)
% x/ B  Z8 w2 llmiterm([2 9 9 mu3],-2,eye(9))
/ ^! R. @+ k9 M  H" slmiterm([2 9 12 P3],barXi',1)
& k$ o; d) j* e: L# ]0 _lmiterm([2 10 10 Q],-1,1)
( z2 e7 ?/ p6 N: v# p) ylmiterm([2 10 10 mu3],-2,eye(9))
9 L5 W5 A' Q) @# l1 Qlmiterm([2 10 13 P3],barXi',1)
lmiterm([2 11 11 P3],barD'*K1',K1*barD)
2 {5 l8 a+ c5 g9 i5 u  B# r) almiterm([2 11 11 P3],barD'*K2',K2*barD)
lmiterm([2 11 11 P3],barD'*K3',K3*barD)
( F2 @" n, N4 L* i- z0 Qlmiterm([2 11 11 0],-(gamma^2)*eye(4))
lmiterm([2 11 12 P3],barB',1)
lmiterm([2 11 13 P3],barB',1)
lmiterm([2 11 13 -Z3],-barD',1)
' h+ A4 }3 s3 l4 e! \lmiterm([2 12 14 P3],1,M)
lmiterm([2 13 14 P3],1,M)
( G5 B9 h  o/ {! v: \: p
9 ^; m3 g$ A/ y" h5 zG=getlmis;
: Z. l0 ]9 k: \8 p# ^[tmin b]=feasp(G);

P11 = dec2mat(G,b,P11);
P21 = dec2mat(G,b,P21);
, d' ?9 T: P/ M- y8 kP31 = dec2mat(G,b,P31);
3 D( R. h7 ]  }( |. \P1 = [kron(eye(3), P11), kron(eye(3), P21); kron(eye(3), P21'), kron(eye(3), P31)]
& \3 ^; {* R+ SP12 = dec2mat(G,b,P12);
P22 = dec2mat(G,b,P22);
2 |3 v+ u/ K6 L3 H! SP32 = dec2mat(G,b,P32);
; X; J$ G; U9 JP2 = [kron(eye(3), P12), kron(eye(3), P22); kron(eye(3), P22'), kron(eye(3), P32)]
, v7 p) h6 s4 f. V) M' eP13 = dec2mat(G,b,P13);
P23 = dec2mat(G,b,P23);
, P. q6 l7 L0 Q9 `P33 = dec2mat(G,b,P33);
9 v; g" s: N9 K9 K" A0 I, C6 d6 z  YP3 = [kron(eye(3), P13), kron(eye(3), P23); kron(eye(3), P23'), kron(eye(3), P33)]
Z(1) = dec2mat(G,b,Z(1));
1 {: q. |  f3 R7 L. BZ(2) = dec2mat(G,b,Z(2));
$ f9 J# {: B0 [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)
mu1 = dec2mat(G,b,mu1)
mu2 = dec2mat(G,b,mu2)
/ h$ @$ Q! O* p3 fmu3 = dec2mat(G,b,mu3)
s = dec2mat(G,b,s)
* ?8 X5 c! i! U# c. i3 }

zaiyiaaaa

zzz.dan

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

IBB-EUT

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