常微分方程数值解法机上实验报告.docxVIP

常微分方程数值解法上机实验报告实验题1.ODE:，采用Euler折线法、梯形法、RK3、RK4、三步三阶Adams方法求解。要求：（1）绘制误差下降曲线；（2）确定每个方法的“数值阶”实验图像和数据：（1）误差曲线图：①为了比较不同的方法在区间固定、所分区间数N不断变大的情况下，在整个区间上数值解与真解的最大误差下降情况，取、，得到如下图像：图一：不同方法随着区间数N的增大，最大误差的下降曲线②为了探究三步三阶Admas方法的初值选取对后续误差的下降速度的影响，分别用Euler方法给初值，梯形方法给初值，和同阶的RK3方法给初值，同样考虑区间数N变大时，最大误差的下降情况，得到如下图二。 图二：使用不同阶数的方法给出初值对Admas方法的影响（2）数值阶：利用如下方法来计算数值阶：取，T=100，对、、、、，分别计算，再利用上述公式，两两相邻计算数值阶，可得如下表格：Euler方法1.0230763111.0111956971.00551595711.002706721梯形方法1.9887729752.0001562161.9984676212.000017355RK33.1010330853.0511392423.0265982583.013430782RK44.0214762184.0112186024.0057134444.002890934三步三阶Admas2.938319