%Ejemplo de código de Matlab/Octave para obtención de diagrama de polos y ceros, cálculo de expansión en fracciones parciales y visualización de la respuesta en frecuencia.

% DIAGRAMA DE POLOS Y CEROS - RESIDUOS
% Realiza la gráfica polos-ceros de la transformada racional en que el
% denominador posee los coeficientes de las potencias en orden decreciente
% dadas en el vector num, lo mismo para el denominador en den. Luego se
% calculan residuos y polos para la expansión en fracciones parciales.

% num=[2 5  12];
% den=[1 4 14 20];
% Hsys=tf(num,den);
% [p,z]=pzmap(Hsys)
% 
% [rs ps ks]=residue(num,den)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Visualización de la respuesta en frecuencia.
% omega=[-4:0.05:4];
% H=freqresp(Hsys,omega);
% Hmag=abs(squeeze(H));
% figure;
% plot(omega,Hmag);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Valores de los zeros, polos y ganancia
%[zeros polos gain]=tf2zp(num,den)