matlab中差分画图,有限的差分法求解偏微分方程MATLAB.doc

实用标准文案

精彩文档

南京理工大学

课程考核论文

课程名称： 高等数值分析

论文题目： 有限差分法求解偏微分方程

姓 名： 罗 晨

学 号： 115104000545

成 绩：

任课教师评语：

任课教师评语：

签名：

年 月 日

有限差分法求解偏微分方程

一、主要内容

1.有限差分法求解偏微分方程，偏微分方程如一般形式的一维抛物线型方程：

具体求解的偏微分方程如下：

2.推导五种差分格式、截断误差并分析其稳定性；

3.编写MATLAB程序实现五种差分格式对偏微分方程的求解及误差分析；

4.结论及完成本次实验报告的感想。

二、推导几种差分格式的过程：

有限差分法(finite-difference methods)是一种数值方法通过有限个微分方程近似求导从而寻求微分方程的近似解。有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。

推导差分方程的过程中需要用到的泰勒展开公式如下：

(2-1)

求解区域的网格划分步长参数如下：

(2-2)

2.1 古典显格式

2.1