数字信号处理第七章.doc
- 文档编号:1644441
- 上传时间:2022-10-23
- 格式:DOC
- 页数:27
- 大小:2.23MB
数字信号处理第七章.doc
《数字信号处理第七章.doc》由会员分享,可在线阅读,更多相关《数字信号处理第七章.doc(27页珍藏版)》请在冰豆网上搜索。
第七章数字滤波器设计
7.1:
无限冲激响应滤波器的阶数的估计
Q7.1用MATTAB确定一个数字无限冲激响应低通滤波器所有四种类型的最低阶数。
指标如下:
40kHz的抽样率,,4kHz的通带边界频率,8kHz的阻带边界频率,0.5dB的通带波纹,40dB的最小阻带衰减。
评论你的结果。
答:
标准通带边缘角频率Wp是:
标准阻带边缘角频率Ws是:
理想通带波纹Rp是0.5dB
理想阻带波纹Rs是40dB
1.使用这些值得到巴特沃斯低通滤波器最低阶数N=8,相应的标准通带边缘频率Wn是0.2469.
2.使用这些值得到切比雪夫1型低通滤波器最低阶数N=5,相应的标准通带边缘频率Wn是0.2000.
3/使用这些值得到切比雪夫2型低通滤波器最低阶数N=5,相应的标准通带边缘频率Wn是0.4000.
4.使用这些值得到椭圆低通滤波器最低阶数N=8,相应的标准通带边缘频率Wn是0.2000.
从以上结果中观察到椭圆滤波器的阶数最低,并且符合要求。
Q7.2用MATLAB确定一个数字无限冲激响应高通滤波器所有四种类型的最低阶数。
指标如下:
3500Hz的抽样率,1050Hz的通带边界频率,600Hz的阻带边界频率,1dB的通带波纹,50dB的最小阻带衰减。
评论你的结果
答:
标准通带边缘角频率Wp是:
标准阻带边缘角频率Ws是:
理想通带波纹Rp是1dB
理想阻带波纹Rs是50dB
1.使用这些值得到巴特沃斯高通滤波器最低阶数N=8,相应的标准通带边缘频率Wn是0.5646.
2.使用这些值得到切比雪夫1型高通滤波器最低阶数N=5,相应的标准通带边缘频率Wn是0.6000.
3.使用这些值得到切比雪夫2型高通滤波器最低阶数N=5,相应的标准通带边缘频率Wn是0.3429.
4.使用这些值得到椭圆低通滤波器最低阶数N=4,相应的标准通带边缘频率Wn是0.6000.
从以上结果中观察到椭圆滤波器的阶数最低,并且符合要求。
Q7.3用MATLAB确定一个数字无限冲激响应带通滤波器所有四种类型的最低阶数。
指标如下:
7kHz的抽样率,1.4kHz和2.1kHz的通带边界频率,1.05kHz和2.45kHz的阻带边界频率,,0.4dB的通带波纹,50dB的最小阻带衰减。
评论你的结果。
答:
标准通带边缘角频率Wp是:
标准阻带边缘角频率Ws是:
理想通带波纹Rp是0.4dB
理想阻带波纹Rs是50dB
1.使用这些值得到巴特沃斯带通滤波器最低阶数2N=18,相应的标准通带边缘频率Wn是[0.38350.6165].
2.使用这些值得到切比雪夫1型带通滤波器最低阶数2N=12,相应的标准通带边缘频率Wn是[0.40000.6000].
3.使用这些值得到切比雪夫2型带通滤波器最低阶数2N=12,相应的标准通带边缘频率Wn是[0.30000.7000].
4.使用这些值得到椭圆带通滤波器最低阶数2N=8,相应的标准通带边缘频率Wn是[0.40000.6000].
从以上结果中观察到椭圆滤波器的阶数最低,并且符合要求。
Q7.4用MATLAB确定一个数字无限冲激响应带阻滤波器所有四种类型的最低阶数。
指标如下:
12kHz的抽样率,2.1kHz和4.5kHz的通带边界频率,2.7kHz和3.9kHz的阻带边界频率,0.6dB的通带波纹,45dB的最小阻带衰减。
评论你的结果。
答:
标准通带边缘角频率Wp是:
标准阻带边缘角频率Ws是:
理想通带波纹Rp是0.6dB
理想阻带波纹Rs是45dB
1.使用这些值得到巴特沃斯带阻滤波器最低阶数2N=18,相应的标准通带边缘频率Wn是[0.38730.7123].
2.使用这些值得到切比雪夫1型带阻滤波器最低阶数2N=10,相应的标准通带边缘频率Wn是[0.35000.7500].
3.使用这些值得到切比雪夫2型带阻滤波器最低阶数2N=10,相应的标准通带边缘频率Wn是[0.45000.6500].
4.使用这些值得到椭圆带阻滤波器最低阶数2N=8,相应的标准通带边缘频率Wn是[0.35000.7500].
从以上结果中观察到椭圆滤波器的阶数最低,并且符合要求。
7.2:
无限冲激响应滤波器设计
程序P7.1说明巴特沃斯带阻滤波器的设计。
%巴特沃斯带阻滤波器的设计
Ws=[0.40.6];Wp=[0.20.8];Rp=0.4;Rs=50;
%估计滤波器的阶数
[N1,Wn1]=buttord(Wp,Ws,Rp,Rs);
%设计滤波器
[num,den]=butter(N1,Wn1,'stop');
%显示传输函数
disp('分子系数是');disp(num);
disp('分母系数是');disp(den);
%计算增益响应
[g,w]=gain(num,den);
%绘制增益响应
plot(w/pi,g);grid
axis([01-605]);
xlabel('\omega/\pi');ylabel('增益,dB');
title('巴特沃斯带阻滤波器的增益响应');
Q7.5通过运行程序P7.1来设计巴特沃兹带阻滤波器。
写出所产生的传输函数的准确表达式。
滤波器的指标是什么,你的设计符合指标吗,使用MATLAB,计算并绘制滤波器的未畸变的相位响应及群延迟响应。
答:
表达式是:
滤波器参数是:
Wp1=0.2π,Ws1=0.4π,Ws2=0.6π,Wp2=0.8π,Rp=0.4dB,Rs=50dB.
设计的滤波器增益响应如下:
从图中可以总结出设计符合指标。
滤波器的未畸变的相位响应及群延迟响应如下:
Q7.6修改程序P7.1来设计符合习题Q7.1所给指标的切比雪夫1型低通滤波器。
写出所产生的传输函数的准确表达式。
你的设计符合指标吗?
使用MATLAB,计算并绘制滤波器的未畸变的相位响应及群延迟响应。
答:
表达式如下:
设计的滤波器增益响应如下:
从图中可以总结出设计符合指标。
滤波器的未畸变的相位响应及群延迟响应如下:
Q7.7修改程序P7.1来设计符合习题Q7.2所给指标的切比雪夫2型高通滤波器。
写出所产生的传输函数的准确表达式。
你的设计符合指标吗?
使用MATLAB,计算并绘制滤波器的未畸变的相位响应及群延迟响应。
答:
表达式如下:
设计的滤波器增益响应如下:
从图中可以总结出设计符合指标。
滤波器的未畸变的相位响应及群延迟响应如下:
Q7.8修改程序P7.1来设计符合习题Q7.3所给指标的椭圆带通滤波器。
写出所产生的传输函数的准确表达式。
你的设计符合指标吗,使用MATLAB,计算井绘制滤波器的未畸变的相位响应及群延迟响应。
答:
表达式如下:
设计的滤波器增益响应如下:
从图中可以总结出设计符合指标。
滤波器的未畸变的相位响应及群延迟响应如下:
7.3:
吉布斯现象
Q7.9使用函数sinc编写一个MATLAB程序,以产生截止频率在Wc=0.4π处、长度分别为81,61,41和21的四个零相位低通滤波器的冲激响应系数,然后计算并画出它们的幅度响应。
使用冒号“:
”运算符从长度为81的滤波器的冲激响应系数中抽出较短长度滤波器的冲激响应系数。
在每一个滤波器的截止频率两边研究频率响应的摆动行为。
波纹的数量与滤波器的长度之间有什么关系?
最大波纹的高度与滤波器的长度之间有什么关系?
你将怎样修改上述程序以产生一个偶数长度的零相位低通滤波器的冲激响应系数?
答:
长度为81时幅度响应如下:
长度分别为61,41和21的幅度响应如下:
从中可以观察到由于吉布斯现象产生的幅度响应的摆动行为。
波纹的数量与滤波器的长度之间的关系——波纹的数量减少与长度成正比。
最大波纹的高度与滤波器的长度之间的关系——最大波纹的高度与长度无关。
Q7.10使用函数sinc编写一个MATLAB程序,以产生一个截止频率在Wc=0.4π处、长度为45的零相位高通滤波器的冲激响应系数,计算并画出其幅度响应。
在每一个滤波器的截止频率两边研究频率响应的摆动行为。
你将怎样修改上述程序以产生一个偶数长度的零相位高通滤波器的冲激响应系数?
答:
长度为45时幅度响应如下:
从中可以观察到由于吉布斯现象产生的幅度响应摆动行为。
在这种情况下你不能改变长度。
原因:
这是一个零相位滤波器,这意味着它也是一个线性相位滤波器,因为零相是一种特殊的线性相位的子集。
现在,理想的有限脉冲响应长度甚至有对称的中点h[n]。
使其成了一个线性相位FIR滤波器。
二型滤波器不可能是高通滤波器,因为必须在z=-1处有零点,意味着w=+-π。
Q7.11编写一个MATLAB程序,以产生长度分别为81,61,41和21的四个零相位微分器的冲激响应系数,计算并画出它们的幅度响应。
下面的代码段显示了怎样产生一个长度为2M+1的微分器。
n=1:
M;
b=cos(pi*n)./n;
num=[-fliplr(b)0b];
对于每种情况,研究微分器的频率响应的摆动行为。
波纹的数量与微分器的长度之间有什么关系,最大波纹的高度与滤波器的长度之间有什么关系?
答:
幅度响应分别如下:
从中可以观察到由于吉布斯现象产生的幅度响应的摆动行为。
波纹的数量与微分器的长度之间的关系——两者成正比。
最大波纹的高度与滤波器的长度之间的关系——两者间没有关系。
Q7.12编写一个MA11AB程序,以产生长度分别为81,61.41和21的四个离散时间希尔伯特变换器的冲激响应系数,计算并画出它们的幅度响应。
下面的代码段显示了怎样产生一个长度为2M十1的希尔伯特变换器。
n=1:
M;
c=sin(pi*n)./2;
b=2*(c.*c)./(pi*n);
num=[-fliplr(b)0b];
对于每种情况,研究希尔伯特变换器的频率响应的摆动行为。
波纹的数量与希尔伯特变换器的长度之间有什么关系?
最大波纹的高度与滤波器的长度之间有什么关系?
答:
幅度响应如下:
从中可以观察到由于吉布斯现像产生的幅度响应的摆动行为。
波纹的数量与希尔伯特变换器的长度之间的关系——两者成正比。
最大波纹的高度与滤波器的长度之间的关系——两者无关系。
7.4:
有限冲激响应滤波器的阶数估计
Q7.13线性相位低通FIR滤波器的阶数估算,参数如下:
wp=2kHz,ws=2.5kHz,dp=0.005,ds=0.005,FT=10kHz
使用kaiord的结果为N=46
使用ceil命令的目的是朝正方向最接近整数方向取整。
使用nargin命令的目的是表明函数M文件体内变量的数目。
Q7.14(a)线性相位FIR滤波器的阶数估算,其中采样频率改为FT=20kHz,则结果为N=91。
(b)线性相位FIR滤波器阶数的估计,其中通带波纹改成dp=0.002和ds=0.002结果为N=57。
(c)线性相位FIR滤波器的阶数估算,其中阻带宽度改成ws=2.3kHz,结果为N=76.
从上述结果和7.13的对比我们可以观察到:
滤波器阶数和采样频率的关系为–对于一个给定的模拟过渡带宽,采样频率的增加导致估算阶数也相应增加,朝下一个整数取整。
其中模拟过渡带宽|Fp-Fs|和Δω的关系:
Δω=2pi*|Fp-Fs|/FT。
因此增加FT会减小Δω。
滤波器阶数和通带波纹宽度的关系为估计的阶数大致和log(底数为10)成比例的扩散。
滤波器阶数和过渡带宽度的关系为在舍入的时候,阶数随着过渡带宽成比例的改变。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 第七