版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

mavisym

金蟲 (小有名氣)

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

[求助] 求助，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 轉(zhuǎn)換為 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))+...

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

回復(fù)此樓

» 猜你喜歡

311求調(diào)劑已經(jīng)有4人回復(fù)
304求調(diào)劑已經(jīng)有4人回復(fù)
354求調(diào)劑已經(jīng)有7人回復(fù)
321求調(diào)劑已經(jīng)有3人回復(fù)
274求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有15人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有8人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有10人回復(fù)
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑已經(jīng)有8人回復(fù)
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有10人回復(fù)

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

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

gczhusainan

新蟲 (小有名氣)

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

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

贊一下

回復(fù)此樓

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