,
0
0
fn
(x)
0
则方程(6.1.1)可改为
(6.1.2)
F(x) 0
当n 1时，方程组 (6.1.1)是一个非线性方程式
f (x) 0
(6.1.3)
定义6.1 设x * 使f (x) 0，则称x *为方程 (6.1.3)的根或零点。g(x)
1F (x(k ) k)) F(
) x
(k
1)
)
F
(
x
(
k
)
)
Bk
1
Bk
Bk
(k 0,1, 2, )
用x 表示迭代近似值，用 x表示新的近似值
迭代方程为： x＝x B1F(x) 拟Newton方程为: B(x x) F(x ) F(x) 令 s x x, y F(x) F(x) 则拟Newton方程为 Bs y
2) Newton下山法
为使 f (xk ) 具有单调下降性质：| f (xk1) || f (xk ) |
将Newton迭代法与下山法结合起来应用：
xk 1
xk
f (xk ) f (xk )
xk1 xk1 (1 )xk
即：Newton下山法计算公式为：
xk 1
xk
f (xk ) f (xk )
f1 f1
( (
x) x)
3 x1
cos(x2 x3
)
1 2
0
x12 81(x2 0.1)2 sin
x3
1.06
0
f3(x)
e x1x2
20x3
10
3
3
0
6.3 拟Newton法
拟Newton法----Newton法的修正过程
x(k 1)＝x(k ) -Bk Bk 1(x(k 1) x(
牛顿法
牛顿法的迭代公式
f(x) y
x
x*
x2 x1
x0
原理：将非线性方程线性化 —— Taylor 展开
取 x0 x*，将 f (x)在 x0 做一阶Taylor展开:
f (x)
f ( x0 )
f ( x0 )(x x0 )
f
(
2!
)
(
x
x0
)，2
在
x0
和
x
之间。
将 (x* x0)2 看成高阶小量，则有：
反复执行步骤2、3，便可得到一系列有根区间： (a, b), (a1, b1), …, (ak, bk), …
4、当 bk1 ak1 时
5、则
xk 1
1 2
(ak
bk
)
即为根的近似
①简单; ② 对f (x) 要求不高(只要连续即可) .
①无法求复根及重根 ② 收敛慢
区间套定理：若一系列闭区间an,bn
k 1, 2,
例: 用割线法求方程 f (x) xex 1 0 于[0.5,0.6]内的根
例：应用Newton法于方程
f (x) xn a 0 和f (x) 1 a 0 (a 0)分别 xn
导出n
a
的迭代公式，求
lim
k
ek 1 ek 2
, ek
n
a
xk
6.2 解非线性方程组的Newton迭代法
y 0 g '(x) 1 p1 p0
y=x y=g(x)
✓
x
x0
x1 x*
y
y=x
y=g(x)
g '(x) 1
p0
p1
x x1 x0 x*
y p0 1 g '(x) 0
y=x
✓
y=g(x) p1
x0
x*
y
y=g(x) p0
x x1
y=x
g '(x) 1
p1
x x0 x* x1
收敛性条件
x (k) i
)的二次以上的项，有
f1(x)
f1( x(k )
)
( x1
x (k) 1
)
f1 ( x(k ) x1
)
( xn
x (k) n
)
f1( x(k ) ) xn
0
,
f
2
(
x
)
f2
(
x(k
)
)
(
x1
x1( k
)
)
f2
( x(k x1
)
)
( xn
x (k) n
)
f2 ( x(k ) ) xn
1 2
(
xk3
3), k
0,1, 2,
求得近似解为： x0=1.900 x1=1.930 x2=2.095 x3=3.098 x4=13.37 得到的近似解是不收敛的，越来越发散。
由此可见迭代函数 g(x) 选取的适当，近似解将会 收敛；选取的不适当，近似解将会发散。
那么选择怎样的g(x)迭代格式才会收敛呢？
x1=1.8945647
x2=1.89352114 …………… x8=1.89328920
x9=1.89328920
由于 x8、x9 相当接近，故可取 x*≈x8=1.89328920。
如果将原方程 x3-2x-3=0 改为得 x 1 ( x3 3)
2
仍取初值 x0=1.9 , 得迭代格式如下：
xk 1
0 f ( x*) f ( x0 ) f ( x0 )( x * x0 )
xk1
xk
f ( xk ) f ( xk )
x*
x0
f ( x0 ) f ( x0 )
(2.3)
只要 f C1，每一步迭代都有 f ’( xk ) 0，而且
lim
k
xk
x *，则
x*就是
f
(x)的根。
2. Newton迭代法的收敛性与收敛阶
1
1
f1 ( x) f２( x)
2 f1(x) n f1(x)
2
f２(
x)
n
f２(
x)
x (k 1) 1
x (k 1) 2
x (k) 1
x (k) 2
0
fn (x)
1 fn (x)
2
fn (x)
n
fn (x) xx(k)
x (k 1) n
x (k) n
矩阵形式为
0
f
n
(
x)
fn
( x(k )
)
( x1
x (k) 1
)
fn
( x(k x1
)
)
( xn
x (k) n
)
fn ( x(k ) xn
)
0
记x *的第k 1次近似值为x(k1) ( x1(k1) , x2(k1) ,
以及
j
fi
(x)
fi (x) x j
, xn(k 1) )T
f1(x(k ) ) f2(x(k) )
点x*，如果迭代误差 k x*满 x足k 渐进关系式
lim
k
k 1 k p
lim k
x* xk1 x* xk p
C 0
则称此迭代格式是p阶收敛的，或称该方法具有p阶敛速。
当p = 1时，称迭代格式为线性（一次）收敛；
当p >1时，称迭代格式为超线性收敛。
当p = 2时，称迭代格式为平方（二次）收敛；
且0 g(x*) ,则称x *为方程 (6.1.3)m重根。当m 1时，x * 单根，
x *满足条件 f (x*) 0，
f '(x*) 0
6.1.1 二分法
定理1：设函数 f (x) 在区间[a, b]上连续，如果f (a) f (b) < 0， 则方程 f (x) = 0 在[a, b]内至少有一实根x*。
F ( x(k ) ) F '( x(k ) )( x(k1) x(k ) ) 0
1 f1(x) 2 f1(x) n f1(x)
其中
F '( x(k ) )
1
f２(
x)
2
f２(
x)
n
f２(
x)
1 fn (x)
2
fn
(
x)
n
fn
(x)
x x( k
)
为F ( x)于点x( k )处的Jacobi矩阵
| (x1)- (x2)|≤L|x1-x2|
例：求方程f (x) xex 1 0 (或x ex )在[1，ln 2]中的解。 2
若要求 x * xk 106 ，迭代次数k至少应为多少？
6.1.3 Newton迭代法
1．迭代过程的收敛速度
定义6.3：设迭代格式xk+1= (xk)收敛于方程x= (x)的不动
定理6.1：假定x = (x) 满足下列条件
(1) (x)在[a, b]上连续
(2) 当x[a, b]时， (x) [a, b]
(3) '(x) 存在，且对任x[a, b]，有
'(x) L 1
(6.1.10)
则方程x = (x)在[a, b]上有唯一的根x*, 且对任意初
值 x0[a, b]时，迭代序列xk+1= (xk) (k = 0, 1, …)所定
第6章 非线性方程(组)和最优化问题 的计算方法
6.1 方程式求根(二分法、迭代法和 Newton迭代法)
给定如下非线性方程组
fi (x1, x2 , , xn ) 0, i 1, 2, , n 引入向量，向量函数记号
(6.1.1)
x1
x
x2
,