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

登錄注冊

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

24小時(shí)熱門版塊排行榜

返回列表

mavisym

金蟲 (小有名氣)

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

[求助] 求助，MATLAB報(bào)錯(cuò)

我的公式明明除的是實(shí)數(shù)，matlab報(bào)錯(cuò)如下：
錯(cuò)誤使用 mupadmex
Error in MuPAD command: Division by zero.
[_power]
  Evaluating: symobj::trysubs

出錯(cuò) sym/subs>mupadsubs (line 139)
G = mupadmex('symobj::fullsubs',F.s,X2,Y2);

出錯(cuò) sym/subs (line 124)
G = mupadsubs(F,X,Y);

出錯(cuò) 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

我注釋的三個(gè)公式本來是一個(gè)長的，分開找錯(cuò)誤的，如果是個(gè)長的就報(bào)別的錯(cuò)
從 sym 轉(zhuǎn)換為 double 時(shí)出現(xiàn)以下錯(cuò)誤:
錯(cuò)誤使用 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.

出錯(cuò) 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))+...

求大神幫忙，實(shí)在不知道怎么搞

回復(fù)此樓

» 猜你喜歡

274求調(diào)劑已經(jīng)有4人回復(fù)
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有10人回復(fù)
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有7人回復(fù)
材料專業(yè)求調(diào)劑已經(jīng)有5人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有10人回復(fù)
344求調(diào)劑已經(jīng)有7人回復(fù)
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑已經(jīng)有7人回復(fù)
271材料工程求調(diào)劑已經(jīng)有4人回復(fù)
材料專碩英一數(shù)二306 已經(jīng)有4人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有4人回復(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小時(shí)
蟲號(hào): 2429796
注冊: 2013-04-22
專業(yè): 巖土與基礎(chǔ)工程

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

贊一下

回復(fù)此樓

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