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(' 输出信号频谱');
