matlab 问题

nkkopd

Error using lmiterm (line 281)
& J3 l0 `6 O' |$ x3 m! {# Q# Ilhs of LMI #1, block (1,1): term dimensions incompatible with other
" k7 p* x9 W1 Q$ T' l2 {7 oterms in row #1 or column #1
& o$ m# l2 {- M3 B
5 C5 v# {( i1 r% ?. O0 P, Z  hError in SCP_03 (line 103)
' r+ o( k1 `8 k0 p! P$ plmiterm([1 1 1 P1],-1,1)

! l! w3 L: ?+ r% [0 oclear
clc
; `/ [# D/ D& A% K. A" @; _
: i: U" b: G% Z+ H, ~# B/ _' lgamma = 0.95;

0 _, H* r# v# R+ Y" T5 b& Tsigma1 = 2;
2 v8 d5 e9 a2 Z! ]sigma2 = 1.98;
sigma3 = 1.86;
9 F5 H  i" {2 g* ?" b; Msigma = sigma1+sigma2+sigma3;
0 k, B) M0 c: h. y+ e% m
+ p9 c5 V4 z9 ?( C4 Y$ Y! Egamma1 = 0.89;
) b* @6 R' t+ j" F* Hgamma2 = 0.88;
gamma3 = 0.85;
1 x$ h/ D2 d3 L  ZGamma = [gamma1, gamma1, gamma1; gamma2, gamma2, gamma2; gamma3, gamma3, gamma3];%3*3

Lambda = diag([2, 2]);%2*2
! W& Z3 k# a6 M, T% c, F: D
W = [kron(Gamma, Lambda), zeros(6,3); zeros(3,6), zeros(3)];%9*9

beta1 = 0.88;
3 i+ o. ?4 ?. `  t. d6 K/ Cbeta2 = 0.85;
# m4 {/ u6 I& l' U& xbeta3 = 0.87;
Xi = diag([beta1, beta2, beta3]);%3*3
+ E6 B6 e  `; |barXi = [kron(Xi, eye(2)), zeros(6,3); zeros(3,6), zeros(3)];%9*9
; Y  F* s4 g2 p( J; J: |4 ~
U = diag([1.68, 1.5, 1.68, 1.5, 1.68, 1.5, 0, 0, 0]);%花U1+U2   9*9
5 \1 C8 q; x5 B" f" @0 g$ iL = diag([0.34, 0.34, 0.34, 0.34, 0.34, 0.34, 0, 0, 0]);%花L1+L2  9*9
/ ^9 G2 g8 Z2 @0 c
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
; ~+ f2 U! w  A# r! LM = [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
: v$ W$ I5 Y, j* D
# O# f% O% S& i! X5 G4 ^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
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
' t5 d% S$ j0 m( T' D
' v6 d* ^) w# z; O" }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

E = [0.3;0.3;0.3];%3*1
0 N0 V5 d0 q: l' x8 d6 o# C
1 @  U; ?# n, PD = diag([0.5, 0.7, 0.4]);%3*3
9 E# [. o+ y7 ?5 Z
5 b; u! q8 I' _+ d* |; C# KA = [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
4 S1 c, A" O* H: f+ `8 \
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
5 k6 {; U7 ]' Z3 O
$ I2 o$ `, E* `  ~: Ntheta(1) = 2;
theta(2) = 3;
  m1 e: l8 W2 r, K' b: btheta(3) = 1;
theta(4) = 3;
( {. s5 D8 }$ T

4 N. p" L; E8 k' [" P# y! n
for k = 1:3
    i = theta(k);
    j = theta(k+1);
   
  Phi = zeros(3);
* P8 D8 x3 s1 c" \7 J6 x( m% `8 F  
# C+ Q3 V, z% G1 {: `: ofor t = 1:3
) t; a- T  {8 E     
    if (i-t) ~= 0;
          delta(t) = 0;
- \! ]& M- {( I7 w7 z+ Q    elseif (i-t) == 0;
1 b! J0 C; a$ d2 z% L          delta(t) = 1;
         Phi(t,t) = delta(t);
9 D" w( M' D4 Y. U' ?1 t4 _, J    end   
  
end

( @( r1 `. m( q/ i2 p; t1 n  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
  barD = [(kron(Phi,eye(1)))*E, kron(Phi,eye(1))*D];%3*4
end
  % N=3  n=2  m=1 p=1 r=1 q=1  M=3
  setlmis([])
7 }. b0 e3 D4 z+ ^  P11 = lmivar(1,[2 1]);
- C; `; f" o: @8 `0 T0 j0 C* @  P21 = lmivar(2,[2 1]);
2 w7 x. }9 N( P/ p& j- T  P31 = lmivar(1,[1 1]);
6 ~! k$ Q5 W* _* O  P1 = lmivar(3,[kron(eye(3),P11),kron(eye(3),P21);kron(eye(3),P21'),kron(eye(3),P31)]);%9*9
# X: i- j) k) w/ e' _4 k% S" B  P12 = lmivar(1,[2 1]);
  P22 = lmivar(2,[2 1]);
( Q4 |5 Z/ ~1 E) y" y9 g  P32 = lmivar(1,[1 1]);
! r* P( |1 ?/ j/ H' B4 i: K/ e7 Y  P2 = lmivar(3,[kron(eye(3),P12),kron(eye(3),P22);kron(eye(3),P22'),kron(eye(3),P32)]);%9*9
5 {2 M* v" q0 K6 h, K9 X& H4 B  P13 = lmivar(1,[2 1]);
  P23 = lmivar(2,[2 1]);
  P33 = lmivar(1,[1 1]);
4 J1 O6 B; F; Q. B) J5 g  P3 = lmivar(3,[kron(eye(3),P13),kron(eye(3),P23);kron(eye(3),P23'),kron(eye(3),P33)]);%9*9
( O* Q7 ~: u) O  Z1 = lmivar(2,[9 3]);
# k- }0 p2 G1 P& q( n  L! y  Z2 = lmivar(2,[9 3]);
  Z3 = lmivar(2,[9 3]);
, j3 a! T0 W- p* y5 W$ x  Q = lmivar(1,[9 1]);
( F; ]$ q* A6 h. s, I  mu1 = lmivar(1,[1,1]);
2 c/ ]; Y* W4 m5 t& t  mu2 = lmivar(1,[1,1]);
' w+ I& r0 e9 L: F: }  mu3 = lmivar(1,[1,1]);
6 L2 [. X. h. Z" i$ T) X  s = lmivar(1,[1,1]);
  barGamma = lmivar(1,[9 1]);
' R9 j: a3 K2 g4 M  E# [  barLambda = lmivar(1,[9 1]);
0 L2 T& G% y7 G4 e- |" m/ Z* V  
' e9 m, H4 D; olmiterm([1 1 1 barLambda],2,1)
lmiterm([1 1 1 P1],-1,1)
& i) h0 P# w4 a& _+ g6 flmiterm([1 1 1 mu1],-1,U+U')
lmiterm([1 1 1 mu2],-1,L+L')
: t9 s+ `  [* h$ f8 S- d. Z" o" Olmiterm([1 1 1 s],1,N'*N)
/ B+ P) a7 j$ K) n7 O" Tlmiterm([1 1 1 0],H'*H)
5 _* Y6 }2 M/ e6 d) B, ?lmiterm([1 1 5 mu1],1,U')
. w# c$ a7 B. \* Qlmiterm([1 1 7 mu2],1,L')
5 p9 w6 K9 k0 hlmiterm([1 1 12 P2],barA'+W',1)
lmiterm([1 2 2 P2],barC'*K1',K1*barC)
% s7 n$ R. H' elmiterm([1 2 2 P2],barC'*K2',K2*barC)
5 I. n+ n$ ?# \9 m. r9 e8 O7 h5 vlmiterm([1 2 2 P2],barC'*K3',K3*barC)
" f! N( L" {7 B7 Blmiterm([1 2 2 P1],-1,1)
lmiterm([1 2 2 mu1],-1,U+U')
, O5 N' V& q- n" M" tlmiterm([1 2 2 mu2],-1,L+L')
8 f: G; a# \+ blmiterm([1 2 6 mu1],1,U')
lmiterm([1 2 8 mu2],1,L')
lmiterm([1 2 11 P2],barC'*K1',K1*barD)
9 b0 `" ?0 R0 Z. ~' Ilmiterm([1 2 11 P2],barC'*K2',K2*barD)
lmiterm([1 2 11 P2],barC'*K3',K3*barD)
# ~5 P- s' b( clmiterm([1 2 13 P2],barA'+W',1)
9 R$ P1 F1 o2 v' \6 r0 m; c: dlmiterm([1 2 13 -Z2],-barC',1)
lmiterm([1 3 3 mu3],-1,L')
8 |' L- b0 s; {1 Z" X* t4 olmiterm([1 3 9 mu3],1,L+L')
* S5 O  S: |, glmiterm([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)
lmiterm([1 6 6 mu1],-2,eye(9))
+ n* \/ G+ k( e+ O/ S) L0 Z' Wlmiterm([1 6 13 P2],1,1)
  q2 B; f* o  c+ w. a2 ulmiterm([1 7 7 Q],1,1)
lmiterm([1 7 7 mu2],-2,eye(9))
, L8 Z% b' d, N. o3 A4 r' plmiterm([1 8 8 Q],1,1)
0 b' B, O) ^9 o' Rlmiterm([1 8 8 mu2],-2,eye(9))
lmiterm([1 9 9 barGamma],2,1)
" d" o% j7 T2 j+ j4 q) Xlmiterm([1 9 9 Q],-1,1)
lmiterm([1 9 9 mu3],-2,eye(9))
1 P8 P( d% E- ^; A" nlmiterm([1 9 12 P2],barXi',1)
lmiterm([1 10 10 Q],-1,1)
lmiterm([1 10 10 mu3],-2,eye(9))
lmiterm([1 10 13 P2],barXi',1)
lmiterm([1 11 11 P2],barD'*K1',K1*barD)
8 }/ t! f! O* Q& w3 v  blmiterm([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)
2 J4 h( A% e% B7 Z1 N4 Rlmiterm([1 11 13 -Z2],-barD',1)
lmiterm([1 12 14 P2],1,M)
& u  u: c! K- [8 }7 ~2 K1 Slmiterm([1 13 14 P2],1,M)

1 b5 H& |+ k7 ?. dlmiterm([-1 1 1 P1],1,1)
' K8 ]. p( g  m, F  o6 u: Almiterm([-2 1 1 P2],1,1)
lmiterm([-3 1 1 P3],1,1)

# O% J$ C, ?, t0 Jlmiterm([2 1 1 barLambda],2,1)
! p9 J1 M6 d( e7 {+ vlmiterm([2 1 1 P2],-1,1)
8 x8 O: S/ k  w( clmiterm([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)
# o$ d( q0 ]3 O$ d2 S) Rlmiterm([2 1 5 mu1],1,U')
+ g, h& m; j8 x5 \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)
1 z8 n+ W% x0 z1 u# g9 @1 olmiterm([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')
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)
5 ?- q) T. L: nlmiterm([2 2 11 P3],barC'*K3',K3*barD)
. M4 o* g2 e' J  c- F: o3 Ulmiterm([2 2 13 P3],barA'+W',1)
2 P* E$ T' m# b: D' dlmiterm([2 2 13 -Z3],-barC',1)
9 T' o' d7 Q7 H* plmiterm([2 3 3 mu3],-1,L')
) t& A5 ~1 q4 B0 i6 H2 r( j  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)
lmiterm([2 6 6 mu1],-2,eye(9))
- o$ x% `1 X+ I2 d& D% \( T* p/ G1 ]lmiterm([2 6 13 P3],1,1)
9 T* I8 `) }! X" y6 elmiterm([2 7 7 Q],1,1)
5 w& E7 h+ r2 W7 C' J7 c4 Q' b3 p. @lmiterm([2 7 7 mu2],-2,eye(9))
4 }$ d* N! u- }4 z9 D: T9 ~lmiterm([2 8 8 Q],1,1)
9 Z% G: v' m: h6 v0 m7 J$ }lmiterm([2 8 8 mu2],-2,eye(9))
lmiterm([2 9 9 barGamma],2,1)
" R% |# V# ~# {- m( Zlmiterm([2 9 9 Q],-1,1)
lmiterm([2 9 9 mu3],-2,eye(9))
5 m6 y! {: [" T5 @0 F" Plmiterm([2 9 12 P3],barXi',1)
lmiterm([2 10 10 Q],-1,1)
" t1 Y/ y( ^4 s5 P  Nlmiterm([2 10 10 mu3],-2,eye(9))
, {9 ]% w. l: }5 m7 e7 Clmiterm([2 10 13 P3],barXi',1)
lmiterm([2 11 11 P3],barD'*K1',K1*barD)
lmiterm([2 11 11 P3],barD'*K2',K2*barD)
lmiterm([2 11 11 P3],barD'*K3',K3*barD)
lmiterm([2 11 11 0],-(gamma^2)*eye(4))
0 d: u6 K4 A8 }% z6 o  O& Almiterm([2 11 12 P3],barB',1)
lmiterm([2 11 13 P3],barB',1)
* g( u+ t, S8 |. W) `8 l9 H" plmiterm([2 11 13 -Z3],-barD',1)
8 t+ H8 ^' O" g. Z6 Wlmiterm([2 12 14 P3],1,M)
7 j+ h4 n* v" ^" j% x1 ]  z8 W* Rlmiterm([2 13 14 P3],1,M)
6 y  g( G& q% y& `* H
6 q7 ]% Q& D6 YG=getlmis;
[tmin b]=feasp(G);
: I. X  {+ F# ^8 u2 A# \0 j
P11 = dec2mat(G,b,P11);
. o4 U9 l0 U4 [/ o4 \4 [P21 = dec2mat(G,b,P21);
P31 = dec2mat(G,b,P31);
% n; K) _1 N2 b$ r8 }0 FP1 = [kron(eye(3), P11), kron(eye(3), P21); kron(eye(3), P21'), kron(eye(3), P31)]
, Y% H% A6 N, [8 GP12 = dec2mat(G,b,P12);
P22 = dec2mat(G,b,P22);
$ W5 y* m$ W* o5 ~$ j! cP32 = dec2mat(G,b,P32);
P2 = [kron(eye(3), P12), kron(eye(3), P22); kron(eye(3), P22'), kron(eye(3), P32)]
P13 = dec2mat(G,b,P13);
; E) X0 _0 v: q. c7 _2 u; SP23 = 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)]
2 ^: [6 @7 i% v! PZ(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)];
5 A6 F8 B  a0 R9 o  b; l+ MbarLambda = [kron(sigma*eye(3), Lambda'*(P(1i).P(i))*Lambda), zeros(6,3); zeros(3,6), zeros(3)];
Q = dec2mat(G,b,Q)
9 M7 q: w( R% _6 v* n' _. ymu1 = dec2mat(G,b,mu1)
mu2 = dec2mat(G,b,mu2)
mu3 = dec2mat(G,b,mu3)
s = dec2mat(G,b,s)

IBB-EUT

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

zzz.dan

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

zaiyiaaaa