%《MATLAB在電子信息課程中的應用》第七章繪圖子程序
% h(n),|H(ejw)|計算與繪圖函數hnwplot
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

function a=hnwplot(hn,k)
N=length(hn);
n=0:N-1;
hw=fft(hn,512);w=0:511;w=2*w/512;
subplot(3,2,k+1)
plot(w,20*log10(abs(hw)));grid;title('幅度特性(dB)');
axis([0,max(w),-100,5]);
xlabel('ω/π');ylabel('Magnitude(dB)')
subplot(3,2,k);stem(n,hn,'.');
axis([0,N-1,min(hn)-0.05,max(hn)+0.1]);xlabel('n');ylabel('h(n)')
line([0,N-1],[0,0])