实验五 FIR 滤波器的设计.docVIP

HUNAN UNIVERSITY 数字信号 课程实验报告 专业班级： 通信工程一班 完 成 日 期 ： 　 2014\06\01 实验五 FIR 滤波器的设计 1、实验目的 认真复习 FIR 数字滤波器的基本概念，线性相位FIR 滤波器的条件和特点、幅度函数 特点、零点位置的基本特点与性质；窗函数设计法的基本概念与方法，各种窗函数的性能和 设计步骤，线性相位FIR 低通、高通、带通和带阻滤波器的设计方法，频率采样设计法的基本概念和线性相位的实现方法。 掌握几种线性相位的特点，熟悉和掌握矩形窗、三角形窗、汉宁窗、海明窗、布莱克曼 窗、凯塞窗设计IIR 数字滤波器的方法，熟悉和掌握频率抽样设计法的线性相位的设计方法， 并对各种线性相位的频率抽样法的设计给出调整和改进。 熟悉利用 MATLAB 进行各类FIR 数字滤波器的设计方法。 2、实验内容 a. 设线性相位FIR 滤波器单位抽样响应分别为 h(n) ={ -4,1, -1, -2,5,6,5, -2, -1,1, -4} ↑ h(n)={ -4,1,- 1, -2,5,6,6,5, -2,-1,1, -4} ↑ h(n) ={ -4,1, -1, -2,5,0, -5,2,1, -1,4} ↑ h(n) { -4,1, -1, -2,5,6, -6, -5,2,1, -1,4} ↑ 分别求出滤波器的幅度频率响应 H(ω)，系统函数H(z)以及零极点分布，并绘制相应的 波形和分布图。 1、 代码： function [H,w,a,L]=H_Type1(h); %求幅度频率响应 M=length(h); L=floor((M-1)/2); a=[h(L+1) 2*h(L:-1:1)]; n=[0:1:L]; w=[0:1:500]*pi/500; H=cos(w*n)*a; 主函数 clc; syms z; %偶对称单位冲激响应h(n)，N为奇数 h1=[-4,1,-1,-2,5,6,5,-2,-1,1,-4]; N1=length(h1); n1=[0:1:N1-1]; [H1,w1,a1,L1]=H_Type1(h1); subplot(221); stem(h1); title(h1(n)); subplot(222); plot(H1); title(幅度函数H1(w)); subplot(223); zplane(h1); title(零极点); HZ1=h1./(z.^n1) %系统函数H(z) 结果： HZ1 = [ -4, 1/z, -1/z^2, -2/z^3, 5/z^4, 6/z^5, 5/z^6, -2/z^7, -1/z^8, 1/z^9, -4/z^10] 2、 代码： function [H,w,b,L]=H_Type2(h); %求幅度频率响应 M=length(h); L=floor(M/2); b=2*h(L:-1:1); n=[1:1:L]; w=[0:1:500]*pi/500; H=cos(w*(n-1/2))*b; 主函数 %偶对称单位冲激响应h(n)，N为偶数 figure; h2=[-4,1,-1,-2,5,6,6,5,-2,-1,1,-4]; N2=length(h2); n2=[0:1:N2-1]; [H2,w2,a2,L2]=H_Type2(h2); subplot(221); stem(h2); title(h2(n)); subplot(222); plot(H2); title(幅度函数H2(w)); subplot(223); zplane(h2); title(零极点); HZ2=h2./(z.^n2) %系统函数H(z) 结果： HZ2 = [ -4, 1/z, -1/z^2, -2/z^3, 5/z^4, 6/z^5, 6/z^6, 5/z^7, -2/z^8, -1/z^9, 1/z^10, -4/z^11] 3、 代码： function [H,w,c,L]=H_Type3(h); %求幅度频率响应 M=length(h); L=floor((M-1)/2); c=2*h(L:-1:1); n=[1:1:L]; w=[0:1:500]*pi/500; H=sin(w*n)*c; 主函数： %奇对称单位冲激响应h(n)，N为奇数