本专栏:数学系的数字信号处理 的前置知识主要有:数学分析(傅立叶级数的部分),泛函分析( L p L^p Lp空间的部分)
连续信号、滤波器与采样定理
我们在数学上粗略地定义信号和滤波器,目的是快速理解相关知识的结构,而不是花费大量的时间在细节上。
本文中的 F \mathscr{F} F 符号均指傅立叶变换,我们也常记 F ( f ) \mathscr{F}(f) F(f) 为 f ^ \hat{f} f^
定义
信号:即把实数域映到复数域的函数 f : R → C f:\mathbb{R}\to\mathbb{C} f:R→C,一般是分段连续的
滤波器:把信号映为信号的变换
定义:时不变性
信号 f f f 的时不变算子是如下定义的 f a f_a fa
f a ( t ) = f ( t − a ) , a ∈ R f_a(t)=f(t-a),a\in\mathbb{R} fa(t)=f(t−a),a∈R滤波器 L L L 称为时不变的,若对任意信号 f f f 和实数 a a a,有 L ( f a ) = ( L f ) a L(f_a)=(Lf)_a L(fa)=(Lf)a
定理:线性时不变滤波器的结构
设 L L L 是分段连续信号 f f f 上的线性时不变滤波器,则存在可积函数 h h h ,使得对任意信号 f f f ,有 L ( f ) = f ∗ h L(f)=f*h L(f)=f∗h;
此时称 h h h 为冲激响应函数,称 h h h 的傅立叶变换 h ~ \tilde{h} h~ 为系统函数
证明思路
引理:(线性时不变滤波器使输入输出同频)
设线性时不变滤波器 L L L, t t t 为自变量,则 ∀ λ ∈ R \forall \lambda\in\mathbb{R} ∀λ∈R,存在函数 h h h ,使得
L ( e i λ t ) = 2 π h ~ ( λ ) e i λ t L(e^{i\lambda t})=\sqrt{2\pi}\tilde{h}(\lambda)e^{i\lambda t} L(eiλt)=2πh~(λ)eiλt
引理的证明思路:(不严格)
设 h λ ( t ) = L ( e i λ t ) h^{\lambda}(t)=L(e^{i\lambda t}) hλ(t)=L(eiλt),由时不变性,
h λ ( t − a ) = L ( e i λ ( t − a ) ) h^{\lambda}(t-a)=L(e^{i\lambda (t-a)}) hλ(t−a)=L(eiλ(t−a))由线性性 L ( e i λ ( t − a ) ) = e − i λ a h λ ( t ) L(e^{i\lambda (t-a)})=e^{-i\lambda a}h^{\lambda}(t) L(eiλ(t−a))=e−iλahλ(t)故 h λ ( t ) = e i λ a h λ ( t − a ) h^{\lambda}(t)=e^{i\lambda a}h^{\lambda}(t-a) hλ(t)=eiλahλ(t−a)取 a = t a=t a=t,得到
h λ ( t ) = e i λ t h λ ( 0 ) h^{\lambda}(t)=e^{i\lambda t}h^{\lambda}(0) hλ(t)=eiλthλ(0)再取 h ~ \tilde{h} h~,使其满足
h ~ ( λ ) = h λ ( 0 ) 2 π \tilde{h}(\lambda)=\frac{h^{\lambda}(0)}{\sqrt{2\pi}} h~(λ)=2πhλ(0)
注:这个证明并没有解释结论形式的来源,也未解释满足所要求 h ~ \tilde{h} h~ 的存在性
定理的证明(不严格)
L ( f ) ( t ) = L F − 1 F ( f ) ( t ) = L [ 1 2 π ∫ − ∞ ∞ f ^ ( λ ) e i λ t d λ ] L(f)(t)=L\mathscr{F}^{-1}\mathscr{F}(f)(t)=L[\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}\hat{f}(\lambda)e^{i\lambda t}\mathrm{d}\lambda] L(f)(t)=LF−1F(f)(t)=L[2π