FIR数字滤波器设计实验.docVIP

附录：实验3 FIR数字滤波器的设计 1、程序一： n=15; wn=[0.3,0.5] b=fir1(n,wn,Hanning(n+1)); [h,w]=freqz(b,1,512); freqzplot(h,w); title(n=15 Hanning); figure; n=45; wn=[0.3,0.5] b=fir1(n,wn,Hanning(n+1)); [h,w]=freqz(b,1,512); freqzplot(h,w); title(n=45 Hanning); 图3-1 （a） N=15时，汉宁窗设计的线性相位带通滤波器 图3-1 （b） N=45时，汉宁窗设计的线性相位带通滤波器 2、程序二： n=15; wn=[0.3,0.5] b=fir1(n,wn,boxcar(n+1)); [h,w]=freqz(b,1,512); freqzplot(h,w); title(n=15 boxcar); figure; n=45; b=fir1(n,wn,boxcar(n+1)); [h,w]=freqz(b,1,512); freqzplot(h,w); title(n=45 boxcar); figure; n=15; b=fir1(n,wn,blackman(n+1)); [h,w]=freqz(b,1,512); freqzplot(h,w); title(n=15 blackman); figure; n=45; b=fir1(n,wn,blackman(n+1)); [h,w]=freqz(b,1,512); freqzplot(h,w); title(n=45 blackman);  图3-2（a）矩形窗设计的线性相位带通滤波器 图3-2（b）汉宁窗设计的线性相位带通滤波器图3-2（c）布莱克曼窗设计的线性相位带通滤波器 3、程序三： n=40;wn=[0.2,0.4,0.6,0.8]; b=fir1(n,wn,kaiser(n+1,4)); [h,w]=freqz(b,1,512); freqzplot(h,w); title(kaiser beta=4);figure; b=fir1(n,wn,kaiser(n+1,6)); [h,w]=freqz(b,1,512); freqzplot(h,w); title(kaiser beta=6);figure; b=fir1(n,wn,kaiser(n+1,10)); [h,w]=freqz(b,1,512); freqzplot(h,w); title(kaiser beta=10); 图3-3（a）β=4 凯塞窗设计的线性相位滤波器 图3-3（b）β=6 凯塞窗设计的线性相位滤波器 图3-3（c）β=10 凯塞窗设计的线性相位滤波器 4、程序四： n=40; w=[0 0.175 0.2 0.225 0.375 0.4 0.425 0.575 0.6 0.625 0.775 0.8 0.825 1]; m=[0 0 0.5 1 1 0.5 0 0 0.5 1 1 0.5 0 0]; b=fir2(n,w,m); [h,w]=freqz(b,1,512); freqzplot(h,w); title(?μ?ê2é?ù·¨); 图3-4 频率采样法设计线性相位滤波器