信号与线性系统实验指导书(4)

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,'实轴') p =

-0.6155 + 0.6674i -0.6155 - 0.6674i 0.4905 + 0.7196i 0.4905 - 0.7196i z = -2 -1

例2：某离散系统的系统函数为：

z?1H(z)?5 43z?z?1求出系统函数的零极点，画出零极点分布图。 a=[3 -1 0 0 0 1]; b=[1 1];

p=roots(a); %求系统极点 z=roots(b); %求系统零点 p=p'; %将极点列向量转置为行向量 z=z'; %将零点列向量转置为行向量 x=max(abs([p z 1])); %确定纵坐标范围 x=x+0.1;

y=x; %确定横坐标范围 figure(1),clf,hold on axis([-x x -y y]) %确定坐标轴显示范围 w=0:pi/300:2*pi; t=exp(i*w); plot(t) %画单位园