% 30.10.00 HA V3-kurssi.
% Matlab-skripti, jolla piirretään siirtofunktion nollakohdat (o) ja navat (x).
% Siirtofunktio (rationaalifunktio) esitetään vektoriparina P, Q.
%   P:ssä on osoittajan kertoimet, korkeimmasta alimpaan asteeseen
%   Q:ssa nimittäjän kertoimet       --------------- " ------------
%
% Modifioi tätä esimerkkiä:

Q=[9 12 13]; P=[2 3];
z=roots(P)  % zeros, nollakohdat
p=roots(Q)  % poles, navat
a=min([real(z(:)'),real(p(:)')]);b=max([real(z(:)'),real(p(:)')])
c=min([imag(z(:)'),imag(p(:)')]);d=max([imag(z(:)'),imag(p(:)')])

clf
plot(z+0.00001*i,'or'); hold on; plot(p+0.000001*i,'xb');grid;
axis([a-0.2, b+0.2, c-0.2, d+0.7]);legend('Nollakohdat','Navat')
shg
