1 题目要求
二分法求函数根的原理为:如果连续函数f(x)在区间[a,b]的两个端点取值异号,即f(a)f(b)<0,则它在这个区间内至少存在1个根r,即f®=0。
二分法的步骤为:
检查区间长度,如果小于给定阈值,则停止,输出区间中点(a+b)/2;否则
如果f(a)f(b)<0,则计算中点的值f((a+b)/2);
如果f((a+b)/2)正好为0,则(a+b)/2就是要求的根;否则
如果f((a+b)/2)与f(a)同号,则说明根在区间[(a+b)/2,b],令a=(a+b)/2,重复循环;
如果f((a+b)/2)与f(b)同号,则说明根在区间[a,(a+b)/2],令b=(a+b)/2,重复循环。
本题目要求编写程序,计算给定3阶多项式
f(x)=a3x^3 +a2x2 +a1x+a0
在给定区间[a,b]内的根。
输入格式:
输入在第1行中顺序给出多项式的4个系数a3 、a2 、a1 、a0 ,
在第2行中顺序给出区间端点a和b。题目保证多项式在给定区间内存在唯一单根。
输出格式:
在一行中输出该多项式在该区间内的根,精确到小数点后2位。
2 样例
输入样例:
3 -1 -3 1
-0.5 0.5
输出样例:
0.33
3 分析
(1)首先去掉题目中的干扰句子——根据输入格式中“题目保证多项式在给定区间内存在唯一单根。”这句话,去掉 二分法的步骤为 中的“否则 如果f(a)f(b)<0,”,不然会出错。
(2)检查区间长度,如果小于给定阈值,则停止——
这句话中有“阈值”, 但样例输入中没有阈值这个变量,在实际操作中任意给一个趋近于0,但不等于0的值,比如0.001
(3)根据步骤中“令 a=(a+b)/2,重复循环;”等类似 表述,采用递归进行实现
4 代码
(1)我的方法——递归形式
#include <stdio.h>
#include <math.h>
double polynomial(double a3,double a2,double a1,double a0,double x)
{
double f=a3*pow(x,3)+a2*pow(x,2)+a1*x+a0;
return f;
}
double binaryFind(double a3,double a2,double a1,double a0,double a,double b)
{
double x=(a+b)/2;
double fa=polynomial(a3,a2,a1,a0,a);
double fb=polynomial(a3,a2,a1,a0,b);
if((b-a)<=1e-3)
return x;
//if(fa*fb<0)
//{
double faplusb_2=polynomial(a3,a2,a1,a0,x);
if (faplusb_2==0)
{
return x;
}
// else
// {
if(fa*faplusb_2>0)
return binaryFind(a3,a2,a1,a0,x,b);//
if(fb*faplusb_2>0)
return binaryFind(a3,a2,a1,a0,a,x);
// }
//}
}
int main(){
double a3,a2,a1,a0;
double a,b;
scanf("%lf %lf %lf %lf",&a3,&a2,&a1,&a0);
scanf("%lf %lf",&a,&b);
double result=binaryFind(a3,a2,a1,a0,a,b);
printf("%.2f",result);
return 0;
}
注意:
递归外不能有while包围
递归要有返回值,且调用递归声明时,前面要return
更可能优先出现的if条件放在前面
(2)非递归形式 : 链接.
#include <iostream>
#include <cmath>
using namespace std;
double a3, a2, a1, a0;
double f(double x)
{
return a3*pow(x,3) + a2*pow(x,2) + a1*x + a0;
}
double Find(double a, double b)
{
while((b-a) >= 0.0001) //注意这个条件,如果区间差距很小默认该区间重合,否则会运行超时
{
double mid = (a+b)/2;
if(f(mid) == 0) return mid;
if(f(a)*f(mid) > 0) a = mid;
else b = mid;
}
return a;
}
int main()
{
double a, b;
scanf("%lf%lf%lf%lf%lf%lf",&a3,&a2,&a1,&a0,&a,&b);
printf("%.2lf",Find(a,b));
return 0;
}
5 总结
(1)递归 链接:
使用递归的三个步骤
递归算法转换为非递归算法的技巧.
(2)double 和 0 的比较 问题
(3)