大学项目：如何将BSpline/PPoly对象转为多项式数组写入文件？

核心思路 #

用scipy.interpolate.PPoly.from_spline()把BSpline转为PPoly对象后，直接通过PPoly的属性获取节点和分段系数，就能得到可写入文件的NumPy数组。

完整代码示例 #

import numpy as np
from scipy.interpolate import make_interp_spline, PPoly

# 1. 从txt读取数据（实际替换为你的读取逻辑）
x_data, y_data = np.loadtxt("input_data.txt", unpack=True)

# 2. 生成三次样条的BSpline对象
bspline = make_interp_spline(x_data, y_data, k=3)

# 3. 转换为PPoly分段多项式对象
ppoly = PPoly.from_spline(bspline)

# 4. 获取节点和系数数组
# ppoly.x: 样条节点数组，长度为n
# ppoly.coeffs: 系数数组，形状为(4, n-1)，每一列对应一个区间的[三次项系数, 二次项, 一次项, 常数项]
nodes = ppoly.x
coeffs = ppoly.coeffs

# 5. 将结果写入txt文件
with open("spline_coeffs.txt", "w") as f:
    # 写入节点
    f.write("Spline Nodes:\n")
    np.savetxt(f, nodes.reshape(1, -1), fmt="%.6f")
    
    # 写入分段系数，每个区间一行，按三次、二次、一次、常数项顺序
    f.write("\nSegment Coefficients (c3, c2, c1, c0):\n")
    np.savetxt(f, coeffs.T, fmt="%.6f")

关键细节说明 #

• 三次样条对应PPoly的系数数组是4行（覆盖0到3阶项），每一列对应一个区间的多项式系数，转置后(coeffs.T)每一行对应一个区间，更符合写入文件后的阅读逻辑。
• ppoly.coeffs本身就是标准NumPy数组，可直接用于保存、计算等操作。
• 若要验证转换正确性，可对比bspline(x)和ppoly(x)的计算结果，两者数值完全一致。

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