24小时热门版块排行榜    

查看: 539  |  回复: 1

yoyo0201

金虫 (著名写手)

[求助] 关于Matalb中LMI工具箱求值问题

(Ai+Bi*Gi)'P+P'(Ai+Bi*Gi)'<0,i=1,2,3;已知A1,A2,A3,B1,B2,B3;根据不等式用LMI工具箱求P,G1,G2,G3.程序如下,运行出错,请问问题出在哪里?
%%%%%赋值
A1=[0 0 0 1 -0.0349 0;0 0 0 0.0349 1 0;0 0 0 0 0 1;-0.0358 0 0 -0.0797 0 0;0 -0.0208 0 0 -0.0818 -0.1224;0 -0.0394 0 0 -0.2254 -0.2468];
A2=[0 0 0 0.0349 -1 0;0 0 0 1 0.0349 0;0 0 0 0 0 1;-0.0358 0 0 -0.0797 0 0;0 -0.0208 0 0 -0.0818 -0.1224;0 -0.0394 0 0 -0.2254 -0.2468];
A3=[0 0 0 0.0349 1 0;0 0 0 -1 0.0349 0;0 0 0 0 0 1;-0.0358 0 0 -0.0797 0 0;0 -0.0208 0 0 -0.0818 -0.1224;0 -0.0394 0 0 -0.2254 -0.2468];
B1=[0 0 0 0.9215 0 0;0 0 0 0 0.7802 1.4811;0 0 0 0 1.4811 7.4562]';
B2=B1;B3=B1;
%%%LMI%%%%%%
%初始化LMI
setlmis([]);
%定义变量
P=lmivar(1,[6,1]);
G1=lmivar(2,[3,6]);
G2=lmivar(2,[3,6]);
G3=lmivar(2,[3,6]);
%添加项
lmiterm([1 1 1 P],1,A1+B1*G1,'s');
lmiterm([2 1 1 P],1,A2+B2*G2,'s');
lmiterm([3 1 1 P],1,A3+B3*G3,'s');
%获取LMI系统描述
lmis=getlmis;
[tmin,xfeas]=feasp(lmis)
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iyangxt

新虫 (小有名气)

【答案】应助回帖

★ ★
xiegangmai: 金币+2, 谢谢参与 2013-01-06 14:44:54
这个问题解决了吗?
lmiterm([1 1 1 P],1,A1+B1*G1,'s');
lmiterm([2 1 1 P],1,A2+B2*G2,'s');
lmiterm([3 1 1 P],1,A3+B3*G3,'s');
有问题
天王盖地虎
2楼2012-12-31 16:51:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yoyo0201 的主题更新
信息提示
请填处理意见