%《數字信號處理教程——MATLAB釋義與實現》第七章例7.3.3程序hc733
% 三角窗累積符幅特性
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
clear
N=51; n=1:N; wd=triang(N);		% 給出三角窗參數，產生窗函數
[A,omega]=amplres(wd');			% 調用amplres，求符幅函數，變wd為行向量
Ad=[fliplr(A),A(2:end)];		% 將符幅函數對稱擴展至負頻率區
omegad=[-fliplr(omega),omega(2:end)]; % 將頻率對稱擴展至負頻率區
Gd=cumtrapz(Ad)*(omega(2)-omega(1))/pi; % 積分求累積符幅函數Gd
AdB=20*log10(abs(Ad));			% 將符幅函數轉換為分貝
GdB=20*log10(Gd); 				% 將累積符幅函數轉換為分貝
subplot(2,2,1),plot(omegad,Ad)  % 繪圖
subplot(2,2,2),plot(omegad,Gd)
subplot(2,2,3),plot(omegad,AdB)
axis([-pi,pi,-50,50]),grid on
subplot(2,2,4),plot(omegad,GdB)       % 繪分貝幅特性
axis([-pi,pi,-50,10]),grid on
set(gcf,'color','w');                % 置圖形背景色為白