SymPy的fourier_series是否存在问题？运行官方示例计算0<x<1上f=x的傅里叶级数结果不符

参数与官方参考示例的计算区间不匹配 #

运行代码时给fourier_series传入的区间参数为(x, 0, 1)，但SymPy官方文档示例图对应的计算区间为(x, -pi, pi)，两者计算得到的傅里叶级数存在本质差异：

• 区间(-pi, pi)对应周期长度为2π，计算得到的是经典斜变锯齿波的傅里叶级数，10项截断后可在连续区间内较好贴合原函数走势
• 区间(0,1)对应周期长度为1，计算得到的是带间断点的锯齿波（在x=1位置函数值从1直接跳变回0），受吉布斯现象影响，即便增加截断项数，在间断点附近也会存在固定比例的幅值超调，视觉上和“完全拟合原函数”的预期存在明显差距

旧版本的计算逻辑bug #

SymPy 1.7及更早的旧版本中，fourier_series对非对称、长度非2π的周期区间存在基函数归一化的逻辑错误：

• 旧版本不会根据实际区间长度计算对应的基函数角频率，默认按2π周期推导三角项，最终生成的截断级数和理论傅里叶级数存在量级偏差，直接表现为0<x<1区间内完全无法匹配f=x的走势

绘图配置的视觉干扰 #

默认绘图参数无法直观展示拟合效果：

• SymPy默认绘图会自动适配所有曲线的坐标范围，当截断级数存在吉布斯超调的极值点时，纵轴范围会被过度拉大，视觉上会显得两条曲线偏差极大
• 未对曲线样式做区分时，很容易混淆原函数和截断级数的对应关系，误判为拟合失败

复现官方示例的操作：将傅里叶级数的区间参数修改为(x, -pi, pi)，重新运行代码即可得到和文档示例一致的效果。如果确实需要计算0<x<1区间上周期为1的f=x的傅里叶级数，需先将SymPy升级至1.10及以上稳定版本，再重新执行计算流程。

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