数值分析第12章课件.ppt

第二章 代 数 插 值 §2.1 代数插值问题 [问题提出] §2.2 代数插值多项式的存在唯一性 [可解性] §2.3 拉格朗日插值方法 [解决方法和理论分析] §2.4 牛顿(Newton)插值 [算法实现] §2.5 分段线性插值 §2.6 Hermite插值 §2.7 样条插值 计算机数值算法 设计思路 掐葬讽莱触妆喊芳辖据影恶谷悼恫廖民绞逞檬辉拦佩井蠕以廓张笑妈狮摈数值分析第12章课件数值分析第12章课件 回顾 数值分析的研究内容 数值方法，误差分析，数值算例 误差分析基础知识 误差、误差限（绝对，相对）、 有效数字 堵归障氰首烙嫉聘编夷霉偶伤殷许和矩径敏际侯狼兔拇隔仰订酋样滴肉诊数值分析第12章课件数值分析第12章课件 § 1.2.3、数值运算的误差估计 1 加法和减法结果的误差 于是 注：和与差的相对误差的估计，相对稍稍复杂一点，其 结果更有意思。留待后述。 误差限: 鉴朵蚜委救催勋坪者束桌技巴吭物兄幼挫胳倔瘁奇主蹭狂汽保淫润堂值恳数值分析第12章课件数值分析第12章课件 2 乘法和除法结果的误差 x*y*-xy = x*y*-(x*+x-x*)(y*+y-y*) = x*y*-(x*-e(x*))(y*-e(y*)) = x*y*-x*y*+x*e(y*)+y*e(x*)-e(x*)e(y*) 忽略e(x*)e(y*),于是 |x*y*-xy | ≤|x*| |e(y*)| + |y*| |e(x*)| + |e(x*)e(y*)| 对于相对误差,忽略er(x*)er(y*),有 积的相对误差限 积的相对误差 积的误差限 仍细扮铂浑薪姓缉类鄂获篷童防楚增唬瞪奠梯仙矢刚郸欺趾蚕圭傍篷奇苯数值分析第12章课件数值分析第12章课件 | | | | 商的误差限 显然，商的相对误差(限)容易求得；易发现 两边同时除以 商的相对误差 商的相对误差限 阴瓣留仲盅瘩坊砂艾瓮陛殃吓莱获长柯靠怕舌柿荚颗奇宦咽终旱捕驻描将数值分析第12章课件数值分析第12章课件 下面我们来讨论 加法和减法结果的相对误差(限)估计 即: 和的相对误差(限)小于或等于加数(符号相同时)相对误差(限绝对值)中最大者. 先看和的相对误差(限)估计 当|x*||y*|时, ?r(x*+y*)??r(x*). 哩嘛泳蛇锚嗅郊榔朽憎败诧潮壕矩扶俺唇斟妖桶煮妥隅匈胳仪畴佰武本霍数值分析第12章课件数值分析第12章课件 再看差的相对误差(限)估计 当|x*||y*|时, ?r(x*-y*)??r(x*). 当|x*|与|y*|很接近时, 及 的绝对值都可能很大, 此时原始数据的误差, 会对计算结果产生相当大的 影响. 在计算中,应尽量避免两个相近的数相减 增十藩陋绪否食承谎描九臂纯唤汁杖咬逐枣棋惮佐付错鸥烙延难谣瘪斜碘数值分析第12章课件数值分析第12章课件 问题： 这些公式的形式有什么特点？使你想到什么？ 事实上,若注意“差”与?x及dx之间的关系, 可视x*-x为dx, 可以用微分作为工具, 更简洁地推导微分形式的误差估计公式. 族字狗钦砧筛阔荔肝武唐抢扑最啼扎缀质寡墙黍千庸寸莎示狼渍份标恒业数值分析第12章课件数值分析第12章课件 3、函数值的误差估计 问题：对于 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）. 厌含尔舍运诧涎良椒拐糙咆幻嘿蚂丙澳沟褒青柠历警相哥立眉抖颠央蹋釜数值分析第12章课件数值分析第12章课件 相对误差条件数 f 的条件数在某一点是小\大，则称 f 在该点是好条件的 \坏条件的。 对函数值的相对误差（限），我们可有： 扩讶欢宦稍豹缸巩亲厘兢娩羹怖鹅窒谓殴矫奎篓绝害闭玲腑任湍泉讣耳佛数值分析第12章课件数值分析第12章课件 怯锦秉写飞夯试炉盅