数值分析第8节答案.doc

第八章 常微分方程初值问题数值解法 1、解：欧拉法公式为 代入上式，计算结果为 2、解：改进的欧拉法为 将代入上式，得 同理，梯形法公式为 将代入上二式，，计算结果见表9—5 表 9—5 改进欧拉 梯形法 0．1 0．2 0．3 0．4 0．5 0．005500 0．021927500 0．050144388 0．090930671 0．144992257 0．005238095 0．021405896 0．049367239 0．089903692 0．143722388 可见梯形方法比改进的欧拉法精确。 3、证明：梯形公式为 代入上式，得 解得 因为，故 对，以h为步长经n步运算可求得的近似值，故代入上式有 4、解：令，则有初值问题 对上述问题应用欧拉法，取h=0.5,计算公式为 由得 5、解： 四阶经典龙格-库塔方法计算公式见式（9.7）。对于问题（1）， ；对于问题（2），。取h=0.2,,分别计算两问题的近似解见表9-6。 表 9-6 （1）的解 （2）的解 0.2 0.4 0.6 0.8 1.0 1.242800000 1.583635920 2.044212913 2.651041652 3.436502273 1.727548209 2.742951299 4.094181355 5.829210728 7.996012143 6、证明： 根据定义9.2，只要证明即可。而 因此只须将和都在x处展开即可得到余项表达式： 所以 故对任意参数t，题中方法是二阶的。 7、解： 因此，中点公式是二阶的。 对模型方程使用中点公式求解，得 易知，当时，中点公式绝对稳定。特别当为实数且时，上不等式的解为 8．解： （1）用欧拉法求解题中初值问题，当满足 时绝对稳定，即当时欧拉法绝对稳定。 （2）当满足不等 时，四阶龙格-库塔法绝对稳定，也即当满足 时绝对稳定。 （3）对于梯形公式，当时，绝对稳定，此条件对都成立，即梯形法对无限制。 9、解： 二阶阿达姆斯显式和隐式方法分别为 将代入上二式，化简得 显式方法 隐式方法 取，计算结果如表9-7所示 表 9-7 显式 隐式 0．4 0．6 0．8 1．0 0．3267 0．44679 0．545423 0．6264751 00．451743801 0．551426746 0 可见，隐式方法比显式方法精确。 10．证明：根据局部截断误差的定义知 故方法是二阶的，局部截断误差的主项为。 11、解 由局部截断误差的定义知 所以当时 方法为二阶；当时 方法为三阶。 12、解： 根据刚性比的定义，若方程组的矩阵的特征值满足条件则 称为刚性比，易知A的两个特征值为 所以刚性比s=20。 当时，数值稳定。因此当时才能保证数值稳定。 1 数值分析第八章