数字信号处理08.docVIP

第八章（有限冲击响应）FIR DF设计 一个离散时间系统H(z)=B(z)/A(z)，若分母多项式A(z)的系数a1=…=aN=0，那么该系统即变成FIR系统： 显然，系数b0，b1，…，bM即是该系统的单位抽样响应h(0)，h(1)，…，h(M)，且当nM时，h(n)≡0。 由于FIR系统只有零点，因此不象IIR系统那样容易取得较好的通带与带阻衰减特性。要取得较好的衰减特性，一般要求H(z)的阶次要高，即M要大。但FIR也有自己突出的优点，系统始终是稳定的，既易于实现线性相位，又允许设计多通带、多阻带的滤波器。而这后两项都是IIR系统不易实现的。 IIR DF设计中借助于模拟滤波器、面向极点系统的设计方法并不适用于仅包含零点的FIR系统，因此，FIR DF的设计主要建立在对理想滤波器频率特性的某种近似基础上，这些近似方法包括窗函数法、频率抽样法和最佳一致逼近法等。 §8.1 FIR DF的线性相位特征 离散时间系统的频率响应包括幅频响应和相频响应，幅频响应反映信号x(n)通过该系统后各频率成分衰减的情况，而相频响应反映信号x(n)中各成分通过该系统后在时间上发生移位的情况。一个理想的离散时间系统，除了具有所希望的幅频特性外，最好还能具有线性相位，这在诸如语音合成、波形传输等应用领域会带来极大的便利。 设一个离散时间系统的幅频特性为1，而相频特性具有如下线性相位： arg[H(ejω)]=-kω （3-1-4） 其中k为常数，该式表明系统的相移与频率成正比。则当信号x(n)通过该系统后，其输出y(n)的频率特性为： 所以，y(n)= x(n-k)，这样，输出y(n)等于输入在时间上的移位，达到了无失真传输的目的。 例：若x(n)=cos(ω0n)+ cos(2ω0n)，令H(ejω)= e-jkω，那么x(n)通过该系统后输出： y(n)= cos[ω0(n-k)]+ cos[2ω0(n-k)]， x(n)，y(n)分别如图所示，二者仅在时间轴上移了k个抽样周期。若再次令|H(ejω)|=1，而令：；则输出y(n)为： y(n)= cos(ω0n-π/4)+ cos(2ω0n-π)，如图c所示，波形明显发生了失真。 若arg[H(ejω)]==kω+β （3-1-5）也称该相频特性为线性相位。 下面将要证明当h(n)=± h(N-1-n)时，FIR滤波器具有线性相位，因为N可以取奇数或偶数，所以共有四种情况： 第一类FIR滤波器，h(n)= h(N-1-n)，即h(n)偶对称 当N为奇数时， ，经推导， ； 令：；（3-1-5） 则：， 因此，H(ejω)具有线性相位，即：arg[H(ejω)]=- (N-1)ω/2 （3-1-7） 其幅频特性具有低通特性。 当N为偶数时， ； 经推导，，（3-1-8） 其中：b(n)=2h(N/2-n) （3-1-9） 相频特性不变，且幅频特性也是低通的。 第二类FIR滤波器，h(n)=- h(N-1-n)，即h(n)奇对称 当N为奇数时，这时因为h(n)以中心(N-1)/2为奇对称，所以必有：h[(N-1)/2]=0，经推导：；式中：c(n)=2h[(N-1)/2-n]， n=1，2，…，(N-1)/2；（3-1-10） 相频特性为：arg[H(ejω)]=- (N-1)ω/2+π/2 （3-1-11） 当N为偶数时， ；式中：d(n)=2h(N/2-n)， n=1，2，…，N/2；（3-1-12）相频特性不变。 由以上讨论可知，当FIR DF的抽样响应满足奇偶对称性时，滤波器具有线性相位。当h(n)奇对称时，通过该滤波器的所有频率成分还会发生90°的相移，这将该信号先通过一个90°的相移器，然后再做滤波，由于这时它的幅频特性是正弦函数的组合，因此第二类FIR滤波器的幅频特性近似于差分器的幅频特性，因此我们设计一般用途的滤波器时，h(n)多取偶对称，而长度N则多取为奇数。 全通滤波器与最小相位滤波器 对于IIR系统，单位抽样响应为无限长，很难满足线性相位的条件，h(n)=± h(N-1-n)，所以在实际工作中，经常用一个全通滤波器与已经设计好的IIR DF相级联，在不改变幅频响应的情况下，对相位做某种校正，尽可能去近似线性相位。 全通滤波器 如果一个因果系统的幅频响应对于所有频率都等于1，或其它某个常数，即有： ；则称该系统Hap(z)为全通滤波器。例如： ，是一个一阶全通滤波器，它的极点在z=λ处，而零点在z=1/λ处，极点与零点是以单位圆为镜像对称的。可以证明，其幅频响应和相频响应分别为： ； 一个二阶全通滤波器的转移函数为： ；式中λ为复数，且|λ|＜1，显然它存在一对单位圆内的共轭极点，一对共轭零点与极点以单位圆为镜像对称。 一般高阶全通滤波器的转移函数