| 查看: 2516 | 回复: 1 | ||
mavisym金虫 (小有名气)
|
[求助]
求助，MATLAB报错
|
|
我的公式明明除的是实数，matlab报错如下： 错误使用 mupadmex Error in MuPAD command: Division by zero. [_power] Evaluating: symobj::trysubs 出错 sym/subs>mupadsubs (line 139) G = mupadmex('symobj::fullsubs',F.s,X2,Y2); 出错 sym/subs (line 124) G = mupadsubs(F,X,Y); 出错 Chebyshev (line 67) K1(alf,alf)=4*D*a*b.*(XGM(3,3).*PSI(1,1)+XGM(1,1).*PSI(3,3)+nju.*XGM(3,1).*PSI(1,3))+... 程序如下 clear clc a=2;b=1;h=0.1;lu=2.71; nju=0.3; D=20; %SSSS % kx0=5e7;Kx0=0;kxa=5e7;Kxa=0; ky0=5e7;Ky0=0;kyb=5e7;Kyb=0; %} N=5;M=5; K1=zeros((M+1)*(N+1),(M+1)*(N+1)); M1=zeros((M+1)*(N+1),(M+1)*(N+1)); Q=zeros((M+1)*(N+1),1); F=zeros((M+1)*(N+1),1); for n=1:N for m=1:M syms z e; Tn=cos(n*acos(z));Tm=cos(m*acos(e)); for tao=1:3 for tao1=1:3 if (tao==1&&tao1==1) X=Tn.*Tn; P=Tm.*Tm; elseif tao==1 X=Tn.*diff(Tn,tao1-1)/a.^(tao1-1); P=Tm.*diff(Tm,tao1-1)/b.^(tao1-1); elseif tao1==1 X=diff(Tn,tao-1)/a.^(tao-1).*Tn; P=diff(Tm,tao-1)/b.^(tao-1).*Tm; else X=diff(Tn,tao-1)/a.^(tao-1).*diff(Tn,tao1-1)/a.^(tao1-1); P=diff(Tm,tao-1)/b.^(tao-1).*diff(Tm,tao1-1)/b.^(tao1-1); end XGM(tao,tao1)=double(int(X,z,-1,1)); PSI(tao,tao1)=double(int(P,e,-1,1)); end end alf=(M+1)*(N+1); Tn1=diff(Tn,'z');Tm1=diff(Tm,'e'); K1(alf,alf)=4*D*a*b.*(XGM(3,3).*PSI(1,1)+XGM(1,1).*PSI(3,3)+nju.*XGM(3,1).*PSI(1,3))+... 4*D*a*b.*(nju.*XGM(1,3).*PSI(3,1)+2*(1-nju).*XGM(2,2).*PSI(2,2))+... (kx0*b.*(eval(subs(Tn1,z,-1))).^2/2).*PSI(1,1)+2*Kx0*b.*(eval(subs(Tn1,z,-1))).^2/a^2; % (kxa*b.*(eval(subs(Tn1,z,1))).^2/2).*PSI(1,1)+2*Kxa*b.*(eval(subs(Tn1,z,1))).^2/a^2+... % (ky0*a.*(eval(subs(Tm1,z,-1))).^2/2).*XGM(1,1)+2*Ky0*a.*(eval(subs(Tm1,z,-1))).^2/b^2+... % (kyb*a.*(eval(subs(Tm1,z,1))).^2/2).*XGM(1,1)+2*Kyb*a.*(eval(subs(Tm1,z,1))).^2/b^2; M1(alf,alf)=lu*h*a*b.*XGM(1,1).*PSI(1,1)/4; end end 我注释的三个公式本来是一个长的，分开找错误的，如果是个长的就报别的错 从 sym 转换为 double 时出现以下错误: 错误使用 mupadmex Error in MuPAD command: DOUBLE cannot convert the input expression into a double array. If the input expression contains a symbolic variable, use the VPA function instead. 出错 Chebyshev (line 67) K1(alf,alf)=4*D*a*b.*(XGM(3,3).*PSI(1,1)+XGM(1,1).*PSI(3,3)+nju.*XGM(3,1).*PSI(1,3))+... 求大神帮忙，实在不知道怎么搞  |
回复此楼» 猜你喜欢
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有15人回复
-
面上的水很深
已经有18人回复
-
国社科通讯评审
已经有5人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
生命口C13面上会评时间
已经有11人回复
-
内心匮乏
已经有16人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
伯胺上甲基求助
已经有7人回复
-
跨出版社商投稿
已经有5人回复
2楼2017-03-31 10:10:53