r1(N-1,wc,hanning(N));Р [h1,w1]=freqz(hn,1);Р subplot(4,3,2);Р plot(w1/pi,20*log10(abs(h1)));grid onР axis([0,1,-80,5])Р subplot(4,3,5);Р plot(w1/pi,20*log10(angle(h1)));grid onР h=conv(y,hn);Р [ft,Wt]=freqz(h,512);Р subplot(4,3,3);Р plot(Wt/pi,abs(ft));Р subplot(4,3,6);Р plot(Wt/pi,angle(ft));Р汉宁窗高通滤波器的callback为:Р只要把hn=fir1(N-1,wc,hanning(N));改为hn=fir1(N-1,wc,'high',hanning(N)); 即可Р其他窗函数只要改变此句(hn=fir1(N-1,wc,hanning(N));)里的窗口类型和ftype即可。Р这6种窗函数的产生函数和调用格式为:Р矩形窗: wn=boxcar(N)Р三角窗函数: wn=bartlett(N)Р汉宁窗函数: wn=hanning(N)Р哈明窗函数: wn=hamming(N)Р布莱克曼窗函数: wn=blackman(N)Р凯塞窗函数: wn=kaiser(N)Р滤波器频谱分析设计: РGUI界面滤波器的幅频、相频特性曲线输出部分设计如下图:Р当ap=3,as=30,fp=0.1,fs=0.3时:Р(1)巴特沃斯低通滤波器的频谱分析为:Р(2)切比雪夫低通滤波器的频谱分析为:Р当ap=3,as=30,fp=0.3,fs=0.1时:Р(1)巴特沃斯高通滤波器的频谱分析为:Р(2)切比雪夫高通滤波器的频谱分析为:Р当fp=0.1,fs=0.3时:Р(1)汉宁窗低通滤波器的频谱分析如下:
全文预览
