编程实现,用二分法求 在区间 内的一个实根,要得到精确到小数点后三位位精确值的近似
注意:编程的循环次数
易证该方程在区间内单调递增,所以以此来确定二分方向
C++实现:
#include <bits/stdc++.h>
/****************
author:WANG_zibi
time:2020/2/25
title:二分法求方程近似值
course:数值计算方法
****************/
using namespace std;
const int N = 1e6+5;
typedef long long ll;
typedef unsigned long long ull;
const int INF= 0x3f3f3f3f;
#define f(i, a, b) for (int i = (a); i <= (b); ++i)
const double eps = 1e-4;//精度
double func(double x) //函数
{
return x*x*x+4*x*x-10.0;
}
int main()
{
double left,right,mid,answer,sum;//左右区间,答案等变量生命
left=1,right=2;
mid=(left+right)/2.0;
sum=func(mid);
answer=mid;
while(abs(sum)>eps)
{
if(sum>0)
{
right=mid;
mid=(right+left)/2.0;
sum=func(mid);
answer=mid;
}
else{
left=mid;
mid=(right+left)/2.0;
sum=func(mid);
answer=mid;
}
}
cout<<"最终的答案为:"<<answer<<'\n';
return 0;
}
MATLAB实现:
%author:WANG_zibi
%date:2020/2/25
%course:数值计算方法
%name:二分法求方程近似解
%方程
l=1.0;
r=2.0;
x=(l+r)/2.0;
eps=0.0001;
f=x^3+4*(x^2)-10;
while abs(f)>eps
if(f>0)
r=x;
x=(l+r)/2.0;
f=x^3+4*(x^2)-10;
end
if(f<0)
l=x;
x=(l+r)/2.0;
f=x^3+4*(x^2)-10;
end
end
ans=x