第一章数值计算方法与误差分析分析.ppt

计算积分In=∫01 xn ex-1dx 的改进方法 根据 In= 1– nIn-1 可知 In-1 =(1– In)/n 再由 e-1/10 I9 1/10 取 I9 ≈ (e-1/10 +1/10)/2 =0.0684 按 I9 ≈ (e-1/10 +1/10)/2 =0.0684 In-1=(1– In)/n （n=9,8,7,…,1) 计算，保留四 位小数依次得I9 , I8 , I7 , I6 , I5 , I4 , I3 , I2 , I1 , I0为： 0.0684, 0.1035, 0.1121, 0.1268, 0.1455, 0.1709, 0.2073, 0.2642, 0.3679, 0.6321 I0 = 0.6321全部为有效数字。为什么会有这么好的计算 结果呢？下面我们分析一下这种方法的误差积累。 计算积分In=∫01 xn ex-1dx 改进方法的误差分析 根据 I9 ≈ (e-1/10 +1/10)/2 =0.0684 In-1=(1– In)/n （n=9,8,7,…,1) 计算I9 , I8 , I7 , I6 , I5 , I4 , I3 , I2 , I1 , I0 ，假定I9与准确值的误差为ε(I9) 。那么， ε(I8)= – ε(I9)/9 ε(I7)= – ε(I8)/8 = –(– ε(I9)/9)/8 = ε(I9)/(9.8) ε(I6)= –ε(I7)/7= –(–ε(I8)/8)/7= –(–(–ε(I9)/9)/8)/7= –ε(I9)/(9.8.7) … … … … … ε(I0)= –ε(I1)/1= –(–ε(I2)/2)/1= … = – ε(I9)/9! （如：I4=(1- I5)/5= (1- (1- I6)/6)/5= (1- (1- (1- I7)/7)/6)/5 = (1- (1- (1- (1- I8)/8)/7)/6)/5 = (1- (1- (1- (1- (1- I9)/9)/8)/7)/6)/5 = 1/5 – … –I9/(9.8.7.6.5) ) 这说明由初始值I9的舍入而产生的误差在计算过程中绝对值会逐渐缩小。若 |ε(I9)|=5×10-5 则 |ε(I0) |= 5×10-5/9! 故这个算法是稳定的。 * 数值计算方法：研究怎样利用计算工具求出数学问题的数值解，并对算法的收敛性、稳定性和误差进行分析计算的全过程。 构建一个完整的数值算法，包含着以下环节： 1. 提出数值问题（即对对象建立数学模型） 2 .构思处理数值问题的基本思想（即提出理论） 3 .列出计算公式 4 .设计程序框图 5 .编制源程序并调试 6 .做出算法的误差分析 从工程实际中抽象出来的数学问题往往很复杂，典型的有： 1、数据点的插值 2 、曲线拟和 3、复杂函数的微积分运算 4、非线性方程f(x)=0的根的求解 5、当n很大时线性方程组AX=B的求解 6、常微分方程的求解 参考书籍的几种名称： 1、数值分析 2、数值计算原理 3、计算方法 4、算法设计 5、计算机数值计算方法与程序设计 数值计算中的误差 1、误差的种类和来源 ① 模型误差 ② 观测误差 ③ 截断误差 ④ 舍入误差 2、误差的有关概念： ① 绝对误差： 真值 近似值 ② 绝对误差限： ③ 相对误差： ④ 相对误差限： 相对误差限的概念 在实际计算中, 由于ε(x)与x都不能准 确地求得,因此相对误差 εr(x)也不可能准确 地得到, 我们只能估计它的大小范围。即指 定一个适当小的正数η，使 |εr(x)|= |ε(x)|/| x* | ≤η 称η为近似值x*的相对误差限。 当|εr(x)|较小时，可以用下式来计算η： η=ξ/|x*| 有效数字 为了既能表示近似数的大小，又能表示近似数的精确程度，我们下面介绍有效数字的概念（注意：有效数字既能表示近似数的大小，又能表示近似数的精确程度）。 半个单位的概念