计算方法实验（四）：牛顿迭代法

Newton迭代法数学原理

求非线性方程$f(x)=0$的根$x^{\ast }$，牛顿迭代法计算公式

$$x_0=\alpha$$

$$x_{n+1}=x_n-\frac{f(x_n)}{f^{{}^{\prime }}(x_n)}$$

$$n=0,1,\cdots$$

一般地，牛顿迭代法具有局部收敛性，为保证迭代收敛，要求，对充分小的$\delta >0$，$\alpha \in O(x^{\ast },\delta )$。如果$f(x)\in C^2[a,b]$，$f(x^{\ast })=0$，$f^{{}^{\prime }}(x^{\ast })\ne 0$，那么，对充分小的$\delta >0$，当$\alpha \in O(x^{\ast },\delta )$时，由牛顿迭代法计算出的$\{x_n\}$收敛于$x^{\ast }$，且收敛速度是2阶的；如果$f(x)\in C^m[a,b]$，$f(x^{\ast })=f^{{}^{\prime }}(x^{\ast })=\cdots =f^{(m-1)}(x^{\ast })=0$，$f^{(m)}(x^{\ast })\ne 0(m>1)$，那么，对充分小的$\delta >0$，当$\alpha \in O(x^{\ast },\delta )$时，由牛顿迭代法计算出的$\{x_n\}$收敛于$x^{\ast }$，且收敛速度是1阶的；

问题

利用牛顿迭代法求$f(x)=0$的根

输入：初值$\alpha$，精度${\epsilon }_1,{\epsilon }_2$，最大迭代次数$N$

输 出：方程$f(x)=0$根$x^{\ast }$的近似值或计算失败标志

程序流程

核心代码

#include 
#include 
#include 
using namespace std;

double x, e1, e2;
int n;

double f(double x) { return cos(x) - x; }
double df(double x) { return -sin(x) - 1; }

int main() {
    scanf("%lf%lf%lf%d", &x, &e1, &e2, &n);
    for (int i = 1; i <= n; i++) {
        double F = f(x), DF = df(x);
        if (fabs(F) < e1) {
            printf("%lf", x);
            return 0;
        }
        if (fabs(DF) < e2) {
            printf("Failed");
            return 0;
        }
        double x1 = x - F / DF;
        double tol = fabs(x - x1);
        if (tol < e1) {
            printf("%lf", x1);
            return 0;
        }
        x = x1;
    }
    printf("Failed");
    return 0;
}

详细报告