| 查看: 510 | 回复: 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; |
回复此楼» 猜你喜欢
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
-
什么时候出结果，有咨询渠道？？？
已经有3人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有26人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有7人回复
100