鱼C论坛

 找回密码
 立即注册
查看: 687|回复: 1

一元三次方程求解用的二分,为什么结果不对

[复制链接]
发表于 2023-11-10 18:17:13 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
// 这题其实得用二分来做,为实数的二分
#include <bits/stdc++.h>
using namespace std;
double a,b,c,d;
double solve(double x)
{
    return a*x*x*x+b*x*x+c*x+d;
}
int main()
{
    cin>>a>>b>>c>>d;
    for(double i=-100;i<=100;i++)  //即在-100到100中一个一个区间去判断
    {
        double x1=i,x2=i+2.0;
        if(solve(x1)==0.0)printf("%.2lf",i); //如果带进去x1==0,则说明就是此方程的根
        if(solve(x1)*solve(x2)<0)//标明此处必有一个根,然后去判断
        {
            while(x2-x1>=0.001)
            {
                double mid=(x1+x2)/2.0;
                if(solve(x1)*solve(mid)<=0.0)x2=mid;
                else x1=mid;


            }
            printf("%.2lf",x1);
        }

    }
    return 0;


}
其中有一个地方不懂x2-x1>=0.001
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

发表于 2023-11-10 18:17:22 | 显示全部楼层
这段代码中的 x2-x1>=0.001 是用来控制二分法的迭代次数的条件。在二分法中,我们希會在每次迭代中减小区间的长度,直到满足精度要求。这里的0.001是一个预先设定的精度值,表示当区间长度小于等于0.001时,就认为已经找到了足够接近的解。因此,当 x2-x1 的值小于等于0.001时,就会退出迭代,输出当前的 x1 值作为方程的解。

然而,在实际应用中,这种固定的迭代终止条件并不总是适用的。在一些情况下,可能需要根据具体问题来调整迭代的终止条件,以确保找到足够精确的解。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2024-6-15 13:32

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表