%《MATLAB及在電子信息課程中的應用》第七章繪圖子程序
% 由h(n)計算其分貝幅特性20log(|H(ejw)|)并繪圖的函數hnwplot
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
% a=hnwplot(hn,i)
% ------------------------------------------
% h(n)為給出的任意序列
% i 為繪制時間序列的子圖位置，全圖取3x2個子圖
% 在i+1子圖上將給出此時間序列的分貝幅頻特性
%
function a=hnwplot(hn,i)
N=length(hn);
n=0:N-1;
hw=fft(hn,512);k=0:511;w=2*k*pi/512;    % 求頻譜及頻率自變量序列
subplot(3,2,i);stem(n,hn,'.');          % 在第i個子圖上畫出時間序列
set(gcf,'color','w')                    % 將圖形背景色設為白色
% axis([0,N-1,min(hn)-0.05,max(hn)+0.1]); 
xlabel('n');ylabel('h(n)')
line([0,N-1],[0,0])
subplot(3,2,i+1)                        % 在第i+1個子圖上
% plot(w,abs(hw));                      % 畫出頻譜特性
plot(w,20*log10(abs(hw)));              % 畫出分貝頻譜特性
grid;title('幅度特性(dB)');
xlabel('ω/π');ylabel('Magnitude(dB)')