Laboratorio de Control Analógico I Práctica No. 4.

- Modelos Matemáticos de Sistemas Lineales

Práctica 4
Modelos Matemáticos de Sistemas Lineales

Objetivo. El objetivo de esta práctica es dar una introducción al uso de algunas funciones de
Scilab orientadas al manejo de modelos matemáticos de sistemas lineales invariantes en el
tiempo, en función de transferencia y en variables de estado, tales como syslin, poly, pfss,
trfmod, tf2ss, ss2tf y abcd.

Introducción.
El toolbox CACSD (Computer Aided Control Systems Design) de Scilab proporciona un enorme
conjunto de funciones que van desde las funciones básicas para introducir modelos de sistemas
lineales invariantes en el tiempo (SLIT's) hasta herramientas poderosas para el análisis y diseño de
sistemas de control tanto de una entrada - una salida (SISO) o de múltiples entradas y múltiples
salidas (MIMO).

El modelo en Función de Transferencia.

Un sistema lineal invariante en el tiempo (SLIT por sus siglas en español y LTI por sus siglas en
inglés) es aquel que cumple con el principio de superposición. Si el sistema tiene una entrada u (t )
y una salida y (t ) , su representación en el dominio del tiempo corresponde a una ecuación
diferencial lineal con coeficientes constantes como la siguiente:

(n) ( n −1) (m) ( m −1)

an y + an −1 y + ... + a1 y
+ a0 y = bm u + bm −1 u + ... + b1u
+ b0u (4.1)

(k )
donde los coeficientes a0 , a1 ,..., an , b0 , bm ,..., bm son constantes, y donde la notación y representa
dk y
la k-ésima derivada , para k=0,1,2,....,n.
dt k

Aplicando Transformada de Laplace a ambos miembros de la ecuación (4.1), y suponiendo

Y ( s)
condiciones iniciales cero, podemos despejar el cociente , para obtener la denominada Función
U ( s)
de Transferencia del sistema:
Y ( s ) bm s m + bm −1s m −1 + ... + b1s + b0
G (s) = = (4.2)
U ( s ) an s n + an −1s n −1 + ... + a1s + a0

Obsérvese que la función de transferencia de un sistema del tipo (4.1) siempre es una función
racional, es decir, una división de polinomios en potencias de la variable de Laplace (s).

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 1

Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales
Una de las virtudes de la representación en función de transferencia (4.2), es que permite escribir
la salida Y ( s ) del sistema como una simple multiplicación de la entrada U ( s ) por G ( s ) , es decir,

Y ( s ) = G ( s)U ( s ) (4.3)

En otras palabras, la función de transferencia G ( s ) viene siendo la ganancia aplicada a la entrada

U ( s) para obtener la salida Y ( s) .

El razonamiento anterior sustentado en la expresión (4.3) permite obtener la regla más simple de la
representación en diagrama de bloques del sistema dado por la ecuación diferencial (4.1) de acuerdo a la
figura siguiente

U(s) G(s) Y(s)=G(s)U(s)

Figura 4.1.- Diagrama de Bloques del sistema dado por (4.1)

Ejemplo: El sistema dado por la ecuación diferencial 2

y + 6 y
+ 4 y = u
+ 3u tiene la representación
en función de transferencia siguiente
s+3
G (s) = 2
(4.4)
2s + 6s + 4

La representación (4.4) no es única, ya que numerador y denominador son polinomios en

potencias de s y se pueden escribir de muchas maneras. Algunas de ellas son las siguientes:

Forma factorizada

Todo polinomio se puede factorizar en binomios de primer orden. Esta forma se prefiere cuando se quieren
evidenciar los polos y los ceros de la función de transferencia, así, por ejemplo, la función de transferencia
(4.4) se puede escribir en su forma factorizada como sigue
s+3
G( s) =
(2 s + 4)( s + 1)

y procurando que todos los factores de primer orden tengan coeficiente 1 en la variable s , se obtiene la
forma estándar:
0.5( s + 3)
G( s) = (4.5)
( s + 2)( s + 1)
De la forma (4.5) resultan evidentes:
los polos: -2 y -1,
los ceros: -3
la ganancia del sistema: 0.5

Observación: Aunque todo polinomio se puede factorizar como producto de binomios de

primer orden, cuando el polinomio tiene raíces complejas esta factorización conduce al uso
de números complejos, por esta razón se prefiere dejar sin factorizar los factores cuadráticos
que tienen raíces complejas conjugadas.

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 2

Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales

Ejemplo 4.1. El polinomio p ( s) = s 3 + 3s 2 + 4s + 4 se puede factorizar como el producto de tres

binomios de primer orden, un binomio por cada una de sus raíces, reales o complejas, es decir,
7 7
p( s ) = ( s + 2)( s + 0.5 + 2 i )( s + 0.5 − 2 i) , pero para evitar el uso de números complejos se preferirá
factorizar como p( s) = ( s + 2)( s 2 + s + 2) . Es decir, los dos binomios con raíces complejas se convierten
en el trinomio cuadrático ( s 2 + s + 2) . Mediante comandos de Scilab esto se puede verificar como sigue

p=poly([ 4 4 3 1],'s','c')//define el polinomio mediante sus coeficientes

p =
2 3
4 + 4s + 3s + s
roots(p) //calcula sus raíces
ans =
- 2.
- 0.5 + 1.3228757i
- 0.5 - 1.3228757i
polfact(p) //obtiene los factores irreducibles sin usar números complejos
ans =
2
1 2 + s 2 + s + s

Forma de fracciones parciales

Usando la técnica de expansión en fracciones parciales, la forma (4.4) también puede expandirse
en una suma de fracciones, cada una con denominador de primer orden (llamadas fracciones
simples o parciales) una por cada polo. Para el ejemplo anterior, la forma de fracciones parciales
queda
1 0.5
G (s) = − (4.6)
s +1 s + 2

La forma de fracciones parciales es de interés en el problema de diseño, ya que permite una

implementación de la función de transferencia como la conexión en paralelo de subsistemas
sencillos de primer orden.

El comando syslin para definir funciones de transferencia (F.T.)

Es el comando básico para introducir modelos de sistemas lineales tanto en función de

transferencia (F.T.) como en variables de estado, para sistemas SISO (una entrada - una salida) y
para sistemas MIMO (múltiples entradas - múltiples salidas), es el comando syslin. Primero se
describirá su uso para la introducción de modelos en función de transferencia.

La sintaxis del comando syslin es G=syslin(dominio,num,den). Este comando define la

función de transferencia G con el numerador especificado por el polinomio num y el denominador
especificados por el polinomio den. El dominio puede ser 'c' para un sistema continuo (analógico)
o 'd' para un sistema discret (digital).

Ejemplo 4.2. Se puede introducir la función de transferencia dada por (4.4) de tres maneras:

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 3

Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales
1) Definiendo los polinomios del numerador y denominador por sus coeficientes
num=poly([3 1],'s','c');//Define numerador a partir de sus coeficientes
den=poly([4 6 2],'s','c');//Define denomin. a partir de sus coeficientes
G=syslin('c',num,den) //Define la F.T. continua
G =
3 + s
-----------
2
4 + 6s + 2s

2) Definiendo primero la variable de Laplace y luego escribiendo la expresión algebraica (4.4)

s=poly(0,'s'); //Define la variable de Laplace
G=(s+3)/(2*s^2+6*s+4) //Define la F.T.
G =
3 + s
-----------
2
4 + 6s + 2s

3) Utilizando la forma factorizada dada por (4.5) se procede como sigue:

K=0.5; //Ganancia del sistema
ceros=poly([-3],'s'); //Numerador definido por su vector de raíces
polos=poly([-2 -1],'s');//Denominador definido por su vector de raíces
G=syslin('c',K*ceros,polos) //Define la F.T.
G =
1.5 + 0.5s
----------
2
2 + 3s + s

El comando trfmod y la forma factorizada de una F.T.

El comando trfmod permite visualizar en una ventana los factores irreducibles del numerador y
denominador de una función de transferencia de un sistema SISO.

Ejemplo 4.3: Para visualizar la misma F.T. dada del ejemplo anterior una vez definida por
cualquiera de los tres métodos descritos, se ejecuta el comando trfmod como sigue y Scilab
despliega la ventana mostrada en la figura 4.2, en la cual se pueden ver claramente: la ganancia y
los factores del numerador y del denominador, los cuales coinciden con la forma factorizada (4.5)

Figura 4.2. Ventana desplegada por el comando trfmod(G,'p')

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 4

Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales

trfmod(G,'p'); //el parámetro 'p' es opcional

El comando pfss y la forma de fracciones parciales.

El comando pfss(G) realiza el cálculo de la expansión en fracciones parciales de la función de

transferencia especificada G.

Ejemplo 4.4. Para obtener la representación de la F. T. (4.4) en su forma de fracciones parciales se

puede proceder como sigue:
num=poly([1 3],'s','c');// Define numerador
den=poly([2 6 4],'s','c');// define denominador
G=syslin('c',num,den); //Define F.T.
frac=pfss(G) //Calcula vector de fracciones parciales
frac =
frac(1)
1
-----
1 + s
frac(2)
- 0.5
-----
2 + s

Es decir, la representación de (4.4) en fracciones parciales es la esperada (4.6). Obsérvese que la

suma de estas fracciones parciales produce la F.T. original como es de esperarse:
frac(1)+frac(2) //Recalcula la función original.
ans =
1.5 + 0.5s
----------
2
2 + 3s + s

Los comandos presentados arriba se pueden complementar con los comandos vistos para manejar
polinomios explicados en la práctica No. 2. A continuación se da un ejemplo.

Ejemplo 4.5: En ocasiones la F.T. no está escrita en ninguna de las tres formas mencionadas arriba,
s2 + 2
sino en una combinación de ellas, por ejemplo: G ( s ) = . En este caso se
s ( s 2 + 3s + 2)( s + 0.5)
puede introducir la F.T. de diferentes maneras:

1) Usando multiplicación de polinomios mediante el comando conv.

num=poly([2 0 1],'s','c'); //Numerador
d1=[0 1]; //coeficientes del factor s
d2=[2 3 1]; //coeficientes del factor s^2+3s+2
d3=[0.5 1]; //coeficientes del factor s+0.5
den=conv(d1,d2);
den=conv(den,d3); //Obtiene coeficientes del denominador
den=poly(den,'s','c'); //convierte coeficientes a polinomio
G=syslin('c',num,den) //define F.T.
G =
2

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 5

Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales
2 + s
-------------------
2 3 4
s + 3.5s + 3.5s + s

2) Factorizando en tres F.T. y luego multiplicándolas.

num1=poly([2 0 1],'s','c');
den1= poly([2 3 1],'s','c');
G1=syslin('c',num1,den1)
G1 =
2
2 + s
---------
2
2 + 3s + s
num2=poly([1],'s','c');
den2= poly([0 1],'s','c');
G2=syslin('c',num2,den2)
G2 =
1
-
s
num3=poly([1],'s','c');
den3= poly([0.5 1],'s','c');
G3=syslin('c',num3,den3)
G3 =

1
-------
0.5 + s
G=G1*G2*G3 //Multiplica los tres factores
G =
2
2 + s
-------------------
2 3 4
s + 3.5s + 3.5s + s

3) Definiendo la variable de Laplace 's' y escribiendo la expresión algebraica de la F.T.

s=poly([0],'s'); //los corchetes del cero son opcionales
G=(s^2+2)/s/(s^2+3*s+2)/(s+0.5)
G =
2
2 + s
-------------------
2 3 4
s + 3.5s + 3.5s + s

Ejemplo de obtención de un modelo en función de transferencia.

Consideremos el circuito RLC serie de la figura 4.3a Para obtener un modelo en función de
transferencia lo primero que se debe saber es cuales son las variables de entrada y de salida. En
este caso se desea obtener la función de transferencia que relaciona el voltaje vc en el capacitor
Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 6
Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales
(salida) al voltaje Vi aplicado por la fuente (entrada).

La manera más sencilla de obtener el modelo en función de transferencia de un sistema es

expresar sus componentes directamente en el dominio de Laplace, es decir, en forma de
impedancias complejas, como se muestra en la figura 4.3b

L ZL(s) ZR
R

C + +
+ i(t) + I(s)
Vi(t) vC(t) Vi(s) ZC(s) VC(s)
- -
- -

(a) (b)
Figura 4.3 a) Circuito RLC serie. b) Circuito RLC serie en términos de impedancias complejas.

El modelo de impedancias complejas permite aplicar la ley de Ohm directamente sobre cada
impedancia para obtener la caída de voltaje entre las terminales de cada componente
directamente en el dominio de Lapalace de acuerdo a la siguiente Tabla

Impedancia
Componente Relación v(t), i(t) Relación V(s), I(s)
Z=V(s)/I(s)
Resistencia (R) v(t ) = R i (t ) V (s) = R I (s) ZR = R
di (t )
Inductancia (L) v(t ) = L V ( s) = Ls I ( s ) Z L = Ls
dt
1 t 1 1
C ∫0
Capacitancia (C) v(t ) = i(t )dt V ( s) = I (s) ZC =
Cs Cs

Por lo tanto, aplicando la regla del divisor de tensión al circuito de la figura 4.3b obtenemos
ZC
Vc ( s ) = Vi ( s ) (4.7)
Z L + Z R + ZC
y sustituyendo las expresiones de la tabla anterior se obtiene
1
Vc ( s) = Cs V ( s)
1 i
sL + R +
Cs
Simplificando
1
Vc ( s ) = Vi ( s )
LCs 2 + RcS + 1

Por lo tanto, la función de transferencia buscada es

Vc ( s ) 1
G(s) = = 2
(4.8)
Vi ( s ) LCs + RcS + 1

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 7

Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales

Modelo en Variables de Estado (o de Espacio de Estado).

Consideremos un SLIT SISO cuya entrada es u (t ) y salida y (t ) , el cual se puede representar por la
ecuación diferencial lineal de orden n dada por (4.1), la cual se reescribe a continuación
(n) ( n −1) (m) ( m −1)
an y + an −1 y + ... + a1 y
+ a0 y = bm u + bm −1 u + ... + b1u
+ b0u (4.9)

Aunque no siempre es sencillo, se puede elegir un conjunto de variables de estado (que en

general pueden ser sumas ponderadas de la entrada, la salida y/o algunas de sus derivadas) para
transformar (4.9) en n ecuaciones diferenciales lineales de primer orden llamadas ecuaciones de
estado, de la forma
x
1 = a11 x1 + a12 x2 + ... + a1n xn + b1u
x
2 = a21 x1 + a22 x2 + ... + a2 n xn + b2 u
(4.10)
...
x
n = an1 x1 + an 2 x2 + ... + ann xn + bn u

Además, la salida del sistema normalmente es una combinación de estados y de la entrada que se
puede escribir como la ecuación de salida siguiente

y = c1 x1 + c2 x2 + ....cn xn + du (4.11)

Es decir, en forma matricial compacta

x
= Ax + bu (4.12)
y = cx + du (4.13)

Donde x es el vector de estados, A es una matriz cuadrada n × n de coeficientes constantes, b es

un vector columna de coeficientes constantes, c es un vector renglón de coeficientes constantes y
d es un escalar constante (en la mayoría de los casos d = 0 ), es decir,

 x1   x
1   a11 a12 ...a1n   b1 

x   x
 a a22 
... a2 n  b 
x =  2  , x
=  2  , A =  21 , b =  2  , c = [ c1 c2 ... cn ] (4.14)
        
       
 xn   x
n   an1 an 2 ... ann  bn 

Las ecuaciones de estado (4.12) y la ecuación de salida (4.13) se pueden transformar de manera
directa a su Función de Transferencia equivalente aplicando Transformada de Laplace suponiendo
condiciones iniciales cero, obteniéndose
sX ( s ) = AX ( s ) + bU ( s )
(4.15)
Y ( s ) = cX ( s ) + dU ( s )

despejando el vector de estado

X ( s ) = ( sI − A) −1 bU ( s ) (4.16)

sustituyendo en la ecuación de salida

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 8
Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales

Y ( s ) = c( sI − A) −1 b + d  U ( s) (4.17)

Es decir, la función de transferencia correspondiente al modelo en variables de estado (4.12),

(4.13) es
Y (s)
G (s) = = c( sI − A) −1 b + d (4.18)
U (s)

Observación: Dado un modelo en variables de estado, le corresponde una única función de

transferencia dada por (4.18), sin embargo, dada una función de transferencia, ésta puede
corresponder a una infinidad de modelos en variables de estado válidos, dependiendo del
conjunto de variables de estado elegido. Esta situación se ilustra en la figura 4.3

Modelo de estado 1 G ( s) = c1[ sI − A1 ]−1 b1 + d1

x
1 = A1 x1 + b1u
y = c1 x1 + d1u

Modelo de estado 2 G ( s) = c2 [ sI − A2 ]−1 b2 + d 2 Modelo en Función

x
2 = A2 x2 + b2u de Transferencia
y = c2 x2 + d 2u Y ( s ) = G ( s)U ( s )

 
Modelo de estado N
x
N = AN xN + bN u
y = cN xN + bN u G ( s ) = cN [ sI − AN ]−1 bN + d N

Figura 4.3. A una sola F. T. le corresponde una infinidad de Modelos en Variables de Estado

Comandos de Scilab para manejar modelos en Espacio de Estados.

Scilab provee los comandos explicados en la siguiente tabla para manejar modelos en el espacio
de estado de SLIT's continuos.

Comando explicación
Crea un sistema tipo SS (Espacio de Estado), es decir, sist contiene la
información del modelo de espacio de estado:
sist=syslin('c',A,B,C,D,x0) dx/dt = Ax(t) + Bu(t)
y(t) = Cx(t) + Du(t)
con las condiciones iniciales x(0)=x0 (opcionales)
Obtiene las matrices A,B,C,D correspondientes al modelo de espacio de
[A,B,C,D]=abcd(sist)
estado de sist.
Convierte el modelo de espacio de estado dado por sist a la matriz de
G=ss2tf(sist) transferencia correspondiente G. Es decir, evalúa la ecuación (4.18) para el
caso MIMO.
Convierte la matriz de funciones de transferencia G a una de sus
sist=tf2ss(G)
representaciones en espacio de estado.

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 9

Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales

Ejemplo 4.6. Conversión entre F.T. y Espacio de Estado.

num=poly([1 1],'s','c');
den=poly([6 5 1],'s','c');
G=syslin('c',num,den) // Define la F.T. G=(s+1)/(s^2+5*s+6)
G =
1 + s
---------
2
6 + 5s + s

sist=tf2ss(G)// Obtiene modelo de espacio de estado correspondiente a G

sist =
sist(1) (state-space system:)
!lss A B C D X0 dt !
sist(2) = A matrix =
- 4.4 - 0.8
4.2 - 0.6
sist(3) = B matrix =
- 1.2649111
0.6324555
sist(4) = C matrix =
- 0.7905694 5.551D-17
sist(5) = D matrix =
0.
sist(6) = X0 (initial state) =
0.
0.
sist(7) = Time domain =
c

[A,B,C,D]=abcd(sist);//Obtiene las Matrices del modelo de estado

s=poly(0,'s');
G1=C*inv(s*eye(2,2)-A)*B+D // Evalúa la fórmula (4.18)
G1 =
1 + s
---------
2
6 + 5s + s

Obsérvese que la F.T. obtenida evaluando la fórmula (4.18) es corresponde perfectamente a la

función de transferencia original. Si se obtiene mediante el comando ss2tf el resultado es el
mismo:

G2=ss2tf(sist)
G2 =
1 + s
---------
2
6 + 5s + s

Ejemplo de obtención de un modelo en variables de estado.

Consideremos el circuito RLC serie de la figura 4.3a, para el cual ya se obtuvo la función de
Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 10
Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales
transferencia (4.8). Se desea ahora obtener un modelo en variables de estado.

Como ya se mencionó, hay una infinidad de modelos posibles, dependiendo de las variables de
estado elegidas. Por esta razón, si no se especifica otra cosa es tradicional elegir como variables
de estado del circuito: voltajes en capacitores y corrientes en inductores, es decir, en este caso las
variables de estado serán:

 x1   vc (t ) 
 =  (4.19)
 x2   i(t ) 

Aplicando la Ley de Voltajes de Kirchoff a la unica malla del circuito, se obtiene

vL (t ) + vR (t ) + vC (t ) = Vi (t ) (4.20)
es decir,
di (t )
L + Ri(t ) + vC (t ) = Vi (t ) (4.21)
dt
Además, se sabe que
1 t
vC (t ) = ∫ i (t )dt (4.22)
C 0
por lo tanto,
dvC (t ) 1
= i (t ) (4.23)
dt C
despejando de (4.21)
di (t ) 1 R 1
= − vC (t ) − i(t ) + Vi (t ) (4.24)
dt L L L

Expresando (4.23) y (4.24) en términos de las variables de estado definidas por (4.19) se obtiene
1
x
1 = x2 (4.25)
C
1 R 1
x
2 = − x1 − x2 + Vi (t ) (4.26)
L L L

Las dos anteriores son las ecuaciones de estado del circuito, sin embargo, es costumbre
expresarlas en forma matricial como sigue

 x
1   0   x
1   0 
1
C
 x
 =  − 1  + V (t ) (4.27)
 2  L −   x
2   L1  i
R
L

y la ecuación de salida es
y = x1 (4.28)

Es decir, las matrices A, b, c y d correspondientes a (4.12), (4.13) son las siguientes

0 1
C 0
A= 1  , b =  1  , c = [1 0] , d = 0
R
(4.29)
− L − L L

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 11

Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales
Caso MIMO:

Cuando un sistema tiene varias entradas y varias salidas, la noción de Función de Transferencia se
puede generalizar, a la de una matriz que contiene como elementos funciones de transferencia
escalares y que se denomina Matriz de Transferencia. En la figura 4.4 se muestra en forma de
bloque un sistema MIMO con p entradas y q salidas.

U1(s) Y1(s)
U2(s) Y2(s)
... G(s) ...
Up(s) Yq(s)
Y(s) = G(s) U(s)
Figura 4.4.- Diagrama de Bloques de un sistema con p entradas y q salidas

Para el sistema MIMO de la figura 4.4 sigue siendo válida la relación

Y ( s ) = G ( s)U ( s ) (4.30)

sin embargo, (4.30) es ahora una expresión matricial donde:

 Y1 ( s )   U1 ( s )   G11 ( s ) G12 ( s ) ... G1 p ( s ) 

Y ( s )  U ( s )  G ( s ) G ( s ) ... G2 p ( s ) 
Y (s) =  2  , U (s) =  2  , G( s) = 
21 22
(4.31)
 ...   ...       
     
Yq ( s )  U p ( s )  Gq1 ( s ) Gq 2 ( s ) ... Gqp ( s ) 

Donde G ( s ) es la matriz de transferencia cuyos elementos se pueden obtener por superposición

como sigue:
Yi ( s )
Gij ( s ) = para k ≠ j , i = 1, 2,..., q, j = 1, 2,..., p (4.32)
U j ( s)
U k ( s ) =0

Tanto el comando syslin como los comandos ss2tf y tf2ss trabajan con sistemas MIMO y
por lo tanto manejan matrices de transferencia de manera directa.

 1 s 
 s2 + s + 1 s +1 
Ejemplo: Para introducir en Scilab la Matriz de Transferencia: G ( s ) =   se
 s+2 s 
 s + 1 s + s + 1 
2

puede proceder como sigue:

//define numeradores:
num11=1; num12=poly([0 1],'s','c');
num21=poly([2 1],'s','c'); num22=poly([0 1],'s','c');
//Define denominadores:
den11=poly([1 1 1],'s','c'); den12=poly([1 1],'s','c');
den21=poly([1 1],'s','c'); den22=poly([1 1 1],'s','c');
Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 12
Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales
//Define Matriz de transferencia
G=syslin('c',[num11,num12;num21,num22],[den11,den12;den21,den22])

O también:
//define numeradores:
num11=1; num12=poly([0 1],'s','c');
num21=poly([2 1],'s','c'); num22=poly([0 1],'s','c');
//Define denominadores:
den11=poly([1 1 1],'s','c'); den12=poly([1 1],'s','c');
den21=poly([1 1],'s','c'); den22=poly([1 1 1],'s','c');
//Define componenentes de la matriz de transferencia
G11= syslin('c',num11,den11); G12= syslin('c',num12,den12);
G21= syslin('c',num21,den21); G22=syslin('c',num22,den22);
//Define la Matriz a partir de sus componentes
G=[G11,G12;G21,G22]

O también, usando la variable de Laplace:

s=poly(0,'s');
G=[1/(s^2+s+1), s/(s+1);(s+2)/(s+1), s/(s^2+s+1)]

Una vez que se ha introducido correctamente la matriz de tranferencia G(s) por cualquiera de los
tres métodos anteriores, Scilab la desplegará como sigue:
G =
1 s
--------- -----
2
1 + s + s 1 + s

2 + s s
----- ---------
2
1 + s 1 + s + s

Ejercicio: Escribir una función de Scilab en un archivo-M que a partir de los coeficientes de los
polinomios del numerador (num) y denominador (den) de una F. T. obtenga la siguiente
información:
• La función de transferencia (G)
• Los polos (p) y los ceros (z).
• Las matrices A,B,C,D del modelo en espacio de estado.

El comando para ejecutar la función deberá ser: [G,p,z,A,B,C,D]=fun_trans(num,den).

Desarrollo de la Práctica.

1. Probar todos los ejemplos propuestos por el profesor conforme los va explicando.
2. Realizar todos los ejercicios propuestos.
3. Contestar el cuestionario de evaluación de la práctica.

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 13

Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 4.- Modelos Matemáticos de Sistemas Lineales
Reportar:

1) El código escrito y comentado para la función fun_trans del ejercicio propuesto acompañado
de dos ejemplos de ejecución.
VC ( s )
2) Obtener la F.T. correspondiente al siguiente circuito y escribir una función que la defina
Vi ( s )
en términos de los valores de R, L y C en Scilab, de la forma G=FuncionTrans(R,L,C).

iL(t)
+
+ R
Vi(t) C VC(t)
-
-

Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 14

Elaboró: José Juan Rincón Pasaye