ex_1.m

来自「信号与系统分析以matlab实现」· M 代码 · 共 30 行

%ex_1
clear all;
a=[8 2 3 1 5];
b=[1 3 2];
%
p=roots(a)              %求系统极点
z=roots(b)             %求系统零点
p=p';			    	%将极点列向量转置为行向量
z=z';						%将零点列向量转置为行向量
x=max(abs([p z]));			%确定纵坐标范围
x=x+0.1;
y=x;						%确定横坐标范围
figure(1)
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(z),imag(z),'o') 	%画零点
title('连续系统零极点图')	%标注标题
text(0.2,x-0.2,'虚轴')
text(y-0.2,0.2,'实轴')
%
[z,p]=tf2zp(b,a)
figure(2)
zplane(z,p),title('连续系统零极点图')