>> fs=10000;
fpl=1000;fph=3000;
fsl=500;fsh=3500;
wpl=2*fpl/fs;wph=2*fph/fs;
wsl=2*fsl/fs;wsh=2*fsh/fs;
wc=[wpl,wph];
B=wpl-wsl;%过渡带的宽度
wdl=(wsl+wpl)/2;wdh=(wph+wsh)/2;%低频截止频率和高频截止频率
wp=wdh-wdl;%通带的宽度
N=ceil(8/B);
n=0:N-1;
window=Hamming(N);
%///////////////////////////////////////////////////////////////////////////
[h1,w]=freqz(window,1);
figure;
stem(window);
xlabel('n');
title('海明窗函数');
%///////////////////////////////////////////////////////////////////////
figure;
plot(w/pi,20*log(abs(h1)/abs(h1(1))));
grid;
xlabel('w/pi');
ylabel('幅度(dB)');
title('海明窗函数的频谱');
hn = fir1(N-1,wc,Hamming(N));
%/////////////////////////////////////////////////////////////////////////
本内容试读结束,登录后可阅读更多
下载后可阅读完整内容,剩余1页未读,立即下载