用窗函数法设计FIR数字滤波器_第1页
用窗函数法设计FIR数字滤波器_第2页
用窗函数法设计FIR数字滤波器_第3页
用窗函数法设计FIR数字滤波器_第4页
用窗函数法设计FIR数字滤波器_第5页
已阅读5页,还剩3页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

实验五用窗函数法设计FIR数字滤波器一、实验目的:掌握用窗函数法设计FIR数字滤波器的原理和方法熟悉线性相位FIR数字滤波器特性。了解各种窗函数对滤波特性的影响。二、实验原理线性相位特点在实际应用中非常重要,如在数据通信、图像处理、语音信号处理等领域,往往要求系统具有线性相位特性,因而常采用容易设计成线性相位的有限冲激响应FIR数字滤波器来实现。常用窗函数:矩形窗(5.21)Hann(汉纳)窗(5.22)Hamming(汉明)窗(5.23)Blackman(布莱克曼)窗(5.24)Kaiser(凯泽)窗(5.25)其中下面介绍用窗函数设计FIR滤波器的步骤:根据技术要求确定待求滤波器的单位取样响应。根据对过渡带和阻带衰减的要求,选择窗函数的形式,并估计窗函数长度,决定于窗口的形式,表示滤波器的过渡带。利用选好的窗函数计算滤波器的单位取样响应。验算技术指标是否满足要求。设计出来的滤波器的频率响应用来计算。Matlab数字信号处理工具箱中常用的FIR数字滤波器设计函数hanning汉纳窗函数。hamming汉明窗函数。blackman布莱克曼窗函数。kaiser凯泽窗函数。kaiserord凯泽窗函数设计法的阶数估计。fir1窗函数法FIR数字滤波器设计:低通、高通、带通、带阻、多频率滤波。fir2频率抽样法FIR数字滤波器设计:任意频率响应。三、实验内容分别用矩形窗、hann窗,hamming窗设计一个阶数为9的FIR低通数字滤波器,截止频率。画出各种方法设计的数字滤波器的单位脉冲响应。画出它们的幅频响应,并比较各滤波器的通带纹波和阻带纹波,有何结论?若当输入,计算各滤波器的输出并画出其波形。程序过程:clc;clearall;%阶数为M=9,数字截止频率为pi/3;设计II型低通线性相位滤波器M=9;Wc=pi/3;k=0:M;hd=Wc*sinc(Wc*(k-0.5*M))/pi;xk=1+2*cos(pi*k/4)+cos(pi*k/2);figure(1);%以下是矩形窗截断wk=ones(1,M+1);hk=hd.*wk;[H,w]=freqz(hk,1);subplot(311);stem(k,hk,'.');title('矩形窗截断的单位脉冲响应');%以下是hann窗截断wk=hanning(M+1);hk=hd.*wk';[H,w]=freqz(hk,1);subplot(312);stem(k,hk,'.');title('hanniing窗截断的单位脉冲响应');%以下是hamming窗截断wk=hamming(M+1);hk=hd.*wk';[H,w]=freqz(hk,1);subplot(313);stem(k,hk,'.');Wp=0.4*pi;Ws=0.6*pi;Ap=0.5;As=45;Wc=(Wp+Ws)/2;%Blackman窗的近似过渡带宽度为11.4pi/N;窗函数的长度NN=ceil(11.4*pi/(Ws-Wp));%N=58,滤波器阶次M=N-1=57可以设计II型低通线性相位系统M=N-1;k=0:M;hd=Wc*sinc(Wc*(k-0.5*M))/pi;wk=blackman(N);hk=hd.*wk';[H,w]=freqz(hk,1);subplot(211);plot(w/pi,20*log10(abs(H)));grid;xlabel('Normalizedfrequency');ylabel('GainindB');title('blackman窗设计的FIR滤波器');%kaiser窗设计subplot(212);f=[Wp/pi,Ws/pi];a=[1,0];dev=[1-10^(-0.05*Ap),10^(-0.05*As)];[M1,Wc1,beta,ftype]=kaiserord(f,a,dev);wk1=kaiser(M1+1,beta);hk1=fir1(M1,Wc1,ftype,wk1);[H1,w1]=freqz(hk1,1);plot(w1/pi,20*log10(abs(H1)));grid;xlabel('Normalizedfrequency');ylabel('GainindB');title('kaiser窗设计的FIR滤波器');比较:kaiser窗的过渡带较长,在阻带的衰减波动逐渐减小;利用blackman窗设计出的低通滤波器阻带衰减最大。用频率取样法设计一个的Ⅰ型线性相位带通FIR滤波器。带通滤波器的通带截止频率分别为。%用频率取样法设计一个M=44的Ⅰ型线性相位带通FIR滤波器。%带通滤波器的通带截止频率分别为clc;clearall;Wp1=0.3*pi;Wp2=0.5*pi;M=44;m=0:M/2;Wm=2*pi.*m/(M+1);%设计理想滤波器的幅度函数Ad[m]mtr1=floor(Wp2*(M+1)/(2*pi))+2;Ad1=double([Wm<=Wp2]);mtr2=ceil(Wp1*(M+1)/(2*pi));Ad2=double([Wp1<=Wm]);Ad=Ad1.*Ad2;Ad(mtr1)=0.38;Ad(mtr2)=0.28;Hd_1=Ad.*exp(-j*Wm*M/2);Hd_2=conj(fliplr(Hd_1(2:M/2)));Hd=[Hd_1,Hd_2];hk=real(ifft(Hd));w=linspace(0,pi,1000);H=freqz(hk,1,w);%归一化频率下的幅频响应plot(w/pi,abs(H));grid;xlabel('Normalizedfrequency');ylabel('GainindB');title('频率取样法设计的FIR滤波器');已知理想低通滤波器为,矩形窗函数求理想低通滤波器的单位脉冲响应,并画出。当时,画出矩形窗函数的幅频响应。,画出加窗处理以后的低通滤波器的幅频响应。实现过程:clc;clearall;OmegaC=0.5*pi;M=15;k=0:M;hd=OmegaC*sinc(OmegaC*(k-0.5*M))/pi;subplot(311);stem(k,hd,'.');grid;title('理想低通滤波器的单位脉冲响应');wk=ones(1,M+1);w=linspace(-pi,pi,1000);Wm=freqz(wk,1,w);subplot(312);plot(w/pi,abs(Wm));grid;title('矩形窗函数的幅频响应N=16');hk=hd.*wk;w=linspace(-pi,pi,1000);H=freqz(hk,1,w);subplot(313);plot(w/pi,abs(H));grid;title('加窗处理以后的低通滤波器的幅频响应');四、思考题FIR滤波器是否需要考虑稳定性问题?为什么?答:不需要;FIR滤波器的单位脉冲响应是有限长的,系统总是稳定的窗函数法和频率抽样法的优缺点是什么?答:窗函数法是利用有限长的单位脉冲响应h[k]逼近无限长的理想滤波器的hd[k],从而使设计的FIR滤波器的频率响应逼近理想滤波器的频率响应频率取样法是使设计的M阶FIR滤波器的频率响应在M+1个取样点上与理想滤波器的频率响应相等,不足的是设计出的FI

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论