第一章误差理论.ppt

一、误差来源、分类 ;在用数值方法解题过程中可能产生的误差归纳起来有如下几类： 1. 模型误差 2. 观测误差 3. 截断误差 4. 舍入误差;用数学方法解决一个具体的实际问题，首先要建立数学模型，这就要对实际问题进行抽象、简化，因而数学模型本身总含有误差，这种误差叫做模型误差 数学模型是指那些利用数学语言模拟现实而建立起来的有关量的描述 数学模型的准确解与实际问题的真解不同;在数学模型中通常包含各种各样的参变量，如温度、长度、电压等，这些参数往往是通过观测得到的，因此也带来了误差，这种误差叫观测误差 数学模型中的参数和原始数据，是由观测和试验得到的 由于测量工具的精度、观测方法或客观条件的限制,使数据含有测量误差,这类误差叫做观测误差或数据误差 根据实际情况可以得到误差上下界 数值方法中需要了解观测误差,以便选择合理的数值方法与之适应;精确公式用近似公式代替时,所产生的误差叫截断误差 例如, 函数f(x)用泰勒(Taylor)多项式 ;在数值计算中只能对有限位字长的数值进行运算 需要对参数、中间结果、最终结果作有限位字长的处理工作，这种处理工作称作舍入处理 用有限位数字代替精确数，这种误差叫做舍入误差，是数值计算中必须考虑的一类误差;误差;二、误差、有效数字;定义2 ;例2 考察三位有效数字重力加速度g,;解1：若取近似值x*=3.1415，绝对误差是0.0000926…，有;定理１ ;例5 指出下列各数具有几位有效数字，及其绝对误差限和相对误差限： －0.002 00 9 000.00 解 因为x1*=－0.002 00， m=－3 绝对误差限0.000 005= 因为m=－3，n=3， x1*= －0.002 00有3位有效数字. a1=2， 相对误差限?r=;三、数值运算的误差估计;;;;误差分析简介 ;一、病态问题与条件数 ;二、算法的数值稳定性;控制递推公式中误差的传播 ; ; ;准确的理论递推式 实际运算的递推式 从而有 ;算法的数值稳定性;设计算法时应遵循的原则 要有数值要稳定性,即能控制误差的传播. 避免大数吃小数,即两数相加时,防止较小的数加不到较大的数上. 避免两相近的数相减,以免有效数字的大量丢失. 避免分母很小(或乘法因子很大),以免产生溢出.;三、避免误??危害的若干原则;减少运算误差原则;减少运算误差原则;减少运算误差若干原则;;例11 仿计算机在3位十进制下， ; 例 求二次方程x2-105x+1=0的根 解：按二次方程求根公式 x1=(105+(1010-4)1/2)/2 x2=(105-(1010-4)1/2)/2 在8位浮点数计算得 x1=(105+105 )/2=105 (正确）, x2=(105-105 )/2=0 (错误） 产生错误的原因 ① 出现大数1010吃掉小数4的情况 ② 分子部分出现两个相近数相减而丧失有 效数位常称为灾难性的抵消;4、绝对值太小的数不宜做除数 当分母为两个相近数相减时,会丧失有效数字;例 计算; 算法2。分成三组因子。每组只取六位小数计算 a=0.0005/0.0003=1.666667(有舍入) b=0.0143/0.0125=1.144000 c=0.0012/0.0135=0.088889 (有舍入) D=a*b*c=1. 666667* 1.144000* 0.088889 =0.169482,准确到小数后5位。;5、简化计算步骤，减少运算次数减少运算次数可以不但节省时间，而且减少舍入误差 例：x255=xx2x4x8x16x32x64x128 原先要做254次乘法现只需14次即可 例 如计算多项式 p(x)=anxn? an-1xn-1 ? … ? a1x ? a0 的值 若直接计算akxk,再逐项相加，一共要做 n+(n-1)+…+2+1=n(n+1)/2次乘法和n次加法 ;如果将前n项提出x，则有 p(x)=（anxn-1? an-1xn-2 ? … ? a1 ）x ? a0 =((anxn-2?an-1xn-3?…? a2)x?a1）x ? a0 =(…(anx ? an-1)x?…?a2)x ? a1)x ? a0 写成递推公式 ;本章小结 ; 规格化浮点数 x= ? 0.a1 a2...at?10c ai?{0,1,2,…,9}, a1?0，L?c ?U ;计算机中数的计算特点: 1. 加法先对阶,