数值分析与方法01.pptVIP

第一章 数值计算方法与误差分析 第一节 数值计算方法 第二节 误差的来源 模型误差 观测误差 截断(方法)误差 舍入误差 例1:函数f(x)用泰勒多项式 Pn(x)=f(0)+ 近似代替, 则截断误差Rn(x)=f(x)- Pn(x) = ( 介于0与x之间) 例2: 用3.14159近似 ,误差R= -3.14159=0.0000026… 为舍入误差. 第三节 近似数的误差表示法 1、绝对误差、绝对误差限 设x为准确值，x*为x的一个近似值，称e*=x*-x为近似值的(绝对)误差。 若|x-x*|≤ε*，则称ε*为绝对误差限。 2、相对误差、相对误差限 称为近似值x*的相对误差，记作er*。 通常取er*= ，相对误差限εr*= 3、有效数字 若近似值x*的误差限是某一位的半个单位，该位到x*的第一位非零数字共有n位，则称x*有n位有效数字。 x*有n位有效数字的标准形式为: x*=±10m×(a1+a2×10-1+…+an×10-(n-1)) (2.1) 其中,a1是1到9中的某个数字,a2,a3,…,an为0到9中的数字,m为整数,且|x-x*|≤ (2.2) 定理1.1 若用(2.1)表示的近似数x*具有n位有效数字,则其相对误差限为 反之,若x*的相对误差限 则x*至少具有n位有效数字。 证明：由(2.1)可得a1×10m≤|x*|≤(a1+1)×10m 当x*有n位有效数字时, 反之,由 所以x*有n位有效数字。 第四节 数值运算误差分析 设f(x)为一元函数，x的近似值为x*，以f(x*)近似f(x)，其误差界 ，可用泰勒展式 其中， 介于x，x*之间， 则 第五节 减小运算误差若干原则 1、避免两个相近数相减 因为两个相近数相减会使有效数字严重损失。解决办法是改变计算公式，或增加有效位数。 一般地，当f(x)≈f(x*)时，可用泰勒展式近似： 取右端的有限项近似左端。 2、防止大数“吃掉”小数 参加运算的数的数量级相差很大时，因为计算机的位数有限，需注意运算次序。 3、简化计算步骤，减少运算次数 可节省计算时间，减少舍入误差。 4、绝对值太小的数不宜做除数 设x，y分别有近似值x*，y*，z= 的近似值 则其绝对误差 |e*(z)|≈ 当|y*|很小时，近似值z*的绝对误差e*(z)可能很大。 5、注意计算过程中误差的传播积累, 防止误差恶性放大 * *