hc962a.m

matlab数字信号处理代码

%《数字信号处理教程——MATLAB释义与实现》第四章例9.6.2程序hc962
% 一阶缓冲滤波器的传递函数
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
K=0.5;alpha0=[0.1,0.2:0.2:0.8];     % 设定一个K 和四个alpha值
for k=1:5
    alpha=alpha0(k);
    b1=[alpha,-1]; a1=[1,-alpha];   % 列出A1(z)分子分母系数
    b = 0.5*(1+K)*a1+0.5*(1-K)*b1;  % 均衡器分子系数
    [H,w]=freqz(b,a1);              % 求均衡器频率特性
    semilogx(w,20*log10(abs(H))),hold on    % 画幅频特性曲线
end
%axis([0.01*pi,pi,0,10])
set(gcf,'color','w')                % 置图形背景色为白