数字信号处理实验IIR滤波器设计.docx
- 文档编号:25633299
- 上传时间:2023-06-10
- 格式:DOCX
- 页数:13
- 大小:188.08KB
数字信号处理实验IIR滤波器设计.docx
《数字信号处理实验IIR滤波器设计.docx》由会员分享,可在线阅读,更多相关《数字信号处理实验IIR滤波器设计.docx(13页珍藏版)》请在冰豆网上搜索。
数字信号处理实验IIR滤波器设计
实验四IIR数字滤波器的设计
(1)
,
,
,
;设计一切比雪夫高通滤波器,观察其通带损耗和阻带衰减是否满足要求。
解:
程序:
clear;
fc=300;fr=200;fs=1000;rp=0.8;rs=20;
wc=2*fs*tan(2*pi*fc/(2*fs));
wt=2*fs*tan(2*pi*fr/(2*fs));
[N,wn]=cheb1ord(wc,wt,rp,rs,'s');
[B,A]=cheby1(N,rp,wn,'high','s');
[bz,az]=bilinear(B,A,fs);
[h,w]=freqz(bz,az);
f=w*fs/(2*pi);
plot(f,20*log10(abs(h)));
axis([0,fs/2,-80,10]);
grid;
xlabel('频率/Hz');
ylabel('幅度/dB');
分析:
f=200Hz时阻带衰减大于30dB,通过修改axis([0,fs/2,-80,10])为axis([200,fs/2,-1,1])
发现通带波动rs满足<0.8。
bz=[0.0262-0.10470.1570-0.10470.0262]
az=[1.00001.52891.65370.94520.2796]
系统函数为:
(2)
,
,
,
;分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字低通滤波器,观察所设计数字滤波器的幅频特性曲线,记录带宽和衰减量,检查是否满足要求。
比较这两种方法的优缺点。
解:
程序:
clear;
fs=1000;fc=200;fr=300;rp=1;rs=25;
%脉冲响应不变法
wp=2*pi*fc;
ws=2*pi*fr;
[N,wn]=buttord(wp,ws,rp,rs,'s');
[b1a1]=butter(N,wn,'s');
[bz1,az1]=impinvar(b1,a1,fs);
[h1,w]=freqz(bz1,az1);
%双线性变换法
wp=2*fs*tan(2*pi*fc/fs/2);
ws=2*fs*tan(2*pi*fr/fs/2);
[N,wn]=buttord(wp,ws,rp,rs,'s');
[b2a2]=butter(N,wn,'s');
[bz2,az2]=bilinear(b2,a2,fs);
[h2,w]=freqz(bz2,az2);
f=w/(2*pi)*fs;
figure;plot(f,abs(h1),'-.r',f,abs(h2),'-b');
grid;xlabel('频率/Hz');ylabel('幅度');
legend('脉冲响应不变法','双线性变换法');
title('巴特沃思低通滤波器,线性幅度谱');
bz1=[0.00000.00020.01530.09950.14440.06110.00750.00020.00000]
az1=[1.0000-1.91992.5324-2.20531.3868-0.63090.2045-0.04500.0060-0.0004]
因此脉冲响应不变法的系统函数为:
bz2=[0.01790.10720.26810.35750.26810.10720.0179]
az2=[1.0000-0.60190.9130-0.29890.1501-0.02080.0025]
因此双线性变换法的系统函数为:
分析:
脉冲响应不变法的N=9,双线性变换法的N=6,由图知它们都满足要求,但脉冲响应的衰减较快,双线性变换的过渡带窄一些,且阶数比脉冲小,容易实现。
(3)利用双线性变换法分别设计满足下列指标的巴特沃思型、切比雪夫型和椭圆型数字低通滤波器,并作图验证设计结果:
,
,
,
。
解:
程序:
clear;
fs=8000;fc=1200;fr=2000;rp=0.5;rs=40;
%巴特沃思低通滤波器
wp=2*fs*tan(2*pi*fc/fs/2);
ws=2*fs*tan(2*pi*fr/fs/2);
[N,wn]=buttord(wp,ws,rp,rs,'s');
[b1a1]=butter(N,wn,'s');
[bz1,az1]=bilinear(b1,a1,fs);
[h1,w]=freqz(bz1,az1);
H1=20*log10(abs(h1));
f=w/(2*pi)*fs;
figure;plot(f,H1);%对数幅度谱
axis([0,fs/2,-100,10]);
grid;xlabel('频率/Hz');ylabel('幅度');
title('巴特沃思低通滤波器,对数幅度谱');
%切比雪夫低通滤波器
wc=2*fs*tan(2*pi*fc/(2*fs));
wt=2*fs*tan(2*pi*fr/(2*fs));
[N,wn]=cheb1ord(wc,wt,rp,rs,'s');
[b2,a2]=cheby1(N,rp,wn,'low','s');
[bz2,az2]=bilinear(b2,a2,fs);
[h2,w]=freqz(bz2,az2);
H2=20*log10(abs(h2));
f=w*fs/(2*pi);
figure;
plot(f,H2);
axis([0,fs/2,-100,10]);
grid;
xlabel('频率/Hz');
ylabel('幅度/dB');
title('切比雪夫低通滤波器,对数幅度谱');
%椭圆型数字低通滤波器
wp=2*fs*tan(2*pi*fc/fs/2);%双线性变换法
ws=2*fs*tan(2*pi*fr/fs/2);
[N,wp]=ellipord(wp,ws,rp,rs,'s');
[b3,a3]=ellip(N,rp,rs,wp,'low','s');
[bz3,az3]=bilinear(b3,a3,fs);
[h3,w]=freqz(bz3,az3);
H3=20*log10(abs(h3));
f=w/(2*pi)*fs;
figure;plot(f,H3);
axis([0,fs/2,-100,10]);
grid;xlabel('频率/Hz');ylabel('幅度/dB');
title('椭圆型数字低通滤波器,对数幅度谱');
巴特沃思低通滤波器:
bz1=[0.00040.00320.01290.03020.04530.04530.03020.01290.00320.0004]
az1=[1.0000-2.79964.4582-4.54123.2404-1.63300.5780-0.13700.0197-0.0013]
系统函数为:
分析:
N=9,为九阶巴特沃思低通滤波器,从图中可以看出通带波动和阻带衰减都满足设计要求。
切比雪夫低通滤波器:
bz2=[0.00260.01320.02640.02640.01320.0026]
az2=[1.0000-2.97754.2932-3.51241.6145-0.3334]
系统函数为:
分析:
N=5,为五阶切比雪夫低通滤波器,从图中可以看出通带波动和阻带衰减都满足设计要求。
椭圆型数字低通滤波器:
bz=[0.03890.03630.06650.03630.0389]
az=[1.0000-2.14442.3658-1.32500.3332]
系统函数为:
分析:
N=4,为四阶椭圆型数字低通滤波器,从图中可以看出通带波动和阻带衰减都满足设计要求。
(4)分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字带通滤波器,已知
,其等效的模拟滤波器指标为
,
,
,
,
,
。
解:
程序:
clear;
fs=30000;fc=[2000,3000];
fr=[1500,6000];rp=3;rs=20;
%脉冲响应不变法
wp=2*pi*fc;
ws=2*pi*fr;
[N,wn]=buttord(wp,ws,rp,rs,'s');
[b1,a1]=butter(N,wn,'s');
[bz1,az1]=impinvar(b1,a1,fs);
[h1,w]=freqz(bz1,az1);
%双线性变换法
wp=2*fs*tan(2*pi*fc/fs/2);
ws=2*fs*tan(2*pi*fr/fs/2);
[N,wn]=buttord(wp,ws,rp,rs,'s');
[b2,a2]=butter(N,wn,'s');
[bz2,az2]=bilinear(b2,a2,fs);
[h2,w]=freqz(bz2,az2);
f=w/(2*pi)*fs;
figure;
plot(f,abs(h1),'-.r',f,abs(h2),'-b');
grid;
xlabel('频率/Hz');
ylabel('幅度');
legend('脉冲响应不变法','双线性变换法');
bz1=[-0.00000.0057-0.01220.00250.0089-0.00490]
az1=[1.0000-4.805610.2376-12.26258.7012-3.47190.6145]
因此脉冲响应不变法的系统函数为:
bz2=[0.0014-0.0000-0.0042-0.00000.0042-0.0000-0.0014]
az2=[1.0000-4.807110.2473-12.28388.7245-3.48490.6176]
因此双线性变换法的系统函数为:
(5)利用双线性变换法设计满足下列指标的切比雪夫型数字带阻滤波器,并作图验证设计结果:
当
时,
,当
以及
时,
;采样频率
。
解:
程序:
clear;
fs=10000;fc=[5003000];
fr=[10002000];rp=3;rs=20;
wp=2*fs*tan(2*pi*fc/fs/2);
ws=2*fs*tan(2*pi*fr/fs/2);
[N,wn]=cheb1ord(wp,ws,rp,rs,'s');
[b,a]=cheby1(N,rp,wp,‘stop’,'s');
[bz,az]=bilinear(b,a,fs);
[h,w]=freqz(bz,az);f=w/(2*pi)*fs;
figure;plot(f,20*log10(abs(h)));
axis([0,2000,-200,10]);
grid;xlabel('频率/Hz');ylabel('幅度/dB');
title('切比雪夫带阻滤波器');
bz=[0.2192-0.54180.7732-0.54180.2192]
az=[1.0000-1.10510.2612-0.42560.4502]
因此数字滤波器的系统函数为:
分析:
N=4,为四阶切比雪夫带阻滤波器,由图知满足设计要求。
Welcome!
!
!
欢迎您的下载,
资料仅供参考!
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 实验 IIR 滤波器 设计