试论《有限差分法在微分方程中的应用》课程论文.docVIP

　　试论《有限差分法在微分方程中的应用》课程论文 试论《有限差分法在微分方程中的应用》课程论文_论文范文导读： 课 程 论 文 有限差分法在微分方程中的应用本学期学习了《微分方程数值解》 ，本书中有限差分法给我留下 的印象比较深刻，下边说说自己在方面的一点理解，请老师指正。 1.有限差分法的基本思想：当系统的数学模型建立后，我们面对的主要问题就是微分积分方程的求解。 基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的 微分方程和定解条件， 并把差分方程组的解作为微分方程定解问题的近似解。 将 原方程及边界条件中的微分用差分来近似， 对于方程中的积分用求和或及机械求 积公式来近似代替，从而把原微分积分方程和边界条件转化成差分方程组。 2.有限差分法求解偏微分方程的步骤： 区域离散，即把所给偏微分方程的求解区域细分成由有限个格点组成的X 格，这些离散点称作X格的节点； 近似替代，即采用有限差分公式替代每一个格点的导数。 逼近求解， 换而言之， 这一过程可以看作是用一个插值多项式及其微分来代 替偏微分方程的解的过程。 从原则上说， 这种方法仍然可以达到任意满意的计算精度。 因为方程的连续数值 解可以通过减小独立变量离散取值的间格， 或者通过离散点上的函数值进行插值 计算来近似得到。 理论上，当X格步长趋近于零时，差分方程组的解应该收敛于 精确解，但由于机器字节的限制，X格步长不可能也没有必要取得无限小，那么 差分法的收敛性或者说算法的稳定性就显得至关重要。 因此，在运用有限差分法 时，除了要保证精度外，还必须要保证其收敛性。 3.构造差分法的几种形式： 主要草用的是泰勒级数展开的方法。 其基本差分表达式主要有三种形式：一 阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等。 其中前两种形式 为一阶计算精度，后一种为二阶计算精度。 4.有限差分法的应用: 4.1 抛物线形的差分法中的一维常系数抛物线型方程 考虑最简单的以为常系数抛物线型方程u u Lu a 2 f ( x) t x2 ( x ,t ) 其中 是（x.t）平面内的给定区域，可以是有节区域或无解区域；agt;0 是常数，L 是微分算 子。 根据定解条件的不同，可以将上述方程分为两类： 1.初值问题 在区域 {( x, t ) | x , t 0} 上求解方程满足初始条件 u( x,0) ( x), x 的解。 2，初边值问题（混合问题） 在区域 {( x, t ) | 0 x l ,0 x T } 内求方程满足初始条件 u( x,0) ( x),0 x l 和下列边界条件之一的解。 第一边届条件 u (0, t ) 1, u (l , t ) 2第二边界条件 0t T u x (0, t ) 1, ux (l , t ) 2第三边界条件 0t T (u 1 (t )ux ) |x 0 r1 (t ) (u 2 (t 试论《有限差分法在微分方程中的应用》课程论文_论文范文(2)导读：f1(a12hb1)u0h2f2g )u x ) |x 0 r2 (t ) i 0 ,i 1 , 2 , 0 t T 用适当的差商代替方程中相应的偏导数，可得到以下几种最简差分格式： 古典显示格式： L u (1) h k j 1 uk uk j j a k k uk j 1 2u j u j 1 h2 f jk 古典隐式格式：k 1 L(2) h uj 1 uk uk j j a 1 k 1 1 uk uk j 1 2u j j 1 h 2 f jk 1 加权六点隐式格式：1 uk uk j j 1 k 1 1 k k uk uk uk j 1 2u j j 1 j 1 2u j u j 1 k 1 k a a (1 ) f j (1 ) f j 2 2 h h L u (3) h k j 4.2 椭圆型方程边值问题的差分法 考虑如下两点边值问题 au bu cu f u (0) u0 , u (1) u1 x (0,1) 其中 u 0 ， u1 为常数，系数 a=a(x),b=b(x),c=c(x),f=f(x)为一致的充分光滑函数，且满足 a(x)gt;0， c(x)gt;0. 首先将区间 离散化， 我们采用剖分部分， 取正整数 M， 将区间 M 等分， 的 M+1 个节点： 0 x0 x1 xM 1其中 x j jh( j 0,1, M ), h 1 。 设 U 为定义在节点 x j jh( j 0,1, M ) 上的X格函 M 数并用 U i 近似 u