数字信号管理方案计划实验报告实验五

.\物理与电子信息工程学院实验报告实验课程名称：数字信号处理实验名称：FIR数字滤波器设计与软件实现班 级：1012341姓 名：严娅学 号：101234153成 绩：_______实验时间：2012年12月20日h(n).\一、实验目的（1）掌握用窗函数法设计FIR数字滤波器的原理和方法。谢谢阅读（2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。谢谢阅读（3）掌握FIR滤波器的快速卷积实现原理。（4）学会调用MATLAB函数设计与实现FIR滤波器。精品文档放心下载二、实验原理1、用窗函数法设计FIR数字滤波器的原理和方法。谢谢阅读如果所希望的滤波器的理想频率响应函数为 H(ej)，则谢谢阅读d其对应的单位脉冲响应为h(n)=1H(ej)ejd（2-1）dd2窗函数设计法的基本原理是用有限长单位脉冲响应序列h(n)逼近h(n)。由于h(n)往往是无限长序列，且是非因果的，谢谢阅读dd截断，并进行加权处理，得到：所以用窗函数(n)将h(n)d(2-2)h(n)＝h(n)(n)dh(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数H(ej)为：谢谢阅读dH(ej)＝h(n)ej(2-3)N1d0中，N为所选窗函数(n)的长度。第七章可知，用窗函数法设计的滤波器性能取决于窗函数(n)的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见第七章。感谢阅读这样选定窗函数类型和长度N后，求出单位脉冲响应＝h(n)·(n)，并按式（2-3）求出H(ej)。H(ej)是否满足要感谢阅读求，要进行验算。一般在h(n)尾部加零使长度满足于2的整数次幂，以便用FFT计算H(ej)。如果要观察细节，补零点数增多即可。如果H(ej)不满足要求，则要重新选择窗函数精品文档放心下载.\类型和长度N，再次验算，直至满足要求。如果要求线性相位特性，则h(n)还必须满足h(n)h(N1n)(2-4)根据上式中的正负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据设计的滤波特性正确选择其中一类。例如，要设计线性低通特征，可选择h(n)h(N1n)一类，而不能选h(n)h(N1n)一类。感谢阅读2、等波纹最佳逼近法等波纹最佳逼近法是一种优化设计法，它克服了窗函数设计法和频率采样法的缺点，使最大误差（即波纹的峰值）最小化，并在整个逼近频段上均匀分布。用等波纹最佳逼近法设计的FIR数字滤波器的幅频响应在通带和阻带都是等波纹的，而且可以分别控制通带和阻带波纹幅度。这就是等波纹的含义。最佳逼近是指在滤波器长度给定的条件下，使加权误差波纹幅度最小化。与窗函数设计法和频率采样法比较，由于这种设计法使滤波器的最大逼近误差均匀分布，所以设计的滤波器性能价格比最高。阶数相同时，这种设计法使滤波器的最大逼近误差最小，即通带最大衰减最小，阻带最小衰减最大；指标相同时，这种设计法使滤波器阶数最低。实现FIR数字滤波器的等波纹最佳逼近法的MATLAB信号处理工具函数为remez和remezord。Remez函数采用数值分析中的remez多重交换迭代算法求解等波纹最佳逼近问题，求的满足等波纹最佳逼近准则的FIR数字滤波器的单位脉冲响应h（n）。感谢阅读三、实验内容及步骤（1）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图10.5.1所示；感谢阅读.\图10.5.1具有加性噪声的信号x(t)及其频谱如图精品文档放心下载（2）设计滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。根据窗函数设计滤波器原理设计滤波器，根据滤波器指标选择合适的窗函数，计算窗函数的长度N，调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序，调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。谢谢阅读（3）改用等波纹最佳逼近法设计滤波器，滤波器指标不变，调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数和滤出来的波形效果。感谢阅读（4）修改程序：将步骤（2）中用窗函数设计的低通滤波器参数改为fp=100;fs=110，单独滤出上边带频谱线，并绘制波形图，谢谢阅读.\（5）修改程序：将步骤（3）中用等波纹设计的低通滤波器改为带通精品文档放心下载滤波器，单独滤出下边带频谱线，根据带通滤波器的参数选取通、阻感谢阅读带截止频率,带通滤波器的设计指标为：fsl=95；fpl=105;fpu=115;fsu=125;精品文档放心下载fb=[fsl,fpl,fpu,fsu];m=[0,1,0]；精品文档放心下载dev=[10^(-As/20),(10^(Rp/20)-1)/(10^(Rp/20)+1),10^(-As/20)]精品文档放心下载并绘制波形图。四、实验程序清单1．信号产生函数xtg程序清单functionxt=xtg(N)%实验五信号x(t)产生,并显示信号的幅频特性曲线感谢阅读%xt=xtg(N)产生一个长度为N,有加性高频噪声的单频调幅信号xt,采样频率Fs=1000Hz精品文档放心下载%载波频率fc=Fs/10=100Hz,调制正弦波频率f0=fc/10=10Hz.精品文档放心下载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);ct=cos(2*pi*fc*t);xt=mt.*ct;nt=2*rand(1,N)-1;精品文档放心下载

%产生单频正弦波调制信号mt，频率为f0%产生载波正弦波信号ct，频率为fc%相乘产生单频调制信号xt%产生随机噪声nt精品文档放心下载.\%设计高通滤波器hn,用于滤除噪声nt中的低频成分,生成高通噪声fp=150;fs=200;Rp=0.1;As=70;%滤波器指标fb=[fp,fs];m=[0,1];%计算remezord函数所需参数f,m,devdev=[10^(-As/20),(10^(Rp/20)-1)/(10^(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;subplot(2,1,1);plot(t,xt);grid;xlabel('t/s');ylabel('x(精品文档放心下载t)');axis([0,Tp/5,min(xt),max(xt)]);title('(a)信号加噪声波形')subplot(2,1,2);plot(f,abs(fst)/max(abs(fst)));grid;titl感谢阅读e('(b)信号加噪声的频谱')axis([0,Fs/2,0,1.2]);xlabel('f/Hz');ylabel('幅度')谢谢阅读2、实验程序清单%《数字信号处理(第三版）学习指导》第10章实验5程序exp5.m精品文档放心下载%FIR数字滤波器设计及软件实现clearall;closeall;.\%调用xtg产生信号xt,xt长度N=1000,并显示xt及其频谱谢谢阅读N=1000;xt=xtg(N);fp=120;fs=150;Rp=0.2;As=60;Fs=1000;%输入给定指标%(1)用窗函数法设计滤波器============================精品文档放心下载wc=(fp+fs)/Fs; %理想低通滤波器截止频率(关于pi归一化）精品文档放心下载B=2*pi*(fs-fp)/Fs;%过渡带宽度指标Nb=ceil(11*pi/B);%blackman窗的长度Nhn=fir1(Nb-1,wc,blackman(Nb));Hw=abs(fft(hn,1024));%求设计的滤波器频率特性ywt=fftfilt(hn,xt,N);%调用函数fftfilt对xt滤波%以下为用窗函数法设计法的绘图部分====subplot(2,1,1);%绘制低通滤波器幅频特性曲线精品文档放心下载myplot(hn,1);title('(a)低通滤波器幅频特性');subplot(2,1,2); %绘制滤除噪声后的信号波形精品文档放心下载yt='y_w(t)';tplot(ywt,T,yt);title('滤除噪声后的信号波形')%(2)用等波纹最佳逼近法设计滤波器=====================fb=[fp,fs];m=[1,0];%确定remezord函数所需参数f,m,devdev=[(10^(Rp/20)-1)/(10^(Rp/20)+1),10^(-As/20)];[Ne,fo,mo,W]=remezord(fb,m,dev,Fs);%确定remez函数所需精品文档放心下载.\参数hn=remez(Ne,fo,mo,W); %调用remez函数进行设计感谢阅读Hw=abs(fft(hn,1024)); %求设计的滤波器频率特性精品文档放心下载yet=fftfilt(hn,xt,N); %调用函数fftfilt对xt滤波精品文档放心下载%以下为用等波纹设计法的绘图部分=================精品文档放心下载subplot(2,1,1); %绘制低通滤波器幅频特性曲线感谢阅读myplot(hn,1);title('(a)低通滤波器幅频特性');subplot(2,1,2); %绘制滤除噪声后的信号波形感谢阅读yt='y_w(t)';tplot(ywt,T,yt);title('滤除噪声后的信号波形')五、实验程序运行结果用窗函数法设计滤波器，滤波器长度Nb=184。滤波器损耗函数和谢谢阅读滤波器输出yw(nT)分别如图10.5.2(a)和（b）所示。感谢阅读用等波纹最佳逼近法设计滤波器，滤波器长度Ne=83。滤波器损精品文档放心下载耗函数和滤波器输出ye(nT)分别如图10.5.2(c)和（d）所示谢谢阅读.\图10.5.2.\图10.5.2用窗函数法设计低通滤波器滤出上边带频谱滤波器损耗函数和滤波感谢阅读器输出yw(nT)分别如图10.5.3(a)和（b）所示。感谢阅读用等波纹最佳逼近法设计带通滤波器滤出下边带频谱滤波器损耗函精品文档放心下载数和滤波器输出yw(nT)分别如图10.5.3(c)和（d）所示。谢谢阅读.\图10.5.3.\六、实验心得及体会1、通过本次实验我学会了调用MATLAB函数设计和实现FIR数字滤波器，理解了窗函数法设计FIR数字滤波器的原理，和等波纹最佳逼近法设计FIR数字滤波器的原理和方法。谢谢阅读2、FIR数字滤波器和IIR数字滤波器的设计方法有很大差别，IIR数字滤波器的设计方法是利用模拟滤波器成熟的理