应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 信息科学 » 程序算法 » MATLAB中LMI程序求解H无穷可行解问题
21/1返回列表
查看: 4231  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

_拧发条的鸟

新虫 (初入文坛)

[求助] MATLAB中LMI程序求解H无穷可行解问题

MATLAB中LMI程序运行报错：Error using lmiterm (line 271)lhs of LMI #4, block (1,1): terms in diagonal blocks must be square；
程序如下：
clc
clear
A=[0.0000000e+00   3.6301310e-21  -9.7026387e+00  -4.9849839e-05   0.0000000e+00  4.9698601e-07   0 0 0 0 0 0 0;
      7.9257028e-24   0.0000000e+00   7.8831201e+03   0.0000000e+00   0.0000000e+00 0 0 0 0 0 0 0 0;
         3.1226519e-07   1.9087724e-10   0.0000000e+00   1.3008627e-08   0.0000000e+00 1.5307304e-11   0.0000000e+00 0 0 0 0 0 0;
           -3.1226519e-07  -1.9087724e-10   0.0000000e+00  -1.3008627e-08   1.0000000e+00 -1.5307304e-11   0 0 0 0 0 0 0;
             -3.1895452e-23   5.2382432e-23   0.0000000e+00  -5.0419522e-07  -4.5526839e-08 0 0 0 0 0 0 0 0;
             0 0 0 0 0 0 1 0 0 0 0 0 0;
             0 0 0 0 0 -0.49 -28 0 0 0 0 0 0;
             0 0 0 0 0 0 0 0 0 0 0 0 0;
             0 0 0 0 0 0 0 0 0 0 0 0 0;
             0 0 0 0 0 0 0 0 0 0 0 1 0;
             0 0 0 0 0 0 0 0 0 0 0 0 1;
             0 0 0 0 0 0 0 1 0 -1 0 -0.16 0;
             0 0 0 0 0 0 0 0 1 0 -1 0 -0.16];
B=[0 0;
    0 0;
    0 0;
    0 0;
    0 1.5905035e-06;
    0 0;
    0.49 0;
    0 0;
    0 0;
    0 0;
    0 0;
    0 0];
G=[1 0 0 0 0;
   0 1 0 0 0;
   0 0 1 0 0;
   0 0 0 1 0;
   0 0 0 0 1;
   0 0 0 0 0;
   0 0 0 0 0;
   0 0 0 0 0;
   0 0 0 0 0;
   0 0 0 0 0;
   0 0 0 0 0;
   0 0 0 0 0;
   0 0 0 0 0];
F=[1 0 0 0 0 0 0 0 0 0 0 0 0
   0 1 0 0 0 0 0 0 0 0 0 0 0
   0 0 1 0 0 0 0 0 0 0 0 0 0
   0 0 0 1 0 0 0 0 0 0 0 0 0
   0 0 0 0 1 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0];
E1=[1 0 0 0 0 0 1 0 0 0 0 0 0
   0 1 0 0 0 0 0 0 0 0 0 0 0
   0 0 1 0 0 0 0 0 0 0 0 0 0
   0 0 0 1 0 0 0 0 0 0 0 0 0
   0 0 0 0 1 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 0 0 0 0 0];
E2=[0 1;
    0 0;
     0 0;
     0 0;
     0 0;
     0 0;
     0 0;
     0 0;
     0 0;
     0 0;
     0 0;
     0 0;
     0 0];
C=[0 0 1 0 0 0 -1 0 0 0 0 0 0;
   0 0 0 1 0 0 0 -1 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 ];
D=[0 0;
    0 0;
    1 0;
    0 1];

setlmis([])
Q=lmivar(1,[13 1]);
Y=lmivar(2,[2 13]);
r=lmivar(1,[1 1]);%
r1=lmivar(1,[1 1]);%
r2=lmivar(1,[1 1]);%

%strictly proper
lmiterm([-1 1 1 Q],1,1);
lmiterm([-2 1 1 r],1,1);
lmiterm([-3 1 1 r2],1,1);
lmiterm([4 1 1 Q],A,1,'s');
lmiterm([4 1 1 Y],B,1,'s');
lmiterm([4 1 1 r],1,F*F');
lmiterm([4 2 1 Q],E1,1);
lmiterm([4 2 1 Y],E2,1);
lmiterm([4 3 1 Q],C,1);
lmiterm([4 3 1 Y],D,1);
lmiterm([4 4 1 0],G');
lmiterm([4 2 2 r],-1,1);
lmiterm([4 3 3 0],-1);
lmiterm([4 4 4 r2],-1,1);
lmisys=getlmis;

[ropt xopt]=feasp(lmisys);%
r=dec2mat(lmisys,xopt,r);
r2=dec2mat(lmisys,xopt,r2);
Q=dec2mat(lmisys,xopt,Q);
Y=dec2mat(lmisys,xopt,Y);
K=Y*inv(Q);

运行时报错，但是矩阵维数应该是没有问题，求大神帮助，LMI不等式如下

MATLAB中LMI程序求解H无穷可行解问题
%@Y4C)Q1IUU4{KC@901WF(R.png
回复此楼

» 猜你喜欢
1楼 2015-12-14 10:59:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

云月禅心

铜虫 (小有名气)
理解，但是爱莫能助呀
赞一下 回复此楼
撸起袖子加油干
2楼2016-05-18 10:57:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 _拧发条的鸟 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +7 枫桥ZL 2026-03-18 9/450 2026-03-21 09:56 by JineShine
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 324求调剂 +5 lucky呀呀呀鸭 2026-03-20 5/250 2026-03-20 22:30 by 促天成
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
信息提示
关闭
请填处理意见
关闭