%《數(shù)字信號處理教程——MATLAB釋義與實(shí)現(xiàn)》第七章例7.2.1程序hc721
% 幅頻特性與符幅特性的區(qū)別
% 電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%
clear
b=[-1,2,4,2,-1];a=1;            %  給定濾波器系數(shù)向量
[H,w]=freqz(b,a);               %  求頻率特性
subplot(2,2,1),plot(w,abs(H))   %  求并畫出幅頻特性
grid on
subplot(2,2,3),plot(w,angle(H)) %  求并畫出幅頻特性
grid on
A=-2*cos(2*w)+4*cos(w)+4;       %  求符幅特性
theta=-2*w;                     %  求線性相位特性
subplot(2,2,2),plot(w,A)        %  畫出符幅特性
line([0,4],[0,0]),grid on
subplot(2,2,4),plot(w,theta)    %  畫出線性相位特性
grid on
set(gcf,'color','w');           % 置圖形背景色為白