常微分方程数值解实验报告.docVIP

《常微分方程数值解》 实验报告 姓名： 张 彬 班级： 信计092 学号： 0909290229 实验目的： 应用前面实验程序求解常微分方程初值问题，通过改变的初值，重点了解Euler方法及其改进方法的稳定性及收敛性 实验内容： ；求满足初值条件x=2，y=2的解；x∈[2，4]，步长h=0.05，并通过数值与图像进行比较。 实验原理： 1.前述几种数值计算方法对于初值及步长的选取都有依赖性。 2.对于无法获得精确解的常微分方程初值问题，更需要分析计算结果的稳定性与收敛性。 实验结果： h=0.05 x=2 y=2 欧拉法 表1 Yend= 自变量 欧拉值 精确解 误差估计 2 2 2 0 2.05 1.965934 1.964661 -0.00065 2.1 1.929429 1.927094 -0.00121 2.15 1.890974 1.887814 -0.00167 2.2 1.85111 1.847387 -0.00201 2.25 1.810418 1.806424 -0.00221 2.3 1.769528 1.765583 -0.00223 2.35 1.729109 1.72