常微分方程实验报告_闵安游_u201310085.pdfVIP

华中科技大学 本科实验报告 课程名称：常微分方程数值解 姓 名：闵安游 学 号：U201310085 专业班级: 信计1301 班 指导老师：黄乘明 实验题目：常微分方程数值解上机 日 期：2015.12.25 实验成绩： 1. 实验目的 通过上机实验，理解利用计算机迭代求解线性方程组的整个过 程，加深对所学计算方法的理论及算法特点的理解。 2. 实验题目 1、分别应用2 步2 阶与3 步3 阶显式Adams 方法求解初值问题 y  2 y 2t 1, t [0, 2], y(0) 1 h 0.2 0.1 t 2 列出当步长 和 时计算到终点 的数值解及其误差，其中 的初值考虑两种方案，一种是使用显式Euler 方法，另一种是直接 2t 采用准确值(其真解为y  e t ) 。要求对所获得结果进行分析，并 附两个算法的MATLAB 程序。 2 、画出2 步3 阶和3 步4 阶隐式Adams 方法的边界轨迹，并确定 方法的绝对稳定区域。 3、应用2 步3 阶和3 步4 阶隐式Adams 方法求解初值问题 u 9u 24u 5cos t 1/ 3sin t 1 1 2 u 24u 51u 9cos t 1/ 3sin t 2 1 2 u (0) 4 / 3,u 2 / 3 1 2 h  0.2, 0.1, 0.05 t 2 列出当步长 时计算到终点 的数值解及其误差，其 中初始值直接采用准确值，该问题的真解是 u (t)  2e3t e39t 1 / 3cos t 1 u (t)  e3t 2e39t 1 / 3cos t 2 要求对所获得结果进行分析，并附两个算法的MATLAB 程序。 4 、应用显示中点方法和4 级4 阶Gill 方法求解初值问题 y y t ,t [0,1],y (0) 3 列出当步长 和 时计算到终点 的数值解及其误差，并附 h 0.2 0.1 t 1 t y t e 算法的MATLAB 程序。（问题的真解为 ( )  2 t 1 ） 5、应用3 级3 阶Heun 方法求解初值问题 y   y  t  t  2 2 1, [0, 2], y(0) 1 h 0.2 0.1 t 2 列出当步长 和 时计算到终点 的数值解及其误差，并附 2t 算法的MATLAB 程序.(该问题的真解为y  e t ，其中方法所需要 的初值直接使用准确值) 3. 程序文本 第 1 题: ⑴、2 步2 阶Adams 方法： ①、初值由真解给定： 主程序： function Twosteps_