数字滤波器实验五.docVIP

实验五FIR 滤波器设计 实验目的 1.认真复习 FIR 数字滤波器的基本概念，线性相位FIR 滤波器的条件和特点、幅度函数特点、零点位置的基本特点与性质；窗函数设计法的基本概念与方法，各种窗函数的性能和设计步骤，线性相位FIR 低通、高通、带通和带阻滤波器的设计方法，频率采样设计法的基本概念和线性相位的实现方法。 2.掌握几种线性相位的特点，熟悉和掌握矩形窗、三角形窗、汉宁窗、海明窗、布莱克曼窗、凯塞窗设计IIR 数字滤波器的方法，熟悉和掌握频率抽样设计法的线性相位的设计方法，并对各种线性相位的频率抽样法的设计给出调整和改进。 3.熟悉利用 MATLAB 进行各类FIR 数字滤波器的设计方法。 实验内容 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)以及零极点分布，并绘制相应的波形和分布图。 实验代码： clc; syms z; %偶对称单位冲激响应h(n)，N为奇数 figure(name,偶对称单位冲激响应h(n)，N为奇数时的幅度函数和零极点分布----xieying,numbertitle,off) 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,r); title(h(n));grid; subplot(222);plot(H1,g); title(幅度函数H(w));grid;text(0,3,\rightarrow a(0)=3) subplot(223);zplane(h1); title(零极点);grid; HZ1=h1./(z.^n1)%求系统函数H(Z) figure(name,偶对称单位冲激响应h(n)，N为偶数时的幅度函数和零极点分布----xieying,numbertitle,off) %偶对称单位冲激响应h(n)，N为偶数 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,r); title(h(n));grid; subplot(222);plot(H2,g); title(幅度函数H(w));grid;%text(0,10,\rightarrow H2(0)=10) subplot(223);zplane(h2); title(零极点);grid; HZ2=h2./(z.^n2) figure(name,奇对称单位冲激响应h(n)，N为奇数时的幅度函数和零极点分布----xieying,numbertitle,off) %奇对称单位冲激响应h(n)，N为奇数 h3=[-4,1,-1,-2,5,0,-5,2,1,-1,4]; N3=length(h3); n3=[0:1:N3-1]; [H3,w3,a3,L3]=H_Type3(h3); subplot(221); stem(h3,r); title(h(n));grid; subplot(222);plot(H3,g); title(幅度函数H(w));grid; subplot(223);zplane(h3); title(零极点);grid; HZ3=h3./(z.^n3) figure(name,奇对称单位冲激响应h(n)，N为偶数时的幅度函数和零极点分布----xieying,numbertitle,off) %奇对称单位冲激响应h(n)，N为偶数 h4=[-4,1,-1,-2,5,6,-6,-5,2,1,-1,4]; N4=length(h4); n4=[0:1:N4-1]; [H4,w4,a4,L4]=H_Type4(h4); subplot(221); stem(h4,r); title(h(n));grid; subplot(222);plot(H4,g); title(幅度函数H(w));grid; subplot(223);zplane(h4); title(零极点); grid; HZ4=h4./(z.^n4)