转载:C/C++牛顿迭代法求非线性方程近似根


最近用到牛顿迭代法,找了很多,感觉下文比较通俗易懂。 

原文地址:https://blog.csdn.net/qq_23908539/article/details/51454562?locationNum=10&fps=1

把f(x)在x0附近展开成泰勒级数f(x) = f(x0)+f'(x0)(x-x0)+f''(x0)/2!*(x-x0)^2+...然后取其线性部分,作为非线性方程f(x) = 0的

近似方程,即泰勒展开的前两项,则有
f(x) = f'(x0)x - x0*f'(x0) + f(x0) = 0
f'(x0)x = x0*f'(x0) - f(x0)
x = x0 - f(x0)/f'(x0)
得到牛顿的一个迭代序列:
->x(n+1) = x(n)-f(x(n))/f'(x(n))

 

例:求方程f(x) = 2*x^3-4*x^2+3*x-6在x = 1.5附近的根,求出总的迭代次数和每次的近似根。

#include <iostream>
#include <cmath>
using namespace std;

double fun1(double x)
{
    return 2*x*x*x-4*x*x+3*x-6;
}

double fun2(double x)
{
    return 6*x*x-8*x+3;
}

int main()
{
    double f1,f2,x,d;
    int cnt = 0;
    x = 100;
    do {
        f1 = fun1(x);//方程的值 
        f2 = fun2(x);//方程的导数 
        d = f1/f2;//"斜率" 
        x -= d;//更新方程的值 
        cnt ++;
        printf("第%d次迭代方程的值为:%lf\n",cnt,fun1(x));
        cout<<"当前的近似根为:";
        printf("%.lf\n",x);
    } while(fabs(d) > 1e-5);
    cout<<"总迭代次数:"<<cnt<<endl;
    cout<<"最终的近似根为:";
    printf("%lf\n",x);
    return 0;
}
转自:https://blog.csdn.net/qq_23908539/article/details/51454562?locationNum=10&fps=1
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值