Problema Repaso (Laboratorio)

Para esta semana en el laboratorio de automatización a todos nos toco resolver un problema mayormente relacionados con transformadas de Laplace.

Mi problema

EL problema se trata de obtener una expresión numerador/ denominador donde tiene:

• Ceros en s=-1, s = -2. 
• Polos en s = 0, s = 4, s = -6
• Ganancia k = 5 

 Esto hacerlo utilizando Matlab, pero como no tengo acceso a ese software yo utilizare Octave.

La respuesta la obtuve basandome en un código para MatLab que viene en el libro citado y el código es el siguiente:

function null =  lab2(null)
  z = [-1, -2];
  p = [0, 4, -6];
  k = 5;
  [num, den] = zp2tf(z,p,k);
  num
  den
endfunction

Donde la función zp2tf se encarga de generar la expresión utilizando los datos dados de parámetro la cuál obtuve de la siguiente página y es la siguiente:

http://users.powernet.co.uk/kienzle/octave/matcompat/scripts/FIXES/zp2tf.m

function [b, a] = zp2tf(z, p, g)
  if nargin != 3 || nargout != 2
    usage("[b, a] = zp2tf(z, p, g)");
  endif
  try cplxpair(z); catch error("zp2tf: could not pair complex zeros"); end
  try cplxpair(p); catch error("zp2tf: could not pair complex poles"); end
  b = g*real(poly(z));
  a = real(poly(p));
endfunction

El programa me arrojo como resultado lo siguiente:

octave:1> lab2.m
num =

    5   15   10

den =

    1    2  -24    0

Lo cual se interpreta como:

El problema fue sacado del libro:
Ogata, Katsuhiko. Modern Control Engineering / Katsuiko Ogata.
Tehran: A'eizh, 1381 = 2002.
964.: ill
Reprint of 4th ed. 2002, Prentice Hall, New Jersey.


Referencias:
http://users.powernet.co.uk/kienzle/octave/matcompat/scripts/FIXES/zp2tf.m