版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (5383)
>考研 (1033)
>導師招生 (945)
>文獻求助 (384)
>蟲友互識 (251)
>基金申請 (185)
>休閑灌水 (170)
>論文投稿 (137)
>碩博家園 (117)
>考博 (112)
>招聘信息布告欄 (92)
>博后之家 (65)
>教師之家 (64)
>綠色求助(高懸賞) (38)
>公派出國 (34)
>留學生活 (32)

返回列表

mavisym

金蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 856.2
散金: 14
帖子: 58
在線: 19.5小時
蟲號: 3144160
注冊: 2014-04-14
性別: MM
專業(yè): 機械結構強度學

[求助] 求助，MATLAB報錯

我的公式明明除的是實數(shù)，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 時出現(xiàn)以下錯誤:
錯誤使用 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))+...

求大神幫忙，實在不知道怎么搞

回復此樓

» 猜你喜歡

一志愿武理材料305分求調劑已經有5人回復
材料專碩英一數(shù)二306 已經有3人回復
085600材料與化工調劑 324分已經有7人回復
311求調劑已經有3人回復
26調劑/材料/英一數(shù)二/總分289/已過A區(qū)線已經有7人回復
295求調劑已經有5人回復
354求調劑已經有5人回復
材料專業(yè)求調劑已經有3人回復
求材料調劑已經有7人回復
一志愿天大材料與化工（085600）總分338 已經有3人回復

1樓 2015-11-30 16:45:07

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

gczhusainan

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1089.9
帖子: 57
在線: 14.7小時
蟲號: 2429796
注冊: 2013-04-22
專業(yè): 巖土與基礎工程

我現(xiàn)在也遇到這個問題了，樓主咋解決的啊

贊一下

回復此樓

2樓2017-03-31 10:10:53