一类Bessel方程在不同边界条件下的差分解.docVIP

一类Bessel方程在不同边界条件下的差分解 信息与计算机科学专业 赵鑫 指导老师：李顺初 教授 【摘要】 本文主要阐述的是使用差分法解一类Bessel方程在不同边界条件下的离散解，以及使用追赶法求解得到的线性方程组.在得到方程的离散解后，利用MATLAB程序，绘制方程的离散解，并分析这些离散解构成的解曲线，及各种参数的变化对解曲线的影响。 关键字：微分方程 Bessel方程 差分法 MATLAB程序 The kind of Bessel equation dissolve in different boundary 【Abstract】 This thesis introduces solving the kind of Bessel equation in different boundary by difference method and gains discrete solution, and solving the linear equation by chase method . Protracting the discrete solution of the Bessel equation by Matlab, and observing the curve which is compose of the discrete solution, discussing the condition of the curve in different parameter. Key word: Bessel function Matlab program Difference method 一 前言 Bessel函数是工程技术中常用的一种特殊函数，主要地来源于圆柱形或圆柱物体的有关物理问题，因此也叫做圆柱函数。这种函数是研究一些物理问题时所归结成的贝塞尔方程的解。 早在贝塞尔之前，欧拉（Euler）、拉格朗日（Lagrange）等人研究过这个方程，并求出了它的通解。贝塞尔在1824年关于天文学问题的研究中，系统地研究了此方程的解，同时编制出这些函数具有十位小数的函数表。此后，称这种函数为贝塞尔函数，并被广泛地应用到物理学和技术科学中。 本文主要阐述的是使用差分法解一类Bessel方程在不同边界条件下的定解问题的离散解，并借助MATLAB程序的绘图工具对其给予直观的表现，分析在不同边界条件下参数的变化对解曲线的影响。 二 函数模型及其求解 2.1 模型的建立 本文要求解的模型如下： 2.2 模型的分析与求解 采用差分法求解上面的模型，首先应该利用导数的定义将导数化为： (1) (2) 先写出差分格式。将（1）（2）代入微分方程得到： (3) 在区间[1，R]上，将这个区间等分成n等份，每一等份的大小为，为了求出解曲线上的这个点的值，故分别取为：.代入(3)式得到一个含个未知数线性方程组: ………………………………………………………………………………..(4) 方程组(4)含有个方程,为了解出这个线性方程组,必须再增加两个条件,即左边界和右边界条件,这里先列出左边界条件的情况,将(1)式代入 得到: 即: (5) 右边界有两种情况,下面将分别讨论. 2.2.1第一种右边界条件 因为 所以 …………(6) 故由(4)(5)(6)可以得到一个阶线性方程组: …………………………………………………………………………………………(7) 2.2.2第二种右边界条件 根据导数的定义得: 即: 也即: …………………(8) 故由(4)(5)(8)可以得到一个n阶线性方程组: ………………………………(9) 在Matlab程序实现中将采用追赶法来解出即为微分方程解曲线上的点. 三 Matlab程序实现 3．1 Matlab工具简介 MATLAB是一种功能非常强大的科学计算软件。它源于MATrix LABoratory一词，愿意为矩阵实验室。一开始它是一种专门用于矩阵数值计算的软件。随着MATLAB逐渐市场化，它不仅具有了数值计算功能，而且具有了数据可视化功能。在MATLAB6.x中，它不仅在数值计算、符号运算