数值第一章基础知识解读.ppt

*/76 郑州大学2014-2015学年硕士研究生课程 数值分析 Numerical Analysis §1.5 数值算法的稳定性 一个算法，如果在执行它的过程中，舍入误差可以在一定条件下得到控制，从而不影响算法得到可靠的结果，则称它是数值稳定的。否则称之是不稳定的。 如果对输入数值的微小误差，经过某个算法的运算，输出结果产生了较大变化，则称算法是病态的；否则称算法是良性的。 */76 郑州大学2014-2015学年硕士研究生课程 数值分析 Numerical Analysis §1.5 数值算法的稳定性 例1.5.1建立一个算法来计算积分 解：注意到 两边在区间[0,1]上积分可得 从而可得迭代递推格式 */76 郑州大学2014-2015学年硕士研究生课程 数值分析 Numerical Analysis §1.5 数值算法的稳定性 按照以上关系编程计算可得 n= 0 积分值I(n)= 1.000500e-003 n= 1 积分值I(n)= 5.001673e-004 n= 2 积分值I(n)= 3.328343e-004 n= 3 积分值I(n)= 8.318346e-004 n= 4 积分值I(n)= -5.810028e-001 n= 5 积分值I(n)= 5.806218e+002 n= 6 积分值I(n)= -5.800410e+005 */76 郑州大学2014-2015学年硕士研究生课程 数值分析 Numerical Analysis §1.5 数值算法的稳定性 计算结果可靠码？首先当 ，所以 因此 是单调递减收敛于0的数列，因此计算结果严重错 误。产生计算谬误的原因是： 如果 ，则计算 时产生999倍 的误差，计算 时产生 倍 的误差，误差的传播像洪水冲破闸门一 样泛滥成灾！可见以上的算法是病态和不稳定的。 */76 郑州大学2014-2015学年硕士研究生课程 数值分析 Numerical Analysis §1.5 数值算法的稳定性 我们可以改变以上格式，取 按照以上关系编程计算可得 n= 0 积分值I(n)= 1.000500e-003 n= 1 积分值I(n)= 5.001668e-004 上限1/n =1.000000 n= 2 积分值I(n)= 3.334167e-004 上限1/n= 0.500000 n= 3 积分值I(n)= 2.500500e-004 上限1/n =0.333333 n= 4 积分值I(n)= 2.000333e-004 上限1/n =0.250000 n= 5 积分值I(n)= 1.666905e-004 上限1/n= 0.200000 n= 6 积分值I(n)= 1.428750e-004 上限1/n =0.166667 算法是稳定的 舍入误差得到 良好的控制 */76 郑州大学2014-2015学年硕士研究生课程 数值分析 Numerical Analysis §1.6 计算工具 本书涉及的数值算法已由前人用各种编程语言实现，它们经历了数十年的考验和无数专业人士的调试、优化，是成熟、可靠的。 在今后的学习和工作中，我们完全可以利用这些已有的算法库减轻编程负担，提高工作效率。可以将主要精力放在对问题的描述和对计算结果的分析上面，而把繁琐的算法实现和优化细节交给现成的科学计算与工程软件/库来处理。 然而，在应用中要想准确的选择及灵活使用已实现的各种算法，必须对该算法的数学理论有透彻的了解，并根据自己的课题做取舍。这也是工科各专业研究生学习本课程的切入点和动因。 */76 郑州大学2014-2015学年硕士研究生课程 数值分析 Numerical Analysis 数学软件包：针对某一类数值计算问题的子程序集，大部分 为自由软件。 常用的数学软件包有： ( ,/petsc) LAPACK(Linear Algebra PACKage),它是用于求解线性代数方程组、线性最小二乘问题和矩阵特征值问题的Fortran77子程序包； ITPACK，它是用于迭代求解大型稀疏线性代数方程组的Fortran77子程序包； §1.6 计算工具 */76 郑州大学2014-2015学年硕士研究生课程 数值分析 Numerical Analysis MINPACK，它