计算方法第六章(迭代法)分析
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
得到
x0 [a, b] , xk 1 ( xk ) ( k 0,1, 2, ) lim xk x , ( x ) x
* * k *
例：在区间[1，）上求解方程
x ( x ) , ( x) x 1 1 1 ( x) 1 2 x 2
( x) ( x) (1 ) x
1 ( ) ( ) (1 ) 0 1 ( )
xk 1 xk 2 xk 1 ( xk 1 ) xk 2 ( xk 2 )
例子：求方程
x 2x 5 0
收敛阶： 定义：设 lim xk
k
, k xk ( xk ) ，如果存在实数 p
和非零常数 c，使：
k 1 k
p
c
则称序列 p 阶收敛，特别，p＝1时，称为线性收敛，p > 1 时，
称为超线性收敛，p＝2时称为平方收敛。 p 越大，序列收敛越快。如果是线性收敛，则 0 < c < 1
设迭代 xk 1 ( xk )收敛到a，充分可导，则迭代p 阶收敛
(a ) (a )
( p 1) (a ) 0, 但 ( p ) (a ) 0
加速收敛技术： 1、松弛法 选择适当的常数 （松弛因子），令
* xk xk (1 ) xk 1 ( xk 1 ) (1 ) xk 1
第六章
迭代法
第一节 非线性方程求根 （ f ( x) 0 ） 1、二分法 利用连续函 数的性质进 行对分。 计算框图为：
压缩映射： 集合 A 上的映射 : A A ， A 上两个点 x1 , x2 之间的距离 记为 d ( x1 , x2 ) ，如映射满足下面条件，称为压缩映射
x1, x2 A, k (0,1) , s.t. d ( ( x1 ), ( x2 )) kd ( x1, x2 ) 例：设函数 f ( x) 满足： f ( x) k 1 ，则该函数为压缩映射
Aitken加速法（适用于线性收敛情况）
xk xk 1 xk 1 xk 2
( xk )( xk 2 ) ( xk 1 ) 2
2
( xk xk 1 ) xk xk 2 xk 1 xk 2 ( xk xk 1 ) x xk xk 2 xk 1 xk 2
定理： ( x) 满足下面条件时，为压缩映射： （1）当 x [a, b] 时， ( x) [a, b] （2）存在正数 L < 1，使得
x ( x) ，可以通过迭代求解。
x [a, b] , ( x) L 1 则方程在区间上有唯一解 x* ，且解可以用下面迭代
3
的根
迭代格式：
xk 1 3 2 xk 1 5
取 ＝1.15，
x xk wk.baidu.com(1 ) xk 1
* k
计算结果要求准确到小数后8位数字
2.154434739312699 2.103612082648350 2.095927464327627 2.094760599916342 2.094583304649520 2.094556363492997 2.094552269550262 2.094551647438705 2.094551552903205 2.094551538537676 2.094551536354704
2.102599958448522 2.094749937881704 2.094556446501749 2.094551657513653 2.094551538972266 2.094551536038016
x=2.5 10 y=x x=(2*y+5)**(1.0/3.0) if (abs(x-y).lt.0.00000001) then goto 15 endif 15 x=2.5 goto 10 20 y=x x=(2*y+5)**(1.0/3.0) x=1.15*x+(1.0-1.15)*y if (abs(x-y).lt.0.00000001) then goto 30 endif goto 20 30 end
* k
( xk 1 xk )2 x xk 1 xk 1 2 xk xk 1
2
( xk 1 xk ) x xk 1 xk 1 2 xk xk 1
2 * k
3、插值加速法
由线性插值公式：
x xk x xk 1 y xk xk 1 xk 1 xk xk xk 1
x xk x xk 1 x xk xk 1 xk 1 xk xk xk 1
xk 1 xk 1 xk 2 x xk 1 2 xk xk 1
可用迭代法求解，迭代序列 x0 100, xn1 ( xn )
xn
误差估计：第k步迭代计算值与精确值误差为
Lk a xk x1 x0 1 L
使用迭代法求解方程值得注意的事项：
1、将要求解的方程化成 x ( x) 的形式。
2、该迭代法第一个条件不易验证。因此，实际使用时，总在 根的附近区间内进行迭代计算，以保证每次迭代的值都在 迭代区间内。 3、L很小时迭代收敛非常快，但如果L与1很接近，则收敛相 当慢。
定理：如果 为闭集合 A 上的压缩映射，则方程 x = (x) 在集 合 A 上有唯一解。且解可以用下面迭代得到：
x0 A , xn 1 ( xn ) xn x , x ( x )
* * *
k xn x x1 x0 1 k
*
n
2、简单迭代：
对于形如的方程
x0 [a, b] , xk 1 ( xk ) ( k 0,1, 2, ) lim xk x , ( x ) x
* * k *
例：在区间[1，）上求解方程
x ( x ) , ( x) x 1 1 1 ( x) 1 2 x 2
( x) ( x) (1 ) x
1 ( ) ( ) (1 ) 0 1 ( )
xk 1 xk 2 xk 1 ( xk 1 ) xk 2 ( xk 2 )
例子：求方程
x 2x 5 0
收敛阶： 定义：设 lim xk
k
, k xk ( xk ) ，如果存在实数 p
和非零常数 c，使：
k 1 k
p
c
则称序列 p 阶收敛，特别，p＝1时，称为线性收敛，p > 1 时，
称为超线性收敛，p＝2时称为平方收敛。 p 越大，序列收敛越快。如果是线性收敛，则 0 < c < 1
设迭代 xk 1 ( xk )收敛到a，充分可导，则迭代p 阶收敛
(a ) (a )
( p 1) (a ) 0, 但 ( p ) (a ) 0
加速收敛技术： 1、松弛法 选择适当的常数 （松弛因子），令
* xk xk (1 ) xk 1 ( xk 1 ) (1 ) xk 1
第六章
迭代法
第一节 非线性方程求根 （ f ( x) 0 ） 1、二分法 利用连续函 数的性质进 行对分。 计算框图为：
压缩映射： 集合 A 上的映射 : A A ， A 上两个点 x1 , x2 之间的距离 记为 d ( x1 , x2 ) ，如映射满足下面条件，称为压缩映射
x1, x2 A, k (0,1) , s.t. d ( ( x1 ), ( x2 )) kd ( x1, x2 ) 例：设函数 f ( x) 满足： f ( x) k 1 ，则该函数为压缩映射
Aitken加速法（适用于线性收敛情况）
xk xk 1 xk 1 xk 2
( xk )( xk 2 ) ( xk 1 ) 2
2
( xk xk 1 ) xk xk 2 xk 1 xk 2 ( xk xk 1 ) x xk xk 2 xk 1 xk 2
定理： ( x) 满足下面条件时，为压缩映射： （1）当 x [a, b] 时， ( x) [a, b] （2）存在正数 L < 1，使得
x ( x) ，可以通过迭代求解。
x [a, b] , ( x) L 1 则方程在区间上有唯一解 x* ，且解可以用下面迭代
3
的根
迭代格式：
xk 1 3 2 xk 1 5
取 ＝1.15，
x xk wk.baidu.com(1 ) xk 1
* k
计算结果要求准确到小数后8位数字
2.154434739312699 2.103612082648350 2.095927464327627 2.094760599916342 2.094583304649520 2.094556363492997 2.094552269550262 2.094551647438705 2.094551552903205 2.094551538537676 2.094551536354704
2.102599958448522 2.094749937881704 2.094556446501749 2.094551657513653 2.094551538972266 2.094551536038016
x=2.5 10 y=x x=(2*y+5)**(1.0/3.0) if (abs(x-y).lt.0.00000001) then goto 15 endif 15 x=2.5 goto 10 20 y=x x=(2*y+5)**(1.0/3.0) x=1.15*x+(1.0-1.15)*y if (abs(x-y).lt.0.00000001) then goto 30 endif goto 20 30 end
* k
( xk 1 xk )2 x xk 1 xk 1 2 xk xk 1
2
( xk 1 xk ) x xk 1 xk 1 2 xk xk 1
2 * k
3、插值加速法
由线性插值公式：
x xk x xk 1 y xk xk 1 xk 1 xk xk xk 1
x xk x xk 1 x xk xk 1 xk 1 xk xk xk 1
xk 1 xk 1 xk 2 x xk 1 2 xk xk 1
可用迭代法求解，迭代序列 x0 100, xn1 ( xn )
xn
误差估计：第k步迭代计算值与精确值误差为
Lk a xk x1 x0 1 L
使用迭代法求解方程值得注意的事项：
1、将要求解的方程化成 x ( x) 的形式。
2、该迭代法第一个条件不易验证。因此，实际使用时，总在 根的附近区间内进行迭代计算，以保证每次迭代的值都在 迭代区间内。 3、L很小时迭代收敛非常快，但如果L与1很接近，则收敛相 当慢。
定理：如果 为闭集合 A 上的压缩映射，则方程 x = (x) 在集 合 A 上有唯一解。且解可以用下面迭代得到：
x0 A , xn 1 ( xn ) xn x , x ( x )
* * *
k xn x x1 x0 1 k
*
n
2、简单迭代：
对于形如的方程