24小时热门版块排行榜    

查看: 462  |  回复: 0

cugbwangyang

新虫 (初入文坛)

[求助] LMI工具箱编程求助

本人初学MATLAB的LMI工具箱,请问哪位大神指导下面的程序为什么没有解呢?是编程出问题了吗?这个是一篇文章里的例子,应该是有解的啊。谢谢大家了!

A1=[0 1.0000 0 0; 1.9600 0 -0.4000 -0.6000; 0 0 0 1.0000; 1.9600 0 -2.4000 -3.6000];
A2=[0 1.0000 0 0; 0 0 -0.4000 -0.6000; 0 0 0 1.0000; 0 0 -2.4000 -3.6000];
B1=[0; 0; 0; -1];
B2=[0; 0; 0; -1];

x1=[1.0000 0 0 0]';

setlmis([]);%初始化
%定义矩阵变量
X=lmivar(2,[4,4]);
M1=lmivar(2,[1,4]);
M2=lmivar(2,[1,4]);
K1=lmivar(2,[1,4]);
K2=lmivar(2,[1,4]);
P=lmivar(2,[4,4]);

%描述第一个矩阵
lmiterm([1 1 1 X],A1,1,'s');
lmiterm([1 1 1 M1],B1,1,'s');
%描述第二个矩阵
lmiterm([2 1 1 X],A2,1,'s');
lmiterm([2 1 1 M2],B2,1,'s');
%描述第三个矩阵
lmiterm([3 1 1 X],A1,1,'s');
lmiterm([3 1 1 M2],B1,1,'s');   
lmiterm([3 1 1 X],A2,1,'s');
lmiterm([3 1 1 M1],B2,1,'s');

lmisys=getlmis;
[r,xfeas]=gevp(lmisys,3,[0, 0, 0, 0, 1]);
X=dec2mat(lmisys, xfeas, X);
M1=dec2mat(lmisys, xfeas, M1);
M2=dec2mat(lmisys, xfeas, M2);
P=inv(X);
K1=M1*P;
K2=M2*P;
回复此楼

» 猜你喜欢

一切都会好起来的,即便不是在今天!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cugbwangyang 的主题更新
信息提示
请填处理意见