fgp835a.m

digital signal processing-MATLAB

%《数字信号处理教程——MATLAB释义与实现》第八章图8.3.5生成程序fgp835a
% 双线性变换中s 平面到z平面的映射（T任设）
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
a=[-2:0.1:0];                       % 设定实部分度
b=[-5:0.1:5];                       % 设定虚部分度
[X,Y]=meshgrid(a,b);                % 建立平面两维自变量网格
S=X+j*Y;                            % 建立平面复数自变量
subplot(2,2,1),plot(X,Y)            % 画等实部线
T=input('T= ');
Z=(1+T*S/2)./(1-T*S/2);             % 进行S到Z的复数映射
C=real(Z);
D=imag(Z);
subplot(2,2,2),plot(C,D),grid on    % 画等实部线所映射的z平面上曲线 
subplot(2,2,3),plot(X',Y')                 % 画等虚部线
subplot(2,2,4),plot(C',D'),grid on   % 画等虚部线所映射的z平面上曲线
set(gcf,'color','w')