1.4数值算法的稳定性.ppt

有效数字 如果|e| = |x* - x|? 1/2 ?10m-n 称近似数x准确到小数点后第n位,从这小数点后第n位数字直到最左边非零数字之间的所有数字都称为有效数字. 截断误差用我们原有的数学思维方式就比较容易理解的，而舍入误差则是本课程引入的一个新概念。要真正理解舍入误差，特别是它在计算中的传播及最终对计算结果的影响，是初步具备科学计算能力的重要标志。希望大家在完成实验后，认真仔细去体会截断误差和舍入误差的含义及对计算结果的影响。 图 0 10-7 10-5 10-5 10-4 10-3 10-2 10-1 10-3 10-1 误差 步长h 实验分析：不论采用怎样的算法，计算结果通常都会有误差。比如算法(1)，由Taylor公式，知： 所以有 利用上式来计算f ’(x0)，其截断误差为： 所以误差是存在的，并且当步长h越来越小时，(1)的近似程度也越来越好。 截断误差 */32 线性代数一④ 设x*是准确数, x是x*的近似数，称e = x* - x 为近似值x的绝对误差，简称误差。 上节课内容回顾 ——反映是近似值与精确值的绝对差值 称 为近似值x的相对误差 ——反映是近似值与精确值的近似程度 通常用百分数来表示，相对误差越小，近似程度越高 绝对误差限，相对误差限 则称?r 为近似值x的相对误差限。 |e| = |x* - x|? ?， 称? 为近似值x的绝对误差限，简称误差限或精度 有效数字越多,误差越小,计算结果越精确. 近似数x=±0.a1a2…an×10m 3.142=0.3142×101， 3.1416=0.31416×101， 3.141=0.3141×101 的近似值如下 例： 用四舍五入得到的数都是有效数字. 4位有效数字 5位有效数字 3位有效数字 m=1, n=4 m=1, n=5 m=1, n=3 |e| = |x* - x|? 1/2 ?10m-n 相对误差与有效数字的关系如下： 定理1.1 设近似数x=±0.a1a2…an×10m有n位有效数字,则其相对误差限为 定理1.2 设近似数x=±0.a1a2…an×10m的相对误差限为 则它至少有n位有效数字。 三、四则运算的误差计算 绝对误差：e = x* - x =?x?dx 相对误差：er =(x*-x)/x* ?dx/x=dlnx 利用这个关系可以讨论四则运算的误差和函数的误差 例如下列式子说明什么误差结果? d(x+y)=dx+dy dln(xy)=dlnx+dlny dln(xn)=ndlnx 思考: 若 则 作业：P12—13 1、3、 4、6 四则运算误差限的公式: 近似计算中误差是不可避免的 ，但不能无限扩大。如何控制误差的传播，是设计算法时必须考虑的问题。本节将要介绍的内容就是如何控制误差的传播，给出一般的原则。 解决一个计算问题往往有多种算法，用不同算法计算的结果其精度往往大不相同。这是因为初始数据的误差或计算中的舍入误差在计算过程中的传播造成的，这种传播因算法不同而相异。 §1.4 算法的数值稳定性 (数值计算中值得注意的问题) 一个算法如果输入数据有误差，而在计算过程中舍入误差不增长，则称此算法是数值稳定的，否则称此算法为不稳定的。 换句话说：若误差传播是可控制的,则称此算法是数值稳定的,否则称此算法为不稳定的。 见教材第2、10—11页 例：计算 (1)In0; (2)In单调递减 In有以下性质 在该例中,用上述公式计算积分的值,I0=ln6-ln5 ≈0.182322的舍入误差在计算过程迅速传播，每次扩大5倍，致使I12= -010-2 严重失真，所以这一公式是不稳定的。 有递推公式 -6354.23 20 2.04388 15 0.0147169 10 0.0284686 5 1270.86 19 -0.395442 14 0.0170566 9 0.0343063 4 -2540161 18 -0.0933742 13 0.0212378 8 0.0431387 3 50.8433 17 -012 0.0188109 7 0.058038 2 -10.1569 16 0.0173247 11 0.0243239 6 0.0883922 1 In n In n In n In n 舍入误差在计算过程迅速传播，每次扩大5倍. 所以此算法是不稳定的。 然后取充分大的m对应的Im的一个估计值为计算初值，再逐步用上式算出Im-1 ，Im-2 ，...，I1。 用