 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1210 | 回復(fù): 36 | ||
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)該是怎么樣的呢？十分感謝！@月只藍(lán)@beefly |
回復(fù)此樓» 猜你喜歡
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有151人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有2人回復(fù)
 |
2樓2017-08-17 17:25:05
 |
3樓2017-08-17 19:47:06
|
|
4樓2017-08-17 22:59:40
|
|
5樓2017-08-18 01:07:01
|
|
6樓2017-08-18 02:06:45
 |
7樓2017-08-18 13:54:45
 |
8樓2017-08-18 15:56:36
|
|
9樓2017-08-18 16:43:44
|
|
10樓2017-08-19 00:30:27
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wgdd 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +12 | junqihahaha 2026-03-26 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +11 | 芋泥小鈴鐺 2026-03-28 | 11/550 |
|
|
|
[考研] 生物技術(shù)與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +12 | 李李不服輸 2026-03-26 | 12/600 |
|
|
|
[考研] 375求調(diào)劑 +6 | 雨夏整夜 2026-03-29 | 6/300 |
|
|
|
[考研] 309求調(diào)劑 +10 | 誰不是少年 2026-03-29 | 10/500 |
|
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +7 | 26要上岸 2026-03-28 | 7/350 |
|
|
|
[考研] 0856材料化工調(diào)劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
|
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[有機(jī)交流]
 有機(jī)合成求助
 20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|