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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 这段matlab程序调试不出来结果，求助。。。。
北京石油化工学院2026年研究生招生接收调剂公告
31/1返回列表
查看: 620  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

孙家宝贝

金虫 (小有名气)

[求助] 这段matlab程序调试不出来结果，求助。。。。

clear
setlmis([]);
E=[-1 0;0 0];
A=[-0.6 0.25;0.8 -0.1];
Ad=[-0.7 0.3;0.25 -0.1];
B=[1;2];
Bw=[-0.6 -0.2;0 -1];
D=[-1;1];
M=[0.2;0.2];
H1=[0.2 0.2];
H2=[0.2 0.2];
H3=[1 1];
C=[-1 2;0 -1];
Cd=[-2 0;8.1 -1];
Dw=[-2.5 0;0 1];

D1=[0.2 0.2];
G1=[0.2 0.2];
N1=[0 0;0 0];
N2=[0 0;0 0];
L=[0 0;0 0];
R=[0;1];
d=9.5;
d0=0;
a1=0.05;
a2=0.5;
a3=0.15;
a4=0.05;
a5=0.7;
a6=0.05;

P=lmivar(1,[2 1]);
Z=lmivar(1,[2 1]);
X=lmivar(1,[2 1]);
Q=lmivar(1,[2 1]);
S=lmivar(2,[2 1]);
Sd=lmivar(2,[2 1]);
W=lmivar(2,[1 2]);
V1=lmivar(1,[2 1]);
V2=lmivar(1,[2 1]);
V3=lmivar(2,[2 1]);
V4=lmivar(2,[2 1]);
K=lmivar(2,[1 2]);
lmiterm([1 1 3 X],1,E'*L);
lmiterm([1 1 3 V3],1,R'*Bw);
lmiterm([1 1 4 X],d,N1');
lmiterm([1 1 5 X],d,A');
lmiterm([1 1 5 -W],d,B');
lmiterm([1 1 5 X],-1,E');
lmiterm([1 1 6 X],1,A');
lmiterm([1 1 6 -W],1,B');
lmiterm([1 1 7 X],1,C');
lmiterm([1 1 7 -W],1,D');
lmiterm([1 1 8 V3],1,R'*B*D1);
lmiterm([1 1 17 X],1,E');
lmiterm([1 1 18 X],1,1);
lmiterm([1 1 19 X],1,G1');
lmiterm([1 1 20 V3],1,R'*M);
lmiterm([1 1 24 X],1,E');
lmiterm([1 1 25 X],1,N1');
lmiterm([1 1 26 -W],1,B');
lmiterm([1 1 26 X],1,A');
lmiterm([1 1 27 V3],1,R');

lmiterm([1 2 3 X],-1,E'*L);
lmiterm([1 2 3 V4],1,R'*Bw);
lmiterm([1 2 3 X],1,a6*(H2')*H3);
lmiterm([1 2 4 X],1,d*N2');
lmiterm([1 2 5 X],1,d*Ad');
lmiterm([1 2 6 X],1,Ad');
lmiterm([1 2 7 X],1,Cd');

lmiterm([1 2 8 V4],1,R'*B*D1);
lmiterm([1 2 9 X],1,H2');
lmiterm([1 2 11 V4],1,R'*M);
lmiterm([1 2 21 X],1,1);
lmiterm([1 2 22 X],1,H2');
lmiterm([1 2 23 V4],1,R'*M);
lmiterm([1 2 24 X],-1,E');
lmiterm([1 2 25 X],1,N2');
lmiterm([1 2 26 X],1,Ad');
lmiterm([1 2 27 V4],1,R');
lmiterm([1 3 3 0],-a5^(2));
lmiterm([1 3 3 0],a6*(H3')*H3);
lmiterm([1 3 4 0],d*L');
lmiterm([1 3 5 0],d*Bw');
lmiterm([1 3 6 0],Bw');
lmiterm([1 3 7 0],Dw');
lmiterm([1 3 9 0],H3');
lmiterm([1 3 13 0],H3')
lmiterm([1 4 4 Z],d,1);
lmiterm([1 5 5 V2],-d,1);
lmiterm([1 5 10 0],B*D1);
lmiterm([1 5 11 0],d*M);
lmiterm([1 5 16 0],d*M);
lmiterm([1 6 6 X],1,1);
lmiterm([1 6 8 0],B*D1);
lmiterm([1 6 11 0],M);
lmiterm([1 6 16 0],M);
lmiterm([1 7 7 0],-1);
lmiterm([1 8 8 0],-a3);
lmiterm([1 9 9 0],-a3);
lmiterm([1 10 10 0],-a3);
lmiterm([1 11 11 0],-a4);
lmiterm([1 12 12 0],-a4);
lmiterm([1 13 13 0],-a4);
lmiterm([1 14 14 0],-a6);
lmiterm([1 15 15 0],-a6);
lmiterm([1 16 16 0],-a6);
lmiterm([1 17 17 X],1,1);
lmiterm([1 18 18 V1],-1,(d-d0+1)^(-1));
lmiterm([1 19 19 0],(a3*(1+d*d))^(-1));
lmiterm([1 20 20 0],-1);
lmiterm([1 21 21 V1],1,1);
lmiterm([1 22 22 0],-a4^(-1));
lmiterm([1 23 23 0],-a4^(-1));
lmiterm([1 24 24 0],-a1^(-1));
lmiterm([1 25 25 0],-a1);
lmiterm([1 26 26 0],-a2^(-1));
lmiterm([1 27 27 0],-a2);



lmiterm([-2 1 1 P],1,1)
lmiterm([2 1 1 0],0)
lmiterm([-3 1 1 Z],1,1)
lmiterm([3 1 1 0],0)
lmiterm([-4 1 1 X],1,1)
lmiterm([4 1 1 0],0)
lmiterm([-5 1 1 Q],1,1)
lmiterm([5 1 1 0],0)

lmis=getlmis
[tmin xfeas]=feasp(lmis)
%P=dec2mat(lmis,xfeas,P)
%Q=dec2mat(lmis,xfeas,Q)
%Z=dec2mat(lmis,xfeas,Z)
X=dec2mat(lmis,xfeas,X)
%Sd=dec2mat(lmis,xfeas,Sd)
%S=dec2mat(lmis,xfeas,S)
W=dec2mat(lmis,xfeas,W)
V1=dec2mat(lmis,xfeas,V1)
V2=dec2mat(lmis,xfeas,V2)
V3=dec2mat(lmis,xfeas,V3)
V4=dec2mat(lmis,xfeas,V4)
%X=inv(P);
%V1=inv(Q);
%V2=inv(Z);
%V3=X*S;
%V4=X*Sd;
%W=K*inv(P);


K=W*inv(X)
问题：运行后，无论如何调节已知参数，tmin都不能变成负数。也就是说，这个线性矩阵不等式没有可行解。这是为什么？求指导。。。。。。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
相信自己就能成功
1楼 2012-12-19 12:28:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙家宝贝

金虫 (小有名气)
clear
setlmis([]);
E=[-1 0;0 0];
A=[-0.6 0.25;0.8 -0.1];
Ad=[-0.7 0.3;0.25 -0.1];
B=[1;2];
Bw=[-0.6 -0.2;0 -1];
D=[-1;1];
M=[0.2;0.2];
H1=[0.2 0.2];
H2=[0.2 0.2];
H3=[1 1];
C=[-1 2;0 -1];
Cd=[-2 0;8.1 -1];
Dw=[-2.5 0;0 1];

D1=[0.2 0.2];
G1=[0.2 0.2];
N1=[0 0;0 0];
N2=[0 0;0 0];
L=[0 0;0 0];
R=[0;1];
d=9.5;
d0=0;
a1=0.05;
a2=0.5;
a3=0.15;
a4=0.05;
a5=0.7;
a6=0.05;

P=lmivar(1,[2 1]);
Z=lmivar(1,[2 1]);
X=lmivar(1,[2 1]);
Q=lmivar(1,[2 1]);
S=lmivar(2,[2 1]);
Sd=lmivar(2,[2 1]);
W=lmivar(2,[1 2]);
V1=lmivar(1,[2 1]);
V2=lmivar(1,[2 1]);
V3=lmivar(2,[2 1]);
V4=lmivar(2,[2 1]);
K=lmivar(2,[1 2]);
lmiterm([1 1 3 X],1,E'*L);
lmiterm([1 1 3 V3],1,R'*Bw);
lmiterm([1 1 4 X],d,N1');
lmiterm([1 1 5 X],d,A');
lmiterm([1 1 5 -W],d,B');
lmiterm([1 1 5 X],-1,E');
lmiterm([1 1 6 X],1,A');
lmiterm([1 1 6 -W],1,B');
lmiterm([1 1 7 X],1,C');
lmiterm([1 1 7 -W],1,D');
lmiterm([1 1 8 V3],1,R'*B*D1);
lmiterm([1 1 17 X],1,E');
lmiterm([1 1 18 X],1,1);
lmiterm([1 1 19 X],1,G1');
lmiterm([1 1 20 V3],1,R'*M);
lmiterm([1 1 24 X],1,E');
lmiterm([1 1 25 X],1,N1');
lmiterm([1 1 26 -W],1,B');
lmiterm([1 1 26 X],1,A');
lmiterm([1 1 27 V3],1,R');

lmiterm([1 2 3 X],-1,E'*L);
lmiterm([1 2 3 V4],1,R'*Bw);
lmiterm([1 2 3 X],1,a6*(H2')*H3);
lmiterm([1 2 4 X],1,d*N2');
lmiterm([1 2 5 X],1,d*Ad');
lmiterm([1 2 6 X],1,Ad');
lmiterm([1 2 7 X],1,Cd');

lmiterm([1 2 8 V4],1,R'*B*D1);
lmiterm([1 2 9 X],1,H2');
lmiterm([1 2 11 V4],1,R'*M);
lmiterm([1 2 21 X],1,1);
lmiterm([1 2 22 X],1,H2');
lmiterm([1 2 23 V4],1,R'*M);
lmiterm([1 2 24 X],-1,E');
lmiterm([1 2 25 X],1,N2');
lmiterm([1 2 26 X],1,Ad');
lmiterm([1 2 27 V4],1,R');
lmiterm([1 3 3 0],-a5^(2));
lmiterm([1 3 3 0],a6*(H3')*H3);
lmiterm([1 3 4 0],d*L');
lmiterm([1 3 5 0],d*Bw');
lmiterm([1 3 6 0],Bw');
lmiterm([1 3 7 0],Dw');
lmiterm([1 3 9 0],H3');
lmiterm([1 3 13 0],H3')
lmiterm([1 4 4 Z],d,1);
lmiterm([1 5 5 V2],-d,1);
lmiterm([1 5 10 0],B*D1);
lmiterm([1 5 11 0],d*M);
lmiterm([1 5 16 0],d*M);
lmiterm([1 6 6 X],1,1);
lmiterm([1 6 8 0],B*D1);
lmiterm([1 6 11 0],M);
lmiterm([1 6 16 0],M);
lmiterm([1 7 7 0],-1);
lmiterm([1 8 8 0],-a3);
lmiterm([1 9 9 0],-a3);
lmiterm([1 10 10 0],-a3);
lmiterm([1 11 11 0],-a4);
lmiterm([1 12 12 0],-a4);
lmiterm([1 13 13 0],-a4);
lmiterm([1 14 14 0],-a6);
lmiterm([1 15 15 0],-a6);
lmiterm([1 16 16 0],-a6);
lmiterm([1 17 17 X],1,1);
lmiterm([1 18 18 V1],-1,(d-d0+1)^(-1));
lmiterm([1 19 19 0],(a3*(1+d*d))^(-1));
lmiterm([1 20 20 0],-1);
lmiterm([1 21 21 V1],1,1);
lmiterm([1 22 22 0],-a4^(-1));
lmiterm([1 23 23 0],-a4^(-1));
lmiterm([1 24 24 0],-a1^(-1));
lmiterm([1 25 25 0],-a1);
lmiterm([1 26 26 0],-a2^(-1));
lmiterm([1 27 27 0],-a2);



lmiterm([-2 1 1 X],1,1)
lmiterm([2 1 1 0],0)
lmiterm([-3 1 1 V2],1,1)
lmiterm([3 1 1 0],0)
lmiterm([-4 1 1 V1],1,1)
lmiterm([4 1 1 0],0)
%lmiterm([-5 1 1 Q],1,1)
%lmiterm([5 1 1 0],0)

lmis=getlmis
[tmin xfeas]=feasp(lmis)
%P=dec2mat(lmis,xfeas,P)
%Q=dec2mat(lmis,xfeas,Q)
%Z=dec2mat(lmis,xfeas,Z)
X=dec2mat(lmis,xfeas,X)
%Sd=dec2mat(lmis,xfeas,Sd)
%S=dec2mat(lmis,xfeas,S)
W=dec2mat(lmis,xfeas,W)
V1=dec2mat(lmis,xfeas,V1)
V2=dec2mat(lmis,xfeas,V2)
V3=dec2mat(lmis,xfeas,V3)
V4=dec2mat(lmis,xfeas,V4)
%X=inv(P);
%V1=inv(Q);
%V2=inv(Z);
%V3=X*S;
%V4=X*Sd;
%W=K*inv(P);


K=W*inv(X)
应该是这个程序，上面的弄错了
赞一下 回复此楼
相信自己就能成功
2楼2012-12-19 12:37:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jaywu003

新虫 (初入文坛)

【答案】应助回帖

★
csgt0: 金币+1, 欢迎交流 2013-01-08 17:00:50
不等式有问题！第一个LMI 不会是负定的!
赞一下 回复此楼
做一只快乐的小虫！
3楼2013-01-08 14:30:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 孙家宝贝 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物与医药调剂 +4 十七sa 2026-04-05 4/200 2026-04-05 20:05 by lys0704
[考研] 307分材料专业求调剂 +7 Hll胡 2026-04-05 7/350 2026-04-05 18:47 by 无际的草原
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[考研] 296材料专硕求调剂 +21 202451007219 2026-04-02 22/1100 2026-04-04 21:48 by hemengdong
[考研] 26调剂 086003 +6 失活的细胞 2026-04-04 6/300 2026-04-04 09:50 by zhangdingwa
[考研] 265求调剂 +20 梁梁校校 2026-04-01 21/1050 2026-04-04 00:38 by userper
[考研] 317分 一志愿江南大学 化学工程学硕 求调剂 +6 YinTai 2026-04-03 6/300 2026-04-03 22:30 by 无际的草原
[考研] 320调剂 +4 农业工程与信息� 2026-04-03 4/200 2026-04-03 21:40 by lbsjt
[考研] 279求调剂 +6 qazplm0852 2026-04-02 6/300 2026-04-03 10:03 by 蓝云思雨
[考研] 一志愿深大085601材料工程专业（专硕）300分可以调剂去哪 +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[考研] 调剂 +7 祉岷. 2026-04-02 7/350 2026-04-03 09:11 by 花呗还欠600
[考研] 366求调剂一志愿东北大学 +8 运气来得若有似� 2026-04-02 8/400 2026-04-02 21:39 by dongzh2009
[论文投稿] chinese chemical letters英文版投稿求助 120+4 Yishengeryi 2026-03-30 6/300 2026-04-02 17:19 by Yishengeryi
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 13:03 by yulian1987
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 285求调剂 +7 AZMK 2026-03-30 13/650 2026-04-01 17:00 by 七度不信任
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 物理学调剂 +4 小羊36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
[考研] 材料专硕 085600求调剂 +7 BBQ233 2026-03-30 7/350 2026-03-30 17:44 by oooqiao
信息提示
关闭
请填处理意见
关闭