设计滤波器.doc

第5章作业 5.1. 用窗函数法设计一个线性相位低通FIRDF，要求通带截止频率为π/4 rad，过渡带宽度为8π/51 rad，阻带最小衰减为45 dB。 （1） 选择合适的窗函数及其长度，求出h(n)的表达式。 ① 根据对阻带衰减及过渡带的指标要求， 选择窗函数的类型， 并估计窗口长度N。 由习题9中教材表7.2.2， 本题应选择哈明窗。 因为过渡带宽度Bt=8π/51， 所以窗口长度N为N≥6.6π/Bt=42.075， 取N=43。 窗函数表达式为 （2*） 用MATLAB画出损耗函数曲线和相频特性曲线。 wp=pi/4; Bt=8*pi/51;wc=wp+Bt/2; N0=ceil(6.6*pi/Bt); hn=fir1(N0-1, wc/pi, low, hamming(N0)); [h1,w1]=freqz(hn,1); f=w1/pi; figure(1); subplot(121); plot(f,20*log10(abs(h1))); grid;title(损耗函数曲线); xlabel(归一化频率/\pi);ylabel(幅度/dB); subplot(122); plot(f,angle(h1));grid; title(相频特性曲线); xlabel(归一化频率/\pi);ylabel(相位/\pi); figure(2); stem([0:N0-1],hn); title(单位脉冲响应);xlabel(n);ylabel(h(n)); 5.2 用频率采样法设计一个线性相位低通滤波器，N=32，Wc=pi/2,边缘上设一点过渡带|H(k)|=0.39,试求个点采样值H(k).并用MATLAB画出该滤波器的损耗曲线图。 因为π/2=（2π/32）×8，所以通带边界频率在k=8处，即以此点作为过渡带采样点。又因N为偶数，属于第二类线性相位滤波器，幅度函数在0到2π上以π为中心呈奇对称分布，应满足Hk=-HN-k ,所以频率采样值的幅度为： 频率采样值的相位为： 频率采样值H(k)为： N=32;Hk=[ones(1,8) 0.39 zeros(1,15) -0.39 -ones(1,7)]; k=0:N-1; hn=real(ifft(Hk.*exp(-j*pi*(N-1)*k/N))); [H w]=freqz(hn,1); f=w/pi; plot(f,20*log10(abs(H)));grid; title(损耗函数曲线);xlabel(归一化频率/\pi); ylabel(幅度/dB); 5.3*. 调用MATLAB工具箱函数fir1设计线性相位带通FIR滤波器。要求通带截止频率为0.55π和0.7π，阻带截止频率为0.45π和0.8π，通带最大衰减为0.15 dB，阻带最小衰减为40 dB， 显示所设计的单位脉冲响应h(n)的数据，并画出损耗函数曲线。 wpl=0.55*pi; wpu=0.7*pi; wsl=0.45*pi; wsu=0.8*pi; rs=40; wc=[(wpl+wsl)/2/pi,(wpu+wsu)/2/pi]; Bt=wpl-wsl; N0=ceil(6.2*pi/Bt); N=N0+mod(N0+1, 2); hn=fir1(N-1, wc, bandpass, hanning(N)); [h1,w1]=freqz(hn,1);f=w1/pi; figure(1); plot(f,20*log10(abs(h1)));grid; title(损耗函数曲线); xlabel(归一化频率/\pi);ylabel(幅度/dB); figure(2); stem(0:N-1,hn); title(单位脉冲响应);xlabel(n);ylabel(h(n)); 运行结果：N=63，h(n)的长度为63 5.4*. 调用remezord和remez函数完成题5.3题所给技术指标的滤波器的设计，并比较设计结果（主要比较滤波器阶数的高低和幅频特性）。 f=[0.45,0.55,0.7,0.8];m=[0,1,0];rp=0.15;rs=40; dat1=(10^(rp/20)-1)/(10^(rp/20)+1);dat2=(10^(-rs/20)); rip=[dat2,dat1,dat2]; [M,fo,mo,w]=remezord(f,m,rip); hn=remez(M,fo,mo,w); [h1,w1]=freqz(hn,1);f=w1/pi; figure(1); plot(f,20*log10(abs(h1)));grid; title(损耗函数曲线); xlabel(归一化频率/\pi);ylabel(幅度/dB);