计算方法 第1章 预备知识与误差分析.pptxVIP

计算方法;1.1绪论; 由于本课程的授课对象主要是计算机学院和软件学院的本科学生，在算法的理论基础方面不追究精巧的数学推导，力图在算法的编程实现方面给出一些作者的体会，在设计算法的各个步骤时，尽量给出便于改写成程序语言的计算公式。 附录给出本课程主要算法的C程序，这些程序仅用于辅助算法过程的学习和数值试验，不一定适用用于商业软件开发，所以程序中没有过多的考虑数值计算中的异常情况，如果将这些程序用于商业软件开发，还需在程序中添加异常情况的保护机制。 ;1. 本课程主要内容;2. 学习计算方法的重要性和必要性;例1.1 考虑一个求解n阶线性方程组的问题：; 通常度量一个算法的时间复杂度是统计该算法中乘除法的运算次数，我们知道每计算一个n阶行列式，需要乘法次数为n!，克莱姆法则求解线性方程组需要计算n+1个行列式，所需总的乘法次数为： ;绪论; 我们再看一个实例，从中可以发现，有时直接使用高等数学中给出的很简单明了的数学表达式进行计算并不一定能够得到我们预期的结果。 例1.2 考虑导数的近似计算问题，根据导数的定义 ;表1.1 差商近似导数的单精度计算结果;表1.2 差商近似导数的双精度计算结果; 从以上两个表格可以看出，用差商近似导数并非Δx愈小愈精确，但是按高等数学中导数的定义Δx愈小愈接近导数值 ，实际计算结果似乎有悖于高等数学中所学的知识！ 但是如果用本课程中的误差分析的理论则很容易解释此类现象为何发生，因此具有这些误差分析的知识对编程人员来说将是很重要的也是必须的。 ;1.2 误差分析; (1)模型误差、根据实际问题建立数学模型的过程中通常会忽略某些次要因素而对问题进行简化，由此产生的误差称为模型误差。例如自由落体运动规律的公式 是忽略了空气阻力的影响，因此产生的误差就是模型误差。 ; (3)截断误差、通常在选定数学模型及参数后还要选择数值算法，很多数值算法往往是将一个无限过程截断为有限过程，此类误差称为截断误差或方法误差。例如前述的差商近似导数的计算方法实质上是在如下展开式中; (4)舍入误差、因为计算机表示浮点数的字长有限(通常4字节或8字节)，除极少量能够准确表示的数据外绝大部分是在超过界定的某位4舍5入所得的近似值，比如圆周率等存入计算机内存时均需要在某位进行4舍5入，此种误差称为舍入误差。 虽然舍入误差在单步运算时或许微不足道，但是在一个较为复杂且连贯的数值算法中，舍入误差可能会积累、传播，所以我们必须加以关注。; 2. 绝对误差 定义1.1 设x是准确值，x*是它的一个近似值，称 为x的绝对误差，简称误差，误差的绝对值的上限ε* 称为x*的绝对误差限，简称误差限。 例如在使用圆周率时通常取为π*=3.14，其绝对误差及误差限为 在工程技术或商品规格中通常用如下方式表示产品的误差限：; 3. 相对误差 度量误差的另一种方式是相对误差，例如从10±0.1和1000±1中我们肯定不会认为绝对误差限小的比绝对误差限大的更精确。 定义1.2 设x为准确值, x*是它的一个近似值，称比值 为近似值x*的相对误差，相对误差的绝对值的上限记为 称为x*的相对误差限，即有 ; 对于x=10±0.1和1000±1，我们可以计算出它们的相对误差分别为： 从绝对误差看，后者误差限更大些，但从相对误差的角度看，后者更为精确。 在(1.9)式中，在可以确定准确值x的前提下也可以取为x作为分母。 ; 4. 有效数字 定义1.3 设x为准确值，其近似值为 这里 为整数，且 则称近似值x*有p位有效数字。 例1.3 圆周率 ，如果取3.1415和 3.14159作为其近似值，问分别有几位有效数字？ 解： (1) 取 则3.1415有4为有效数字，最后一位数字5是无效数字。 ; (2) 取 则取3.14159有6为有效数字。 我们在近似计算中四舍五入保留小数点以后n位数字时有效数字的个数，如同样保留小数点后4位数字 1.414213… 近似值1.4142有5位有效数字 0.012353… 近似值0.0124有3位有效数字 28.08918… 近似值28.0892有6位有效数字 总之：对精确值进行四舍五入所得结果后所有非零数字均是有效数字。 ; 定理1.1 设近似值x*具有(1.11)式的形式，如果有n位有效数字或是由其准确值4舍5入后得到的近似值，则x*的相对误差限为： 证明：如果 有n位有效数字，则其绝对误差限为 相对误差限为 如果x*是由其准确值4舍5入后得到的近似值，则此近似值必然有n位有效数字，从而结论也必然成立。 ; 5. 近似计算中减少误差的几个策略 1) 避免两个相近的数相减 我