matlab滤波器设计.概要.pptx

数字滤波器设计专题研讨;1．IIR 数字滤波器设计 设计一个IIR数字低通滤波器，其能取代下列指标的模拟低通滤波器(系统的抽样频率为44.1kHz) fp=2kHz， fs=10kHz , Ap=0.5dB, As=50dB 分别用双线性变换和冲激响应不变法设计一个BW型数字低通滤波器，并进行比较。;设计步骤;（1）冲激响应不变法设计一个BW型数字低通滤波器 仿真程序： Fs=44100; fs=10000;fp=2000; ws=fs*2*pi;wp=fp*2*pi; %确定模拟滤波器指标 Ws=ws/Fs;Wp=wp/Fs; Ap=0.5;As=50; N=buttord(wp,ws,Ap,As,‘s’);%确定模拟滤波器阶数 fprintf(N=%.0f\n,N); wc=wp/10^(0.1*Ap-1)^(1/N/2);%确定3dB截频 [numa,dena]=butter(N,wc,‘s’);%求模拟滤波器系统函数 [numd,dend]=impinvar(numa,dena,Fs); w=linspace(0,pi,512); h=freqz(numd,dend,w);%冲激响应不变法模拟到数字滤波器的转换 norm=max(abs(h)); numd=numd/norm; plot(w/pi,20*log10(abs(h)/norm)); xlabel(‘Normalized frequency ’); ylabel(‘Gain in dB’);%计算频谱 w=[Wp Ws]; h=freqz(numd,dend,w); fprintf(Ap=%.4f\n,-20*log10(abs(h(1)))); fprintf(‘As=%.4f\n’,-20*log10(abs(h(2))));%计算通带阻带衰减 ;仿真结果;（2）双线性变换法设计一个BW型数字低通滤波器 仿真程序： Fs=44100;fp=2000;fs=10000; Ap=0.5;As=50; wp=fp*2*pi; ws=fs*2*pi;%确定模拟滤波器指标 Wp=wp/Fs;Ws=ws/Fs; Omegap=2*Fs*tan(Wp/2); Omegas=2*Fs*tan(Ws/2); [N,wc]=buttord(Omegap,Omegas,Ap,As,‘s’);%确定模拟滤波器阶数和3dB截频 [bt,at]=butter(N,wc,‘s’);%求模拟滤波器系统函数 [bz,az]=bilinear(bt,at,Fs); w=linspace(0,pi,512); h=freqz(bz,az,w);%双线性变换法模拟到数字转换 norm=max(abs(h)); bz=bz/norm; plot((w/pi),20*log10(abs(h)/norm)); xlabel(‘Normalized frequency’);ylabel(‘Gain in dB’);%计算频谱 w=[Wp Ws]; h=freqz(bz,az,w); fprintf(Ap=%.4f\n,-20*log10(abs(h(1)))); fprintf(‘As=%.4f\n’,-20*log10(abs(h(2))));%计算通带阻带衰减 ;仿真结果：;结果分析：;2． 窗函数法设计FIR 数字滤波器 分别用Blackman窗和Kaiser窗法设计一个满足下列指标的线性相位的FIR低通滤波器 ?Wp=0.4p rad, Ap=0.5 dB, Ws=0.6p rad, As=55dB 并讨论两种窗函数的结果;设计步骤：;（1）.用Blackman窗设计线性相位的FIR低通滤波器 仿真程序：;仿真结果：;（2）.用Kaiser窗设计线性相位FIR低通滤波器;仿真结果：;Blackman窗;1.比较两种窗的设计结果可得 (1) 用Blackman窗设计的FIR低通滤波器N=53，通带和阻带衰减分别为Ap≈0dB，As≈74dB。 (2) 用Kaiser窗函数法设计的线性相位FIR数字滤波器长度N=27，Kaiser窗的参数β=3.3953.滤波器通带和阻带衰减分别为Ap≈0dB，As≈46dB 。 2.通过比较用Blackman窗和Kaiser 窗设计的两种结果可知，Kaiser窗设计出的滤波器更接近设计指标。相比于常用窗函数，用Kaiser窗设计出的滤波器阶数较低，但滤波器的阻带波纹衰减较慢。;问题探究：如何控制滤波器的阻带衰减