解一元三次方程

浅谈一元三次方程

一元三次方程的标准形式（即所有一元三次方程经整理都能得到的形式）是ax³+bx²+cx+d=0（a，b，c，d为常数，x为未知数，且a≠0）。一元三次方程的公式解法有卡尔丹公式法与盛金公式法。两种公式法都可以解标准型的一元三次方程。由于用卡尔丹公式解题存在复杂性，相比之下，盛金公式解题更为直观，效率更高。

盛金公式法：

三次方程应用广泛。用根号解一元三次方程，虽然有著名的卡尔丹公式，并有相应的判别法，但使用卡尔丹公式解题比较复杂，缺乏直观性。范盛金推导出一套直接用a、b、c、d表达的较简明形式的一元三次方程的一般式新求根公式——盛金公式，并建立了新判别法——盛金判别法。

用Python代码实现一元三次方程求解

第一步：用户输入

print("整理格式：ax³+bx²+cx+d=0（a，b，c，d为常数，x为未知数，且a≠0）")
a = float(input("请输入三次项系数："))
b = float(input("请输入二次项系数："))
c = float(input("请输入一次项系数："))
d = float(input("请输入常数项："))

整理格式后依次输入系数，用 float 强转。

第二步：定义重根判别式

A = b**2 - 3*a*c
B = b*c - 9*a*d
C = c**2 - 3*b*d

第三步：定义总判别式

D = B**2 - 4*A*C

盛金判别法

当A=B=0时，方程有一个三重实根。
当Δ=B²-4AC>0时，方程有一个实根和一对共轭虚根。
当Δ=B²-4AC=0时，方程有三个实根，其中有一个二重根。
当Δ=B²-4AC<0时，方程有三个不相等的实根。

第四步：根据总判别式分情况讨论

if A == B == 0:
    #盛金公式1
    print("方程有三个相同的实数根：",-1*c/b)
elif D > 0:
    # 盛金公式2
    Y1 = A * b + 3 * a * (1 / 2 * (-B + (D ** (1 / 2))))
    Y2 = A * b + 3 * a * (1 / 2 * (-B - (D ** (1 / 2))))
    print("X1=", 1 / (a * 3) * (-b - Y1 ** (1 / 3) - Y2 ** (1 / 3)))
    print("X2=", 1 / (a * 6) * (
    -2 * b + Y1 ** (1 / 3) + Y2 ** (1 / 3) - 3 ** (1 / 2) * (Y1 ** (1 / 3) - Y2 ** (1 / 3)) * ((-1) ** (1 / 2))))
    print("X3=", 1 / (a * 6) * (
    -2 * b + Y1 ** (1 / 3) + Y2 ** (1 / 3) + 3 ** (1 / 2) * (Y1 ** (1 / 3) - Y2 ** (1 / 3)) * ((-1) ** (1 / 2))))
elif D == 0:
    # 盛金公式3
    K = B / A
    print("X1=",-b/a + K)
    print("X2=X3=",-K/2)
elif D < 0:
    # 盛金公式4
    import math
    T = (2*A*b - 3 * a * B) / (2 * (A ** (3 / 2)))
    p = math.acos(T)
    g3 = math.cos(p / 3) - math.sin(p / 3) * (3 ** (1 / 2))
    g2 = math.cos(p / 3) + math.sin(p / 3) * (3 ** (1 / 2))
    print("X1=",1 / (3 * a) * (-b - 2 * (A ** (1/2) * math.cos(p / 3))))
    print("X2=",1 / (3 * a) * (-b + (A ** (1/2) * g2 )))
    print("X3=", 1 / (3 * a) * (-b + (A ** (1 / 2) * g3)))

代码示例：

博主推荐

n阶行列式求解
解一元二次方程
3秒解一元二次不等式