### Python 实现多项式零点式展开并提取各次项系数 在 Python 中可以利用 `numpy` 和 `sympy` 这两个强大的科学计算库来完成多项式的零点式展开以及提取特定次幂的系数的任务。以下是具体的实现方法： #### 使用 SymPy 展开多项式并提取系数 SymPy 是一个用于符号数学的强大库，它可以直接表示和操作多项式表达式。 ```python from sympy import symbols, expand def polynomial_from_zeros(zeros): x = symbols('x') # 构建多项式 (x - zero1)(x - zero2)...(x - zeron) polynomial = 1 for zero in zeros: polynomial *= (x - zero) expanded_poly = expand(polynomial) # 展开多项式 return expanded_poly def extract_coefficients(expanded_poly): coeffs_dict = {} degree = expanded_poly.as_poly().degree() # 获取最高次幂 for i in range(degree + 1): # 遍历从高到低的所有次幂 coeff = expanded_poly.coeff('x', i) # 提取指定次幂的系数 if i == 2: coeffs_dict['x_squared'] = coeff elif i == 1: coeffs_dict['x_linear'] = coeff elif i == 0: coeffs_dict['constant_term'] = coeff return coeffs_dict # 示例：给定零点 [-1, 2, 3] zeros = [-1, 2, 3] expanded_poly = polynomial_from_zeros(zeros) coeffs = extract_coefficients(expanded_poly) print(f"Expanded Polynomial: {expanded_poly}") print(f"Coefficient of x^2: {coeffs.get('x_squared')}") print(f"Coefficient of x: {coeffs.get('x_linear')}") print(f"Constant Term: {coeffs.get('constant_term')}") ``` #### 输出解释 上述代码实现了以下功能： 1. **构建多项式**：通过输入一组零点，程序会自动构造对应的多项式 `(x - zero1)(x - zero2)...(x - zeron)` 并将其展开[^1]。 2. **提取系数**：对于展开后的多项式，分别提取 \(x^2\) 的系数、\(x\) 的线性系数以及常数项[^1]。 运行以上代码的结果如下（假设零点为 \([-1, 2, 3]\)）： ``` Expanded Polynomial: x**3 - 4*x**2 + x + 6 Coefficient of x^2: -4 Coefficient of x: 1 Constant Term: 6 ``` --- #### 使用 NumPy 多项式模块 NumPy 提供了一个专门处理多项式的子模块 `numpy.polynomial.Polynomial`，也可以用来解决这个问题。 ```python import numpy as np from numpy.polynomial.polynomial import Polynomial def polynomial_from_zeros_numpy(zeros): roots = list(map(lambda r: (-r), zeros)) # 转换根的形式以便于使用Polynomial.fromroots() poly = Polynomial.fromroots(roots).convert(domain=[-1, 1]) # 创建多项式对象 return poly.coef[::-1] # 返回按降序排列的系数列表 # 示例：给定零点 [-1, 2, 3] zeros = [-1, 2, 3] coef_list = polynomial_from_zeros_numpy(zeros) # 手动解析系数 if len(coef_list) >= 3: coef_x_squared = coef_list[2] else: coef_x_squared = 0 if len(coef_list) >= 2: coef_x_linear = coef_list[1] else: coef_x_linear = 0 if len(coef_list) >= 1: constant_term = coef_list[0] else: constant_term = 0 print(f"Coefficient of x^2: {coef_x_squared}") print(f"Coefficient of x: {coef_x_linear}") print(f"Constant Term: {constant_term}") ``` #### 输出解释 此代码片段同样能够生成多项式并提取所需的系数。注意，在 NumPy 中，多项式的系数是以升序存储的，因此需要反转数组以匹配标准形式[^1]。 --- ### 总结 两种方法各有优劣： - 如果更倾向于符号运算或者希望获得完整的多项式表达式，推荐使用 SymPy； - 若仅需数值结果而无需显示具体表达式，则 NumPy 更加高效。