zhang_neu
2009-05-26, 12:29
求解如下线性矩阵不等式,变量为α ,β ,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, α+β)
因为手中材料有限,代码最后几句自己也不是很懂,希望得到高人指点,非常非常之感激!!!
自己编写的代码如下:
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, α+β)
因为手中材料有限,代码最后几句自己也不是很懂,希望得到高人指点,非常非常之感激!!!