第一类线性相位条件约束——数字图像FIR滤波器设计对单位脉冲响应的要求

如果 θ ( ω ) \theta(\omega) θ(ω)为过原点的直线,则称为第一类线性相位,即
θ ( ω ) = − τ ω \theta(\omega) = -\tau \omega θ(ω)=τω
其中, τ \tau τ表示滤波器的群延时。

在第一类线性相位约束条件下,对 FIR 数字滤波器单位脉冲响应的要求。

H ( e j ω ) = ∑ n = 0 N − 1 h ( n ) e − j ω n = H ( ω ) e j θ ( ω ) = H ( ω ) e − j τ ω H({\rm e}^{{\rm j}\omega}) = \sum_{n=0}^{N-1} h(n) {\rm e}^{-{\rm j}\omega n} = H(\omega) {\rm e}^{{\rm j}\theta(\omega)} = H(\omega) {\rm e}^{-{\rm j}\tau\omega} H(ejω)=n=0N1h(n)ejωn=H(ω)ejθ(ω)=H(ω)ejτω

则有

∑ n = 0 N − 1 h ( n ) ( cos ⁡ n ω − j sin ⁡ n ω ) = H ( ω ) ( cos ⁡ τ ω − j sin ⁡ τ ω ) \sum_{n=0}^{N-1} h(n) (\cos n\omega - {\rm j} \sin n\omega) = H(\omega) (\cos \tau\omega - {\rm j} \sin \tau\omega) n=0N1h(n)(cosjsin)=H(ω)(cosτωjsinτω)

式中等号两边实部与实部相等,虚部与虚部相等,可得

∑ n = 0 N − 1 h ( n ) cos ⁡ n ω = H ( ω ) cos ⁡ τ ω (1) \sum_{n=0}^{N-1} h(n) \cos n\omega = H(\omega) \cos \tau\omega \tag{1} n=0N1h(n)cos=H(ω)cosτω(1)

∑ n = 0 N − 1 h ( n ) sin ⁡ n ω = H ( ω ) sin ⁡ τ ω (2) \sum_{n=0}^{N-1} h(n) \sin n\omega = H(\omega) \sin \tau\omega \tag{2} n=0N1h(n)sin=H(ω)sinτω(2)

将式(1)与式(2)相除,可得

∑ n = 0 N − 1 h ( n ) sin ⁡ n ω ∑ n = 0 N − 1 h ( n ) cos ⁡ n ω = sin ⁡ τ ω cos ⁡ τ ω (7.1.9) \frac{\sum\limits_{n=0}^{N-1} h(n) \sin n\omega}{\sum\limits_{n=0}^{N-1} h(n) \cos n\omega} = \frac{\sin \tau\omega}{\cos \tau\omega} \tag{7.1.9} n=0N1h(n)cosn=0N1h(n)sin=cosτωsinτω(7.1.9)

即可得

∑ n = 0 N − 1 h ( n ) sin ⁡ n ω cos ⁡ τ ω = ∑ n = 0 N − 1 h ( n ) cos ⁡ n ω sin ⁡ τ ω \sum\limits_{n=0}^{N-1} h(n) \sin n\omega \cos \tau\omega = \sum\limits_{n=0}^{N-1} h(n) \cos n\omega \sin \tau\omega n=0N1h(n)sincosτω=n=0N1h(n)cossinτω

利用三角函数公式可得

∑ n = 0 N − 1 h ( n ) sin ⁡ [ ( τ − n ) ω ] = 0 (3) \sum_{n=0}^{N-1} h(n) \sin[(\tau-n)\omega] = 0 \tag{3} n=0N1h(n)sin[(τn)ω]=0(3)

如果式(3)成立,则要求 h ( n ) sin ⁡ [ ( τ − n ) ω ] h(n) \sin[(\tau-n)\omega] h(n)sin[(τn)ω]对于求和区间中心位置呈奇对称。由于 sin ⁡ [ ( τ − n ) ω ] \sin[(\tau-n)\omega] sin[(τn)ω]对于 n = τ n=\tau n=τ呈奇对称关系,则要求 h ( n ) h(n) h(n)对于求和区间中心位置呈偶对称。因此,如果希望 FIR 数字滤波器满足第一类线性相位要求,则 h ( n ) h(n) h(n) τ \tau τ需要满足以下条件,

{ h ( n ) = h ( N − 1 − n ) , 0 ⩽ n ⩽ N − 1 (圆周偶对称) τ = ( N − 1 ) / 2 (4) \begin{cases} h(n) = h(N-1-n), & 0 \leqslant n \leqslant N-1(圆周偶对称) \\ \tau = (N-1)/2 \end{cases} \tag{4} {h(n)=h(N1n),τ=(N1)/20nN1(圆周偶对称)(4)

需要注意的是,式(4)是保证 FIR 数字滤波器具有第一类严格线性相位特性的充分条件,但不是必要条件。

线性相位的条件约束

类型 h ( n ) h(n) h(n)对称性 N N N的奇偶性 幅度函数对称性 相位函数 适用范围
第1类 偶对称 奇数 关于 ω = 0 , π , 2 π \omega=0, \pi, 2\pi ω=0,π,2π偶对称 第一类线性相位 低通、高通、带通、带阻
第2类 偶对称 偶数 关于 ω = π \omega=\pi ω=π为奇对称,关于 ω = 0 , 2 π \omega=0, 2\pi ω=0,2π为偶对称,且在 ω = π \omega=\pi ω=π H ( ω ) = 0 H(\omega)=0 H(ω)=0 θ ( ω ) = − N − 1 2 ω \theta(\omega)=-\frac{N-1}{2}\omega θ(ω)=2N1ω 低通、带通
第3类 奇对称 奇数 关于 ω = 0 , π , 2 π \omega=0, \pi, 2\pi ω=0,π,2π奇对称,且在 ω = 0 , π , 2 π \omega=0, \pi, 2\pi ω=0,π,2π H ( ω ) = 0 H(\omega)=0 H(ω)=0 第二类线性相位 带通
第4类 奇对称 偶数 关于 ω = π \omega=\pi ω=π为奇对称,关于 ω = 0 , 2 π \omega=0, 2\pi ω=0,2π为偶对称,且在 ω = 0 , π , 2 π \omega=0, \pi, 2\pi ω=0,π,2π H ( ω ) = 0 H(\omega)=0 H(ω)=0 θ ( ω ) = − N − 1 2 ω \theta(\omega)=-\frac{N-1}{2}\omega θ(ω)=2N1ω 低通、高通

你可能感兴趣的:(数字信号处理,图像处理)