收敛加速的方法ppt课件

n
|
x0
x*
|
0
( 0<L<1 )
所以,
lim
k
xk
x* ，故迭代格式收敛
| xk x* || xk xk1 xk1 x* |
| xk xk1 | | xk1 x* |
| xk xk1 | L | xk x* |
(1 L) | xk x* || xk xk1 |
|
x*
xk
|
此区间有唯一的根 x.
[a,b]称为有根区间.
二分算法
原理:
把区间[a,b]二等分,分点为 x0
1 (a 2
b),计算函
数值 f ( x0 ),则可能出现三种情况:
1 如果 f ( x0 ) 0，则 x0是 f ( x)在区间[a,b]的零点.
则
x a b 2
3
2 如果 f ( x0 ) f (a) 0,则 f ( x)在区间[a, x0 ]内有零点;
f (x)= 0 x ( x)
(1)
若存在 x*,使得x* ( x*) ,则称x*为不动点。
在根x*的附近取一点x0作为x*的预测值,也叫迭代初 值。
13
2 按迭代格式进行计算
把x0代入(1)的右端，得
x1 ( x0 )
如果 x1 x0 ，则 x1 x* 。
如果 x1 x0 ，把 x1作为根的新的预测值代入(1)，得
x2 ( x1 ) 如果 x2 x1 ，则 x2 x* 。
如果 x2 x1 ，把 x2作为根的新的预测值代入(1)......
如此重复上述步骤，则有迭代公式
xk1 ( xk ) ( k = 0, 1, 2, ···)
14
其中，( x) ：迭代函数，得到迭代序列
x0 , x1 , , xk ,
2. 利用迭代法计算满足一定精度的根的近似值.在隔 根区间找到一个(或者多个)出发值 x0,按某种方法产 生一个序列 x0 ,..., xn ,...此序列在某种条件下收敛于 方程的根 x.
2
2.1 二分法
设函数 f ( x)在区间[a,b]上连续,且 f (a) f (b) 0,根据 连续函数的性质知 f ( x) 0在[a,b]内有根.不妨设它在
4
1
对任意 k 0,设第 k 个区间为[ak,bk],取中点 xk 2 (ak bk ) , 计算 f(xk).
若 f(xk)=0，则 xk 就是方程的根; 若 f(xk) f(ak)<0, 取 ak+1=ak, bk+1= xk; 若 f(xk) f(bk)<0, 取 ak+1= xk, bk+1= bk.
k 0,1,2,
x1 2.375, x2 12.39
当 k增大时，xk随之增大而不趋于任何极限，此时迭 代过程发散。
通过此例说明，迭代过程只有在一定条件下才可能
收敛。一个发散的过程没有任何意义。
20
定理2.3 如果 (x) C1[a, b] ,满足条件:
(1) a (x) b ; (2) | (x) | L 1
21
唯一性：
设在[a, b] 上存在两个根x1*和x2*，则 ( x1*) x1 * ( x1*) x1 *
由微分中值定理
| x1* x2* || ( x1* ) ( x2* ) |
| ( ) || x1* x2* |
L 1 ，必有
L | x1* x2* |
|
x1*
x
▪将一个计算过程反复进行
▪一种常见常用的计算技术
▪构造有效的迭代格式
▪选取合适的迭代初值
▪对迭代格式进行收敛性分析
12
迭代原理 迭代是一种逐次逼近过程求解问题的方法. 已知方程 f ( x) 0的一个近似根后，通常使用某个固定公式反 复校正根的近似值
作法：迭代的三个主要部分
1 选取初值
把给定的方程 f ( x) 0 改写成等价形式
1 1
L
|
xk 1
xk
|
24
估计迭代次数 满足精度 要求所需的迭代次数
只需要 1
1 L
xk1
xk
成立即可.
注意 当L越小时,序列{ xk }收敛越快.
在实际算法设计中,我们常用条件 xk1 xk 作 为迭代结束;也可以用相对误差 xk1 xk 作为
xk1
迭代终止条件。
为避免出现死循环,可以设置一个最大的迭代次数.
Th2.4 给的是一个全局收敛的迭代法.条件不易检验， L不 易得,实际应用时通常只在不动点 x*的邻近考察其收敛性。
25
三 局部收敛性和收敛阶
定义 2.2 如果存在( x)的不动点 x的某个邻域
U( x, ) [x , x ],
对任意的 x0 U ( x , ), 迭代过程 xk1 ( xk )产生的序 列{ xk }均收敛于 x,则称迭代公式是局部收敛的. 定理 2.5（局部收敛性的充分性条件）设( x)在 x ( x)
的根 x邻近有连续的一阶导数,且
'(x) 1
则迭代过程局部收敛.
证明 因'( x)连续,所以存在 x的一个邻域 U( x, ) [x , x ],
当 x U( x, )时,有 '( x) L 1,且
p1
o x* x2
x1
x0
x
如果 p1, p2 , , pk , 逐渐逼近p*，－－－迭代过程收敛
16
y (x)
y
y=x
p*
p0
p0
p1
p1
p2
o
x
x2
x1 x0 x*
如果 p1, p2 , , pk , 逐渐远离p*，－－－迭代过程发散
（无意义）
17
例2.2 求方程 f (x)=x3 – x – 1 = 0
在x =1.5附近的根 x*。
解 设将方程改写成下列形式
由此得迭代公式
x 3 x1
xk1 3 xk 1, k 0,1,2,
迭代初值取x0 =1.5，计算值用6位数字表示。
迭代结果如下表
18
k
xk
k
xk
0
1.5
5
1.32476
1
1.35721
6
1.32473
2
1.33086
7
1.32472
3
3 判别收敛 如果迭代序列的极限存在，则迭代过程收敛，显然有
存在，则称迭代过程发散。
上述迭代过程也称不动点迭代法。
15
几何意义
方程 x ( x) 求根，在几何上就是确定曲线 y ( x)
与直线 y x 的交点 p*
y
y=x
y (x)
p0
p0
p1 p* p2
第四步: 若 b a 1,则转第二步;否则输出 x0,结束.
作业：1、用二分法求方程 f ( x) x3 4x2 10 0
在区间[1,2］内的一个实根，要求误差不超过0.005。
11
2.2 迭代法
111
一种圆周率的计算方案: 初值: x0=1
4
(1)n
1
3
5
7
迭代格式: xn xn1 2n 1 ( n=1,2,3,······)
x0 , x1 , , xk ,
由于
6
定理 2.2 设 x 为方程 f ( x) 0 在区间[a,b]内的唯一根,
f ( x)满足 f (a) f (b) 0,则二分法计算过程中第 k 个区间
[a,b]的中点 xk 满足不等式
x
xk
1 2 (bk
ak )
1 2k1
(b
a
).
证明 因为ak x bk ,所以
1.32588
8
1.32472
4
1.32494
从表中可看到 x7与 x8完全相同，这时可认为x8已 满足方程， x8 即为所求根的近似值。
x* x8 1.32472
上述迭代过程是收敛的。
19
如果将方程改写成下列形式
据此有迭代公式
x x3 1
xk1
x
3 k
1,
迭代初值仍取 x0 =1.5，则有
[1, 1.5］内的一个实根，要求误差不超过0.005。 解 由公式估计所要 二分的次数
k ln(1.5 1) ln 0.01 5.644 ln 2
即只要二分6次，便能达到所要求的精度。
9
计算结果 f (1) 1 0, f (1.5) 0
k
ak
bk
xk
f(xk)
0
1.0
1.5
1.25
x xk
1 2
(ak
bk
)
x
=
1 2
(ak
x )
1 2
(bk
x
)
1 2 ( ak
x
bk
x
)
1 2
(( x
ak
)
bk
x )),
7
因此
x xk
1 2 (bk
ak )
ba 2k1 .
因此,二分法产生的序列收敛.
利用二分法收敛定理,对于事先给定的精度要求,可以估 计出二分法结果满足要求的迭代次数.
* 2
|
0 即x1*
x2*
x*
22
定理2.4 如果 (x) C1[a, b] ,满足条件:
(1) a ( x) b ; (2) | ( x) | L 1
则对任意的 x0∈ [a, b] , 迭代格式 xk1 ( xk )
产生的序列 { xk }收敛到不动点x*,且有事后误
差估计式
|
x*