数字信号处理之大作业.docx

数字信号处理教程之大作业郭航（2014212596）解：参考书p386解：参考书p396数字滤波器 通带截止频率π/5rad 通带最大衰减3dB 阻带截止频率3*π/5rad 阻带最小衰减20dB2.频率预畸变（T=2）设计其系统函数（p355-（7.5.24）） 所以 N=2 查表7.4得求数字滤波器的系统函数根据 表7.9 得附录一T=2; %设置采样周期为2 fs=1/T; %采样频率为周期倒数 Wp=0.2*pi/T; Ws=0.6*pi/T; %设置归一化通带和阻带截止频率 Ap=3; As=20; %设置通带最大和最小衰减 [N,Wc]=buttord(Wp,Ws,Ap,As,s); %调用butter函数确定巴特沃斯滤波器阶数 [B,A]=butter(N,Wc,s); %调用butter函数设计巴特沃斯滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应subplot(2,1,1); plot(W/pi,abs(hf)/abs(hf(1))); %绘出巴特沃斯模拟滤波器的幅频特性曲线 grid on; title(巴特沃斯模拟滤波器); xlabel(Frequency/Hz); ylabel(Magnitude); [D,C]= bilinear (B,A,fs); %调用双线性变换法Hz=freqz(D,C,W); %返回频率响应subplot(2,1,2); plot(W/pi,abs(Hz)/abs(Hz(1))); %绘出巴特沃斯数字低通滤波器的幅频特性曲线 grid on; title(巴特沃斯数字滤波器); xlabel(Frequency/Hz);ylabel(Magnitude);附录二wp=0.2*pi; ws=0.6*pi; Rp=3;As=20; ripple=10^(-Rp/20);Attn=10^(-As/20); Fs=0.5;T=1/Fs; Omgp=(2/T)*tan(wp/2); Omgs=(2/T)*tan(ws/2); [n,Omgc]=buttord(Omgp,Omgs,Rp,As,s) [ba1,aa1]=butter(n,Omgc,s); [bd,ad]=bilinear(ba1,aa1,Fs) [sos,g]=tf2sos(bd,ad) [H,w]=freqz(bd,ad); dbH=20*log10((abs(H)+eps)/max(abs(H)));subplot(2,2,1),plot(w/pi,abs(H)); ylabel(|H|);title(幅度响应);axis([0,1,0,1.1]); set(gca,XTickMode,manual,XTick,[0,0.25,0.4,1]); set(gca,YTickMode,manual,YTick,[0,Attn,ripple,1]);grid subplot(2,2,2),plot(w/pi,angle(H)/pi); ylabel(＼phi);title(相位响应);axis([0,1,-1,1]); set(gca,XTickMode,manual,XTick,[0,0.25,0.4,1]); set(gca,YTickMode,manual,YTick,[-1,0,1]);grid subplot(2,2,3),plot(w/pi,dbH);title(幅度响应(dB)); ylabel(dB);xlabel(频率(＼pi));axis([0,1,-40,5]); set(gca,XTickMode,manual,XTick,[0,0.25,0.4,1]); set(gca,YTickMode,manual,YTick,[-50,-15,-1,0]);gridsubplot(2,2,4),zplane(bd,ad); axis([-1.1,1.1,-1.1,1.1]);title(零极图);n = 2Omgc = 0.4363bd = 0.1053 0.2107 0.1053ad = 1.0000 -0.8958 0.3172sos = 1.0000 2.0000 1.0000 1.0000