版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

孙家宝贝

金虫 (小有名气)

应助: 4 (幼儿园)
金币: 471.4
散金: 455
帖子: 89
在线: 62.2小时
虫号: 1439436
注册: 2011-10-12
性别: MM
专业: 控制论中的数学方法

[求助] 这段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都不能变成负数。也就是说，这个线性矩阵不等式没有可行解。这是为什么？求指导。。。。。。

回复此楼

» 猜你喜欢

帮我的英语口语老师找学生已经有4人回复
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值已经有7人回复
物理学I论文润色/翻译怎么收费? 已经有69人回复
2026年循环经济功能材料国际会议（ICFMCE 2026）已经有0人回复
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）已经有0人回复
哈尔滨理工大学物理系招收物理学考研调剂已经有5人回复
0702一志愿吉大B区求调剂已经有5人回复
求调剂已经有0人回复
0702一志愿吉大B区求调剂有论文已经有0人回复
请问还有没有用Latex写文章的小伙伴们？已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab程序出错调试已经有7人回复
matlab画图求助已经有11人回复
求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟已经有10人回复
matlab调试问题已经有3人回复
求助matlab编程已经有13人回复
【求助】matlab 中，几个m文件调用已经有14人回复
【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助已经有16人回复
【求助】用matlab编程，需要将数据自检分类（高手进来指点下）已经有5人回复
【求助】请教一个matlab程序，谢谢帮忙！！已经有15人回复
【求助】matlab中如何将程序的结果输出到一个文件中已经有11人回复
【求助】VS2005数值计算程序调试【急求！！！！加金币】已经有9人回复
【求助】请教把java程序转化成matlab程序已经有9人回复
【求助】初学matlab 已经有11人回复
【求助】请教虫友关于matlab差分方程的求解和画图程序已经有4人回复
【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】已经有15人回复
【求助】求MAtlab求解程序！！！已经有20人回复
【求助】正态分布如何用MATLAB检验已经有6人回复
【求助】帮我编个程啊~C语言，C++，Matlab，VB...都行~ 已经有14人回复
【求助】求HMM的Matlab程序已经有4人回复

相信自己就能成功

1楼 2012-12-19 12:28:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

孙家宝贝

金虫 (小有名气)

应助: 4 (幼儿园)
金币: 471.4
散金: 455
帖子: 89
在线: 62.2小时
虫号: 1439436
注册: 2011-10-12
性别: MM
专业: 控制论中的数学方法

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

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 84.7
帖子: 23
在线: 9小时
虫号: 1401831
注册: 2011-09-14
专业: 控制理论与方法

【答案】应助回帖

★
csgt0: 金币+1, 欢迎交流 2013-01-08 17:00:50

不等式有问题！第一个LMI 不会是负定的!

赞一下

回复此楼

做一只快乐的小虫！

3楼2013-01-08 14:30:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主孙家宝贝的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 326求调剂 +4	诺贝尔化学奖觊� 2026-03-15	4/200	2026-03-16 16:47 by 我的船我的海
[考研] 环境工程调剂 +5	大可digkids 2026-03-16	5/250	2026-03-16 16:45 by houyaoxu
[考研] 0703化学调剂，六级已过，有科研经历 +7	曦熙兮 2026-03-15	7/350	2026-03-16 16:34 by houyaoxu
[考研] 285求调剂 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 304求调剂 +6	小熊joy 2026-03-14	6/300	2026-03-16 12:59 by Iveryant
[考研] 326求调剂 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 材料与化工（0856）304求B区调剂 +7	邱gl 2026-03-10	11/550	2026-03-14 12:18 by 邱gl
[考研] 一志愿郑大070303，338分，求调剂 +4	dadawaf 2026-03-10	5/250	2026-03-14 01:20 by lsw010101
[考研] 一志愿华中农业大学071010，总分三百二,求调剂 +3	困困困困坤坤 2026-03-10	3/150	2026-03-14 00:35 by JourneyLucky
[考研] 318求调剂 +3	李新光 2026-03-10	3/150	2026-03-14 00:21 by JourneyLucky
[考研] 329求调剂 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 求调剂 +7	18880831720 2026-03-11	7/350	2026-03-13 16:10 by JourneyLucky
[考研] 研究生招生 +3	徐海涛11 2026-03-10	7/350	2026-03-12 14:26 by 徐海涛11
[考研] 085600 材料与化工 295 求调剂 +10	dream…… 2026-03-10	12/600	2026-03-12 13:46 by dream……
[考博] 读博申请 +5	感dd 2026-03-10	7/350	2026-03-11 17:02 by QGZDSYS
[考研] 哈工大材料324求调剂 +6	闫旭东 2026-03-10	8/400	2026-03-10 22:49 by 星空星月
[考博] 26申博求助 +3	跳跃饼干 2026-03-10	4/200	2026-03-10 21:15 by Tntcnn
[考研] 0703化学调剂 +3	三dd. 2026-03-10	3/150	2026-03-10 15:45 by peike
[考研] 一志愿：武汉理工，材料工程，英二数二总分314 +3	2202020125 2026-03-10	4/200	2026-03-10 13:54 by xiongyaxuan
[硕博家园] 木虫好像不热闹了，是不是？ +4	偏振片 2026-03-10	4/200	2026-03-10 09:51 by longwave