 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1565 | 回復(fù): 34 | |||
wgdd木蟲 (正式寫手)
|
[求助]
求助一個(gè)數(shù)值積分問題，用matlab的quadgk函數(shù)來計(jì)算，謝謝！
|
|
現(xiàn)求解一個(gè)半無限震蕩積分問題： 被積函數(shù)為s*besselj(0,s*R)./(s.^2-k0^2), 積分區(qū)間為[0,inf], 其中R=6, k0=10, 均為已知常數(shù)。我打算用兩種方法來計(jì)算這個(gè)問題。(解析解和數(shù)值解) 方法一，解析解 對(duì)于這樣一個(gè)積分，可以查到數(shù)學(xué)手冊(cè)中積分結(jié)果為pi*i/2*besselh(0,1,k0*R). 經(jīng)過計(jì)算可以得到結(jié)果為-0.0744。 計(jì)算程序 % analytic method k0=10; R=6; p_inc=real(pi*1i/2*besselh(0,1,k0*R)); 方法二，采用matlab中的quadgk函數(shù)進(jìn)行計(jì)算，由于k0為奇點(diǎn)，所以將積分區(qū)間分為兩塊。[0,k0]和[k0,inf]，計(jì)算結(jié)果為0.0390 計(jì)算程序 % 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ù)此樓» 猜你喜歡
309求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(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ù)
 |
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ā)表 | |
|---|---|---|---|---|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +4 | 大饞小子 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +10 | leaflight 2026-03-24 | 10/500 |
|
|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +5 | 真名有冰 2026-03-29 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +7 | 小蠟新筆 2026-03-29 | 7/350 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +11 | 北道巷 2026-03-24 | 11/550 |
|
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
|
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑！�。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|