第1课 数值分析NumericalAnalysis.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 解法之一：将 作Taylor展开后再积分 S4 R4 /* Remainder */ 取 则 称为截断误差 /* Truncation Error */ | 舍入误差 /* Roundoff Error */ | = 0.747… … 由截去部分 /* excluded terms */ 引起 由留下部分 /* included terms */ 引起 积分不等式 2.误差的传播与积累 /* Spread Accumulation */ 例：蝴蝶效应 —— 纽约的一只蝴蝶翅膀一拍，风和日丽的桂林就下起狂风暴雨了？！ NY GL 以上是一个病态问题 /* ill-posed problem*/ 病态：自变量的微小变化就可能引起函数值的巨大变化的情况；反之称良态。 例：计算 ?公式一： 注意此公式精确成立 记为 则初始误差 ? ?? ? ! ! ! What happened?! 考察第n步的误差 造成这种情况的是不稳定的算法 /* unstable algorithm */。 迅速积累，舍入误差递增。 可见初始的小扰动 ? 公式二 注意此公式与公式一 在理论上等价。 方法：先估计一个IN ,再反推要求的In ( n N )。 可取 取 We just got lucky? 考察反推一步的误差： 以此类推，对 n N 有： 误差逐步递减, 这样的算法称为稳定的算法 /* stable algorithm */ 在我们今后的讨论中，误差将不可回避， 算法的稳定性是一个非常重要的问题。 3.数值算法的稳定性 /* Algorithm Stability */ 对于某个给定的算法，如果输入数据的误差随着算法的运算过程不断增长而得不到控制，那么我们 就说该算法是数值不稳定的；否则就是数值稳定的。 误差线性增长，算法稳定，不可避免； 误差指数增长，算法不稳定，必须避免。 假设输入数据的误差为 ，经 次运算后，结果的