哈工程数字信号试验六.doc

数字信号处理实验试验六：离散时间滤波器设计班级：学号：姓名：指导老师：日期：一、 IIR数字滤波器设计1. 实验原理2. 程序试验1clearwp=6*pi*103;ws=9*pi*103;ap=1;as=15;%Fs=30*103;%wp1=wp/Fs;ws1=ws/Fs;%N,Wc=cheb1ord(wp,ws,ap,as,s);%N3dBWcb,a=cheby1(N,ap,Wc,low,s);%1bz,az=impinvar(b,a,Fs);%w0=wp1,ws1;hx=freqz(bz,az,w0);H,W=freqz(bz,az);%dbHx=-20*log10(abs(hx)/max(abs(H);%wp1ws1plot(W,abs(H);xlabel();ylabel();%grid试验2clearwp=1.4*pi*103;ws=pi*103;ap=1;as=32;%Fs=2*103;%wp1=wp/Fs;ws1=ws/Fs;%omp1=2*Fs*tan(wp1/2);omp2=2*Fs*tan(ws1/2);%N,Wc=cheb1ord(omp1,omp2,ap,as,s);%N3dBWcb,a=cheby1(N,ap,Wc,high,s);%1bz,az=bilinear(b,a,Fs);%w0=wp1,ws1;hx=freqz(bz,az,w0);H,W=freqz(bz,az);%dbHx=-20*log10(abs(hx)/max(abs(H);%wp1ws1plot(W,abs(H);xlabel();ylabel();%grid3. 实验结果试验1试验2二、 窗函数法设计FIR数字滤波器1. 实验原理滤波器的理想频率响应函数为Hd(ej)，则其对应的单位脉冲响应为窗函数设计法的基本原理是用有限长单位脉冲响应序列h(n)逼hd(n)。由于hd(n)往往是无限长序列，且是非因果的，所以用窗函数。w(n)将hd(n)截断，并进行加权处理：h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数H(ej)为。用窗函数法设计的滤波器性能取决于窗函数w(n)的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。2. 程序试验1clearN=26;wc=0.4;h1=fir1(25,wc,boxcar(N)H,W=freqz(h1,1);subplot(221);plot(W/pi,abs(H);title(/dB);xlabel();ylabel(H(w);h2=fir1(25,wc,hanning(N)H,W=freqz(h2,1);subplot(222);plot(W/pi,abs(H);title(/dB);xlabel();ylabel(H(w);h3=fir1(25,wc,blackman(N)H,W=freqz(h3,1);subplot(223);plot(W/pi,abs(H);title(/dB);xlabel(); ylabel(H(w);clearN1=26;wc=0.4;h1=fir1(25,wc,boxcar(N1)H,W=freqz(h1,1);subplot(221);plot(W/pi,abs(H);title(/dB);xlabel();ylabel(H(w);N2=51;h2=fir1(50,wc,boxcar(N2)H,W=freqz(h2,1);subplot(222);plot(W/pi,abs(H);title(/dB);xlabel();ylabel(H(w);N3=101;h3=fir1(100,wc,boxcar(N3)H,W=freqz(h3,1);subplot(223);plot(W/pi,abs(H);title(/dB);xlabel();ylabel(H(w);试验2clearN=51;wc=0.6;h=fir1(50,wc,high,kaiser(N,8)H,W=freqz(h,1)plot(W/pi,abs(H);title(/dB);xlabel();ylable(H(w)试验3clearN=66;wc=0.25;h=fir1(65