大学课件 数值分析Numerical Analysis.pptVIP

* 计算函数值的误差估计 当 为多元函数时，我们有如下结果： 　其中　在 与 之间，如果 ， 和 的比值不太大，可得： * 当 为多元函数时，我们有如下结果： 因此： * 的相对误差限为： 因此： * 的相对误差限为： * 例题 已测得某场地长 的值为 ，宽 的值为 ，已知 ，试求面积 的绝对误差限和相对误差限。 解：因 ，所以 ， ，故： * 例题 其中 而 ， 解：因 ，所以 ， ，故： * 其中 而 ， 于是绝对误差限为： 相对误差限为 * 其中 而 ， 于是绝对误差限为： 相对误差限为 * 1.3 误差定性分析与避免误差危害 前面讨论的误差限的方法是最坏情况 对于千万次的运算可用概率统计的方法 20世纪60年代后提出 向后误差分析法 区间分析法 目前尚无有效方法对误差做出定量分析 本节讨论： 病态问题与条件数 算法的数值稳定性 避免误差危害的若干原则 * 1.3 误差定性分析与避免误差危害 输入数据的微小变动导致输出数据的较大误差，就被称为病态问题。 衡量是否病态的标准：条件数 对于函数值计算问题，条件数定义为： 目前尚无有效方法对误差做出定量分析 本节讨论： 病态问题与条件数 算法的数值稳定性 避免误差危害的若干原则 * 1.3.1 病态问题与条件数 输入数据的微小变动导致输出数据的较大误差，就被称为病态问题。 衡量是否病态的标准：条件数 对于函数值计算问题，条件数定义为： 不同的问题，条件数具体定义不同。 * 条件数的例子 设 ，则 ，即对于该函数，误差会被放大 倍。 一般情况下，条件数大于10，就认为问题病态。 不同的问题，条件数具体定义不同。 * 条件数的例子 设 ，则 ，即对于该函数，误差会被放大 倍。 一般情况下，条件数大于10，就认为问题病态。 * 1.3.2 算法的数值稳定性 由于初始数据的误差在计算过程中的传播使得计算结果的误差很快增长，就称该问题是数值不稳定的。 计算 并估计误差 由分部积分公式可得： * 1.3.2 算法的数值稳定性 按公式 可以计算出下表 由分部积分公式可得： * 按公式 可以计算出下表 n n 0 0.6321 5 0.1480 1 0.3679 6 0.1120 2 0.2642 7 0.2160 3 0.2074 8 -0.7180 4 0.1704 9 7.552 * 现在改变方法，原积分有估计式： n n 0 0.6321 5 0.1480 1 0.3679 6 0.1120 2 0.2642 7 0.2160 3 0.2074 8 -0.7180 4 0.1704 9 7.552 * 现在改变方法，原积分有估计式： 取n=9，可得： 我们取 接着用公式 得到下面的计算结果 * 接着用公式 得到下面的计算结果 n n 0 0.6321 5 0.1455 1 0.3679 6 0.1268 2 0.2643 7 0.1121 3 0.2073 8 0.1035 4 0.1708 9 0.0684 * 可见第二种方法计算稳定。 n n 0 0.6321 5 0.1455 1 0.3679 6 0.1268 2 0.2643 7 0.1121 3 0.2073 8 0.1035 4 0.1708 9 0.0684 * 误差分析 对于第一种计算方法 设 ，则有 因此，误差是逐次放大的。 对于第二种计算方法 设 ，则有 因此，误差是逐次缩小的。 * 数值稳定的定义 一个算法如果输入数据有误差，而在计算过程中舍入误差不增长，则称此算法是数值稳定的，否则称此算法为不稳定的。 对于第二种计算方法 设 ，则有 因此，