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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 这段matlab程序调试不出来结果，求助。。。。
北京石油化工学院2026年研究生招生接收调剂公告
31/1返回列表
查看: 621  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂，一志愿郑州大学材料与化工专硕，英二数二342分，求老师收留 +19 v12abo 2026-04-02 21/1050 2026-04-06 09:29 by 蓝云思雨
[考研] 302分求调剂 一志愿安徽大学085601 +8 zyx上岸！ 2026-04-04 8/400 2026-04-05 23:08 by lbsjt
[考研] 085600，320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 工科08-机械专硕-求调剂 +3 雷欧飞踢 2026-04-02 3/150 2026-04-05 18:49 by 蓝云思雨
[考研] 298求调剂 +7 manman511 2026-04-05 7/350 2026-04-05 10:29 by 唐沐儿
[考研] 359求调剂22408 +3 123456789qw 2026-03-31 3/150 2026-04-05 10:09 by zhq0425
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 一志愿双非085502，267分，过四级求调剂 +3 再忙也要吃饭啊 2026-04-03 3/150 2026-04-04 05:03 by gswylq
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 生物学308分求调剂（一志愿华东师大） +7 相信必会光芒万� 2026-04-02 7/350 2026-04-03 16:48 by rzh123456
[考研] 334求调剂 +9 Trying] 2026-03-31 9/450 2026-04-03 15:18 by 琢珥丶
[考研] 283求调剂 +3 jiouuu 2026-04-03 4/200 2026-04-03 13:28 by jiouuu
[考研] 一志愿武汉理工0856，初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
[考研] 一志愿郑大材料工程290求调剂 +20 Youth_ 2026-03-30 20/1000 2026-04-02 14:48 by 5896
[考研] 0710生物学，325求调剂 +3 mkkkkkl 2026-04-01 3/150 2026-04-02 09:48 by Jaylen.
[考研] 08生物与医药专硕初试346找调剂 +6 dianeeee 2026-04-01 7/350 2026-04-02 08:23 by guoweigw
[考研] 350求调剂 +7 阿佳～ 2026-03-31 7/350 2026-04-01 16:12 by yanflower7133
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
信息提示
关闭
请填处理意见
关闭