学习笔记-频率域滤波(1)-基本概念

基本概念

复数

由欧拉公式 $e^{j\theta }=\cos \theta +j\sin \theta$ 可推出：
$$C=R+jI$$

$$C=|C|(\cos \theta +j\sin \theta )$$

$$C=|C|e^{j\theta }$$

傅里叶级数

周期为T的周期函数$f(t)$可表示为傅里叶级数：
$$f(t)=\sum _{n=-\infty }^{\infty }{c}_n{e}^{j\frac{2\pi n}{T}t}$$

$$c_n=\frac{1}{T}{\int }_{-T/2}^{T/2}f(t)e^{-j\frac{2\pi n}{T}t}\mathrm{d}t$$

离散冲激

令 $x$ 表示一个离散变量。单位离散冲激 $\delta (x)$ 定义为：
$$\delta (x)=\{\begin{array}{l}1,x=0\\ 0,x\ne 0\end{array}$$
其取样性质的形式为：
$$\sum _{x=-\infty }^{\infty }f(x)\delta (x)=f(0)$$
或更一般地使用 $x=x_0$ 处的离散冲激形式：
$$\sum _{x=-\infty }^{\infty }f(x)\delta (x-x_0)=f(x_0)$$

傅里叶变换

连续单变量 $t$ 的连续函数 $f(t)$ 的傅里叶变换由 ℑ { f ( t ) } \Im\{f(t)\} ℑ{f(t)} 表示：
F ( μ ) = ℑ { f ( t ) } = ∫ − ∞ ∞ f ( t ) e − j 2 π μ t d t (1) F(\mu)=\Im\{f(t)\}=\int_{-\infty}^\infty f(t)e^{-j2\pi\mu t}\mathrm dt \tag{1} F(μ)=ℑ{f(t)}=∫−∞∞​f(t)e−j2πμtdt(1)

f ( t ) = ℑ − 1 { F ( μ ) } = ∫ − ∞ ∞ F ( μ ) e j 2 π μ t d μ (2) f(t)=\Im^{-1}\{F(\mu)\}=\int_{-\infty}^\infty F(\mu)e^{j2\pi\mu t}\mathrm d\mu \tag{2} f(t)=ℑ−1{F(μ)}=∫−∞∞​F(μ)ej2πμtdμ(2)
于是，$t=t_0$ 处的冲激的傅里叶变换为：
F ( μ ) = ℑ { δ ( t − t 0 ) } = ∫ − ∞ ∞ δ ( t − t 0 ) e − j 2 π μ t d t = 仅 在 t 0 处 有 值 e − j 2 π μ t 0 F(\mu)=\Im\{\delta(t-t_0)\}=\int_{-\infty}^\infty\delta(t-t_0)e^{-j2\pi\mu t}\mathrm dt\xlongequal{仅在t_0处有值}e^{-j2\pi\mu t_0} F(μ)=ℑ{δ(t−t0​)}=∫−∞∞​δ(t−t0​)e−j2πμtdt仅在t0​处有值 e−j2πμt0​
而对 $e^{j2\pi at}$ 做傅里叶变换为：
F ( μ ) = ℑ { e j 2 π a t } = ∫ − ∞ ∞ e j 2 π a t e − j 2 π μ t d t = ∫ − ∞ ∞ e j 2 π ( a − μ ) t d t = δ ( μ − a ) F(\mu)=\Im\{e^{j2\pi at}\}=\int_{-\infty}^\infty e^{j2\pi at}e^{-j2\pi\mu t}\mathrm dt=\int_{-\infty}^\infty e^{j2\pi(a-\mu)t}\mathrm dt=\delta(\mu-a) F(μ)=ℑ{ej2πat}=∫−∞∞​ej2πate−j2πμtdt=∫−∞∞​ej2π(a−μ)tdt=δ(μ−a)
因为当 $\mu \ne a$ 时上述积分式 $e^{j\theta }$ 为复平面内的单位圆 积分为零
由此可求取冲激串 $s_{\Delta T}(t)$ 的傅里叶变换 冲激串定义为一串冲激的和：
$$s_{\Delta T}(t)=\sum _{k=-\infty }^{\infty }\delta (t-k\Delta T)$$
冲激串是周期为 $\Delta T$ 的函数 可以表示为傅里叶级数：
$$s_{\Delta T}(t)=\sum _{n=-\infty }^{\infty }{c}_n{e}^{j\frac{2\pi n}{\Delta T}t}$$
式中 $c_n$ 代入 $s_{\Delta T}(t)$ ：
$$c_n=\frac{1}{\Delta T}{\int }_{-\Delta T/2}^{\Delta T/2}{s}_{\Delta T}(t)e^{-j\frac{2\pi n}{T}t}\mathrm{d}t$$
可发现在 $[-\frac{1}{\Delta T},\frac{1}{\Delta T}]$ 区间内仅有一个冲激 因此：
$$c_n=\frac{1}{\Delta T}{\int }_{-\Delta T/2}^{\Delta T/2}\delta (t)e^{-j\frac{2\pi n}{T}t}\mathrm{d}t=\frac{1}{\Delta T}{e}^0=\frac{1}{\Delta T}$$
于是傅里叶级数变为：
$$s_{\Delta T}(t)=\frac{1}{\Delta T}\sum _{n=-\infty }^{\infty }{e}^{j\frac{2\pi n}{\Delta T}t}$$
对其求傅里叶变换：
S ( μ ) = ℑ { s Δ T ( t ) } = ℑ { 1 Δ T ∑ n = − ∞ ∞ e j 2 π n Δ T t } = 1 Δ T ∑ n = − ∞ ∞ ℑ { e j 2 π n Δ T t } S(\mu)=\Im\{s_{\Delta T}(t)\}=\Im\left\{\frac1{\Delta T}\sum_{n=-\infty}^\infty e^{j\frac{2\pi n}{\Delta T}t}\right\}=\frac1{\Delta T}\sum_{n=-\infty}^\infty\Im\{e^{j\frac{2\pi n}{\Delta T}t}\} S(μ)=ℑ{sΔT​(t)}=ℑ{ΔT1​n=−∞∑∞​ejΔT2πn​t}=ΔT1​n=−∞∑∞​ℑ{ejΔT2πn​t}
在之前已得到 ℑ { e j 2 π a t } = δ ( μ − a ) \Im\{e^{j2\pi at}\}=\delta(\mu-a) ℑ{ej2πat}=δ(μ−a) 故有 ℑ { e j 2 π n Δ T t } = δ ( μ − n Δ T ) \Im\{e^{j\frac{2\pi n}{\Delta T}t}\}=\delta(\mu-\frac{n}{\Delta T}) ℑ{ejΔT2πn​t}=δ(μ−ΔTn​) 代入:
S ( μ ) = 1 Δ T ∑ n = − ∞ ∞ ℑ { e j 2 π n Δ T t } = 1 Δ T ∑ n = − ∞ ∞ δ ( μ − n Δ T ) S(\mu)=\frac1{\Delta T}\sum_{n=-\infty}^\infty\Im\{e^{j\frac{2\pi n}{\Delta T}t}\}=\frac1{\Delta T}\sum_{n=-\infty}^\infty\delta(\mu-\frac{n}{\Delta T}) S(μ)=ΔT1​n=−∞∑∞​ℑ{ejΔT2πn​t}=ΔT1​n=−∞∑∞​δ(μ−ΔTn​)
这表示周期为 $T$ 的冲激串 其傅里叶变换仍然为冲激串 其周期为 $\frac{1}{T}$

卷积

卷积算子由 $\star$ 表示，定义为：
$$(f\star h)(t)={\int }_{-\infty }^{\infty }f(\tau )h(t-\tau )\mathrm{d}\tau \text{\hspace{1em}(3)}$$
给定 $f(t)$ 的傅里叶变换为 $F(\mu )$， $h(t)$ 的傅里叶变换为 $H(\mu )$，对式 (3) 求傅里叶变换：
ℑ { ( f ⋆ h ) ( t ) } = ∫ − ∞ ∞ [ ∫ − ∞ ∞ f ( τ ) h ( t − τ ) d τ ] e − j 2 π μ t d t = 交 换 积 分 顺 序 ∫ − ∞ ∞ f ( τ ) [ ∫ − ∞ ∞ h ( t − τ ) e − j 2 π μ t d t ] d τ = ∫ − ∞ ∞ f ( τ ) [ ∫ − ∞ ∞ h ( t − τ ) e − j 2 π μ ( t − τ ) e − j 2 π μ τ d ( t − τ ) ] d τ = ∫ − ∞ ∞ f ( τ ) [ e − j 2 π μ τ ∫ − ∞ ∞ h ( x ) e − j 2 π μ ( x ) d x ] d τ = ∫ − ∞ ∞ f ( τ ) [   e − j 2 π μ τ H ( μ )   ] d τ = H ( μ ) ∫ − ∞ ∞ f ( τ ) e − j 2 π μ τ d τ = H ( μ ) F ( μ ) = ( H ⋅ F ) ( μ ) \begin{aligned} \Im\{(f\star h)(t)\}&=\int_{-\infty}^\infty\left[\int_{-\infty}^\infty f(\tau)h(t-\tau)\mathrm d\tau\right]e^{-j2\pi\mu t}\mathrm dt \\[4ex] &\xlongequal{交换积分顺序}\int_{-\infty}^\infty f(\tau)\left[\int_{-\infty}^\infty h(t-\tau)e^{-j2\pi\mu t}\mathrm dt\right]\mathrm d\tau \\[4ex] &=\int_{-\infty}^\infty f(\tau)\left[\int_{-\infty}^\infty h(t-\tau)e^{-j2\pi\mu(t-\tau)}e^{-j2\pi\mu\tau}\mathrm d(t-\tau)\right]\mathrm d\tau\\[4ex] &=\int_{-\infty}^\infty f(\tau)\left[e^{-j2\pi\mu\tau}\int_{-\infty}^\infty h(x)e^{-j2\pi\mu(x)}\mathrm dx\right]\mathrm d\tau\\[4ex] &=\int_{-\infty}^\infty f(\tau)[\ e^{-j2\pi\mu\tau}H(\mu)\ ]\mathrm d\tau\\[4ex] &=H(\mu)\int_{-\infty}^\infty f(\tau)e^{-j2\pi\mu\tau}\mathrm d\tau\\[4ex] &=H(\mu)F(\mu)=(H\cdot F)(\mu) \end{aligned} ℑ{(f⋆h)(t)}​=∫−∞∞​[∫−∞∞​f(τ)h(t−τ)dτ]e−j2πμtdt交换积分顺序 ∫−∞∞​f(τ)[∫−∞∞​h(t−τ)e−j2πμtdt]dτ=∫−∞∞​f(τ)[∫−∞∞​h(t−τ)e−j2πμ(t−τ)e−j2πμτd(t−τ)]dτ=∫−∞∞​f(τ)[e−j2πμτ∫−∞∞​h(x)e−j2πμ(x)dx]dτ=∫−∞∞​f(τ)[ e−j2πμτH(μ) ]dτ=H(μ)∫−∞∞​f(τ)e−j2πμτdτ=H(μ)F(μ)=(H⋅F)(μ)​
于是得到了卷积定理的一半：
$$(f\star h)(t)\iff (H\cdot F)(\mu )$$
类似地，可以推出另一半：
$$(f\cdot h)(t)\iff (H\star F)(\mu )$$
这说明频域的卷积类似于空间域的乘积，空间域的卷积类似于频域的乘积。