有限冲击响应数字滤波器.pptVIP

有限冲激响应数字滤波器 数字滤波器主要有无限冲击响应（ Infinite Impulse Response，IIR）滤波器和有限冲击响应（Finite Impulse Response，FIR）滤波器。其中，IIR滤波器是指单位冲击响应长度为无限长的滤波器，FIR滤波器是指单位冲击响应长度为有限长的滤波器。 FIR滤波器的单位冲击响应h(n)定义为： 如果希望得到的理想滤波器的频率响应为 那么FIR 滤波器的设计就在于寻找一个传递函数： 根据具体的滤波要求，计算出FIR滤波器的系数是FIR滤波器设计的主要内容之一。计算滤波器系数的方法多种多样，常用的方法有窗函数法、频率采样法及最优化方法。 利用最优化方法计算滤波器系数 1、根据给定的滤波器性能参数，确定对应的理想滤波器的频率响应 * * 式中，N称为滤波器的长度。N-1是滤波器的介数，因为长度为N的FIR滤波器零点的个数为N-1. FIR滤波器的传递函数为h(n)的Z变换，用数学公式表示为： 对于FIR滤波器来说，单位冲击响应h(n) (n=0,1,2,..N-1)也即滤波器的系数。 FIR滤波器的输入/输出关系可用如下的差分方程描述： 去逼近 ，逼近的方法有以下三种： 1、窗口设计法(时域逼近) 2、频率采样法(频域逼近) 3、最优化设计法(等波纹逼近) 窗函数法：选取某一种合适的理想频率滤波器，然后对理想的冲击响应加窗进行截断来计算滤波器的系数。这种方法主要是选择恰当的窗函数和合适的理想滤波器。 频率采样法：直接在频域对理想滤波器进行采样，然后利用IDFT的方法计算滤波器系数。 最优化方法：利用了线性规划技术，使得实际滤波器和理想滤波器的逼近误差在全频带区间上是均匀分布 的，能在满足相同技术指标的情况下，所需要的滤波系数最少。 一个截止频率为 的线性理想低通滤波器的频率响应为： 为滤波器的延时常数 2、确定滤波器的长度N 设计滤波器之前，滤波器的长度是不知道的。但是要求：通带越平，阻带衰减越大，过渡带越窄，滤波器的长度越大。下面的公式可以用来简单估算滤波器的长度N 3、计算加权函数 若用 表示理想滤波器的频率响应 , 表示实际滤波器的响应，则理想滤波器和实际滤波器的误差可表示为： 表示误差加权函数，它的取值根据通带或阻带的逼近精度 要求不同而不同。通常，在逼近精度高的频带， 取值大；要求逼近 精度低的频带， 取值小。误差加权函数的取值通常由通带偏差 和阻带偏差 确定，用数学公式表示如下： 4、计算滤波器系数 根据前面步骤的结果，在Matlab里调用最优方法计算函数，得到一组滤波器系数h(n)。 在Matlab信号处工具箱中，最优化方法计算滤波器系数的函数是remez函数，其标准的调用格式为： 其中： h为待计算的滤波器系数； N为滤波器的长度； f为滤波器给定的边沿频率，其单位为pi； A为各指定的频率上理想的幅度响应。