 |
 |
 |
| 24小時熱門版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1563 | 回復(fù): 34 | ||
wgdd木蟲 (正式寫手)
|
[求助]
求助一個數(shù)值積分問題，用matlab的quadgk函數(shù)來計算，謝謝！
|
|
現(xiàn)求解一個半無限震蕩積分問題： 被積函數(shù)為s*besselj(0,s*R)./(s.^2-k0^2), 積分區(qū)間為[0,inf], 其中R=6, k0=10, 均為已知常數(shù)。我打算用兩種方法來計算這個問題。(解析解和數(shù)值解) 方法一，解析解 對于這樣一個積分，可以查到數(shù)學(xué)手冊中積分結(jié)果為pi*i/2*besselh(0,1,k0*R). 經(jīng)過計算可以得到結(jié)果為-0.0744。 計算程序 % analytic method k0=10; R=6; p_inc=real(pi*1i/2*besselh(0,1,k0*R)); 方法二，采用matlab中的quadgk函數(shù)進行計算，由于k0為奇點，所以將積分區(qū)間分為兩塊。[0,k0]和[k0,inf]，計算結(jié)果為0.0390 計算程序 % numerical method k0=10; R=6; p_f=@(s)(s).*besselj(0,s*R)./(s.^2-k0^2); p1=quadgk(p_f,0,k0); p2=quadgk(p_f,k0,inf); p_incN=p1+p2; 求各位大牛支招，為啥兩種方法差別這么大呀~正確的數(shù)值積分方法應(yīng)該是怎么樣的呢？十分感謝！@wuming524 |
回復(fù)此樓» 猜你喜歡
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復(fù)
-
290求調(diào)劑
已經(jīng)有3人回復(fù)
-
327求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600 286分 材料求調(diào)劑
已經(jīng)有11人回復(fù)
 |
2樓2017-08-17 17:25:11
 |
3樓2017-08-17 19:47:18
 |
4樓2017-08-18 01:07:09
 |
5樓2017-08-18 02:06:54
 |
6樓2017-08-18 13:54:55
 |
7樓2017-08-18 15:56:44
|
|
8樓2017-08-18 16:43:55
 |
9樓2017-08-18 22:33:48
|
|
10樓2017-08-19 00:41:39
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wgdd 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工328分調(diào)劑 +7 | 。，。，。，。i 2026-03-23 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 337求調(diào)劑 +4 | 《樹》 2026-03-29 | 4/200 |
|
|
[考研] 277跪求調(diào)劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +4 | 篤行致遠(yuǎn). 2026-03-28 | 4/200 |
|
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個楠 2026-03-25 | 4/200 |
|
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研]
 材料考研調(diào)劑生
+3
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|