NA1绪论.PPTVIP

00000.02 1 00.008 00.004096 算　式 例 利用下列恒等式计算　　　　　，结果见下表 ： 2.3 Spread Estimation for Error 问题：对于 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*)|为放大因子 /* amplification factor */ 或 绝对条件数 /* absolute condition number */. 2.3 Spread Estimation for Error 相对误差条件数 /* relative condition number*/ f 的条件数在某一点是小\大，则称 f 在该点是好条件的 /* well-conditioned */ \坏条件的 /* ill-conditioned */。 2.3 Spread Estimation for Error 1. Introduction 例(P.12 4.) 测得某地纬度φ = 45o56’23” ，测量误差不超过0.5s。问计算sinφ的误差界是多少，算出sinφ=0.718608590有几位数字准确? 解： 设 y=sinφ, e*(y)= y’(φ) e*(φ) = cosφe*(φ) φ=(45+56/60+23/3600)o =45 =0(rad) e*(φ) =0.5”=010-5　(rad) e*(y)=1.6857330×10-60.17×10-50.5×10-5. sinφ=0.718608590具有 有效数字。 5位 ?多元函数 y = f (x1,x2,…,xn)的误差估计 公式中系数绝对值 便是绝对条件数和相对条件数。 若用近似值 x* 取代 x，则函数的绝对误差 2.3 Spread Estimation for Error 相对误差： 2.3 Spread Estimation for Error 例 (P.13 6.)证明下列命题: (1)设函数y=f(x)可导，f(a)=0 ，则在x≈a时计算f(x)的值，按相对误差是病态问题。 (2) 设a是方程p(x)=c0xn+ c1xn-1+…+cn-1x+cn的根，则当接近重根时(即p’(a)≈ 0时)按照绝对误差求根a的问题是病态问题。 证(1)： 因为函数y=f(x)可导，则y的相对误差条件数为 由于f(a)=0，当x≈a时有 故是病态问题。 2.3 Spread Estimation for Error 例 (P.13 6.)证明下列命题: (2) 设a是方程p(x)=c0xn+ c1xn-1+…+cn-1x+cn的根，则当接近重根时(即p’(a)≈ 0时)按照绝对误差求根a的问题是病态问题。 证(2)：因为方程p(x)=0的根取决于系数ci, i=0,1,…,n，所以有 当 时， 是坏条件数。 因此求根的条件数是 根据隐函数求导法有 设参数x1，x2的绝对误差和相对误差分别是Δx1,Δx2,δx1,δx2 ?算术运算 y = f (x1,x2)的误差估计 １．加减运算. 2．乘法运算. 3．除法运算. 4．幂运算. 2.3 Spread Estimation for Error 例 　　考察算式　　　　，　　　　，　　　　， 　　　　　　的计算精度。 解：设 绝对误差限最小 绝对误差限最大 2.3 Spread Estimation for Error 2.4 几点注意事项 /* Remarks */ 2.4.1 避免相近二数相减 取Δx=0.1，有： 例 用中心差商公式 求 在x=2处的导数值。（准确值f’(x)x=2=0.353553…） 取Δx=0.0001，有： 解:当x=2时有： 有效位数损失殆尽 对策： ? 几种经验性避免方法： 当