24小时热门版块排行榜

返回列表

cugbwangyang

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 93.4
帖子: 23
在线: 12.9小时
虫号: 3512846
注册: 2014-11-01
性别: GG
专业: 机构学与机器人

[求助] 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;

回复此楼

» 猜你喜欢

一切都会好起来的，即便不是在今天！

1楼 2015-08-24 19:42:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 cugbwangyang 的主题更新

返回列表