| 5 | 1/1 | 返回列表 |
| 查看: 1657 | 回復(fù): 4 | |||
| 本帖產(chǎn)生 1 個(gè) 數(shù)學(xué)EPI ,點(diǎn)擊這里進(jìn)行查看 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
zk19900227木蟲(chóng) (著名寫(xiě)手)
|
[求助]
求大神來(lái)救,關(guān)于貝塞爾函數(shù)的理論和編程的問(wèn)題。。。。 已有1人參與
|
||
|
J = besselj(nu,Z,1) computes besselj(nu,Z).*exp(-abs(imag(Z))). 以上是matlab中第一類貝塞爾函數(shù)中的一個(gè),特點(diǎn)就是將貝塞爾函數(shù)中的指數(shù)項(xiàng)提取出來(lái),其理論和編程方法是什么? 最好講一講 啊,或是提供一份明確便于學(xué)習(xí)的資料。。 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 。。。。。。。。。。。以下是matlab中bessel函數(shù)的鏈接。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 http://cn.mathworks.com/help/matlab/ref/besselj.html |

木蟲(chóng) (著名寫(xiě)手)

榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
鐵桿木蟲(chóng) (著名寫(xiě)手)
工程和科學(xué)數(shù)值計(jì)算咨詢
|
猜測(cè)lz意圖解釋一下。 興趣點(diǎn)似乎在 besselj(nu,Z,1),即把 besselj(nu,Z)指數(shù)增長(zhǎng)系數(shù)消掉后所剩部分的有關(guān)原理和背景。一般第一類Bessel函數(shù)的使用都限制在實(shí)數(shù)范圍內(nèi),在z很大時(shí)漸進(jìn)展開(kāi)首項(xiàng)是 sqrt(2/(pi z))*cos(z-nu *pi/4-pi/4) ,這時(shí)besselj(nu,Z,1)和 besselj(nu,Z) 是一樣的。但是如果 z變量是復(fù)數(shù),那么besselj(nu,Z) 有一個(gè) exp (|imag(z)|)/z^(3/2) 的附加項(xiàng),這是指數(shù)增長(zhǎng)的。besselj(nu,Z,1) = besselj(nu,Z) * exp(-abs(imag(Z))) 的意思就是把 這個(gè)指數(shù)系數(shù)去掉。這個(gè)在估計(jì)指數(shù)增長(zhǎng)量級(jí)或漸進(jìn)展開(kāi)是很有用途的,否則浮點(diǎn)計(jì)算很快會(huì)上溢,就什么都看不到了。關(guān)于besselj(nu,Z,1)的編程我想可以根據(jù) 以下兩點(diǎn)就可以做了:(1) J(nu,z)和modified Bessel function I(nu,z)的關(guān)系 (2) I(nu,z) 的帶exp(z)的漸進(jìn)展開(kāi)。細(xì)節(jié)繁復(fù),這里就不具體展開(kāi)了。 其實(shí)2樓給出的鏈接 http://en.wikipedia.org/wiki/Bes ... _I.CE.B1.2C_K.CE.B1 就是一份最好的學(xué)習(xí)資料,已經(jīng)能解釋大部分問(wèn)題了;癁閷(shí)用當(dāng)然還需要一些數(shù)值經(jīng)驗(yàn)。 |
木蟲(chóng) (著名寫(xiě)手)
|
謝謝回答!我是做計(jì)算電磁學(xué)的,最近用fortran做了一個(gè)三層介質(zhì)磁偶極子源激發(fā)的電磁場(chǎng)的程序。。。里面有對(duì)貝塞爾函數(shù)J0(實(shí)數(shù))*exp(ix)樣式的公式的積分,出現(xiàn)了算不準(zhǔn)的幾種情況,得到的積分值的虛部過(guò)大。想解決這個(gè)問(wèn)題,看看能不能從貝塞爾函數(shù)上入手,但是看你的回答,實(shí)數(shù)范圍內(nèi)兩種貝塞爾函數(shù)是一樣的,這條路估計(jì)走不通了。。。。我用的貝塞爾函數(shù)是copy的,是使用插值法得到的。。。。 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 0856,材料與化工321分求調(diào)劑 +7 | 大饞小子 2026-03-27 | 8/400 |
|
|
[考研] 291求調(diào)劑 +6 | 孅華 2026-03-22 | 6/300 |
|
|
[考研] 286求調(diào)劑 +8 | PolarBear11 2026-03-26 | 8/400 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 材料求調(diào)劑 +8 | @taotao 2026-03-21 | 8/400 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研] 349求調(diào)劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 考研調(diào)劑 +6 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-21 | 7/350 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研] 化學(xué)調(diào)劑 +6 | yzysaa 2026-03-21 | 6/300 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問(wèn)津 2026-03-23 | 5/250 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|