计算方法作业chapter8.docVIP

习题八 取，用欧拉法、后退欧拉法、中点法、梯形法、改进欧拉法求解下列初值问题： （3）. 解：求精确解，设则原式化为，可得：， 故，已知，所以， 求得各解如下表所示： 根据题意，按欧拉法公式， 后退欧拉法公式：， 中点法公式：，（其中取后退欧拉法值）， 梯形法公式：， 改进欧拉法公式： 计算结果如下表所列： 精确解 欧拉法 后退欧拉法 中点法 梯形法 改进欧拉法 0.1 1.095445115 1.10000 1.09074 1.10000 1.09566 1.09591 0.2 1.183215957 1.19182 1.17408 1.18364 1.18359 1.18410 0.3 1.264911064 1.27744 1.25125 1.26914 1.26544 1.26620 0.4 1.341640786 1.35821 1.32309 1.34291 1.34232 1.34336 0.5 1.414213562 1.43513 1.39018 1.41858 1.41506 1.41640 0.6 1.483239697 1.50897 1.45287 1.48564 1.48427 1.48596 0.7 1.549193338 1.50834 1.51138 1.55416 1.55043 1.55251 0.8 1.612451550 1.64978 1.56577 1.61631 1.61393 1.61647 0.9 1.673320053 1.71778 1.61598 1.67944 1.67508 1.67817 1 1.732050848 1.78477 1.66181 1.73784 1.73415 1.73787 5、利用R-K法求解第1题. 解： Excel的计算表格如下： Xi Yi K1 Xi+h/2 Yi+K1/2 K2 Yi+K2/2 K3 Xi+h Yi+K3 K4 0 1 0.1 0.05 1.05 0.095476 1.047738 0.095229 0.1 1.095229 0.091262 0.1 1.095446 0.091287 0.15 1.141089 0.087818 1.139355 0.087605 0.2 1.18305 0.084494 0.2 1.183217 0.084516 0.25 1.225475 0.081747 1.22409 0.081562 0.3 1.264779 0.079039 0.3 1.264912 0.079057 0.35 1.304441 0.076781 1.303303 0.076621 0.4 1.341533 0.07452 0.4 1.341642 0.074536 0.45 1.37891 0.072622 1.377953 0.072481 0.5 1.414123 0.070697 0.5 1.414216 0.070711 0.55 1.449571 0.069073 1.448752 0.068948 0.6 1.483163 0.067408 0.6 1.483242 0.06742 0.65 1.516952 0.065997 1.516241 0.065886 0.7 1.549128 0.064539 0.7 1.549196 0.06455 0.75 1.581472 0.063299 1.580846 0.063199 0.8 1.612395 0.062008 0.8 1.612455 0.062018 0.85 1.643464 0.060906 1.642909 0.060816 0.9 1.673271 0.059753 0.9 1.673325 0.059762 0.95 1.703206 0.058766 1.702708 0.058684 1 1.732008 0.057728 1 1.732056 7.利用待定系数法确定如下求解公式的系数，使其阶数尽可能高，并导出局部截断误差的单表达式。 (1) (3) 解：(1)、其局部截断误差为：（在处泰勒展开） 要使应取的系数为0。 所以简单局部截断误差为： （3）、其局部截断误差为：（在处泰勒展开） 可解得： 又 简单局部截断误差为： 8 .仿照修正哈明预测-校正公式的推导，对下列预测-校正公式进行修正： （1）预测：来尔恩公式（8—36） 校正：辛普森公式（8—20） 解：上述公式为： 故由两式的局部截断误差表示式知，真解与计算解之差为： 两式相减，假定 得： 代回两式，便得局部截断