第一章数值计算方法绪论

1、*)()(xxxxxExErbAx ni, 3 ,2 中国石油大学（华东）理学院制作:刘新海)(*yEr*2*2*2*1*1*1rrxfyxxfyx&教材教材 (Text Book) 数值计算方法数值计算方法（第二版）（第二版） 李维国、同登科李维国、同登科主编主编 （中国石油大学出版社，（中国石油大学出版社，2009年）年）& 参考书目参考书目 (Reference) 数值计算方法数值计算方法 (上、下册上、下册) 林成森林成森 编著编著 （科学出版社（科学出版社 1998年）年） Principle of Numerical Analysis 数值分析原理数值分析原理 封建

2、湖、车刚明、聂玉峰封建湖、车刚明、聂玉峰 编著编著 （科学出版社，（科学出版社，2001年）年） 提问：提问：数值计算方法是做什么用的？数值计算方法是做什么用的？研究对象研究对象：数值问题数值问题有限个输入数据（问题的自有限个输入数据（问题的自变量、原始数据）与有限个输出数据（待求解数据）之变量、原始数据）与有限个输出数据（待求解数据）之间函数关系的一个明确无歧义的描述。间函数关系的一个明确无歧义的描述。如一阶微分方程初值问题如一阶微分方程初值问题2(0)1dyxdxy求函数解析表达式求函数解析表达式( )yy x数学问题数学问题求函数求函数 在某些点在某些点( )yy x 1niix的近似函

3、数值的近似函数值数值问题数值问题程序程序设计设计上机上机计算计算设计高效、可设计高效、可靠的数值方法靠的数值方法数值数值问题问题求解求解近似结果近似结果输出输出重点讨论重点讨论数值问题的来源：数值问题的来源：实际实际问题问题建立数学模型建立数学模型数值数值问题问题数值方法的设计原数值方法的设计原则则收敛性：方法的可行性收敛性：方法的可行性稳定性：初始数据等产生的误差对结果的影响稳定性：初始数据等产生的误差对结果的影响便于编程实现：逻辑复杂度要小便于编程实现：逻辑复杂度要小计算量要小：时间复杂度要小，运行时间要短计算量要小：时间复杂度要小，运行时间要短存贮量要尽量小：空间复杂度要小存贮量要尽量小

4、：空间复杂度要小可可靠靠性性分分析析计算复杂性计算复杂性误差估计：运算结果不能产生太大的偏差且误差估计：运算结果不能产生太大的偏差且能够控制误差能够控制误差1 误误 差差 /* Error */一、一、 误差的来源与分类误差的来源与分类 /* Source & Classification */ 1、从实际问题中抽象出数学模型从实际问题中抽象出数学模型 模型误差模型误差 /* Modeling Error */ 2、通过观测得到模型中某些参数（或物理量）的值通过观测得到模型中某些参数（或物理量）的值 观测误差观测误差 /* Measurement Error */ 3、数学模型与数值算

5、法之间的误差求近似解数学模型与数值算法之间的误差求近似解 方法误差方法误差 (截断误差截断误差 /* Truncation Error */ ) 4、由于机器字长有限，原始数据和计算过程会产生新的误差由于机器字长有限，原始数据和计算过程会产生新的误差 舍入误差舍入误差 /* Roundoff Error */二、二、 误差分析的基本概念误差分析的基本概念 /* Basic Concepts */ 设设 为真值（精确值），为真值（精确值）， 为为 的一个近似值的一个近似值 称称 为近似值为近似值 的绝对误差，简称误差。的绝对误差，简称误差。 xxx x exx注：注：误差可正可负，常常是无限位的

6、误差可正可负，常常是无限位的绝对绝对误差限误差限/* accuracy */ 绝对值的上界绝对值的上界 exx 如：如：53 141591103 14159262.(.) 绝对绝对误差还不能完全表示近似值的好坏误差还不能完全表示近似值的好坏（绝对误差（绝对误差/* absolute error */）1 1 .Def近似值近似值 的误差的误差 与准确值与准确值 的比值：的比值：x e xexxxx 称为近似值称为近似值 的相对误差，记作的相对误差，记作reex x 注：注：实际计算时，相对误差通常取实际计算时，相对误差通常取rexxexx 因为因为221()()()()eeeexxexexxx

7、xxxex 1 2 .Def（相对误差（相对误差/* relative error */ ）相对相对误差也可正可负误差也可正可负rexxxx 1 3 .Def（有效数字有效数字/*Significant Digits*/ ）相对相对误差限误差限相对相对误差的绝对值的上界误差的绝对值的上界r /* relative accuracy */如：如：3 1415926. 3 14. 3 141592. 3位位21102e6位位51102e若近似值若近似值 与准确值的误差与准确值的误差绝对绝对值不超过某一位的值不超过某一位的半个单位，该位到半个单位，该位到 的第一位非零数字共有的第一位非零数字共有 位

8、，则位，则x x nx n称称 有有 位位有效数字有效数字 有效数字有效数字(另外一种定义形式另外一种定义形式)用科学计数法，记用科学计数法，记 其中其中 ， 若若 （即（即 的截取按四舍五入规则），则称的截取按四舍五入规则），则称 为有为有 位有效位有效数字，精确到数字，精确到 。12010mnx.a aa 01 a| 0 5 10m nxx.nax nm 10 120 19, ,nmZ a aan3.1415926535897932;3.1415 例例1：问：问： 有几位有效数字？请证明你的结论。有几位有效数字？请证明你的结论。 *10501050*and103141504131 .| ,

9、.*证明：证明：有有 位有效数字，精确到小数点后第位有效数字，精确到小数点后第 位。位。43规格化形式规格化形式若若 的每一位都是有效数字，则的每一位都是有效数字，则 称是称是有效数有效数x x 特别，经特别，经“四舍五入四舍五入”得到的数均为有效数得到的数均为有效数 一个算法如果输入数据有扰动（即误差），而计算一个算法如果输入数据有扰动（即误差），而计算过程中舍入误差不增长过程中舍入误差不增长, ,则称此算法是数值稳定的，否则则称此算法是数值稳定的，否则此算法就称为不稳定的。此算法就称为不稳定的。 1 4 .Def（数值稳定性（数值稳定性/ /* Numerical Stability */

10、） 对数学问题本身如果输入数据有微小扰动，引起对数学问题本身如果输入数据有微小扰动，引起输出数据（即问题真解）的很大扰动，这就是病态问输出数据（即问题真解）的很大扰动，这就是病态问题。题。1 5 .Def（病态问题（病态问题/ /* ill-posed problem */） 它是数学问题本身性质所决定的，与算法无关，它是数学问题本身性质所决定的，与算法无关，也就是说对病态问题，用任何算法（或方法）直接计也就是说对病态问题，用任何算法（或方法）直接计算都将产生不稳定性。算都将产生不稳定性。 三、数值算法及稳定性三、数值算法及稳定性 /* Numerical Algorithm and Stab

11、ility */此公式此公式精确精确成立成立80001050 .IIE记为记为*0I632120560111100.edxeeIx 则初始误差则初始误差111111110010 nI)e(ndxexeIdxexennnn10109111012111312141315141 10 36787944.1 100 088128001 110 030592001 120 632896001 137 22764801 1494 9594241 151423 3914II.II.II.II.II.II.II. What happened?!例例2 2 计算计算1010 1 2, , ,.nxnIx e d

12、xne 11101011nxnxnnIx enxe dxn Ie 公式一：公式一：考察第考察第n步的误差步的误差nE11|(1)(1)|nnnnnEIInInI| !01En|Enn 我们有责任改变。我们有责任改变。造成这种情况的是造成这种情况的是不稳定的算法不稳定的算法 /* unstable algorithm */迅速积累，误差快速递增。迅速积累，误差快速递增。初始的小扰动初始的小扰动801050| .E)1(1111nnnnInIInI 公式二：公式二：注意此公式与公式一注意此公式与公式一在理论上等价。在理论上等价。方法：先估计一个方法：先估计一个IN , ,再反推要求的再反推要求的I

13、n ( n N )。11)1(1 NINeN1112(1)1NNIIe NN 可取可取0* NNNIIEN, ,时时当当()()()()()()1514151314121311121011121110 04274623321616110 06381691815110 06687022014110 07177921413110 07735173212110 08387711511110 367879442I.eII.II.II.II.II.II. ()01110 632120561II.取取考察反推一步的误差：考察反推一步的误差：()()1111|11|NNNNEIIENNN 以此类推，对以此类

14、推，对 n N 有：有：|)1(.)1(1|NnEnNNE 误差逐步递减误差逐步递减, 这样的算法称为这样的算法称为稳定的算法稳定的算法 /* stable algorithm */ 在我们今后的讨论中，在我们今后的讨论中，误差误差将不可回避，将不可回避， 算法的算法的稳定性稳定性将会是一个非常重要的话题。将会是一个非常重要的话题。例例3：蝴蝶效应：蝴蝶效应 纽约的一只蝴蝶翅膀一拍，风和日纽约的一只蝴蝶翅膀一拍，风和日 丽的北京就刮起台风来了？！丽的北京就刮起台风来了？！纽约纽约北京北京这是一个这是一个病态问题病态问题蝴蝶效应蝴蝶效应2 误差分析的方法和原则误差分析的方法和原则 /* Erro

15、r Analysis */一、误差分析的方法一、误差分析的方法1、向前误差分析法向前误差分析法：利用误差限，：利用误差限，随着计算过程随着计算过程 逐步向前进行分析，直至估计出最后的结果。逐步向前进行分析，直至估计出最后的结果。1212()()()xxxx121221()()()xxxxxx12211222()()()xxxxxxx 注：注：两个近似数两个近似数 ，四则运算得到的误差限分别为，四则运算得到的误差限分别为12,xx（1 1）（2 2）对于函数对于函数 y = f (x)，若用，若用 x* 取代取代 x，将对，将对y 产生什产生什 么影响么影响？分析：分析：e*(y) = f (x

16、*) f (x) e*(x) = x* xMean Value Theorem( )()fxx x* 与与 x 非常接近时，可认为非常接近时，可认为 ，则有：，则有：( )()ffx ( )()( )eyfxex 即：即： 产生的误差经过产生的误差经过 作用后被作用后被放大放大/ /缩小缩小了了 倍。倍。故称故称 为为放大放大/缩小缩小因子因子 /* amplification factor */ 或或 绝对绝对条件数条件数 /* absolute condition number */.x f()f x ()f x ( )|( )|()reyeyf x ( )|( )|rexexx ()(

17、)( )()()( )()rrf xf xxxxeyxxf xxxfxexf x 相对误差条件数相对误差条件数 /* relative condition number*/ f 的条件数在某一点是的条件数在某一点是小小 大大，则称则称 f 在该点是在该点是好条好条件的件的 /* well-conditioned */ 坏条件的坏条件的 /* ill-conditioned */。注：注：关于多元函数关于多元函数 可类似讨论可类似讨论，理论工具：理论工具：Taylor公式公式 (教材第教材第7页页).,(21nx,x,xfy ()( )()fxexf x 例例4 4105%x 设设 ,试求函数试

18、求函数 的相对误差限的相对误差限. ( )nf xx 解：解：由题设知由题设知:近似值为近似值为 ,绝对误差限为绝对误差限为 10 x ()5%x 1111()()nnfxxxnnx ()() ()()0.005()()()re ffx e xe xeff xf xnxn 2、向后误差分析法向后误差分析法：把舍入误差的累积与导出：把舍入误差的累积与导出 的已的已知知 量量 的某种摄动（微小误差）等价起来，的某种摄动（微小误差）等价起来， A 12,nx xx即令即令1122(,)nnAf xxx 利用摄动理论，由利用摄动理论，由 的界估计出最后的舍入误差界。的界估计出最后的舍入误差界。i 3、

19、区间分析法区间分析法：把参加运算的数都看成区间量，：把参加运算的数都看成区间量，根根据区间运算规则求得最后结果的近似值和误差限。据区间运算规则求得最后结果的近似值和误差限。4、概率分析法概率分析法：利用概率统计方法，将：利用概率统计方法，将数据和运算中数据和运算中 的舍入误差视为适合某种分布的随机变量，然后确定计的舍入误差视为适合某种分布的随机变量，然后确定计算结果的误差分布。算结果的误差分布。 二、几点注意事项二、几点注意事项 /* Remarks */1、 避免相近二数相减避免相近二数相减例：例：a1 = 0.12345，a2 = 0.12346，各有，各有5位有效数字。位有效数字。 而而

20、 a2 a1 = 0.00001，只剩下，只剩下1位有效数字位有效数字。 几种经验性避免方法：几种经验性避免方法：;xxxx ;1lnlnln xxx当当 | x | 1 时：时：;2sin2cos12xx .6121112xxxex2、 避免小分母避免小分母 : 分母小会造成浮点溢出分母小会造成浮点溢出 /* over flow */3、避免大数避免大数吃吃小数小数例：例：用单精度计算用单精度计算 的根。的根。010)110(992 xx精确解为精确解为110291 x,x 算法算法1 1：利用求根公式利用求根公式aacbbx242 在计算机内，在计算机内，109存为存为0.1 1010，1存为存为0.1 101。做加法时，做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。即两加数的指数先向大指数对齐，再将浮点部分相加。即1 的指数部分须变为的指数部分须变为1010，则：，则：1 = 0.0000000001 1010，取取单精度时就成为：单精度时就成为： 109+1=0.10000000 1010+0.00000000 1010=