求助关于LMI的求解问题

ganea

求助一个LMI的解法问题，原问题在附件里，刚刚接触这东西，很菜，希望高手指点，不胜感激
a=[-10 10 0;28 -1 -30;0 30 -8/3];
c=[1 0 0];
d=[-5 5 0;56 -1 -30;0 30 -4/3];
e=[0.1;0;0];
k=[5 -5 0;28 0 0;0 0 4/3];
I=eye(3);

setlmis([])
p = lmivar(1,[3 1])
m = lmivar(2,[3 1])
w = lmivar(1,[1 1])
x = lmivar(1,[1 1])
lmiterm([-1 1 1 1],-1,a,'s')
lmiterm([-1 1 1 2],1,c,'s')
lmiterm([-1 2 2 0],1)
lmiterm([-1 2 1 0],c)
lmiterm([-2 1 1 3],1,1)
lmiterm([-2 2 1 0],e)
lmiterm([-2 2 2 inv(p)],1,1)
lmiterm([3 1 1 3],1,1)
lmiterm([-3 1 1 0],0.25)
lmiterm([-4 2 1 1],a,1)
lmiterm([-4 2 1 1],d,-1)
lmiterm([-4 2 1 1],k,1)
lmiterm([-4 1 1 x],1,I)
lmiterm([-4 2 2 0],I)
lmis = getlmis
c=mat2dec(lmis,0,0,0,eye(1))
options=[1e-5,0,0,0,0]
[copt,xopt]=mincx(lmis,c,options)
P=dec2mat(lmis,xopt,p)
W=dec2mat(lmis,xopt,w)
M=dec2mat(lmis,xopt,m)

[ 本帖最后由 eight 于 2008-5-8 22:08 编辑 ]