 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1108 | 回復(fù): 1 | |||||
[求助]
求助：用matlab 循環(huán)計(jì)算橢圓積分 已有1人參與
|
|
大家好， 請(qǐng)教大家一個(gè)在循環(huán)中計(jì)算橢圓積分的問(wèn)題。我要計(jì)算的基本方程如附圖。方程有三個(gè)變量，y,f,v. 我想畫(huà)出y跟f之間的函數(shù)圖像。 我的code如下： syms fi;v = [0:0.001:pi/2]; L = length(v); y = zeros(1,L); k = zeros(1,L); for i = 1:L t1 = 1./sqrt(sin(v(i))-sin(fi)); m1 = int(t1,0,v(i)); k(i) = vpa(m1) t2 = sin(fi)/(sin(v(i))-sin(fi)); m2 = 1/k(i)*int(t1,0,v(i)); y(i) = vpa(m2) f(i) = k(i)^2./2; end plot(f,y,'b'); 系統(tǒng)提示錯(cuò)誤如下： The following error occurred converting from sym to double: Error using 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. Error in NonlinearStiffness (line 10) k(i) = vpa(m1) 但是我將積分范圍換成具體數(shù)值而不是數(shù)組元素時(shí)，是可以得到一個(gè)結(jié)果： >> t1 = 1./sqrt(sin(0.2)-sin(fi)); m1 = int(t1,0,0.2); >> m1 m1 = (7217745006463825^(1/2)*9007199254740992^(1/2)*ellipticF(pi/4, 18014398509481984/7217745006463825)*2*i)/7217745006463825 - (2*7217745006463825^(1/2)*ellipticF(pi/4 - 1/10, 18014398509481984/7217745006463825)*(9007199254740992*sin(1/5) - 1789454248277167)^(1/2))/(7217745006463825*(1789454248277167/9007199254740992 - sin(1/5))^(1/2)) >> vpa(m1) ans = 0.901 + 1.145e-36*i 我有兩個(gè)疑問(wèn)： 1. 為什么語(yǔ)句在循環(huán)中就運(yùn)行不通呢？ 2. 為什么積分結(jié)果會(huì)有虛部呢？ 橢圓積分.png |
回復(fù)此樓» 收錄本帖的淘帖專(zhuān)輯推薦
 matlab編程繪圖 |
» 猜你喜歡
0856調(diào)劑，是學(xué)校就去
已經(jīng)有9人回復(fù)
-
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有9人回復(fù)
-
265求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
已經(jīng)有7人回復(fù)
-
材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦
已經(jīng)有6人回復(fù)
-
295復(fù)試調(diào)劑
已經(jīng)有7人回復(fù)
-
本人考085602 化學(xué)工程 專(zhuān)碩
已經(jīng)有20人回復(fù)
-
一志愿中南化學(xué)337求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有23人回復(fù)
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專(zhuān)家經(jīng)驗(yàn): +1059 - 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲(chóng)號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專(zhuān)業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
iyearn: 金幣+10, ★★★很有幫助 2015-10-02 09:47:23
感謝參與，應(yīng)助指數(shù) +1
iyearn: 金幣+10, ★★★很有幫助 2015-10-02 09:47:23
|
k y f 1.0e+003 * 0 0 0 0.0004 0.0018 0.0001 0.0006 0.0036 0.0002 0.0008 0.0054 0.0003 0.0009 0.0074 0.0004 0.0010 0.0092 0.0005 0.0011 0.0112 0.0006 0.0012 0.0130 0.0007 0.0013 0.0153 0.0008 0.0014 0.0173 0.0010 0.0015 0.0197 0.0011 0.0016 0.0219 0.0012 0.0017 0.0243 0.0014 0.0018 0.0275 0.0015 0.0018 0.0302 0.0017 0.0019 0.0335 0.0019 0.0020 0.0371 0.0021 0.0021 0.0410 0.0023 0.0022 0.0447 0.0025 0.0024 0.0495 0.0028 0.0025 0.0557 0.0031 0.0026 0.0620 0.0034 0.0028 0.0704 0.0038 0.0029 0.0807 0.0043 0.0031 0.0916 0.0048 0.0033 0.1076 0.0055 0.0035 0.1290 0.0063 0.0038 0.1594 0.0074 0.0042 0.2067 0.0088 0.0047 0.2905 0.0110 0.0054 0.4838 0.0148 0.0072 1.4962 0.0257 附圖1.png |
2樓2015-08-21 11:41:24
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 iyearn 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿中南化學(xué)337求調(diào)劑 +4 | niko- 2026-03-19 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +8 | S.H1 2026-03-18 | 8/400 |
|
|
[考研] 材料080500調(diào)劑求收留 +6 | 一顆meteor 2026-03-13 | 6/300 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +8 | 步川酷紫123 2026-03-13 | 8/400 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|