貝塞爾函數(shù)3-文檔資料

1、1貝塞爾函數(shù)貝塞爾函數(shù)32貝塞爾函數(shù)第二次課內(nèi)容總結(jié)貝塞爾函數(shù)第二次課內(nèi)容總結(jié)貝塞爾函數(shù)的遞推公式函數(shù)展成貝塞爾函數(shù)的級數(shù)貝塞爾函數(shù)應(yīng)用舉例貝塞爾函數(shù)應(yīng)用舉例3 1(1)nnnndx Jxx Jxdx 1(2)nnnndxJxxJxdx 上面兩式左邊的導(dǎo)數(shù)求出來上面兩式左邊的導(dǎo)數(shù)求出來, , 并經(jīng)過化簡并經(jīng)過化簡 1(3)nnnxJxnJxxJx 1(4)nnnxJxnJxxJx 4兩式相加減分別消去兩式相加減分別消去 和和 nJx nJx 11(25)nnnnJxJxJxx 112 (6)nnnJxJxJx5 1221221sin1mmmmdxJxxx dxx 1221221cosmmmdx

2、Jxxx dxx 這里微分算子這里微分算子 表示算子表示算子 連續(xù)連續(xù)作用作用 m 次的縮寫次的縮寫. . 1mdx dx1 dx dx半奇數(shù)級貝塞爾函數(shù)的表達(dá)式半奇數(shù)級貝塞爾函數(shù)的表達(dá)式可見，半奇數(shù)階的貝塞爾函數(shù)都是初等函數(shù)。6方程的通解為方程的通解為 nnP rCJrDYr 令令 ,xxry xP 方程轉(zhuǎn)化為方程轉(zhuǎn)化為 222220d ydyxxxny xdxdx ,nP rCJr 從而從而由于由于 , , 由條件由條件 知知 , |(0)|,P 0D (0)nY 7由由 可得可得: : 0P R 0.nJR 求特征問題求特征問題 2220r Prr PrrnP r 0,0.P RP 因此

3、，必須判明因此，必須判明 的零點(diǎn)的零點(diǎn)是否存在；如果存是否存在；如果存在，則需要研究其分布情形。在，則需要研究其分布情形。 nJx0rR貝塞爾函數(shù)的零點(diǎn)貝塞爾函數(shù)的零點(diǎn)8關(guān)于貝塞爾函數(shù)零點(diǎn)的結(jié)論：關(guān)于貝塞爾函數(shù)零點(diǎn)的結(jié)論：l 有無窮多個(gè)單重實(shí)零點(diǎn)有無窮多個(gè)單重實(shí)零點(diǎn), , 這些零點(diǎn)在這些零點(diǎn)在 x軸上關(guān)于原點(diǎn)對稱分布軸上關(guān)于原點(diǎn)對稱分布, , 因而因而 有無窮多個(gè)有無窮多個(gè)正的零點(diǎn)；正的零點(diǎn)； ( )nJx( )nJx24681012o1.00.5-0.50( )Jx1( )Jx9l 的零點(diǎn)和的零點(diǎn)和 的零點(diǎn)是彼此相間分的零點(diǎn)是彼此相間分布，布，即即 的任意兩個(gè)相鄰零點(diǎn)之間有且僅有的任意兩個(gè)相

4、鄰零點(diǎn)之間有且僅有一個(gè)一個(gè) 的零點(diǎn)，反之亦然；的零點(diǎn)，反之亦然； 1( )nJx 24681012o1.00.5-0.5( )nJx( )nJx1( )nJx 0( )Jx1( )Jx1024681012o1.00.5-0.5l以以 表示表示 的非負(fù)零點(diǎn)的非負(fù)零點(diǎn), , 則則 ( ) (1,2,)nmm( )( )1lim.nnmmm0( )Jx1( )Jx( )nJx函數(shù)以函數(shù)以 為周期振蕩為周期振蕩11與這些特征值相應(yīng)的與這些特征值相應(yīng)的特征函數(shù)特征函數(shù)為為 1,2,nmmnPrJrRm 方程方程 的解為：的解為： 0nJR ,1,2,nmRm 即貝塞爾方程相應(yīng)定解問題的即貝塞爾方程相應(yīng)定

5、解問題的特征值特征值為為 )2(,1,2,nnmmmR 12貝塞爾函數(shù)的正交性貝塞爾函數(shù)的正交性結(jié)論：結(jié)論：n 階貝塞爾特征函數(shù)系階貝塞爾特征函數(shù)系 ,1,2,nmmnPrJrmR 在區(qū)間在區(qū)間 (0, R) 上上帶權(quán)帶權(quán) r 正交正交， 2nmnJrR 模值模值的平方的平方 0dnnRmknnJrr JrrRR 2222110,22，nnnmnmmkRRJJmk 即即13結(jié)論結(jié)論2：在區(qū)間在區(qū)間0, R上具有一階連續(xù)導(dǎo)數(shù)以及上具有一階連續(xù)導(dǎo)數(shù)以及分段連續(xù)的二階導(dǎo)數(shù)的函數(shù)分段連續(xù)的二階導(dǎo)數(shù)的函數(shù) f (r) 如果在如果在 r=0 處處有界有界, , 在在 r =R 處等于零處等于零, , 則它

6、必可以展開為如則它必可以展開為如下形式的絕對且一致收斂的級數(shù)：下形式的絕對且一致收斂的級數(shù)： 1nmmnmf rA JrR 其中其中 22112mnnmARJ 0dRnmnrf r JrrR 模的平方模的平方權(quán)權(quán)函數(shù)函數(shù)貝塞爾函數(shù)貝塞爾函數(shù)14例例 設(shè)有半徑為設(shè)有半徑為 1 的薄均勻圓盤，其側(cè)面絕緣，的薄均勻圓盤，其側(cè)面絕緣，邊界上的溫度始終保持為零度，初始圓盤內(nèi)溫度邊界上的溫度始終保持為零度，初始圓盤內(nèi)溫度分布為分布為 其中其中 r 為圓盤內(nèi)任一點(diǎn)的極半徑，為圓盤內(nèi)任一點(diǎn)的極半徑，求圓盤的溫度分布規(guī)律。求圓盤的溫度分布規(guī)律。21,r 分析分析: : 由于是在圓域內(nèi)求解問題由于是在圓域內(nèi)求解問

7、題, , 故采用極坐故采用極坐標(biāo)標(biāo). . 考慮到定解條件和考慮到定解條件和 無關(guān)無關(guān), , 所以溫度所以溫度 u 只只能是能是 t 和和 r 的函數(shù)的函數(shù). . 15解解 根據(jù)問題的要求根據(jù)問題的要求, , 即可歸結(jié)為求下列方程即可歸結(jié)為求下列方程的定解問題：的定解問題：22222211,01uuuuartrrrr 由于由于 u 和和 無關(guān)無關(guān), , 可以化簡為問題可以化簡為問題 0,u 2221, 01,0,uuuarttrrr 201,01,turr 10,0.rut 16此外此外, , 由問題的物理意義由問題的物理意義, , 還有條件還有條件 ,u 且且 時(shí)，時(shí)， t 0.u 令令 ,u

8、 r tF r T t 代入到上述方程代入到上述方程, , 有有21,FFTrFa T 由此得由此得20,(1)TaT 220,(2)r FrFr F 解解(1)(1)得：得： 2,atT tCe 時(shí)，時(shí)，t 0,u 0 17(2) 為零階非標(biāo)準(zhǔn)的貝塞爾方程為零階非標(biāo)準(zhǔn)的貝塞爾方程, ,由由 u(r, t) 的有界性的有界性, 可以知道可以知道20.C 再由條件再由條件10,ru 知：知： 00,J 即即 是是 的零點(diǎn)的零點(diǎn). . 0( )Jx用用 ( (n =1,2) 表示表示 的正零點(diǎn)的正零點(diǎn), , 綜合綜合以上結(jié)果可得以上結(jié)果可得: :n 0( )Jx方程方程 (1) 的解為的解為 22

9、.atT tCe 2, 令令 則則 102 0F rC JrC Yr它的通解為：它的通解為：18220r FrFr F 方程方程 的特征值為：的特征值為： 2:1,2,nn 相應(yīng)的特征函數(shù)為：相應(yīng)的特征函數(shù)為： 0,nnFrJr 這時(shí)方程這時(shí)方程 的解為：的解為：20TaT 22natnnTtC e 從而從而 220,natnnnur tC eJr 19由疊加原理由疊加原理, , 可得原問題的解為可得原問題的解為 2201,natnnnu r tC eJr 由初始條件得由初始條件得 2011nnnrC Jr 其中其中 12200112nnnrrJrCdrJ 20 1130022001122nn

10、nnnCrJr drr Jr drJJ 因?yàn)橐驗(yàn)?1nnnnd x Jt xtx Jt x dx 令令 1,n 10nnnrJrdrJr dr 所以所以 1111000d.nnnnnrJrJrJrr 21 11132000nnnrJrr Jr drr d 131121002nnnnr Jrr Jr dr 1212202()nnnnJr Jr 1222()nnnnJJ 1nnnnd x Jt xtx Jt x dx 22從而從而 22214nnnnJCJ 所求定解問題的解為：所求定解問題的解為： 222022114( , )natnnnnnJu r tJr eJ 其中其中 是是 的正零點(diǎn)。的正零

11、點(diǎn)。 n 0( )Jr23例例hbU（圓柱形域內(nèi)的電勢分布）由導(dǎo)體壁構(gòu)成的空圓柱，其高為，半徑為，設(shè)上底的電勢為，側(cè)面和下底的電勢為零，試求圓柱體內(nèi)部的電勢分布。 oxyz xo 24解解由于區(qū)域?yàn)閳A柱形，所以采用柱坐標(biāo)系。2222211()0, 02, 0uuuRzz 一、建立方程一、建立方程2222222110uuuuz252222010 (0, 02 , 0) (3)0 ,(4)0(5)zz hbuuubzhzuuUu “翻譯翻譯”邊界邊界條件條件一、建立方程一、建立方程uz由于邊界條件與角度無關(guān)，因此所求的電勢只是、兩個(gè)自變量的函數(shù)，于是歸結(jié)為求下列定解問題：U為常數(shù)，為上底的電勢。2

12、6( , )( ) ( )uzRZ z應(yīng)用分離變量法，令，代入方程得1RRzRz 由此得0 (6)zz220 (7)RRR222( )( )() ( )0 0r P rrP rrnP rrR( )()()nnP rAJrBYr的通解為一、建立方程一、建立方程我們知道我們知道277方程（ ）為零階貝塞爾方程，其通解為00( )()()RCJDYu 由問題的物理意義知，由此推出(0)(8)R 一、建立方程一、建立方程5再由邊界條件（ ）得( )0 (9)R b 28于是，構(gòu)成具體問題的新的方程組為220 (7)RRR(0)(8)R ( )0 (9)R b NoImage一、建立方程一、建立方程29

13、7方程（ ）為零階貝塞爾方程，其通解為00( )()()RCJDY80 .D 由條件（ ）知二、求本征值、本征函數(shù)二、求本征值、本征函數(shù)3000(0)0()0( )( )mJbbJxJx即，由此可知是的零點(diǎn)。以表示的正零點(diǎn)，有(0)0()0mJ789從而，得到方程（ ）在條件（ ）、（ ）下的本征值以及本征函數(shù)(0)()mmb(0)0( )()mmRJb(1,2)m 二、求本征值、本征函數(shù)二、求本征值、本征函數(shù)9,再由條件（ ）得0( )()0R bCJb31三、由疊加原理三、由疊加原理 寫出解。寫出解。0 (6)zz6將上述的值代入方程（ ），可得(0)(0)( )mmzzeemmmZzCD

14、從而(0)(0)(0)0( , )()()mmzzmeemmmuzCDJb35由疊加原理，方程（ ）滿足（ ）的解為(0)(0)(0)01( , )()() (10)mmzzmeemmmuzCDJb32四、確定常數(shù)四、確定常數(shù)(0)(0)(0)01( , )()() (10)mmzzmeemmmuzCDJb00 ,(4)zz huuU4由條件( ），得(0)01( ,0)()()0mmmmuCDJb于是得0 (m1,2,) (11)mmCD5再由條件（ ）得(0)(0)(0)01( , )()()mmhhmeemmmuhCDJUb0(5)bu33(0)(0)(0)00(0)(0)(0)10 1()2(12)()()2mmbmhheemmmmmUJdbUCDbJJ四、確定常數(shù)四、確定常數(shù)依據(jù)貝塞爾級數(shù)展開系數(shù)公式，得)(2)()()(12)(0nknnknRkJRrdrRJrfrA 340 (m1,2,) (11)mmCD5.115.12將（）與先前得到的（）聯(lián)立，解得(0)(0)(0)1()mmmmUCshhJb(0)(0)(0)1()mmmmUDshhJb 將上述結(jié)果代入（10），得到原問題的解(0)(0)(0)01( , )()()m