关于线性矩阵不等式的出错问题
哪位高手能帮我看看这个问题!clear
close
clc
A=[-2.73 4.26;2.08 -3.12];
B=;
C=;
G=;
E1=;
E2=;
F1=;
F2=;
L=;
I=;
I1=[1 0 0 0 0 0;
0 1 0 0 0 0;
0 0 1 0 0 0;
0 0 0 1 0 0;
0 0 0 0 1 0;
0 0 0 0 0 1];
setlmis([])
X=lmivar(1,);
Q=lmivar(1,);
r=lmivar(1,);
e1=lmivar(1,);
e2=lmivar(1,);
lmiterm(,1,A','s');
lmiterm(,B*L,1,'s');
lmiterm(,1,E1*E1');
lmiterm(,G);
lmiterm(,-1,I);
lmiterm(,C,1);
lmiterm(,0);
lmiterm(,-1,I);
lmiterm(,1,E2*E2');
lmiterm(,F1,1);
lmiterm(,F2*L,1);
lmiterm(,0);
lmiterm(,0);
lmiterm(,-1,I1);
lmiterm(,F1,1);
lmiterm(,F2*L,1);
lmiterm(,0);
lmiterm(,0);
lmiterm(,0);
lmiterm(,-1,I1);
Xpos=newlmi;
lmiterm([-Xpos 1 1 X],1,1);
rpos=newlmi;
lmiterm([-rpos 1 1 r],1,1);
e1pos=newlmi;
lmiterm([-rpos 1 1 e1],1,1);
e2pos=newlmi;
lmiterm([-rpos 1 1 e2],1,1);
uncertainHinf2=getlmis;
c=mat2dec(uncertainHinf2,eye(2),eye(2),1,1,1);
options=;
=mincx(uncertainHinf2,c,options);
X=dec2mat(uncertainHinf2,xopt,X)
Q=dec2mat(uncertainHinf2,xopt,Q)
r=dec2mat(uncertainHinf2,xopt,r)
e1=dec2mat(uncertainHinf2,xopt,e1)
e2=dec2mat(uncertainHinf2,xopt,e2)
K=Q*inv(X)
最后就是显示不了结果!显示
Result:could not establish feasibility nor infeasibility
The LMI constraints were found INFEASIBLE
??? Error using ==> dec2mat
DECVARS must be a vector of length 9
Error in ==> D:\Matlab\work\uncertainHinf2.m
On line 94==> X=dec2mat(uncertainHinf2,xopt,X)
不知道什么原因,用feasp的时候就好用为什么换成mincx的时候就出现错误!请高手指点指点,多谢
[ 本帖最后由 eight 于 2007-9-2 21:14 编辑 ] 原帖由 fang_317 于 2007-9-2 11:02 发表 http://www.chinavib.com/forum/images/common/back.gif
哪位高手能帮我看看这个问题!
clear
close
clc
A=[-2.73 4.26;2.08 -3.12];
B=;
C=;
G=;
E1=;
E2=
根据出错地方自己调试，详细描述请看看置顶帖：聚宝盆 我也遇到同样的问题
不知道你解决了没有.
请与我联系
QQ:15709027 楼上你好 我也遇到了相同的问题 请问怎么解决的？ 特别着急加我qq好吗？649013663 矩阵不等式没有可行解
相关可能原因及解决办法如下：
1）、先检查以下LMI有没有问题，比如推导过程、输入到MATLAB，也许是打字打错了。如果还不行，可能是理论有问题。不是所有正确的推导都有解出现；
2）、程序可以运行，说明语法没有问题。有可能需要检查一下LMI，检查矩阵的主对角线元素，若对角线元素有严格大于0存在，肯定是无可行解的；
3）、一般可行性问题的程序都有多个LMI，可逐个检查，一个LMI一个的检查是否可行，找到到底是哪个LMI不可行，再想办法解决。
页:
[1]