|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
各位高人,看到篇关于LMI的论文,想仿真结果,发现和论文的结果相差很多,请教高人我的程序中是否出现什么问题
,多谢!
已知:q=[ p*k1+a -0.5*(p*a-a+c+Mz) -0.5My
-0.5*(p*a-a+c+Mz) k2-c 0
-0.5My 0 k3+b ] >0
其中 a=35 b=3,c=28,Mz=53,My=32, p=My^2/(a*b), 求解k1,k2,k3
程序如下:
clear,clc
a=35;
b=3;
c=28;
mcy=32;
mcz=53;
p=mcy^2/(a*b);
setlmis([])
K1=lmivar(1,[1 1]);
K2=lmivar(1,[1 1]);
K3=lmivar(1,[1 1]);
lmiterm([-1 1 1 K1],p,1);
lmiterm([-1 1 1 0],p*a);
lmiterm([-1 1 2 0],-1/2*(p*a-a+c+mcz));
lmiterm([-1 1 3 0],-1/2*mcy);
lmiterm([-1 2 2 K2],1,1);
lmiterm([-1 2 2 0],-c);
lmiterm([-1 3 3 K3],1,1);
lmiterm([-1 3 3 0],b);
lmisys=getlmis;
[tmin,xfeas]=feasp(lmisys);
KK1=dec2mat(lmisys,xfeas,K1)
KK2=dec2mat(lmisys,xfeas,K2)
KK3=dec2mat(lmisys,xfeas,K3)
程序运行的结果和原论文有很多的差别,竟然相差10^3左右,望高人指导!!!
多谢多谢!!
[ 本帖最后由 xmwhit 于 2007-9-28 08:42 编辑 ] |
|