数值分析简介.pptVIP

理学院 信计教研室 例15 : 计算 y = ln x。若 x ? 20，则取 x 的几位有效数字可保证 y 的相对误差 0.1% ? 解：设截取 n 位有效数字后得 x* ? x，则 估计 x 和 y 的相对误差上限满足近似关系 不知道怎么办啊？ x 可能是20.#，也可能是19.#，取最坏情况，即a1 = 1。 ? n ? 4 例：计算 ，取 4 位有效，即 , 则相对误差 §4 Error Estimation for Functions 理学院 信计教研室 特例----算术运算的误差和绝对误差限分析: §4 Error Estimation for Functions 理学院 信计教研室 §4 Error Estimation for Functions 理学院 信计教研室 §4 Error Estimation for Functions 理学院 信计教研室 特例----算术运算的相对误差和相对误差限分析: §4 Error Estimation for Functions 理学院 信计教研室 §4 Error Estimation for Functions 理学院 信计教研室 §4 Error Estimation for Functions 理学院 信计教研室 §4 Error Estimation for Functions 理学院 信计教研室 1.3 算法的数值稳定性 用一个算法进行计算，如果初始数据误差在计算中传播 使计算结果的误差增长很快，这个算法就是数值不稳定的. 计算 并估计误差. 由分部积分可得计算 的递推公式 （3.2） 例 问题：由递推公式（3.2）可以构造几种计算积分序列{In}(n=1,2, …,N)的算法？假设计算中只有初始值带有误差，计算过程中不 产生新的误差，请判断哪种方法能给出更精确的结果。 理学院 信计教研室 例 计算 ( n =0,1,···, 20 ) 5/18 理学院 信计教研室 例6：计算 ? 公式一： 注意此公式精确成立 记为 则初始误差 ? ?? ? ! ! ! What happened?! §2 Introduction of error 理学院 信计教研室 考察第n步的误差 我们有责任改变。 造成这种情况的是不稳定的算法 /* unstable algorithm */ 迅速积累，误差呈递增走势。 可见初始的小扰动 ? 公式二： 注意此公式与公式一 在理论上等价。 方法：先估计一个IN ,再反推要求的In ( n N )。 可取 §2 Introduction of error 理学院 信计教研室 取 We just got lucky? §2 Introduction of error 理学院 信计教研室 考察反推一步的误差： 以此类推，对 n N 有： 误差逐步递减, 这样的算法称为稳定的算法 /* stable algorithm */ 在我们今后的讨论中，误差将不可回避， 算法的稳定性会是一个非常重要的话题。 §2 Introduction of error 理学院 信计教研室 传播与积累 /* Spread Accumulation */ 例：蝴蝶效应 —— 纽约的一只蝴蝶翅膀一拍，风和日丽的北京就刮起台风来了？！ NY BJ 以上是一个病态问题 /* ill-posed problem*/ 关于本身是病态的问题，我们还是留给数学家去头痛吧！ §2 Introduction of error 理学院 信计教研室 1.2 误差 /* Error */ §1 误差的背景介绍 /* Introduction */ 1. 来源与分类 /* Source Classification */ 从实际问题中抽象出数学模型 —— 模型误差 /* Modeling Error */ 通过测量得到模型中参数的值 —— 观测误差 /* Measurement Error */ 求近似解 —— 方法误差 (截断误差 Truncation Error） 机器字长有限 —— 舍入误差 /* Roundoff Error */ 理学院 信计教研室 §1.2.4 误