主题:[讨论]我是新手,不怎么会用,刚才发的好像不对,重新发了一次
自己编了一个解线性矩阵不等式的程序,但是运行完得到的矩阵L为零阵,拜托帮我看一下。谢谢!
E = [1 0;0 0];
A = [-2 0;0 1];
B = [0 0;0 1];
A1 = [0.3 0;0 0.7];
Q = [1 0;0 0];
R = 1;
% Initial a LMI system
setlmis([]);
% Define Variables
X = lmivar(2, [2 2]);
L = lmivar(2, [1 2]);
% pos in (1, 1)
lmiterm([1 1 1 X], A, 1,'s');
lmiterm([1 1 1 0], 1);
lmiterm([1 1 1 0], B);
% pos (1, 2)
lmiterm([1 1 2 X], 1, A1);
% pos(1, 3)
lmiterm([1 1 3 -L], 1, 1);
% pos(1, 4)
lmiterm([1 1 4 -L], 0.707, 1);
% pos(1, 5)
lmiterm([1 1 5 -L], 0.707, R);
% pos(1, 6)
lmiterm([1 1 6 X], 1, Q);
% pos(2, 2)
lmiterm([1 2 2 0], -0.5);
% pos(3, 3)
lmiterm([1 3 3 0], -1);
% pos(4, 4)
lmiterm([1 4 4 0], -1);
% pos(5, 5)
lmiterm([1 5 5 0], -1);
% pos(6, 6)
lmiterm([1 6 6 0],-1);
lmis=getlmis;
[tmin, feas]=feasp(lmis)