验证改进欧拉公式求数值微分的准确性以及验证复合辛普森公式.doc

实验四报告 学院：计信学院 专业：计科 班级： 姓名 学号 实验组 4 实验时间 指导教师 成绩 实验项目名称 验证改进欧拉公式求数值微分的准确性以及验证复合辛普森公式 实验目的 1、编写算法，比较利用改进欧拉公式求出数值微分的近似值与精确值的误差； 2、验证复合辛普森公式。 实验要求 能利用改进欧拉公式算法求解数值微分的近似值，能够将近似值与精确值的曲线变化图在MATLAB上绘出，并能根据图形比较两者的误差。 实验原理 改进欧拉公式算法、复合辛普森算法 实验仪器 MATLAB 7.0运行环境 实验步骤 运用MATLAB的实验环境，根据改进欧拉公式算法的基本原理，编写改进欧拉公式常微分方程的算法,具体编写代码如下所示： 2、在上述代码编译成功并无任何错误的前提下，编写Ｍ文件调用maeuler函数，计算代求的常微分方程。 此处，为了代入方便，在此编写了的函数，如下所示 此函数的目的即是若以后的验证中需要用到形如此类问题的的解时，只需调用 即可，不必再例举出该形式方程。 为了验证复合辛普森公式的精确性，同样此处编写该算法： 同样，为了计算方便，同样定义的函数， 实验内容 验证课本167页的例1 例1：取h=0.1，求解初值问题 验证课本171页的例3 例3：用Euler法、梯形法以及改进Euler法求解初值问题 取h=0.1，计算至x=0.5。 验证课本133页的第2题 计算定积分 实验数据 1、验证例1：输入以下所示的代码，并使用ode45函数求出精确解，并绘出利用改进欧拉公式求出 的近似值，然后利用plot函数绘出两者的曲线 结果如下所示： 绘出的图如下所示： 验证例2：输入以下所示的代码，并使用ode45函数求出精确解，并绘出利用改进欧拉公式求出的近似值，然后利用plot函数绘出两者的曲线 结果如下所示： 绘出的图如下所示： 3、验证复合辛普森公式：当输入n=1时， 得s=0.946145882273587，此时误差很大，所以将n=2代入， 得此时值与精确值相比，还是存在误差，所以将n=3代入 此时可得出精确值。 实验总结 从改进欧拉公式算法得出的图形可得出，改进欧拉公式得出的近似值与精确值 非常的接近，可得知用改进欧拉公式得出的值误差非常小。 指导教师意见 签名： 年 月 日