牛顿法在m重根情况下的线性收敛性通用证明方法问询
嘿,你这个思路完全找对方向了!把有m重根的函数写成 ( f(x) = (x-r)^m g(x) ) 的形式就是证明的核心,我来一步步给你把推导过程理清楚:
第一步:明确m重根的函数形式
既然r是f(x)的m重根,根据定义我们有 ( f(r) = f'(r) = \dots = f^{(m-1)}(r) = 0 ),且 ( f^{(m)}(r) \neq 0 )。这时候我们可以把f(x)分解为 ( f(x) = (x-r)^m g(x) ),这里的g(x)满足两个关键条件:一是 ( g(r) \neq 0 )(否则根的重数就超过m了),二是g(x)和f(x)一样是(m+1)阶连续可导的(因为f(x)是)。第二步:计算f(x)的导数
用乘积法则对f(x)求导:
[
f'(x) = m(x-r)^{m-1}g(x) + (x-r)^m g'(x) = (x-r)^{m-1}\left[ m g(x) + (x-r)g'(x) \right]
]第三步:代入牛顿迭代公式
牛顿法的迭代公式是 ( x_{i+1} = x_i - \frac{f(x_i)}{f'(x_i)} ),把刚才得到的f(x_i)和f'(x_i)代入进去:
( f(x_i) = (x_i - r)^m g(x_i) ),( f'(x_i) = (x_i - r)^{m-1}\left[ m g(x_i) + (x_i - r)g'(x_i) \right] )
两者相除时,( (x_i - r)^{m-1} ) 可以直接约掉(因为x_i靠近r但不等于r,所以这个项不为0),得到:
[
\frac{f(x_i)}{f'(x_i)} = \frac{(x_i - r)g(x_i)}{m g(x_i) + (x_i - r)g'(x_i)}
]第四步:用误差项简化表达式
定义误差 ( e_i = x_i - r ),那么 ( x_i = r + e_i ),代入迭代公式得到 ( e_{i+1} = x_{i+1} - r ),整理后:
[
e_{i+1} = e_i - \frac{e_i g(r + e_i)}{m g(r + e_i) + e_i g'(r + e_i)}
]
给右边通分并提取公因子 ( e_i ),可以得到:
[
e_{i+1} = e_i \cdot \frac{(m-1)g(r + e_i) + e_i g'(r + e_i)}{m g(r + e_i) + e_i g'(r + e_i)}
]第五步:取极限分析收敛性
因为牛顿法是局部收敛到r的,当 ( i \to \infty ) 时,( e_i \to 0 ),所以 ( r + e_i \to r )。根据g(x)的连续性,( g(r + e_i) \to g(r) ),( g'(r + e_i) \to g'(r) )。把这些极限代入上面的式子:
分子的极限是 ( (m-1)g(r) + 0 = (m-1)g(r) ),分母的极限是 ( m g(r) + 0 = m g(r) )。
由于 ( g(r) \neq 0 ),可以直接约掉,最终得到:
[
\lim_{i \to \infty} \frac{e_{i+1}}{e_i} = \frac{m-1}{m}
]
这样就完成了任意m重根函数的收敛性证明,你看,整个过程就是基于你想到的那个因式分解展开的,是不是很顺畅?
备注:内容来源于stack exchange,提问作者Erik Eriksson

