数值分析李庆杨-第一章 引论.ppt

数值分析 五、教学参考书 《数值计算引论》 白峰杉 高等教育出版社 《科学和工程计算基础》 施妙根等 清华大学出版社 《数值分析》，易大义等编，浙江科学技术出版社 《数值方法教程》，刘钦圣等编 , 冶金出版社,1998 《计算方法》，秦林祥等编, 兵器工业出版社,1992 《数值分析基础》，关治等编 , 高教出版社,1998 在用数值方法解题过程中可能产生的误差归纳起来有如下几类： 1. 模型误差 2. 观测误差 3. 截断误差 4. 舍入误差 用数学方法解决一个具体的实际问题，首先要建立数学模型，这就要对实际问题进行抽象、简化，因而数学模型本身总含有误差，这种误差叫做模型误差 数学模型是指那些利用数学语言模拟现实而建立起来的有关量的描述 数学模型的准确解与实际问题的真解不同 在数学模型中通常包含各种各样的参变量，如温度、长度、电压等，这些参数往往是通过观测得到的，因此也带来了误差，这种误差叫观测误差 数学模型中的参数和原始数据，是由观测和试验得到的 由于测量工具的精度、观测方法或客观条件的限制,使数据含有测量误差,这类误差叫做观测误差或数据误差 根据实际情况可以得到误差上下界 数值方法中需要了解观测误差,以便选择合理的数值方法与之适应 精确公式用近似公式代替时,所产生的误差叫截断误差 例如, 函数f(x)用泰勒(Taylor)多项式 在数值计算中只能对有限位字长的数值进行运算 需要对参数、中间结果、最终结果作有限位字长的处理工作，这种处理工作称作舍入处理 用有限位数字代替精确数，这种误差叫做舍入误差，是数值计算中必须考虑的一类误差 算法的数值稳定性 算法优劣的标准 从截断误差观点看,算法必须是截断误差小,收敛敛速要快。即运算量小,机器用时少. 从舍入误差观点看,舍入误差在计算过程中要能控制,即算法的数值要稳定. 从实现算法的观点看,算法的逻辑结构不宜太复杂，便于程序编制和上机实现. 设计算法时应遵循的原则 要有数值要稳定性,即能控制误差的传播. 避免大数吃小数,即两数相加时,防止较小的数加不到较大的数上. 避免两相近的数相减,以免有效数字的大量丢失. 避免分母很小(或乘法因子很大),以免产生溢出. 例 求二次方程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位。 三、避免误差危害的若干原则 除了分清问题是否病态和算法是否数值稳定外,还要考虑避免误差危害和防止有效数字损失的如下原则. 1.避免‘大数’除以‘小数’ 例8 仿计算机，采用3位十进制，用消元法求解方程组 解： 错.为什么,怎么办？ 减少运算误差原则 2、两个相近的数相减，会严重损失有效数字 例如x =1958.75，y =1958.32都具有五位 有效数字，但x-y=0.43只有两位有效数字 通常采用的方法是改变计算公式,例如当与 很接近时,由于 用右端代替左端公式计算,有效数字就不会损失 减少运算误差原则 当x很大时可作相应的变换 则用右端来代替左端。 减少运算误差若干原则 当x接近0时 一般情况，当f(x)≈f(x*)时，可用泰勒展开 取右端的有限项近似左端。 如果计算公式不能改变，则可采用增加有效位数的方法保证精度 例11 仿计算机在3位十进制下， ３、防止‘大数’吃‘小数’ 这里分子的误差被扩大104倍,再如 若将分母变为0.0011,即分母只有0.0001的变化时,计算结果却有了很大变化 减少运算误差若干原则 解: 分子分母分别计算后相除(取9位小数) A=0.0005*0.0143*0.0012=00.0012