请问关于LMI工具箱中直接求矩阵>0的算法

superstones

本人对编程很水，最近碰到一个关于LMI的问题，希望各位高手能帮助一下，多谢！！
问题描述如下：
已知一个对称、正定矩阵P>0,(比如P=(x,1,2; 1,y,3; 2,3,z)>0    )
问如何利用LMI方法,在Matlab中编程实现，求出变量x,y,z？
查看了一些关于LMI的论文，发现都没有直接求已知矩阵>0或者<0的
希望各位高手能不吝指教！
   谢谢！

[ 本帖最后由 eight 于 2007-4-17 23:14 编辑 ]

superstones

有人知道吗？谢谢啊

ff630

a=[0 0 0 0;0 -0.088317 0.625 0;0 0 0 1;0 -0.002725 0.3283625 0];
b=[0;-0.883167;0;2.356553];
c=[1 0 0 0;0 0 1 0];
setlmis([]);
x=lmivar(1,[4 1]);
y=lmivar(1,[1 1]);
ew=lmivar(2,[1 1]);

lmiterm([1 1 1 x],1,a','s'); % LMI #1: x*a'+a*x
lmiterm([1 2 1 0],b'); % LMI #1: b'
lmiterm([1 2 2 y],1,-1); % LMI #1: -y
lmiterm([1 3 1 x],c,1); % LMI #1: c*x
lmiterm([1 3 3 0],-inv(ew)); % LMI #1: -inv(ew)

lmiterm([-2 1 1 x],1,1); % LMI #2: x 即令正定x>0，默认为<0，所以用-2表示第二个LMI大于0
lmiterm([-3 1 1 ew],1,1); % LMI #3: ew 即令α>0,是正实数，此处将α用ew表示，-3表示第三个LMI大于0
s1=getlmis;

[lp kp]=feasp(s1) %此处用来求s1的可行解，返回的lp<0，则该LMI可行。

x=dec2mat(s1,kp,x)