fig72_1_2.m

高西全2009年最新版<数字信号处理>源码程序

%《数字信号处理--原理、实现及应用》第六章 图7.2.1绘图程序:fig721.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
close all;
wc=1/4;N=30;alph=(N-1)/2;
n=-10:40;
hdn=sin(wc*pi*(n-alph))./(pi*(n-alph));
subplot(3,1,1);stem(n,hdn,'.')
ylabel('h_d(n)');axis([-10,40,-0.1,0.3])
wrn=boxcar(N);
hn=fir1(N-1,wc,boxcar(N));
n=0:N-1;
subplot(3,1,2);stem(n,wrn,'.')
ylabel('R_N(n)');axis([-10,40,0,1.3])
subplot(3,1,3);stem(n,hn,'.');ylabel('h(n)')
line([0,50],[0,0])
axis([-10,40,-0.1,0.3]);xlabel('n')
Hk=fft(hn,512);
w=(0:511)*2*pi/512;
figure(2);subplot(2,2,2),plot(w/pi,20*log10(abs(Hk)))
xlabel('ω/π');ylabel('幅度/dB')
title('(b) 损耗函数曲线')
axis([0,1,-60,10])
grid on;
subplot(2,2,1),plot(w/pi,abs(Hk))
xlabel('ω/π');ylabel('幅度')
title('(a) 幅频响应函数曲线')
axis([0,1,0,1.2])
grid on;