%《數字信號處理教程——MATLAB釋義與實現》第八章例8.4.1程序hc841b
% 由低通模擬變換為其它類型模擬頻帶的變換示例
% 電子工業出版社出版  陳懷琛編著 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.1;                                      % 采樣周期（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')                                    % 置圖形背景色為白