二分法-Newton法-割线法求解教学教材.docxVIP

二分法-Newton法-割线法求解 非线性方程数值求解的算法实现 二分法： 1.理论基础： 对于方程f（x）=0，如果在区间[a,b]上至少有一个根，就称[a,b]是方程的一个有根区间。例如，若f(x)连续，且f(x)*f(b)0，由介值定理可知[a,b]是一个有根区间。如果在[a.b]上方程有且只有一个根，那就把方程的根隔离出来了，这时若把有根区间不断缩小，便可逐步得到根的近似值。 设方程f(x)=0的一个有根区间为[a,b]，且在区间[a,b]上只有一个根，满足f(x)*f(b)0。可以用区间对分的方法形成有限区间的序列。令，对于区间，区间中点。检验的符号，若，则取新的有限区间，有根区间向左压缩。若，则取为，有根区间向右压缩。 以上产生的有根区间的序列，满足，其中各区间的长度等于上一个区间的一半，区间中点的序列就是方程的根的近似解序列。不难分析，对有。而是的中点，满足。这保证了。 代码： function y=f(x) y=x^3-x-1; end 主程序： clear all clc a=input(输入下界 a=); b=input(输入上界 b=); eps=input(输入误差范围 eps=); distance=b-a; num=0; format long while(distanceeps) x0=(a+b)/2; if f(a)*f(x0)0 b=x0; else if f(a)*f(x0)0 a=x0; else a=x0; b=x0; end distance=distance/2; num=num+1; end x=(a+b)/2; end 3.运行结果： Newton法： 1.理论基础： 为了求解方程组的根。设已有一个近似值，如果存在且连续，由Taylor展开式得，其中在和之间。因为f(x*)=0,如果，可得。如果把右端含的项略去，剩下的两项就作为x*新的一个近似值，记为，即。 2.代码 ： function Newton(x0,nums) x=x0; for i=1:nums x1=x-(x^3+4*x^2-10)/(3*x^2+8*x); x=x1; end x 3.运行结果 割线法： 理论基础： Newton法每步不但要计算函数值，还要计算导数值。有时的计算比较麻烦，可以用点上的差商近似导数，即，Newton法迭代变换为，。这就是割线法的计算公式。 代码： function GeXian(x0,x1,nums) x(1)=x0; x(2)=x1; for i=3:nums x(i)=x(i-1)-(x(i-1)^3+4*x(i-1)^2-10)*(x(i-1)-x(i-2))/((x(i-1)^3+4*x(i-1)^2-10)-(x(i-2)^3+4*x(i-2)^2-10)); end x 3.运行结果： 比较： 在论文中求解的例子可以看出，二分法需要14次迭代才能精确到小数点后四位有效数字。而Newton法和割线法分别需要3次和4次能达到如此的效果。很明显Newton法和割线法优于二分法。这是因为二分法是线性收敛的，而Newton法至少是二阶收敛的。割线法是超线性收敛的，而且收敛速度小于2.所以，在本例中，Newton法是优于割线法的。 二分法的收敛速度虽然较慢，但它是全局收敛的，而Newton法和割线法是局部收敛的。