fir滤波器之系统函数

1、fir滤波器设计的目标是找出系统的传递函数或单位冲击响应。常见的设计方法有窗函数法。

2、fir最有用的特点是它的线性相位。线性相位能够保证一个由多个频率组成的信号在通过滤波器后，信号的波形不发生变化。

3、fir滤波器的数学表示
差分方程：
$$y(n)=b_{0}x(n)+b_{1}x(n-1)+....+b_{k}x(n-k)=\sum _{k=0}^{M-1}{b}_{k}x(n-k)公式一$$
单位冲击响应：
可以将输出序列写成系统的单位冲击响应h(n)和输入信号的卷积形式：
$$y(n)=\sum _{k=0}^{M-1}h(k)x(n-k)公式二$$
fir滤波器同样能够用它的系统函数来表征：
$$H(z)=\sum _{k=0}^{M-1}h(k)z^{-k}公式三$$
由上述公式可看出，求得M个单位冲击响应系数h(k)与输入信号做卷积就能够得到系统输出。系数可以由指定期望的系统频率响应通过公式三求得。常用的是窗函数设计法。

4、使用窗函数设计线性相位fir滤波器
这种方法是先指定期望的频率响应 $H_d(w)$，然后求出相应的单位冲击响应$h_d(n)$。实际上 $H_d(w)$和$h_d(n)$是傅里叶变换对关系：
$$\begin{gathered} H_d(w)=\sum _{n=0}^{\infty }{h}_d(n)e^{-jwn} \\ {h}_d(n)=\frac{1}{2\pi }{\int }_{-\pi }^{\pi }{H}_d(w)e^{jwn}dw \end{gathered}$$
当指定$H_d(w)$的取值后，就可求出$h_d(n)$。
理想情况下$h_d(n)$是无限长的，需要将它截取一段长度M，称之为M点fir滤波器。截取M点相当于给$h_d(n)$乘一个矩形窗：
$$w(n)=\{\begin{array}{ll}1& \text{n}=0,1,...,M-1\\ 0& \text{o}thers\end{array}$$
于是单位冲击响应就变成：
$$h(n)=h_d(n)w(n)$$
由于时域的乘积等价于频域的卷积，所以可以得到$h(n)$的频域响应函数：
$$H(w)=\frac{1}{2\pi }{\int }_{-\pi }^{\pi }{H}_d(v)W(w-v)dv$$
先看窗函数的频域响应。
矩形窗：
$$W(w)=\sum _{n=0}^{M-1}1e^{-jwn}=e^{-jw\frac{(M-1)}{2}}\frac{\sin (\frac{wM}{2})}{\sin (\frac{w}{2})}$$

幅度响应：
$$|W(w)|=|\frac{\sin (\frac{wM}{2})}{\sin (\frac{w}{2})}|$$
相位响应：
$$\Theta (w)=-w(\frac{M-1}{2})$$
当M=20时

当M=100时

随着M的增大，主瓣的宽度减小。旁瓣随着M的增大而变窄，峰值变高，旁瓣的面积不随M变化。

blackman（M=20）：

Hamming（M=20）：

从上述描述可以看出不同窗函数对频域响应是有影响的。

5、窗函数的特点
1）因为频域响应是由窗函数卷积得到的，所以具有平滑作用。
2）矩形窗的截取会产生振铃现象，导致旁瓣增大。可以使用其他窗截取进行抑制。
3）其他窗的使用可以有效降低旁瓣，使整个频域响应更加平滑，但是主瓣会展宽，峰值降低，过渡带增加。

6、窗系数的求法：
常用求取滤波器系统的工具是使用matlab设计。