数值计算方法作业-二分法

编程实现,用二分法求 f ( x ) = x 3 + 4 x 2 10 f(x)=x^3+4x^2-10 在区间 [ 1 , 2 ] [1,2] 内的一个实根,要得到精确到小数点后三位位精确值的近似
注意:编程的循环次数

易证该方程在区间内单调递增,所以以此来确定二分方向

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

发布了84 篇原创文章 · 获赞 12 · 访问量 2909

猜你喜欢

转载自blog.csdn.net/qq_43294914/article/details/104541131