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

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有165人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有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)
应该是这个程序，上面的弄错了