第一章计算方法素材.pptVIP

板书说明有效数字， 定义结束看习题3，例1结束习题2 有效数字与误差限之间的关系 有效数字刻画的是精确度 * 习题5 * 板书证明 * * 桥梁 习题4 作业 习题4（2），5（2) * 例题精选1.2 p19 * * 习题六（1）（3） * 周三 习题六（2） 条件数太大则是病态问题 周五下午 * * * 习题十 * 习题7 * * 习题9 * 问题：对于 y = f (x)，若用 x* 取代 x，将对y 产生什么影响？ 分析：e*(y) = f (x*) ? f (x) e*(x) = x* ? x Mean Value Theorem = f ’(? )(x* ? x) x* 与 x 非常接近时，可认为 f ’(? ) ? f ’(x*) ，则有： |e*(y)| ? | f ’(x*)|·|e*(x)| 即：x*产生的误差经过 f 作用后被放大/缩小了| f ’(x*)|倍。故称| f ’(x*)|为放大因子或绝对条件数。 一 函数的误差估计 相对误差条件数 f 的条件数在某一点是小\大，则称 f 在该点是好条件的 \坏条件的. 二、算法的数值稳定性 考虑初始数据误差在计算中的传播问题. P4，例题1.1 控制递推公式中误差的传播 对于一个数学问题的求解往往有多种数值方法在选择数值方法时，要注意所用的数值方法不应将计算过程中难以避免的误差放大的较快，造成计算结果完全失真。 例 计算积分 并估计误差 解 容易得到递推公式 即 为 则准确的理论递推式 实际运算的递推式 两式相减有 这就是说,若 与 的误差为 = - ,即 ，则误差的递推规律为 于是 计算 时的误差被扩大了 倍,显然算法是数值不稳定的。 如果将递推公式 变换一种形式 准确的理论递推式 实际运算的递推式 从而有 即 于是有 则这个算法的误差传递规律为 即每计算一步的误差的绝对值是上一步的十分之一，误差的传播逐步缩小，得到很好的控制，这个算法是数值稳定的 算法的数值稳定性 算法优劣的标准 从截断误差观点看,算法必须是截断误差小,收敛敛速要快。即运算量小,机器用时少. 从舍入误差观点看,舍入误差在计算过程中要能控制,即算法的数值要稳定. 从实现算法的观点看,算法的逻辑结构不宜太复杂，便于程序编制和上机实现. 设计算法时应遵循的原则 要有数值要稳定性,即能控制误差的传播. 避免大数吃小数,即两数相加时,防止较小的数加不到较大的数上. 避免两相近的数相减,以免有效数字的大量丢失. 避免分母很小(或乘法因子很大),以免产生溢出. 三、避免误差危害的若干原则 除了分清问题是否病态和算法是否数值稳定外,还要考虑避免误差危害和防止有效数字损失的如下原则. 1.避免‘大数’除以‘小数’ 例8 仿计算机，采用3位十进制，用消元法求解方程组 解： 错.为什么 减少运算误差原则 2、两个相近的数相减，会严重损失有效数字 例如x =1958.75，y =1958.32都具有五位 有效数字，但x-y=0.43只有两位有效数字 通常采用的方法是改变计算公式,例如当与 很接近时,由于 用右端代替左端公式计算,有效数字就不会损失 减少运算误差原则 当x很大时可作相应的变换 则用右端来代替左端。 减少运算误差若干原则 当x接近0时 一般情况，当f(x)≈f(x*)时，可用泰勒展开 取右端的有限项近似左端。 如果计算公式不能改变，则可采用增加有效位数的方法保证精度 例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、绝对值太小的数不宜做除数 当分母为两个相近数相减时,会丧失有效数字 这里分子的误差被扩大104倍,再如 若将分母变为0.0011,即分母只有0.0001的变化时,计算结果却有了很大变化 减少运算误差若干原则 例 计算 解: 分子分母分别计算后相除