数字信号处理笔记09:数字滤波器设计方法

一、模拟滤波器设计

模拟滤波器指标

数字信号处理笔记09:数字滤波器设计方法_第1张图片

\Omega_P:通带截止频率

\Omega_S:阻带起始频率

\delta_P:通带内最大衰减

\delta_S:阻带内最小衰减

由幅度平方函数确定系统函数

|H(j\Omega)|^2 = H(j\Omega) H^*(j\Omega)

h(t)\leftrightarrow H(j\Omega) \leftrightarrow H(s)

h^*(-t) \leftrightarrow H^*(j\Omega) \leftrightarrow H^*(-s^*)

数字信号处理笔记09:数字滤波器设计方法_第2张图片

确定|H(j\Omega)|^2后,选择一半的关于虚轴对称的零极点则可组成H(j\Omega)

低通巴特沃斯滤波器

|H(j\Omega)|^2 = \frac{1}{1+\left( \frac{\Omega}{\Omega_c} \right )^{2N}}

数字信号处理笔记09:数字滤波器设计方法_第3张图片

设计低通滤波器,要求如下:

 解:

0.89125^2 = \frac{1}{1+\left( \frac{2\pi 400}{\Omega_c} \right )^{2N}}

0.17783^2 = \frac{1}{1+\left( \frac{2\pi 600}{\Omega_c} \right )^{2N}}

化简得

\left( \frac{2\pi 400}{\Omega_c} \right )^{2N} = \frac{1}{0.89125^2}-1
\left( \frac{2\pi 600}{\Omega_c} \right )^{2N} = \frac{1}{0.17783^2} -1

计算可得

N = \frac{1}{2} \frac{ln\left( \frac{1-\frac{1}{0.89125^2}}{1-\frac{1}{0.17783^2}} \right )}{ln(\frac{2}{3})} = 5.88574

取N为6,计算可得到\Omega_c = 2813

因此,|H(j\Omega)|^2 = \frac{1}{1+\left( \frac{\Omega}{2813} \right )^{12}} = \frac{1}{1+(\frac{j\Omega}{j2813})^{12}}

那么,|H(s)|^2 = \frac{1}{1+\left( \frac{s}{j2813} \right)^{12}}

极点:

1+\left[\frac{s_p}{j2813}\right]^{12} = 0

可得:s_p = (-1)^{\frac{1}{12}}j2813 = \left( e^{j2k\pi}e^{j\pi} \right )^{\frac{1}{12}} e^{j\frac{\pi}{2}} 2813 =2813 e^{j\frac{2k+7}{12}\pi}

数字信号处理笔记09:数字滤波器设计方法_第4张图片

取左半平面的极点

H(s) = \frac{2813^6}{(s-2813e^{j\frac{7}{12}\pi})(s-2813e^{-j\frac{7}{12}\pi})(s-2813e^{j\frac{3}{4}\pi})(s-2813e^{-j\frac{3}{4}\pi})(s-2813e^{j\frac{11}{12}\pi})(s-2813e^{-j\frac{11}{12}\pi})}

数字信号处理笔记09:数字滤波器设计方法_第5张图片

数字信号处理笔记09:数字滤波器设计方法_第6张图片

频率转换与高通、带通、带阻滤波器设计

数字信号处理笔记09:数字滤波器设计方法_第7张图片

 二、离散时间系统对连续时间信号进行滤波

数字信号处理笔记09:数字滤波器设计方法_第8张图片

 X(e^{jw}) = \frac{1}{T_s} \sum_{k = -\infty}^{\infty} X_c\left( j \frac{w}{T_s} - j \frac{2\pi k}{T_s} \right )

Y(e^{jw} ) = X(e^{jw})H(e^{jw}) = \frac{1}{T_s} \sum_{k = -\infty}^{\infty} X_c\left( j \frac{w}{T_s} - j \frac{2\pi k}{T_s} \right ) H(e^{jw})

y[n]转化为连续时间信号y_s(t)

Y_s(j\Omega) = Y(e^{jw})|_{w = \Omega T_s} = Y(e^{j\Omega T_s}) \\= \frac{1}{T_s} \sum_{k=-\infty}^{\infty} X_c \left( j\Omega - j\frac{2\pi k}{T_s} \right) H(e^{j\Omega T_s})

连续时间信号y_s(t)经增益为T_s,通带为\left( -\frac{\Omega_s}{2}, \frac{\Omega_s}{2} \right )的重构滤波器h_r(t)后得到重建的连续时间信号y_r(t)

Y_r(j\Omega) = H_r(j\Omega) Y_s(j\Omega) \\= H_r(j\Omega) \frac{1}{T_s} \sum_{k=-\infty}^{\infty} X_c \left( j\Omega - j\frac{2\pi k}{T_s} \right) H(e^{j\Omega T_s}) \\= X_c(j\Omega)H(e^{j\Omega T_s})

|\Omega| < \frac{\Omega_s}{2}

数字滤波器H(e^{jw})的模拟等效滤波器:H_{eff}(j\Omega) = H(e^{j\Omega T_s}),|\Omega| < \frac{\pi}{T_s}

三、冲激响应不变法设计IIR滤波器 

数字滤波器设计是根据幅频响应|H(e^{jw})|,设计出满足要求的系统函数H(z)。根据H(z)可以画出信号流图,在数字系统中编程实现。

IIR数字滤波器的设计:

(1)将待设计的数字滤波器的性能指标w转换为模拟滤波器的性能指标\Omega

(2)设计满足要求的模拟滤波器H(s)

(3)将模拟滤波器转换为数字滤波器H(z)

冲击响应不变法把H(s)s平面映射到z平面的H(z)z = e^{sT_s}

 对于有理分式

H(s) = \sum_{k=1}^{N} \frac{A_k}{s-s_k}

h(t) = \sum_{k=1}^{N} A_k e^{s_k t} u(t)

h[n] = h(nT_s) =\sum_{k=1}^{N} A_k e^{s_k nT_s} u(nT_s) = \sum_{k=1}^{N} A_k e^{s_k T_s n} u(n)

数字信号处理笔记09:数字滤波器设计方法_第9张图片

四、双线性变换法设计IIR滤波器

数字信号处理笔记09:数字滤波器设计方法_第10张图片

 五、窗函数法设计FIR滤波器

理想数字低通滤波器:

 H_d(e^{jw}) = 1, |w| < w_c

h_d[n] = \frac{1}{2\pi} \int_{-\pi}^{\pi} H_d(e^{jw}) e^{jwn} dw \\ = \frac{1}{2\pi} \int_{-w_c}^{w_c} e^{jwn} dw \\= \frac{1}{2\pi} \left[ \frac{e^{jwn}}{jn} \right ]_{-w_c}^{w_c} \\= \frac{1}{2\pi} \frac{e^{jw_c n} - e^{-jw_c n}}{jn} \\= \frac{1}{2\pi} \frac{2jsin(w_c n)}{jn} \\= \frac{sin(w_c n)}{\pi n}

通过移位和加窗的方法,把理想滤波器的单位脉冲响应h_d[n]变成一个有限长的因果序列h[n]

h[n] = h_d[n-M/2] w[n]

其中,w[n] = 1, 0\leq n \leq M

对于w[n],其DTFT为

W(e^{jw}) = \sum_{n=-\infty}^{\infty} w[n] e^{-jwn} = \sum_{n=0}^{M} e^{-jwn} = \frac{1-e^{-jw(M+1)}}{1-e^{-jw}} \\ = \frac{e^{-jw(M+1)/2}}{e^{-jw/2}} \frac{e^{jw(M+1)/2} - e^{-jw(M+1)/2}}{e^{jw/2} - e^{-jw/2}} \\= e^{-jwM/2} \frac{sin\left[w \frac{M+1}{2} \right ]}{sin\left[ \frac{w}{2} \right ]}

对于h_d[n-M/2],其DTFT为H_d(e^{jw}) e^{-jwM/2}

因为h[n] = h_d[n-M/2] w[n]

所以,

H(e^{jw}) = \frac{1}{2\pi} \int_{-\pi}^{\pi} W(e^{j(w-\theta)}) H_d(e^{j\theta}) e^{-j\theta M/2} d\theta \\= \frac{1}{2\pi} \int_{-w_c}^{w_c} e^{-j(w-\theta)M/2} \frac{sin\left[(w-\theta) \frac{M+1}{2} \right ]}{sin\left[ \frac{w-\theta}{2} \right ]} e^{-j\theta M/2} d\theta \\= \frac{e^{-jwM/2}}{2\pi} \int_{-w_c}^{w_c} \frac{sin\left[(w-\theta) \frac{M+1}{2} \right ]}{sin\left[ \frac{w-\theta}{2} \right ]} d\theta

数字信号处理笔记09:数字滤波器设计方法_第11张图片

数字信号处理笔记09:数字滤波器设计方法_第12张图片

 常用窗函数

  • 矩形窗

w[n] = \left\{\begin{matrix} 1, & 0\leq n \leq M \\ 0, & others \end{matrix}\right.

数字信号处理笔记09:数字滤波器设计方法_第13张图片

  •  Bartlett窗

w[n] = \left\{\begin{matrix} 2n/M & 0\leq n \leq M/2 \\ 2-2n/M & M/2 < n \leq M\\ 0 & others \end{matrix}\right.

数字信号处理笔记09:数字滤波器设计方法_第14张图片

数字信号处理笔记09:数字滤波器设计方法_第15张图片

  •  Hanning窗

w[n] \left\{\begin{matrix} 0.5 - 0.5 cos(2\pi n/M) & 0 \leq n \leq M \\ 0 & others \end{matrix}\right.

数字信号处理笔记09:数字滤波器设计方法_第16张图片

  •  Hamming窗

w[n] = \left\{\begin{matrix} 0.54 - 0.46 cos(2\pi n /M) & 0\leq n \leq M \\ 0 & others \end{matrix}\right.

数字信号处理笔记09:数字滤波器设计方法_第17张图片

数字信号处理笔记09:数字滤波器设计方法_第18张图片

  •  Blackman窗

w[n] = \left\{\begin{matrix} 0.42 - 0.5 cos(2\pi n /M) + 0.08 cos(4\pi n /M) & 0 \leq n \leq M \\ 0 & {\rm others} \end{matrix}\right.

主瓣宽度 旁瓣峰值 过渡带近似宽度 阻带最小衰减 等效Kaiser窗\beta
矩形 \frac{4\pi}{M+1} -13dB \frac{1.8\pi}{M+1} -21dB 0
Bartlett \frac{8\pi}{M+1} -25dB \frac{4.1\pi}{M+1} -25dB 1.33
Hanning \frac{8\pi}{M+1} -31dB \frac{6.2\pi}{M+1} -44dB 3.86
Hamming \frac{8\pi}{M+1} -41dB \frac{6.6\pi}{M+1} -53dB 4.86
Blackman \frac{12\pi}{M+1} -57dB \frac{11\pi}{M+1} -74dB 7.04

设计以通带截止频率w_p= 0.4\piw_s = 0.6\pi,通带纹波\delta_1 = 0.01,阻带纹波\delta_2 = 0.001的低通滤波器。

最终的数字低通滤波器单位冲击响应函数为

h[n] = h_d[n-M/2] w[n] = \frac{sin \left[ w_c (n-M/2) \right ]}{\pi (n-M/2)} w[n]

上式中需要确定的参数为w_cMw[n]

(1) 确定窗函数形状w[n]

A = -20lg \delta_2 = 60dB

选择Blackman窗

w[n] = \left\{\begin{matrix} 0.42 - 0.5 cos(2\pi n /M) + 0.08 cos(4\pi n /M) & 0 \leq n \leq M \\ 0 & {\rm others} \end{matrix}\right.

(2)确定M

过渡带宽度\frac{11\pi}{M+1} \leq w_s - w_p,计算得到M \geq 54

(3) 确定wc

计算参数w_c = \frac{w_p+w_s}{2} = 0.5\pi

由此可得到滤波器的单位脉冲响应

h[n] = \frac{sin \left[ 0.5\pi (n-27) \right ]}{\pi (n-27)} w[n]

数字信号处理笔记09:数字滤波器设计方法_第19张图片

数字信号处理笔记09:数字滤波器设计方法_第20张图片

数字信号处理笔记09:数字滤波器设计方法_第21张图片

  •  Kaiser窗

w[n] = \left\{\begin{matrix} \frac{I_0\left[ \beta \left( 1 - (n-\alpha)^2/\alpha^2 \right )^{0.5} \right ]}{I_0(\beta)} & 0\leq n \leq M\\ 0& {\rm others} \end{matrix}\right.

\beta = \left\{\begin{matrix} 0.1102 (A-8.7) & A>50 \\ 0.5842(A-21)^{0.4} + 0.07886(A-21) & 21\leq A \leq 50 \\ 0 & A<21 \end{matrix}\right.

\alpha = \frac{M}{2}

M = \frac{A-8}{2.285 (w_s - w_p)}

A = -20lg \delta

I_0(\cdot)表示第一类零阶修正的贝塞尔函数

六、频率采样法设计FIR滤波器

假设某FIR系统的单位脉冲响应为h[n]h[n]的DTFT为H(e^{jw})h[n]的DFT为H[k]H[k] = H(e^{jw})|_{w = \frac{2\pi k}{ N}}

根据给定的滤波器的系统性能可确定系统的频率响应H(e^{jw}),进而计算出其采样值H[k],然后由H[k]的IDFT得到所设计滤波器的单位脉冲响应h[n]

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