有限冲激响应数字滤波器设计.doc
文本预览下载声明
实验5 有限冲激响应数字滤波器设计
一、实验目的：
1、加深对数字滤波器的常用指标理解。
2、学习数字滤波器的设计方法。
二、实验原理：
低通滤波器的常用指标：
（1）通带边缘频率；
（2）阻带边缘频率；
（3）通带起伏；
（4）通带峰值起伏，
（5）阻带起伏，最小阻带衰减。
三、预习要求
在MATLAB中，熟悉函数fir1、kaiserord 、remezord、remez的使用； B = fir1(n,Wn,high,noscale)设计滤波器； [n,Wn,beta,ftype] = kaiserord(f,a,dev)估计滤波器数； [n,fo,ao,w] = remezord (f,a,dev,fs)计算等波纹滤波器阶数n和加权函数w(ω)； B=remez(n,f,a)进行等波纹滤波器的设计。
阅读扩展练习中的实例，学习FIR滤波器的设计方法及其在MATLAB中的实现；
给出FIR数字滤波器的冲激响应，绘出它们的幅度和相位频响曲线，讨论它们各自的实现形式和特点。
数字滤波器有IIR和FIR两种类型，它们的特点和设计方法不同。
四、实验内容：
利用MATLAB编程,分别用窗函数法和等波纹滤波器法设计两种FIR数字滤波器，指标要求如下：
通带边缘频率：，通带峰值起伏：。阻带边缘频率：，最小阻带衰减：。 N0=ceil(4.46*pi/Bt(2));
B=fir1(N0-1,wc,bandpass,kaiser(N0,3.384));
figure(5)
freqz(B,1)
title(kaiser window)
figure(6)
stem(B)
title(kaiser window)
（2）用等波纹法设计：
调用函数[n,fpts,mag,wt]=remezord(f，a，dev)
f=[0.3 0.45 0.65 0.8]
a=[0 1 0]
dev=[0.01 0.1087 0.01]
其含义同函数[n,wn,bta,ftype]=kaiserord（f，a，dev，fs）中的参数相同。
f=[0.3,0.45,0.65,0.75];
m=[0 1 0];
rp=1;rs=40;
dat1=(10^(rp/20)-1)/(10^(rp/20)+1);
dat2=10^(-rs/20);
rip=[dat1,dat2,dat1];
[M,fo,mo,w]=remezord(f,m,rip);M=M+1;
hn=remez(M,fo,mo,w);
figure(9)
freqz(hn,1)
title(remez)
figure(10)
stem(hn)
title(remez)
实验总结
本实验主要是加深对数字滤波器的常用指标理解学习数字滤波器的设计方法。在MATLAB中，函B = fir1(n,Wn,high,noscale)来设计滤波器； [n,Wn,beta,ftype] = kaiserord(f,a,dev)估计滤波器数； [n,fo,ao,w] = remezord (f,a,dev,fs)计算等波纹滤波器阶数n和加权函数w(ω)； B=remez(n,f,a)进行等波纹滤波器的设计。
Personal loans in these measures refers to eligible natural persons paid for personal consumption, production and operation of legitimate uses, such as loans. 3.2 abbreviations No 3.3 classification No 4. responsibility and authority departments/sectors/positions explain duties incompatible loan investigate gang
Personal loans in these measures refers to eligible natural persons paid for personal consumption, production and operation of legitimate uses, such as loans. 3.2 abbreviations No 3.3 classification No 4. responsibility and authority departments/sectors/positions explain duties incompatible loan investigate g
显示全部