Home > other >  Matlab LMI problem for help
Matlab LMI problem for help

Time:09-23

About the solution of the theorem 1, below is my code, seek out no feasible solution, can you tell me what is going wrong?
CLC
The clear
A=[0, 1, 0, 0, 0, 1, 1-5-3].
Ad=0.6 * eye (3);
G=0.1 * eye (3);
B=[0.9 0.1 0.5 0.3 0, 0; - 0.5-0.2 to 0.1];
E=eye (3);
A1=[zeros (6, 3), [E; A]].
Ad1=[zeros (6, 3), [zeros (3), Ad]].
% Ad1=0.6 * eye (6);
G1=[zeros (3), G];
B1=[zeros (3, 3), B].
H1=[- eye (3), zeros (3)];
KN=[1.4443, 0.4449 to 0.3640-2.2578-0.6375-0.9807
0.4155 1.0127-0.1189-0.7128-0.2795-0.9251
- 0.5218-1.4436-1.1944-0.8794-0.0251 to 0.6915];
KN=0.5 * ones (3, 6);
Mu=0.4;
Z1=0.1 * eye (3);
Z1=[0.1 0 0 0 0; 0.2 0 0 0.3];
Z2=0.1 * eye (3);
=Z2 [0.3 0 0 0 0; 0.4 0 0 0.5];
Z=[Z1, zeros (3), zeros (3), Z2];
R=0.1 * eye (3);
% x1=[1 1 1 1 1 1) ';
Gamma=1;
% %
Setlmis ([])
Kappa=lmivar (1, 1, 1);
P=lmivar (2, 6, 6);
Q=lmivar (2, 6, 6);

Lmiterm (1 P] [- 2 1, 1, 1);
Lmiterm ([1 1-3 Q], 1, 1);
Lmiterm (1 1 kappa] [- 4, 1, 1);

Lmiterm (1 P] [1, 1, 1, A1 + B1 * KN);
Lmiterm (1 - P] [1, 1, + B1 * KN (A1) ", 1);
Lmiterm (1 1 Q] [1, 1, 1);
Lmiterm ([1 1 1 0], Z + KN '* R * KN);
Lmiterm (1 1 kappa] [1, 1, 1);
Lmiterm (1 2 P] [1, 1, Ad1);
Lmiterm (1 3 P] [1, 1, H1);
Lmiterm (2 Q] [1, 2, - (1 - mu), 1);
Lmiterm ([1 2 3 0], 0).
Lmiterm ([1 3 0] 3, - gamma);
%
Lmis=getlmis;
[tmin, xfea]=feasp (lmis).
P=dec2mat (lmis, xfea, P)
Q=dec2mat (lmis, xfea, Q)
Kappa=dec2mat (lmis, xfea, kappa)
  • Related