非线性方程的数值解法——二分法
考虑单变量线性方程:
f
(
x
)
=
0
f(x)=0
f(x)=0
第一步:先通过隔离法选取初始的含根区间。若已知f(x)在(a,∞)上有根,取步长h>0,计算f(a+ih)*f(a+(i+1)h),i=0,1… 直到出现j,使得f(a+jh)*f(a+(j+1)h)<0,含根区间为(a+jh,a+(j+1)h),下文记为(a,b)。
第二步:二分法计算根。
二分法总是收敛的。实际计算中,对于给定精度epsilon>0,当区间长度|bk-ak|<2epsilon时终止运算,取x*≈x_k+1=1/2(ak+bk),此时误差小于半个区间长度:
∣
x
k
+
1
−
x
∗
∣
<
=
1
2
(
b
k
−
a
k
)
=
1
2
k
+
1
(
b
−
a
)
<
ϵ
|x_{k+1}-x*|<=\frac{1}{2}(b_k-a_k)=\frac{1}{2^{k+1}}(b-a)<\epsilon
∣xk+1−x∗∣<=21(bk−ak)=2k+11(b−a)<ϵ
所以可以用上式求出二分的次数k,使算法在有限步终止。
matlab实现:
与上面不同的是,给定已知的求根区间(a,b);结束计算的条件是函数值的差小于误差而非根的误差。
funciton[r,k] = half(fun,a,b,e)
%Inputs -fun 需要求根的函数
-a,b 含根区间的左右端点
-e 容忍误差
%outputs -r 根
-k 二分次数
while abs(fun(b)-fun(a))>e %函数值的差小于等于容忍误差时停止计算
c = 0.5*(a+b);
if fun(c)==0
k = k+1;
break;
end
if fun(a)*fun(c)<0
b = c;
elseif fun(b)*fun(c)<0
a = c;
end
k = k+1; % 二分次数
end
r = c; % 根
%函数
fun = inline('x^3-2*x-5','x');
%调用half
[r,k] = half(fun,2,3,10^-3);
二分法的优点是算法简单,总是收敛。
缺点是收敛速度慢,一般用于求近似根。
参考资料:《数值计算方法》;倪勤等;高等教育出版社