跳到内容

课程笔记

DSP 期末考试速查

结合 2022–2025 四套真题整理的开卷速查:高频公式、做题流程、易错检查和题型-章节映射。

DSP 期末考试速查

本页是把 11 章笔记按考试使用频率重新排列的精简版。考场翻到这一页,先确认题型属于哪一类,再按流程写答案。

1. 题型 → 章节速查表

题干关键词章节核心操作
real sequence DFT、conjugate symmetry、odd/even DFTch5X[k]=X[kN]X[k]=X^*[\langle-k\rangle_N]
DIT-FFT、bit reversal、butterfly、复乘复加ch5+ch11N/2log2NN/2\cdot\log_2 N 复乘,DIT 输入比特反转
DFT length change、zero padding、repeat、upsamplech5补零→更密采样DTFT;重复→奇数bin为零;插零→频谱压缩+镜像
sampling、aliasing、reconstruction、(1)n(-1)^nch3ω=ΩT=2πf/Fs\omega=\Omega T=2\pi f/F_s(1)n=ejπn(-1)^n=e^{j\pi n} 平移 π\pi
linear、time-invariant、causal、stablech4按定义证明
cascade、parallel、overall responsech4级联 H=H1H2H=H_1H_2;并联 H=H1+H2H=H_1+H_2
input/output spectrum、system identificationch4H=Y/XH=Y/X 只在 X0X\ne0 频率可确定;δ[n]\delta[n] 最佳
H(z)H(z)、ROC、pole-zero、difference equationch6差分方程 H(z)h[n]\leftrightarrow H(z) \leftrightarrow h[n];因果 ROC 在最外极点外
Direct Form II、canonical structurech8分母反馈 ak-a_k、分子前馈 bkb_k 共用延时链
linear phase FIR Type I–IVch7对称/反对称 + 长度奇偶
minimum phase、maximum phase、allpass、equalizerch7单位圆外零点倒共轭 1/z01/z_0^* 进单位圆
Butterworth、bilinear transform、prewarpingch9Ω=2Ttan(ω/2)\Omega=\frac{2}{T}\tan(\omega/2);阶数向上取整
window method、Hamming/Hann/Blackmanch10指标→ωc\omega_c→选窗→算阶数→h[n]=hd[nα]w[n]h[n]=h_d[n-\alpha]w[n]
fewest multipliers、linear phase structurech8+ch10对称先加再乘

2. 频率换算三件套

ω=ΩT=2πfFs\boxed{\omega=\Omega T=\frac{2\pi f}{F_s}}
已知公式
Hz ff数字频率 ω\omegaω=2πf/Fs\omega=2\pi f/F_s
数字频率 ω\omegaHz fff=ωFs/(2π)f=\omega F_s/(2\pi)
模拟角频率 Ω\Omega数字频率 ω\omegaω=ΩT\omega=\Omega T
DFT bin kkHz fkf_kfk=kFs/Nf_k=kF_s/NkN/2k\le N/2);fk=(kN)Fs/Nf_k=(k-N)F_s/Nk>N/2k>N/2
f=Fs/2f=F_s/2ω\omegaω=π\omega=\pi(Nyquist)

⚠️ 考场第一件事: 看题目给的频率单位。给 Hz 先换 ω\omega,给 kHz 先除 1000。

3. DFT / FFT 核心公式

3.1 DFT 定义

X[k]=n=0N1x[n]WNkn,WN=ej2π/NX[k]=\sum_{n=0}^{N-1}x[n]W_N^{kn},\quad W_N=e^{-j2\pi/N} x[n]=1Nk=0N1X[k]WNknx[n]=\frac{1}{N}\sum_{k=0}^{N-1}X[k]W_N^{-kn}

⚠️ IDFT 有 1/N1/N。考试常见丢分点。

3.2 快速检查公式

检查项公式
X[0]X[0](直流)n=0N1x[n]\sum_{n=0}^{N-1}x[n]
x[0]x[0]1Nk=0N1X[k]\frac{1}{N}\sum_{k=0}^{N-1}X[k]
能量(Parseval)x[n]2=1NX[k]2\sum\|x[n]\|^2=\frac{1}{N}\sum\|X[k]\|^2

3.3 实序列 DFT 对称

X[k]=X[kN]X[k]=X^*[\langle -k\rangle_N]
  • X[0]X[0] 一定实数
  • NN 偶数时 X[N/2]X[N/2] 一定实数
  • 实部偶对称、虚部奇对称、幅度偶对称、相位奇对称

3.4 DFT 性质速查

时域频域
x[nn0N]x[\langle n-n_0\rangle_N]WNkn0X[k]W_N^{kn_0}X[k]
WNk0nx[n]W_N^{-k_0n}x[n]X[kk0N]X[\langle k-k_0\rangle_N]
(1)nx[n](-1)^n x[n]NN 偶)X[kN/2N]X[\langle k-N/2\rangle_N]
Y[k]=(1)kX[k]Y[k]=(-1)^kX[k]y[n]=x[nN/2N]y[n]=x[\langle n-N/2\rangle_N]

3.5 FFT 复杂度

  • 级数:log2N\log_2 N
  • 复乘:N2log2N\frac{N}{2}\log_2 N
  • 复加:Nlog2NN\log_2 N
  • DIT:输入 bit-reversed,输出 normal
  • DIF:输入 normal,输出 bit-reversed

8 点 DIT 输入顺序:0,4,2,6,1,5,3,70,4,2,6,1,5,3,7

3.6 DFT 长度变化

操作频域效果
末尾补零到 NN'DTFT 更密采样,不增加真实分辨率
时域重复成 2N2NY[2k]=2X[k]Y[2k]=2X[k]Y[2k+1]=0Y[2k+1]=0
时域插零(y[2n]=x[n]y[2n]=x[n]Y[k]=X[kmodN]Y[k]=X[k\bmod N](频谱压缩+镜像)

4. 采样与恢复

4.1 采样频谱

Gp(jΩ)=1Tk=Ga(j(ΩkΩs))G_p(j\Omega)=\frac{1}{T}\sum_{k=-\infty}^{\infty}G_a(j(\Omega-k\Omega_s))
  • 幅度有 1/T1/T 因子
  • 副本间隔 Ωs=2π/T=2πFs\Omega_s=2\pi/T=2\pi F_s

4.2 Nyquist 条件

Fs>2fmΩs>2ΩmF_s>2f_m\quad\Longleftrightarrow\quad\Omega_s>2\Omega_m

4.3 混叠折回

频率 ff 采样后等效为距最近 kFskF_s 的差:

falias=fkFsf_{alias}=|f-kF_s|

数字频率 >π>\pi 折回:cos(ωn)=cos((ω2π)n)\cos(\omega n)=\cos((\omega-2\pi)n)ω>π\omega>\pi 时)

4.4 重构滤波器

通带增益 =T=T(抵消采样时 1/T1/T

4.5 (1)n(-1)^n 频谱搬移

(1)n=ejπn,(1)nx[n]X(ej(ωπ))(-1)^n=e^{j\pi n},\quad (-1)^n x[n]\leftrightarrow X(e^{j(\omega-\pi)})

频谱复用 y[n]=x1[n]+(1)nx2[n]y[n]=x_1[n]+(-1)^n x_2[n]:低频 x1x_1 + 高频搬移的 x2x_2,可用低通/高通分离。

5. Z 变换与系统函数

5.1 常用变换对

序列Z 变换ROC
anu[n]a^n u[n]11az1\frac{1}{1-az^{-1}}z>a\|z\|>\|a\|
anu[n1]-a^n u[-n-1]11az1\frac{1}{1-az^{-1}}z<a\|z\|<\|a\|
δ[n]\delta[n]11全平面

5.2 ROC 判断

序列类型ROC
右边序列(因果)最外极点外
左边序列(反因果)最内极点内
双边序列极点之间圆环

5.3 因果稳定

  • 因果:ROC 在最外极点外
  • 稳定:ROC 包含单位圆
  • 因果且稳定有理系统:所有极点在单位圆内

5.4 差分方程 → H(z)H(z)

aky[nk]=bkx[nk]H(z)=bkzkakzk\sum a_k y[n-k]=\sum b_k x[n-k]\quad\Longleftrightarrow\quad H(z)=\frac{\sum b_k z^{-k}}{\sum a_k z^{-k}}

⚠️ 反馈项符号:分母 1+a1z11+a_1 z^{-1} 对应 y[n]=a1y[n1]+y[n]=-a_1 y[n-1]+\cdots

5.5 频率响应

H(ejω)=H(z)z=ejωH(e^{j\omega})=H(z)\big|_{z=e^{j\omega}}

前提是单位圆在 ROC 中。

6. 系统性质证明模板

线性

  1. x1y1x_1\to y_1x2y2x_2\to y_2
  2. x=ax1+bx2x=ax_1+bx_2
  3. 算输出,看是否等于 ay1+by2ay_1+by_2

时不变

  1. 原输入 x1y1x_1\to y_1
  2. 延迟输入 x2[n]=x1[nn0]x_2[n]=x_1[n-n_0]
  3. 新输出 y2[n]y_2[n] vs 原输出延迟 y1[nn0]y_1[n-n_0]

因果

y[n0]y[n_0] 是否依赖 x[n]x[n]n>n0n>n_0

稳定(LTI)

nh[n]<\sum_n |h[n]|<\infty

7. 线性相位 FIR 四类型

类型对称/反对称长度强制零点适合
Type I对称奇数低通(首选)
Type II对称偶数ω=π\omega=\pi低通(可行);不做普通高通
Type III反对称奇数ω=0,π\omega=0,\piHilbert、微分器
Type IV反对称偶数ω=0\omega=0Hilbert、微分器

最小相位构造

零点 z0z_0 在单位圆外 → 替换为 1/z01/z_0^*(倒共轭到单位圆内) → 原系统 = 最小相位 × 全通

均衡器

原系统有单位圆外零点 → 逆系统有单位圆外极点 → 不存在稳定因果逆

8. IIR 滤波器设计(Butterworth + 双线性)

考场完整流程

Hzω=2πf/Fs数字频率Ω=2Ttan(ω/2)预畸变ϵ2,A2阶数BLTH(z)\boxed{ \text{Hz}\xrightarrow{\omega=2\pi f/F_s}\text{数字频率}\xrightarrow{\Omega=\frac{2}{T}\tan(\omega/2)}\text{预畸变}\xrightarrow{\epsilon^2,A^2}\text{阶数}\xrightarrow{\text{BLT}}H(z) }

步骤 1: 统一数字频率 ωp=2πfp/Fs\omega_p=2\pi f_p/F_sωs=2πfs/Fs\omega_s=2\pi f_s/F_s

步骤 2: 预畸变 Ωp=2Ttan(ωp/2)\Omega_p=\frac{2}{T}\tan(\omega_p/2)Ωs=2Ttan(ωs/2)\Omega_s=\frac{2}{T}\tan(\omega_s/2)

步骤 3: dB → 线性

ϵ2=10αp/101,A2=10αs/10\epsilon^2=10^{\alpha_p/10}-1,\quad A^2=10^{\alpha_s/10}

步骤 4: Butterworth 阶数(向上取整

Nlog10[(A21)/ϵ2]2log10(Ωs/Ωp)N\ge\frac{\log_{10}[(A^2-1)/\epsilon^2]}{2\log_{10}(\Omega_s/\Omega_p)}

步骤 5: 截止频率 Ωc=Ωp/ϵ1/N\Omega_c=\Omega_p/\epsilon^{1/N}

步骤 6: 双线性变换 s=2T1z11+z1s=\frac{2}{T}\frac{1-z^{-1}}{1+z^{-1}}

⚠️ 常见错:不做预畸变、阶数四舍五入、dB 用 10log10\log 代替 20log20\log

9. FIR 滤波器设计(窗函数法)

考场完整流程

HzωΔωωc选窗M,Nhd[m]h[n]=hd[nα]w[n]\boxed{ \text{Hz}\to\omega\to\Delta\omega\to\omega_c\to\text{选窗}\to M,N\to h_d[m]\to h[n]=h_d[n-\alpha]w[n] }

步骤 1: ωp=2πfp/Fs\omega_p=2\pi f_p/F_sωs=2πfs/Fs\omega_s=2\pi f_s/F_s

步骤 2: Δω=ωsωp\Delta\omega=\omega_s-\omega_pωc=(ωp+ωs)/2\omega_c=(\omega_p+\omega_s)/2

步骤 3: 选窗(阻带衰减 A=20log10δA=-20\log_{10}\deltaδ=min(δp,δs)\delta=\min(\delta_p,\delta_s)

衰减过渡带宽
Rectangular~21 dB4π/N4\pi/N
Hann~44 dB8π/N8\pi/N
Hamming~53 dB8π/N8\pi/N
Blackman~74 dB12π/N12\pi/N

步骤 4: 估算 NCπ/ΔωN\approx C\pi/\Delta\omega,Type I 取奇数 NNM=N1M=N-1α=M/2\alpha=M/2

步骤 5: 理想低通 hd[m]=sin(ωcm)πmh_d[m]=\frac{\sin(\omega_c m)}{\pi m}m0m\ne0),hd[0]=ωc/πh_d[0]=\omega_c/\pi

步骤 6: 加窗 h[n]=hd[nα]w[n]h[n]=h_d[n-\alpha]w[n]

步骤 7: 最省乘法结构:对称先加再乘,乘法器数 N/2\lceil N/2\rceil

高通 / 带通 / 带阻

  • 高通:hHP[m]=δ[m]hLP[m]h_{HP}[m]=\delta[m]-h_{LP}[m]
  • 带通:hBP[m]=sin(ω2m)sin(ω1m)πmh_{BP}[m]=\frac{\sin(\omega_2 m)-\sin(\omega_1 m)}{\pi m}
  • 带阻:hBS[m]=δ[m]hBP[m]h_{BS}[m]=\delta[m]-h_{BP}[m]

10. 滤波器结构

基本构件

构件符号作用
延时器z1z^{-1}存储一个采样点
乘法器系数 aka_k, bkb_k乘以常数
加法器++求和

DF-I vs DF-II

特性DF-IDF-II
延时链两条(前馈 + 反馈)一条(共用 w[n]w[n]
延时器数M+NM+Nmax(M,N)\max(M,N)
别名canonical structure
优先级概念清晰考试首选

DF-II 核心方程

w[n]=x[n]k=1Nakw[nk]w[n]=x[n]-\sum_{k=1}^{N}a_k w[n-k] y[n]=k=0Mbkw[nk]y[n]=\sum_{k=0}^{M}b_k w[n-k]

延时器数:max(M,N)\max(M,N)(canonical)

⚠️ 反馈来自 w[nk]w[n-k],不是 y[nk]y[n-k]

级联 vs 并联

特性级联 (Cascade)并联 (Parallel)
系统函数H=HiH=\prod H_i(相乘)H=C+HiH=C+\sum H_i(相加)
实现方式前一节输出接后一节输入所有支路共用输入,输出相加
展开方法因式分解部分分式展开
优点量化误差可控各支路独立,易于并行

级联结构步骤:

  1. 因式分解 H(z)=GHi(z)H(z)=G\prod H_i(z)
  2. 每个 Hi(z)H_i(z) 用一阶或二阶节实现
  3. 复共轭极/零点放同一二阶节

并联结构步骤:

  1. 部分分式展开 H(z)=C+Hi(z)H(z)=C+\sum H_i(z)
  2. 不要漏掉常数直通项 CC
  3. 每个支路独立实现

FIR 对称省乘法

h[n]=h[Mn]h[n]{x[kn]+x[kM+n]}h[n]=h[M-n]\Rightarrow h[n]\{x[k-n]+x[k-M+n]\}

先加后乘,乘法器数从 NN 降到 N/2\lceil N/2\rceil

11. 考场策略

2022–2025 真题分布

年份Q1Q2Q3Q4Q5Q6
2022线性相位FIR(ch7)DFT/FFT(ch5)采样恢复(ch3)FIR设计(ch10)
2023H(z)零极点(ch6)DFT对称(ch5)线性相位FIR设计(ch7)IIR设计(ch9)
2024DFT性质(ch5)零极点+hn最小相位(ch7)FIR设计(ch10)
2025采样+(1)n(-1)^n(ch3)级联系统(ch4)系统辨识(ch4)H(z)+DF-II(ch6+ch8)DFT长度(ch5)自选设计(ch9/10)

必考公式 Top 10

  1. ω=2πf/Fs\omega=2\pi f/F_s(频率换算)
  2. X[k]=x[n]WNknX[k]=\sum x[n]W_N^{kn}(DFT 定义)
  3. X[k]=X[kN]X[k]=X^*[\langle -k\rangle_N](实序列对称)
  4. x[n]2=1NX[k]2\sum|x[n]|^2=\frac{1}{N}\sum|X[k]|^2(Parseval)
  5. (1)n=ejπn(-1)^n=e^{j\pi n}(频谱搬移 π\pi
  6. anu[n]1/(1az1)a^n u[n]\leftrightarrow 1/(1-az^{-1})z>a|z|>|a|(Z 变换)
  7. Ω=2Ttan(ω/2)\Omega=\frac{2}{T}\tan(\omega/2)(预畸变)
  8. Nlog[(A21)/ϵ2]2log(Ωs/Ωp)N\ge\frac{\log[(A^2-1)/\epsilon^2]}{2\log(\Omega_s/\Omega_p)}(Butterworth 阶数)
  9. hd[m]=sin(ωcm)/(πm)h_d[m]=\sin(\omega_c m)/(\pi m)(理想低通)
  10. h[n]=hd[nα]w[n]h[n]=h_d[n-\alpha]w[n](窗函数法)

时间分配建议(2 小时 / 100 分)

  • 每个 Q 约 20–25 分钟
  • 先做最熟的题型,拿稳基础分
  • 设计题(FIR/IIR)步骤多但套路固定,留足 25 分钟
  • 最后检查:频率单位、1/N1/N、ROC、阶数取整、对称类型