hc841a.m

这是数字信号处理——matlab释义与实现一书的源代码

%《数字信号处理教程——MATLAB释义与实现》第八章例8.4.1程序hc841a
% 由低通模拟变换为其它类型模拟频带的变换示例
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
num=1;den=[1,2,2,1];                        % 给定归一化巴特沃斯模拟滤波器分子分母系数向量
OmegaC=10*2*pi; B=2*2*pi;                   % 边界频率和带宽
[numT1,denT1] = lp2lp(num,den,OmegaC)	    % 低通模拟变低通模拟	
[numT2,denT2] = lp2hp(num,den,OmegaC)	    % 低通模拟变高通模拟
[numT3,denT3] = lp2bp(num,den,OmegaC,B)	    % 低通模拟变带通模拟
[numT4,denT4] = lp2bs(num,den,OmegaC,B)	    % 低通模拟变带阻模拟
[Ha1,Omega1]=freqs(numT1,denT1);            % 低通模拟频率特性
[Ha2,Omega2]=freqs(numT2,denT2);            % 高通模拟频率特性
[Ha3,Omega3]=freqs(numT3,denT3);            % 带通模拟频率特性
[Ha4,Omega4]=freqs(numT4,denT4);            % 带阻模拟频率特性
figure(1)
subplot(2,2,1),plot(Omega1,abs(Ha1)),axis([0,100,0,1.1]),grid on        % 画模拟幅频特性图
title('变换后的模拟低通特性'),xlabel('\Omega'),ylabel('|Ha1(\Omega)|')
subplot(2,2,2),plot(Omega2,abs(Ha2)),axis([0,100,0,1.1]),grid on
title('变换后的模拟高通特性'),xlabel('\Omega'),ylabel('|Ha1(\Omega)|')
subplot(2,2,3),plot(Omega3,abs(Ha3)),axis([0,100,0,1.1]),grid on
title('变换后的模拟带通特性'),xlabel('\Omega'),ylabel('|Ha1(\Omega)|')
subplot(2,2,4),plot(Omega4,abs(Ha4)),axis([0,100,0,1.1]),grid on
title('变换后的模拟带阻特性'),xlabel('\Omega'),ylabel('|Ha1(\Omega)|')
set(gcf,'color','w')                                    % 置图形背景色为白
T=0.001;                                     % 采样周期（hc841a和hc841a只改了这个参数）
[bd1,ad1]=bilinear(numT1,denT1,1/T)         % 低通模拟变换为低通数字
[bd2,ad2]=bilinear(numT2,denT2,1/T)         % 高通模拟变换为高通数字
[bd3,ad3]=bilinear(numT3,denT3,1/T)         % 带通模拟变换为带通数字
[bd4,ad4]=bilinear(numT4,denT4,1/T)         % 带阻模拟变换为带阻数字
[H1,w1]=freqz(bd1,ad1);                     % 低通数字频率特性
[H2,w2]=freqz(bd2,ad2);                     % 高通数字频率特性
[H3,w3]=freqz(bd3,ad3);                     % 带通数字频率特性
[H4,w4]=freqz(bd4,ad4);                     % 带阻数字频率特性
figure(2)
subplot(2,2,1),plot(w1,abs(H1)),axis([0,5,0,1.1])                        % 画数字幅频特性图
title('变换后的数字低通特性'),xlabel('\omega'),ylabel('|H1(\omega)|'),grid on
subplot(2,2,2),plot(w2,abs(H2)),axis([0,5,0,1.1])
title('变换后的数字高通特性'),xlabel('\omega'),ylabel('|H2(\omega)|'),grid on
subplot(2,2,3),plot(w3,abs(H3)),axis([0,5,0,1.1])
title('变换后的数字带通特性'),xlabel('\omega'),ylabel('|H3(\omega)|'),grid on
subplot(2,2,4),plot(w4,abs(H4)),axis([0,5,0,1.1])
title('变换后的数字带阻特性'),xlabel('\omega'),ylabel('|H4(\omega)|'),grid on
set(gcf,'color','w')                                    % 置图形背景色为白