数值分析总复习提纲教材
数值分析总复习提纲
数值分析课程学习的内容看上去比较庞杂,不同的教程也给出了不同的概括,但总的来说无非是误差分析与算法分析、基本计算与基本算法、数值计算与数值分析三个基本内容。在实际的分析计算中,所采用的方法也无非是递推与迭代、泰勒展开、待定系数法、基函数法等几个基本方法。
一、误差分析与算法分析
误差分析与算法设计包括这样几个方面: (一)误差计算 1、截断误差的计算
截断误差根据泰勒余项进行计算。
基本的问题是
(1)1
()(01)(1)!
n n f x x n θεθ++<<<+,已知ε求n 。
例1.1:计算e 的近似值,使其误差不超过10-6。
解:令f(x)=e x ,而f (k)(x)=e x ,f (k)(0)=e 0=1。由麦克劳林公式,可知
211(01)2!!(1)!
n x x
n x x e e x x n n θθ+=+++++<<+
当x=1时,1
111(01)2!
!(1)!
e e n n θθ=+++
++
<<+
故3
(1)(1)!(1)!
n e R n n θ=<++。
当n =9时,R n (1)<10-6,符合要求。此时, e≈2.718 285。
2、绝对误差、相对误差及误差限计算
绝对误差、相对误差和误差限的计算直接利用公式即可。 基本的计算公式是:
①e(x)=x *-x =△x =dx
② *()()()ln r e x e x dx
e x d x x x x
====
③(())()()()e f x f x dx f x e x ''== ④(())(ln ())r e f x d f x =
⑤121212121122121122((,))(,)(,)(,)()(,)()x x x x e f x x f x x dx f x x dx f x x e x f x x e x ''''=+=+ ⑥121212((,))
((,))(,)
f x x f x x f x x εδ=
⑦ x
ε
δ=
注意:求和差积商或函数的相对误差和相对误差限一般不是根据误差的关系而是直接从定义计算,即求出绝对误差或绝对误差限,求出近似值,直接套用定义式
()()r e x e x x =
或x
ε
δ=, 这样计算简单。
例1.2:测得圆环的外径d 1=10±0.05(cm),内径d 2=5±0.1(cm)。求其面积的近似值和相应的绝对误差限、相对误差限。
解:圆环的面积公式为: 22
12()4
S d d π
=
-
所以,圆环面积的近似值为 222(105)58.905()4
S cm π
=
-≈
由上述讨论,面积近似值的绝对误差限为
112211222()(2()2())(()())
4
2
(100.0550.1)2
1.57()
S d d d d d d d d cm π
π
εεεεεπ
≤
+=
+=
?+?≈
相对误差为
() 1.57
()100% 2.7%58.905
S S S εδ==?≈
相对误差要化成百分数。
3、绝对误差、相对误差、有效数字的关系计算
绝对误差、相对误差、有效数字的关系依据如下结论讨论: ①如果一个数
*123
11
10.(0)n n n x a a a a a a a -+=±≠
其近似值
123
10.n n x a a a a a -=±
是对x*的第n+1位进行四舍五入后得到的,则x 有n 位有效数字,且其绝对误差不超过
1
102
n -? ,即 1
*102
n x x --≤? 。
②如果一个数
*123
1110.10(0)m n n n x a a a a a a a -+=±?≠
的近似值
123
10.10m n n x a a a a a -=±?
是对x*的第n+1位进行四舍五入后得到的,则x 有n 位有效数字,且其绝对误差不超过
1
102
m n -? ,即 1
*102
m n x x --≤?。
③设12310.10m n n x a a a a a -=±?是x*的具有n 位有效数字的近似值,则其相对误差限为
11
1102n a δ-=?
反之,若x 的相对误差限111
102(1)
n a δ-=
?+
则x 至少具有n 位有效数字。
例1.3
的近似值,使其绝对误差不超过31
102
-?。
解:因为12<<
所以,化成12310.10m n n x a a a a a -=±?的形式,有11,1a m ==。
而31411
101022
ε--=?=?,
所以,由定理2,n=4,
所以近似值应保留4位有效数字。
1.732≈。
例1.4
的近似值的相对误差不超过410-,应取几位有效数字?
(5%)
解:设取n 个有效数字可使相对误差小于410-,则 1411
10102n a --?<,
而34≤≤,显然13a =,此时,
114111*********n n a ---?=?, 即141
10106
n --?<, 也即561010n ?> 所以,n=5。
例1.5:已知近似数x 的相对误差限为0.3%,问x 至少有几个有效数字? 解:设x 有n 位有效数字,其第一位有效数字按最不利情况取为9,则
1131111
0.3%10101010002(91)2102210n n n n
δ---===?=?=?=+??由上可得
6101000n ?=,
n ≈2.2,
所以取n=2。 指出:
也可以按首位为1,9分别计算,取较小者。
4、计算方法的余项计算
各种计算方法的余项的计算根据相应的余项定理进行。 (二)误差分析
精度水平的分析主要依据两个结论: 相对误差越小,近似数的精确度越高。
一个近似数的有效数字越多,它的相对误差越小,也就越精确。反之亦然。
例1.6: 测量一个长度a 为400米,其绝对误差不超过0.5米,测量另一长度b 为20米,其绝对误差不超过0.05米。问,哪一个测量的更精确些?
解:
0.50.125%
400
0.05
0.25%
20
a a b
b a b εδεδ======
显然,δa < δb 所以测值a 更准确一些。 答:测值a 更准确一些。 指出:
衡量测量工作的好坏用相对误差。 解决这样的题目就是三个步骤: 第一,求出两个相对误差。 第二,比较两个相对误差的大小。 第三,结论。
(三)算法分析 1、稳定性分析
算法的稳定性通过对计算的误差的扩缩情况进行分析。
例1.7:设近似值T 0=S 0=35.70具有四位有效数字,计算中无舍入误差,试分析分别用递推式
15142.8i i T T +=-和11
142.85
i i S S +=-
计算T 20和S 20所得结果是否可靠。
解:设计算T i 的绝对误差为e(T i )=T i *-T i ,其中计算T 0的误差为ε,那么计算T 20的误差为
e(T 20)=T 20*-T 20=(5T 19*-142.8)-(5T 19-142.8)=5(T 19*-T 19) =5e(T 19)=52e(T 18)=……=520e(T 0) 显然误差被放大,结果不可靠。
同理,20
2001()()5e S e S ??
= ???,误差缩小,结果可靠。
指出:
注意理论分析,因此初始近似值本身是不必要的。 2、收敛性分析
算法的收敛性分析主要是迭代法解方程的收敛性分析和迭代法解方程组的收敛性分析,其他计算方法的收敛性分析一般在具体计算过程中体现。
(1)迭代法收敛性判定的基本结论是:
定理(迭代法基本定理):对于任意的f ∈R n ,和任意的初始向量x (0)∈R n ,迭代法
x (k+1)=Bx (k)+f(k=0,1,2,…)
收敛的充分必要条件是迭代矩阵B 的谱半径ρ(B)<1。
推论:若1B <,则迭代格式x (k+1)=Bx (k)+f(k=0,1,2,…)收敛。 (2)判定雅可比迭代法、高斯—赛德尔迭代法收敛的基本依据是: 定理: 设线性方程组Ax=b ,其系数矩阵为
1112121
22
212(0)n n ii n n nn a a a a a a A a a
a a ?? ? ?=≠ ? ? ???
则雅可比迭代法迭代矩阵的特征值满足如下条件:
11121212221
2
0n n n n nn a a a a a a a a a λλλ
=;
高斯-赛德尔迭代法迭代矩阵的特征值满足如下条件:
1112
12122212
0n n
n n nn
a a a a a a a a a λλλλλλ=。
(3)系数矩阵为严格对角占优矩阵的方程组的迭代法收敛性:
定理:系数矩阵为严格对角占优的线性方程组,它的雅可比迭代和高斯-赛德尔迭代都是收敛的。
指出:
迭代法基本定理是一般结论,对任意迭代法的收敛性都能分析。限定雅可比迭代法和高斯-赛德尔迭代法则不必应用基本定理,以回避求迭代矩阵。
例1.8:已知线性方程组
1231231
232211221
x x x x x x x x x +-=??
++=??++=? 求解这个方程组的雅可比迭代法和高斯—赛德尔迭代法是否收敛?
解:
122111221A -?? ?= ? ???
,
令
2
2
1102
2
λ
λ
λ
-=,
则312300λλλλ=?===, 所以ρ(B J )=0<1
所以雅可比迭代法收敛。 而
212322
10(2)00,222λλλλλλλλλλλ
-=?-=?===, 所以ρ(B G-S )=2>1
所以高斯—赛德尔迭代法发散。
二、基本计算与基本算法 (一)秦九韶算法
秦九韶算法是一种求多项式的值的计算方法。 对任意给定的x ,计算代数多项式
1110()n n n n n P x a x a x a x a --=++
++
的值,可以利用下面的方法计算:
1210()((()))n n n n P x a x a x a x a x a --=+++++
这种算法就是著名的秦九韶算法。是我国宋朝伟大的数学家秦九韶的伟大发现。
秦九韶算法可以写成递推的形式:
10(1,
2,1,0)()n n k k k n
s a s xs a k n p x s
+=??
=+=-??=?
具体计算式,递推格式是采用如下表格形式进行计算: