数值与算法_第1讲_引论概览.pptVIP

* * 2、函数误差 当自变量有误差时计算函数值也产生误差，可以利用函数的泰勒展开式进行估计。 * * * * 一、病态问题与条件数 1、病态问题：对一个数值问题本身如果输入数据有微小扰动（即误差），引起输出数据（即问题解）相对误差很大，就是病态问题。 误差估计定性分析与避免误差危害 * * 二、算法的稳定性 1、用一个算法进行计算，由于初始数据误差在计算中传播使计算结果误差增长很快使得数值不稳定。 * * 计算结果： n 法一 (A) 法二 (B) 0 1 2 3 4 5 6 7 8 9 0.6321 0.3679 0.2642 0.2074 0.1704 0.1480 0.1120 0.2160 -0.7280 7.552 0.6321 0.3679 0.2643 0.2073 0.1708 0.1455 0.1268 0.1121 0.1035 0.0684 例(对舍入误差不敏感算法)：有长度为100的数组，每个数组元素都是仅有两位有效数字的实数，并且假设执行两位精度的浮点运算，要求这些元素之和。 常规的算法是按数组存储顺序依次累加数值，但它可能导致结果有很大误差。例如数组中的数依次为:1.0, 0.01, 0.01，……，0.01(共99个相同的0.01),由于浮点计算的舍入误差，按常规算法将得到结构为1.0,大大偏离准确率。 怎样提高准确率？ * * 三、避免误差危害的若干原则 1、要避免除数绝对值远远小于被除数绝对值的除法。 用绝对值小的数作除数舍入误差会增大，如计算 x/y， 若0|y||x|，则可能对计算结果带来严重影响，应尽量避免。 2、要避免两相近数相减 在数值中两相近数相减有效数字会严重损失。 例如，x=532.65，y=532.52都具有五位有效数字，但 x - y=0.13只有两位有效数字。通过改变算法可以避免两相近 数相减。 * * 3、要防止“大数”吃掉小数 数值运算中参加运算的数有时数量级相差很大，而计算机位数有限，如不注意运算次序就可能出现大数“吃掉”小数的现象，影响计算结果的可靠性。 如用六位浮点数计算某市的工业总产值，原始数据是各企业的工业产值，当加法进行到一定程度，部分和超过100亿元 （0.1×1011），再加产值不足10万元的小企业产值，将再也加不进去。而这部分企业可能为数不少，合计产值相当大.这种情况应将小数先分别加成大数，然后相加，结果才比较正确。这个例子告诉我们，在计算机数系中，加法的交换律和结合律可能不成立，这是在大规模数据处理时应注意的问题。 * 4、注意简化计算步骤，减少运算次数 减少算术运算的次数不但可计算机的计算时间，还能减少误差的积累效应。使参加运算的数字精度应尽量保持一致，否则那些较高精度的量的精度没有太大意义。 教材与参考书 * ? Numerical Analysis:Mathematics of Scientific Computing (Third Edition) 数值分析 （英文版 第3版 ） David Kincaid Ward Cheney（机械工业出版社) ? Numerical Analysis (Seventh Edition) 数值分析 （第七版 影印版） Richard L. Burden J. Douglas Faires （高等教育出版社） 工程数值分析 王立秋等编著（山东大学出版社） 工程数值分析题解 现代数值分析 蔺小林、蒋耀林 编著（国防工业出版社） 参考书 * 先修课程 微积分 线性代数 会用Matlab * 学习方法 1.注意掌握各种方法的基本原理 2.注意各种方法的构造手法 3.重视各种方法的误差分析 4.做一定量的习题 5.注意与实际问题相联系 * * 主要教学内容 一.绪论 (数值分析与科学计算引论) 误差的基本概念(有效数字)、数值运算的误差限、问题的敏感性与条件数、数值算法的稳定性、“抵消”现象、减小舍入误差的建议 二.非线性方程解法 三.线性方程组的直接解法 四.线性方程组的迭代解法 五.矩阵特征值计算 六-1.函数逼近与线性最小二乘 六-2.函数插值 七.数值积分与微分 八.常微分方程初值问题 附加.Matlab数值计算与应用 * * * * SIAM:Society for Industry and Applied Mathematics工业与应用数学学会成立于二十世纪五十年代前期，是一个促进应用和计算数学的研究、发展、应用为目的的协会 * SIAM:Society for Industr