实验4用窗函数法设计FIR数字滤波器.docVIP

用窗函数法设计FIR数字滤波器 实验目的? 加深对窗函数法设计FIR数字滤波器的基本原理的理解 学习用MATLAB语言的窗函数法编写设计FIR数字滤波器的程序 3、了解MATLAB有关窗函数法设计的常用子函数? ? 二、实验涉及的MATLAB子函数? 1、boxcar：矩形窗? 2、triang：三角窗? 3、bartlett：巴特利特窗? 4、hamming：哈明窗? 5、hanning：汉宁窗 ?6、blackman：布莱克曼窗? chebwin：切比雪夫窗? 8、kaiser：凯瑟窗? 9、firl：基于窗函数的FIR数字滤波器设计——标准频率响应，以经典方法实现加窗线性相位FIR滤波器设计，可设计出标准的低通、带通、高通和带阻滤波器。? ? 三、实验原理? 1、运用窗函数法设计FIR数字滤波器? 分母a0为1，其余ak全都为0时的一个特例。由于极点全部集中在零点，稳定和线性相位特性是FIR滤波器的突出优点，因此在实际中广泛使用。? FIR滤波器的设计任务是选择有限长度的h(n)，使传输函数H(ejw)满足技术要求。?用窗函数法设计FIR数字滤波器的基本步骤是：? 根据过渡带和阻带衰减设计指标选择窗函数类型，估算滤波器的阶数N； 2)由数字滤波器的理想频率响应H(ejw)求出其单位冲击响应hd(n)。 2、各种窗函数特性的比较? 3、用窗函数设计FIR数字低通滤波器?4、用窗函数法设计FIR数字高通滤波器?5、用窗函数法设计FIR数字带通滤波器?6、用窗函数法设计FIR数字带阻滤波器? ? 四、实验内容? 选择合适的窗函数设计FIR数字低通滤波器，要求：通带ωp=0.2π，Rp=0.05dB；阻带ωs=0.3π，As=40dB。描绘实际滤波器的脉冲响应、窗函数及滤波器的幅频响应曲线和相拼相应曲线。? ? 五、实验步骤? 1、打开MATLAB软件，选择“File/New”创建一个新的文件；? 2、编程： wp=0.2*pi; ws=0.3*pi; deltaw=ws-wp; N0=ceil(6.6*pi/deltaw); N=N0+mod(N0+1,2); windows=(hamming(N)); wc=(ws+wp)/2; hd=ideal_lp(wc,N); b=hd.*windows; [db,mag,pha,grd,w]=freqz_m(b,1); n=0:N-1;dw=2*pi/1000; Rp=-(min(db(1:wp/dw+1))); As=-round(max(db(ws/dw+1:501))); subplot(2,2,1),stem(n,b); axis([0,N,1.1*min(b),1.1*max(b)]);title(实际脉冲响应)； xlabel(n);ylabel(h(n)); subplot(2,2,2),stem(n,windows); axis([0,N,0,1.1]);title(窗函数特性); xlabel(n);ylabel(wd(n)); subplot(2,2,3),stem(w/pi,db); axis([0,1,-80,10]);title(幅度频率响应); xlabel(频率（单位：\pi）);ylabel(H(e^{j\omega})); subplot(2,2,4),stem(w/pi,pha); axis([0,1,-80,10]);title(相位频率响应); 3）将上述程序在MATLAB中运行，并对实验结果进行分析?? 六、实验结果? 实际滤波器的脉冲响应、窗函数及滤波器的幅频响应曲线和相拼相应曲线如图所示：