1、I新疆大学 毕业论 文 (设计 )题 目: 求解热传导方程的高精度隐式差分格式 所属院系： 数学与系统科学学院 专 业： 信息与计算科学 II声 明本人郑重声明该毕业论文（设计）是本人在开依沙尔老师指导下独立完成的，本人拥有自主知识产权，没有抄袭、剽窃他人成果，由此造成的知识产权纠纷由本人负责。声明人（签名）： 年 月 日亚库甫江.买买提同学在指导老师的指导下，按照任务书的内容，独立完成了该毕业论文（设计） ，指导教师已经详细审阅该毕业论文（设计） 。指导教师（签名）： 年 月 日III新 疆 大 学毕业论文（设计）任务书班 级：信计 07-2 姓 名：亚库甫江.买买提论文（设计）题目：求解热
2、传导方程的高精度隐式差分格式 专 题： 毕业设计 论文（设计）来源： 教师自拟 要 求 完 成 的 内 容 ： 学习和掌握一维热传导方程已有的各种差分格式的基础上，扩散方程对空间变量应用紧致格式离散，对时间变量应用梯形方法，构造热传导方程的精度为 数值格式，讨24Oh论格式的稳定性，最后数值例子来验证。发题日期：2012 年 12 月 25 日 完成日期：2012 年 5 月 28 日实习实训单位： 数 学 学 院 地点： 数 学 学 院 论文页数： 19 页； 图纸张数： 4 指 导 教 师 ： 开依沙尔 老 师 教研室主任 院 长 （系 主 任 ） IV摘要本文首先对热传导方程经典差分格式
3、进行复习和讨论，然后热传导方程对空间变量四阶紧致格式进行离散,时间变量保持不变，把一维热传导方程转化为常微分方程组的初值问题, 再利用梯形方法构造热传导方程方程的时间二阶空间四阶精度的一种差分格式，并稳定性进行分析，数值结果与 Crank-Nicholson 格式进行比较，数值结果表明, 该方法是有效求解热传导方程的数值计算 .关键词: 热传导方程,高精度紧致格式 ; 梯形方法；两层隐格式; Crank-Nicolson 格式ABSTRACTThis paper first study on some classical finite difference for the heat condu
4、ction equation, secondely secondely we apply compact finite difference approximation of fourth order for discretizing spatial derivatives but leave the time variable Continuous. This approach results in a system of ODEs, which can then be used trapezodial formula derived fourth order in space and se
5、cond order in time unconditionally stable implicit scheme .the stability and local truncation error of the obtained method are analysied. Numerical experiments shows that this method Useful, efficient method for solving diffusion equationKeywords: Heat conduction eqution;Higher- oder compact scheme;
6、 Trapezodial formula ;Two- level implict scheme; Crank- Nicolson schemeV目 录引言.1预备知识.21.扩散方程的经典差分格式.31.1 显式差分格.31.1.1 显式的截断误差. .41.1.2 显式差分格式的稳定性.41.2 隐式差分格式.51.2.1 隐式差分格式的截断误差.51.2.2 隐式差分格式的稳定性.61.3 Crank-Nicolson 格式.61.3.1 Crank-Nicolson 差分格式的截断误差.71.3.2 Crank-Nicolson 差分格式的稳定性.82.高精度格式的构造.92.1 梯形方
7、法.92.2 本文格式的构造.102.3 稳定性分析.113.数值实验.13结论.17致谢.18参考文献.19VI1引言热传导方程是一类描述物理量随时间的扩散和衰减规律的抛物型微分方程自然环境、工程设备及生物机体中的许多物理现象，诸如气体的扩散、液体的渗透、热的传导、以及半导体材料中杂质的扩散等都可用热传导方程来描述由于物理问题本身的复杂性，其精确解往往不容易求得，因此研究其数值求解方法无疑具有非常重要的理论意义和工程应用价值 【1】 求解该问题的数值方法主要有 差分法、有限元法、边界元法等，其中有限差分方法数值求解扩散方程的应用广泛的有效地方法之一。目前求解该问题的主要的差分格式有显式格式，
8、隐式格式，Crank-Nicolson 格式等 1,2,4。 虽然显式格式计算简单，但是稳定性有所限制，一般隐式格式和Crank-Nicolson 格式分别为一阶和二阶精度的绝对稳定的隐式格式，还显得误差阶不够高， 得到的结果也往往不能令人满意, 考虑到这些不足文7中半离散方法构造格式结果Crank-Nicolson 格式进行比较，在文 10待定参数法构造2Oh精度 的显式格式但是稳定性条件比较苛刻，它文的稳定性条件为36，本文热传导方程对空间变量应用紧致格式离散，对时间变量应用梯21arh形方法，构造热传导方程的精度为 的绝对 稳定的隐式差分格式，并24Oh讨论了稳定性，数值值结果与经典Cr
9、ank-Nicholson 格式进行比较，数值结果表明,该方法是有效求解扩散方程的数值计算.本文分为三大部分，第一部分简单介绍热传导方程的经典差分格式，第二部分主要介绍热传导方程的高精度格式的构造和稳定性，第三部分给出具体的数值算例，结果与Crank-Nicolson格式，准确值进行比较，最后给出结论。2预备知识利用下面的各种数值微分公式得到不同的差分格式111 21111 21 12,2,njnjnjnjnjnjnjnjnjnjnjnjnjnjnjjnjnjnuxtuxtuOttttuxtuxtuhhxttOuxtuxtuhhxtttu 22njhx截断误差：一般说来，微分方程的解不会精确地
10、满足差分方程。将差分方程中的各个项同时用微分方程的解在相应点的值代入，利用泰勒展开，就会得到一个误差项，这个误差项就是截断误差。相容性：若时间步长 以及空间步长 同时趋于 ，截断误差 ，就h00njR说差分格式与微分方程是相容的。一个差分格式与一个微分方程相容，则表明当 时，差分算子与微分算子对任一光滑函数的作用是相同的，所以可用0,h相容的差分格式近似相应的微分方程，而截断误差则是对这一近似程度的一个度量。收敛性：考察差分格式在理论上的准确解能否任意逼近微分方程的解。如果当时间步长 以及空间步长 趋于 时， ，我们称差分h00),(njjnj utxe格式是收敛的，即时间步长 以及空间步长
11、趋于 时，差分格式的解逼近于微h0分方程的解。稳定性：差分格式的计算是逐层计算的，计算第 层上的 时，要用1n1nj到第 层上计算出来的结果。计算 时的舍入误差，必然会影响到 的值，nnjuu从而就要分析这种误差传播的情况。因此，一个有实用价值的数值方法应该具有能够控制这种误差影响的性能，这就是数值方法的稳定性。精度：如果一个差分格式的截断误差 ，就说差分格式对时间)(pqhOE是 阶精度的，对空间 是 阶精度的。tqxpLax 等价定理 ：给定一个适定的线性初值问题以及与其相容的差分5式，则差分格式的稳定性是差分格式收敛性的充分必要条件。3定理1（von Neumann条件） 微分方程的差分
12、格式稳定的必要条件是当 ，0,对所有 有TnRk， MkGj 1),( pj,21其中 为增长因子（或增长矩阵） ， 表示 的特征值，),(G),(kGj)(k为常数。M定理2 如果差分格式的增长矩阵 是正规矩阵，则 von Neumann 条),(件是差分格式稳定的必要且充分条件。推论2.1 当 为实对称矩阵，酉矩阵，Hermite矩阵时，von Neumann ),(k条件是差分格式稳定的充分必要条件。推论2.2 当 时，即 只有一个元素，则von Neumann 条件是差1p),(kG分格式稳定的充要条件。定理3 如果存在常数 使得0,K， ，k)( 0则差分格式是稳定的。1. 热传导方
13、程的经典差分格式考虑一维热传导方程的初边界问题： 20, ,0(,)(,uaxbttxfat tbt 1.1 显式差分格式我们可以对 用向前差分ut1,njnjnjuxtuxtuOt用二阶差商2x21 1 22,njnjnjnjuxtuxtxtuhhx 得到差分格式为(1.1.1)11120nnnjjjjjuuuh41.1.1 显式差分格式的截断误差证：（用 taylor 展开）(,)njjnuxt 121 1(,)(,),(0)!nn nj jnjjuuxtuxttt )10(!4!3!2),(),( )(!, 31 24321 njnjnjnjnjnj njnjnjjnjnj xuhxuh
14、xuhxutxutxutt把上述代入差分格式中，得截断误差为： 1 2 322244(,)nnn nnnj j j j j juuhuhukxtttxxx )1,0().( 242 3212 442 321 hotxt njnjnjnjnj从上述可知，截断误差为 ，它对空间方向为一阶截断误差而2(,)njkxtoh对时间方向为二阶截断误差。1.1.2 显式差分格式的稳定性证：先把差分格式公式(1.1.1)改写为：其中111(2)nnnnjjjjjuuu 2h利用稳定性的 Fourier 方法，令 ，并将它代入上式就得到ikjhjve1()nikjhnikjhijikjnve v 消去共因子有 1(2)nikjikjne由此得到增长因子 (,)(2)(1)()1)cos2cosikjikj ikjikjGkeehh 2sin012sin41 1|4|)(|,( 2kkk
