数字信号处理课程设计报告精要.docVIP

基本原理： 双线性变换法是将整个s平面映射到整个z平面，其映射关系为 　　或 双线性变换法克服了脉冲响应不变法从s平面到z平面的多值映射的缺点，消除了频谱混叠现象。但其在变换过程中产生了非线性的畸变，在设计IIR数字滤波器的过程中需要进行一定的预修正。 题目1.1涉及的MATLAB子函数： 1.Bilinear：调用格式[numd,dend]=bilinear(num,den,Fs);将模拟域传递函数变换为数字域传递函数,Fs为取样频率。 2.buttord 3. butter 设计流程： (1)输入给定的数字滤波器设计指标； (2)根据公式W＝(2/T)*tan(wp/2)进行预修正，将数字滤波器指标转换成模拟滤波器设计指标； (3)确定模拟滤波器的阶数和3db截止频率； (4)计算巴特沃斯模拟低通滤波器的系统传递函数； (5) 用双线性变换法将模拟滤波器转换为数字滤波器 (6) 计算滤波器的频率响应，绘制频率特性曲线。 程序清单： wp=0.25*pi; %通带边界频率 ws=0.4*pi;%阻带边界频率 Rp=1;As=15;%通带最大衰减 阻带最小衰减(db) ripple=10^(-Rp/20);Attn=10^(-As/20);%计算通带最大衰减 阻带最小衰减(非db) Fs=100;T=1/Fs;%采样频率及采样时间间隔 Omgp=(2/T)*tan(wp/2); Omgs=(2/T)*tan(ws/2);%由数字角频率计算模拟角频率 [n,Omgc]=buttord(Omgp,Omgs,Rp,As,s)%计算巴特沃斯模拟滤波器的阶数n及3db截止频率Omgc [ba1,aa1]=butter(n,Omgc,s);%计算巴特沃斯模拟滤波器的系统函数的系数向量ba1,aa1 [bd,ad]=bilinear(ba1,aa1,Fs)%将模拟域传递函数变换为数字域传递函数 [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,XTick,[0,0.25,0.4,1]); set(gca,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,XTick,[0,0.25,0.4,1]); set(gca,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,XTick,[0,0.25,0.4,1]); set(gca,YTick,[-50,-15,-1,0]);grid subplot(2,2,4),zplane(bd,ad); axis([-1.1,1.1,-1.1,1.1]);title(零极图); 运行结果： n =5 Omgc =103.2016 bd = 0.0072 0.0362 0.0725 0.0725 0.0362 0.0072 ad = 1.0000 -1.9434 1.9680 -1.0702 0.3166 -0.0392 sos = 1.0000 0.9962