如何使用Mathematica求解一般方程及任意次数多项式方程

一、求解给定的具体方程 #

你给出的方程x^2 + 2x^2 + 2x + 1 == 0可直接用Mathematica处理，无需手动合并同类项：

1. 直接调用Solve函数 #

在Mathematica中输入以下代码并运行：

Solve[x^2 + 2x^2 + 2x + 1 == 0, x]

系统会自动合并同类项得到3x² + 2x + 1 == 0，并输出精确复数解：

{{x -> (-1 - I Sqrt[2])/3}, {x -> (-1 + I Sqrt[2])/3}}

2. 验证解的正确性 #

可以用两种方式验证：

• 转换为数值解查看近似值：

  N[Solve[x^2 + 2x^2 + 2x + 1 == 0, x]]
  

• 将解代入原方程验证是否等于0：

  (x^2 + 2x^2 + 2x + 1) /. x -> (-1 - I Sqrt[2])/3
  

运行后会得到0，确认解的正确性。

二、求解任意次数的一般多项式方程 #

Mathematica针对不同次数的多项式，有不同的处理方案：

1. 低次多项式（5次及以下） #

对于次数≤5的多项式，Solve函数可直接给出精确代数解。比如三次方程：

Solve[x^3 - 6x^2 + 11x - 6 == 0, x]

会输出整数解{{x -> 1}, {x -> 2}, {x -> 3}}；无有理根的低次方程则会给出根式形式的解。

2. 高次多项式（6次及以上） #

根据阿贝尔-鲁菲尼定理，一般6次及以上多项式没有根式解，可通过以下方式处理：

• 求数值解：用NSolve直接获取高精度数值解，示例：

  NSolve[x^6 - 5x^5 + 3x^4 + 10x^3 - 12x^2 - 8x + 16 == 0, x]
  

• 提取有理根：先用Factor函数因式分解多项式，拆分出有理根对应的因子后分别求解：

  Factor[x^6 - 5x^5 + 3x^4 + 10x^3 - 12x^2 - 8x + 16]
  

• 求区间内的实根：如果只需要某区间内的实根，用FindRoot并提供初始猜测值：

  FindRoot[x^6 - 5x^5 + 3x^4 + 10x^3 - 12x^2 - 8x + 16 == 0, {x, 1}]
  

3. 含参数的一般多项式 #

如果多项式包含未知参数（如a x^2 + b x + c == 0），Solve会给出参数化的通解：

Solve[a x^2 + b x + c == 0, x]

输出为{{x -> (-b - Sqrt[b^2 - 4 a c])/(2 a)}, {x -> (-b + Sqrt[b^2 - 4 a c])/(2 a)}}（默认假设a≠0）。

内容的提问来源于stack exchange，提问作者lokesh