计算方法 第一章 数值计算方法.pptVIP

本课程的基本要求 掌握数值方法的基本原理 掌握常用的科学与工程计算的基本方法 能用所学方法在计算机上算出正确结果  数学问题的数值解法例示 解 数学问题的数值解法例示 现取h=0.05,其结果见下表: 本章内容 §1 引言 §2 误差的来源及分类 §3 误差的度量 §4 误差的传播 §5 减少运算误差的原则 要求掌握的内容 1.1 引言 数值分析又称计算方法, 它是研究各种数 学问题的数值解法及其理论的一门学科。 1. 对于要解决的问题建立数学模型 2. 研究用于求解该数学问题近似解的算法和过程 3. 按照2进行计算，得到计算结果 在任何科学计算中其解的精确性总是相对的,而误差则是绝对的.我们从下面这个例子就可以了解误差产生的原因. 试求摆长为L的单摆运动周期. 在用数值方法解题过程中可能产生的误差归纳起来有如下几类： 1. 模型误差 2. 观测误差 3. 截断误差 4. 舍入误差 用数学方法解决一个具体的实际问题，首先要建立数学模型，这就要对实际问题进行抽象、简化，因而数学模型本身总含有误差，这种误差叫做模型误差 数学模型是指那些利用数学语言模拟现实而建立起来的有关量的描述 数学模型的准确解与实际问题的真解不同 在数学模型中通常包含各种各样的参变量，如温度、长度、电压等，这些参数往往是通过观测得到的，因此也带来了误差，这种误差叫观测误差 数学模型中的参数和原始数据，是由观测和试验得到的 由于测量工具的精度、观测方法或客观条件的限制,使数据含有测量误差,这类误差叫做观测误差或数据误差 根据实际情况可以得到误差上下界 数值方法中需要了解观测误差,以便选择合理的数值方法与之适应 精确公式用近似公式代替时,所产生的误差叫截断误差 例如, 函数f(x)用泰勒(Taylor)多项式 在数值计算中只能对有限位字长的数值进行运算 需要对参数、中间结果、最终结果作有限位字长的处理工作，这种处理工作称作舍入处理 用有限位数字代替精确数，这种误差叫做舍入误差，是数值计算中必须考虑的一类误差 多元函数误差估计 例题 （2）防止大数“吃掉”小数 例 求二次方程x2-105x+1=0的根 解：按二次方程求根公式 x1=(105+(1010-4)1/2)/2 x2=(105-(1010-4)1/2)/2 在8位浮点数计算得 x1=(105+105 )/2=105 (正确）, x2=(105-105 )/2=0 (错误） （3）绝对值太小的数不宜做除数 当分母为两个相近数相减时,会丧失有效数字 例1.8 计算 算法2。分成三组因子。每组只取六位小数计算 a=0.0005/0.0003=1.666667(有舍入) b=0.0143/0.0125=1.144000 c=0.0012/0.0135=0.088889 (有舍入) D=a*b*c=1. 666667* 1.144000* 0.088889 =0.169482,准确到小数后5位。 例题 1.3.4 有效数字与相对误差 定理1.1 若近似数x*=?0.x1x2…xn?10m具有 n 位 有效数字，则其相对误差 证: ∵ x* = ?0.x1x2…xn?10m ∴ ?x*? ≥x1?10 m-1 又 ∵ x*具有n位有效数字,则?x- x*?≤1/2?10 m - n ∴ 一般应用中可以取?r*=1/2x1 ?10-(n-1),n越大,?r*越小, ∴有效数字越多，相对误差就越小 1.3.4 有效数字与相对误差 例7 取3.14作为?的四舍五入的近似值时，求其 相对误差 解：3.14=0.314 ?101 x1=3 m=1 ∵ 四舍五入的近似值,其各位都是有效数字 ∴ n=3 ?r*=1/2x1 ?10-(n-1)=1/2*3 ?10-2=17% 例8 已知近似数x*有两位有效数字，试求其相 对误差限 1.3.4 有效数字与相对误差 解：已知 n=2 代入公式 ?r*=1/2x1 ?10-(n-1)得 ?r*=1/2x1 ?10-1 x*的第一位有效数字x1没有给出，可进行如下讨论：当 x1=1 ?r*=1/2x1 ?10-1=1/2*1 ?10-1=5% x1=9 ?r*=1/2x1 ?10-1=1/2*9 ?10-1=0.56% 取 x1=1 时相对误差为最大，即 5%