实验八 MATLAB语言的基本使用方法.docVIP

实验八 MATLAB语言的基本使用方法 实验目的： 1、掌握用窗函数法设计FIR数字滤波器的内容。 2、熟悉线性相位FIR数字滤波器特性。 3、了解各种窗对滤波器的影响。 实验原理： 有限长单位冲激响应（FIR）数字滤波器的特点： 1、可以实现严格的线性相位，可具有任意的幅度特性。 2、FIR DF的h(n)是有限长的，因而滤波器一定是稳定的。 3、任何非因果的有限长序列，经过一定的延时，都能成为因果的有限长序列，因而非因果系统。 实验结果： 理想线性相位低通滤波器的单位冲击响应hd(n) function hd=ideallp(wc,N) tao=(N-1)/2; n=[0:1:(N-1)]; m=n-tao+eps; hd=sin(wc*m)./(pi*m); 检验函数 function [db,mag,pha,grd,w]=freq_z(b,a) [H,w]=freqz(b,a,1000,whole); H=(H(1:1:501)); w=(w(1:1:501)); mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(H); grd=grpdelay(pha); M文件： clc;clear all wp=0.2*pi; ws=0.3*pi; deltaw=ws-wp; N0=ceil(6.6*pi/deltaw); N=N0+mod(N0+1,2); n=[0:N-1]; wd=(hamming(N)); wc=(wp+ws)/2; hd=ideallp(wc,N); h=hd.*wd; [db,mag,pha,grd,w]=freqz_m(h,[1]); dw=2*pi/1000; Rp=-min(db(1:wp/dw+1:501)); subplot(1,2,1);plot(n,wd);grid; title(海明窗); xlabel(n); ylabel(w(n)); axis([0,N,0,1]); subplot(1,2,2);plot(w/pi,db); title(幅度响应（db）); xlabel(\omega/\pi); ylabel(20log|H(e^j^\omega/\pi)); axis([0,1,-80,5]);grid; N=45;Rp=0.0428;As=50.3822; Sy8_1 M文件： clc;clear all wp1=0.35*pi; ws1=0.2*pi; wp2=0.65*pi; ws2=0.8*pi; AS=60; deltaw=min((wp1-ws1),(ws2-wp2)); N=ceil(11*pi/deltaw); M=N-1; n=[0:N-1]; wc1=(ws1+wp1)/2; wc2=(ws2+wp2)/2; wn=[wc1/pi,wc2/pi]; h=fir1(M,wn,blackman(N)); [db,mag,pha,grd,w]=freqz_m(h,[1]); dw=2*pi/1000; Rp=-min(db(wp1/dw+1:1:wp2/dw+1)); As=-min(db(ws2/dw+1:1:501)); subplot(1,2,1); stem(n,h,.);grid; title(布莱克曼窗); xlabel(n); ylabel(w(n)); axis([0,N,0,0.4]); subplot(1,2,2);plot(w/pi,db); title(幅度响应（db）); xlabel(\omega/\pi); ylabel(20log|H(e^j^\omega）|（db)); axis([0,1,-120,20]);grid; N=110;Rp=0.0036;As=72.3801; Sy8_2