声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2095|回复: 4

[控制理论] [求助]请问LIM(线性矩阵不等式)的一个程序问题

[复制链接]
发表于 2007-9-28 00:03 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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 编辑 ]
回复
分享到:

使用道具 举报

 楼主| 发表于 2007-9-28 09:33 | 显示全部楼层
晕,是LMI,竟然打成LIM,抱歉抱歉
发表于 2007-9-29 11:11 | 显示全部楼层
第二个lmiterm项怎么是lmiterm([-1 1 1 0],p*a),而不是lmiterm([-1 1 1 0],a)
 楼主| 发表于 2007-9-29 23:10 | 显示全部楼层
不好意思啊,原矩阵应该是
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


麻烦你们看看那里出了问题,怎么结果相差很大啊。
多谢多谢!!
发表于 2007-12-13 20:12 | 显示全部楼层
,结果是
KK1 =

  5.2473e+007


KK2 =

  5.1174e+008


KK3 =

  5.1174e+008
不对吗?:@)
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-1 07:47 , Processed in 0.067345 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表