matlab 问题

nkkopd

Error using lmiterm (line 281)
2 g5 u& p0 X" Y% r: _lhs of LMI #1, block (1,1): term dimensions incompatible with other
terms in row #1 or column #1
/ X, W0 _/ J* K0 j) t4 n" b
6 d# e) h' i+ f  R( m8 S& e8 x5 FError in SCP_03 (line 103)
3 A$ {2 z0 P( I) Z; F- ylmiterm([1 1 1 P1],-1,1)
' K* {1 O9 Z, d9 ]
- }* v, r. R0 Bclear
& w) N% a& [, |3 ~4 Nclc

, l" L% o* C3 I% D- J" a+ hgamma = 0.95;

$ y( L& c+ K% o! z$ i% v, Y' {' Hsigma1 = 2;
1 o& D9 b& L/ S$ j2 vsigma2 = 1.98;
sigma3 = 1.86;
' Y  I7 Y# D  `0 `& Usigma = sigma1+sigma2+sigma3;
2 M, m% D9 A% r
gamma1 = 0.89;
gamma2 = 0.88;
gamma3 = 0.85;
Gamma = [gamma1, gamma1, gamma1; gamma2, gamma2, gamma2; gamma3, gamma3, gamma3];%3*3

Lambda = diag([2, 2]);%2*2

/ A1 Y6 v# Y$ h7 E+ s; qW = [kron(Gamma, Lambda), zeros(6,3); zeros(3,6), zeros(3)];%9*9
, F; C1 A& D$ U; N, _/ T
beta1 = 0.88;
: B0 G! d5 |. {; O* Gbeta2 = 0.85;
% W* S: F9 m- V; _beta3 = 0.87;
7 }2 Z' Q/ W$ eXi = diag([beta1, beta2, beta3]);%3*3
barXi = [kron(Xi, eye(2)), zeros(6,3); zeros(3,6), zeros(3)];%9*9

7 ^- }1 ~# U, P3 V$ y; I* KU = 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
: b: n; I: f' f! I: y5 W/ ^+ J
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
, |; I: U: E1 q! TM = [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

6 N  E! F1 m5 oK1 = [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
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
& G8 ^8 X8 C% F
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

7 D# u- w: @. KE = [0.3;0.3;0.3];%3*1

D = diag([0.5, 0.7, 0.4]);%3*3

% x1 e7 q$ b0 R$ z1 PA = [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
# s2 o# G4 x6 z# R3 s/ t
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
$ P& r% ^. b) ^* b. q9 z& i
theta(1) = 2;
theta(2) = 3;
1 k4 ?/ j1 [9 C; w3 s* A+ Qtheta(3) = 1;
theta(4) = 3;

1 h3 ^& ~! W, m1 ~& Q  {6 i
: l' P- n: m0 l4 v- a+ {4 j% z
( N' B  |+ u$ @. @+ Y2 Lfor k = 1:3
    i = theta(k);
+ I9 U8 t1 ~: J: l    j = theta(k+1);
% |! R) z- m4 s; C   
6 v) n' r( G) W$ ?( g$ P, K  Phi = zeros(3);
  
0 f* J2 T# p3 _) P6 Xfor t = 1:3
! Z: i/ z  n8 e: w& L6 Q9 G: z/ h     
: J$ Y$ ]  R' @$ p    if (i-t) ~= 0;
          delta(t) = 0;
  {; m8 Q1 e# a6 @) x    elseif (i-t) == 0;
. ^3 W* X2 H- g5 i; h7 ^* h- v1 r          delta(t) = 1;
         Phi(t,t) = delta(t);
; T* ]0 s+ ^5 O, v    end   
9 C* O, l7 y; [/ u4 R  M  
end

  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
5 s& X1 |; [, Y5 X. }8 B6 V  z  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
: k! ]7 g: E) Eend
$ ?! P" a% u) e- ^3 s  % N=3  n=2  m=1 p=1 r=1 q=1  M=3
8 R" M; g- \2 R( Q: A9 h) g- r1 k  setlmis([])
  P11 = lmivar(1,[2 1]);
  ^. Z6 A% F6 X  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
  P12 = lmivar(1,[2 1]);
1 B; J, m4 H& ^5 K% {$ N  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
9 l) `8 }; E$ D; O  P13 = lmivar(1,[2 1]);
. M& t8 a  B  [" E2 W/ _. y  P23 = lmivar(2,[2 1]);
6 i; K; ^9 C+ U$ k  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
; J4 `; j: T* J+ N. l! z  Z1 = lmivar(2,[9 3]);
$ T: \3 I, M6 H8 M  Z2 = lmivar(2,[9 3]);
  Z3 = lmivar(2,[9 3]);
  Q = lmivar(1,[9 1]);
  mu1 = lmivar(1,[1,1]);
( i  {/ n, ~. @' A. u  mu2 = lmivar(1,[1,1]);
; i3 V; o* O9 k, ?; A7 [/ P$ n  mu3 = lmivar(1,[1,1]);
9 Q2 B4 _! m+ M  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)
lmiterm([1 1 1 mu1],-1,U+U')
lmiterm([1 1 1 mu2],-1,L+L')
$ ]3 ]* @9 z) ^( V0 g( blmiterm([1 1 1 s],1,N'*N)
lmiterm([1 1 1 0],H'*H)
' E- g: |. [+ U1 w6 Flmiterm([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)
7 F. @/ Q4 \6 k2 elmiterm([1 2 2 P2],barC'*K3',K3*barC)
% S. ~* A1 P8 E' X  K3 mlmiterm([1 2 2 P1],-1,1)
lmiterm([1 2 2 mu1],-1,U+U')
# c2 L2 ~. ^+ U& ?' }lmiterm([1 2 2 mu2],-1,L+L')
. u" Q4 [+ M0 {: Ilmiterm([1 2 6 mu1],1,U')
6 ?/ X& e" r1 Q% |- vlmiterm([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)
. p- D" ~8 l7 [7 N! k3 J; Nlmiterm([1 2 13 P2],barA'+W',1)
lmiterm([1 2 13 -Z2],-barC',1)
9 t# |! }8 G4 F' j2 z2 D( w  klmiterm([1 3 3 mu3],-1,L')
, Y) J1 ]0 o; S. ^; G/ Tlmiterm([1 3 9 mu3],1,L+L')
. i* K& P. c4 ]& y2 W! Z* S7 Klmiterm([1 4 4 mu3],-1,L+L')
lmiterm([1 4 10 mu3],1,L')
$ H" h) @/ O, b1 Q0 Ylmiterm([1 5 5 mu1],-2,eye(9))
2 s' x5 c$ \2 q& w+ \' rlmiterm([1 5 12 P2],1,1)
lmiterm([1 6 6 mu1],-2,eye(9))
+ ~) {9 J" P, N2 F- e5 _# ilmiterm([1 6 13 P2],1,1)
lmiterm([1 7 7 Q],1,1)
lmiterm([1 7 7 mu2],-2,eye(9))
0 F2 v' P' ~# I+ ^+ }% |$ {lmiterm([1 8 8 Q],1,1)
lmiterm([1 8 8 mu2],-2,eye(9))
( E8 R. x/ O( L; V$ p- q% p! flmiterm([1 9 9 barGamma],2,1)
lmiterm([1 9 9 Q],-1,1)
3 u* x* v0 l6 C: tlmiterm([1 9 9 mu3],-2,eye(9))
8 M# e0 E- e# N- ilmiterm([1 9 12 P2],barXi',1)
lmiterm([1 10 10 Q],-1,1)
1 x3 O- e: K$ k* q. Z8 n+ xlmiterm([1 10 10 mu3],-2,eye(9))
lmiterm([1 10 13 P2],barXi',1)
' ?' S# W$ D8 G- \+ L$ C* G$ [& g- J5 ^lmiterm([1 11 11 P2],barD'*K1',K1*barD)
' }) i3 m* G! g, N, Vlmiterm([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)
; U) {$ P! d4 ?lmiterm([1 11 13 P2],barB',1)
( \5 {" ^( o+ s: J( `lmiterm([1 11 13 -Z2],-barD',1)
4 Y4 B5 f6 g4 j+ blmiterm([1 12 14 P2],1,M)
lmiterm([1 13 14 P2],1,M)
3 F9 A, k" Q1 ]9 D
lmiterm([-1 1 1 P1],1,1)
lmiterm([-2 1 1 P2],1,1)
; A" ^0 f) n  Y+ e& v- clmiterm([-3 1 1 P3],1,1)

lmiterm([2 1 1 barLambda],2,1)
+ Z& r' ~+ g' K+ e+ p3 qlmiterm([2 1 1 P2],-1,1)
lmiterm([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)
lmiterm([2 2 2 P3],barC'*K1',K1*barC)
lmiterm([2 2 2 P3],barC'*K2',K2*barC)
: H7 U: Z2 v/ `6 r! H$ x5 t/ `lmiterm([2 2 2 P3],barC'*K3',K3*barC)
1 b9 z2 S4 d; {6 A7 n# M: Ylmiterm([2 2 2 P2],-1,1)
lmiterm([2 2 2 mu1],-1,U+U')
lmiterm([2 2 2 mu2],-1,L+L')
lmiterm([2 2 6 mu1],1,U')
; g' N1 a: t  n+ E9 tlmiterm([2 2 8 mu2],1,L')
' \6 y4 J4 R6 ?lmiterm([2 2 11 P3],barC'*K1',K1*barD)
lmiterm([2 2 11 P3],barC'*K2',K2*barD)
. g1 E: l+ Q& w6 Tlmiterm([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')
  Y/ v" ^# N4 w  e' dlmiterm([2 4 4 mu3],-1,L+L')
6 m1 j) c! l. ?, m3 zlmiterm([2 4 10 mu3],1,L')
1 I6 |1 h1 d) j& U3 almiterm([2 5 5 mu1],-2,eye(9))
$ ^: c; N. A9 z; @  wlmiterm([2 5 12 P3],1,1)
6 _8 F& i4 Z. N* rlmiterm([2 6 6 mu1],-2,eye(9))
lmiterm([2 6 13 P3],1,1)
lmiterm([2 7 7 Q],1,1)
  Y8 M  x1 l. m1 v5 ?6 [( G1 ^lmiterm([2 7 7 mu2],-2,eye(9))
lmiterm([2 8 8 Q],1,1)
0 v: a8 i: c  Vlmiterm([2 8 8 mu2],-2,eye(9))
lmiterm([2 9 9 barGamma],2,1)
5 I# N6 L6 v& J$ ]+ tlmiterm([2 9 9 Q],-1,1)
lmiterm([2 9 9 mu3],-2,eye(9))
4 v# ?: r# \) C9 L. A4 Rlmiterm([2 9 12 P3],barXi',1)
# u, `0 R+ ?& W$ xlmiterm([2 10 10 Q],-1,1)
lmiterm([2 10 10 mu3],-2,eye(9))
7 c1 A' [" Z+ ]/ }7 \lmiterm([2 10 13 P3],barXi',1)
- T+ M* ?5 z/ U0 `8 slmiterm([2 11 11 P3],barD'*K1',K1*barD)
& D4 S2 g& o. M* X& @3 d1 C1 Jlmiterm([2 11 11 P3],barD'*K2',K2*barD)
* ^5 p% [' ~/ v5 B6 }: ]lmiterm([2 11 11 P3],barD'*K3',K3*barD)
lmiterm([2 11 11 0],-(gamma^2)*eye(4))
lmiterm([2 11 12 P3],barB',1)
4 A% Q9 i/ e! M+ Mlmiterm([2 11 13 P3],barB',1)
lmiterm([2 11 13 -Z3],-barD',1)
lmiterm([2 12 14 P3],1,M)
* R* L( h' }  F  Y. n+ W7 almiterm([2 13 14 P3],1,M)
: G. \, R5 k- Q- G- M3 U) A
G=getlmis;
$ p% i( B3 D* ~7 n; I& y[tmin b]=feasp(G);
* G* H; ~* ~! j
P11 = dec2mat(G,b,P11);
P21 = dec2mat(G,b,P21);
P31 = dec2mat(G,b,P31);
8 I3 O( |$ \  h, {$ ]+ iP1 = [kron(eye(3), P11), kron(eye(3), P21); kron(eye(3), P21'), kron(eye(3), P31)]
- n9 H8 ^/ s( t  SP12 = dec2mat(G,b,P12);
P22 = dec2mat(G,b,P22);
* I7 G! O& z+ D9 T1 {6 d' N$ ZP32 = dec2mat(G,b,P32);
9 Y! `. k+ U8 U1 [% l# F' M5 [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);
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));
$ F) R% r- @8 ~$ A- B* i  k' y8 Z* TZ(3) = dec2mat(G,b,Z(3));
1 X3 W$ v5 v( ^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)
mu3 = dec2mat(G,b,mu3)
s = dec2mat(G,b,s)
5 e; x4 V7 Z. P! s/ m$ z

zaiyiaaaa

zzz.dan

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

IBB-EUT

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