第1章误差教案.pptVIP

(4) 数值字长已定，待定近似式项数 初步选定一个项数，根据所定的项数和数值的舍入误差估计出?值、R值，如果它们的值比较接近则选定的项数即为所求；否则另选项数，重复上述过程，直到选出符合要求的项数为止。 1.5 误差分配原则和处理方法 例1.17 对于0?x?1,计算ex 的台劳展开式的部分和 解：根据台劳公式的展开式的截断误差公式: 近似ex ,若在计算机中,x的计算公式及以上公式中各项结果均截取至小数后5位,试确定上式中应取几项为好? 1.5 误差分配原则和处理方法 初取n=3,则有R3(x)3/4!=0.125 ?3=3?(0.5 ? 10-5)=0.000015 因为R3(x)?3,再取n=6,有 R6(x)3/7!=0.00059， ?6=6 ? (0.5 ? 10-5)=0.00003 因为R6(x)?6,增大至n=7,有 R7(x)3/8!=0.000074， ?7=7 ?(0.5 ? 10-5)=0.000035 这时R7(x) ≈?7，所以应取8项进行计算为宜。 1.5 误差分配原则和处理方法 * §1.3 算术运算中的误差 数值计算中求得的解与参量（原始数据）x1, x2,…,xn 有关，计为：y=f(x1, x2,…,xn) xi,yi为准确值, xi*,yi*分别为其近似值; y*=f(x1*, x2*,…,xn*) 由微分学：当自变量改变量（误差）很小时，函数的微分作为函数改变量的主要线性部分可以近似函数的改变量，故利用微分运算公式可导出误差运算公式。 假设： §1.3 算术运算中的误差 绝对误差： §1.3 算术运算中的误差 二元函数绝对误差： §1.3 算术运算中的误差 相对误差： §1.3 算术运算中的误差 二元函数相对误差： §1.3 算术运算中的误差 和、差、积、商之误差公式： 和差的绝对误差限 和差的绝对误差 §1.3 算术运算中的误差 和、差、积、商之误差公式： §1.3 算术运算中的误差 积、商之相对误差限公式： 结论: 和、差的（绝对）误差为各项误差之和、差； 和、差的（绝对）误差限为各项误差限之和; 积、商的相对误差为各项相对误差之和、差； 积、商的相对误差限为各项之相对误差限之和; §1.3 算术运算中的误差 则其相对误差限为: 4. c=xp(p1) 则其绝对误差为: 结论:x的p次幂的相对误差等于x本身的相对误差的p倍. §1.3 算术运算中的误差 则其相对误差限为: 5. c=x1/q(1p=1/q) 结论:x的q次根的相对误差等于x本身的相对误差的1/q倍. ?例1.10 已知球体的直径D=3.7cm,按v=?D3/6计算体积，求其绝对误差限与相对误差限. 1.3 算术运算中的误差 解决方案： （1）改变计算公式：更有效 （2）多保留几位有效数字 1. 避免两个相近的数相减 影响：容易造成有效数字的严重丢失。 §1.4 数值计算中应该注意的问题 取8位有效数字： §1.4 数值计算中应该注意的问题 例如： 改变计算公式： 结果只有1位有效数字 结果有8位有效数字 §1.4 数值计算中应该注意的问题 改变公式的常见情况： 当x接近0时，应作变换： 1-cosx=2sin2(x/2) (1-cosx)/sinx=sinx/(1+sinx) 当x充分大时，应作变换： §1.4 数值计算中应该注意的问题 2. 避免大数吃小数 例：用单精度计算x2-(109+1)x+109=0的根。 精确解为:x1=10 9 ,x2=1 ? 算法1：利用求根公式 §1.4 数值计算中应该注意的问题 在计算机内，109存为0.1?1010，1存为0.1?101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。即1 的指数部分须变为1010，则： 1 = 0.0000000001 ? 1010， 取单精度时就成为： 109+1=01010+0?1010 =0?1010 解决方案：改变计算公式 §1.4 数值计算中应该注意的问题 4. 先化简再计算，减少运算次数，避免误差积累。 5. 选用稳定的算法。 舍入误差对计算结果影响小的算法 3. 避免除数的绝对值远小于被除数的绝对值 故 时，舍入误差可能会很大。 §1.4 数值计算中应该注意的问题 例1.14 计算In=∫01xnex-1dx,n=0,…,7 解：算法1。 用分部积分法可以推知In满足以下递推公式 In=1-nIn-1 取I0=∫01ex-1dx=ex-1|01 =1-e-1=0.6321 逐次递推得I1,I2,…,I7.