数值计算方法上机实习题答案.docxVIP

数值计算方法上机作业热能工程 数值计算方法上机作业 热能工程 设 I ? ? 1 xn dx ， n 0 5 ? x 由递推公式 I n ? ?5I  n?1 1 ，从 I n 0 的几个近似值出发，计算I ； 20 解：易得： I ? ln6-ln5=0.1823, 0 程序为： I=0.182; for n=1:20 I=(-5)*I+1/n; end I 输出结果为： I = -3.0666e+010 20 粗糙估计 I 20 ，用 I  n?1 ? ? 1 I 5  n?1 1 ，计算 I ； 5n 0 因为 0.0079 ? ?1 x 20 dx ? I ? ?1 x 20 dx ? 0.0095 0 6 20 0 5 所以取 I 20 ? 1 (0.0079 ? 0.0095) ? 0.0087 2 程序为：I=0.0087; for n=1:20 I=(-1/5)*I+1/(5*n); end I I = 0.0083 0 分析结果的可靠性及产生此现象的原因(重点分析原因)。 首先分析两种递推式的误差；设第一递推式中开始时的误差为E ? I 0 0  I ? ，递推过程的舍入误差不计。 0 并记 E n ? I ? I ? ，则有 E n n n ? ?5E  n?1 ? ? ? (?5)n E 0 。因为 E ? 20 1 1