非线性现象广泛存在物质时间与社会生活,由于自然现象和实际问题的复杂性,无法求解出一般方程的准确解,因此很有必要进行数值解的研究,本文介绍非线性方程求根问题的数值方法。包括二分法、牛顿迭代法以及他们的收敛性判断和收敛速度的求解。
一.方程求根的理论依据
(1)代数基本定理:任何一个一元复系数方程式都至少有一个复数根。
(2)零点定理,f(x)在(a,b)区间连续,且f(a)·f(b)<0。则称f(x)在区间[a,b]内有根。
二.直接法
2.1二分法
2.1.1 二分法的理论依据:
设函数f(x)在[a,b]上连续、严格单调,且f(a)·f(b)<0,则在[a,b]有且仅有一根。
2.1.2二分法的基本思想:
反复对分区间,从而逐步缩小有根区间,直至满足精度为止。
2.1.3二分法的步骤
取区间[a,b]的中点=(a+b)/2,计算f(
),若 f(
)=0,即为根,否则:
(1)若f(a)·f()<0,则根在区间[a,
]中, [a,
]为新区间,并对新区间进行二等分。 (2)若 f(b)·f(
)<0,则根在区间[b,
]中, [b,
]为新区间,并对新区间进行二等分。
如此反复下去,无限逼近根所在的区间,确定二分的次数,得到一定精确值的解。
2.1.4二分法的终止条件
已知区间[a,b],二分的次数k,给定的精度 。解下面不等式得二分次数k
三.迭代法
3.1不动点迭代
迭代法即逐次逼近方法,是数值计算中一类典型方法,被用于数值计算得各个方面中。
3.1.1迭代过程及基本思想
将方程f(x)=0,写成等价的迭代形式x= 。由此确定了迭代法:
例如:
能够写成:
构造迭代算法:
3.1.2收敛性判断
设方程f(x)=0在[a,b]上存在唯一解,x=(x) 是方程的等价形式
如果满足映内压缩性,则迭代公式对于区间任意取初值收敛:
(1)映内性:对任意的x[a,b],对应
[a,b]。
(2)压缩性:,其中L是
在a点的取值。
局部收敛性与收敛阶
局部收敛,在x*处邻近连续,且满足:
,则迭代过程
在邻近具有局部的收敛性。
收敛速度
收敛阶,
(1) 若压缩因子0<L<1,则迭代过程为线性收敛.
(2) 若压缩因子0=L,则迭代过程为超线性收敛.
(3) 一般地,收敛阶p越大,迭代过程收敛越快.
三.牛顿迭代法及其收敛性
牛顿迭代法是对f(x)进行一阶的泰勒展开,则有牛顿公式:
收敛性,若f(x)在单根x*附近存在连续的二阶导数,且初值充分接近x*,则牛顿迭代法收敛。