|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
哪位高手能帮我看看这个问题!
clear
close
clc
A=[-2.73 4.26;2.08 -3.12];
B=[1.27 0.72;0.45 1.18];
C=[1.20 2.50;0.80 1.80];
G=[1 2;0 5];
E1=[0.16 0 0.1 0 0 0;0 0.16 0 0.1 0 0];
E2=[0 0 0 0 0.12 0;0 0 0 0 0 0.15];
F1=[1 0;0 1;0 0;0 0;1 0;0 1];
F2=[0 0;0 0;1 0;0 1;0 0;0 0];
L=[0 0;0 1];
I=[1 0;0 1];
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,[2 1]);
Q=lmivar(1,[2 1]);
r=lmivar(1,[1 0]);
e1=lmivar(1,[1 0]);
e2=lmivar(1,[1 0]);
lmiterm([1 1 1 X],1,A','s');
lmiterm([1 1 1 Q],B*L,1,'s');
lmiterm([1 1 1 e1],1,E1*E1');
lmiterm([1 1 2 0],G);
lmiterm([1 2 2 r],-1,I);
lmiterm([1 3 1 X],C,1);
lmiterm([1 3 2 0],0);
lmiterm([1 3 3 r],-1,I);
lmiterm([1 3 3 e2],1,E2*E2');
lmiterm([1 4 1 X],F1,1);
lmiterm([1 4 1 Q],F2*L,1);
lmiterm([1 4 2 0],0);
lmiterm([1 4 3 0],0);
lmiterm([1 4 4 e2],-1,I1);
lmiterm([1 5 1 X],F1,1);
lmiterm([1 5 1 Q],F2*L,1);
lmiterm([1 5 2 0],0);
lmiterm([1 5 3 0],0);
lmiterm([1 5 4 0],0);
lmiterm([1 5 5 e1],-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=[1e-5,0,0,0,0];
[copt,xopt]=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 编辑 ] |
|