fir滤波器设计与实现

《fir滤波器设计与实现》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、FIR滤波器设计与实现实验报告实验时间：2014年5月29日专业：微电子学实验人员：12微电子许一飞、祁响、刘秦华一．实验目的:（1）通过实验巩固FIR滤波器的认识和理解。（2）熟悉掌握FIR低通滤波器的窗函数设计方法。（3）理解FIR的具体应用。二．实验内容:在通信、信息处理以及信号检测等应用领域广泛使用滤波器进行去噪和信号的增强。FIR滤波器由于可实现线性相位特性以及固有的稳定特征而等到广泛应用，其典型的设计方法是窗函数设计法。设计流程如下： （1）设定指标：截止频率fc，过渡带宽度△f，阻带衰减A。 （2）求理想低通滤波器（LPF）的时域响应hd（n）。 （3

2、）选择窗函数w（n），确定窗长N。 （4）将hd（n）右移（N-1）/2点并加窗获取线性相位FIR滤波器的单位脉冲响应h（n）。 （5）求FIR的频域响应H（ejw），分析是否满足指标。如不满足，转（3）重新选择，否则继续。 （6）求FIR的系统函数H（z）。 （7）依据差分方程由软件实现FIR滤波器或依据系统函数由硬件实现。  实验要求采用哈明窗设计一个FIR低通滤波器并由软件实现。哈明窗函数如下： w（n）=0.54-0.46cos（2πn/（N-1）），0≤n≤N-1； （公式 1） 设采样频率为fs=10kHz。实验中，窗长度N和截止频率fc应该都能调节。具

3、体实验内容如下： （1）设计FIR低通滤波器（FIR_LPF）（书面进行）。 （2）依据差分方程编程实现FIR低通滤波器。 （3）输入信号x（n）=3.0sin（0.16πn）+cos（0.8πn）到fc=2000Hz，N=65的FIR_LPF，求输出信号y（n），理论计算并画出0≤f≤fs范围输入信号x（n）和输出信号y（n）的幅度谱，标出峰值频率，观察滤波器的实际输出结果，分析其正确性。 （4）输入信号x（n）=1.5sin（0.2nπ）-cos（0.4nπ）+1.2sin（0.9nπ）到fc=1100Hz，N=65的FIR_LPF，求输出信号y（n），理论计算

4、并画出0≤f≤fs范围输入信号x（n）和输出信号y（n）的幅度谱，标出峰值频率，观察滤波器的实际输出结果，分析其正确性。 （5）输入信号x（n）=1.5sin（0.2nπ）-cos（0.4nπ）+1.2sin（0.9nπ）到fc=2100Hz，N=65的FIR_LPF，求输出信号y（n），理论计算并画出0≤f≤fs范围输入信号x（n）和输出信号y（n）的幅度谱，标出峰值频率，观察滤波器的实际输出结果，分析其正确性。 （6）输入信号x（n）=1.5sin（0.2nπ）-cos（0.4nπ）+5.0sin（0.9nπ）到fc=1100Hz，N=65的FIR_LPF，求输

5、出信号y（n），理论计算并画出0≤f≤fs范围输入信号x（n）和输出信号y（n）的幅度谱，标出峰值频率，观察滤波器的实际输出结果，分析其正确性。 （7）输入信号x（n）=1.5sin（0.2nπ）-cos（0.4nπ）+1.2sin（0.9nπ）到fc=1990Hz，N=65的FIR_LPF，求输出信号y（n），理论计算并画出0≤f≤fs范围输入信号x（n）和输出信号y（n）的幅度谱，标出峰值频率，观察滤波器的实际输出结果，分析其正确性。三．实验程序functionf=FIR(A,a,B,b,C,c,fc,N,S)fs=10000;wc=fc*2*pi/fs;n=0

6、:2*N-1;x=A*sin(a*pi*n)+B*cos(b*pi*n)+C*sin(c*pi*n);n=0:(N-1);w=0.54-0.46*cos(2*pi*n/(N-1));m=n-(N-1)/2+eps;hd=(sin(wc*m))./(pi*m);h=hd.*w;form=0:N-1y(m+1)=0;fork=0:N-1y(m+1)=y(m+1)+h(k+1)*x(m-k+1+N);endendx=x(n+1+N);s=0:S-1;f=fs*s/S;X=abs(fft(x,S))/(N/2);H=abs(fft(h,S));Y=X.*H;maxim=0;

7、maxi=0;maxk=0;maxl=0;fork=0:S-1ifY(k+1)>maximmaxim=Y(k+1);maxk=k;endendforl=0:S-1ifX(l+1)>maximaxi=X(l+1);maxl=l;endendfa=maxl/S*fs;fb=maxk/S*fs;Ma=max(X);Mb=max(Y);figure(1);subplot(6,1,1);plot(x,'.-')title('x(n)');subplot(6,1,2);plot(n,h,'.-')title('h(n)');subplot(6,1,3);plot(y,'.-