| 查看: 2456 | 回复: 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))+... 求大神帮忙，实在不知道怎么搞  |
回复此楼» 猜你喜欢
296求调剂
已经有5人回复
-
0703化学调剂
已经有13人回复
-
085600材料与化工 求调剂
已经有15人回复
-
085700资源与环境308求调剂
已经有4人回复
-
生物学调剂招人！！！
已经有4人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有17人回复
-
一志愿武汉理工材料工程专硕调剂
已经有5人回复
-
0856调剂，是学校就去
已经有7人回复
-
申博26年
已经有3人回复
-
一志愿吉林大学材料学硕321求调剂
已经有10人回复
2楼2017-03-31 10:10:53
