 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 2450 | 回復(fù): 1 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
mavisym金蟲(chóng) (小有名氣)
|
[求助]
求助，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è)公式本來(lái)是一個(gè)長(zhǎng)的，分開(kāi)找錯(cuò)誤的，如果是個(gè)長(zhǎng)的就報(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ù)此樓» 猜你喜歡
一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑
已經(jīng)有3人回復(fù)
-
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
【同濟(jì)軟件】軟件（085405）考研求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有4人回復(fù)
-
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿武理材料305分求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專(zhuān)碩英一數(shù)二306
已經(jīng)有3人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有7人回復(fù)
gczhusainan
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1089.9
- 帖子: 57
- 在線(xiàn): 14.7小時(shí)
- 蟲(chóng)號(hào): 2429796
- 注冊(cè): 2013-04-22
- 專(zhuān)業(yè): 巖土與基礎(chǔ)工程
2樓2017-03-31 10:10:53
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 mavisym 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線(xiàn) +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 材料工程專(zhuān)碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研]
 機(jī)械專(zhuān)碩調(diào)劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|
|
|
[考研] 中科大材料專(zhuān)碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 328化工專(zhuān)碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|