24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1730  |  回复: 7

yue_s008

至尊木虫 (职业作家)

[求助] 为什么 LMI 的解不满足设定的不等式条件(附程序)已有1人参与

A=[0.8,1.2,-4;0,0.8,0.8;0,0,1];
E=[-1,0,-40;0,1,-4];

setlmis([])
P=lmivar(1,[3,1]);U=lmivar(1,[2,1]);

lmiterm([-1,1,1,P],1,1);
lmiterm([-1,1,1,U],E',-E)


lmis=getlmis;
[tmin,feas]=feasp(lmis);

我这样得出的解为什么不满足0<P-E'*U*E的条件
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kuikui9802

金虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
你两个表达式里面都是负1,应该有一个正一
lmiterm([-1,1,1,P],1,1);
lmiterm([1,1,1,U],E',-E)
2楼2018-07-16 08:19:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yue_s008

至尊木虫 (职业作家)

引用回帖:
2楼: Originally posted by kuikui9802 at 2018-07-16 08:19:00
你两个表达式里面都是负1,应该有一个正一
lmiterm(,1,1);
lmiterm(,E',-E)

你再好好看看,你这个和我原来的不等式都不是同一个了,减号变加号了。

发自小木虫IOS客户端
3楼2018-07-16 08:56:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yue_s008

至尊木虫 (职业作家)

引用回帖:
2楼: Originally posted by kuikui9802 at 2018-07-16 08:19:00
你两个表达式里面都是负1,应该有一个正一
lmiterm(,1,1);
lmiterm(,E',-E)

你这个是 -E'*U*E<P
4楼2018-07-16 08:57:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kuikui9802

金虫 (著名写手)

哦,没注意,那应该没错啊
5楼2018-07-16 09:30:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yue_s008

至尊木虫 (职业作家)

引用回帖:
5楼: Originally posted by kuikui9802 at 2018-07-16 09:30:30
哦,没注意,那应该没错啊

我就用了两个类似的lmi,就得到不严格可行的答案了。

发自小木虫IOS客户端
6楼2018-07-16 12:56:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
7楼2018-10-15 12:20:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1124127307

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by yue_s008 at 2018-07-16 12:56:44
我就用了两个类似的lmi,就得到不严格可行的答案了。
...

不严格可行解,确实不满足约束条件。所以必须求解到可行解,即那个时间 t<0. (t越小,解越严格,越好)
8楼2019-04-12 15:37:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yue_s008 的主题更新
信息提示
请填处理意见