把用户的数学/物理概念转化为 Manim CE 动画视频（零…

# 傅里叶级数逼近方波动画制作指南

一、数学原理（课堂讲解用）

• 方波函数：f(x) = 1（0<x<π），f(x) = -1（π<x<2π），周期为 2π
• 傅里叶展开只含奇次谐波：

f(x) ≈ (4/π) × Σ [sin((2k-1)x) / (2k-1)]，k = 1, 2, 3, …

• 项数 N 越大，逼近越好，但在间断点附近会出现吉布斯现象（约9%的过冲），这是很好的课堂亮点

二、Python 实现代码

```python

import numpy as np

import matplotlib.pyplot as plt

from matplotlib.animation import FuncAnimation

x = np.linspace(0, 2*np.pi, 1000)

def partial_sum(n): # n 为谐波个数

s = np.zeros_like(x)

for k in range(1, n+1):

s += np.sin((2*k-1)*x) / (2*k-1)

return 4/np.pi * s

fig, ax = plt.subplots(figsize=(12.8, 7.2), dpi=100) # 1280x720，1080p可改dpi

ax.set_xlim(0, 2*np.pi); ax.set_ylim(-1.6, 1.6)

ax.axhline(0, color='gray', lw=0.5)

sq = np.sign(np.sin(x)) # 目标方波

ax.plot(x, sq, 'k--', label='方波')

line, = ax.plot([], [], 'r-', lw=2)

text = ax.text(0.1, 1.3, '', fontsize=14)

def update(frame):

n = [1, 3, 5, 10, 50][frame]

line.set_data(x, partial_sum(n))

text.set_text(f'叠加前 {n} 项谐波')

return line, text

ani = FuncAnimation(fig, update, frames=5, interval=1500, repeat=True)

ani.save('fourier_square.mp4', fps=1, dpi=150) # dpi=150时约1080p

```

三、导出 1080p 视频要点

• 分辨率：figsize=(16, 9) 且 dpi=120，即 1920×1080
• 依赖：需安装 ffmpeg（pip install matplotlib 后配置 ffmpeg 路径），否则改用 ani.save('xxx.gif')
• 每帧停留 1.5–2 秒便于讲解；也可额外做逐项累加的连续动画版本

四、课堂展示建议

• 对比图中用不同颜色保留前一帧曲线，直观展示收敛过程
• 重点讲解吉布斯现象：50 项时边缘过冲依然存在，只是位置向间断点靠拢
• 可补充提问：“为什么方波只有 sin 奇次项？”（奇函数对称性 + 半波对称）

五、备选工具

• 不写代码：GeoGebra / Desmos 手动输入级数表达式也可动态展示
• 进阶版：用 manim 生成数学动画，观感更专业