[MATLAB工具箱] 请教高手MATLAB软件的LMI工具包中的mincx命令相关问题

[zhang_neu]

求解如下线性矩阵不等式，变量为α ，β ，Q，X

自己编写的代码如下：

A=[-0.0451,0.0361,0,-0.183;-0.37,-2.024,1,0;0.336,-6.94,-2.985,0;0,0,1,0];
B=[0,1.0;-0.16,0;-1.1,0;0,0];
D=eye(4);
E1=[-0.0045,0.0036,0,-0.0183;-0.0370,-0.2024,1,0;0.0336,-0.6940,-0.2985,0;0,0,1,0];
r=1.2;
q=1.3;

setlmis([]);

Q = lmivar(2, [2 4]);
X = lmivar(1, [4 1]);
α = lmivar(1, [1 0]);
β = lmivar(1, [1 0]);

lmiterm([1 1 1 0], D);
lmiterm([1 1 1 X], -r, 1);
lmiterm([1 1 2 X], A, 1);
lmiterm([1 1 2 Q], B, 1);
lmiterm([1 1 2 X], -q, 1);
lmiterm([1 2 2 X], -r, 1);
lmiterm([1 2 3 -X], 1, E1');
lmiterm([1 3 3 0], -1);

lmiterm([2 1 1 α], -1,1);
lmiterm([2 1 2 -Q], 1,1);
lmiterm([2 2 2 0], -1);

lmiterm([-3 1 1 X], 1,1);
lmiterm([-3 1 2 0], 1);
lmiterm([-3 2 2 β], 1,1);

lmis = getlmis;
c = mat2dec(lmis,eye(4));
options = [le-5,0,0,0,0];
[copt,xopt] = mincx(lmis,c,options)



Xopt=dec2mat(lmis, xopt, α+β)



因为手中材料有限，代码最后几句自己也不是很懂，希望得到高人指点，非常非常之感激！！！