基于matlab的海明窗带阻滤波器的设计

fs=15000; wc1=0.22*pi; wc2=0.72*pi; N=31; n=[0:1:30]; %窗函数的宽度 r=(31-1)/2; %线性相位的移位 hd=(sin(pi*(n-r))+sin(wc1*(n-r))-sin(wc2*(n-r)))./(pi*(n-r+eps)); %理想带阻滤波器单位冲击响应 h=fir1(30,[wc1/pi wc2/pi],'stop',hanning(31)); %调用汉宁窗 [h1,w1]=freqz(h,1); %求数字滤波器的频响 figure(1) subplot(2,1,1); plot(w1/pi,20*log10(abs(h1))); %做滤波器的特性曲线 axis([0,1,-110,0]); grid xlabel('归一化频率w'); ylabel('幅度 '); title('幅频特性'); subplot(2,1,2) plot(w1,180/pi*(angle(h1))); %画出频率响应伯德图 xlabel('频率/hz') ; ylabel('相位/^o'); grid on; %显示网格线 figure(2) N=100; t=(0:N)/fs; %采样间隔和时间序列 s = sin(2*pi*t*800)+sin(2*pi*t*3000)+sin(2*pi*t*7000); %多频率输入信号 sf = filter(h,1,s); %对信号进行滤波 subplot(2,1,1); stem(t,s,'.') grid xlabel('时间/s '); title('输入信号 ') s1=fft(s,128); %对输入信号进行快速福利叶变换 f1=(0:length(s1)-1)'*fs/length(s1); %进行对应的频率转换 g1=abs(s1); subplot(2,1,2); stem(f1,g1,'.') %输入信号的离散信号 grid xlabel('w');ylabel('幅度'); title(' 输入信号频谱'); figure(3) subplot(2,1,1); stem(t,sf,'.') grid xlabel('时间/s '); title(' 输出信号'); s2=fft(sf,128); f2=(0:length(s2)-1)'*fs/length(s2); g2=abs(s2); subplot(2,1,2); stem(f2,g2,'.') grid xlabel('w');ylabel('幅度 '); title(' 输出信号频谱');