一元三次方程的实数根

一元三次方程的实数根

一元三次方程的一般形式

一元三次方程的一般形式为

$$ax^3+bx^2+cx+d=0,\quad a\neq 0.$$
令 $x=y-\dfrac{b}{3a}$，则原式变成
$$\begin{equation*} \begin{aligned} &a\left(y-\dfrac{b}{3a}\right)^3 +b\left(y-\dfrac{b}{3a}\right)^2 +c\left(y-\dfrac{b}{3a}\right) +d=0, \\ &a\left(y^3-\dfrac{by^2}{a}+\dfrac{b^2y}{3a^2}-\dfrac{b^3}{27a^3}\right) +b\left(y^2-\dfrac{2by}{3a}+\dfrac{b^2}{9a^2}\right) +c\left(y-\dfrac{b}{3a}\right) +d =0, \\ &ay^3-by^2+\dfrac{b^2}{3a}y-\dfrac{b^3}{27a^2} +by^2-\dfrac{2b^2}{3a}y+\dfrac{b^3}{9a^2} +cy-\dfrac{bc}{3a} +d=0, \\ &ay^3+\left(c-\dfrac{b^2}{3a}\right)y+\left(d+\dfrac{2b^3}{27a^2}-\dfrac{bc}{3a}\right) =0, \\ &y^3+\left(\dfrac{c}{a}-\dfrac{b^2}{3a^2}\right)y+\left(\dfrac{d}{a}+\dfrac{2b^3}{27a^3}-\dfrac{bc}{3a^2}\right) =0. \end{aligned} \end{equation*}$$
如此一来二次项就不见了，化成 $y^3+3py+2q=0$，其中 $p=\dfrac{c}{3a}-\dfrac{b^2}{9a^2}$， $q=\dfrac{d}{2a}+\dfrac{b^3}{27a^3}-\dfrac{bc}{6a^2}$。

方程$y^3+3py+2q=0$的实数根

令$\Delta=q^2+p^3$，则

$(i)$ 当$\Delta>0$时，方程只有一个实根。

$$y=\sqrt[3]{-q+\sqrt{\Delta}}+\sqrt[3]{-q-\sqrt{\Delta}}.$$
$(ii)$ 当 $\Delta=0$时，方程有三个实根，其中至少有两个相等的实根。
$$\begin{equation*} \begin{aligned} &y_1=-2\sqrt[3]{q}, \\ &y_2=y_3=\sqrt[3]{q}. \end{aligned} \end{equation*}$$
$(iii)$ 当 $\Delta<0$时，方程有三个不等实根。

令$\alpha=\dfrac{1}{3}\arccos\dfrac{-q\sqrt{-p}}{p^2}$，则

$$\begin{equation*} \begin{aligned} &y_1=2\sqrt{-p}\cos\alpha, \\ &y_2=2\sqrt{-p}\cos\left(\alpha+120^\circ\right), \\ &y_3=2\sqrt{-p}\cos\left(\alpha+240^\circ\right). \end{aligned} \end{equation*}$$

C++代码实现

// π
const double PIE = 3.1415926535897932384626433832795;
// 求一个实数的立方根
static double SolveCubicRoot(double value);
// 求一元一次方程的实根：ax+b=0
static std::vector<double> SolveLinearEquation(double a, double b);
// 求一元二次方程的实根：ax^2+bx+c=0
static std::vector<double> SolveQuadraticEquation(double a, double b, double c);

// 求一元三次方程的实根：ax^3+bx^2+cx+d=0
std::vector<double> SolveCubicEquation(double a, double b, double c, double d)
{
    // 判断三次项系数是否为零
    if (a == 0)
    {
        return SolveQuadraticEquation(b, c, d);
    }

    std::vector<double> root;

    // 系数
    double p = (c / a - b * b / (3 * a * a)) / 3;
    double q = (d / a + 2 * b * b * b / (27 * a * a * a) - b * c / (3 * a * a)) / 2;
    double diff = -b / (3 * a);

    // 判别式
    double delta = p * p * p + q * q;
    if (delta > 0)
    {
        // 方程只有一个实根
        double sqrtDelta = sqrt(delta);
        root.push_back(SolveCubicRoot(-q + sqrtDelta) + SolveCubicRoot(-q - sqrtDelta));
    }
    else if (delta < 0)
    {
        // 方程有三个不等的实根
        double angle = acos(-q * sqrt(-p) / (p * p)) / 3;
        root.push_back(2.0 * sqrt(-p) * cos(angle));
        root.push_back(2.0 * sqrt(-p) * cos(angle + 2 * PIE / 3));
        root.push_back(2.0 * sqrt(-p) * cos(angle + 4 * PIE / 3));
    }
    else
    {
        // 方程有三个实根，其中至少有两个相等的实根
        if (q == 0)
        {
            root.push_back(0);
        }
        else
        {
            root.push_back(SolveCubicRoot(q));
            root.push_back(-2 * SolveCubicRoot(q));
        }
    }

    // 将结果加上余量
    for (int i = 0, maxSize = root.size(); i < maxSize; ++i)
    {
        root[i] += diff;
    }

    return root;
}

double SolveCubicRoot(double value)
{
    return value > 0 ? pow(value, 1.0 / 3) : -pow(-value, 1.0 / 3);
}

std::vector<double> SolveLinearEquation(double a, double b)
{
    std::vector<double> root;

    // 判断一次项系数是否为零
    if (a != 0)
    {
        root.push_back(-b / a);
    }

    return root;
}

std::vector<double> SolveQuadraticEquation(double a, double b, double c)
{
    // 判断二次项系数是否为零
    if (a == 0)
    {
        return SolveLinearEquation(b, c);
    }

    std::vector<double> root;

    // 计算判别式
    double delta = b * b - 4 * a * c;
    if (delta > 0)
    {
        // 方程有两个不等的实根
        root.push_back((-b + sqrt(delta)) / (2 * a));
        root.push_back((-b - sqrt(delta)) / (2 * a));
    }
    else if (delta == 0)
    {
        // 方程有两个相等的实根
        root.push_back(-b / (2 * a));
    }

    return root;
}