数值分析课件 第1章 绪论_0

------------------------------------------- --------------------最新资料推荐------------------------ ------------------------------

数值分析课件 第1章 绪论

第一章 绪论 上世纪中叶诞生的计算机给科学、 工程技术和人
类的社会生活带来一场新的革命。
它使科学计算平行于理论分析和实验研究， 成为人类探索未知
科学领域和进行大型工程设计的第三种方法和手段。
在独创性工作的先行性研究中， 科学计算更有突出的作用。
在今天， 熟练地运用电子计算机进行科学计算， 已成为科学工
作者的一项基本技能。
然而， 科学计算并不是计算机本身的自然产物， 而是数学与计
算机结合的结果， 它的核心内容是以现代化的计算机及数学软件为
工具， 以数学模型为基础进行模拟研究。
近年来， 它同时也成为数学科学本身发展的源泉和途径之一。
1 数值分析的研究对象与特点 数值分析是计算数学的一个主
要部分， 计算数学是数学科学的一个分支， 它研究用计算机求解各
种数学问题的数值计算方法及其理论与软件实现。
一般地说， 用计算机解决科学计算问题， 首先需要针对实际问
题提炼出相应的数学模型， 然后为解决数学模型设计出数值计算方
法， 经过程序设计之后上机计算， 求出数值结果， 再由实验来检
验。
概括为 由实际问题的提出到上机求得问题的解答的整个过程
都可看作是应用数学的任务。



1 21

如果细分的话， 由实际问题应用有关科学知识和数学理论建立
数学模型这一过程， 通常作为应用数学的任务， 而根据数学模型提
出求解的数值计算方法直到编出程序上机计算出结果， 这一过程则
是计算数学的任务， 即数值分析研究的对象。
因此， 数值分析是寻求数学问题近似解的方法、 过程及其理论
的一个数学分支。
它以纯数学作为基础， 但却不完全像纯数学那样只研究数学本
身的理论， 而是着重研究数学问题求解的数值方法及与此有关的理
论， 包括方法的收敛性， 稳定性及误差分析； 还要根据计算机的
特点研究计算时间最省（或计算费用最省）的计算方法。
有的方法在理论上虽然还不够完善与严密， 但通过对比分析，
实际计算和实践检验等手段， 被证明是行之有效的方法也可采用。
因此数值分析既有纯数学高度抽象性与严实际问题数学模型计
算方法程序设计计算结果密科学性的特点， 又有应用的广泛性与实
际试验的高度技术性的特点， 是一门与使用计算机密切结合的实用
性很强的数学课程。
在电子计算机成为数值计算机的主要工具以后， 则要求研究适
合计算机使用的， 满足精确要求， 计算时间省的有效算法及其相关
的理论。
在实现这些算法时往往还要根据计算机的容量、 字长、 速度等
指标， 研究具体的求解步骤和程序设计技巧。
有的方法在理论上虽还不够严格， 但通过实际计算、 对比分析

----------------------------------- ----------------------------最新资料推荐---------------- --------------------------------------

等手段，证明是行之有效的方法， 也应采用。
这些就是数值分析具有的特点， 概括起来有四点：
第一， 面向计算机， 要根据计算机特点提供切实可行的有效
算法。
即算法只能包括加、 减、乘、 除运算和逻辑运算， 这些运算
是计算机能直接处理的运算。
.第二， 有可靠的理论分析， 能任意逼近并达到精确要求， 对
近似算法要保证收敛性和数值稳定性， 还要对误差进行分析。
这些都建立在相应数学理论的基础上。
第三， 要有好的计算复杂性， 时间复杂性好是指节省时间， 空
间复杂性好是指节省存储量， 这也是建立算法要研究的问题， 他关
系到算法能否在计算机上实现。
第四， 要有数值试验， 即任何一个算法除了从理论上要满足
上述三点外， 还要通过数值试验证明是行之有效的。
根据数值分析 课程的特点， 学习是我们首先要注意掌握方法
的基本原理和思想， 要注意方法处理的技巧及其与计算机的结合，
要重视误差分析、 收敛性及稳定性的基本理论； 其次， 要通过例
子， 学习使用各种数值方法解决实际计算问题； 最后， 为了掌握
本课的内容，还应作一定数量的理论分析与计算练习。
由于本课内容包括了微积分、 代数、 常微分方程的数值方法，
读者必须掌握这几门课的基本内容才能学好这门课。

3 21

2 数值计算的误差 2.1 误差来源与分类 我们算出数学模型
的近似解和一个物理量的真的值往往不相等， 它们之差称为误差。
用电子计算机进行解决实际问题的数值计算，误差是不可避免的。
引起数值结果中的误差的原因是多方面， 通常来自 固有误差和
计算误差， 如下面所示：
模型误差：
用数学方法解决实际问题， 首先必须建立该问题的数学模型。
即把实际问题经过抽象， 忽略一些次要的因素， 简化成一个确
定的数学问题。
数学模型只是对实际问题的一种近似、 一种粗糙的描述， 因而
它与实际问题或客观现象之间必然存在误差， 这种误差称为模型误
差。
这样的误差常常是可以忽略不计的。
例如在经典力学问题中， 我们常常忽略相对论效应。
但是， 如果这种误差不可忽略，说明数学模型选择得不好。
那么不论数值计算多么精确， 其结果都将存在不可忽略的误差。
观测误差：
数学问题中总包含一些参量（或物理量， 如电压、 电流、 温
度、 长度等）， 它们的值（输入数据） 往往是由观测得到的。
而模型误差 观测误差 误差来源 计算误差 固有误差 截断误差
舍入误差 观测的误差是难以避免的， 由此产生的误差称为观测误
差。

--------------------------------------- ------------------------最新资料推荐-------------------- ----------------------------------

由于观测误差通常具有随机的性质， 所以想用分析的方法来估
计它们的影响常常是一件非常困难的事。
由于固有误差的产生往往涉及各专业知识及实验手段， 所以它
不是数值分析所研究的内容。
截断误差：
数值计算的本质就是用有限的过程、 离散的数据来近似（刻画）
无限的过程、连续的量。
这样对一些具有连续量、 无限计算量的数学问题的求解的过程
中， 我们不得不对一些连续的量进行有限的离散化近似， 不得截去
无限的计算过程， 只进行有限次的计算。
如求一个收敛的无穷级数之和， 必须截去该级数后面的无穷多
项， 而用前面有限项的部分来近似代替， 于是产生了有限过程代替
无限过程的误差， 称为截断误差， 这是计算方法本身所出现的误差，
所以也称为方法误差。
例如
小时， 可以用212x作为cos x近似值。
由交错级数判敛的莱布尼兹（Leibniz） 准则，它的截断误差的
绝对值不超过424x。
有限过程代替无限过程的误差和计算量取决于其方法的收敛性
及收敛速度。
舍入误差：

5 21

， 当 x 很

计算中遇到的数据可能位数很多或是无穷小数， 如
21.41421356=长的限制， 无穷小数和位数很多的数必须舍入成一定
的位数（机器字长）。
舍入方法：
a) 舍入机。
采用四舍五入的办法， 如将
截断机。
采用截尾的办法。
如 2在八位字长的截断机里取成1.4142135。
这样产生的误差称为舍入误差。
少量的舍入误差是微不足道的， 但在计算机是作了成千上万次
运算后， 舍入误差的累积有时可能是十分惊人的。
它取决于方法的稳定性。
如果方法能够累积大量的误差， 此算法是不稳定的， 反之为稳
定算法。
如算法 1：
受机器字1nnyany=+。
舍为1.4142136； b)
显然是不稳定的； 算法 2：
11nnyayn=+， 研究计算结果的误差是否满足精度要求就是误差
估计问题， 本书主要讨论算法的截断误差与舍入误差， 而截断误差
将结合具体算法讨论。
为分析数值运算的舍入误差， 先要对误差基本概念做简单介绍。

--------- -------------------------------------------------- ----最新资料推荐---------------------------------------- --------------

2.2 绝对误差与相对误差 定义 1 设x为准确值，称exx= 为
近似值的绝对误差， 简称误差。
注意这样定义的误差差为正时近似值偏大， 叫强（赢） 近似
值； 当绝对误差为负时近似值偏小， 叫弱（亏） 近似值。
通常我们不能算出准确值x， 当然也不能算出误差情况估计出
误差的绝对值不超过某正数就是误差绝对值的一个上界。
误差限， 它总是正数。
一般情形xx， 即*x 为 x的一个近似值，***e 可正可负， 当
绝对误*e 的准确值， 只能根据测量工具或计算* ， 也* 叫做近似
值的******xxx*x+。
这个不等式有时也表示为我们把近似值的误差*x=。
*e 与准确值x的比值 *xx *exx=称为近似值*x 的相对误差，
记作*re 。
在实际计算中， 由于真值x总是不知道的，
*****rexx***rex****()eee xx==通常取exx==作为*x 的相对误差，
条件是e=较小， 此时 * 2** 2)*******( )e(()1()exxxx xx xeex=
是*r 的平方项级， 故可忽略不计。
相对误差也可正可负， 它的绝对值上界叫做相对误差限， 记
作*r ， 即**r*x =。
注 1：
1） 绝对误差是有量纲的量， 它与所研究问题的背景有关， 因

7 21

此我们不能单单从绝对误差值的大小来判断计算结果的精度， 还必
须考虑到实际问题的应用背景； 2） 通过引入相对误差， 我们可
比较不同算法， 不同应用问题的计算精度， 这一点是无法通过比较
绝对误差做到的， 因为不同值的东西， 在量上是无法比较的。
正如我们不能通过比较两个商品的使用价值来判定商品的真贵，
只能通过抽象的价值。
2.3 有效数字 有效数字是近似值的一种表示法， 它既能表示
近似值的大小， 又能表示其精确程度。
定义 2 若近似值个单位， 该位到则称近似值有n位有效数字。
在科学记数法中， 将近似值式为
10,aa i间的整数。
按照定义 2， 近似值字当且仅当 1xx*x 的误差限是某一位的半
*x 的第一位非零数字共有n位，*x 写成规格化形
9 之*x 有n位有效数
为 0 到
其中m为整数，
-1) (1,2,)i=*102m n
(1-2) 因此在m相同的情形下， n越大则误差越小，亦即一个近似
值的有效位数越多其误差限越小。
例 1 按四舍五入原则写出下列各数具有5 位有 效 数 字 的
近 似 数 :187.9325,0.03785551, 8.000033, 2.7182818。
按定义， 上述各数具有 5 位有效数字的近似数分别是 187.93,
0.037856, 8.0000, 2.7183 8.000033x =的 5 位有效数字的近似数
是 8.0000 而不是 8， 因 为 8 只有 1 位有效数字。
例 2 重力 常 数 g ， 如果以120.980 10gm s；若 以220.980

------ -------------------------------------------------- -------最新资料推荐------------------------------------- -----------------

10gkm s， 它们都具有 3 位有效数字， 因为按第一种写法
19.80102按第二种写法 10.00980102他们虽然写法不同， 但都具有
3 位有效数字。
至于绝对误差限， 由于单位不同结果也不同，1110,22差都是
0.0059.80=例 2 说明有效位数与小数点后有多少位数有关。
然而， 从(2-2)可以得到具有n位有效数注意2m s 为单位，
2km s 为 单 位 ，21 31102g=， 52 3 1102g=， *122*25210m skm
s==， 而相对误*0.0000050.00980r =。
字的近似数*x ， 其绝对误差限为*1102m n=，在m相同的情况
下， n越大则10m n效位数越多， 绝对误差限越小。
越小， 故有关于一个近似数的有效位数与其相对误差的关系，
有下面的定理 定理 1.1 设近似数(1)若为 1ra(2) 如果 1ra+则证
明 由(1-1)可得到 1*110ax*x 有n位有效数字时 *0.5 10rax*x 具
有规格化形式(1-1)， *x 具有n位有效数字， 则其相对误差限
*11102n + (1-3) *11102(1)n + (1-4) *x 至少具有
n位有效数字。
11(1) 10mma+ 当*11*11110102m nnmxxa +== 即证(1)。
反之， 由 1***r11110(1) 102(1)1102nmm nxxxaa +=++= 即说
明定理说明， 有效位数越多， 相对误差限越少。
例 3 要使 20 的近似值的相对误差限小于0.1%， 要取几位有
效数字？ 设 取 n 位 有 效 数 字 ，1102a4n = ， 就有 *30.125

9 21

10即只要对 20的近似值取 4 位有效数字， 其相对 误 差 就 小
于 0.1% 。
此 时 由 开 方 表 得204.472。
3 数值计算中误差的传播 *x 具有n位有效数字， (2)得证。
由 定 理 1.1 ，*r11n +。
由于 ， 知14a = ，故只要取310r0.1%= 3.1 基本运
算中的误差估计 数值运算中误差传播情况比较复杂， 估计起来比较
困难。
本节所讨论的运算是四则运算与一些常用函数的计算。
由微分学， 当自变量改变（误差） 很小时，函数的微分作为
函数的改变量的主要线性部分可以近似函数的改变量， 故可以利用
微分运算公式可导出误差运算公式。
设数值计算中求得的解与参量（原始数据）
为 ( ,yf x x=参量的误差必然引起解的误差。
设的近似值分别为
较小时， 解的绝对误差为
， 相应的解为
*1*2*n12*1*2*n*1*2*n*i1*1*2*n**i1(,,,)(
可微， 则当数据误差
,f x
(yf x x=假设 f 在点12(,,,x
有关， 记
x,,)(,,,)(,,,)()(,,,)()nniiiniif x xxxdf x xxf x xxxxxf x xxe
-5) 其相对误差为
***r ***1*2*n**i*1*2*n1*1*2*n*i*r*i*1*2*n1()()(ln )(,,,)()(,
,,)(,,,)()(,,,)niiniie ye ydfyf x xxe xxf x xxf x xxxe xxf x

------------ -------------------------------------------------- -最新资料推荐------------------------------------------- -----------

-6) 注 2：
函数值的绝对误差等于函数的全微分，自变量的微分即为自变量
的误差； 函数值的相对误差等于函数的对数的全微分。
将式(1-5)及(1-6)中的差限 和特别地， 由式(1-5)及(1-6)可
得和、 差、 积、商之误差及相对误差公式
和*( )re 分别换成误r ， 求和的
各项变成绝对值。
1212122 1x e1 2x e11221222()()1()e xxeee x xxxeeexxx=+
(1-7) 2121122()()()rrrrrrrrrxxe xxeexxxxe x
-8) 由以上的式(1-7)及(1-8)
可得出 (e xx2112122)()()rrrrrrrrrreeeee x
xeeeexeeeeex=+=++=+ (1-9) 因此， 和、 差的误差限不超过各数
的误差限之和， 积、 商的相对误差限不超过各数的相对误差限之和。
例 4 设yx=， 求 y的相对误差与 x的相对误差之间的关系。
解 由式(1-9)得 ( )(ln)re ydx==所以地，x的相对误差是x
的相对误差的一半。
例 5 设0x ， x的相对误差为 ， 求ln x的绝对误差。
n(ln )( )nrndxne x= nx 的相对误差是x的相对误差的n倍 ，
特别解 由于( )e x( )re xx ==， 即 ( )e xx=， 所以 ( )e
x(ln )(ln )( )rexdxe xx=== 例 6 已测得某场地长为l的值为d 的
值 为d=*|| 0.1ddm， 试求面积 s相对误差。

11 21

)80 ,(dl( )于是绝对误差限 ( )相对误
差限 ( )*|s 3.2 误差定性分析 *100 ,| 0.2lmm=宽，*80 ,m已 知
*|ll的绝对误差与ld=解 因,,,sdld===由(1-5)知
其
**)110smlmd====， 而 ()**0.2 ,0.1
中
,ldm==
()()()2*800.21100.127sm+= ( )*****( )s270.31%|8800rssl d===
数值运算中的误差分析是一个很重要而复杂的问题， 上节讨论了不
精确数据运算结果的误差限， 他只适用于简单情形， 然而一个工程
或科学计算问题往往要运算千万次， 由于每步运算都有误差， 如果
每步都做误差分析是不可能的， 也不科学， 因为误差积累有正有负，
绝对值有大有小， 都按最坏情况估计误差限得到的结果比实际误差
大得多， 这种保守的误差估计不反映实际误差积累。
考虑到误差分布的随机性， 有人用概率统计方法， 将数据和运
算中的舍入误差视为适合某种分布的随机变量，然后确定计算结果的
误差分布， 这样得到的误差估计更能接近实际， 这种方法称为概率
分析法。
20 世纪 60 年代以后对舍入误差分析提出了一些新的方法，
较重要的有以下两种：
1.向后误差分析法是把新算出的量由某个公式表达， 它仅含基
本的算术运算， 如假定
算出量 (xg a=， ，若
出结果为
是前面已算出的量或原始数据， 新
的摄动为i ， 使得由浮点运算得
则可根据xxi 的界由摄动理论估

------------------------------------------------ ---------------最新资料推荐----------------------------- -------------------------

计最后舍入误差的界， 威克逊（Wilkinson） 将这种方法fl应用于
数值代数（矩阵运算） 的误差分析， 取得较好的效果。
2.区间分析法是把参加运算的数 , , ,看成区间量
根据区间运算规则求得最后结果的近似值及误差限。
例如，似数为, ， 由 于x[,],xX y +=若计算*zx y=（*为运
算符号）， 由[ , ][,]z zzz zz=+,则 z为所求近似值， 而 z则为
误差限。
上面简略介绍了误差分析的几种方法， 但都不是十分有效的，
目前尚无有效的方法对误差作出定量估计。
为了确保数值计算结果的正确性， 首先需对数值计算问题作出
定性分析， 为此本节讨论以下三个问题。
3.3 病态问题与条件数 对一个数值问题本身 如果输入数据
又微小扰动（即误差）， 导致输出数据（即问题解）都, x y
的近]+*X Y=,,y[,则Y= , =Z相对误差很大， 这就是病态问题， 例
如计算函数值 ( )f x 时， 若x有扰动x?，函 数 值f x*xxx? =，
其相对误差 为xf x*()的 相 对 误 差 为*( )f x()( )f x。
相对误差比值 *( )f x()( )( )f x( )f xpf xxxf xCx?= (1-10)
pC 称为计算函数值问题的条件数。
自变量相对误差一般不会太大， 如果条件数pC 很大， 将引起
函数值相对误差很大， 出现这种情况的问题就是病态问题。
例如，( )f xx=则有差 可 能 放 大 n 倍 。

13 21

如(1.02)1.24f， 若取x=误差为 2%， 函数值相对误差为 24%，
这是问题可以认为是病态的。
一般情况条件数C 越大病态越严重。
npCn= ， 他表示相对误10n =，有*1,1.02x=自变量相对
(1)1f= ，10pC 就认为是病态的，p其他计算问题也要分析是否病态。
例如解线性方程组， 如果输入数据有微小误差引起解的巨大误
差， 就认为是病态方程组， 我们将在第5 章用矩阵的条件数来分
析这种现象。
3.4 算法的数值稳定性 数值运算的误差分析是个很重要而复
杂的问题， 对于一个工程或科学计算问题往往要运算千万次， 由于
几乎每一步运算都会产生误差， 如果每步都做误差分析是不可能的，
也不科学， 因为误差累积有正有负， 绝对值有大有小， 很可能有
一部分相互抵消。
如果每一步都按最坏情况估计误差不反映实际情况。
近年来， 虽然发展起来一些新的误差估计办法， 如概论统计法，
事后误差估计法及区间分析法等等。
但这些方法目前也仅仅在特定的一些问题上取得了较为满意的
结果， 还正在发展中。
对于计算步骤较多的算法， 一般不宜采用上述的一些定量估计
的办法。
而是采用定性分析的办法， 即讨论算法的数值稳定性。
定义 3 一个算法如果输入数据有误差， 而在计算过程中舍入

--------------- ------------------------------------------------最新 资料推荐---------------------------------------------- --------

误差得到控制， 则称此算法是数值稳定的， 否则称此算法是不稳定
的。
在一种算法中， 如果某一步有了 绝对值为的误差， 而以后各
步计算都准确地进行， 仅由 所引起的误差的绝对值， 始终不超过 ，
就说算法是稳定的。
对于数值稳定性的算法， 不用做具体的误差估计， 就认为其结
果是可靠的。
而数值不稳定的算法尽量不要使用。
1nIe=计误差。
由分部积分可得计算若计算出值。
要算出展开部分和 ( 1)1( 1)2!并取7k = ， 用 4 位小数计算，
则得e=例 7 计算并估nI 的递推公式
-11) 0I ， 代 入 (1-11)， 可
逐次求出的0I 就要先计算1e， 若用泰勒多项式
， 10.3679，截断误差14710.3679108!4R。
计 算过程中小数点后第 5 位的数字按四舍五入原则舍入， 由
此产生的舍入误差这里先不讨论。
当初始值取为
（A）
计算结果见表 1-1 的误差
列。

15 21

时， 用(1-11)递推的计算公式为 0I
。
结果是有影响的。

用近似0I 产生的0I就是初始误差， 它对后面计算表 1-1
用(A)算) (用(B)算) (用(A)算) (用(B)
算) 0 1 2 3 4 0.6321 0.3679 0.2642 0.2074 0.1704 0.6321 0.3679
0.2643 0.2073 0.1708 5 6 7 8 9 0.1480 0.1120 0.2160 -0.7280
7.552 0.1455 0.1268 0.1121 0.1035 0.0684 从表中可以看到矛盾。
实际上， 有积分估值得 出现负值， 这与一切
相100111001(min )11(max )1xnnxxnxex dxIeex dxn==+ (1-12) 因
此， 当n较大时， 用的。
这里的计算公式与每步计算都是正确的，那么是什么原因使计算
结果错误呢？ 主要就
差有误差
近似nI 显然是不正确是初值步计算的误
， 由此引起以后各满足关系 1,2,)n
容易推得 0( 1)E 导致第n步的误差扩大 !
n!nnEn E= ， 由此看出：
误差倍， 当n较大时， 误差将淹没真值， 因此用似
nI 显然是不正确的， 这种递推公式不宜采用。
1E=例如，8n = ， 若40102， 则 808!2EE= 这就说明是数值
不稳定的。
现在换一种计算方法。
由(1-12)取完全不能近似8I 了。
近
它表明公式(A)9n = ， 有 1911010I 我们粗略取
19911()0.06842 1010eII+==， 然后将公式(1-11)倒过来算， 即由
公式为 *9I 算出（B）

--- -------------------------------------------------- ----------最新资料推荐---------------------------------- --------------------

计算结果见
表 1-1 的误 差 不 超 过1nEEn*9 E 较大， 但由于误差逐步缩小，
故可用*nI 列。
我们发现。
记E*0I 与*nI0I 的则410*nnI=，*0*!=，*0E 比*nE 缩小了 ! n
倍， 因此， 尽管*nI 近似nI 。
反之， 当用方案(A)计算时， 尽管初值当准确， 由于误差传播
是逐步扩大的， 因而计算结果不可靠。
此例说明， 数值不稳定的算法是不能使用的。
在例 7 中， 算法(B)使数值稳定的， 而算法(A)是不稳定的。
数值不稳定现象属于误差危害现象， 如何防止误差危害下面将
进一步讨论。
4 数值计算中应注意的几个问题 相数值计算中首先要
分清问题是否病态和算法是否数值稳定， 计算时还应尽量避免误差
危害， 防止有效数字的损失， 下面给出若干原则：
1． 要避免除数绝对值远远小于被除数绝对值的除法。
因为 2( )( )( )yxye xxe yey= 故当 y例 8 线性方程组 x时，
舍入误差可能增大很多。
的准确解为
201900x=121999980.50000125,0.99999539999 9199999x=== 现在四
位浮点十进制数（仿机器实际计算） 下用消去法求解， 上述方程写

17 21

成 程， 则
出现用小的数除大的数， 得到 411100.100010
由此解出 87 0, 10 0.1000
若用（除第一方程减第二方1166210
0.1000,1显然严重失真。
若反过来用第二个方程消去第一个方程中含66210 0.100010
由此求的相当好的近似解
0.5000,xx=2． 要避免两相近数相减 两数之差uxy= 的相对误差为
1x 的项， 则避免了大数被小数除， 得到 1111210 0.201910
。
( )e x( )e y( )()rre ue xyxy== 当 x与 y很接近时， u的
相对误差会很大， 有效数字位数将严重丢失。
例如，532.52y =都具有五位有效数字， 但只有两位有效数字。
这说明必须尽量避免出现这类运算。
最好是改变计算方法， 防止这种现象产生。
先举例说明。
例 9 求1610xx+ = 的根。
解 1863x = +， 8630.8 100.794 10x= *2x 只有一位有效数
字。
若改用 532.65x =，0.13xy=2*220.006 10x== 12118630.627
1015.94863x= =+ 具有 3 位有效数字。

------- -------------------------------------------------- ------最新资料推荐-------------------------------------- ----------------

例 10 计算表 ） 。
由于cos2710 (1A=710 (1cos2 )A=（用四位数学用
0.9994cos2 )=， 直接计算 710 (10.9994)36 10= =。
x只有一位有效数字， 若利用21cos2sin2x=，则 A具有三位
有效数字（这里sin1此例说明， 可通过改变计算公式避免或减少有
效数字的损失。
类似的如果则 727310 (1cos2 )2 (sin1 )= 106.13 10）。
=== 0.01751x 和2x 很接近时，1122lglglgxxxx= 用右边算式
有效数字就不损失。
当 x很大时， 111xxxx+ =+ + 都 用 右 边 算 式 代 替 左
边 。
一 般 情 况 ， 当( )( )*f xf x时， 就用泰勒展开 ( )f
取右端的有限项近似左端。
如果无法改变算式， 则采取增加有效位数进行运算； 在计算机
上则采用双倍字长进行运算， 但是增加计算时间和多占内存单位。
3． 要防止大数吃掉 小数 在运算中参加运算的数有时数量级
相差很大， 而计算机位数有限， 如不注意运算次序就可能出现大数
吃掉 小数的现象， 影响计算结果的可靠性。
例 11 在五位十进制计算机上， 计算 100052492ii=把运算的
数写成规格化形式 1A=+， 其中0.10.9i 1000515.2492 10iiA==+ 由
于计算机内计算时要对阶， 若取阶时0.000009 10示为机器零， 因

19 21

此 0.9i =， 对5i =， 在五位的计算机中表55555.2492 100.000009
100.000009 105.2492 10A?==+++ （符号的。
这是由于运算中大数吃掉 小数造成的。
例 12 求二次方程根。
解 利用因式分解容易求出此方程的两个根为1210 ,1xx== ，
但若用求根公式， 则得 9101(10x=?=表示机器中相等）。
结果显然是不可靠0.9i =299(101)100xx++= 的99291)4 10 2+ +
若用 8 位小数的计算机运算， 由于对阶， 有 9101+
9929910(1010 ,=1)4 10 0= ， 结果显然是错误的。
10+x 这样求得为避免这种情形出现， 也可采用改变计算公式的
方法。
如将式 9101(10x=912x92921)4 10 2+ + 改变成 9299292
10101(101)4 10 x=+ ++ 则有 92992 10+11010x== 此结果是精确的。
4 要用简化计算， 减少运算次数， 提高效率 求一个问题的数
值解法有多种算法， 不同的算法需要不同的计算量， 如果能减少运
算次数， 不但可节省计算机的计算时间， 还能减少舍入误差累积。
这是数值计算必须遵从的原则， 也是数值分析 要研究的重要内
容。
例如， 若选用级数 111ln2 1234的前n项部分和来计算ln2的
近似值， 截断误差1。
如果要求误差不超过
即要前十万项求和。
为1n +510， 则510n ，

--------------------------------------------- ------------------最新资料推荐-------------------------- ----------------------------

这样做不仅计算量大， 而且舍入误差的累积将使有效数字丢失
严重。
如果利用级数 224111ln2 (1x)13521n
计算， 当13x = 时有
来
取
前 5 项之和作为近似值， 产生的截断误差为 211(3 11 913
9211(1311 992113 11 919显然， 第二种算法比第一种算法有效。
又如计算多项式 ( )nnnP xa xax=+567525451)15 91)91112 11
的值。
若直接按上式计算， 共需作(2次乘法与n次加法。
若按秦九韶算法 10 (1,2,,1,0)( )nnkkknuauxuaknnP
计算， 即将前式改写成如下形式 ( ){nP xax
a=+则只需作n次乘法和n次加法。






21 21