脉冲串函数在数字信号处理中的核心应用与价值

发布于:2025-09-15 ⋅ 阅读:(14) ⋅ 点赞:(0)

在数字信号处理(DSP)领域,脉冲串函数是连接信号维度转换、描述离散域关键操作的基础性数学工具。它以单位冲激信号 / 序列为核心,通过 “选通” 或 “加权” 机制,将复杂的信号处理操作(如采样、抽取、插值、频谱计算)转化为可借助傅里叶变换(FT/DTFT/DFT)分析的形式,解决了 “连续与离散信号衔接” 及 “离散操作数学建模” 两大核心难题,是 DSP 理论体系与工程实践中不可或缺的关键组件。

一、脉冲串函数的定义与核心形式

脉冲串函数本质是由单位冲激信号 / 序列构成的周期性或非周期性集合,其核心特征是 “在特定时刻 / 索引处保留信号值,其余时刻 / 索引处为零”,从而实现对信号的 “筛选” 或 “采样”。根据作用对象的不同,DSP 中常用的脉冲串函数主要分为两类:

1. 连续时间脉冲串(采样脉冲)

连续时间脉冲串主要用于将连续信号转化为离散信号,其数学表达式为:
p(t)=∑n=−∞∞δ(t−nTs) p(t) = \sum_{n=-\infty}^{\infty} \delta(t - nT_s) p(t)=n=δ(tnTs)

其中,TsT_sTs为采样周期(对应采样频率 fs=1/Tsf_s = 1/T_sfs=1/Ts),δ(t)\delta(t)δ(t) 为连续时间单位冲激信号(满足 ∫−∞∞δ(t)dt=1\int_{-\infty}^{\infty} \delta(t)dt = 1δ(t)dt=1,仅在 t=0t=0t=0 处有非零值)。该脉冲串的核心特征是周期性,仅在 t=nTst = nT_st=nTsnnn 为整数)处产生冲激,恰好对应连续信号的采样时刻。

2. 离散时间脉冲串(选通脉冲)

离散时间脉冲串主要用于离散信号的多速率处理(如抽取、插值),其数学表达式为:
pd[n]=∑k=−∞∞δ[n−kM] p_d[n] = \sum_{k=-\infty}^{\infty} \delta[n - kM] pd[n]=k=δ[nkM]

其中,MMM 为周期(对应抽取 / 插值间隔),δ[n]\delta[n]δ[n] 为离散时间单位脉冲序列(满足 δ[0]=1\delta[0] = 1δ[0]=1δ[n]=0\delta[n] = 0δ[n]=0n≠0n \neq 0n=0))。该脉冲串同样具有周期性,仅在 n=kMn = kMn=kMkkk 为整数)处为 1,其余索引处为 0,可实现对离散信号特定索引值的 “筛选”。

二、脉冲串函数的四大核心应用场景

脉冲串函数的应用贯穿 DSP 的关键环节,从信号的初始采样到后续的速率转换、频谱计算,均依赖其实现数学建模与频域分析。以下结合具体原理与实例,详细阐述其核心应用:

1. 连续信号采样:实现 “连续→离散” 的维度转换

采样是 DSP 的起点,需将模拟域的连续信号 xa(t)x_a(t)xa(t) 转化为数字域的离散序列 x[n]=xa(nTs)x[n] = x_a(nT_s)x[n]=xa(nTs)。脉冲串函数在此过程中扮演 “采样模板” 的角色,通过与连续信号相乘,精准 “选通” 采样时刻的信号值,同时为频域分析提供数学基础。

(1)数学建模

连续信号 xa(t)x_a(t)xa(t) 与采样脉冲串 p(t)p(t)p(t) 相乘,得到 “采样后的连续信号”(本质是离散序列的连续域表示):
xa∗(t)=xa(t)⋅p(t)=xa(t)⋅∑n=−∞∞δ(t−nTs) x_a^*(t) = x_a(t) \cdot p(t) = x_a(t) \cdot \sum_{n=-\infty}^{\infty} \delta(t - nT_s) xa(t)=xa(t)p(t)=xa(t)n=δ(tnTs)

根据连续冲激信号的 “筛选性质”(f(t)δ(t−t0)=f(t0)δ(t−t0)f(t)\delta(t-t_0) = f(t_0)\delta(t-t_0)f(t)δ(tt0)=f(t0)δ(tt0)),上式可简化为:
xa∗(t)=∑n=−∞∞xa(nTs)δ(t−nTs) x_a^*(t) = \sum_{n=-\infty}^{\infty} x_a(nT_s) \delta(t - nT_s) xa(t)=n=xa(nTs)δ(tnTs)

该式清晰体现:采样后的信号仅在 t=nTst = nT_st=nTs 处保留 $ x_a(nT_s)$(即离散序列 x[n]x[n]x[n] 的幅值),其余时刻为零,完美实现了 “连续信号离散化” 的数学描述。

(2)频域分析价值

若直接定义 x[n]=xa(nTs)x[n] = x_a(nT_s)x[n]=xa(nTs),仅能描述离散序列的幅值,无法分析采样过程的频域特性(如混叠、采样定理)。而采样脉冲串 p(t)p(t)p(t) 的傅里叶变换为:
P(jΩ)=2πTs∑k=−∞∞δ(Ω−kΩs) P(j\Omega) = \frac{2\pi}{T_s} \sum_{k=-\infty}^{\infty} \delta(\Omega - k\Omega_s) P(jΩ)=Ts2πk=δ(ΩkΩs)

(其中 Ωs=2π/Ts\Omega_s = 2\pi/T_sΩs=2π/Ts 为采样角频率)。结合傅里叶变换的 “乘积→卷积” 性质,可得到采样信号的频域表达式:
Xa∗(jΩ)=12πXa(jΩ)∗P(jΩ)=1Ts∑k=−∞∞Xa(j(Ω−kΩs)) X_a^*(j\Omega) = \frac{1}{2\pi} X_a(j\Omega) * P(j\Omega) = \frac{1}{T_s} \sum_{k=-\infty}^{\infty} X_a(j(\Omega - k\Omega_s)) Xa(jΩ)=2π1Xa(jΩ)P(jΩ)=Ts1k=Xa(j(ΩkΩs))

该式正是奈奎斯特采样定理的核心数学依据:当 Ωs>2Ωm\Omega_s > 2\Omega_mΩs>2ΩmΩm\Omega_mΩm 为信号最高角频率)时,各周期延拓的频谱不重叠,可通过低通滤波无失真恢复原连续信号;若采样率不足,则会出现频谱混叠。若无脉冲串函数,这一关键频域关系将无法建立。

2. 离散信号抽取:保障 “降采样” 的频域分析严谨性

抽取(Downsampling)是多速率信号处理中 “降低采样率” 的操作,定义为 “对离散序列 x[n]x[n]x[n] 每隔 M−1M-1M1 个点取 1 个值”,即抽取后序列 y[n]=x[Mn]y[n] = x[Mn]y[n]=x[Mn]($ M $ 为抽取因子)。脉冲串函数在此处的核心作用是将 “间隔选点” 这一离散操作转化为 “序列相乘”,从而关联抽取前后信号的 DTFT,分析频谱混叠问题。

(1)数学建模

构造周期为 MMM 的离散脉冲串 pd[n]=∑k=−∞∞δ[n−kM]p_d[n] = \sum_{k=-\infty}^{\infty} \delta[n - kM]pd[n]=k=δ[nkM],抽取操作可通过 “原始序列与脉冲串相乘” 实现:
xp[n]=x[n]⋅pd[n]=x[n]⋅∑k=−∞∞δ[n−kM]=∑k=−∞∞x[kM]δ[n−kM] x_p[n] = x[n] \cdot p_d[n] = x[n] \cdot \sum_{k=-\infty}^{\infty} \delta[n - kM] = \sum_{k=-\infty}^{\infty} x[kM] \delta[n - kM] xp[n]=x[n]pd[n]=x[n]k=δ[nkM]=k=x[kM]δ[nkM]

该式中,xp[n]x_p[n]xp[n] 是 “完整索引的序列”,但仅在 n=kMn = kMn=kM 处保留 x[kM]x[kM]x[kM](即抽取序列的幅值),其余索引处为零。通过 “重索引”(令 n=Mn′n = Mn'n=Mn),可直接得到抽取序列 y[n′]=xp[Mn′]=x[Mn′]y[n'] = x_p[Mn'] = x[Mn']y[n]=xp[Mn]=x[Mn],完成了抽取操作的数学具象化。

(2)频域分析价值

若不引入脉冲串,直接对 y[n]=x[Mn]y[n] = x[Mn]y[n]=x[Mn] 求 DTFT,无法关联 Y(ejω)Y(e^{j\omega})Y(e)X(ejω)X(e^{j\omega})X(e)。而通过脉冲串 pd[n]p_d[n]pd[n],可先计算 xp[n]x_p[n]xp[n] 的 DTFT:
Xp(ejω)=DTFT{x[n]⋅pd[n]}=12πX(ejω)∗Pd(ejω) X_p(e^{j\omega}) = \text{DTFT}\{x[n] \cdot p_d[n]\} = \frac{1}{2\pi} X(e^{j\omega}) * P_d(e^{j\omega}) Xp(e)=DTFT{x[n]pd[n]}=2π1X(e)Pd(e)

其中,离散脉冲串的 DTFT 为 Pd(ejω)=2πM∑k=0M−1δ(ω−2πkM)P_d(e^{j\omega}) = \frac{2\pi}{M} \sum_{k=0}^{M-1} \delta\left(\omega - \frac{2\pi k}{M}\right)Pd(e)=M2πk=0M1δ(ωM2πk)。代入卷积性质展开后可得:
Xp(ejω)=1M∑k=0M−1X(ej(ω−2πkM)) X_p(e^{j\omega}) = \frac{1}{M} \sum_{k=0}^{M-1} X\left(e^{j\left(\omega - \frac{2\pi k}{M}\right)}\right) Xp(e)=M1k=0M1X(ej(ωM2πk))

由于 y[n]=xp[Mn]y[n] = x_p[Mn]y[n]=xp[Mn],根据 DTFT 的 “时移缩放性质”,抽取序列的 DTFT 为:
Y(ejω)=Xp(ejω/M)=1M∑k=0M−1X(ej(ω−2πkM)) Y(e^{j\omega}) = X_p(e^{j\omega/M}) = \frac{1}{M} \sum_{k=0}^{M-1} X\left(e^{j\left(\frac{\omega - 2\pi k}{M}\right)}\right) Y(e)=Xp(e/M)=M1k=0M1X(ej(Mω2πk))

该式揭示:抽取后的频谱是原频谱 “展宽 MMM 倍” 并 “周期延拓叠加” 的结果,若原信号带宽超过 π/M\pi/Mπ/M,会出现频谱混叠 —— 这一结论为抽取前需加低通滤波器(抗混叠滤波)提供了理论依据,而脉冲串函数是推导该关系的唯一数学桥梁。

3. 离散信号插值:支撑 “升采样” 的频谱重构

插值(Upsampling)是 “提高采样率” 的操作,定义为 “在离散序列相邻点间插入 L−1L-1L1 个零”(LLL 为插值因子),即插值后序列 y[n]=x[n/L]y[n] = x[n/L]y[n]=x[n/L]nnnLLL 的整数倍时取 x[n/L]x[n/L]x[n/L],否则为零)。脉冲串函数在此处的作用是描述 “插入零” 的操作,同时揭示频谱扩展规律,为后续滤波重构提供指导。

(1)数学建模

构造周期为 LLL 的离散脉冲串 pL[n]=∑k=−∞∞δ[n−kL]p_L[n] = \sum_{k=-\infty}^{\infty} \delta[n - kL]pL[n]=k=δ[nkL],插值序列可表示为:
y[n]=x[⌊nL⌋]⋅pL[n] y[n] = x\left[\left\lfloor \frac{n}{L} \right\rfloor\right] \cdot p_L[n] y[n]=x[Ln]pL[n]

该式表明:插值序列 y[n]y[n]y[n] 仅在 n=kLn = kLn=kL 处保留原序列 x[k]x[k]x[k] 的值,其余索引处填充零,与抽取操作的 “选通” 逻辑形成互补,精准描述了 “升采样” 的核心过程。

(2)频域分析价值

对插值序列 y[n]y[n]y[n] 求 DTFT,结合脉冲串的频域特性可得:
Y(ejω)=1L∑k=0L−1X(ej(ω−2πkL)) Y(e^{j\omega}) = \frac{1}{L} \sum_{k=0}^{L-1} X\left(e^{j\left(\omega - \frac{2\pi k}{L}\right)}\right) Y(e)=L1k=0L1X(ej(ωL2πk))

该式揭示:插值后的频谱是原频谱 “压缩 LLL 倍” 并沿 ω=2πk/L\omega = 2\pi k/Lω=2πk/L 周期延拓的结果(即频谱扩展)。若要恢复无失真的高采样率信号,需通过低通滤波(截止频率 π/L\pi/Lπ/L)滤除延拓产生的镜像频谱 —— 这一 “滤波需求” 的理论依据,正是通过脉冲串函数关联原序列与插值序列后推导得出的。

4. DFT 的周期延拓:连接 DTFT 与 DFT 的理论纽带

DFT 是 DSP 中实际计算信号频谱的工具,但 DFT 的本质是 “对离散序列的周期延拓求傅里叶级数(DFS)”。脉冲串函数通过描述 “周期延拓” 过程,揭示了 DFT 与 DTFT 的内在联系,为频谱计算提供了理论支撑。

(1)数学建模

离散序列 x[n]x[n]x[n] 的 DTFT 是 ω\omegaω 的连续周期函数(周期 2π2\pi2π),而 DFT 是 DTFT 在 ω=2πk/N\omega = 2\pi k/Nω=2πk/NNNN 为 DFT 点数)处的采样值。DFT 的计算基于 “周期延拓序列” x~[n]=∑k=−∞∞x[n−kN]\tilde{x}[n] = \sum_{k=-\infty}^{\infty} x[n - kN]x~[n]=k=x[nkN],该延拓过程可通过脉冲串卷积实现:
x~[n]=x[n]∗pN[n] \tilde{x}[n] = x[n] * p_N[n] x~[n]=x[n]pN[n]

其中,pN[n]=∑k=−∞∞δ[n−kN]p_N[n] = \sum_{k=-\infty}^{\infty} \delta[n - kN]pN[n]=k=δ[nkN] 是周期为 NNN 的离散脉冲串。卷积操作的物理意义是:将原序列 x[n]x[n]x[n] 沿 n=kNn = kNn=kN 平移并叠加,形成周期为 NNN 的序列 x~[n]\tilde{x}[n]x~[n]

(2)理论关联价值

对周期延拓序列 x~[n]\tilde{x}[n]x~[n] 求 DFS,其系数 X~[k]\tilde{X}[k]X~[k] 满足:
X~[k]=∑n=0N−1x~[n]e−j2πkn/N=∑n=0N−1x[n]e−j2πkn/N=X[k] \tilde{X}[k] = \sum_{n=0}^{N-1} \tilde{x}[n] e^{-j2\pi kn/N} = \sum_{n=0}^{N-1} x[n] e^{-j2\pi kn/N} = X[k] X~[k]=n=0N1x~[n]ej2πkn/N=n=0N1x[n]ej2πkn/N=X[k]

(即 DFS 系数等于 DFT 值)。同时,通过脉冲串的卷积性质,x~[n]=x[n]∗pN[n]\tilde{x}[n] = x[n] * p_N[n]x~[n]=x[n]pN[n] 的 DTFT 为:
X~(ejω)=X(ejω)⋅PN(ejω)=X(ejω)⋅2πN∑k=−∞∞δ(ω−2πkN) \tilde{X}(e^{j\omega}) = X(e^{j\omega}) \cdot P_N(e^{j\omega}) = X(e^{j\omega}) \cdot \frac{2\pi}{N} \sum_{k=-\infty}^{\infty} \delta\left(\omega - \frac{2\pi k}{N}\right) X~(e)=X(e)PN(e)=X(e)N2πk=δ(ωN2πk)

该式表明:周期序列的 DTFT 是原序列 DTFT 在 ω=2πk/N\omega = 2\pi k/Nω=2πk/N 处的冲激采样,而 DFT 值正是这些冲激的强度 —— 脉冲串函数在此处完美连接了 “非周期序列的连续频谱(DTFT)” 与 “周期序列的离散频谱(DFT)”,解释了 DFT 的物理意义,为工程中通过 DFT 计算频谱提供了理论依据。

三、脉冲串函数的核心价值总结

脉冲串函数在 DSP 中的应用,本质是解决 “信号处理操作的数学建模” 与 “频域分析可行性” 两大问题,其核心价值可归纳为三点:

  1. 维度转换桥梁:通过连续时间脉冲串,实现连续信号到离散信号的采样转换,同时保留频域可分析性,为采样定理的推导与应用奠定基础;

  2. 离散操作具象化:将抽取(间隔选点)、插值(插入零)等离散操作转化为 “序列 × 脉冲串” 的数学形式,借助傅里叶变换的卷积性质,分析操作对频谱的影响(如混叠、扩展),为多速率处理的滤波器设计提供指导;

  3. 理论关联纽带:通过离散脉冲串的卷积操作,实现离散序列的周期延拓,连接 DTFT(连续频谱)与 DFT(离散频谱),揭示 DFT 的本质,为实际频谱计算提供理论支撑。

综上,脉冲串函数是 DSP 从 “直观操作” 走向 “数学严谨性” 的关键工具,其应用贯穿信号处理的全流程,是理解与掌握 DSP 理论及工程实践的核心基础。