非线性方程的单根、重根:
\,\,\,\,\,\,\,\,\,\, 给定非线性方程 f ( x ) = 0 f(x)=0 f(x)=0:
- 设有正整数 m 使得 f ( x ) = ( x − a ) m g ( x ) f(x)=(x-a)^mg(x) f(x)=(x−a)mg(x)且 g ( a ) ≠ 0 g(a) \neq0 g(a)=0,则当 m≥2 时,称 α 为 f ( x ) = 0 f(x)=0 f(x)=0的 m重根;当 m =1 时,称 α 为 f ( x ) = 0 f(x)=0 f(x)=0的单根。
- 若 α 为 f ( x ) = 0 f(x)=0 f(x)=0的 m重根的充要条件: f ( a ) = f ′ ( a ) = . . . = f m − 1 ( a ) = 0 , f m ( a ) ≠ 0 f(a)=f'(a)=...=f^{m-1}(a)=0,f^m(a)\neq 0 f(a)=f′(a)=...=fm−1(a)=0,fm(a)=0。
数值方法解非线性方程的根:
- 二分法
- 迭代法:
- 单点迭代法(不动点迭代, Newton 迭代法等)
- 多点迭代法(弦截法等)
(一)迭代法的基本理论:
\,\,\,\,\,\,\,\,\,\, 迭代法是一种逐次逼近的方法,它的基本思想是通过构造一个递推关系式 迭代格式 ) ,计算出根的近似值序列 并要求该序列收敛于方程的根。
- (一)单点迭代法的基本形式:(以Picard不动点迭代法为例)
\,\,\,\,\,\,\,\,\,\, 1.建立迭代公式:将方程f(x)=0改写成等价形式x=φ(x),进而得 x k + 1 = φ ( x k ) x_{k+1}=φ(x_k) xk+1=φ(xk)。
\,\,\,\,\,\,\,\,\,\, 2.判断迭代序列是否收敛:由迭代公式 x k + 1 = φ ( x k ) x_{k+1}=φ(x_k) xk+1=φ(xk)可得序列 { x k } k = 0 ∞ \{x_k\}^\infty_{k=0} { xk}k=0∞。需要首先判断该序列是否收敛,即是否存在α使得 lim k − > ∞ x k = α \lim \limits_{k->\infty} x_k=α k−>∞limxk=α。
\,\,\,\,\,\,\,\,\,\, 3.求出根的数值解:因为有 x k + 1 = φ ( x k ) x_{k+1}=φ(x_k) xk+1=φ(xk),则方程两边同时对k趋近于无穷取极限,有 lim k − > ∞ x k = α \lim \limits_{k->\infty} x_k=α k−>∞limxk=α,而且 lim k − > ∞ φ ( x k ) = φ ( α ) \lim\limits_{k->\infty} φ(x_k)=φ(α) k−>∞limφ(xk)=φ(α)。进而可知α满足等式x=φ(x),从而α是迭代函数φ(x)的不动点(不断迭代的过程中,函数值在此处保持不动)。因此α也是方程f(x)=0的一个根。 - (二)多点迭代法的基本形式:
\,\,\,\,\,\,\,\,\,\, 建立迭代公式时,迭代函数是一个n元函数。即 x k + 1 = φ ( x k , x k − 1 , x k − 2 , . . . , x k − n + 1 ) x_{k+1}=φ(x_k,x_{k-1},x_{k-2},...,x_{k-n+1}) xk+1=φ(xk,xk−1,xk−2,...,xk−n+1)。其余与单点迭代法相似。 - (三)迭代法的主要性能指标:
- 1.收敛:
- 全局收敛性:设α为 f ( x ) = 0 f(x)=0 f(x)=0的根,如果 ∀ x 0 ∈ [ a , b ] , [ a , b ] 是 f ( x ) 的 定 义 域 , x 0 是 迭 代 法 的 初 始 值 \forall x_0\in[a,b],[a,b]是f(x)的定义域,x_0是迭代法的初始值 ∀x0∈[a,b],[a,b]是f(x)的定义域,x0是迭代法的初始值,由迭代法产生的序列都收敛于根α则称该迭代法是全局收敛的。
- 局部收敛性: 设α为 f ( x ) = 0 f(x)=0 f(x)=0的根 , 如果存在α 的某个邻域Δ:|x-α|≤δ,对任意初值 x 0
- 1.收敛: