计算方法误差度量.ppt

? 初值误差传播：假设每一步都是准确 计算，即不考虑截断误差和由运算进 一步引入的舍入误差， 仅介绍初始数 据的误差传播规律。 – 研究方法： ? 泰勒（ Taylor ）方法 – n 元函数 复习泰勒公式 ? ? ? ? ? ? ? ? ? ) )( ( ) )( ( ! 1 1 ) ( ) ( * * * 1 1 * 1 * n n n x x p f x x p f p f p f ? 记点 ) , , , ( * * 2 * 1 n x x x ? 为 p * , 点 ) , , , ( 2 1 n x x x ? 为 p ， n 元泰勒公式 ： ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 2 1 1 1 2 2 2 1 1 2 2 21 1 1 1 2 1 1 11 2 1 ) x x )( p ( f ) x x )( x x )( p ( f ) x x )( x x )( p ( f ) x x )( x x )( p ( f ) x x )( x x )( p ( f ) x x )( p ( f ! * n n * nn * * n n * n * n n * * n * * * * n n * * n * * 泰勒公式分析初值误差传播 设 n 元可微函数 ) , , , ( n 2 1 x x x f y ? ? 中的自变量 x 1 、 x 2 、 … 、 x n 是相互独立的。 用自变量的近似值进行准确计算，得 y * = ) , , , ( * * * n 2 1 x x x f ? 。 当 * 1 x 、 * 2 x 、…、 * n x 很好地近似了相应真值时，利用多元函 数一阶 Taylor 公式求得 y * 的绝对误差： ? ? ? ? ? ? n 1 i * i i * n * 1 i * * ) )( , , ( ) ( x x x x f y y y e ? ? ? ? n 1 i * i * n * 1 i ) ( ) , , ( x e x x f ? ? 相对误差： ? ? ? ? n 1 i * i * i * n * 1 i * * i * * * ) ( ) , , ( ) ( ) ( x x e x x f y x y y e y e r ? ? ? ? n 1 i * i * n * 1 i * * i ) ( ) , , ( x e x x f y x r ? ? 进而得到如下绝对误差限和相对误差限传播关系： ? ? ? ? n 1 i * i * n * 1 i * ) ( ) , , ( ) ( x x x f y ? ? ? ? ? ? ? n 1 i * i * n * 1 i * * i * ) ( ) , , ( ) ( x x x f y x y r r ? ? ? 对于一元函数，有如下初值误差传播近 似计算公式： ) ( ) ( ) ( * * * x e x f y e ? ? ) ( ) ( ) ( * * * * * x e y x x f y e r r ? ? ? 二元函数算术运算误差传播规律 ) ( ) ( ) ( * 2 * 1 * 2 * 1 x x x x ? ? ? ? ? ? ) ( ) ( ) ( * 2 * 1 * 1 * 2 * 2 * 1 x x x x x x ? ? ? ? ? ) 0 ( ) ( ) ( ) ( * 2 2 * 2 * 2 * 1 * 1 * 2 * 2 * 1 ? ? ? x x x x x x x x ? ? ? ? 绝对误差限 ? 相对误差限 )} ( ), ( { max ) ( * 2 * 1 * 2 * 1 x x x x r r r ? ? ? ? ? ) ( * * 0 2 1 ? x x ) ( ) ( ) ( * 2 * 1 * 2 * 1 x x x x r r r ? ? ? ? ? ) ( * * 0 2 1 ? x x ) ( ) ( ) ( * 2 * 1 * 2 * 1 x x x x r r r ? ? ? ? ? ) ( * * 0 2 1 ? x x ? 尽量避免相近的数相减 ? 例 x=52.127 x * =52.129 四位有效数字 y=52.123 y * =52.121 四位有效数字 A=x-y=0.004 A * =x * -y * =0.008 零位有效数字 ? 结论：避免相近