求解bessel方程

1、1一類一類BesselBessel方程在不同邊界條件下的差分解方程在不同邊界條件下的差分解信息與計算機科學專業(yè)趙鑫指導老師：李順初教授【摘要】本文主要闡述的是使用差分法解一類Bessel方程在不同邊界條件下的離散解，以及使用追趕法求解得到的線性方程組.在得到方程的離散解后，利用MATLAB程序，繪制方程的離散解，并分析這些離散解構成的解曲線，及各種參數的變化對解曲線的影響。關鍵字：微分方程Bessel方程差分法MATLAB程序Theki

2、ndofBesselequationdissolveindifferentboundary【Abstract】ThisthesisintroducessolvingthekindofBesselequationindifferentboundarybydifferencemethodgainsdiscretesolutionsolvingthelinearequationbychasemethod.Protractingthediscr

3、etesolutionoftheBesselequationbyMatlabobservingthecurvewhichiscomposeofthediscretesolutiondiscussingtheconditionofthecurveindifferentparameter.Keywd:BesselfunctionMatlabprogramDifferencemethod一前言前言Bessel函數是工程技術中常用的一種特殊函數

4、，主要地來源于圓柱形或圓柱物體的有關物理問題，因此也叫做圓柱函數。這種函數是研究一些物理問題時所歸結成的貝塞爾方程的解。????????0222???????xyxxyxxyx?早在貝塞爾之前，歐拉（Euler）、拉格朗日（Lagrange）等人研究過這個方程，并求出了它的通解。貝塞爾在1824年關于天文學問題的研究中，系統地研究了此方程的解，同時編制出這些函數具有十位小數的函數表。此后，稱這種函數為貝塞爾函數，并被廣泛地應用到物理學和

5、技術科學中。本文主要闡述的是使用差分法解一類Bessel方程在不同邊界條件下的定解問題的離散解，并借助MATLAB程序的繪圖工具對其給予直觀的表現，分析在不同邊界條件下參數的變化對解曲線的影響。二函數模型函數模型及其求解其求解2.1模型的建立本文要求解的模型如下：??????????????????????????????????????.1000)(01122RuRbaRyRyuQdxdybxyayudxdyxdxydxx?均為常數；

6、且其中：右邊界條件左邊界條件3右邊界有兩種情況下面將分別討論.2.2.1第一種右邊界條件0)(?Ry因為所以Rnh??1…………(6)0)1(??nhy故由(4)(5)(6)可以得到一個階線性方程組:n??????????????????????????????????????????????????????????????????????????????????????????????0))1(1()])2(1(2[))2(1(]))

7、2(1)(1[(.0))1(1())3(1())2(1()])3(1(2[))3(1(]))3(1)(1[(.........................................0)51()31()41()52()31(])31)(1[(0)41()21()31()32()21(])21)(1[(0)31()1()21()2()1(])1)(1[(0)21()1()2()1(])1[()1()1()(222222hnfhnh

8、hnfhhnuhhnfhnhnfhnhhnfhhnuhhfhhfhhfhhuhhfhhfhhfhhuhhfhhfhhfhhuhhfhfhfhuhqhhbffbah…………………………………………………………………………………………(7)2.2.2第二種右邊界條件0)(?Ry根據導數的定義得:0)()(|)()(|??????????hRfhRfxxfhxfdxdyRxRx即:)()(RfhRf??也即:…………………(8))1())1(

9、1(nhfhnf????故由(4)(5)(8)可以得到一個n階線性方程組:??????????????????????????????????????????????????????????????????????????????????????????????0))1(1()])2(1([))2(1(]))2(1)(1[(.0))1(1())3(1())2(1()])3(1(2[))3(1(]))3(1)(1[(...........

10、..............................0)51()31()41()52()31(])31)(1[(0)41()21()31()32()21(])21)(1[(0)31()1()21()2()1(])1)(1[(0)21()1()2()1(])1[()1()1()(222222hnfhnhhnfhhnuhhnfhnhnfhnhhnfhhnuhhfhhfhhfhhuhhfhhfhhfhhuhhfhhfhhfhhuhhf