基于MATLAB的FIR数字滤波器设计资料.docVIP

基于MATLAB的FIR数字 滤波器设计 系别： 班级： 姓名： 学号： 1.实验目的： 加深对低通滤波器的设计过程的理解； 进一步学习用窗函数法和等波纹法实现数字滤波器设计的方法； (3) 比较两种方法所得出的结果的不同之处。 实验内容： 利用MATLAB编程设计一个数字带通滤波器，指标要求如下： 通带边缘频率：，，通带峰值起伏：; 阻带边缘频率：，，最小阻带衰减： ; 分别用窗函数法和等波纹滤波器法设计两种FIR数字滤波器。 3.实验原理： （1）用窗函数法实现： 调用函数[n,wn,bta,ftype]=kaiserord（f，a，dev，fs） 参数：f=[0.3 0.45 0.65 0.8]为对应数字频率；，，，， a=[0 1 0]为由f指定的各个频带上的幅值向量， 一般只有0和1表示；和f长度关系为（2*a的长度）—2=（f的长度）,devs=[0.01 0.1087 0.01]用于指定各个频带输出滤波器的频率响应与其期望幅值之间的最大输出误差或偏差，长度与a相等，计算公式：阻带衰减误差=通带波动衰减误差=，fs缺省值为2HZ。 （2）用等波纹法实现： 调用函数: [n,fpts,mag,wt]=remezord(f，a，dev); f=[0.3 0.45 0.65 0.8]; a=[0 1 0]; dev=[0.01 0.1087 0.01]; 其含义同函数[n,wn,bta,ftype]=kaiserord（f，a，dev，fs)中的参数相同。 4.实验要求： 给出FIR数字滤波器的冲激响应，绘出它们的幅度和相位频响曲线。 实验程序清单： （1）用窗函数法实现： clc; clear; clear all; [n,wn,bta,ftype]=kaiserord([0.3 0.45 0.65 0.8],[0 1 0],[0.01 0.1087 0.01]);%用kaiserord函数估计出滤波器阶数n和beta参数 h1=fir1(n,wn,ftype,kaiser(n+1,bta),noscale); [hh1,w1]=freqz(h1,1,256); figure(1) subplot(2,1,1) plot(w1/pi,20*log10(abs(hh1))) Grid; xlabel(归一化频率w);ylabel(幅度/db); subplot(2,1,2) plot(w1/pi,angle(hh1)) grid xlabel(归一化频率w);ylabel(相位/rad); h1; 所得曲线： 若用freqz(h1,1,256)，则为： 运行结果： h1 = Columns 1 through 8 0.0041 0.0055 -0.0091 -0.0018 -0.0056 -0.0000 0.0391 -0.0152 Columns 9 through 16 -0.0381 0.0077 -0.0293 0.0940 0.0907 -0.2630 -0.0517 0.3500 Columns 17 through 24 -0.0517 -0.2630 0.0907 0.0940 -0.0293 0.0077 -0.0381 -0.0152 Columns 25 through 31 0.0391 -0.0000 -0.0056 -0.0018 -0.0091 0.0055 0.0041 (2)用等波纹法设计： clc; clear; clear all; [n,fpts,mag,wt]=remezord([0.3 0.45 0.65 0.8],[0 1 0],[0.01 0.1087 0.01]); h2=remez(n,fpts,mag,wt);%设计出等波纹滤波器 [hh2,w2]=freqz(h2,1,256); figure(2) subplot(2,1,1) plot(w2/pi,20*log10(abs(hh2))) grid xlabel(