fgp844.m

digital signal processing-MATLAB

%《数字信号处理教程——MATLAB释义与实现》第八章图8.4.4生成程序fgp844
% 模拟滤波器频带变换（低通变带通）
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
x0=logspace(-2,1,1000);         % 低通原型的频率设定
x=[-fliplr(x0),0,x0];           % 变换为正负对称频率范围
y1=x/2+sqrt((x/2).^2+1);        % 对频率按(8.4.7a)式进行变换
y2=x/2-sqrt((x/2).^2+1);
plot(x,[y1;y2]),hold on         % 画出全频域的映射曲线
set(gcf,'color','w')            % 置图形背景色为白
x1=x(find(x<=1&x>=-1));         % 低通原型的通带范围（设带宽为1）
y1=x1/2+sqrt((x1/2).^2+1);      % 对通带频率范围按(8.4.7a)式进行变换
y2=x1/2-sqrt((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])