重庆分公司,新征程启航
为企业提供网站建设、域名注册、服务器等服务
二分法计算函数f(x)=x*x*x*x+2*x*x*x-x-1;
成都创新互联主营新林网站建设的网络公司,主营网站建设方案,重庆APP开发,新林h5重庆小程序开发搭建,新林网站营销推广欢迎新林等地区企业咨询
本程序在turbo c或c++下编译
#include "stdio.h"
#include math.h
float f(float x)
{float y;
y=x*x*x*x+2*x*x*x-x-1;
return y;
}
void main()
{float a=0,b=0,h,y,x;
int k,n0;
printf("please input qujian a and b");
scanf("%f%f%d",a,b,n0); /*输入含根区间a,b,循环次数n0 */
for(k=0;k=n0;k++)
{ x=(a+b)/2;
h=(b-a)/2;
y=f(x);
if(h10e-6||fabs(y)10e-6)
{ printf("k=%d,x=%f,y=%f",k,x,y);
break; } /*输出分半次数k,函数的根x,及x对应的函数值.*/
else
{if(f(a)*f(x)0)
b=x;
else a=x;
}
}
}
举个例子:
//二分查找法//
# includestdio.h
void main()
{
int a[16],i,num,flag=0,top,bottom,mid;
//定义一个一维数组a[16]用来存放供查找用的数据,但只用a[1]——a[15]//
//num用来放要查找的数据,flag是表示是否找到的开关变量,top表示查找的起始位置,bottom表示查找的终止位置,mid表示top与bottom的中间位置//
char goon;
//变量goon为'y'或'Y'时表示继续下一轮查找,否则终止程序//
printf("请输入第1个数字:\n");
scanf(" %d",a[1]);
//依次输入第二到第十五个数,并要求输入的数递减//
for(i=2;i=15;i++)
{
printf("请输入第%d个数字:\n",i);
scanf(" %d",a[i]);
if(a[i]=a[i-1])
{
printf("请再次输入,它应该比上一个数小:\n");
scanf(" %d",a[i]);
}
}
//输出刚才输入的数//
printf("你刚才输入的数是:\n");
for(i=1;i=15;i++)
printf(" %d",a[i]);
printf("\n");
//查找循环开始//
do
{
printf("现在请输入你要查找的数:\n");//输入想要查找的数//
scanf(" %d",num);
top=15;
bottom=1;
mid=15/2+1;
if(numa[1] || numa[15])//如果要查找的数据不在规定范围内,令flag=1,输出超出范围的信息//
{
flag=1;
printf("你所要查找的数字不在范围内!\n");
}
while(flag==0 (top-bottom)0)//如果在规定的范围内,开始二分法查找//
{
if(num==a[mid])//找到所需数据,退出本层循环//
{
printf("你所要查找的数字是第%d个。\n",mid);
flag=1;
}
else if(numa[mid])//如果要查找的数据比a[mid]大,在前半数组查找//
{
top=mid+1;
mid=(top+bottom)/2;
}
else //如果要查找的数据比a[mid]小,在后半数组查找//
{
bottom=mid-1;
mid=(top+bottom)/2;
}
}
if(flag==0)//如果未找到数据,输出找不到的信息//
printf("无法找到你要找的数字!\n");
printf("是否继续查找?(Y/N):\n");//询问是否开始下一轮查找//
scanf(" %c",goon);
}while(goon=='y' || goon=='Y');
}
这段代码是求解方程f(x)=0在区间[-10,10]上的根的数值解。
方法的思想就是:一直选取区间中间的数值,如果发现中间的函数值与一侧函数值,异号,那么说明解在这个更小的区间中,采用eps=1e-5作为区间的极限大小,通过迭代的方法求解这个方程的数值解。
所以了解了上述思想,那么else if(f(a)*f(c)0) b=c; 说明的是 f(a)和f(c)异号,那么使用b=(a+b)/2缩小迭代区间,继续迭代;同理else a=c;说明f(a)和f(c)同号,那么使用a(a+b)/2缩小迭代区间,继续迭代!