數值分析Numerical Analysis.pptVIP

数 值 分 析Numerical Analysis 教材 (Text Book) 数值分析 李庆扬 王能超 易大义 编 （清华大学出版社） 参考书目 (Reference) ? Numerical Analysis (Seventh Edition) 数值分析 （第七版 影印版） Richard L. Burden J. Douglas Faires （高等教育出版社） ? Numerical Analysis 数值分析 Rainer Kress （Springer-Verlag 世界图书出版公司） ? 数值分析 李红 （华中科技大学出版社） ? 计算方法简明教程 王能超 （高等教育出版社） ? MATLAB数值分析 周品 （机械工业出版社） 第一章 数值分析与科学计算引论 /* Chapter 1 Introduction */ 什么是数值分析 数值分析（或计算方法）：用计算机求解数学问题的数值计算方法和理论与软件实现 实际问题数学模型计算方法算法程序 算法（数值分析的根本任务）：由基本运算及运算顺序所构成的完整的解决问题的步骤 中国古代数学-算术 中国古代将数学称为——“算术” “凡算之法，先识其位。一纵十横，百立千僵；千十相望，万百相当。” ——《孙子算经》 “算法之术，是用智矣。” ——《周髀算经》 现代数值分析 数值逼近 数值微分 数值积分 数值代数 非线性方程的数值解 微分方程数值解 研究方法 理论分析 算法分析 误差分析 稳定性分析 收敛性分析 好的算法： 面向计算机 有可靠的理论分析 有好的计算复杂性 有数值实验验证 误差分析 从实际问题中抽象出数学模型 —— 模型误差 /* Modeling Error */ 通过测量得到模型中参数的值 —— 观测误差 /* Measurement Error */ 求近似解 ——方法误差 (截断误差 /* Truncation Error */ ) 计算机字长有限 —— 舍入误差 /* Roundoff Error */ 1.误差来源与分类 /* Source Classification */ 1 1 / e = 0.747… … 由截去部分 /* excluded terms */ 引起 由留下部分 /* included terms */ 引起 积分不等式 2.误差的传播与积累 /* Spread Accumulation */ 例：蝴蝶效应 —— 纽约的一只蝴蝶翅膀一拍，风和日丽的桂林就下起狂风暴雨了？！ NY GL 以上是一个病态问题 /* ill-posed problem*/ 病态：自变量的微小变化就可能引起函数值的巨大变化的情况；反之称良态。 注意此公式精确成立 ? ?? ? ! ! ! What happened?! 注意此公式与公式一 在理论上等价。 方法：先估计一个IN ,再反推要求的In ( n N )。 考察反推一步的误差： 以此类推，对 n N 有： 误差逐步递减, 这样的算法称为稳定的算法 /* stable algorithm */ 在我们今后的讨论中，误差将不可回避， 算法的稳定性是一个非常重要的问题。 3.数值算法的稳定性 /* Algorithm Stability */ 对于某个给定的算法，如果输入数据的误差随着算法的运算过程不断增长而得不到控制，那么我们 就说该算法是数值不稳定的；否则就是数值稳定的。 误差线性增长，算法稳定，不可避免； 误差指数增长，算法不稳定，必须避免。 问题：对于 y = f (x)，若用 x* 取代 x，将对y 产生什么影响？ 分析：e*(y) = f (x) ? f (x*) e*(x) = x ? x* = 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 */. 4. 函数的误差估计与条件数 相对误差条件数 /* r