数值分析习题六解答.docVIP

习 题 六 解 答 1、在区间［0，1］上用欧拉法求解下列的初值问题，取步长h=0.1。 (1)(2) 解：取h=0.，本初值问题的欧拉公式具体形式为 由初值y0=y(0)=出发计算，所数值结果如下： x0=0，y0=； x1=0.，x2=0.2， 取h=0.，本初值问题的欧拉公式具体形式为 由初值y0=y(0)=出发计算，所数值结果如下： x0=0，y0=； x1=0.， x2=0.2， 指出： 本小题的求解过程中，函数值计算需要用到计算器。 2、用欧拉法和改进的欧拉法(预测－校正法)求解初值问题，取步长h=0.1。 解：(1) 取h=0.，本初值问题的欧拉公式具体形式为 由初值y0=y(0)=出发计算，所数值结果如下： x0=0，y0=； x1=0.，x2=0.2， ， 取h=0. 由初值y0=y(0)=出发计算，所数值结果如下：x0=0，y0=； x=0.1， x2=0.2， 3、试导出解一阶常微分方程初值问题 的隐式欧拉格式 并估计其局部截断误差。 解：在区间[xn,xn+1]上对常微分方程y /(x)=f(x,y)两端同时积分，得 由右矩形公式得 所以有差分格式 对于隐式欧拉法 假定yn＝y(xn)，上式右边的yn＋1＝y(xn＋1)，则 将y／(xn＋1) 按泰勒公式展开，上式为 将y(xn＋1)按泰勒公式展开，得 两式相减，得 即 所以， 指出： 可以用多种方法导出，其中差商法、数值积分方法是简单的方法。 4、验证改进的欧拉公式对任何不超过二次的多项式 准确成立，并说明理由。 分析：①本题所说的改进的欧拉法，是指梯形公式 。 ②在初值问题中，y是解函数。 ③本题要证明的是，如果解函数是，则用梯形公式求出的数值解等于相应的解函数的函数值，而，即要证明 。 ④为了证明结论成立，先建立求解格式。 ⑤注意，，所以。 解：因为 所以。 记，设 改进的欧拉公式为 将上式对i从0到n－1求和并利用初值条件得 则 所以，改进的欧拉法对任何不超过二次的多项式 准确成立。 指出： 通过累加，把递推关系变成了函数关系。 5、对于初值问题 试用（1）欧拉法；（2）改进的欧拉法；（3）四阶经典龙格-库塔法分别求解，并比较之，取。 解：（1）取，本初值问题的欧拉公式具体形式为 由初值y0=y(0)=出发计算，所数值结果如下： x0=0，y0=； x1=0.，x2=0.4，x2=0.6，， 取 由初值y0=y(0)=出发计算，所数值结果如下：x0=0，y0=； x=0.2， x2=0.4， (3) 四阶经典龙格－库塔公式为 在本题中，， 取，计算得 6、用经典四阶龙格－库塔方法求下列初值问题的数值解。 （1）（2） 解（1）四阶经典龙格－库塔公式为 在本题中，， 取，计算得 i 1 1.242 800 000 1.727 548 209 2 1.583 635 920 2.742 051 299 3 2.044 212 913 4.049 181 355 4 2.651 041 652 5.829 210 728 5 2.436 502 273 7.996 012 143 7、选取参数p、q，使得下列公式 具有二阶精度。 解： 所以 则 而y(xi+1)泰勒展开得 …… 指出： 显然，两个式子不能逐项对比。实际上，不形成递推而形成循环，本题应为错题。 7*、选取参数p、q，使得下列公式 具有二阶精度。 解：因为 则 而y(xi+1)泰勒展开得 比较上面两个关系式，如果三项相等，此方法是二阶数值方法，此时 。 8、用亚当姆斯预报－校正系统求解初值问题 取步长h=0.1计算。 解：四阶的亚当姆斯预测—校正系统为 解：将，代入预报－校正系统得 预报： 校正： 用龙格-库塔法先求出，再用亚当姆斯预测-校正法计算得 i 龙格-库塔法 预测-校正法 1 0.1 0.095 162 5 2 0.2 0.181 269 1 3 0.3 0.259 181 6 4 0.4 0.329 680 1 5 0.5 0.393 469 7 6 0.6 0.451 188 9 7 0.7 0.503 415 4 8 0.8 0.550 671 9 9 0.9 0.593 431 3 10 1.0 0.632 121 6 补充题（一） 1、用欧拉公式求解初值问题 当x取步长为h=0.02，用欧拉公式解初值问题0,0.02,0.04,…,0.10时的解。 2、取步长为h=0.2，用欧拉公式解初值问题 。 答 案 1. 解：将代入欧拉公式，得本初值问题的欧拉公式的具体形式为： ，（） 取由初值y0=y(0)=0出发计算，所得数值结果如下： 用欧拉公式求解的计算结果