fig2_44.m

来自「现代控制系统（matlab程序）」· M 代码 · 共 22 行

% Chapter 2: Figure 2.44, p. 85
%    Transfer function example for G(s) and H(s).
%
numg=[6 0 1]; deng=[1 3 3 1];
z=roots(numg)
%
p=roots(deng)
%
n1=[1 1]; n2=[1 2]; d1=[1 2*i]; d2=[1 -2*i]; d3=[1 3];
numh=conv(n1,n2);
denh=conv(d1,conv(d2,d3));
printsys(numh,denh)
%
num=conv(numg,denh);
den=conv(deng,numh);
printsys(num,den)
%
% Generate pole-zero map in Figure 2.43, p. 84
%
pzmap(num,den)
title('Pole-Zero Map')