贝塞尔方程与贝塞尔函数学习笔记
《數(shù)學(xué)物理方法(顧樵)》第13章學(xué)習(xí)筆記
第一節(jié) 幾個(gè)微分方程的引入
- 三維波動(dòng)方程:
?2v?t2=a2(?2v?x2+?2v?y2+?2v?z2)≡a2?2v\frac{\partial^2 v}{\partial t^2} = a^2 (\frac{\partial^2 v}{\partial x^2} + \frac{\partial^2 v}{\partial y^2} + \frac{\partial^2 v}{\partial z^2}) \equiv a^2 \nabla^2 v?t2?2v?=a2(?x2?2v?+?y2?2v?+?z2?2v?)≡a2?2v - 三維熱傳導(dǎo)方程:
?v?t=a2(?2v?x2+?2v?y2+?2v?z2)≡a2?2v\frac{\partial v}{\partial t} = a^2 (\frac{\partial^2 v}{\partial x^2} + \frac{\partial^2 v}{\partial y^2} + \frac{\partial^2 v}{\partial z^2}) \equiv a^2 \nabla^2 v?t?v?=a2(?x2?2v?+?y2?2v?+?z2?2v?)≡a2?2v
對三維波動(dòng)方程與三維熱傳導(dǎo)方程使用分離變量法,得到時(shí)間上的方程,以及空間上的名為亥姆霍茲方程的方程。
- 亥姆霍茲方程:
?2u(r)+k2u(r)=0\nabla^2 u(\pmb r) + k^2u(\pmb r) = 0?2u(rrr)+k2u(rrr)=0
將亥姆霍茲方程變換到球坐標(biāo)上,再次應(yīng)用分離變量法,得到以半徑為自變量的球貝塞爾方程,以及以半徑與 zzz 軸夾角為自變量再經(jīng)變量代換得到的連帶勒讓德方程。
球貝塞爾方程中設(shè)定特殊值,可以得到歐拉方程。
連帶勒讓德方程中設(shè)定特殊值,可以得到勒讓德方程。
- 球貝塞爾方程:
ddr(r2dRdr)+(k2r2?w2)R=0\fracozvdkddzhkzd{dr}(r^2 \frac{dR}{dr})+(k^2r^2-w^2)R = 0drd?(r2drdR?)+(k2r2?w2)R=0 - 連帶勒讓德方程:
ddr[(1?x2)dydx]+(w2?m21?x2)y=0\fracozvdkddzhkzd{dr}[(1-x^2)\frac{dy}{dx}]+(w^2-\frac{m^2}{1-x^2})y = 0drd?[(1?x2)dxdy?]+(w2?1?x2m2?)y=0 - 歐拉方程:
ddr(r2dRdr)?w2R=0\fracozvdkddzhkzd{dr}(r^2 \frac{dR}{dr})-w^2R = 0drd?(r2drdR?)?w2R=0 - 勒讓德方程:
ddr[(1?x2)dydx]+w2y=0\fracozvdkddzhkzd{dr}[(1-x^2)\frac{dy}{dx}]+w^2y = 0drd?[(1?x2)dxdy?]+w2y=0
將亥姆霍茲方程變換到柱坐標(biāo)上,再次應(yīng)用分離變量法,得到以半徑為自變量的方程,進(jìn)一步應(yīng)用變量代換,得到貝塞爾方程。
- 貝塞爾方程:
x2d2ydx2+xdydx+(x2?m2)y=0x^2\frac{d^2 y}{dx^2} + x\frac{dy}{dx} + (x^2 -m^2)y = 0x2dx2d2y?+xdxdy?+(x2?m2)y=0
這些方程也可以直接通過施圖姆-劉維爾型方程引入。
ddx[k(x)dydx]?q(x)y+λρ(x)y=0\fracozvdkddzhkzd{dx}[k(x)\frac{dy}{dx}]-q(x)y+\lambda \rho(x)y = 0dxd?[k(x)dxdy?]?q(x)y+λρ(x)y=0
所以對于這些函數(shù)的本征函數(shù)集,可以通過施圖姆-劉維爾型方程的結(jié)論驗(yàn)證正交性。
第二節(jié) 伽馬函數(shù)的基本知識
- 定義
Γ(x)=∫0∞e?ttx?1dt(x>0)\Gamma(x)=\int ^\infty _0 e^{-t}t^{x-1}dt\ \ \ \ (x>0)Γ(x)=∫0∞?e?ttx?1dt????(x>0) - 基本性質(zhì)
| Γ(1)=1\Gamma(1)=1Γ(1)=1 | Γ(12)=π\(zhòng)Gamma(\frac{1}{2})=\sqrt{\pi}Γ(21?)=π? |
| Γ(x+1)=xΓ(x)\Gamma(x+1)=x\Gamma(x)Γ(x+1)=xΓ(x) | Γ(n+1)=n!(n=0,1,2,...)\Gamma(n+1)=n!\ \ \ \ (n=0,1,2,...)Γ(n+1)=n!????(n=0,1,2,...) |
| Γ(n+12)=(2n)!22nn!π\(zhòng)Gamma(n+\frac{1}{2})=\frac{(2n)!}{2^{2n}n!}\sqrt{\pi}Γ(n+21?)=22nn!(2n)!?π? | Γ(n+12+1)=(2n+1)!22n+1n!π\(zhòng)Gamma(n+\frac{1}{2}+1)=\frac{(2n+1)!}{2^{2n+1}n!}\sqrt{\pi}Γ(n+21?+1)=22n+1n!(2n+1)!?π? |
當(dāng)n比較大的時(shí)候,使用變量代換,可得到斯特林公式:
n!≈2πnnne?nn!\approx \sqrt{2\pi n}n^n e^{-n}n!≈2πn?nne?n
第三節(jié) 求解貝塞爾方程
使用 Frobenius方法 得到級數(shù)形式的解的系數(shù)的方程,進(jìn)而得到第一類貝塞爾函數(shù)。
貝塞爾方程的通解有兩種形式。
在討論貝塞爾方程通解的第二種形式的時(shí)候,利用第一類貝塞爾方程構(gòu)造得到第二類vvv階貝塞爾函數(shù)(也稱 諾依曼函數(shù) )。
- vvv階第一類貝塞爾函數(shù):
Jv(x)=∑m=0∞(?1)mm!Γ(m+v+1)(x2)2m+vJ_v(x)=\sum^\infty _{m=0}\frac{(-1)^m}{m!\Gamma(m+v+1)}(\frac{x}{2})^{2m+v}Jv?(x)=m=0∑∞?m!Γ(m+v+1)(?1)m?(2x?)2m+v - 諾依曼函數(shù)
Yv(x)={Jv(x)cos?vπ?J?v(x)sin?vπv?Zlim?α→vJα(x)cos?απ?J?α(x)sin?απv∈ZY_v(x)=\begin{dcases} \frac{J_v(x)\cos v\pi-J_{-v}(x)}{\sin v\pi}\ \ \ \ \ \ \ \ \ \ \ \ v \notin \mathbb{Z} \\ \\ \lim_{\alpha \to v}\frac{J_\alpha(x)\cos \alpha\pi-J_{-\alpha}(x)}{\sin \alpha\pi}\ \ \ \ v \in \mathbb{Z}\\ \end{dcases} Yv?(x)=??????????sinvπJv?(x)cosvπ?J?v?(x)?????????????v∈/?Zα→vlim?sinαπJα?(x)cosαπ?J?α?(x)?????v∈Z?
兩個(gè)補(bǔ)充:
Yn(x)=2π(ln?x2+γ)Jn(x)?1π∑k=0n?1(n?l?1)!k!(x2)2k?n?1π∑k=0n?1(?1)kk!(n+k)2k+n\begin{aligned}Y_n(x)=&\frac{2}{\pi}(\ln\frac{x}{2}+\gamma)J_n(x)\\ &-\frac{1}{\pi}\sum^{n-1}_{k=0}\frac{(n-l-1)!}{k!}(\frac{x}{2})^{2k-n}\\&-\frac{1}{\pi}\sum^{n-1}_{k=0}\frac{(-1)^k}{k!(n+k)!}[\Phi(k)+\Phi(n+k)](\frac{x}{2})^{2k+n}\end{aligned}Yn?(x)=?π2?(ln2x?+γ)Jn?(x)?π1?k=0∑n?1?k!(n?l?1)!?(2x?)2k?n?π1?k=0∑n?1?k!(n+k)!(?1)k?[Φ(k)+Φ(n+k)](2x?)2k+n?
Yn(x)=1π∫0πsin?(xsin?θ?nθ)dθ?1π∫0∞[ent+(?1)ne?nt]e?xsinh?tdt\begin{aligned}Y_n(x)=&\frac{1}{\pi}\int ^{\pi} _{0} \sin (x\sin \theta-n\theta)d\theta \\ &-\frac{1}{\pi}\int ^{\infty} _{0}[e^{nt}+(-1)^n e^{-nt}]e^{-x\sinh t}dt\end{aligned}Yn?(x)=?π1?∫0π?sin(xsinθ?nθ)dθ?π1?∫0∞?[ent+(?1)ne?nt]e?xsinhtdt?
其中:
Φ(n)=1+12+13+...+1n,Φ(0)=0γ=lim?n→∞(Φ(n)?ln?n)=0.577\begin{aligned}\Phi(n)&=1+\frac{1}{2}+\frac{1}{3}+...+\frac{1}{n}, \ \ \ \ \Phi(0)=0\\ \gamma&=\lim _{n\to \infty} (\Phi(n)-\ln n)=0.577\end{aligned}Φ(n)γ?=1+21?+31?+...+n1?,????Φ(0)=0=n→∞lim?(Φ(n)?lnn)=0.577?
第四節(jié) 貝塞爾函數(shù)的基本性質(zhì)
- 生成函數(shù):該函數(shù)的級數(shù)展開式的系數(shù)是貝塞爾函數(shù)。
整數(shù)階貝塞爾函數(shù) Jn(x)J_n(x)Jn?(x) 的生成函數(shù):
exp?[x2(r?1r)]=∑n=?∞∞Jn(x)rn\exp [\frac{x}{2}(r-\frac{1}{r})]=\sum^{\infty}_{n=-\infty}J_n(x)r^nexp[2x?(r?r1?)]=n=?∞∑∞?Jn?(x)rn - 性質(zhì)
| ddx[xvJv(x)]=xvJv?1(x)\fracozvdkddzhkzd{dx}[x^vJ_v(x)]=x^vJ_{v-1}(x)dxd?[xvJv?(x)]=xvJv?1?(x) | ddx[xvYv(x)]=xvYv?1(x)\fracozvdkddzhkzd{dx}[x^vY_v(x)]=x^vY_{v-1}(x)dxd?[xvYv?(x)]=xvYv?1?(x) |
| ddx[x?vJv(x)]=?x?vJv+1(x)\fracozvdkddzhkzd{dx}[x^{-v}J_v(x)]=-x^{-v}J_{v+1}(x)dxd?[x?vJv?(x)]=?x?vJv+1?(x) | ddx[x?vYv(x)]=?x?vYv+1(x)\fracozvdkddzhkzd{dx}[x^{-v}Y_v(x)]=-x^{-v}Y_{v+1}(x)dxd?[x?vYv?(x)]=?x?vYv+1?(x) |
| Jv′(x)=12[Jv?1(x)?Jv+1(x)]J'_v(x)=\frac{1}{2}[J_{v-1}(x)-J_{v+1}(x)]Jv′?(x)=21?[Jv?1?(x)?Jv+1?(x)] | Jv′(x)=12[Jv?1(x)?Jv+1(x)]J'_v(x)=\frac{1}{2}[J_{v-1}(x)-J_{v+1}(x)]Jv′?(x)=21?[Jv?1?(x)?Jv+1?(x)] |
| Jv?1(x)+Jv+1(x)=2vxJv(x)J_{v-1}(x)+J_{v+1}(x)=\frac{2v}{x}J_v(x)Jv?1?(x)+Jv+1?(x)=x2v?Jv?(x) | Yv?1(x)+Yv+1(x)=2vxYv(x)Y_{v-1}(x)+Y_{v+1}(x)=\frac{2v}{x}Y_v(x)Yv?1?(x)+Yv+1?(x)=x2v?Yv?(x) |
| xJv?1(x)=vJv(x)+xJv′(x)xJ_{v-1}(x)=vJ_v(x)+xJ'_v(x)xJv?1?(x)=vJv?(x)+xJv′?(x) | xYv?1(x)=vYv(x)+xYv′(x)xY_{v-1}(x)=vY_v(x)+xY'_v(x)xYv?1?(x)=vYv?(x)+xYv′?(x) |
| xJv+1(x)=vJv(x)?xJv′(x)xJ_{v+1}(x)=vJ_v(x)-xJ'_v(x)xJv+1?(x)=vJv?(x)?xJv′?(x) | xYv+1(x)=vYv(x)?xYv′(x)xY_{v+1}(x)=vY_v(x)-xY'_v(x)xYv+1?(x)=vYv?(x)?xYv′?(x) |
- 整數(shù)階貝塞爾函數(shù)積分形式
有兩種方法得到其積分形式。
一是根據(jù)生成函數(shù)在復(fù)數(shù)域上的解析函數(shù),由其洛朗級數(shù)系數(shù)在特殊閉合回路上得到。
二是由同樣的解析函數(shù)出發(fā),在某個(gè)特殊閉合回路上將函數(shù)展開,通過比較等號左右兩邊的形式,結(jié)合三角函數(shù)的正交性,再通過三角函數(shù)公式得到積分形式。
Jn=1π∫0πcos?(xsin?θ?nθ)dθ(n=0,±1,±2,...)J_n=\frac{1}{\pi}\int ^\pi _0 \cos (x\sin \theta-n\theta)d\theta\ \ \ \ (n=0,\pm1, \pm2,...)Jn?=π1?∫0π?cos(xsinθ?nθ)dθ????(n=0,±1,±2,...) - 整數(shù)階貝塞爾函數(shù)漸進(jìn)公式
使用穩(wěn)定相方法獲取其漸進(jìn)公式:
Jn(x)≈2πxcos?(x?π4?nπ2)(n=0,1,2,...)J_n(x)\approx \sqrt{\frac{2}{\pi x}}\cos(x-\frac{\pi}{4}-\frac{n\pi}{2})\ \ \ \ (n=0,1,2,...)Jn?(x)≈πx2??cos(x?4π??2nπ?)????(n=0,1,2,...)
第五節(jié) 貝塞爾函數(shù)的正交完備性
∫0axJv(λvmx)Jv(λvkx)dx=0(m≠k)\int_0^a xJ_v(\lambda_{vm}x)J_v(\lambda_{vk}x)dx=0\ \ \ \ (m\ne k)∫0a?xJv?(λvm?x)Jv?(λvk?x)dx=0????(m?=k)
并計(jì)算模值:
∫0axJv2(λvmx)dx=a22Jv+12(μvm)(m=1,2,...)\int_0^a xJ^2_v(\lambda_{vm}x)dx=\frac{a^2}{2}J^2_{v+1}(\mu_{vm})\ \ \ \ (m=1,2,...)∫0a?xJv2?(λvm?x)dx=2a2?Jv+12?(μvm?)????(m=1,2,...)
f(x)=∑m=1∞AmJv(λvmx)f(x)=\sum^\infty _{m=1} A_mJ_v(\lambda_{vm}x)f(x)=m=1∑∞?Am?Jv?(λvm?x)
Am=2[aJv+1(μvm)]2∫0axJv(λvmx)f(x)dxA_m=\frac{2}{[aJ_{v+1}(\mu_{vm})]^2}\int^a_0 xJ_v(\lambda_{vm}x)f(x)dxAm?=[aJv+1?(μvm?)]22?∫0a?xJv?(λvm?x)f(x)dx
貝塞爾級數(shù)在間斷點(diǎn)處的收斂性由狄利克雷定理確定。
總結(jié)
以上是生活随笔為你收集整理的贝塞尔方程与贝塞尔函数学习笔记的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: correct ways to defi
- 下一篇: c语言标准库 swap,swap