二分法求多项式的一个根

最新推荐文章于 2023-01-03 20:18:03 发布

Durant_kevin

最新推荐文章于 2023-01-03 20:18:03 发布

阅读量1.6k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： C++ 数据结构与算法文章标签：二分法多项式的根

本文链接：https://blog.csdn.net/u010367506/article/details/46310459

数据结构与算法同时被 2 个专栏收录

23 篇文章

订阅专栏

C++

6 篇文章

订阅专栏

本文介绍了二分法求解多项式根的数学原理，详细阐述了求根过程，并提供了相应的代码实现，帮助理解如何通过二分法找到连续函数在一个区间内的根。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

数学原理

二分法求根的数学原理：如果连续函数f(x)在区间[a,b]的两个端点上取值异号，则在该该函数在该区间上必有一个根。

解步骤

二分法求解步骤与二分查找非常相似。具体如下：

1.检查区间的长度，如果小于阈值，则返回中间值，mid=(a+b)/2。

2.求中间值对应的函数值，f(mid)。

3.如果f(mid)==0,返回mid。

4.如果f(mid)与f(a)同号，即f(mid)*f(a)>0,则令a=mid。返回步骤一。

5.如果f(mid)与f(b)同号，即f(mid)*f(b)>0,则令b=mid。返回步骤一。

代码

#include<iostream>
using namespace std;
double THRESHOLD=0.001;
int coefficients[100];
int n=2;//系数个数-1 
double f(double x){
       double res=0;
       double temp=1;
       for(int i=0;i<n;i++){
           res+=res+coefficients[i]*temp;
           temp*=x;
       }
       return res;
}
double root(double a, double b){
       if(f(a)==0)
          return a;
       if(f(b)==0)
          return b;
       double mid;
       while(b-a>=THRESHOLD){
          mid=(a+b)/2;
          if(f(mid)==0)
             return mid;
          if(f(a)*f(mid)>0)
              a=mid;
          else if(f(b)*f(mid)>0)
              b=mid;                             
       }
       return (a+b)/2;
}
int main()
{
     double a, b;
     for(int i=0; i<n; ++i ) {
         cin>>coefficients[i];
     }
     cin>>a>>b;
     
     printf( "%.2lf", root( a, b ) );
     return 0;
}

算法流程很简单，要做到bug-free。