有限差分法求解偏微分方程matlab.doc

《有限差分法求解偏微分方程matlab.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、.南京理工大学课程考核论文课程名称：高等数值分析论文题目：有限差分法求解偏微分方程姓名：罗晨学号：115104000545成绩：任课教师评语：签名：年月日有限差分法求解偏微分方程word范文.一、主要内容1.有限差分法求解偏微分方程，偏微分方程如一般形式的一维抛物线型方程：具体求解的偏微分方程如下：2.推导五种差分格式、截断误差并分析其稳定性；3.编写MATLAB程序实现五种差分格式对偏微分方程的求解及误差分析；4.结论及完成本次实验报告的感想。二、推导几种差分格式的过程：有限差分法（finite-differencemethods）

2、是一种数值方法通过有限个微分方程近似求导从而寻求微分方程的近似解。有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。推导差分方程的过程中需要用到的泰勒展开公式如下：(2-1)求解区域的网格划分步长参数如下：(2-2)word范文.2.1古典显格式2.1.1古典显格式的推导由泰勒

3、展开公式将对时间展开得(2-3)当时有(2-4)得到对时间的一阶偏导数(2-5)由泰勒展开公式将对位置展开得(2-6)当时，代入式(2-6)得(2-7)因为，代入上式得(2-8)得到对位置的二阶偏导数word范文.(2-9)将式(2-5)、(2-9)代入一般形式的抛物线型偏微分方程得(2-10)为了方便我们可以将式(2-10)写成(2-11)(2-12)最后得到古典显格式的差分格式为(2-13)，古典显格式的差分格式的截断误差是。2.1.2古典显格式稳定性分析古典显格式(2-13)写成矩阵形式为(2-14)上面的C矩阵的特征值是：wo

4、rd范文.(2-15)使，即结论：当时，所以古典显格式是稳定的。2.2古典隐格式2.2.1古典隐格式的推导将代入式(2-3)得(2-16)(2-17)得到对时间的一阶偏导数(2-18)将式(2-9)、(2-18)原方程得到(2-19)为了方便把(2-19)写成(2-20)word范文.(2-21)最后得到古典隐格式的差分格式为(2-22)，古典隐格式的差分格式的截断误差是。2.2.2古典隐格式稳定性分析将古典隐格式(2-22)写成矩阵形式如下(2-23)误差传播方程(2-24)所以误差方程的系数矩阵为使，显然恒成立。结论：对于，即任意

5、网格比下，古典隐格式是绝对稳定的。word范文.2.3Richardson格式2.3.1Richardson格式的推导将，代入式(2-3)得(2-25)即(2-26)由此得到可得(2-27)将式(2-9)、(2-27)代入原方程得到下式(2-28)为了方便可以把式(2-28)写成(2-29)即(2-30)最后得到Richardson显格式的差分格式为(2-31)word范文.，古典显格式的差分格式的截断误差是。2.3.2Richardson稳定性分析将Richardson显格式(2-31)写成如下矩阵形式(2-32)误差传播方程矩阵形

6、式(2-33)再将上面的方程组写成矩阵形式(2-34)系数矩阵的特征值是(2-35)解得特征值为(2-36)(恒成立)(2-37)结论：上式对任意的网比都恒成立，即Richardson格式是绝对不稳定的。4.Crank-Nicholson格式3.4.1Crank-Nicholson格式的推导将代入式(2-9)得word范文.(2-40)即(2-41)得到如下方程(2-42)所以处的一阶偏导数可以用下式表示：(2-43)将,代入式(2-6)可以得到式(2-9)；同理,代入式(2-6)可以得到(2-44)所以，处的二阶偏导数用式(2-6)

7、、(2-44)表示：word范文.(2-45)所以，处的函数值可用下式表示：(2-46)原方程变为：(2-47)将差分格式代入上式得：(2-48)为了方便写成：(2-49)最后得到Crank-Nicholson格式的差分格式为(2-50)，Crank-Nicholson格式的差分格式的截断误差是。3.4.1Crank-Nicholson稳定性分析Crank-Nicholson格式写成矩阵形式如下：word范文.(2-51)误差传播方程是：(2-52)可以得到：(2-53)(2-54)使即(2-55)(2-56)(2-57)(2-58)

8、上式恒成立。结论：Crank-Nicholson格式对任意网格比也是绝对稳定的。word范文.5.DuFortFrankle格式（Richardson格式的改进）将代入式(2-31)并化简得到DuFortFrankle：