杜拉瓜耶 2020-04-25 23:09 采纳率: 0%
浏览 482

怎么用二分法求在一定区间内有多个根的一元三次函数?

1、本人小白,刚接触c语言(用的是vs2019)。练习题:用二分法求方程在(-10,10)范围内的根:y=x^3+6x^2+3x-6,收敛精度为10^-5。二分法好像只能计算出一个零点,就想着用循环的方式,把第一次循环得到的x=0.7465代替b的位置,再进行二分法。但是好像不太行,请问我这问题出在什么地方,要怎么改,或者有没有更简单的方法?注:已知这个方程在(-10,10)内有三个根。

图片说明

#include
#include
int main()
{
double a, b, x, fa, fb, fx;
a = -10, b = 10, x = 0;
fa = pow(a, 3) + 6 * pow(a, 2) + 3 * a - 6;
fb = pow(b, 3) + 6 * pow(b, 2) + 3 * b - 6;
do
{
x = (x > 0) ? b : a;
do
{
x = (a + b) / 2;
fx = pow(x, 3) + 6 * pow(x, 2) + 3 * x - 6;
if (fx * fa < 0)
{
b = x;
fb = pow(b, 3) + 6 * pow(b, 2) + 3 * b - 6;
}
else
{
a = x;
fa = pow(a, 3) + 6 * pow(a, 2) + 3 * a - 6;
}
} while (fabs(fx) >= 1e-5);
printf("x:%f\n", x);
} while (x > -10 && x < 10);
return 0;

}

  • 写回答

1条回答 默认 最新

  • threenewbee 2020-04-26 10:50
    关注

    这里不能用二分法。

    因为函数不是单调增/减的

    但是,你可以换一个思路,首先用二分法求这个函数的导函数的根。再对每个区间再用二分法

    不过最简单高效的是使用牛顿迭代。

    评论

报告相同问题?

悬赏问题

  • ¥15 如何在scanpy上做差异基因和通路富集?
  • ¥20 关于#硬件工程#的问题,请各位专家解答!
  • ¥15 关于#matlab#的问题:期望的系统闭环传递函数为G(s)=wn^2/s^2+2¢wn+wn^2阻尼系数¢=0.707,使系统具有较小的超调量
  • ¥15 FLUENT如何实现在堆积颗粒的上表面加载高斯热源
  • ¥30 截图中的mathematics程序转换成matlab
  • ¥15 动力学代码报错,维度不匹配
  • ¥15 Power query添加列问题
  • ¥50 Kubernetes&Fission&Eleasticsearch
  • ¥15 報錯:Person is not mapped,如何解決?
  • ¥15 c++头文件不能识别CDialog