数值分析11讲义课件

数值分析NumericalAnalysis

数值分析

《数值分析》(第2版)

朱晓临主编中国科学技术大学出版社教材《数值分析》(第2版)教材《数值分析》（第5版）

李庆阳，王能超，易大义编著

清华大学大学出版社《数值分析》（第3版）

颜庆津著，

北京航空航天大学

出版社

《NumericalAnalysis》(Ninthed.)RichardBurden,DouglaFaires,

Brooks/Cole,CengageLearning,2011参考书目《数值分析》（第5版）参考书目

微积分

线性代数

常微分方程

算法语言预备知识微积分预备知识一、为什么要学习数值分析？课程简介现实世界的问题可以归结为各种各样的数学问题

—方程求根问题

—解线性方程组的问题

—定积分问题

—常微分方程初值问题

—…….一、为什么要学习数值分析？课方程求根问题在科学计算中常要遇到求解各种方程，例如：高次代数方程x5－3x＋7＝0超越方程方程求根问题在科学计算中常要遇到求解各种方程，例如：高次代数

高次线性方程和超越方程看似简单，但难于求其精确解。对于高次代数方程，由代数基本定理知多项式根的数目和方程的阶相同，但对超越方程就复杂的多，如果有解，其解可能是一个或几个，也可能是无穷多个。高次线性方程和超越方程看似简单，但难于求其精确解。对解线性方程组的问题线性方程组的一般形式解线性方程组的问题线性方程组的一般形式(1)当b≠0时称为非齐次线性方程组，其可能有唯一解、无解或者无穷多个解。当b=0时称为线性齐次方程组，必有零解。(2)由线性代数知识可知：当系数矩阵A非奇异(即detA≠0)时，方程组有唯一解，可用克莱姆法则求解，但它只适合于n很小的情况。(1)当b≠0时称为非齐次线性方程组，其可能有唯一解、无解或克莱姆法则其中克莱姆法则其中例如用克莱姆法则求解一个n阶方程组，要算n+1个n阶行列式的值，总共需要

n!(n-1)(n+1)次乘法。当n充分大时，计算量是相当惊人的。例如用克莱姆法则求解一个n阶方程组，要算n+1个n阶行列比如一个20阶不算太大的方程组，大约要做1021次乘法，这项计算即使每秒1亿次浮点数乘法计算的计算机去做，也要连续工作30万年才能完成。当然这是完全没有实际意义的，故需要寻找有效算法!比如一个20阶不算太大的方程组，大约要做1021次乘法，这项定积分问题对于积分由微积分知识可知：只要找到被积函数f(x)的原函数F(x)，便有下列牛顿—莱布尼兹公式定积分问题对于积分由微积分知识可知：只要找到被积函数f(x)为何要进行数值积分？原因之一：许多形式上很简单的函数，例如已证明它们的原函数不能用初等函数表示成有限形式。为何要进行数值积分？原因之一：许多形式上很简单的函数，例如已原因之二：有些被积函数的原函数过于复杂，例如的一个原函数是原因之二：有些被积函数的原函数过于复杂，例如的一个原函数是要计算f(x)定积分的近似值，上式就不见得方便。要计算f(x)定积分的近似值，上式就不见得方便。原因之三：f(x)以离散数据点形式给出xix0x1…xnyi=f(xi)y0y1…yn原因之三：f(x)以离散数据点形式给出xix0x1…xny常微分方程初值问题一阶常微分方程的初值问题，即例常微分方程初值问题一阶常微分方程的初值问题，即例常微分方程的一般解(解析解)

对一些典型的微分方程(可分离变量方程，一阶线性方程等等)，有可能找出它们的一般解表达式，然后用初始条件确定表达式中的任意常数，这样解即能确定。常微分方程的一般解(解析解)例如

求解解：分离变量得dy=2xdx

积分得y=x2+c

由初值得c=0

故解为y=x2

例如求解解：分离变量得dy=2xdx但是对于求解无法求出一般解！但是对于求解无法求出一般解！二、如何学习数值分析？课程简介1.注意掌握各种方法的基本原理2.注意各种方法的构造手法3.重视各种方法的误差分析4.做一定量的习题5.注意与实际问题相联系二、如何学习数值分析？课1.注意掌握各种方法的基本原理第一章

绪论数值分析第一章

绪论§1数值分析研究的对象与特点一、数值分析研究的对象二、数值分析研究的特点一、数值分析研究的对象二、数值分析研究的特点1.数值分析研究的对象数值分析是计算数学的一个主要部分，计算数学是数学科学的一个分支，它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现.1.数值分析研究的对象数值分析是计算数学的一个主要部分，计算用计算机解决科学计算问题通常经历以下过程实际问题数学模型数值计算方法程序设计上机计算结果应用数学计算数学用计算机解决科学计算问题通常经历以下过程实际问题数学模型数值2.数值分析研究的内容—函数的数值逼近(插值与拟合)—数值积分与数值微分—非线性方程数值解—数值线性代数—常微和偏微数值解，……2.数值分析研究的内容—函数的数值逼近(插值与拟合)数值分析实质上是以数学问题为研究对象，不像纯数学那样只研究数学本身的理论，而是把理论与计算紧密结合，着重研究数学问题的数值方法及理论。数值分析实质上是以数学问题为研究对象，不像纯数学那样只研究数3.数值分析的特点面向计算机，要根据计算机特点设计切实可行的有效算法.

有可靠的理论分析，能任意逼近并达到精度要求，对近似计算要保证收敛性和数值稳定性.3.数值分析的特点面向计算机，要根据计算机特点设计切实可行(3)要有好的计算复杂性，时间复杂性好是指节省时间，空间复杂性好是指节省存贮量，这也是建立算法要研究的问题.(4)要有数值试验，即任何一个算法除了从理论上要满足上述三点外，还要通过数值试验证明是行之有效的.(3)要有好的计算复杂性，时间复杂性好是指节省时间，空间复§2数值计算的误差一、误差来源的分类二、误差分析的重要性三、绝对误差和绝对误差限四、相对误差和相对误差限五、有效数字§2一、误差来源的分类二、误差分析的重要性三、绝对误差和绝对一、误差来源的类型1.模型误差

从实际问题中抽象出数学模型

——模型误差

/*ModelingError*/一、误差来源的类型1.模型误差从实际问题中抽象出数学模2.观测误差

通过测量得到模型中参数的值

——观测误差

/*MeasurementError*/注：通常根据测量工具的精度，可以知道这类误差的上限值。2.观测误差通过测量得到模型中参数的值注：通常根

当数学模型得不到精确解时，要用数值计算方法求它的近似解，由此产生的误差称为截断误差或方法误差3.截断误差

求近似解——方法误差(截断误差)

/*TruncationError*/

当数学模型得不到精确解时，要用数值计算方法求它的近似例如：在微积分中sinx可展开成

但在计算机中计算时，常用前几项来代替，即抛弃了无穷级数的后段，这样就产生了截断误差。

当|x|很小时，常用x代替sinx，其截断误差大约为x3/6。例如：在微积分中sinx可展开成但在计算机中计算时，

由于计算机字长有限，原始数据的输入及浮点数运算过程中都有可能产生误差，这样产生的误差称为舍入误差4.舍入误差

机器字长有限——舍入误差

/*RoundoffError*/由于计算机字长有限，原始数据的输入及浮点数运算过程中在数值分析课程中，主要研究

截断误差

舍入误差

在数值分析课程中，主要研究截断误差二、误差分析的重要性考察如下两个方程组试思考这两个方程组的解的关系？二、误差分析的重要性考察如下两个方程组试思考这两个方程组容易看出系数矩阵完全相同，而常数项矩阵有微小差别，右端系数1.9999变成2.0001，其误差为2.0001-1.9999=0.0002=0.02%容易看出系数矩阵完全相同，而常数项矩阵有微小差别，右端系数1但对应的解为由此看出系数矩阵完全相同，而常数项矩阵有微小差别的方程组，其解竟然相差得很大！解的最大误差=2=200%但对应的解为由此看出系数矩阵完全相同，而常数项矩阵有微小差别据说，美军1910年的一次部队的命令传递是这样的:营长对值班军官:明晚大约8点钟左右，哈雷彗星将可能在这个地区看到，这种彗星每隔76年才能看见一次。命令所有士兵着野战服在操场上集合，我将向他们解释这一罕见的现象。如果下雨的话，就在礼堂集合，我为他们放一部有关彗星的影片。值班军官对连长:根据营长的命令，明晚8点哈雷彗星将在操场上空出现。如果下雨的话，就让士兵穿着野战服列队前往礼堂，这一罕见的现象将在那里出现。连长对排长:根据营长的命令，明晚8点，非凡的哈雷彗星将身穿野战服在礼堂中出现。如果操场上下雨，营长将下达另一个命令，这种命令每隔76年才会出现一次。排长对班长:明晚8点，营长将带着哈雷彗星在礼堂中出现，这是每隔76年才有的事。如果下雨的话，营长将命令彗星穿上野战服到操场上去。班长对士兵:在明晚8点下雨的时候，著名的76岁哈雷将军将在营长的陪同下身着野战服，开着他那“彗星”牌汽车，经过操场前往礼堂。据说，美军1910年的一次部队的命令传递是这样的:三、绝对误差和绝对误差限

定义

设某一量的准确值为x，近似值为x*，则x*与x之差叫做近似值x*的绝对误差(简称误差)，记为？判断题：绝对误差是误差的绝对值．三、绝对误差和绝对误差限定义设某一量的准确值为x，绝对误差的性质(1)绝对误差e(x*)可正可负(2)|e(x*)|的大小标志着x*的精确度(3)绝对误差e(x*)未知绝对误差的性质(1)绝对误差e(x*)可正可负定义

若指定一个适当小的正数

，使

有时用

表示近似值x*的精度或准确值的所在范围。则

称为近似值

x*的绝对误差限。定义有时用表示绝对误差限的性质(1)在实际问题中，绝对误差一般是有量纲的，绝对误差限也是有量纲的。例如，测得某物体的长度为5m，其误差限为0.01m绝对误差限的性质(1)在实际问题中，绝对误差一般是有量纲的，(2)绝对误差限是正的，有无穷多个。若已知

是绝对误差限，由于则比大的任意正数均是绝对误差限。(2)绝对误差限是正的，有无穷多个。若已知是绝对误差

思考题：设有两个温度计，其一测量1000℃时的绝对误差限为5℃，而另一个测量100℃时的绝对误差限为1℃。问：哪一个温度计更精确？四、相对误差和相对误差限思考题：设有两个温度计，其一测量1000℃时的绝对误答：虽然后者绝对误差限的数值较小，但第一种温度计更为精确。

决定一个量的近似值的精确度除了要看绝对误差的大小外，还要考虑到该量本身的大小。答：虽然后者绝对误差限的数值较小，但第一种温度计更为精确。定义:

绝对误差与准确值之比称为x*的相对误差.定义:绝对误差与准确值之比称为x*的相对误差.(2)由于准确值x未知，故实际问题中，当||较小时，常取注(1)相对误差是个无量纲量,值小者精度高。(2)由于准确值x未知，故实际问题中，当|当||较小时，可用下式计算定义若指定一个适当小的正数，使则称

为近似值x*的相对误差限当||较小时，可用下式计算定义若指定一个适

当x有很多位数字时，常按照“四舍五入”原则取前几位数字作为x的近似值例：设x=π=3.1415926···取x1*=3作为π的近似值，则五、有效数字当x有很多位数字时，常按照“四舍五入”原则取前几位数取x2*=3.14作为π的近似值，则取x3*=3.1416作为π的近似值，则它们的误差都不超过某一数位的半个单位。取x2*=3.14作为π的近似值，则取x3*=若近似值x*的绝对误差限是某一位的半个单位，该位到x*的第一位非零数字一共有n位，则称近似值x*有n位有效数字，或说x*精确到该位。定义：有效数字若近似值x*的绝对误差限是某一位的半个单位，该位到x*的第一若将准确值x

的近似值x*表示成标准形式而其误差限则说近似值x*具有n

位有效数字。这里n

为正整数，m

为整数，每个均为0,1，…，9中的一个数字，注：有效数字的等价定义若将准确值x的近似值x*表示成标准形式而其误差限则说近似值定理1.1

若

x*具有n位有效数字，则其相对误差限满足反之，若x*的相对误差限满足则近似值x*至少具有n位有效数字。定理1.1若x*具有n位有效数字，则其相对反之，若x*例：为使的相对误差小于0.001%,至少应取几位有效数字？解：假设*取到n位有效数字，则其相对误差上限为要保证其相对误差小于0.001%，只要保证其上限满足已知a1=3，则从以上不等式可解得n>6log6，即n6，应取*=3.14159。例：为使的相对误差小于0.001%,至少应取几位有效数§3避免误差危害的若干原则(一)要避免相近两数相减(二)要防止大数“吃掉”小数(三)注意简化计算步骤，减少运算次数(四)要避免绝对值小的数

作除数(五)设法控制误差的传播§3(一)要避免相近两数相减(二)要防止大数“吃掉”小数(一)要避免相近两数相减的值。例:求当x=1000，y的准确值为0.01580.

(1)直接相减(2)若将原式改写为则y=0.01581(一)要避免相近两数相减的值。例:求当x=1000，几种经验性避免方法：很小，几种经验性避免方法：很小，(二)要防止大数“吃掉”小数例：用单精度计算的根。精确解为算法1：利用求根公式在计算机内，109存为0.11010，1存为0.1101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。即1的指数部分须变为1010，则：1=0.00000000011010，取单精度时就成为：109+1=0.100000001010+0.000000001010=0.100000001010(二)要防止大数“吃掉”小数例：用单精度计算算法2：先解出

再利用注：求和时从小到大相加，可使和的误差减小。算法2：先解出注：求和时从小到大相加，可使和的误差减小。例：在5位浮点十进制计算机上，计算y=54321+0.4+0.3+0.4解：若按从左到右的顺序进行计算，后三位在对阶过程变为例：在5位浮点十进制计算机上，计算y=54321

后三个数都在对阶过程中变为零，得出含有较大误差的结果y=54321。但若按从右到左的顺序进行计算，后三位在对阶过程变为

这种算法避免了大数“吃掉”小数！后三个数都在对阶过程中变为零，得出含有较大误差的结果一般地，有如下原则若干数相加，采用绝对值较小者先加的算法，结果的相对误差限较小一般地，有如下原则(三)注意简化计算步骤，减少运算次数，避免误差积累例：计算多项式的值解：如果先计算各项然后相加，则乘法次数=4+3+2+1=10，加法次数=4(三)注意简化计算步骤，减少运算次数，避免误差积累例：计算但如改用下式计算

则只需做4次乘法和4次加法。计算量大大减少！注：第二种方法称为“秦九韶算法”（

Horner算法）但如改用下式计算则只需做4次乘法和4次加法。计算量大通常，计算如下n次多项式的值如果先计算各项然后相加，则乘法次数=n+(n-1)+···+2+1=n(n+1)/2

加法次数=n通常，计算如下n次多项式的值如果先计算各项然后相加，则若采用“秦九韶算法”，则乘法次数=n

加法次数=n两种算法的乘法运算次数随n的变化见下表：n=2n=3n=4n=5方法1361015方法22345若采用“秦九韶算法”，则两种算法的乘法运算次数随n的变化见下例：当x接近于0时，的分子、分母都接近0，为避免绝对值小的数作除数，可将原式变形为(四)要避免绝对值小的数作除数例：当x接近于0时，的分子、分母都接近0，为避免绝对值小的数例：当x很大时，的分母接近0，为避免绝对值小的数作除数，可将原式变形为例：当x很大时，的分母接近0，为避免绝对值小的数作除数，可将(五)设法控制误差的传播许多算法具有递推性。递推法运算过程较规律，但多次递推必然导致误差的积累。(五)设法控制误差的传播许多算法具有递推性。递推法运算过程例：求定积分的值.解：直接积分可产生递推公式若取初值例：求定积分的值.解：直接积分可产生递推公式若取初值可得递推公式按公式就可以逐步算出Whathappened?!不稳定的算法！这就是误差传播所引起的危害!注意此公式精确成立，且可得递推公式按公式就可以逐步算出Whathappened?由题设中的递推公式（1）可看出，

的误差扩大了5倍后传给

，因而初值

的误差对以后各步这就造成的计算结果严重失真。计算结果的影响，随着

的增大愈来愈严重。要怎么做才能解决这个问题呢?由题设中的递推公式（1）可看出，的误差扩大了5倍后传给可求得I90.017，按改写后的公式可逐次求得不妨设I9I10，于是由将公式变为可求得I90.017，按改写后的公式可逐次求得不妨设I

I80.019I70.021 I60.024I50.028 I40.034I30.043 I20.058I10.088 I00.182稳定的算法！

在我们今后的讨论中，误差将不可回避，算法的稳定性会是一个非常重要的话题。 I80.019I70.021稳定例：计算积分解：利用分部积分公式例：计算积分解：利用分部积分公式有有从而有递推公式又因为从而有递推公式又因为这说明由初始值E1的误差在计算过程中绝对值会迅速扩大。因此这个算法不可靠！则误差传递规律是这说明由初始值E1的误差在计算过程中绝对值会迅速扩大。因此这但若将公式改为则误差按规律逐步缩小，从而但若将公式改为则误差按规律逐步缩小，从而

所以只要适当选择初值E9，再由公式依次计算E8,E7,···,E1,便可以得到比较精确的结果。因此这个算法是稳定的！所以只要适当选择初值E9，再由公式依次计算E8,E人有了知识，就会具备各种分析能力，明辨是非的能力。所以我们要勤恳读书，广泛阅读，古人说“书中自有黄金屋。”通过阅读科技书籍，我们能丰富知识，培养逻辑思维能力；通过阅读文学作品，我们能提高文学鉴赏水平，培养文学情趣；通过阅读报刊，我们能增长见识，扩大自己的知识面。有许多书籍还能培养我们的道德情操，给我们巨大的精神力量，鼓舞我们前进。人有了知识，就会具备各种分析能力，数值分析11讲义课件数值分析NumericalAnalysis

数值分析

《数值分析》(第2版)

朱晓临主编中国科学技术大学出版社教材《数值分析》(第2版)教材《数值分析》（第5版）

李庆阳，王能超，易大义编著

清华大学大学出版社《数值分析》（第3版）

颜庆津著，

北京航空航天大学

出版社

《NumericalAnalysis》(Ninthed.)RichardBurden,DouglaFaires,

Brooks/Cole,CengageLearning,2011参考书目《数值分析》（第5版）参考书目

微积分

线性代数

常微分方程

算法语言预备知识微积分预备知识一、为什么要学习数值分析？课程简介现实世界的问题可以归结为各种各样的数学问题

—方程求根问题

—解线性方程组的问题

—定积分问题

—常微分方程初值问题

—…….一、为什么要学习数值分析？课方程求根问题在科学计算中常要遇到求解各种方程，例如：高次代数方程x5－3x＋7＝0超越方程方程求根问题在科学计算中常要遇到求解各种方程，例如：高次代数

高次线性方程和超越方程看似简单，但难于求其精确解。对于高次代数方程，由代数基本定理知多项式根的数目和方程的阶相同，但对超越方程就复杂的多，如果有解，其解可能是一个或几个，也可能是无穷多个。高次线性方程和超越方程看似简单，但难于求其精确解。对解线性方程组的问题线性方程组的一般形式解线性方程组的问题线性方程组的一般形式(1)当b≠0时称为非齐次线性方程组，其可能有唯一解、无解或者无穷多个解。当b=0时称为线性齐次方程组，必有零解。(2)由线性代数知识可知：当系数矩阵A非奇异(即detA≠0)时，方程组有唯一解，可用克莱姆法则求解，但它只适合于n很小的情况。(1)当b≠0时称为非齐次线性方程组，其可能有唯一解、无解或克莱姆法则其中克莱姆法则其中例如用克莱姆法则求解一个n阶方程组，要算n+1个n阶行列式的值，总共需要

n!(n-1)(n+1)次乘法。当n充分大时，计算量是相当惊人的。例如用克莱姆法则求解一个n阶方程组，要算n+1个n阶行列比如一个20阶不算太大的方程组，大约要做1021次乘法，这项计算即使每秒1亿次浮点数乘法计算的计算机去做，也要连续工作30万年才能完成。当然这是完全没有实际意义的，故需要寻找有效算法!比如一个20阶不算太大的方程组，大约要做1021次乘法，这项定积分问题对于积分由微积分知识可知：只要找到被积函数f(x)的原函数F(x)，便有下列牛顿—莱布尼兹公式定积分问题对于积分由微积分知识可知：只要找到被积函数f(x)为何要进行数值积分？原因之一：许多形式上很简单的函数，例如已证明它们的原函数不能用初等函数表示成有限形式。为何要进行数值积分？原因之一：许多形式上很简单的函数，例如已原因之二：有些被积函数的原函数过于复杂，例如的一个原函数是原因之二：有些被积函数的原函数过于复杂，例如的一个原函数是要计算f(x)定积分的近似值，上式就不见得方便。要计算f(x)定积分的近似值，上式就不见得方便。原因之三：f(x)以离散数据点形式给出xix0x1…xnyi=f(xi)y0y1…yn原因之三：f(x)以离散数据点形式给出xix0x1…xny常微分方程初值问题一阶常微分方程的初值问题，即例常微分方程初值问题一阶常微分方程的初值问题，即例常微分方程的一般解(解析解)

对一些典型的微分方程(可分离变量方程，一阶线性方程等等)，有可能找出它们的一般解表达式，然后用初始条件确定表达式中的任意常数，这样解即能确定。常微分方程的一般解(解析解)例如

求解解：分离变量得dy=2xdx

积分得y=x2+c

由初值得c=0

故解为y=x2

例如求解解：分离变量得dy=2xdx但是对于求解无法求出一般解！但是对于求解无法求出一般解！二、如何学习数值分析？课程简介1.注意掌握各种方法的基本原理2.注意各种方法的构造手法3.重视各种方法的误差分析4.做一定量的习题5.注意与实际问题相联系二、如何学习数值分析？课1.注意掌握各种方法的基本原理第一章

绪论数值分析第一章

绪论§1数值分析研究的对象与特点一、数值分析研究的对象二、数值分析研究的特点一、数值分析研究的对象二、数值分析研究的特点1.数值分析研究的对象数值分析是计算数学的一个主要部分，计算数学是数学科学的一个分支，它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现.1.数值分析研究的对象数值分析是计算数学的一个主要部分，计算用计算机解决科学计算问题通常经历以下过程实际问题数学模型数值计算方法程序设计上机计算结果应用数学计算数学用计算机解决科学计算问题通常经历以下过程实际问题数学模型数值2.数值分析研究的内容—函数的数值逼近(插值与拟合)—数值积分与数值微分—非线性方程数值解—数值线性代数—常微和偏微数值解，……2.数值分析研究的内容—函数的数值逼近(插值与拟合)数值分析实质上是以数学问题为研究对象，不像纯数学那样只研究数学本身的理论，而是把理论与计算紧密结合，着重研究数学问题的数值方法及理论。数值分析实质上是以数学问题为研究对象，不像纯数学那样只研究数3.数值分析的特点面向计算机，要根据计算机特点设计切实可行的有效算法.

有可靠的理论分析，能任意逼近并达到精度要求，对近似计算要保证收敛性和数值稳定性.3.数值分析的特点面向计算机，要根据计算机特点设计切实可行(3)要有好的计算复杂性，时间复杂性好是指节省时间，空间复杂性好是指节省存贮量，这也是建立算法要研究的问题.(4)要有数值试验，即任何一个算法除了从理论上要满足上述三点外，还要通过数值试验证明是行之有效的.(3)要有好的计算复杂性，时间复杂性好是指节省时间，空间复§2数值计算的误差一、误差来源的分类二、误差分析的重要性三、绝对误差和绝对误差限四、相对误差和相对误差限五、有效数字§2一、误差来源的分类二、误差分析的重要性三、绝对误差和绝对一、误差来源的类型1.模型误差

从实际问题中抽象出数学模型

——模型误差

/*ModelingError*/一、误差来源的类型1.模型误差从实际问题中抽象出数学模2.观测误差

通过测量得到模型中参数的值

——观测误差

/*MeasurementError*/注：通常根据测量工具的精度，可以知道这类误差的上限值。2.观测误差通过测量得到模型中参数的值注：通常根

当数学模型得不到精确解时，要用数值计算方法求它的近似解，由此产生的误差称为截断误差或方法误差3.截断误差

求近似解——方法误差(截断误差)

/*TruncationError*/

当数学模型得不到精确解时，要用数值计算方法求它的近似例如：在微积分中sinx可展开成

但在计算机中计算时，常用前几项来代替，即抛弃了无穷级数的后段，这样就产生了截断误差。

当|x|很小时，常用x代替sinx，其截断误差大约为x3/6。例如：在微积分中sinx可展开成但在计算机中计算时，

由于计算机字长有限，原始数据的输入及浮点数运算过程中都有可能产生误差，这样产生的误差称为舍入误差4.舍入误差

机器字长有限——舍入误差

/*RoundoffError*/由于计算机字长有限，原始数据的输入及浮点数运算过程中在数值分析课程中，主要研究

截断误差

舍入误差

在数值分析课程中，主要研究截断误差二、误差分析的重要性考察如下两个方程组试思考这两个方程组的解的关系？二、误差分析的重要性考察如下两个方程组试思考这两个方程组容易看出系数矩阵完全相同，而常数项矩阵有微小差别，右端系数1.9999变成2.0001，其误差为2.0001-1.9999=0.0002=0.02%容易看出系数矩阵完全相同，而常数项矩阵有微小差别，右端系数1但对应的解为由此看出系数矩阵完全相同，而常数项矩阵有微小差别的方程组，其解竟然相差得很大！解的最大误差=2=200%但对应的解为由此看出系数矩阵完全相同，而常数项矩阵有微小差别据说，美军1910年的一次部队的命令传递是这样的:营长对值班军官:明晚大约8点钟左右，哈雷彗星将可能在这个地区看到，这种彗星每隔76年才能看见一次。命令所有士兵着野战服在操场上集合，我将向他们解释这一罕见的现象。如果下雨的话，就在礼堂集合，我为他们放一部有关彗星的影片。值班军官对连长:根据营长的命令，明晚8点哈雷彗星将在操场上空出现。如果下雨的话，就让士兵穿着野战服列队前往礼堂，这一罕见的现象将在那里出现。连长对排长:根据营长的命令，明晚8点，非凡的哈雷彗星将身穿野战服在礼堂中出现。如果操场上下雨，营长将下达另一个命令，这种命令每隔76年才会出现一次。排长对班长:明晚8点，营长将带着哈雷彗星在礼堂中出现，这是每隔76年才有的事。如果下雨的话，营长将命令彗星穿上野战服到操场上去。班长对士兵:在明晚8点下雨的时候，著名的76岁哈雷将军将在营长的陪同下身着野战服，开着他那“彗星”牌汽车，经过操场前往礼堂。据说，美军1910年的一次部队的命令传递是这样的:三、绝对误差和绝对误差限

定义

设某一量的准确值为x，近似值为x*，则x*与x之差叫做近似值x*的绝对误差(简称误差)，记为？判断题：绝对误差是误差的绝对值．三、绝对误差和绝对误差限定义设某一量的准确值为x，绝对误差的性质(1)绝对误差e(x*)可正可负(2)|e(x*)|的大小标志着x*的精确度(3)绝对误差e(x*)未知绝对误差的性质(1)绝对误差e(x*)可正可负定义

若指定一个适当小的正数

，使

有时用

表示近似值x*的精度或准确值的所在范围。则

称为近似值

x*的绝对误差限。定义有时用表示绝对误差限的性质(1)在实际问题中，绝对误差一般是有量纲的，绝对误差限也是有量纲的。例如，测得某物体的长度为5m，其误差限为0.01m绝对误差限的性质(1)在实际问题中，绝对误差一般是有量纲的，(2)绝对误差限是正的，有无穷多个。若已知

是绝对误差限，由于则比大的任意正数均是绝对误差限。(2)绝对误差限是正的，有无穷多个。若已知是绝对误差

思考题：设有两个温度计，其一测量1000℃时的绝对误差限为5℃，而另一个测量100℃时的绝对误差限为1℃。问：哪一个温度计更精确？四、相对误差和相对误差限思考题：设有两个温度计，其一测量1000℃时的绝对误答：虽然后者绝对误差限的数值较小，但第一种温度计更为精确。

决定一个量的近似值的精确度除了要看绝对误差的大小外，还要考虑到该量本身的大小。答：虽然后者绝对误差限的数值较小，但第一种温度计更为精确。定义:

绝对误差与准确值之比称为x*的相对误差.定义:绝对误差与准确值之比称为x*的相对误差.(2)由于准确值x未知，故实际问题中，当||较小时，常取注(1)相对误差是个无量纲量,值小者精度高。(2)由于准确值x未知，故实际问题中，当|当||较小时，可用下式计算定义若指定一个适当小的正数，使则称

为近似值x*的相对误差限当||较小时，可用下式计算定义若指定一个适

当x有很多位数字时，常按照“四舍五入”原则取前几位数字作为x的近似值例：设x=π=3.1415926···取x1*=3作为π的近似值，则五、有效数字当x有很多位数字时，常按照“四舍五入”原则取前几位数取x2*=3.14作为π的近似值，则取x3*=3.1416作为π的近似值，则它们的误差都不超过某一数位的半个单位。取x2*=3.14作为π的近似值，则取x3*=若近似值x*的绝对误差限是某一位的半个单位，该位到x*的第一位非零数字一共有n位，则称近似值x*有n位有效数字，或说x*精确到该位。定义：有效数字若近似值x*的绝对误差限是某一位的半个单位，该位到x*的第一若将准确值x

的近似值x*表示成标准形式而其误差限则说近似值x*具有n

位有效数字。这里n

为正整数，m

为整数，每个均为0,1，…，9中的一个数字，注：有效数字的等价定义若将准确值x的近似值x*表示成标准形式而其误差限则说近似值定理1.1

若

x*具有n位有效数字，则其相对误差限满足反之，若x*的相对误差限满足则近似值x*至少具有n位有效数字。定理1.1若x*具有n位有效数字，则其相对反之，若x*例：为使的相对误差小于0.001%,至少应取几位有效数字？解：假设*取到n位有效数字，则其相对误差上限为要保证其相对误差小于0.001%，只要保证其上限满足已知a1=3，则从以上不等式可解得n>6log6，即n6，应取*=3.14159。例：为使的相对误差小于0.001%,至少应取几位有效数§3避免误差危害的若干原则(一)要避免相近两数相减(二)要防止大数“吃掉”小数(三)注意简化计算步骤，减少运算次数(四)要避免绝对值小的数

作除数(五)设法控制误差的传播§3(一)要避免相近两数相减(二)要防止大数“吃掉”小数(一)要避免相近两数相减的值。例:求当x=1000，y的准确值为0.01580.

(1)直接相减(2)若将原式改写为则y=0.01581(一)要避免相近两数相减的值。例:求当x=1000，几种经验性避免方法：很小，几种经验性避免方法：很小，(二)要防止大数“吃掉”小数例：用单精度计算的根。精确解为算法1：利用求根公式在计算机内，109存为0.11010，1存为0.1101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。即1的指数部分须变为1010，则：1=0.00000000011010，取单精度时就成为：109+1=0.100000001010+0.000000001010=0.100000001010(二)要防止大数“吃掉”小数例：用单精度计算算法2：先解出

再利用注：求和时从小到大相加，可使和的误差减小。算法2：先解出注：求和时从小到大相加，可使和的误差减小。例：在5位浮点十进制计算机上，计算y=54321+0.4+0.3+0.4解：若按从左到右的顺序进行计算，后三位在对阶过程变为例：在5位浮点十进制计算机上，计算y=54321

后三个数都在对阶过程中变为零，得出含有较大误差的结果y=54321。但若按从右到左的顺序进行计算，后三位在对阶过程变为

这种算法避免了大数“吃掉”小数！后三个数都在对阶过程中变为零，得出含有较大误差的结果一般地，有如下原则若干数相加，采用绝对值较小者先加的算法，结果的相对误差限较小一般地，有如下原则(三)注意简化计算步骤，减少运算次数，避免误差积