《实验四：FIR数字滤波器设计与软件实现.docx》由会员分享，可在线阅读，更多相关《实验四：FIR数字滤波器设计与软件实现.docx（4页珍藏版）》请在三一文库上搜索。
1、实验四：FIR数字滤波器设计与软件实现一、实验目的:(1)掌握用窗函数法设计 FIR数字滤波器的原理和方法。(2)掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。(3)掌握FIR滤波器的快速卷积实现原理。(4)学会调用MATLAB函数设计与实现 FIR滤波器。二、实验内容及步骤:(1)认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理；(2)调用信号产生函数 xtg产生具有加性噪声的信号 xt,并自动显示xt及其频谱，如 图1所示；图1具有加性噪声的信号 x(t)及其频谱如图(3)请设计低通滤波器，从高频噪声中提取 xt中的单频调幅信号，要求信号幅频失真 小于0.1
2、dB,将噪声频谱衰减 60dB。先观察xt的频谱，确定滤波器指标参数。(4)根据滤波器指标选择合适的窗函数，计算窗函数的长度N,调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序，调用 MATLAB快速卷积函数fftfilt实现对xt 的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。(4)重复(3),滤波器指标不变，但改用等波纹最佳逼近法， 调用MATLAB函数remezord 和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。提示:1J MATLAB函数力门和fftfilt的功能及其调用格式请查阅本课本;采样频率 Fs=1000
3、Hz,采样周期 T=1/Fs；fp=120Hz ,阻带截3)根据图10.6.1(b)和实验要求，可选择滤波器指标参数：通带截止频率至频率fs=150Hz,换算成数字频率，通带截止频率 ,=2nfpT=0.24n,通带最大衰为0.1dB,阻带截至频率 供=2nfsT=0.3n ,阻带最小衰为 60dB。实验程序框图如图2所示。图2实验程序框图三、实验程序：1、信号产生函数 xtg程序清单：function xt=xtg%实验五信号x(t)产生，并显示信号的幅频特性曲线%xt=xtg(N)产生一个长度为 N,有加性高频噪声的单频调幅信号xt,采样频率Fs=1000Hz%载波频率fc=Fs/10=1
4、00Hz,调制正弦波频率 f0=fc/10=10Hz.N=1000;Fs=1000;T=1/Fs;Tp=N*T;t=0:T:(N-1)*T;fc=Fs/10;f0=fc/10; % 载波频率fc=Fs/10,单频调制信号频率为f0=Fc/10;mt=cos(2*pi*f0*t);%产生单频正弦波调制信号mt,频率为f0ct=cos(2*pi*fc*t);%产生载波正弦波信号ct,频率为fcxt=mt.*ct;%相乘产生单频调制信号xtnt=2*rand(1,N)-1;% 产生随机噪声 nt%=设计高通滤波器hn,用于滤除噪声nt中的低频成分，生成高通噪声= fp=150; fs=200;Rp=
5、0.1;As=70;% 滤波器指标fb=fp,fs;m=0,1;% 计算 remezord 函数所需参数 f,m,devdev=10A(-As/20),(10A(Rp/20)-1)/(10A(Rp/20)+1);n,fo,mo,W=remezord(fb,m,dev,Fs); % 确定 remez 函数所需参数hn=remez(n,fo,mo,W); %调用remez函数进行设计，用于滤除噪声 nt中的低频成分yt=filter(hn,1,10*nt);%滤除随机噪声中低频成分，生成高通噪声yt%=xt=xt+yt;%噪声加信号fst=fft(xt,N);k=0:N-1;f=k/Tp;subp
6、lot(2,1,1);plot(t,xt);grid;xlabel(t/s);ylabel(x(t);axis(0,Tp/5,min(xt),max(xt);title(信号加噪声波形)subplot(2,1,2);plot(f,abs(fst)/max(abs(fst);grid;title(b) 信号加噪声的频谱)axis(0,Fs/2,0,1.2);xlabel(f/Hz);ylabel(幅度) 2、主程序清单：clear all;close all;xt=xtg;%M 用 xtg 产生 xtfp=120; fs=150;Rp=0.1;As=60;Fs=1000;N=1000;%B 入给
7、定指标wc=(fp+fs)/Fs;%三一化B=2*pi*(fs-fp)/Fs;验渡带宽度N0=ceil(11*pi/B);%blackman 窗的长度hn=fir1(N0-1,wc,blackman(N0);Hw=abs(fft(hn,1024);湖波器频率特性yt=fftfilt(hn,xt,N);%寸 xt 滤波subplot(2,1,1);就图plot(20*log10(Hw/max(Hw);grid; axis(0,Fs/2,-120,20); xlabel(f/Hz);ylabel(幅度)； title(a)低通滤波器的幅频特性)； t=0:N-1/Fs;Tp=N/Fs;subplo
8、t(2,1,2);plot(t,yt);axis(0,Tp/2,-1,1);xlabel(t/s);ylabel(y(t);title(b)滤除噪声后的输出信号 y(t);四、思考题：(1)如果给定通带截止频率和阻带截止频率以及阻带最小衰减，如何用窗函数法设计线性相位低通滤波器？请写出设计步骤。答：用窗函数法设计线性相位低通滤波器的设计步骤教材中有详细的介绍。(2)如果要求用窗函数法设计带通滤波器,且给定通带上、下截止频率为 Wpl和Wpu阻带上、下截止频率为 Wsl和Wsu试求理想带通滤波器的截止频率 Wcl和Wcu答：希望逼近的理想带通滤波器的截止频率Wcl和Wcu分别为：-(；:;)/
9、2- (；:; )/ 2。cl ( sl pl ) j cu ( SU pu ) J(3)解释为什么对同样的技术指标，用等波纹最佳逼近法设计的滤波器阶数低？用窗函数法设计的滤波器，如果在阻带截止频率附近刚好满足，则离开阻带截止频率越远，阻带衰减富裕量越大，即存在资源浪费；几种常用的典型窗函数的通带最大衰减和阻带最小衰减固定，且差别较大，又不能分别控制。所以设计的滤波器的通带最大衰减和阻带最小衰减通常都存在较大富裕。如本实验所选的blackman窗函数，其阻带最小衰减为 74dB,而指标仅为60dB。 用等波纹最佳逼近法设计的滤波器，其通带和阻带均为等波纹特性，且通带最大衰减和阻带最小衰减可以分别控制，所以其指标均匀分布，没有资源浪费，所以期阶数低得多。