数值分析课件 第1章.ppt

数 值 分 析 李庆扬 王能超 易大义 编 第1章 绪 论 1.1 数值分析研究对象与特点 1.1 数值分析研究对象与特点 1.2 数值计算的误差 1.2.2 误差与有效数字 1.2.3 数值运算的误差估计 1.3 误差定性分析与避免误差危害 1.3.1 病态问题与条件数 1.3.2 算法的数值稳定性 1.3.3 避免误差危害的若干原则 设 是一元函数， 的近似值为 ，以 近 似 ，其误差界记作 ， 一般情况下，当自变量有误差时函数值也产生误差， 取绝对值得 其误差限可利用函数的泰勒展开式进行估计. 利用泰勒展开 当 为多元函数，如计算 时, 的近似值为 ， 则 的近似值为 于是由泰勒展开， 函数值 的误差 为 假定 与 的比值不太大，可忽略 的高阶项， 于是可得计算函数的误差限 如果 于是误差限 （2.3） 而 的相对误差限为 （2.4） 已测得某场地长 的值为 ， 宽 的值 为 ， 已知 . 试求面积 的绝对误差限与相对误差限. 因 知 例4 解 由 其中 而 于是绝对误差限 相对误差限 一个工程或科学计算问题往往要运算千万次，由于每 步运算都有误差，如果每步都做误差分析是不可能的，也 不科学. 因为误差积累有正有负，绝对值有大有小，都按最坏 情况估计误差限得到的结果比实际误差大得多，这种保守 的误差估计不反映实际误差积累. 考虑到误差分布的随机性，有人用概率统计方法，将 数据和运算中的舍入误差视为适合某种分布的随机变量， 20世纪60年代以后对舍入误差分析提出了一些新方法， 然后确定计算结果的误差分布，这样得到的误差估计更接近实际， 这种方法称为概率分析法. 较重要的有向后误差分析法和区间分析法两种. 1. 向后误差分析法是把新算出的量由某个公式表达， 若 的摄动为 ， 它仅含基本算术运算，如假定 是前面已算出的量 或原始数据，新算出量 使得由浮点运算得出结果为 则可根据 的界， 的界，威克逊(Wilkinson)将这种方法应用于数值代数(矩阵运算)的误差分析，取得较好效果. 由摄动理论估计最后舍入误差 2. 区间分析法是把参加运算的数 都看成区 间量 ，根据区间运算规则求得最后结果的近似 值及误差限. 例如， 的近似数为 ， 则 由于 若计算 (*为运算符号)， 而 则为误差限. 则 为所求近似值， 由 对一个数值问题本身, 如果输入数据有微小扰动（即误 差），引起输出数据（即问题解）相对误差很大，这就是病 态问题. 例如计算函数值 时， 若 有扰动 ，其 相对误差为 ， 函数值 的相对误差为 （3.1） 称为计算函数值问题的条件数. 相对误差比值 自变量相对误差一般不会太大，如果条件数 很大， 将引起函数值相对误差很大，出现这种情况的问题就是病态问题. 例如， ， 它表示相对误差可能放大 倍. 如 ， 有 ， 自变量相对误差为 ， 函数值相对误差为 ， 一般情况下,条件数 就认为是病态， 越大 病态越严重. 则有 若取 这时问题可以认为是病态的. 其他计算问题也要分析是否病态. 例如解线性方程组，如果输入数据有微小误差引起解 的巨大误差，就认为是病态方程组，第5章将用矩阵的条件 数来分析这种现象. 用一个算法进行计算，如果初始数据误差在计算中传播 使计算结果的误差增长很快，这个算法就是数值不稳定的. 计算 并估计误差. 由分部积分可得计算 的递推公式 若计算出 ， 代入（3.2），可逐次求出 的值. （3.2） 例5 而要算出 就要先计算 . 并取 ， 则得 ， 计算过程中小数点后第5位的数字按