数值分析——Newton迭代法求解方程附Matlab程序

Newton迭代法应用

1. 给定初值 $x_0$及容许误差 $\epsilon$，编制Newton法解方程 $f(x)=0$根的通用程序；
Newton 迭代格式的迭代函数为：
$\phi (x)=x-\frac{f(x)}{f^{\prime }(x)}$即为$x_{k+1}=x_k-\frac{f(x_k)}{f^{\prime }(x_k)}$

Newton迭代法程序如下：
文件名：Newton.m

x0=1;%初值大小
newton=@(x)x - (x^3/3-x)/(x^2-1);  %此处的f(x)=x^3/3-x=0
k=0 ; err=1.000;
phi=10^(-4);%允许误差大小
disp('k         x_k             err');
while(err>phi)
    x1=newton(x0);
    err=abs(x1-x0);
    x0=x1;
    k=k+1;
  end
fprintf('%d  %1.15f   %1.15f\n',k,x1,err);%每次迭代的结果和误差
%输出对于不同初值迭代次数、迭代结果、误差

存在问题：特别对于x=0，时牛顿迭代法可能会失去其作用，并不能通过较小的 $\delta$ 来找到收敛于0的牛顿序列，导数是否在x=0这一点可导

2. 给定方程$f(x)=x^3-x=0$。容易知道其有三个根 $x_1^{\ast }=-\sqrt{3}$, $x_2^{\ast }=0$ ,$x_3^{\ast }=\sqrt{3}$。

1. 由Newton方法的局部收敛法可知：存在 $\delta >0$, 当 $\delta >0$时， Newton迭代序列收敛于根 $\delta >0$,， 是确定尽可能大的 $\delta >0$,；
   由牛顿的局部收敛性，即在最大区间$(x^{\ast }-\delta ,x^{\ast }+\delta )$内，通过牛顿迭代，存在收敛序列使得$x_n收敛到x^{\ast }$，$x^{\ast }$=$0$,由区间的对称性可以在 $[0,1]$内找到最大的 $\delta$ ，使得$x_n能够收敛到x^{\ast }$
   程序如下：

文件名：Converge.m：
max=0.000000;
for i=0:0.0001:1			%将区间按0.0001分割
    x1=i;			%初值大小
    newton=@(x)x - (x^3/3-x)/(x^2-1);%f(x)=x^3/3-x
    k=0;err=1.000;
    phi=10^(-4);					%允许误差大小
    while(err>phi)
        x2=newton(x1);
        err=abs(x2-x1);
        x1=x2;
        k=k+1;
        %fprintf('%d  %1.15d   %1.15f\n',k,x1,err);
        if x1==0
            max=i;        %记录的最大值
        end
    end
end
fprintf('max= %1.15f\n',max)

运行程序：
max= 0.774300000000000
由结果可知当误差取$phi=10^{-4}时$，确定的最大的 $\delta =0.77430$，当$x_0\in (-0.77430,0.77430)时$，Newton迭代序列收敛到根$x_2^{\ast }$

2. 试取若干初始值， 观察当 $x_0\in (-\infty ,-1)$、$(-1,-\delta )、(-\delta ,\delta )$、$（\delta ,1）、(1,+\infty )$ 时，$Newton$序列是否收敛以及收敛于那一个根。
   分别在$x_0\in (-\infty ,-1)$、$(-1,-\delta )$、$(-\delta ,\delta )$、$（\delta ,1）$、$(1,+\infty )$ 时，$Newton$，取初值(注意：在取初值时需要修改源程序中的$x_0$可能在调试中较麻烦，可以建立函数实现输入初值$x_0$)

A.$在x_0\in (-\infty ,-1)$内取初值x=-10，利用（1）的Newton.m程序,运行结果：

即初值x=-10时，通过9次牛顿迭代可以得到在误差范围内的根，误差为err=0.000000241956186，牛顿序列收敛到$x_1^{\ast }=-\sqrt{3}$

B.在$(-1,-\delta )$内取初值x=-0.9，利用（1）的Newton.m程序,运行结果：

即初值x=-0.9时，通过6次牛顿迭代可以得到在误差范围内的根，误差err=0.000003448013598牛顿序列收敛到$x_3^{\ast }=\sqrt{3}$

C.在 $(-\delta ,\delta )$ 内取初值x=-0.6，利用（1）的Newton.m程序,运行结果：

即初值x=-0.6时，通过4次牛顿迭代可以得到在误差范围内的根，误差err=0.000000341186658。牛顿序列收敛到 $x_2^{\ast }=0$

D.在$（\delta ,1）$内取初值x=0.9，利用（1）的Newton.m程序,运行结果：

即初值x=0.6时，通过6次牛顿迭代可以得到在误差范围内的根，误差err=0.000003448013598。牛顿序列收敛到$x_3^{\ast }=-\sqrt{3}$

E.在$(1,+\infty )$ 内取初值x=100，利用（1）的Newton.m程序,运行结果：

即初值x=0.6时，通过14次牛顿迭代可以得到在误差范围内的根，误差err= 0.000088593081718。牛顿序列收敛到$x_3^{\ast }=\sqrt{3}$