离散时间滤波器设计.docVIP

实验六 离散时间滤波器设计 实验原理 IIR数字滤波器的设计 脉冲响应不变法将模拟滤波器的s平面变换成数字滤波器的z平面，从而将模拟滤波器映射成数字滤波器。 为克服脉冲响应不变法产生频率响应的混叠失真，可以采用非线性频率压缩方法，使s平面与z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，这就是双线性变换法。 窗函数法设计FIR数字滤波 FIR滤波器的设计问题，就是要使所设计的FIR滤波器的频率响应逼近所要求的理想滤波器的频率响应。逼近可在时域进行，也可在频域进行。窗函数设计FIR数字滤波器是在时域进行的，用窗函数截取无限长的hd(n),，这样得到的频率响应逼近于理想的频率响应。 实验内容 IIR数字滤波器的设计 要求通带截止频率fp=3KHz,，通带最大衰减ap=1dB，阻带截止频率fs=4.5KHz,，阻带最小衰减as=15dB，采样频率fc=30KHz，用脉冲响应不变法设计一个切比雪夫数字低通滤波器，并图示滤波器的振幅特性，检验wp,ws对应的衰减。 （1）切比雪夫Ⅰ型滤波器的设计 程序： wp=6000*pi;ws=9000*pi; %通带截止频率和阻带截止频率 ap=1;as=15;Fs=30000; %通带最大衰减和阻带最小衰减 wp1=wp/Fs;ws1=ws/Fs; [N,WC]=cheb1ord(wp,ws,ap,as,s); [b,a]=cheby1(N,ap,WC,s); %切比雪夫Ⅰ型 [bz,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))) plot(W/pi,abs(H)); xlabel(相对频率);ylabel(振幅); grid %加网格 结果： （2）切比雪夫Ⅱ型滤波器的设计 程序： clear all clc wp=6000*pi;ws=9000*pi; ap=1;as=15; Fs=30000; wp2=wp/Fs;ws2=ws/Fs; [N,Wc]=cheb2ord(wp,ws,ap,as,s); [b,a]=cheby2(N,as,Wc,s); %切比雪夫Ⅱ型 [bz,az]=impinvar(b,a,Fs) w0=[wp2,ws2]; Hx=freqz(bz,az,w0); [H,W]=freqz(bz,az); dbHx=20*log10(abs(Hx)/max(abs(H))); plot(W,abs(H)) xlabel(相对频率);ylabel(幅频); grid 结果： 用双线性变换法设计一个切比雪夫Ⅰ数字高通滤波器。技术指标为：采样频率fc=2KHz,通带截止频率fp=700Hz, 通带最大衰减ap=1dB, 阻带边缘频率fs=500Hz, 阻带最小衰减as=32dB。 程序： clear all Clc wp=1400*pi;ws=1000*pi; %通带截止频率和阻带截止频率 ap=1;as=32;Fs=2000*pi; %通带最大衰减和阻带最小衰减 wp1=wp/Fs;ws1=ws/Fs; omp1=2*Fs*tan(wp1/2); omps=2*Fs*tan(ws1/2); [N,WC]=cheb1ord(omp1,omps,ap,as,s); [b,a]=cheby1(N,ap,WC,s); [bz,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))); plot(W,abs(H)) xlabel(相对频率);ylabel(幅频); grid 结果： 窗函数法设计FIR数字滤波 窗函数法设计低通数字滤波器 N=26,分别利用矩形窗，汉宁窗和布莱克曼窗设计该滤波器，且滤波器具有线性相位。绘出脉冲响应h(n)及滤波器的频率响应。 M文件： function y=chuangdt(N) wc=0.4; h=fir1(N-1,wc,boxcar(N)); [H,W]=freqz(h,1); subplot(3,3,1); plot(h); title(脉冲响应); subplot(3,3,2); plot(W/pi,abs(H));