3.6脉冲响应不变法­数字信号处理总结课件.pptVIP

%高通滤波器的设计 wp=1/(2*pi*5000);ws=1/(2*pi*1000);Ap=1;As=40; [N,Wc]=buttord(wp,ws,Ap,As,s); [num,den] = butter(N,Wc,s); disp(LP 分子多项式); fprintf(%.4e\n,num); disp(LP 分母多项式); fprintf(%.4e\n,den); [numt,dent] = lp2hp(num,den,1); disp(HP 分子多项式); fprintf(%.4e\n,numt); disp(HP 分母多项式); fprintf(‘%.4e\n’,dent); 例： 设计满足下列条件的模拟BW型高通滤波器 fp=5kHz, fs=1kHz, Ap?1dB, As ? 40dB。 Ap= 40.0000 As= 0.1098 0 1000 2000 3000 4000 5000 6000 -70 -60 -50 -40 -30 -20 -10 0 Frequency in Hz Gain in dB 例： 设计满足下列条件的模拟BW型高通滤波器 fp=5kHz, fs=1kHz, Ap?1dB, As ? 40dB 补充5. 模拟带阻滤波器的设计 模拟带阻滤波器MATLAB实现 [numt,dent] = lp2bs(num,den,W0,B) 模拟带通滤波器MATLAB实现 [numt,dent] = lp2bp(num,den,W0,B) 例：试设计一个满足下列指标的BW型带阻滤波器 Ap=1dB;As=20dB;wp1=10;wp2=30;ws1=19;ws2=21。 Ap=1;As=20;wp1=10;wp2=30;ws1=19;ws2=21; B=ws2-ws1;w0=sqrt(ws1*ws2); wLp1=B*wp1/(w0*w0-wp1*wp1); wLp2=B*wp2/(w0*w0-wp2*wp2); wLp=max(abs(wLp1),abs(wLp2)); [N,Wc]=buttord(wLp,1,Ap,As,s) [num,den] = butter(N,Wc,s ); [numt,dent]=lp2bs(num,den,w0,B); w=linspace(5,35,1000); h=freqs(numt,dent,w); plot(w,20*log10(abs(h))); w=[wp1 ws1 ws2 wp2]; set(gca,xtick,w);grid; h=freqs(numt,dent,w);A=-20*log10(abs(h)) 例：试设计一个满足下列指标的BW型带阻滤波器 Ap=1dB;As=20dB;wp1=10;wp2=30;ws1=19;ws2=21。 例：试设计一个满足下列指标的BW型带阻滤波器 Ap=1dB;As=20dB;wp1=10;wp2=30;ws1=19;ws2=21。 10 19 21 30 -100 -80 -60 -40 -20 0 补充6. 双线性变换法设计DF [numd,dend] = bilinear(num,den,Fs) num,den：AF分子、分母多项式的系数向量 Fs：抽样频率 numd,dend：DF分子、分母多项式的系数向量 例：利用AF-BW filter及双线性变换法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB %determine the DF filter [numd,dend]=bilinear(numa,dena,Fs) %plot the frequency response w=linspace(0,pi,1024); h=freqz(numd,dend,w); plot(w/pi,20*log10(abs(h))); axis([0 1 -50 0]);grid; xlabel(Normalized frequency); ylabel(Gain,dB); %computer Ap As of the designed filter 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)))); * 如何将模拟滤波器转变为数字滤波器? 1. 脉冲响应不变法 2. 双线性变换法 Wp,Ws wp,ws H(s) H(z) 频率 变换 设计模拟 滤波