hc721.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第七章例7.2.1程序hc721
% 幅频特性与符幅特性的区别
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear
b=[-1,2,4,2,-1];a=1;            %  给定滤波器系数向量
[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');           % 置图形背景色为白