FIR数字滤波器设计和实现幻灯片.pptVIP

;主题概述 ;5.1 概述：IIR 和 FIR 比较;5.1 概述：IIR 和 FIR 比较;5.1 概述：FIR DF 设计方法;5.1 概述：FIR DF 零极点; FIR DF 的频率响应为：;信号通过线性滤波器时，其幅度和相位可能会发生改变，滤波器增益 |H(ω)|和相位 θ(ω) 可能会随频率的变化而改变。 如：输入正弦信号 Acos(nω0) 则：输出为 |H(ω0)| Acos(nω0＋θ)，其中相移θ＝θ(ω0) 输出频率和输入频率相同，但幅度和相位都发生了变化 输出信号比输入信号滞后的样点数 n (位移) 可由下式求得: 设：nω0＋θ＝0 ; 对不同的频率有恒定的相移，会产生相位失真. ;主题概述 ;于是: ;可以证明，当 ;于是有：;可以证明，当 ;FIR滤波器单独满足恒定群延时的必要条件为： 冲激响应 h(n) 对中心点 (N-1)/2 成奇对称。此时，无论 N 为奇数或偶数，滤波器的相频特性均为线性，并包含有π/2 的固定相移： 因此，信号通过此类滤波器时不仅产生 (N-1)/2 个取样点的延迟，还将产生 90o 的相移，通常这类滤波器又被称为 90o 移相器，并具有很好的应用价值。; 奇对称：θ(ω) 对所有的频率成分都有一个 90°相移。;5.2.2 线性相移 FIR DF 频率响应：Type I;5.2.2 线性相移 FIR DF 频率响应：Type I; 由此可以看出其线性相位特性。由于 cos(nω) 对于 ω=0、π、2π都是偶对称，所以振度响应 Hr(ω) 对ω=0、π、2π也是偶对称。;h(n) 偶对称，N 为偶数（恒相时延、恒群时延〕 由于h(n) 序列的长度为偶数，因此滤波器的频率响应函数可拆分成如下两部分（前后对称部分，中心点处无值）：;5.2.2 线性相移 FIR DF 频率响应：Type II;5.2.2 线性相移 FIR DF 频率响应：Type II;h(n) 奇对称，N 为奇数（恒群时延〕 h(n) 长度为奇数，拆分成前后两部分：;Hr (w);与 c(n) 或 h(n) 的值无关，因此，这种类型的滤波器不适用于低通、带阻或高通滤波器设计，而且 ，这说明 jHr(w) 是纯虚数，对于逼近理想数字希尔伯特变换和微分器，它是很有用的。理想的希尔伯特变换是一个全通滤波器，它对输入信号产生 90 度的相移，它频繁用于通信系统中的调制。微分器广泛用于模拟和数字系统中对信号求导。 2）由于 sin(nω) 对于 ω=0、π、2 π 都是奇对称，所以，Hr(w) 以ω =0、 π 、2π为奇对称。 ;5.2.2 线性相移 FIR DF 频率响应：Type IV;5.2.2 线性相移 FIR DF 频率响应：Type IV;一般形式： ; 一般的 FIR DF 的零、极点：;令：m=N-1-n;上面提到 Zi 肯定是共轭的，故 Zi* 亦必为其零点 于是零点有：;主题概述 ;思路：;设所要求的 DF 的频率响应是 Hd(ejw)，需要注意：它可能是低通、高通、带通和带阻 FIR DF，没有特指某种类型的数字滤波器。 不管是何种 FIR DF，它的频率响应是频域中的周期函数，周期为 2π，所以它可以展开为傅氏级数形式：;因此，所要求的 DF 的系统函数便可求得： 显然，Hd(z) 是非因果的，且 hd(n) 的持续时间为 -∞～ +∞，物理上不可实现。 我们可以采用逼近 Hd(ejw) 的方法 首先把 hd(n) 先截短为有限项，把 hd(n) 截为2M+1项，得：;然后把截短后的 hd(n) 右移，使之变成因果性的序列。 令 H(z) 等于 H1(z) 乘以 z-M 得： 令 h(n)= hd (n-M), n=0, 1, 2, ......, 2M，则 ;对 hd(n) 的截短必然产生误差，即以 |H(ejw)| 近似 |Hd(ejw)| 。; 因为 |H(ejw)| 是对 hd(n) 截短而产生的，假定： ;说明：;将 hd(n) 截短： ;窗函数的频谱： ; 截短，根据时域相乘映射为频域卷积，得：;θ;现在分析几个特殊频率点的滤波器性能： ω= 0 时： 由于一般情况下都满足 ωc 2π / N，因此，H(0) 的值近似等于窗谱函数 WR(ejw) 与θ轴围出的整个面积。 ;5.3.2 窗口法：Gibbs 效应; 理想滤波器的不连续点演化为过渡带;Gibbs现象；;设计FIR DF时，窗函数不仅可以影响过渡带宽度，还能影响肩峰和波动的大小，因此，选择窗函数应使其频谱： 主瓣宽度尽量小，以使过渡带尽量陡。 旁瓣相对于主瓣