fgp842.m

来自「matlab数字信号处理代码」· M 代码 · 共 18 行

%《数字信号处理教程——MATLAB释义与实现》第八章图8.4.2生成程序fgp842
% 模拟滤波器频带变换（低通变高通）
% 电子工业出版社出版  陈怀琛编著 2004年8月
%
x0=logspace(-2,1,100);          % 低通原型的频率设定
x=[-fliplr(x0),x0];             % 变换为正负对称频率范围
y=1./x;                         % 对频率做倒数变换
plot(x,y),hold on           % 画出全频域的映射曲线
set(gcf,'color','w')
grid on
axis([-3,3,-3,3])
x1=x(find(x<=1&x>=-1));         % 低通原型的通带范围（设带宽为1）
y1=1./x1;                       % 对通带范围频率做倒数变换
plot(x1,y1,'xr')                % 画出通带频域的映射曲线
plot(x1,-3*ones(size(x1)),'+k') % 加深低通原型通带范围的点型和颜色
plot(3*ones(size(y1)),y1,'og')  % 加深高通滤波器通带范围的点型和颜色