有限冲激响应数字滤波器的设计(4)

2018-11-18 21:17

当N为奇数时(如图所示),只有一个实根z=r,对应于一个一阶网络H0(z)。这时的

H(z)为

H(z)?(1?rzN?N1)N(N?1)/2??H(z)?H(z)?k?0? (7.21)

k?1??显然,N等于奇数时的频率采样修正结构由一个一阶网络结构和(N-1)/2个二阶网络结构组

成。

一般来说,当采样点数N较大时,频率采样结构比较复杂,所需的乘法器和延时器比较多。但在以下两种情况下,使用频率采样结构比较经济。 ?

(1)对于窄带滤波器,其多数采样值H(k)为零,谐振器柜中只剩下几个所需要的谐振器。这时采用频率采样结构比直接型结构所用的乘法器少,当然存储器还是要比直接型用得多一些。 (2)在需要同时使用很多并列的滤波器的情况下,这些并列的滤波器可以采用频率采样结构,并且可以共用梳状滤波器和谐振柜,只要将各谐振器的输出适当加权组合就能组成各个并列的滤波器。

7.2.3.4 快速卷积型结构

根据循环卷积和线性卷积的关系可知,两个长度分别为N1和N2的序列的线性卷积,可以用这两个序列的L(?N1?N2?1)点的循环卷积来实现。

由FIR滤波器的直接型结构:滤波器的输出信号y(n)是输入信号x(n)和滤波器单位冲激响应h(n)的线性卷积。所以,对有限长序列x(n), 我们可以通过补零的方法延长x(n)和

h(n)序列,然后计算它们的循环卷积,从而得到FIR系统的输出y(n)。

实际上,循环卷积通常不在时域上实现,而是利用循环卷积定理,采用DFT或FFT实现有限长序列x(n)和h(n)的线性卷积,如图7.12所示。图中L?N1?N2?1。若采用FFT来实现,L还需满足一定条件,如基2FFT则需L?2,M是整数。显然,由于采用FFT的方法来实现卷积,因此当N1和N2足够长时,处理速度将大大提高,因此,称为快速卷积结构。

对x(n)为无限长的一般情况,可用重叠相加法或重叠保留法实现FIR滤波器的快速卷积结构。

M16

x(n)L 点FFTX(k)X(k)·H(k)H(k)L 点FFTh(n)L 点IFFTy(n)图 7.12 FIR的快速卷积型结构

具体过程如下:

1、将x(n)和h(n)补零变成长为L的序列,L?N1?N2?1,L?2

M?h(n)0?n?N2?1?x(n)0?n?N1?1h(n)?, x(n)???0N?n?L?10N?n?L?1?2?12、将x(n)和h(n)分别做L点DFT(或FFT)得X(k),H(k),k?0,1,?,L?1 3、将X(k)和H(k)相乘得Y(k)?X(k)H(k),k?0,1,?,L?1 4、 对Y(k)做L点IDFT(或IFFT)得y(n),n?0,1,?,L?1

7.2.3.5 线性相位FIR滤波器结构

FIR结构可设计线性相位的数字滤波器,这种滤波器在数据传输、图像处理等系统中是重要的,对于这种形式的滤波器,由于对滤波器系数h(n)有特殊要求,因此结构上也有所不同。详见下节。

17

§7.3线性相位FIR数字滤波器的条件、特点及结构

一、常用FIR滤波器的特点及线性相位条件的引出

在讨论线性相位FIR数字滤波器的条件之前,我们先讨论几种常用FIR滤波器的特点。这些滤波器既包括了一般意义下的高通、低通、带通等数字滤波器,也包括希尔波特(Hilbert)变换器和微分器。

1、低通滤波器的时域特性

理想低通滤波器的幅度特性如图7.13所示

H(j?)?????c0?c?

图7.13 理想低通滤波器的幅度特性

不妨假设相位特性是?(j?)?0,则,理想滤波器的频率特性是

?1??c????c (7.22) H(j?)??其他?0可求得其单位冲激响应为

1hd(n)?2?如图7.14所示。

??e?c?cj?nd?=sin(?cn)/n? (7.23)

图7.14 理想低通滤波器的单位冲激响应(?c??/2)

18

显然,理想低通滤波器的冲激响应为hd(n)?sin(?c0n)/n?,是偶对称的,无限长的。 同样,理想高通、带通等一般滤波器的冲激响应hd(n)也是偶对称的。将hd(n)延迟

??(N?1)/2后,以?为中心对称截取N点而得h(n)。则满足关系式h(n)=h(N-1-n),既是

h(n)关于对称中心??(N?1)/2是偶对称的。

(为什么要位移(N-1)/2 ?请思考。因为FIR的单位冲激响应是有限长的,假设是N点,则FIR系数的一种选取方法是从h(n)中选取N点。而h(n)的较大值一般在n=0附近,故应在n=0附近截取N点。为保证对称性及因果性,则要位移(N-1)/2点)

位移后,滤波器的的频率相应近似为

?e-j??H(j?)???0可见,其相位特性是线性的。

??c????c其他 (7.24)

同理,可以证明具有线性相位特性?(?)????的理想高通、理想带通滤波器的时域相应也满足:h(n)=h(N-1-n)。

2、希尔波特变换器的时域特性

理想希尔波特变换器的频率特性为:

??j0???? (7.25) HD(j?)??j?????0?其幅频特性及相频特性如图7.15所示

φ (jω) π/2 1 |H(jω)| -π -π π

(a)幅频特性 (b)相频特性

图7.15 希尔伯特变换器的频率特性

其冲激响应为

ω -π/2 π ω 1hd(n)?2?

??je?0j?n1d??2???0jej?nd??2sin2[?n?n]2 (7.26)

19

显然它是奇对称的。将其右移??(N?1)/2后,频率特性为

??je?j(N?1?)HD(j?)???j(N?1?)?jesin2[/2/20?????????0 (7.27)

相应的冲激响应为

?(n??)2hd(n)?2]0?n?N?1 (7.28)

?n??以??(N?1)/2为中心,对称截取N点得

h(n)?2sin2[?(n??)2]0?n?N?1 (7.29)

?n??满足关系式h(n)??h(n?1?n),即是奇对称的。相应的频率特性如式(7.27)所示,可见,其相频特性为

???/2?(N?1)?/20?????D(j?)?? (7.30)

?/2?(N?1)?/2?????0?显然是线性的。

3、理想微分器

理想微分器的频率特性如图7.16所示,

|H(jω)| π 0 φ (jω) π/2 -π -π π ω -π

π ω

(a)幅频特性 (b)相频特性

图7.16 理想微分器的频率特性

理想微分器的冲激响应位移??(N?1)/2为后的频率特性为

HD(j?)?j?e?j(N?1)?/2其相位特性是线性的:

?????? (7.31)

20


有限冲激响应数字滤波器的设计(4).doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:读《美国人眼中的朝鲜战争》一书有感

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: