(完整版)数值计算方法上机实习题答案.pdfVIP

(完整版)数值计算方法上机实习题答案--第1页 数值计算方法上机作业 热能工程 xn 1．设I 1 dx， n 0 5 x 1 （1） 由递推公式I 5I  ，从 的几个近似值出发，计算 ；I I n n1 n 0 20 解：易得：I ln6-ln5=0.1823, 0 程序为： I=0.182; for n=1:20 I=(-5)*I+1/n; en I I 输出结果为： = -3.0666e+010 20 1 1 （2 ） 粗糙估计 ，用I I  I  ，计算 ；I 20 n1 5 n1 5n 0 x20 x20 因为 0.0079 1 dxI 1 dx0.0095 0 6 20 0 5 1 所以取I  (0.0079 0.0095) 0.0087 20 2 程序为：I=0.0087; for n=1:20 I=(-1/5)*I+1/(5*n); en I I = 0.0083 0 （3 ） 分析结果的可靠性及产生此现象的原因(重点分析原因) 。 首先分析两种递推式的误差；设第一递推式中开始时的误差为E I I，递推过程的舍入误差不计。 0 0 0 n  E I I  E  并记 ，则有E 5E  (5) E 。因为 n n n n n1 0 20 1 1 (5)20 E I ，所此递推式不可靠。而在第二种递推式中E  E ( )nE ，误差在缩小， 0 20 0 5 1 5 n 所以此递推式是可靠的。出现以上运行结果