- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
算法的特点,描述
第1章 绪论
例.x 2.718281828L , x1 2.71828325,则x1的有 6 位有效位数 若fl(x) 2.71828225,则有 7 位有效位数
例.在F(10,5, -2,3)中有多少个数?
3
例.下列各式均与
33+
8 8
等价,在浮点数系F(10,5,-10,10)中
F
)
I
D1
A
Gauss-Seidel:G d
(D (D
E)-1 F E)1b
收敛性判定定理
TH2.6 G 1,则迭代格式收敛
TH2.7 A为严格对角占优, Jacobi格式收敛 TH2.8 A为严格对角占优,Gauss - Seidel格式收敛 TH2.9 A对称正定,Gauss - Seidel收敛; 2D - A对称正定, Jacobi收敛 TH2.10 迭代格式x(k1) Gx(k) d收敛的充要条件为(G) 1 TH2.11 迭代格式的误差估计
哪个公式能获得最准确的结果:(17-6 8)3, 1 , (17+6 8)3
6
1
1
3
8
, 3
6 ,19601 6930 8
8, 19601 6930
8
第1章 绪论
例.证明在浮点数系F ( ,t, L,U )中, 浮点数的相对误差
(x) x - fl(x) 满足 (x) 1 1-t
GG分解, 追赶法
迭代解法
Jacobi迭代法 Gauss-Seidel迭代法
Gauss消去法 :消去的时间复杂度o(n3),回代 : o(n2 )
列主元Gauss消去法 :消去的时间复杂度o(n3),回代 : o(n2 )
LU分解:L : 单位下三角阵,U : 上三角阵,时间复杂度o(n3)
LDU分解:L : 单位下三角阵, D : 对角阵,U : 单位上三角阵,时间复杂度o(
x* x%
x% A
Cond(A)
x*
x* A
(3)当系数矩阵有扰动A, 右端向量有扰动b
x x%
k
x
A 1 k
b b
A A
A
其中k cond (A) A1 A
第2章 线性代数方程组
迭代算法:构造x(k1) Gx(k) d,判断收敛
Jacobi:G d
D -1 ( E D1b
第2章 线性代数方程组
1
-
1
1
2
例:矩阵A= 1
3
01
,则A1 _______, A ___, A ___
1
1 4
0 0 1
P36, P37
2 1
例:
若矩阵A12a来自可以分解为GGT的形式,
其中G为下三角阵,
a 1
且对角元均为正,问a的取值范围,并请按此要求将此a分解
《计算方法》总结
目录
第1章 绪论 第2章 线性代数方程组 第3章 数据近似 第4章 数值微积分 第5章 非线性方程求解 第6章 常微分方程数值解法 第7章 最优化方法简介
(误差分析基础)
(基本工具)
(计算方法应用)
第1章 绪论
1.误差:近似值与真正值之差
分为模型误差、数据误差、截断误差、舍入误差
GG分解:针对对称正定矩阵,o(n3 / 6),加n个开方运算
带状矩阵分解:三对角阵分解,追赶法
第2章 线性代数方程组
范数: 定义,性质.向量与矩阵范数的相容性,等价性
方程组的条件数: m cond(A) A A-1
(1)当右端向量有扰动b
x x* A A1
b
x*
b
(2)当系数矩阵有扰动A
能少,应该式如何计算:_______
例.在浮点数系下,计算x2 16x 1 0的两个根,应如何 计算才能使精度较高?
例.对于函数f (x)在某个区间上连续可微,则求f (x)的近似条件数
第2章 线性代数方程组
Gauss解法
线性方程组解法
数值解法
列主元Gauss解法 矩阵分解法:LU分解, LDU分解,
x
2
证明: 设x
( d1
d2
2
d3
3
d t 1
t1
)l,
其中1 d1 ,0 d j ( j 2,3,...)
若d t 1
1 2
有fl( x)
(d1
d2
2
d3
3
dt
t
)
l
此时x
fl( x)
d t 1
t1
l
1 2
1
t
l
1 lt
2
由于1
d1
,有
x
1
l
x fl( x) 1 1t
良态问题:输入数据相对小的扰动不会引起解的相对大的变化
条件数:当输入数据具有 x的误差,引起问题的结果误差为 f (x) 则cond ( f ) sup f (x) x
5.方法的稳定性
数值稳定:若初始误差导致最终解的误差能被有效地控制 数值不稳定:若初始误差导致最终解的误差不能被有效地控制
6.算法 由有限个无二义性法则组成的一个计算过程
2.数制表示
(1)实数x可以表示以下形式的 进制t位有效数字
x
( d1
d2
2
L
dt
t
) l ,1
d1
,0
dj
,
j 2,3,L ,t
(2)有效数字: 指一个近似数的有意义的数字的位数
若 x 0.d1d2 L dt L 10l , x% 0.d1d2 L d%t 10l , 如果 x x% 0.510lt ,则称x%有t位有效数字
fl
(
x y
)
(1
-
3
)(
x y
)
浮点运算的注意事项
(1)避免产生大结果的运算,尤其是避免小数作为除数 参加运算;
(2)避免“大”“小”数相加减; (3)避免相近数相减,防止大量有效数字损失; (4)尽可能简化运算步骤,减少运算次数。
第1章 绪论
4.问题的性态:问题的解对原始数据扰动的敏感性
病态问题:输入数据相对小的扰动引起解的相对大的变化
x
2
第1章 绪论
同理,
若d t 1
1 2
有fl( x)
(d1
d2
2
d3
3
dt
t
1)
l
x
fl( x)
dt1 t1
l
1 2
1
t
l
1 lt
2
由于1
d1
,有
x
1
l
x fl( x) 1 1t
x
2
第1章 绪论
例.为了使计算y 10 3 4 6 的乘除法次数尽可
x -1 x -12 x -13
(3)浮点数系:表示为F( ,t, L,U ) 数的个数:2( 1) t1(U L 1) 1
上溢:l U 下溢:l L
第1章 绪论
3.舍入误差:对数进行舍入,得到有t位尾数的浮点数 fl( x)
相对舍入误差: (x) x fl(x)
x
(x) 1 1t
2
性质 : fl(x y) (1-1)(x y) fl(xy) (1-2)(xy)