zero_pole.m

来自「经典《信号与系统》教程的matlab例程,对深入理解信号与系统相关概念有很大帮助」· M 代码 · 共 32 行

%===========================================================
%  程序名：zero_pole 
%  本程序绘制系统函数的零点、极点在s平面上的分布图。
%  A:系统函数分母多项式系数向量，实验者从键盘上输入
%  B:系统函数分子多项式系数向量，实验者从键盘上输入
%  p:函数返回的系统函数极点位置行向量
%  q:函数返回的系统函数零点位置行向量
%===========================================================
b=input('请输入系统函数分子多项式系数：');
a=input('请输入系统函数分母多项式系数：');
p=roots(a);                       %求系统极点
q=roots(b);                       %求系统零点
p=p';					          %将极点列向量转置为行向量
q=q';					          %将零点列向量转置为行向量
x=max(abs([p q]));				  %确定纵坐标范围
x=x+1;
y=x;						      %确定横坐标范围
clf
%subplot(2,1,1);
axis([-x x -y y]);		          %确定坐标轴显示范围
axis('square')

%===========================================================
plot([-x x],[0 0]);hold on; 	       %画横坐标轴
plot([0 0],[-y y]);hold on;		       %画纵坐标轴
plot(real(p),imag(p),'x');hold on;     %画极点
plot(real(q),imag(q),'o');hold on;     %画零点
grid on;
title('系统函数零极点分布');
ylabel('虚轴');
xlabel('实轴');