查看: 107|回复: 0|关注: 0

[未答复] 求解LMI问题时计算出错

[复制链接]

新手

5 麦片

财富积分


050


7

主题

12

帖子

0

最佳答案
发表于 5 天前 | 显示全部楼层 |阅读模式
miu=0.1;
h=0.3796;
A=[0,1;-8.883,-0.314];B=[0,0;-0.987,0];
System_Dimension = length(A); %系统维数
setlmis([]);
P   = lmivar(1, [System_Dimension 1]);
Q  = lmivar(1, [System_Dimension 1]);
R   = lmivar(1, [System_Dimension 1]);
W  = lmivar(1, [System_Dimension 1]);
Z   = lmivar(2, [System_Dimension System_Dimension]);
Power1 = newlmi;
lmiterm([Power1 1 1 P  ], 1    , A);
lmiterm([Power1 1 1 P  ], A'   , 1);
lmiterm([Power1 1 1 Q  ], 1    , 1);
lmiterm([Power1 1 1 W  ], 1    , 1);
lmiterm([Power1 1 1 R  ], -1   , 1);
lmiterm([Power1 1 2 P  ], 1    , B);
lmiterm([Power1 1 2 R  ], 1    , 1);
lmiterm([Power1 1 2 Z  ], 1    , 1);
lmiterm([Power1 1 3 Z  ], -1   , 1);
lmiterm([Power1 1 4 R  ], h*A'    , 1);
lmiterm([Power1 2 2 W  ], (-1+miu)  , B);
lmiterm([Power1 2 2 R  ], -2   , 1);
lmiterm([Power1 2 2 Z  ], -1    , 1);
lmiterm([Power1 2 2 Z'  ], -1    , 1);
lmiterm([Power1 2 3 R  ], 1   , 1);
lmiterm([Power1 2 3 Z  ], 1    , 1);
lmiterm([Power1 2 4 R  ], h*B'    , 1);
lmiterm([Power1 3 3 Q  ], -1   , 1      );
lmiterm([Power1 3 3 R  ],  -1   , 1     );

lmiterm([Power1 4 4 R  ],  -1       , 1);
Power2 = newlmi;
lmiterm([-Power2 1 1 R  ], 1    , 1);
lmiterm([-Power2 1 2 Z  ], 1    , 1);
lmiterm([-Power2 2 2 R  ], 1    , 1);
PLmi = newlmi;
      
lmiterm([-PLmi 1 1 P], 1, 1);
QLmi = newlmi;
lmiterm([-QLmi 1 1 Q], 1, 1);
RLmi = newlmi;
lmiterm([-RLmi 1 1 R], 1, 1);
WLmi = newlmi;
lmiterm([-WLmi 1 1 W], 1, 1);
LMI_System = getlmis;
Options(5) = 1;
Options(2) = 100;
Options(3) = 1e9;
Options(4) = 10;
[Goal, Feasible_Solution] = feasp(LMI_System, Options);
if(Goal<0)
    Result = 1
else
    Result = 0
end

当miu=0.1的时候,求解图片中的LMI,我编的程序如上。正确答案应该是h=0.3796,但我始终算不出来,有没有大神能帮我看一下程序呀?
微信截图_20190112105833.png
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /2 下一条

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