西安石油大学现代数值计算方法第1章

1、课程名称：f数值分析f或计算方法(考试时要带上计算器(有函数的)第一章第一章 数值计算中的误差分析数值计算中的误差分析（绪论）（绪论）内容提要内容提要 1.1 1.1 计算方法的任务与特点计算方法的任务与特点 1.2 1.2 误差与误差估计误差与误差估计( (误差知识误差知识) ) 1.3 1.3 选用算法时应遵循的原则选用算法时应遵循的原则1.11.1计算方法的任务与特点计算方法的任务与特点f科学与工程计算过程：科学与工程计算过程：f实际问题实际问题数学模型数学模型数值问题数值问题算法算法程程序序调试调试结果结果计算机用途分类：计算机用途分类：科学计算、数据处理科学计算、数据处理计算方法的特

2、点f严密的科学性、实验的技术性、高度的严密的科学性、实验的技术性、高度的抽象性、应用的广泛性、提供算法、算抽象性、应用的广泛性、提供算法、算法分析、兼顾计算机的特点：有效数字法分析、兼顾计算机的特点：有效数字（精度）、运算量、存储量等。（精度）、运算量、存储量等。第一章第一章 绪论绪论第二章第二章 线性方程组求解线性方程组求解第三章第三章 非线性方程求解非线性方程求解第四章第四章 矩阵特征值问题（不讲）矩阵特征值问题（不讲）第五章第五章 函数的插值函数的插值第六章第六章 曲线拟合曲线拟合第七章第七章 数值积分和数值微分数值积分和数值微分第八章第八章 常微分方程数值解法常微分方程数值解法 本课程

3、主要内容本课程主要内容1.2 1.2 误差知识误差知识（误差与数值计算中的误差估计）（误差与数值计算中的误差估计） 内容提要：内容提要：1.1.误差的来源及其分类误差的来源及其分类2.2.误差的度量（误差与有效数字）误差的度量（误差与有效数字）3.3.数值计算的误差估计数值计算的误差估计一、误差来源及其分类一、误差来源及其分类1)1)模型误差（描述误差）模型误差（描述误差） 反映实际问题有关量之间的计算公式（反映实际问题有关量之间的计算公式（数学模型）通常是近似的。数学模型）通常是近似的。2 2）观测误差）观测误差3 3）截断误差（方法误差）截断误差（方法误差） 数值方法精确解与待求解模型的理

4、论分数值方法精确解与待求解模型的理论分析解之间的差异。析解之间的差异。 这是由于我们需要将无穷过程截断为有这是由于我们需要将无穷过程截断为有限过程，而使得算法必须在有限步内执行结限过程，而使得算法必须在有限步内执行结束而导致的。束而导致的。neeee,!n1!21! 111,!21! 111n例如：例如：4 4）舍入误差）舍入误差 以四舍五入为例（也可以五舍六入等）最多舍去或添加最后一位的半个单位。注意：与截断误差不同！二、误差的度量二、误差的度量1)绝对误差2)相对误差3)有效数字4)度量间的关系1. 绝对误差f绝对误差定义：准确值绝对误差定义：准确值 x 减近似值减近似值 x*)(*xex

5、x 绝对误差限：绝对误差限：书上有错！改正2.相对误差 绝对误差限虽然能够刻划对同一真值不同近似绝对误差限虽然能够刻划对同一真值不同近似的好坏，但它不能刻划对不同真值近似程度的好的好坏，但它不能刻划对不同真值近似程度的好坏坏 。有效数：当有效数：当x*为为四舍五入得到的近似数，四舍五入得到的近似数，则称则称x*为有效数。为有效数。有效数的绝对误差限、相对误差限，有效数字位数举例：例若有效数x*=1.02,则绝对误差限： (x*)=0.005相对误差限： (x*)= (x*)/| x*|= 0.0049，x*具有3位有效数字.若有效数x*= 2500，则绝对误差限 (x*)=0.5，相对误差限

6、(x*)= (x*)/| x*|= 0.0002x*具有4位有效数字.f若有效数x*= 25102，则 (x*)=50 ，(x*)= (x*)/| x*|= 0.02，x*具有2位有效数字.注意：有效数 2500=25.00102但与 25102的误差不同；当然与一个准确数2500也不同。关于某一位的半个单位关于某一位的半个单位3.3.有效数字（教材第有效数字（教材第6 6页中间）：页中间）：f若近似数若近似数x x* *的绝对误差限是（不超过）的绝对误差限是（不超过）某一位的半个单位，则称其精确到这一某一位的半个单位，则称其精确到这一位，且从该位到位，且从该位到x x* *的第一位非零数字共

7、的第一位非零数字共有有n n位，则称近似数位，则称近似数x x* *具有具有n n位有效数字位有效数字.f（黑板举例）分两种情况：有效数的有效数字位数。若既知道x*（非有效数），又知道x时，如何求有效数字位数？见下例：f举例：举例：x=3.1415926x=3.1415926, ,f近似数近似数 x x1 1* *=3.14102, x=3.14102, x2 2* *=3.142=3.1422*11021005. 000057. 0 xx3*210210005. 000040. 0 xx3 3位有效数字，非有效数！位有效数字，非有效数！4 4位有效数字，有效数位有效数字，有效数注意：一个有效

8、数，若知道对应的准确数，此时所求误差限可能不同。而有效数字位数怎样求都一样。 数的标准形式fremark1remark1: : 有效数的误差限是末位数单位的一半，有效数的误差限是末位数单位的一半，可见有效数本身就体现了误差界。可见有效数本身就体现了误差界。fremark2remark2: : 对真值进行四舍五入得到有效数。对真值进行四舍五入得到有效数。fremark3remark3: :准确数字有无穷多位有效数字。准确数字有无穷多位有效数字。fremark4remark4: : 从实验仪器所读的近似数（最后一为是从实验仪器所读的近似数（最后一为是估计位）不是有效数，估计最后一位是为了确保对估计

9、位）不是有效数，估计最后一位是为了确保对最后一位进行四舍五入得到有效数。最后一位进行四舍五入得到有效数。 例例 从最小刻度为厘米的标尺读得的数据从最小刻度为厘米的标尺读得的数据123.4cm123.4cm是为了得到有效数是为了得到有效数123.cm,123.cm,读得数据读得数据156.7cm156.7cm是为了得到有效数是为了得到有效数157.cm157.cm。4.4.误差度量间的联系误差度量间的联系f绝对误差与相对误差绝对误差与相对误差nmxe10)(21*f绝对误差与有效数字（教材第绝对误差与有效数字（教材第7页页1.2.2式）式）f相对误差与有效数字（相对误差与有效数字（教材教材第第8

10、页页1.2.3式）式）)(/ )(*xexxer定理证明o1o2证毕,10) 1(1011*11mmxnmxe1021)(*remarkf1 1、该定理实质上给出了一种求相对误差限的方法。、该定理实质上给出了一种求相对误差限的方法。f2 2、仅从、仅从 并不能保证并不能保证x x* *一定具有一定具有n n位有效数字。如位有效数字。如 设其近似值设其近似值a=0.484a=0.484，其相对误差为：，其相对误差为：我们并不能由此断定我们并不能由此断定a a有两位有效数字，因为有两位有效数字，因为nrxe111021*）（4900. 00229sina21104210125.0012397.04

11、84.0484.04900.0201021005.00600.0484.04900.0 aa例题例题解：解：令01.01041n1， 解得2.4n . 取 n=3 位有效数字。 # 三、数值运算的误差估计三、数值运算的误差估计f由于自变量的误差，引起函数值的误差。也称为误差传播。这种函数值的误差往往很难精确描述，只能大体估计。f例如一元函数，二元函数，等。y*的绝对误差e(y*)= y-y*,相对误差er(y*)=( y-y*)/y*,复习泰勒公式复习泰勒公式)()()()(! 11)()(*111*nnnxxxpfxxxpfpfpf2*2*11*1*2*222*2*11*2212*2*111

12、*22*1111*2)()()()( )()()()( )()()()(! 21nnnnnnnnnnnnnxxxxpfxxxxxxpfxxxxxxpfxxxxxxpfxxxxxxpfxxxxpf泰勒公式分析初值误差传播泰勒公式分析初值误差传播相对误差（教材第相对误差（教材第9页页1.2.6式）：式）：进而得到如下绝对误差限和相对误差限传播关系：进而得到如下绝对误差限和相对误差限传播关系： 例题：参见教材第例题：参见教材第10页例页例3求导数复习：基本初等函数的导数；四则运算的导数；复合函数求导。v选用数值稳定性好的算法。v定义定义：一个算法：一个算法, , 如果在运算过程中舍入误差如果在运算过

13、程中舍入误差在一定条件下能够得到控制在一定条件下能够得到控制, , 或者舍入误差的或者舍入误差的增长不影响产生可靠的结果增长不影响产生可靠的结果, , 则称该算法是数则称该算法是数值稳定的值稳定的, , 否则称其为数值不稳定否则称其为数值不稳定. .dxxxinn105v例：计算如下积分近似值的两种方案比较),(21511nininn方法方法1：1823. 056ln51100dxxi1.3 选用算法应遵循的原则方法方法1 1计算结果计算结果方法一结果分析方法一结果分析n方法一分析：计算结果表明方法一分析：计算结果表明, , 舍入误差的传播近舍入误差的传播近似依似依5 5的幂次进行增长的幂次进

14、行增长, , 因而是一种不稳定的方因而是一种不稳定的方法。法。 方法二：方法二：由此分析知，该方法是稳定的。关于初值的由此分析知，该方法是稳定的。关于初值的近似可由下面式子得到：近似可由下面式子得到：) 1(5156) 1(611010ndxxidxxnnnn151nnini*1*51nnini0221555eeeennnn5511nnini551*1nnininnee511方法方法2 2计算结果计算结果016706615512110.*iv简化计算步骤以减少运算次数。简化计算步骤以减少运算次数。842842284488163*3*3*3*33*3*3*33*3*33*33例例1 1例2 秦九

15、韶算法 01223344axaxaxaxa01234)( ( (axaxaxaxa) 10(98765432110000001kkk例v尽量避免相近的数相减尽量避免相近的数相减 例例 x=52.127 xx=52.127 x* *=52.129 =52.129 四位有效数字四位有效数字 y=52.123 yy=52.123 y* *=52.121 =52.121 四位有效数字四位有效数字 a=x-y=0.004 aa=x-y=0.004 a* *=x=x* *-y-y* *=0.008 =0.008 零位有效数字零位有效数字结论：避免相近数相减v合理安排量级相差很大的数之间的运算合理安排量级相

16、差很大的数之间的运算次序次序, , 尽可能避免大数尽可能避免大数“吃掉吃掉”小数。小数。)(11111xxxxxxxx1lnln) 1ln(一些避免相近数相减示例一些避免相近数相减示例当当|x|1|x|1时时xxxx111)1ln()1ln(22xxxx当当|x|1|x|1时时2221111xxx53arctan53xxxx!5!3sin53xxxxv尽可能避免绝对值很小的数做分母，防尽可能避免绝对值很小的数做分母，防止出现溢出。止出现溢出。当当a,ba,b中有近似值时，由中有近似值时，由)0()()()(2bbaebbeabae若 ，则 可能很大。当a,b都是准确值时，由于 很大，会使其它较小的数加不到 中而引起严重误差，或者会发生计算机“溢出”，导致计算无法进行下去。ab