《微分方程数值解》实验指导书2010.docVIP

微分方程数值解 实验指导书 李光云 数学与计算科学学院 二O一O 十月 一、概述 本课程实验指导书是根据陆金甫，关治编著的《偏微分方程数值解法(第二版)》编写的。 通过上机实验，可帮助学生理解偏微分方程数值计算方法的基本思想，巩固学生所学过的算法，同时让学生进一步学习和掌握Matlab软件在数值计算方面中的应用。 二、实验环境 本书选择的实验环境是计算机以及软件Matlab(版本6.5以上)一套。 三、实验课时安排 共7次实验，16课时，试验五为综合设计性试验4课时，其它实验2个课时。 四、实验要求 上机完成实验指导书中所规定的内容，自行按实验指导书要求完成程序设计和调试，并提交每次实验的实验报告，附带算法程序清单和算法输出结果。 五、实验考核要求 上机完成试验内容，并提交一份算法程序清单和数值结果。  实验一 双曲型方程的迎风格式 一、实验目的 1、掌握双曲型方程的迎风格式的算法思想,格式稳定的条件； 2、培养编程与上机调试能力。 二、实验课时：2个课时 三、实验准备 1、了解偏微分方程的分类，以及双曲型方程的特点； 2、熟悉求解偏微分方程有限差分法的基本原理，双曲型方程迎风格式的推导和特点，熟悉用迎风格式求解双曲型方程的算法步骤； 3、熟悉Matlab软件的基本命令及其数值计算中的基本命令和函数。 4、了解双曲型方程迎风格式稳定的条件。 四、实验内容 用双曲型方程的迎风格式 其中，分别为时间步长和空间步长。求解初值问题 其中 取.计算区间的值至,并画出解的图形。 五、基本思想及主要步骤 考虑对流方程 其为给定常数。其迎风格式的基本思想是在方程中关于空间的偏导数用在特征线方向一侧的单边差商代替，得到对流方程的迎风格式 用Fourier方法讨论格式的稳定性可知，以上两个格式都是条件稳定的，而且都在时稳定。 程序设计步骤如下： (1)首先由的正负确定用哪个格式，再把格式改写为显示格式； (2)利用空间步长定义一个零向量用于存储值； (3)将初值离散存于向量内； (4)利用差分格式沿着坐标逐层计算值，每计算一次覆盖存储于向量内，直到算到层，然后输出结果； (5)利用得到的向量画出在时的函数图像。 六、实验提示： 此算法可以编写M程序，使得该程序可以实现对不同的，时间步长和空间步长求解，由于不用符号的用到的格式不同，迎风格式在条件时稳定，故可以在程序里加一个判断部分，当不稳定时及时报警跳出。 七、实验结果记录和要求 1、实验结果输出，画出解的图形； 2、把得到的解和初值问题的解析解比较，观察比较结果； 3、递交实验报告。 八、思考： 若考虑对流方程的差分格式 是否也能得到差分格式稳定？   实验二 常系数扩散方程的经典差分格式 一、实验目的 1、了解抛物型方程的经典差分格式,格式稳定的条件； 2、掌握常系数扩散方程初边值问题的加权隐式格式的求解方法； 3、培养编程与上机调试能力。 二、实验课时：2个课时 三、实验准备 1、了解偏微分方程的分类方法，抛物型方程的特点； 2、了解常见的抛物型方程的差分格式； 3、熟悉抛物型方程的加权隐式格式及其推导过程； 4、熟悉常系数抛物型方程初边值问题的初值和边界离散方法。 四、实验内容 考虑常系数扩散方程的初边值问题 其解析解为 用加权隐式格式近似求解 其中，取为时间步长，为网格比，对不同的时间步长()，计算当时初边值问题的解，并且与精确解比较，分析比较结果。 五、基本思想及主要步骤 用有限差分法解常系数扩散方程 有加权隐式差分格式 其中，当时为Crank-Nicolson格式，当时为向后差分格式，当时为向前差分格式。 加权隐式格式稳定的条件是 加权隐式格式是两层隐式格式，用第n层计算第n+1层节点值的时候，要解线性方程组。 实验步骤如下： (1)输入，确定加权隐式格式的参数； (2)定义向量v，把初边值条件离散，得到的值存入向量v； (3)利用差分格式由第n层计算第n+1层，建立相应线性方程组，求解并且存入向量v； (4)计算到，输出。 六、实验提示： 求解程序可编写M文件，以便改变，得到不用参数下的值，然后比较与真实值的误差，判断格式的稳定性。 七、实验结果记录和要求 1、输出不同参数下求得的节点的近似值； 2、比较近似值和真实值的误差，印证格式稳定条件； 3、实验结束后，递交实验报告。 八、思考： 抛物型方程和双曲型方程的初边值问题的边界处理有什么不同，用隐式格式计算与显示格式相比，有什么优劣。   实验三 椭圆型方程的五点格式 一、实验目的 1、掌握poisson方程定解问题的五点格式的求解方法； 2、培养编程与上机调试能力。 二、实验课时：2个课时 三、实验准备 1、了解偏微分方程的分类方法，椭圆型方程的特点； 2、了解常见的poisson方程的差分格式； 3