日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

贝塞尔方程与贝塞尔函数学习笔记

發(fā)布時(shí)間:2023/12/10 编程问答 42 豆豆
生活随笔 收集整理的這篇文章主要介紹了 贝塞尔方程与贝塞尔函数学习笔记 小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.

《數(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)?????vZ?

兩個(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)![Φ(k)+Φ(n+k)](x2)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=0n?1?k!(n?l?1)!?(2x?)2k?n?π1?k=0n?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=nlim?(Φ(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ì)
第一類貝塞爾函數(shù)第二類貝塞爾函數(shù)
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ù)的正交完備性

  • 類比 正弦函數(shù)集 Sm(x)=sin?(mπx)S_m(x)=\sin(m\pi x)Sm?(x)=sin(mπx),構(gòu)建正交的貝塞爾函數(shù)集。
  • 結(jié)合參數(shù)形式的貝塞爾函數(shù)集 Jv(λx)J_v(\lambda x)Jv?(λx) 的性質(zhì),在區(qū)間 [0,a][0,a][0,a] 上證明 正交性
    ∫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)容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。