1. 信号与系统分析作业解析：傅里叶级数实战指南

作为信号与系统课程的核心内容，傅里叶级数分析是理解周期信号频域特性的关键工具。本次作业涵盖了从基础计算到综合应用的完整知识体系，我将结合多年教学经验，带你系统掌握每个知识点的解题要点。

1.1 傅里叶级数分解与频谱绘制

1.1.1 三角形式系数计算要点

对于给定的周期信号，计算傅里叶级数系数需要遵循标准流程：

1. 确定信号周期T：首先观察信号波形，测量或计算其完整周期长度。例如第一小题的三角波周期明显为2π。

2. 计算直流分量a₀：

   a_0 = \frac{1}{T}\int_{T} f(t)dt

   对于对称波形（如奇函数），a₀通常为0。

3. 求解余弦系数aₙ：

   a_n = \frac{2}{T}\int_{T} f(t)\cos(nω_1t)dt

   积分区间选择一个完整周期即可，优先选择对称区间简化计算。

4. 求解正弦系数bₙ：

   b_n = \frac{2}{T}\int_{T} f(t)\sin(nω_1t)dt

实操技巧：遇到分段函数时，将积分拆分为各线性段分别计算。例如三角波可分上升段和下降段积分。

1.1.2 复指数形式转换

得到三角系数后，转换为复指数形式的Fₙ：

F_n = \frac{a_n - jb_n}{2}, \quad F_{-n} = \frac{a_n + jb_n}{2}

特别地，直流分量F₀ = a₀。

典型错误警示：相位计算时注意象限判断，建议同时保留实部和虚部形式避免出错。

1.1.3 单边频谱绘制规范

单边频谱需呈现两个子图：

1. 幅度谱：纵轴|Cₙ| = √(aₙ² + bₙ²)，横轴nω₁
2. 相位谱：纵轴θₙ = arctan(-bₙ/aₙ)，横轴nω₁

绘图要点：

• 使用stem函数绘制离散谱线
• 标注关键频率点的精确值
• 对于衰减谱线，建议用对数坐标展示

1.2 波形与频谱关系深度解析

1.2.1 对称性对谐波的影响

作业中第二题考察了信号对称性与谐波分量的关系：

• 奇函数：仅含正弦项（aₙ=0），且只有奇次谐波
• 偶函数：仅含余弦项（bₙ=0），且只有偶次谐波
• 奇谐函数：f(t) = -f(t+T/2)，仅含奇次谐波

重建波形实操：根据给定1/4周期波形：

1. 奇函数情形：先镜像对称得到半周期，再负向复制
2. 偶函数情形：直接轴对称延拓

1.2.2 占空比与谐波零点

第三题揭示了方波占空比与谐波消失点的关系：

• 占空比τ/T = 1/3时，3k次谐波消失
• 通用规律：当τ/T = m/n（最简分数）时，n/m整数倍谐波为零

推导过程： 方波频谱包络为sinc函数，零点满足：

\frac{nπτ}{T} = kπ \Rightarrow n = k\frac{T}{τ}

当T/τ=3时，n=3k谐波消失。

1.3 傅里叶级数性质证明

1.3.1 时移特性证明

设原信号f(t)的FS系数为Fₙ，则时移后信号g(t)=f(t-t₀)的系数：

G_n = \frac{1}{T}\int_T f(t-t_0)e^{-jnω_1t}dt

令τ=t-t₀，得：

G_n = e^{-jnω_1t_0}\frac{1}{T}\int_T f(τ)e^{-jnω_1τ}dτ = F_n e^{-jnω_1t_0}

物理意义：时移引入线性相位偏移，幅度谱不变。

1.3.2 微分特性证明

对f(t)的FS展开式逐项微分：

f'(t) = \sum_{n=1}^∞ [-nω_1a_n\sin(nω_1t) + nω_1b_n\cos(nω_1t)]

对比f'(t)的FS表达式可得：

a'_n = nω_1b_n, \quad b'_n = -nω_1a_n

应用场景：简化含微分运算的系统分析，如通过输出频谱反推系统特性。

1.4 综合应用题精解

1.4.1 谐波抑制波形设计

设计要求：

1. 无偶次谐波 → 需奇谐对称
2. 无3k次谐波 → 占空比1/3
3. 1/n²衰减 → 信号连续（无跳变）

实现方案：

1. 构造占空比1/3的梯形波（上升/下降沿占周期1/6）
2. 平移T/2并取反，与原信号相加
3. 验证：合成信号满足f(t)=-f(t+T/2)

参数计算： 设幅值A，上升时间τ=T/6：

C_n \propto \frac{\sin(nπ/3)}{nπ/3} \cdot \frac{\sin(nπ/6)}{nπ/6} \sim O(1/n^2)

1.4.2 信号参数确定

根据题目条件：

1. 仅有sin项 → 奇函数
2. 仅有b₁,b₂ → 谐波不超过2次
3. 半周期反对称 → 仅含奇次谐波 ⇒ b₂=0
4. 功率P=(b₁²+b₂²)/2=0.5 ⇒ b₁=1
5. 最终形式：x(t) = sin(πt/2)

排查技巧：遇到矛盾条件时（如既有偶次谐波要求又要求半周期反对称），需重新审视条件解读是否正确。

2. 实验作业：Gibbs现象与Dirichlet条件

2.1 Gibbs现象仿真分析

2.1.1 矩形波合成实验

MATLAB实现要点：

E = 1; T = 1; tao = 0.5; % 设定参数 t = linspace(-0.7,0.7,10000); % 高密度采样 data = ones(size(t))*E*tao/T; % 直流分量 for n = 1:200 data = data + 2*E*tao/T*sinc(n*tao/T)*cos(2*pi*n*t/T); end plot(t,data); % 观察9%的过冲

现象解释：

• 不连续点处出现约9%的振荡
• 增加N仅缩小振荡区域，不减小过冲量
• 本质原因是sinc函数的非一致收敛

2.1.2 锯齿波合成改进

原代码存在幅值误差，修正为：

data = zeros(size(t)); for n = 1:100 data = data + (2*E/pi)*(-1)^(n+1)/n * sin(2*pi*n*t/T); end

收敛对比：

• 锯齿波收敛速度慢于矩形波（1/n vs sinc衰减）
• 连续点处振荡幅度随N增大而减小

2.2 Dirichlet条件验证实验

2.2.1 面积无限信号

典型示例：f(t)=1/t (t∈(0,1])

import numpy as np import matplotlib.pyplot as plt t = np.linspace(0.01,1,1000) f = 1/t # FFT合成演示 N = 20000 F = np.fft.fft(f, 2*N+1) recon = np.zeros_like(f) for n in range(1, N+1): recon += 2/N * np.abs(F[n]) * np.cos(2*np.pi*n*t + np.angle(F[n])) plt.plot(t, recon) # 观察发散现象

2.2.2 无限振荡信号

构造Weierstrass函数：

t = 0:0.001:1; f = zeros(size(t)); for k = 1:100 f = f + cos(3^k*t)/2^k; end

收敛特性：

• 傅里叶系数按指数衰减
• 但时域波形处处不可导

2.2.3 无限间断信号

有理数点脉冲函数：

def f(t): return 1 if t in np.rat else 0 # 数值实现使用稠密采样 t = np.linspace(0,1,10000) f = np.isclose(np.round(t*100), t*100).astype(float)

分析结论：

• 间断点导致高频分量丰富
• 合成信号在连续点收敛，间断点不收敛

3. 傅里叶级数工程应用扩展

3.1 实际信号分析技巧

1. 周期识别：

   • 对实测信号先做自相关分析确定基频
   • 谐波幅度应满足|Cₙ|≤1/n²（连续信号）

2. 频谱泄露处理：

   • 严格整周期采样
   • 使用平顶窗提高幅值精度

3. 噪声环境下的谐波检测：

   [pxx,f] = pwelch(signal, hann(1024),512,1024,fs); harm_idx = find(peaks(pxx) > 3*median(pxx));

3.2 典型问题排查指南

问题1：谐波幅度计算值与理论不符

• 检查信号周期设定是否正确
• 确认积分区间覆盖完整周期
• 验证归一化系数（1/T或2/T）

问题2：相位谱出现π跳变

• 使用atan2函数替代atan
• 检查实部符号：real(Fₙ)>0时θₙ∈(-π/2,π/2)

问题3：合成信号边缘不匹配

• 增加谐波数量N
• 检查Gibbs现象影响范围
• 考虑使用Lanczos sigma因子平滑

通过系统完成这些作业题目，不仅能掌握傅里叶级数的核心计算方法，更能深入理解频域分析的物理意义，为后续学习滤波器设计和信号处理奠定坚实基础。在实际操作中遇到具体问题时，建议先用简单信号（如方波、三角波）验证算法正确性，��逐步过渡到复杂信号分析。