DSP数字信号处理实验报告-FIR滤波器设计与实现.doc

数字信号处理基于哈明窗的fir数字滤波器的设计dsp实验报告 fir滤波器设计与实现一、实验目的（1） 通过实验巩固fir滤波器的认识和理解。（2） 熟练掌握fir低通滤波器的窗函数设计方法。（3） 理解fir的具体应用。二、实验内容在通信、信息处理以及信号检测等应用领域广泛使用滤波器进行去噪和信号的增强。fir滤波器由于可实现线性相位特性以及固有的稳定特征而等到广泛应用，其典型的设计方法是窗函数设计法。设计流程如下：（1）设定指标：截止频率fc，过渡带宽度f，阻带衰减a。（2）求理想低通滤波器（lpf）的时域响应hd（n）。（3）选择窗函数w（n），确定窗长n。（4）将hd（n）右移（n-1）/2点并加窗获取线性相位fir滤波器的单位脉冲响应h（n）。（5）求fir的频域响应h（e），分析是否满足指标。如不满足，转（3）重新选择，否则继续。（6）求fir的系统函数h（z）。（7）依据差分方程由软件实现fir滤波器或依据系统函数由硬件实现。实验要求采用哈明窗设计一个fir低通滤波器并由软件实现。哈明窗函数如下：w（n）=0.54-0.46cos（），0nn-1；设采样频率为fs=10khz。实验中，窗长度n和截止频率fc应该都能调节。具体实验内容如下：（1）设计fir低通滤波器（fir_lpf）（书面进行）。（2）依据差分方程编程实现fir低通滤波器。（3）输入信号x（n）=3.0sin（0.16）+cos（0.8）到fc=2000hz，n=65的fir_lpf，求输出信号y（n），理论计算并画出0ffs范围输入信号x（n）和输出信号y（n）的幅度谱，标出峰值频率，观察滤波器的实际输出结果，分析其正确性。（4）输入信号x（n）=1.5sin（0.2）-cos（0.4）+1.2sin（0.9）到fc=1100hz，n=65的fir_lpf，求输出信号y（n），理论计算并画出0ffs范围输入信号x（n）和输出信号y（n）的幅度谱，标出峰值频率，观察滤波器的实际输出结果，分析其正确性。（5）输入信号x（n）=1.5sin（0.2）-cos（0.4）+1.2sin（0.9）到fc=2100hz，n=65的fir_lpf，求输出信号y（n），理论计算并画出0ffs范围输入信号x（n）和输出信号y（n）的幅度谱，标出峰值频率，观察滤波器的实际输出结果，分析其正确性。（6）输入信号x（n）=1.5sin（0.2）-cos（0.4）+5.0sin（0.9）到fc=1100hz，n=65的fir_lpf，求输出信号y（n），理论计算并画出0ffs范围输入信号x（n）和输出信号y（n）的幅度谱，标出峰值频率，观察滤波器的实际输出结果，分析其正确性。（7）输入信号x（n）=1.5sin（0.2）-cos（0.4）+1.2sin（0.9）到fc=1990hz，n=65的fir_lpf，求输出信号y（n），理论计算并画出0ffs范围输入信号x（n）和输出信号y（n）的幅度谱，标出峰值频率，观察滤波器的实际输出结果，分析其正确性。三、实验分析如果采用哈明窗函数设计因果线性相位低通滤波器fir，所设计的fir低通滤波器单位脉冲响应h（n）如式。h（n）=hd（n）w（n）=，0nn-1如果输入函数为x（n），则y（n）=x（n）*h(n)。四、实验结果（1）matlab 程序如下：function fir = fir(fc,n,a,a,b,b,c,c)%窗函数设计n = 0 : n-1; %窗函数、lpf的横坐标fs = 10000; %采样频率10khzwc = 2 * pi * fc/fs;hd = n;for k = 1 : (n-1)/2 hd(k) = sin(wc * (k-1-(n-1)/2) ./ ( pi *(k-1-(n-1)/2); %lpf纵坐标endfor k =(n-1)/2 + 2 : n hd(k) = sin(wc * (k-1-(n-1)/2) ./ ( pi *(k-1-(n-1)/2); %lpf纵坐标endhd(n-1)/2+1) = wc/pi;wn = 0.54 - 0.46*cos(2*pi*n/(n-1); %哈明窗窗函数hn = hd .* wn; %fir滤波器时域表达式%窗函数频谱w = 0 : 0.001*pi :0.999* pi;w = 0:999;for k = 1 : 1000; w(k) = 0; n = 1; while(n n+1) w(k) = hd(n) * exp(-j*k*pi/1000*n)+ w(k); n = n + 1; endend%输入频谱n = 0 : 999;x = a * sin ( a * pi * n ) + b * cos ( b * pi * n) + c* sin ( c * pi * n );x = x;s = 0;for k = 1 : 1000 x(k) = 0; n = 1; while( n 1001) x(k) = x(k) + x(n) * exp(-j * 2* pi *(k-1)*(n-1)/1000); n = n + 1; end if k 500 if s abs(x(k)& (s = abs(x(k) s = abs(x(k); m = k - 1; end endendy = x .* w;t = 0 : n-1; %过渡变量k = 0 : 0.002 * pi : 1.998 *pi;subplot(4,1,1); %图1plot(t,hn); %窗函数时域图subplot(4,1,2); %图2plot(k/2,20*log10(abs(w); %窗函数频谱图subplot(4,1,3); %图3plot(k,abs(x),.-); %输入函数频谱图subplot(4,1,4); plot(k,abs(y),.-);f = m /1000 * fssend (3) x（n）=3.0sin（0.16）+cos（0.8）当f=800hz时，峰值为1500(4) x（n）=1.5sin(0.2)-cos(0.4)+1.2sin(0.9)当f=1000hz时，峰值为750(5) x（n）=1.5sin(0.2)-cos(0.4)+1.2sin(0.9)当f=1000hz时，峰值为750(6) x（n）=1.5sin(0.2)-cos(0.4)+5.0sin(0.9)当f=4500hz时，峰值为2500(7) x（n）=1.5sin(0.2)-cos(0.4)+1.2sin(0.9)当f=1000hz时，峰值为750五、思考题1、当哈明窗长度n比65小（32）或大（129）的话，实验结果如何变化？（1）当n=32时a：x（n）=3.0sin（0.16）+cos（0.8） b：x（n）=1.5sin(0.2)-cos(0.4)+1.2sin(0.9)c：x（n）=1.5sin(0.2)-cos(0.4)+1.2sin(0.9)d：x（n）=1.5sin(0.2)-cos(0.4)+5.0sin(0.9)e x（n）=1.5sin(0.2)-cos(0.4)+1.2sin(0.9)（2）当n=129时a：x（n）=3.0sin（0.16）+cos（0.8） b：x（n）=1.5sin(0.2)-cos(0.4)+1.2sin(0.9)c：x（n）=1.5sin(0.2)-cos(0.4)+1.2sin(0.9)d：x（n）=1.5sin(0.2)-cos(0.4)+5.0sin(0.9)e x（n）=1.5sin(0.2)-cos(0.4)+1.2sin(0.9)2、当采用矩形窗的话，实验（3）（4）的结果是怎样的？矩形窗