数字信号处理—基于计算机的方法-FIR数字滤波器的设计和应用.docVIP

电子科技大学生命科学与技术学院 标 准 实 验 报 告 （实验）课程名称 数字信号处理 2013-2014-第2学期 电子科技大学教务处制表 电 子 科 技 大 学 实 验 报 告 学生姓名 学 号 指导教师 实验时间 2014年 月 日 一、实验室名称：清水河校区，实验大楼 A406 机房 二、实验名称：FIR数字滤波器的设计与应用 三、实验学时：2学时 四、实验原理： 五、实验目的：（详细填写） 利用加窗傅立叶级数（即对时间域的无限长冲激响应进行加窗截短）的方法，设计FIR数字滤波器： 1）掌握其基本思想和步骤； 2）学会使用有关的MATLAB函数。 六、实验内容：（详细填写） 加载有用信号(EEGdata2.txt)；分析信号的幅度谱，确定低通FIR数字滤波器的指标；分别利用3种窗函数（Rectangular, Hamming, Kaiser）设计此FIR滤波器；对信号加随机噪声，并用设计的滤波器对含噪声信号进行滤波。 要求： 1) 画出滤波器幅度谱及其增益图；分析设计的滤波器是否达到指标要求； 2) 画出滤波器的相位谱，是否是线性相位？ 3) 分析滤波器的阶数； 4) 利用设计的三种滤波器对加载的信号分别进行滤波，对比分析滤波结果(时间域及频谱)。 七、实验器材（设备、元器件）： 八、实验步骤： 九、实验数据及结果分析：（详细填写）（包括程序、图、结果等） 程序 %%载入EEGdata2.txt并分析其信号谱和幅度谱，并分析基本指标 data=load(EEGdata2.txt); n=length(data); x0=1:n; x0=x0/n*2; f0=abs(fft(data)); figure(1); subplot(2,1,1),plot(x0,data,r);title(原信号时间域信号谱); subplot(2,1,2),plot(x0,f0,b);title(原信号幅度谱); %% %噪声信号添加 noise=20*rand(n,1)-10; yy=data+noise; f1=abs(fft(yy)); figure(2); subplot(2,1,1),plot(x0,yy,r);title(叠加信号谱); subplot(2,1,2),plot(x0,f1,b);title(叠加信号幅度谱); %% %输入数字指标，确定FIR滤波器阶数 wp=0.23*pi;ws=0.28*pi;Rs=50;Rsm=10^(-Rs/20); transWidth = ws-wp; M=[0.92*pi 3.32*pi]; M=M/transWidth;%[Rectangular Hamming ] N=ceil(2*M+1); fpts=[0.23 0.27]; mag=[1 0]; dev=[Rsm Rsm]; [N(3),Wn0,beta,ftype]=kaiserord(fpts,mag,dev); disp(N); %% %设计给定阶数的窗函数 L=N+1; w1=rectwin(L(1)); w2=hamming(L(2)); w3=kaiser(L(3),beta); %%确定滤波器的传输函数,矩形窗滤波器幅度谱相位谱分析 Wn=(wp+ws)/2/pi; b1=fir1(N(1),Wn,low,w1); [h1,ww1]=freqz(b1,1,512); ww1=ww1/pi; ff1=abs(h1);db1=20*log10(ff1); pha1=angle(h1);pha1=unwrap(pha1); figure(3);subplot(1,3,1),plot(ww1,ff1,b);title(Rectangular滤波器幅度谱); subplot(1,3,2),plot(ww1,db1,b);title(Rectangular滤波器增益谱);%幅度谱和增益谱 subplot(1,3,3),plot(ww1,pha1,b);title(Rectangular滤波器相位谱);%相位谱 %% %Hamming滤波器幅度谱相位谱分析 b2=fir1(N(2),Wn,low,w2); [h2,ww2]=freqz(b2,1,512);ww2=ww2/pi; ff2=abs(h2);db2=20*log10(ff2); pha2=angle(h2);pha2=unwrap(pha2); figure(4);subplot(1,3,1),plot(ww2,ff2,b);title(Hamming滤波器幅度谱); subplot(1,3,2),plot(ww2,db2,b);title(Hamming滤波