实验五fir数字滤波器设计与软件实现.docxVIP

实验五 FIR数字滤波器设计与软件实现 1实验目的 掌握用窗函数法设计 FIR数字滤波器的原理和方法。 掌握用等波纹最佳逼近法设计 FIR数字滤波器的原理和方法。 掌握FIR滤波器的快速卷积实现原理。 学会调用MATLAB函数设计与实现 FIR滤波器。 实验内容及步骤 认真复习第七章中用窗函数法和等波纹最佳逼近法设计 FIR数字滤波器的原理； 调用信号产生函数 xtg产生具有加性噪声的信号 xt,并自动显示xt及其频谱，如 图10.5.1所示； ⑻信号加噪声波形100-100.05 0.10.15 0.2 0.25 0.3 0.35 0.4t/s(b)信号加噪声的频谱10.5050150200250300350400450500f/Hz ⑻信号加噪声波形 10 0 -10 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 t/s (b)信号加噪声的频谱 1 0.5 0 50 150 200 250 300 350 400 450 500 f/Hz 1 L dm」 ...,Li： 1 山LU .1 I J 100 图10.5.1具有加性噪声的信号 x(t)及其频谱如图 请设计低通滤波器，从高频噪声中提取 xt中的单频调幅信号，要求信号幅频失真 小于0.1dB，将噪声频谱衰减 60dB。先观察xt的频谱，确定滤波器指标参数。 根据滤波器指标选择合适的窗函数，计算窗函数的长度 N，调用MATLAB函数 fir1设计一个FIR低通滤波器。并编写程序，调用 MATLAB快速卷积函数fftfilt实现对xt 的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 (4)重复(3),滤波器指标不变，但改用等波纹最佳逼近法， 调用MATLAB函数remezord 和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。 7章和第8章;提示：MATLAB函数fir1和 7章和第8章; 采样频率 Fs=1000Hz，采样周期 T=1/Fs ; 根据图1061(b)和实验要求，可选择滤波器指标参数：通带截止频率 fp=120Hz，阻 带截至频率fs=150Hz，换算成数字频率，通带截止频率 ? .p =2二fp 一二0.24二，通带最大衰为 0.1dB，阻带截至频率 飞=2二fs 一 =0.3二，阻带最小衰为 60dB。］ ④实验程序框图如图10.5.2所示，供读者参考。 信号产生函数xtg程序清单 function xt=xtg(N) %实验五信号x(t)产生，并显示信号的幅频特性曲线 %Xt=xtg(N) 产生一个长度为N,有加性高频噪声的单频调幅信号 xt,采样频率Fs=1000Hz %载波频率fc=Fs/10=100Hz, 调制正弦波频率 f0=fc/10=10Hz. N=2000;Fs=1000;T=1/Fs;Tp=N*T; t=0:T:(N-1)*T; fc=Fs/10;f0=fc/10; % mt=cos(2*pi*f0*t); % ct=cos(2*pi*fc*t); % xt=mt.*ct; % 载波频率fc=Fs/10 ，单频调制信号频率为 f0=Fc/10; 产生单频正弦波调制信号 mt，频率为f0 产生载波正弦波信号 ct，频率为fc 相乘产生单频调制信号 xt  nt=2*ra nd(1,N)-1; % 产生随机噪声 nt %=======设计高通滤波器 hn,用于滤除噪声 nt中的低频成分 ,生成高通噪声 fp=150; fs=200;Rp=0.1;As=70; % 滤波器指标 fb=[fp,fs];m=[0,1]; % 计算 remezord 函数所需参数 f,m,dev dev=[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(h n,1,10* nt); % 滤除随机噪声中低频成分，生成高通噪声 yt %============================================================== xt=xt+yt; % 噪声加信号 fst=fft(xt,N);k=O:N-1;f=k/Tp; subplot(3,1,1);plot(t,xt);grid;xlabel(t/s);ylabel(x(t)); axis([0,Tp/5,min(xt),max(xt)]);title((a) 信号加噪声波形) subplot(3,1,2);plot(f