滤波器的性能指标及MATLAB仿真.docVIP

滤波器的性能指标及MATLAB仿真 常见的FIR数字滤波器有带阻、带通、低通和高通四种，这里选择低通滤波器为例，具体性能指标包括：通带截止频率，阻带截止频率，实际通带波动，最小阻带衰减,抽样频率。 根据上述的性能指标，我们采样窗函数设计法，常见的窗函数有很多，例如汉明窗、布莱克曼窗及凯塞窗等等，MATLAB提供了几个子程序用于实现这些窗函数，以下为子程序： w=boxcar(M)在数组w中产生M点的矩形窗函数； w=triang(M)在数组w中产生M点的Bartlett(三角形)窗函数； w=hanning(M)在数组w中产生M点的Hanning窗函数； w=hamming(M)在数组w中产生M点的Hamming窗函数； w=Blackman(M)在数组w中产生M点的Blackman窗函数； w=Kaiser(M,beta)在数组w中产生beta值的M点矩形窗函数； 利用这些子程序就能用MATLAB设计基于窗函数法的FIR数字滤波器，当然也要求一个理想低通脉冲响应。因此，有一个单一的子程序创建是方便的，该子程序如下给出： Function hd=ideal_lp(wc,M); %Ideal LowPass filter computation %--------------------------------- %[hd]=ideal_lp(wc,M) %hd=ideal impulse response between 0 to M-1 %wc=cutoff frequency in radians %M=length of the ideal filter % alpha=(M-1)/2; n=[0:1:(M-1)]; m=n-alpha+eps; %add smallest number to avoid divide by zero hd=sin(wc*m)./ (pi*m); 接下来进行滤波器的程序设计，在窗函数的选择方面，选择凯塞窗作为对信号进行截断的窗函数，凯塞窗是一种适应性比较强的窗，其窗函数的表示式为： ，0≤n≤N-1 (3-1) 其中是第一类变形零阶贝塞尔函数，是一个可自由选择的参数，它可以同时调整主瓣宽度与旁瓣电平，越大，则窗越窄，而频谱的旁瓣越小，但主瓣宽度也相应增加。 采用凯塞窗进行FIR滤波器设计的程序如下： wp=0.2*pi; ws=0.3*pi; As=50; tr_width=ws-wp; M=ceil((As-7.95)/(14.36*tr_width/(2*pi))+1)+1 M=61 n=[0:1:M-1]; beta=0.1102*(As-8.7) Beta=4.5513 wc=(ws+wp)/2; hd=ideal_lp(wc,M); w_kai=(Kaiser(M,beta))’; h=hd.*w_kai; [db,mag,pha,grd,w]=freqz_m(h,[1]); Delta_w=2*pi/1000 As=-round(max(db(ws/delta_w+1:1:501))) %Min stopband Attenuation As=52 %Plots subplot(1,1,1)； subplot(2,2,1); stem(n,hd); title(‘Ideal Impulse Response’) axis([0 M-1 -0.1 0.3]); xlabel(‘n’); ylabel(‘hd(n)’) subplot(2,2,2); stem (n,w_kai); title(‘Kaiser Window’) axis([0 M-1 0 1.1]); xlabel(‘n’); ylabel(‘w(n)’) subplot(2,2,3); plot(w/pi,db); title(‘Magnitude Response in dB’); grid axis([0 1 -100 10]); xlabel(‘frequency in pi units’); ylabel(‘Decibels’) 此时，Kaiser窗的参数是和，实际的阻带衰减是。 通过MATLAB仿真后,可以得到低通滤波器的图如下： 图3-1 低通滤波器幅度响应(dB) 图3-2 低通滤波器相位图形 图3-3 低通滤波器实际脉冲响应 图3-4 Kaiser窗波形 3.3 对于采用不同窗函数的研究 依然按照上面的性能指标，这次采用汉明窗为截断信号的窗函数，汉明窗函数给出了比较小的过渡带，因此有较低的阶。这里，再给出汉明