牛顿迭代法（Newton's method）又称为牛顿-拉夫逊（拉弗森）方法（Newton-Raphson method），它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。
牛顿迭代公式
设r是\(f(x)=0\)的根，选取\(x_0\)作为r的初始近似值，过点\((x_0,f(x_0))\)
,做曲线
\(y=f(x)\)的切线L，L的方程为\(y=f(x_0)+f’(x_0)(x-x_0)\)
，求出L与x轴交点的横坐标
\[x_1=x_0-\frac{f(x_0)}{f’(x_0)}\]
称\(x_1\)为r的一次近似值。过点\((x_1,f(x_1))\)
做曲线
\(y=f(x)\)的切线，并求该切线与x轴交点的横坐标
\[x_2=x_1-\frac{f(x_1)}{f’(x_1)}\]
称\(x_2\)为r的二次近似值。重复以上过程，得r的近似值序列，其中，
\[x_{n+1}=x_n-\frac{f(x_n)}{f’(x_n)}\]
称为r的\(n+1\)次近似值，上式称为牛顿迭代公式。
用牛顿迭代法解非线性方程，是把非线性方程\(f(x)=0\)线性化的一种近似方法。把
\(f(x)\)在点
\(x_0\)的某邻域内展开成泰勒级数
\[f(x)=f(x_0)+f’(x_0)(x-x_0)+\frac{f’’(x_0)(x-x_0)^2}{2!}+…+\frac{f^{(n)}(x_0)(x-x_0)^n}{n!}+R_n(x)\]
取其线性部分（即泰勒展开的前两项），并令其等于0，即
\(f(x_0)+f’(x_0)(x-x_0)=0\)，以此作为非线性方程
\(f(x)=0\)的近似方程，若\(f’(x_0)\neq0\)，则其解为
\[x_1=x_0-\frac{f(x_0)}{f’(x_0)}\]
这样，得到牛顿迭代法的一个迭代关系式：
\[x_{n+1}=x_n-\frac{f(x_n)}{f’(x_n)}\]
从下面的图中，我们可以看到牛顿迭代的几何意义，每次迭代，都会更加逼近\(f(x)=0\)的解。
下面用牛顿迭代法在matlib解出方程\(x^2+2xe^x+e^{2x}=0\)的根，首先画出函数的图像，猜测根的大致位置。
函数图像如下图所示：
近似结果：
x=-1:0.01:1;
y= x.^2+2*x.*exp(x)+exp(2*x)
plot(x,y);
grid on;
clc;clear;
%syms x;
%diff(x^2+2*x*exp(x)+exp(2*x),x,1)
%clear;
x=0.0
for i=1:100
x1=x-(x^2+2*x*exp(x)+exp(2*x))/(2*x + 2*exp(2*x) + 2*exp(x) + 2*x*exp(x))
if(abs(x1-x)>0.0001)
x=x1;
else
break;
end
end
i
View Code