课程论文有限差分法在微分方程中的应用本学期学习了《微分方程数值解》 ,本书中有限差分法给我留下的印象比较深刻,下边说说自己在方面的一点理解,请老师指正
有限差分法的基本思想:当系统的数学模型建立后, 我们面对的主要问题就是微分积分方程的求解
基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的微分方程和定解条件, 并把差分方程组的解作为微分方程定解问题的近似解
将原方程及边界条件中的微分用差分来近似,对于方程中的积分用求和或及机械求积公式来近似代替,从而把原微分积分方程和边界条件转化成差分方程组
有限差分法求解偏微分方程的步骤:区域离散,即把所给偏微分方程的求解区域细分成由有限个格点组成的网格,这些离散点称作网格的节点;近似替代,即采用有限差分公式替代每一个格点的导数
逼近求解,换而言之,这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程
从原则上说, 这种方法仍然可以达到任意满意的计算精度
因为方程的连续数值解可以通过减小独立变量离散取值的间格,或者通过离散点上的函数值进行插值计算来近似得到
理论上, 当网格步长趋近于零时, 差分方程组的解应该收敛于精确解,但由于机器字节的限制, 网格步长不可能也没有必要取得无限小,那么差分法的收敛性或者说算法的稳定性就显得至关重要
因此,在运用有限差分法时,除了要保证精度外,还必须要保证其收敛性
构造差分法的几种形式:主要草用的是泰勒级数展开的方法
其基本差分表达式主要有三种形式:一阶向前差分、 一阶向后差分、 一阶中心差分和二阶中心差分等
其中前两种形式为一阶计算精度,后一种为二阶计算精度
有限差分法的应用 : 4
1 抛物线形的差分法中的一维常系数抛物线型方程考虑最简单的以为常系数抛物线型方程22( )uuLuaf xtx(, )x t其中是( x
t)平面内的给定区域,可以是有节