第十章非线性方程及非线性方程组解法
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
可先用二分法或经验确定迭代初值x0 0.5,再按牛
顿公式进行迭代。
Newton法具有收敛快,稳定性好,精度高等优点,是求 解非线性方程的有效方法之一。但它每次迭代均需计算函 数值与导数值,故计算量较大。而且当导数值提供有困难 时, Newton法无法进行。
牛顿法应用举例
对于给定正数C,应用牛顿法解二次方程 x2 C 0
f(x) f( x0 ) (x- x0 ) f ' ( x0 )
xx0 2
f ('' x0 )
2!
取线性部分近似代替
f(x) 0 : f(x ) f ' (x )xx 0
0
0
0
解出x作为近似根x : 1 x1 x0 f( x0 ) f (' x0 ) ( f (' x0 ) 0)
xk
C 2
q2k C 1 q2k .
对任意x0 0,总有 q 1,故当k 时,xk C .
(b
a)
n
n,
x 2 x 误差: x
n
ba
n1
n
2
计算简便,容易估计误差,但收敛较慢。
§2. 迭代法
改写方程:f (x) 0 x (x) 且 连续。
x x x 建立迭代格式: ( ),得到序列 { }
n 1
n
n
则
若
{x
}
n
收敛必收敛到
f (x) 0
的根:
limxn1 lim
n
n
证:由于 ' (x*) 0 1,故xk1 (xk )具有局部收敛性。
将 (xk )在根x*处展开,由条件有
(xk
)
( x* )
( p) ( )
p!
( xk
x*) p
xk 1
x*
wk.baidu.com p) ( )
p!
( xk
x*) p
ek1 ( p) (x* )
ekp
p!
迭代过程的收敛速度依赖于迭代函数 (x)的选
L xy
即 ( x)满足上一定理的条件(2),故结论成立。
误差估计式
x* xn
Ln 1 L
x1 x0 表明,常数L越小,
简单迭代法收敛越快。因而构造迭代函数 (x)的原则
是使 ' (x) 在有根区间[a, b]上有尽可能小的上界。
对任意正整数p有
x x x x x x x x
L
n p
于是,对任意正整数n, p,有
x x x x x x x x
L
n p
n
n p
n p1
n p1
n p2
n1
n
L x x L x x L x x n p1 n p2 L n
1
0
1
0
1
0
L L L x x n( p1 p2 L 1)
1
0
1 Lp
1 L
Ln
x1
x0
令p ,得
x x x x* Ln
取。如果当 x [a, b]时 ' (x) 0,则该迭代过程只
可能是线性收敛。
对牛顿公式
xk 1 xk
f (xk ) f ' (xk )
其迭代函数为
(x)
x
f (x) f '(x)
由于
'(x)
f (x) f '' (x)
f
'
2
( x)
假定x*是f (x)的一个单根,即f (x* ) 0, f ' (x*) 0,则
(
x
)
n
lim
n
x
n
若
{x
}
n
收敛,即
lim xn x*,则:
n
x* (x*) f (x*) 0
迭代过程的几何表示
x (x) :
y x 交点即真根。
y (x)
yx
y
Q1
Q2
P* P2
O x* x2
P1
x1
y (x)
P0
x0
x
例:求方程 f (x) x3 x 1 0 在x0 1.5附近的根x*. 解:(1) 将方程改写为 x 3 x 1
实际用迭代法计算时,先用对分区间法求较好的初值, 然后再进行迭代。
迭代法加速(埃特金方法)(1)
设x0是根x*的某个预测值,通过两次迭代校正
有
x1 ( x0 )
x2 ( x1)
由微分中值定理,有
x1 x* (x0 ) (x* ) ' (1)(x0 x* ) x2 x* (x1) (x* ) ' (2 )(x1 x* ) 假定 ' (x)改变不大,近似地取某个近似值L,
n 1 L 1 0
收敛充分性定理(二、1)
定理. 设函数 (x)在区间[a, b]上满足条件 (1)对任意x [a, b],都有a (x) b;
(2)存在常数0 L 1, 使得对一切x [a, b],都有
'(x) L
则方程x (x)在[a, b]内有唯一的根x*,且对任何
初值x0 [a, b], 迭代序列
依此类推,得
xn x* Ln x0 x*
因 0 L 1, 所以
lim
n
xn
x*
即对任意初值x0 [a, b], 迭代序列xn均收
敛到方程的根x*。
类似地,对任意正整数k,有
xk 1 xk (xk ) ( xk 1) L xk xk 1
L Lk x1 x0
收敛充分性定理(一、4)
依次产生迭代格式,称 Newton 法:
x x x f ( )
n1
n
n
f ' (xn),
n 0, 1, 2,
Newton 法的几何解释
当 x0 在取后(在真根 附近),过 x0, f (x0)
作 f (x) 的切线,则切线方程:
y
f (x0)
f
'
(x0)(x x0)
(0,f(0))
则称该迭代过程是p阶收敛的。 p 1为线性收敛,p 1为超线性收敛, p 2为平方收敛。
定理:对于迭代过程xk1 (xk ),如果 ( p) (x)在所
求根x*的邻近连续,并且
' (x*) "(x*) L ( p1) (x*) 0; p (x*) 0
则该迭代过程在点x*邻近是p阶收敛的。
假设方程x (x)在[a, b]内有两个根x1* x2*,
由条件(2),有
x1* x2* (x1* ) (x2* ) L x1* x2* x1* x2*
导出矛盾,唯一性得证。
收敛充分性定理(一、3)
对任意x0 [a, b],由迭代公式有
xn x* (xn1) (x* ) L xn1 x*
由上式知 ' (x* ) 0,由上述定理知,牛顿法在根x*的
邻近式平方收敛的。
例:用牛顿法解方程 xex 1 0.
解:f (x) xex 1, f '(x) ex (1 x)
牛顿法迭代函数为 (x) x
f (x) f '(x)
x
x ex 1 x
牛顿公式为
xk 1
xk
xk 1
e xk xk
改进
x%k1 ( xk )
xk1 ( x%k1)
xk 1
xk 1
(xk 1 x%k 1)2 xk 1 2x%k 1 xk
§3. Newton 法
非线性问题的最简单解法是线性近似.
将非线性方程线性化,以线性方程的解逐步逼
近非线性方程的解,这就是Newton法的基本思 想
f ( x ) 0 在真根附近 x 点展开成 Taylor 级数:
1 2.25
根据定理,任取x0 [1.5, 2],由这两种等价方程所构造的 简单迭代方法都收敛,且第一种所产生的迭代序列收敛较快。
收敛充分性定理(三、1)
定理:如果函数 (x)在x*的一邻域O(x*, *) 内连续可微,x*为方程x (x)的根,且 ' (x) 1,则存在正数 , *,使得对任意 x0 [x* , x* ], 迭代序列
即导出求开方值
C的计算程序
xk 1
1 2
( xk
C xk
).
现证此迭代公式对于初值x0 0都是收敛的。由迭代公式得
xk 1
C
1 2 xk
( xk
C )2
xk 1
C
1 2xk
( xk
C )2
2
xk1 xk1
C C
xk xk
C C
.
2k
xk xk
C C
x0 x0
C C
记 q x0 C , x0 C
则由x1 x* L(x0 x* ) x2 x* L(x1 x* )
x1 x* x2 x*
x0 x* x1 x*
x*
x2
( x2 x1)2 x0 2x1 x2
此种加速需用两次迭代值进行加工。
迭代法的加速(埃特金方法)(2)
如果将一次改进值作为一步,
则计算公式如下:
校正 再校正
0 1 ξ
与x轴交点( y 0 )为第二个近似根 x1 : x1 x0 f ( x1 ) f (' x0 )'
迭代法收敛定义
定义:设迭代过程xk1 (xk )收敛于方程 x (x)的根x*,如果迭代误差ek xk x*
当 k 时成立下列渐进关系式 ek1 C (C 0为常数) ekp
f (2) 1 0
[1.5, 2] 为有根区间。
(1)x x 1 1(x) 因 1.5 1.5 1 1(x) 2 1 2
且
1' (x) 2
1 1 1 x 1 2 2.5 3.162
(2)
x 1 1 x
2(x)
因
1.5
1
1 2
2
(
x)
1
1 1.5
2
且
' 2
(
x)
1 x2
1 1.52
收敛充分性定理(一、1)
定理. 设函数 (x)在区间[a, b]上满足条件 (1)对任意x [a, b],都有a (x) b;
(2)存在常数0 L 1, 使得对一切x, y [a, b],都有
(x) ( y) L x y
则方程x (x)在[a, b]内有唯一的根x*,且对任何
初值x0 [a, b], 迭代序列
n
n p
n p1
n p1
n p2
n1
n
(L L x x x x
p1 p2 L 1) 1
n1
n 1 L n1
n
令p ,得
x x x x * 1
n 1 L n1 n
x x 可通过检查 来判断迭代过程应否终止。
n1
n
例:用简单迭代法求方程 f (x) x2 x 1 0 的根。
解:因 f (1.5) 0.25 0,
第十章 方程求根
求解非线性方程
f (x) 0 f 是非线性函数,
例:代数方程
a x a x a x a f (x) n
n1 L
0, n 1。
n
n1
1
0
例: 超越方程
f (x) ex sin x 0
§1. 非线性方程实根的对分法(二分法)
设 f (x) 在[a,b] 上连续且 [a,b] 有且仅有一个根又
xn1 ( xn ) (n 0,1,L )
均收敛于x*,并有
x* xn
Ln 1 L
x1 x0
收敛充分性定理(二、2)
证:设x, y为[a,b]上任意两点,由微分中值
定理,在x, y之间至少存在一点,使得 (x) ( y) '( )(x y)
(x) ( y) '( )(x y)
f (a) f (b) 0。则可用对分法:
不妨设 f (a) 0, f (b) 0
1), 若
f
a
2
b
0
输出根
x
a b , 否则:若 2
f
a
2
b
0,
令
a1
ab 2
,b1
b
反之
ab b1 2 , a1 a.
2 ),对 [ a1 ,b1 ] 区间重复1)的计算,并产生 [ a2 ,b2 ],
由此建立迭代公式 xk1 3 xk 1 (k 0,1, 2L )
k0 1
2L 7
8
xk 1.5 1.35721 1.33086 L 1.32472 1.32472 迭代收敛。
(2) 若将方程改写为 x x3 1
建立迭代公式 k0
xk1 xk3 1.
1
2L
xk 1.5 2.375 12.39 L 迭代不收敛。
3),
若
f
ai
2
bi
0,则得到根
x ai bi . 2
二分法的收敛性
f (x)
二分法产生一个有根区间:
[a, b]
[a1,b1] L
[a
,
n
b
]
n
[a
,
n
b
]
n
区间长度:
a1 a x*
b1
x0
b
b a b a 1 ( ) L
n
n2
n 1
n 1
a b 当n 足够大时,取近似值
1
2n
xn1 (xn ) (n 0,1, 2,L )
收敛于x* .
收敛充分性定理(三、2)
证:因 ' (x)在O(x*, *)内连续,且 ' (x) 1,故存 在正数L 1, *, 使得对x [x* , x* ],有
'(x) L 1 另一方面,由 (x*) x*, 又有
(x) x* (x) (x*) L x x* 即 (x) [x* , x* ]。由上面定理知,迭代序列 xn1 (xn )收敛于x*。
xn1 (xn ) (n 0,1,L )
均收敛于x*,并有
x* xn
Ln 1 L
x1 x0
收敛充分性定理(一、2)
证:由条件(2)知(x)在[a, b]上连续。 令 (x) x (x),则 (x)在[a,b]上连续,且
(a) a (a) 0, (b) b (b) 0 故存在 [a,b],使得() 0,即 (), 所以方程x (x)在[a,b]内有根。