资源描述:
《实验指导书 7-8章》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、实验六FIR滤波器设计一.实验目的1.学习FIR数字滤波器窗口函数设计法.2.学习FIR数字滤波器频率抽样设计法3.进行FIR、IIR滤波器性能比较.二.实验内容1.设计一个线性相位滤波器,其理想特性如下图所示。通带内幅度为1,阻带内幅度为0。数字截止频率为,求。理想低通特性图2.设计如下图所示理想特性FIR线性相位数字滤波器,求。理想高通特性图3.用频率抽样法设计一个FIR数字低通滤波器,截至频率=,取。三.实验说明1.设数字截止频率为=,的长度N=21,求解实验内容1的MATLAB实现程序如下:--------------------------------------------
2、-------------------------------------------------------------------clear;%清空所有数closeall;%关闭所有窗口N=21;%N定义为21wc=pi/4;n=0:N-1;%定义N的范围且步长为一r=(N-1)/2;hdn=sin(wc*(n-r))/pi./(n-r);%理想滤波器ifrem(N,2)~=0%如果N为奇数,处理不定式hdn(r+1)=wc/pi;endwn1=boxcar(N);%矩形窗hn1=hdn.*wn1';%为得到的序列加窗wn2=hamming(N);%hamming窗hn2=hdn.
3、*wn2';%为得到的序列加窗subplot(2,2,1)%画2*2个图第一个图为hn2stem(n,hn1,'.')%绘制hn1波形line([0,20],[0,0]);%画线title('矩形窗设计的h(n)');%设图的标号名称xlabel('n');ylabel('h(n)');subplot(2,2,3)stem(n,hn2,'.')%绘制单位脉冲响应h(n)波形line([0,20],[0,0]);title('hamming窗设计的h(n)');xlabel('n');ylabel('h(n)');hn11=fft(hn1,512);%进行离散傅里叶变换w=2*[0:51
4、1]/512;subplot(2,2,2)plot(w,20*log10(abs(hn11)))grid;%添加网格线axis([0,2,-80,5]);%坐标范围title('幅度特性');xlabel('w/pi');ylabel('幅度(dB)');hn22=fft(hn2,512);subplot(2,2,4)plot(w,20*log10(abs(hn22)))grid;%添加网格线axis([0,2,-80,5]);title('幅度特性');xlabel('w/pi');ylabel('幅度(dB)');----------------------------------
5、-----------------------------------------------------------------------------程序运行结果如下图所示。用矩形窗和hamming窗设计的FIR低通滤波器