工程应用数学案例.docVIP

实验题（实验题1必做，在实验题2、3、4、5中选做一题） 1.(20) 求解微分方程初值问题。 自选方法，在步长，时求解，并与精确解进行比较；(2) 在步长和求解的基础上，运用外推方法提高求解精度。 解：用Euler法求解初值问题 yi+1=yi + h(-yi-xi-2) 其中xi=kh (0x1) 现分别选择h(1)=0.05, h(2)=0.05/2进行计算，然后运用外推法提高求解精度 yi (外推)=(20yh(1)(x)-yh(2)(x))/(20-1) 下表中，第一列表示区间[0,1]进行20等分后得到的节点，第二列表示在这些节点处解的精确值，第四列和第六列分别是步长h(1), h(2)时在节点处得数值解，可以看出，步长h(2)时的误差要小于步长h(1)时的误差。第七列是根据步长h(1), h(2)的计算结果得到的外推近似解，与精确解相比，外推解的近似解较好。 xi y(xi) i(h(1)) yi i(h(2)) yi yi (外推) 0.00 1.0000 0 1.0000 0 1.00 1.0000 0.05 0.8525 1 0.8500 2 0.8513 0.8499 0.10 0.7097 2 0.7050 4 0.7075 0.7049 0.15 0.5714 3 0.5648 6 0.5682 0.5646 0.20 0.4375 4 0