牛顿迭代法求方程组的根_直接迭代法

(26) 2023-08-29 15:12

Hi,大家好,我是编程小6,很荣幸遇见你,我把这些年在开发过程中遇到的问题或想法写出来,今天说一说牛顿迭代法求方程组的根_直接迭代法,希望能够帮助你!!!。

1.牛顿迭代法求方程根
牛顿迭代法求方程组的根_直接迭代法_https://bianchenghao6.com/blog__第1张

[图片摘自知乎如何通俗易懂地讲解牛顿迭代法求开方(数值分析)]牛顿迭代法求方程组的根_直接迭代法_https://bianchenghao6.com/blog__第2张
y=x3+2x2+3x+4 在1附近的实根
x0 是上一次求出的近似根,在开始是根据题设 x0 =1(题目希望求x在1附近的一个实根,因此第一次的近似值可以设定为1)。

#include <iostream>
#include <stdio.h>
#include <math.h>
int main()
{ 
   
    double x0;
    double x=1,f,f1;
    do
    { 
   
        x0 = x;
        f = x0 * x0* x0 + 2* x0* x0 + 3* x0 + 4;
        f1 = 3 * x0* x0 + 4* x0 +3;
        x = x0 - f / f1;//公式x=x0-f(x0)/f’(x0).
    }while(fabs(x - x0)>= 1e-3);
    printf("%lf\n",x);
    return 0;
}

2.用迭代法求平方根----即求方程y=x^2-a的根
求平方根的迭代公式为:X(n+1)=(Xn+a/Xn) /2。x=(x0+a/x0)/2;
公式推导方法同上

#include <iostream>
#include <stdio.h>
#include <math.h>
int main()
{ 
   
    double x0,a;
    double x=1,f,f1;
    scanf("%lf",&a);
    do
	{ 
   
		x0 = x;
		x=(x0+a/x0)/2.0;
	}while(fabs(x - x0)>1e-5);
    printf("%lf",x);

    return 0;
}

3.二分法求方程根
x1,x2代表区间左右边界,f(x1)*f(x2)异号即根在x1,x2之间,取x1x2中间点xx 再判断根在x1~xx之间,还是在xx ~x2之间

 //伪代码
 do
    { 
   
        xx=(x1+x2)/2;
        if(f(xx)*f(x1)>0)
            x1=xx;
        else
            x2=xx;
    }
    while(fabs(f(xx))>=1e-7);

今天的分享到此就结束了,感谢您的阅读,如果确实帮到您,您可以动动手指转发给其他人。

上一篇

已是最后文章

下一篇

已是最新文章

发表回复