%《數字信號處理教程——MATLAB釋義與實現》第八章圖8.4.6生成程序fgp846
% 模擬濾波器頻帶變換（低通變帶阻）
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
x0=logspace(-2,2,1000);         % 低通原型的頻率設定
x=[-fliplr(x0),x0];             % 變換為正負對稱頻率范圍
y1=1./x/2+sqrt((1./x/2).^2+1);  % 對頻率按(8.4.12a)式進行變換
y2=1./x/2-sqrt((1./x/2).^2+1);
plot(x,y1,x,y2,'-.'),hold on    % 畫出全頻域的映射曲線
set(gcf,'color','w')            % 置圖形背景色為白
x1=x(find(x<=1&x>=-1));         % 低通原型的通帶范圍（設帶寬為1）
y1=1./x1/2+sqrt((1./x1/2).^2+1);% 對通帶頻率范圍按(8.4.12a)式進行變換
y2=1./x1/2-sqrt((1./x1/2).^2+1);
plot(x1,[y1;y2],'xr')           % 加深畫出帶阻通帶頻域的映射曲線
plot(x1,-3*ones(size(x1)),'+k') % 加深低通原型通帶范圍的點型和顏色
plot(3*ones(size(y1)),y1,'og')  % 加深帶阻濾波器通帶范圍的點型和顏色
plot(3*ones(size(y2)),y2,'og')
grid on
axis([-3,3,-3,3])