数字电路滤波器MATLAB课程设计.doc

数字信号处理课程设计报告 一、 数字信号处理课程设计目的 通过课程设计，主要达到以下目的： 1、掌握Matlab语言基础及使用入门知识，掌握使用Matlab语言及其工具箱进行基本信号的分析与处理。加深对数字信号处理理论方面的理解。 3、了解和掌握用MATLAB实现IIR和FIR滤波器的设计方法、过程，为以后的设计打下良好基础。 二、设计工具 1、数字信号处理课程设计所需要的硬件工具 PC机。 2、数字信号处理课程设计所需要的软件 MATLAB软件。 三、设计题目程序及结论 IIR模拟滤波器设计题目：（巴特沃斯带通模拟滤波器） 设计一带通滤波器，采样频率为1000HZ，通带为60HZ到200HZ，通带波纹Rp≤3dB,截止带宽50Hz,且在通带两边衰减为40dB （1）设计原理： 输入信号和输出信号均为连续时间信号，冲击响应也是连续的滤波器，称为模拟滤波器，由于IIR数字滤波器是在已知的低通模拟滤波器基础上设计的，主要包括巴特沃斯低通滤波器、切比雪夫低通滤波器、椭圆低通滤波器，因此把这些模拟低通滤波器称为滤波器原型。巴特沃斯滤波器具有单调下降的幅频特性，通带和阻带幅频都比较平坦。 （2）IIR模拟滤波器设计源程序： Wp=[60 200] /500; Ws=[50 250]/500; Rp=3; Rs=40; [n, Wn]=buttord(Wp, Ws, Rp, Rs) [b,a]=butter(n, Wn); freqz(b,a,128,1000) title( n=16 Butterworth Bandpass Filter ) （3）IIR模拟滤波器设计结果及波形（模拟带通滤波器巴特沃斯） n = 16 Wn =0.1198 0.4005 IIR数字滤波器设计题目：（巴特沃斯低通数字滤波器） 采用双线性变换法设计一个巴特沃斯低通数字滤波器，要求： Wp=0.25*pi,Rp=1dB；Ws=0.4*pi,As=15dB,滤波器采样频率Fs=100Hz。 （1）设计原理：双线性变换法是将整个s平面映射到整个z平面，它克服了脉冲响应不变法从s平面到z平面的多值映射的缺点，消除了频谱混叠现象。但其在变换过程中产生了非线性的畸变，在设计IIR数字滤波器的过程中需要进行一定的预修正。 （2）IIR数字滤波器设计源程序： wp=0.25*pi; ws=0.4*pi; Rp=1;As=15; Fs=100;T=1/Fs; Wp=(2/T)*tan(wp/2); Ws=(2/T)*tan(ws/2); [n,Wc]=buttord(Wp,Ws,Rp,As,s); [z0,p0,k0]=buttap(n); ba=k0*real(poly(z0)); aa=real(poly(p0)); [ba1,aa1]=lp2lp(ba,aa,Wc); [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,1,1),plot(w/pi,abs(H)); ylabel(|H|);title(幅度响应);axis([0,1,0,1.1]); subplot(2,1,2),plot(w/pi,dbH);title(幅度响应（dB）) ylabel(dB);xlabel(频率(\pi));axis([0,1,-40,5]) （3）IIR巴特沃斯型数字低通滤波器设计波形 FIR数字滤波器设计题目：（汉明窗低通FIR数字滤波器） 利用汉明窗设计一个数字低通滤波器，其中滤波器的技术指标为? Wp=0.25? Ws=0.35?。 （1）设计原理：FIR滤波器的冲激响应就是系统函数各次项的系数，所以设计FIR滤波器的方法之一就是：从时域出发，截取有限长的一段冲激响应作为H（z）系数，冲激响应长度N就是系统函数H（z）的阶数。只要N足够长，并且截取的方法合理，总能够满足频域的要求，这就是FIR滤波器的窗函数设计法。 （2）FIR数字滤波器设计源程序 wp=0.25*pi;ws=0.35*pi; wdelta=ws-wp; N=ceil(6.6*pi/wdelta) Nw=N; wc=(wp+ws)/2; n=0:N-1; alpha=(N-1)/2; m=n-alpha+eps;