数字信号处理实验IIR滤波器设计说明Word格式文档下载.docx

上传人:b****6 文档编号:20450739 上传时间:2023-01-23 格式:DOCX 页数:13 大小:188.47KB
下载 相关 举报
数字信号处理实验IIR滤波器设计说明Word格式文档下载.docx_第1页
第1页 / 共13页
数字信号处理实验IIR滤波器设计说明Word格式文档下载.docx_第2页
第2页 / 共13页
数字信号处理实验IIR滤波器设计说明Word格式文档下载.docx_第3页
第3页 / 共13页
数字信号处理实验IIR滤波器设计说明Word格式文档下载.docx_第4页
第4页 / 共13页
数字信号处理实验IIR滤波器设计说明Word格式文档下载.docx_第5页
第5页 / 共13页
点击查看更多>>
下载资源
资源描述

数字信号处理实验IIR滤波器设计说明Word格式文档下载.docx

《数字信号处理实验IIR滤波器设计说明Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《数字信号处理实验IIR滤波器设计说明Word格式文档下载.docx(13页珍藏版)》请在冰豆网上搜索。

数字信号处理实验IIR滤波器设计说明Word格式文档下载.docx

wt=2*fs*tan(2*pi*fr/(2*fs));

[N,wn]=cheb1ord(wc,wt,rp,rs,'

s'

);

[B,A]=cheby1(N,rp,wn,'

high'

'

[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)

分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字低通滤波器,观察所设计数字滤波器的幅频特性曲线,记录带宽和衰减量,检查是否满足要求。

比较这两种方法的优缺点。

fc=200;

fr=300;

rp=1;

rs=25;

%脉冲响应不变法

wp=2*pi*fc;

ws=2*pi*fr;

[N,wn]=buttord(wp,ws,rp,rs,'

[b1a1]=butter(N,wn,'

[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);

[b2a2]=butter(N,wn,'

[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'

xlabel('

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)利用双线性变换法分别设计满足下列指标的巴特沃思型、切比雪夫型和椭圆型数字低通滤波器,并作图验证设计结果:

fs=8000;

fc=1200;

fr=2000;

rp=0.5;

rs=40;

%巴特沃思低通滤波器

[bz1,az1]=bilinear(b1,a1,fs);

H1=20*log10(abs(h1));

plot(f,H1);

%对数幅度谱

axis([0,fs/2,-100,10]);

巴特沃思低通滤波器,对数幅度谱'

%切比雪夫低通滤波器

[b2,a2]=cheby1(N,rp,wn,'

low'

H2=20*log10(abs(h2));

plot(f,H2);

切比雪夫低通滤波器,对数幅度谱'

%椭圆型数字低通滤波器

%双线性变换法

[N,wp]=ellipord(wp,ws,rp,rs,'

[b3,a3]=ellip(N,rp,rs,wp,'

[bz3,az3]=bilinear(b3,a3,fs);

[h3,w]=freqz(bz3,az3);

H3=20*log10(abs(h3));

plot(f,H3);

椭圆型数字低通滤波器,对数幅度谱'

巴特沃思低通滤波器:

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)分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字带通滤波器,已知

,其等效的模拟滤波器指标为

fs=30000;

fc=[2000,3000];

fr=[1500,6000];

rp=3;

[N,wn]=buttord(wp,ws,rp,rs,'

[b1,a1]=butter(N,wn,'

[b2,a2]=butter(N,wn,'

plot(f,abs(h1),'

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)利用双线性变换法设计满足下列指标的切比雪夫型数字带阻滤波器,并作图验证设计结果:

时,

,当

以及

采样频率

解:

程序:

fs=10000;

fc=[5003000];

fr=[10002000];

[N,wn]=cheb1ord(wp,ws,rp,rs,'

[b,a]=cheby1(N,rp,wp,‘stop’,'

[bz,az]=bilinear(b,a,fs);

f=w/(2*pi)*fs;

axis([0,2000,-200,10]);

切比雪夫带阻滤波器'

 

bz=[0.2192-0.54180.7732-0.54180.2192]

az=[1.0000-1.10510.2612-0.42560.4502]

因此数字滤波器的系统函数为:

N=4,为四阶切比雪夫带阻滤波器,由图知满足设计要求。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 总结汇报 > 其它

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1