数字信号处理原理配套习题答案第8章FIR数字滤波器的设计习题解答.docx

8.1：

FIR低通滤波器的脉冲响应h(n)是由理想低通滤波器的脉冲响应?d

h

当n=

?

对于给定的N=21,滤波器的阶数是20，即n=0,1,2,...,20，中心位置为n=10,理想脉冲响应需要乘以矩形窗,矩形窗函数ω(n)的值为1，因此滤波器的脉冲响应h

频率响应H(ejω

Step1:计算理想低通滤波器的脉冲响应

h

?

矩形窗函数：ω

滤波器的脉冲响应:h(n)=

Step2:计算频率响应

通过傅里叶变换(DTFT)，求解H(

H(

8.2：

与上题类似：

Step1:计算理想低通滤波器的脉冲响应

h

?

三角形窗函数ω

h(n)=

Step2:计算频率响应

通过傅里叶变换(DTFT)，求解H(

H(

8.3:

理想高通滤波器的脉冲响应hd

h

h

Step1:N=51，α

h

h

Step2:汉宁窗的公式为：

ω

h(n)=

8.4:

理想带通滤波器的脉冲响应hd

h

Step1:

h

Step2:

海明窗的表达式为：

ω

h(n)=

8.5：

h

Step1：

其中，M=

h

Step2：

布莱克曼窗

ω

ω

8.6:

理想低通滤波器的频率响应为：

H

其对应的时域冲激响应为：

?

凯泽窗:

ω

其中I0

若ω

8.7：

通常，阻带衰减要求与过渡带宽度一起决定了窗函数的类型：

1、汉明窗：阻带衰减约为44dB

2、凯泽窗：通过调整窗系数β可以满足更大的阻带衰减要求

为满足-45dB的阻带衰减，选择凯泽窗并根据阻带衰减和过渡带宽度来确定

Step1：滤波器长度N计算：

凯泽窗法中，N和过渡带宽度?ω的

N

其中，A是是阻带衰减的绝对值，即A=-45dB，带入得到

N

Step2：理想低通滤波器的冲激响应计算

h

Step3：凯泽窗

ω

8.8：

(1)

?

则

H

由上式可以看出

|

因为?1(n)是偶对称序列，根据线性相位的定义，它是一个线性相位滤波器，其时延为N?12=3.5，

令

要知两个滤波器的性能,必须求出它们各自的频率响应的幅度函数,根据它们的通带起伏以及阻带衰减的情况,来加以比较。由于N=8是偶数，又是线性相位，则

H(

可以令

?

?

及

?

?

代入可得

H

H

从阻带看,H1(ω)的阻带衰减大,而H2(ω)的阻带衰减小,这一点

8.9

由H(z)

z

其中z01和z02、z03和z04互为共轭零点。为了使新系统满足单位抽样响应为实数,必须满足新系统的零点保持共轭分布,因此在移动零点过程中应将z01和z

因此总共可得到的新系统的个数为

M=

总共可得到的系统个数为M+1=16

由系统零点分布可知,系统有3个零点分布在单位圆外,即

z

用全通系统将这些零点映射到单位圆内可得到最小相位系统。

H

(1?0.7

其中最小相位系统为

H

(1?

要获得最大相位系统,需要将原H(z)级联一个全通系统,把单位圆内的零点映射到单位圆外。满足要求的全通系统为

H

从而可得最大相位系统为

H

×

8.10

由于HBP(e

H

且φ(

H

所以带阻滤波器可以表示为：

H

(2)

?

考虑到φ(

φ

?

8.11

因滤波器的冲激响应具有反对称性质，即

h(n)=?h(N?1?n)

故当N为奇数时，有

h(

H(

上式中n用N?12

H(

由于滤波器的频率响应为

H(

令

c(n)=2h(

得到

H(

8.12

(1)

H(

首先计算?min

H

再计算?(n)的幅频响应

H(

进行变量替换,令m=N?1?n,则n=N?1?m,当n=0时,m=N?1,当n=N?1时,m=0,因此幅频响应可以改写为

H(

利用e?j

H(

由于幅频响应只关注幅度，不关注相位，因此两者具有相同的幅频响应

(2)

一个系统是最大相位系统，当其单位冲激响应h(n)的所有非零系数的时延响应都是正的。由于h(n)=?min(N?1?n)。可以看到,h(n)的系数是反向排列的，这意味着系统在频域中的相位响应为负的，因为每一个

8.13

题8.1的MATLAB实现：

%参数设置

N=21;%滤波器阶数

wc=0.5*pi;%截止频率

M=(N-1)/2;%中心位置

%计算理想脉冲响应h_d(n)

n=0:(N-1);

hd=sinc((n-M)*0.5);%0.5=wc/pi

%矩形窗w(n)，所有值为1

w=ones(1,N);

%计算加窗后的脉冲响应h(n)

h=hd.*w;

%计算频率响应H(e^jw)

H=fft(h,1024);%使用FFT计算频率响应

H=