hc851.m

书籍中包含的第三部分实例的matlab代码

%《数字信号处理教程——MATLAB释义与实现》第八章例8.5.1程序hc851
% 由低通数字变换为高通数字的频带变换示例
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear
wc=0.3*pi; OmegaClp=10*2*pi;    % 给出原始参数
T=0.01;                         % 
wclp=2*atan(OmegaClp*T/2)      % 求出原低通数字截止频率
alpha = -(cos((wclp+wc)/2))/(cos((wclp-wc)/2))  % 低通-高通频带变换参数
Nz = -[alpha,1]; Dz = [1,alpha];        % 低通高通变换分子分母系数向量
bdlp=[0.0167, 0.0501, 0.0501,  0.0167]; % 原低通数字滤波器分子分母系数向量
adlp=[1, -1.7977, 1.2211, -0.2898];
[bhp,ahp] = zmapping(bdlp,adlp,Nz,Dz)   % 用zmapping求变换后的数字高通系数向量
[Hlp,wlp]=freqz(bdlp,adlp);             % 求低通滤波器的频率特性
subplot(1,2,1),plot(wlp,abs(Hlp))   %  绘制低通滤波器幅频特性
title('(a)低通数字滤波器幅频特性')
xlabel('\omegalp'),ylabel('|Hlp|')
grid on,
axis([0,pi,0,1.1])
set(gcf,'color','w')                 % 置图形背景色为白
[H,w]=freqz(bhp,ahp);                % 求高通滤波器的频率特性
subplot(1,2,2),plot(w,abs(H)),   %  绘制高通滤波器幅频特性
title('(b)变换后高通的幅频特性')
xlabel('\omegahp'),ylabel('|Hhp|')
grid on,axis([0,pi,0,1.1])