sjdt.m

来自「1． 该光盘包含了“信号与系统分析及MATLAB实现”一书的第十一章所有程序； 」· M 代码 · 共 25 行

function [p,q]=sjdt(A,B)
%绘制连续系统零极点图程序 
%A:系统函数分母多项式系数向量
%B:系统函数分子多项式系数向量
%p:函数返回的系统函数极点位置行向量
%q:函数返回的系统函数零点位置行向量
p=roots(A);                                   　%求系统极点
q=roots(B);;                                   	　%求系统零点
p=p';													%将极点列向量转置为行向量
q=q';													％将零点列向量转置为行向量
x=max(abs([p q]));								%确定纵坐标范围
x=x+0.1;
y=x;													%确定横坐标范围
clf
hold on
axis([-x x -y y]);								%确定坐标轴显示范围
axis('square')
plot([-x x],[0 0]) 								%画横坐标轴
plot([0 0],[-y y])								 %画纵坐标轴
plot(real(p),imag(p),'x')						 %画极点
plot(real(q),imag(q),'o') 					%画零点
title('连续系统零极点图')						 %标注标题
text(0.2,x-0.2,'虚轴')
text(y-0.2,0.2,'实轴')