GUIAS DE SIMULINK:

SISTEMAS MULTIVARIABLES

Cerpa Kelly, Hurtado Maylen

UNIVERSIDAD DEL ATLANTICO.

4 De Junio de 2012.

OBJETIVO

 Ilustrar el manejo de Simulink, mediante la aplicación del método de simulación de

espacio de los estados, para obtener de forma grafica el comportamiento de un
proceso según el tipo de cambio aplicado.

EJEMPLO 8.5.1: “Dinámica de una columna de absorción”.

En una columna de absorción la corriente de gas entrante a la columna circula en

contracorriente con el líquido. El gas asciende como consecuencia de la diferencia de presión
entre la entrada y la salida de la columna. El contacto entre las dos fases produce la
transferencia del soluto de la fase gaseosa a la fase líquida, debido a que el soluto presenta
una mayor afinidad por el disolvente. Las columnas de absorción a menudo contienen platos
con una capa líquida fluyendo a través de éstos; estos platos son a menudo modelados como
etapas de equilibrio.

V L
y1 xf

1
2
3

N-2
N-1
N

V L
yN+1 xN

Figura 1. Columna de absorción de N etapas.

Consideremos una columna de absorción de 5 etapas, en la cual el flujo de líquido de alimento
es L=4/3 Kg mol de aceite inerte/min, el flujo de vapores V=5/3 Kg mol de aire/min. La
composición del líquido de alimento es xf=0.0 Kg mol de benceno/Kg mol de aceite inerte y la
composición del vapor de alimento y6 = 0.1 kg mol de benceno/kg mol de aire. En cada etapa
se asume que la retención de líquido es M = 20/3 kg mol.

Para el modelamiento matemático de la dinámica de la columna se tienen en cuenta las

siguientes consideraciones:

 El mayor componente de la corriente de líquido es inerte y no se absorbe en la corriente

de gas.

 El mayor componente de la corriente de gas es inerte y no se absorbe en la corriente de

líquido.

 Cada etapa del proceso es una etapa de equilibrio, es decir, el vapor que deja una etapa
esta en equilibrio termodinámico con el líquido que entra.

 La relación en cada etapa entre las composiciones de vapor y líquido en equilibrio es lineal
y el valor de la constante de proporcionalidad es k = 0.5, es decir que:

yi  kxi

El concepto de etapa de equilibrio es importante para el desarrollo de un modelo dinámico de

una columna de absorción. Una etapa de equilibrio está representada esquemáticamente en la
siguiente figura:

L, xi -1 V, yi

Q i

L, xi V, yi +1

Figura 2. Etapa de equilibrio i.

Balance de Materia en la etapa i (2, 3 y 4)

La cantidad total de soluto en la etapa i es la suma del soluto en la fase liquida y la fase
gaseosa (es decir, Mxi  Wyi ). La velocidad de cambio de la cantidad de soluto es
d ( Mxi  Wyi ) / dt . El balance de materia de los componentes en la etapa i pueden ser
escritos (Acumulación = Entrada - Salida):

d ( Mxi  Wyi )
 Lxi 1  Vyi 1  Lxi  Vyi
dt

Como el líquido es mucho más denso que el vapor, podemos asumir que la mayor contribución
al termino de acumulación es Mxi . Entonces:
 d ( Mxi )
 Lxi 1  Vyi 1  Lxi  Vyi
dt

Siendo la retención de líquido M constante tenemos:

dxi L V L V
 xi 1  kxi 1  xi  kxi
dt M M M M

Luego,

dxi L  L  Vk  Vk
 xi 1    xi  xi 1
dt M  M  M

Donde i representa cada una de las etapas intermedias de la columna.

Balance de Materia en la etapa de tope (1)

El balance de materia en el tope (etapa 1) es:

dx1 L  L  Vk  Vk
 xf    x1  x2
dt M  M  M

Donde x f es conocido (composición del liquido de alimento).

Balance de Materia en la etapa de fondo (N)

El balance de materia en el fondo (etapa N) es:

dx N L  L  Vk  Vk
 x N 1    xN  y N 1
dt M  M  M

Donde y N 1 es conocido (composición del vapor de alimento).

Vemos que para nuestra columna de absorción de 5 etapas obtuvimos un sistema de 5

ecuaciones diferenciales lineales de primer orden, las cuales son:
dx1  L  Vk  Vk L
   x1  x2  xf
dt  M  M M
dx2 L  L  Vk  Vk
 x1    x2  x3
dt M  M  M
dx3 L  L  Vk  Vk
 x2    x3  x4
dt M  M  M
dx4 L  L  Vk  Vk
 x3    x4  x5
dt M  M  M
dx5 L  L  Vk  V
 x4    x5  y6
dt M  M  M

Las variables de estado son xi (i = 1 a 5), y las variables de entrada son x f (composición del
líquido de alimento) y y 6 (composición del vapor de alimento). Se asume que las velocidades
de flujo de líquido y vapor son constantes.

 SIMULACIÓN EN EL ESPACIO DE LOS ESTADOS.

En el espacio de los estados podemos escribir nuestro sistema de ecuaciones de una columna
de absorción de 5 etapas, en forma de matrices, asi:

x  Ax  Bu
y  Cx  Du
Donde:

𝑥̇ = representa el vector de las derivadas de la variable de salida con respecto al tiempo.

A= Es una matriz cuadrática de dimensiones iguales al número de variables de salida, en este

caso 5 x 5 cuyos elementos son los coeficientes de cada uno de los términos correspondientes
a las variables de salida en cada una de las ecuaciones.

B= Es una matriz de dimensiones 5 x 2 cuyos elementos son los coeficientes de cada uno de los
términos correspondientes a las variables de entrada en cada una de las ecuaciones del
sistema.

x= representa el vector de las variables de salida.

  ( L  Vk ) Vk 
 M 0 0 0 
  x1   L 
M
 1 
x ( L  Vk ) 0
   M
L Vk
 x    0 0
   M M M   x2   0 0
 x    0 x f 
2
( L  Vk )
 x3   
L Vk

 3  
0 0 0  
   M M M  y6 
 4 
x L ( L  Vk ) Vk   x4   0 0
 x5   0 0  x   V 
 M M M   5   0 M

 L ( L  Vk ) 
0 0 0 
 M M 

u= representa el vector de las variables de entrada.

C= matriz

D=

 0.325 0.125 0 0 0 
 0.2  0.325 0.125 0 0 

A 0 0.2  0.325 0.125 0 
 
 0 0 0.2  0.325 0.125 
 0 0 0 0.2  0.325

0.2 0 
0 0 

B 0 0 
 
0 0 
 0 0.25

- CAMBIO PASO EN LA COMPOSICIÓN DEL VAPOR DE ALIMENTO.

X  AX  BU
Análisis en el estado Dinámico.
Y  CX  DU
En el siguiente desarrollo se considera que en un instante de tiempo despreciable (t = 0
minutos) la composición de la corriente de vapor que entra a la columna cambia de y 6 = 0.1 kg
mol de benceno/kg mol de aire a y6 = 0.15 kg mol de benceno/kg mol de aire, es decir que
experimenta un cambio paso de 0.05 kg mol de benceno/kg mol de aire y se propone el
determinar gráficamente:

 El cambio en el tiempo de la concentración del líquido que sale por el fondo de la columna
y del vapor que sale por el tope

 El cambio en el tiempo de la concentración de líquido en cada etapa de la columna

 Para determinar el perfil gráfico de la concentración del líquido que sale por el fondo y del
vapor que sale por el tope a través del tiempo, la matriz C del modelo en la forma del
espacio de los estados es

 0 0 0 0 1
C 
0.5 0 0 0 0

El producto de la anterior matriz C por el vector X cuyos elementos son las composiciones del
líquido en los platos 1 a 5 da un vector Y cuyos elementos son las variables de salida X5 y Y1.
La matriz D es:

0 0 
D 
0 0 

- Para realizar nuestra simulación, primero, inicializamos SIMULINK en el programa

MATLAB y elegimos New model de la barra de herramientas.
- Sabemos que realizaremos un cambio paso en la composición del vapor de alimento
de 0.05, entonces de la librería sources de simulink escogemos step y arrastramos
hasta nuestro nuevo modelo.
- Le damos doble clic al step instalado y en la ventana Source block parameters: Step
establecemos las condiciones iniciales y finales de nuestro sistema como se muestra
en la figura 3.

Figura 3. Parámetros del cambio paso.

- Luego, de la librería Continuous de simulink escogemos State-Space y lo llevamos a
nuestra ventana de simulación.
- Conectamos el step con el State-Space, luego, hacemos doble click al State-Space
instalado y en la ventana Function Block Parameters: State-Space, colocamos en cada
una de las casillas los elementos correspondientes a cada una de las matrices.
Recordemos que debemos escribirlos entre corchetes y por filas separando cada
elemento con un espacio, y para establecer que los elementos de cada fila han
terminado los separamos con un punto y coma (;) para empezar a escribir los
elementos correspondientes a la siguiente fila; y escribimos las condiciones iniciales
como se muestra en la figura 4.

Figura 4. Parámetros de State-Space.

- Para mostrar las graficas obtenidas instalamos un Demux de la librería de Signal

Routing y dos Scoope de la librería Sinks de simulink para obtener el perfil gráfico de la
concentración del líquido que sale por el fondo y del vapor que sale por el tope.
 Figura 5. Perfil de la composición del liquido en el Fondo de la columna.

Figura 6. perfil de la composición del vapor en el tope de la columna.

- Para determinar el perfil gráfico de la concentración del líquido que sale por el fondo y
del vapor que sale por el tope a través del tiempo, la matriz C del modelo en la forma
del espacio de los estados es
 1 0 0 0 0
0 1 0 0 0

C  0 0 1 0 0
 
0 0 0 1 0
0 0 0 0 1

El producto de la anterior matriz C por el vector X cuyos elementos son las composiciones del
líquido en los platos 1 a 5 da un vector Y cuyos elementos son las variables de salida X1, X2, X3,
X4, X5. La matriz D es una matriz de ceros de dimensiones 5 x 2:

0 0
0 0

D  0 0
 
0 0
0 0

- Realizamos un nuevo sistema en el cual las matrices C y D son las expresadas

anteriormente.

Figura 7. Esquema de la simulación de la columna de absorción con cambio paso en la

composicipon del vapor de alimento.

- La función del Display es mostrar la respuesta ultima para cada una de las etapas.
- Se instalo un botón Gain de la librería Math Operations y se completó con las
siguientes especificaciones:
 Figura 8. Parámetros de la función Gain.

Figura 9. Perfil de la composición del líquido en cada etapa.

La línea amarilla representa la etapa 1, la siguiente la 2 y asi sucesivamente. Esto nos muestra
que la magnitud del cambio en la composición de líquido tiene el valor más grande en la etapa
de fondo y que la velocidad relativa de respuesta es más rápida en cuanto la etapa esté más
cerca al fondo de la columna.

Figura 10.Respuesta de la composiciones de las etapas a un cambio paso de 0.05 de

incremento en la composición del vapor de alimento a t = 0. Las variables desviación
normalizadas son graficadas.

La Figura muestra los perfiles en cada etapa de equilibrio de las concentraciones de la fase
líquida representando en el eje de las ordenadas la variable desviación escalada de dichas
concentraciones como Xi(t)/Xi(100) para mayor facilidad de comparación

- CAMBIO PASO EN LA COMPOSICIÓN DEL LIQUIDO DE ALIMENTO.

En el siguiente desarrollo se considera que en un instante de tiempo despreciable (t = 0

minutos) la composición de la corriente de líquido que entra a la columna cambia de xf = 0.0
kg mol de benceno/kg mol de aceite inerte a xf = 0.025 kg mol de benceno/kg mol de aceite
inerte, es decir que experimenta un cambio paso de 0.025 kg mol de benceno/kg mol de aceite
inerte y se propone el determinar gráficamente:

 El cambio en el tiempo de la concentración del líquido que sale por el fondo de la columna
y del vapor que sale por el tope.

 El cambio en el tiempo de la concentración de líquido en cada etapa de la columna.

En este caso cambian algunas condiciones con respecto al desarrollado anteriormente, en este
caso tendremos que:
 Figura 11. Parámetros de la función Step.

Para obtener la respuesta de la composición del líquido de fondo a un cambio paso de 0.025
de incremento en la composición del liquido de alimento las matrices A, B, C y D permanecen
iguales.

Figura 12. Parámetros de la función State-Space.

De lo cual obtenemos los siguientes perfiles:

Figura 13.Respuesta de la composición del líquido de fondo a un cambio paso de 0.025 de
incremento en la composición del liquido de alimento a t = 0.

Figura 14.Respuesta de la composición del vapor de tope a un cambio paso de 0.025 de

incremento en la composición del líquido de alimento a t = 0
 - Para determinar el perfil gráfico de la concentración del líquido en cada etapa de la
columna a través del tiempo, la matriz C del modelo en la forma del espacio de los
estados es

1 0 0 0 0
0 1 0 0 0

C  0 0 1 0 0
 
0 0 0 1 0
0 0 0 0 1

El producto de la anterior matriz C por el vector X cuyos elementos son las composiciones del
líquido en los platos 1 a 5 da un vector Y cuyos elementos son las variables de salida X1, X2, X3,
X4, X5. La matriz D es una matriz de ceros de dimensiones 5 x 2:

0 0
0 0

D  0 0
 
0 0
0 0

Realizamos un nuevo sistema en el cual las matrices C y D son las expresadas anteriormente.

Figura 15. Esquema de la simulación de una columna de absorción con cambio paso en la
composición del líquido de alimento.
Figura 16.Respuesta de la composiciones de las etapas a un cambio paso de 0.025 de
incremento en la composición del liquido de alimento a t = 0

Figura 17. Respuesta de la composiciones de las etapas a un cambio paso de 0.025 de

incremento en la composición del liquido de alimento a t = 0. Las variables desviación
normalizadas son graficadas.
EJEMPLO 8.5.2: Dinámica de un tren de reactores en serie

En un sistema de tres reactores de mezcla completa conectados en serie no interactuantes

como se observan en la Figura 8.11, un compuesto B es producido a partir de un reaccionante
A mediante una reacción de primer orden en fase líquida del tipo:

AB

Fo V1 F1 F2 F3
V2 V3
CAo K1 CA1 CA2 CA3
K2 K3

Figura 8.11. Reactores CSTR en serie

Un flujo volumétrico constante Fo = 1 pie3/min se mantiene a través del sistema, sin embargo
el volumen en cada reactor es diferente pero constante y sus valores son V1 = 10 pie3, V2 = 8
pie3 y V3 = 5 pie3. Como la temperatura es diferente en cada reactor pero constante, el valor
de la constante de velocidad de reacción en cada reactor es diferente en cada reactor, k1 = 1/4
min-1, k2 = 1/5 min-1, k3 = 1/2 min-1. La concentración molar de la corriente de entrada es Co = 1
lb mol/ pie3 y puede cambiar con el tiempo. Se asume que la densidad de las corrientes se
mantiene constante (independiente de la concentración). La reacción es de primer orden, de
tal manera que la velocidad de reacción de A es dada por la ecuación r  kCA

 Escribir el modelo dinámico del sistema

 Escribir el modelo dinámico en la forma del espacio de los estados
 Determine los valores característicos y evalúe la estabilidad del sistema
 Determinar las concentraciones de A en el estado estacionario en cada reactor
 Mostrar el perfil gráfico de la concentración de A en cada reactor para un cambio paso
en la concentración de la corriente de entrada de 1 lb mol/ pie3

Modelo Matemático

Se lleva a cabo un balance de materia para cada reactor:

 PRIMER REACTOR

El balance molar del componente A:

dC A1 (t )
V1  Fo C Ao (t )  k1V1C A1 (t )  F1C A1 (t ) (1)
dt

El balance molar global de materia:

 d [ 1V1 ]
  o Fo  1 F1 (2)
dt

Siendo el volumen constante y las densidades iguales y constantes, entonces el balance global
se reduce a una igualdad entre los flujos volumétricos de entrada y salida, es decir que
Fo  F1 . Por lo tanto, la ecuación diferencial que expresa el modelo matemático para el
primer reactor se puede escribir como:

dC A1 (t ) 1 K
  C A1 (t )  1 C Ao (t ) (3)
dt 1 1

Siendo la constante de tiempo y la ganancia parámetros dinámicos del reactor número uno
que se calculan con las siguientes expresiones:

V1 10
1    2.86 min
k1V1  Fo (1 / 4) *10  1
Fo 1
K1    0.286
k1V1  Fo (1 / 4) *10  1

 SEGUNDO REACTOR

Los balances de materia para el segundo y tercer tanque corresponden a ecuaciones similares
a las del primer tanque, razón por lo cual se escriben a continuación las ecuaciones finales para
cada uno de ellos y el cálculo de sus respectivos parámetros dinámicos así:

dC A 2 (t ) 1 K
  C A2 (t )  2 C A1 (t ) (4)
dt 2 2

V2 8
2    3.08 min
k2V2  Fo (1 / 5) * 8  1
Fo 1
K2    0.385
k2V2  Fo (1 / 5) * 8  1

 TERCER REACTOR

dC A3 (t ) 1 K
  C A3 (t )  3 C A2 (t )
dt 3 3
V3 5
3    1.43 min
k3V3  Fo (1 / 2) * 5  1
Fo 1
K2    0.286
k3V3  Fo (1 / 2) * 5  1

El modelo matemático en estado dinámico es el conjunto de las ecuaciones diferenciales (3),

(4) y (5) que expresan las variaciones de las concentraciones de A en cada uno de los tanques
con el tiempo.
SIMULACIÓN EN EL ESPACIO DE LOS ESTADOS

En el espacio de los estados podemos escribir nuestro sistema de ecuaciones de una columna
de absorción de 5 etapas, en forma de matrices, asi:

x  Ax  Bu
y  Cx  Du
Donde:

𝑥̇ = representa el vector de las derivadas de la variable de salida con respecto al tiempo.

A= Es una matriz cuadrática de dimensiones iguales al número de variables de salida, en este

caso 5 x 5 cuyos elementos son los coeficientes de cada uno de los términos correspondientes
a las variables de salida en cada una de las ecuaciones.

B= Es una matriz de dimensiones 5 x 2 cuyos elementos son los coeficientes de cada uno de los
términos correspondientes a las variables de entrada en cada una de las ecuaciones del
sistema.

x= representa el vector de las variables de salida.

De acuerdo al sistema de ecuaciones diferenciales (3), (4) y (5), el vector columna x contiene
como elementos las concentraciones en cada uno de los tanques y el vector columna x
contiene los elementos correspondientes a las velocidades de cambio de las concentraciones
en cada tanque. La variable de entrada es “u”, es decir u = CAo

 dC A1 
 dt 
 C A1    dC 
x  C A 2  x   A2 
 dt 
C A3   dC A3 
 dt 
 
La matriz A de los coeficientes de los términos lineales de las variables de salida en las
ecuaciones (3), (4) y (5) y la matriz B de los coeficientes de los términos lineales de las variables
de entrada son las siguientes

 1 / 1 0 0   0.350 0 0 
  
A  K 2 /  2 1 /  2 0    0.065  0.352 0 
 0 K3 /  3 1 /  3   0 0.35  0.70

 K 1 /  1  0.09
B   0    0 
 0   0 
La matriz, C, en la ecuación lineal del modelo, es una matriz diagonal de 3x3 si el modelo solo
expresa como variables de salida las tres concentraciones en cada uno de los tanques. En este
caso la matriz D es igual a cero. Es decir que

1 0 0 0 
C  0 1 0 D  0
0 0 1 0

Entonces, el modelo escrito en la forma del espacio de los estados es:

 0.350 0 0   C A1  0.09


x   0.065  0.352 0  C A2    0 C Ao
 0 0.35  0.70 C A3   0 
1 0 0  C A1  0
y  0 1 0 C A2   0C Ao
0 0 1 C A3  0

- Para realizar nuestra simulación, primero, inicializamos SIMULINK en el programa

MATLAB y elegimos New model de la barra de herramientas.
- Se sabe que se va a realizar un cambio paso al primer reactor, entonces de la librería
sources de Simulink escogemos step y arrastramos hasta el modelo nuevo.
- Le damos doble clic al step instalado y en la ventana Source block parameters: Step
establecemos las condiciones iniciales y finales de nuestro sistema como se muestra
en la figura 18.

Figura 18: Parámetros de la función step.

- Luego, de la librería Continuous de simulink escogemos el bloque State-Space y lo
llevamos a nuestra ventana de simulación.
- Conectamos el step con el State-Space, luego, hacemos doble click al State-Space
instalado y en la ventana Function Block Parameters: State-Space, colocamos en cada
una de las casillas los elementos correspondientes a cada una de las matrices.
Recordemos que debemos escribirlos entre corchetes y por filas separando cada
elemento con un espacio, y para establecer que los elementos de cada fila han
terminado los separamos con un punto y coma (;) para empezar a escribir los
elementos correspondientes a la siguiente fila; y escribimos las condiciones iniciales
como se muestra en la figura 19. Como la matriz C, es una matriz indentidad de 3x3,
podemos representarle de forma abreviada por la palabra eye y entre paréntesis () se
escribe el tamaño de la matriz.

Figura 19: Parámetros de la función state space

- Para obtener el perfil grafico de la función, instalamos un bloque scope, que se

encuentra en la librería commonly. La función state – space se conecta al bloque
scope. El cual nos da el siguiente perfil grafico.
 Figura 20: Respuesta paso para un sistema de tres reactores en serie.

- Realizamos un nuevo sistema en el cual las matrices C y D son las expresadas

anteriormente.

Figura 21: Esquema de la simulación de reactores CSTR en serie

- La función del Display se utiliza para mostrar la respuesta ultimaen cada uno de los
reactores
- Se instalo un botón Gain de la librería Math Operations y se completó con las
siguientes especificaciones:

Figura 22: Parámetros de la función Gain.

Figura 23: Concentración de A en cada reactor

EJEMPLO 8.5.2: Dinámica de un reactor por lotes – Sistema de reacciones en serie

Considere un reactor que opera por lotes, donde no hay corriente de entrada ni de salida. El
reactor se carga, inicialmente, con un volumen de líquido y una concentración inicial de A, en
mol/litro, CAo. Considere un sistema de reacciones en serie donde el componente A reacciona
para formar el componente deseado B y este reacciona para formar un compuesto no deseado
C. Las reacciones son irreversibles y de cinéticas de primer orden, es decir que:

K1 K2

A B

rA  K1C A

rB  K1C A  K 2CB

rC  K 2C B

Siendo K1 y K2, las constantes específicas de velocidad de la reacción de las conversiones de A

en B y B en C, respectivamente. Para la solución del modelo se emplea como concentración
inicial de A, el valor de 0.5 moles/ litro y se asignarán valores de K1 = 1 min-1 y K2 = 0.1 min-1.
Como el producto deseado es B, se quiere simular la variación de su concentración para
observar su máxima producción porque si el tiempo de reacción sobrepasa este óptimo es
posible que todo el producto B sea convertido a C.

MODELO MATEMATICO

Como en el sistema hay tres componentes, se requiere el planteamiento de un balance por

cada componente de la siguiente manera:

Balance de materia de A:

Como hay flujos de entrada ni de salida, entonces la variación de A en el reactor es lo debido a

la descomposición del mismo, lo que se puede escribir como

d [VC A ]
  K1C AV (1)
dt

Como el volumen del sistema es constante, el balance molar del componente A se reduce a:

dC A
  K1C A (2)
dt

Balance de materia de B:

Como en la primera reacción se produce B y en la segunda se descompone y la relación

estequiométrica es de uno a uno, el balance de moles de B es dado por:
 dC B
 K1C A  K 2 C B
dt
(3)

Balance de materia de C:

El componente C se produce en la segunda reacción, por lo tanto la velocidad de su producción

está dada por:

dCC
 K 2CB
dt
(4)

Las ecuaciones diferenciales (2), (3) y (4) constituyen el modelo matemático del sistema de
reacciones en serie en un reactor por lotes.

La ecuación (2) es diferencial lineal cuya respuesta para la concentración de A será de un perfil
exponencial decreciente. Aun cuando es preferible simular la ecuación diferencial de (3)
habiendo resuelto la ecuación (2), la combinación de ellas permite obtener una ecuación
diferencial lineal de segundo orden para la variación en el tiempo de la concentración de B, y
que corresponde a la siguiente:

d 2CB dC
2
 ( K1  K 2 ) B  K1 K 2 C B  0
dt dt
(5)

En forma similar, una combinación de las ecuaciones (4) y (5) permite mostrar que la variación
de la concentración de C con el tiempo es de acuerdo a una ecuación diferencial lineal de
tercer orden y que corresponde a la siguiente:

d 3 CC d 2 CC dCC
3
 ( K1  K 2 ) 2
 K1 K 2 0
dt dt dt

SIMULACION EN EL ESPACIO DE LOS ESTADOS

De acuerdo al sistema de ecuaciones diferenciales (2), (3) y (4), el vector columna x contiene
como elementos las concentraciones de A, B y C y el vector columna x contiene los elementos
correspondientes a las velocidades de cambio de las concentraciones de A, B y C. El
planteamiento es un ejemplo de un sistema con entrada cero.

C A   dC A 
 dt 
x  CB    dC 
x B
CC   dt 
 dCC 
 dt 

La matriz A de los coeficientes de los términos lineales de las variables de salida en las
ecuaciones (2), (3) y (4) es la siguiente
   K 1 0 0   1 0 0
A   K1  K 2 0   1  0.1 0
 0 K 2 0  0 0.1 0

Las matrices B y D no existen y la matriz, C, en la ecuación lineal del modelo, es una matriz
identidad de 3x3 si el modelo solo expresa como variables de salida las tres concentraciones
en cada uno de los tanques. Es decir que

1 0 0
C  0 1 0
0 0 1

Entonces, el modelo escrito en la forma del espacio de los estados es:

 1 0 0  C A1 

x  Ax   1  0.1 0 C A 2 
 0 0.1 0 C A3 
1 0 0  C A1 
y  Cx  0 1 0 C A2 
0 0 1 C A3 

Aplicando transformada de Laplace al modelo del espacio de los estados en términos de las
variables desviación se tiene que:

X (s)  (sI  A) 1 x(0)

Y (s)  C (sI  A) 1 x(0)

Siendo I la matriz identidad de 3x3 y x(0) el vector de las concentraciones iniciales para A, B y
C. La respuesta del reactor se muestra gráficamente en la Figura 8.12. Se observa, que se
alcanza un máximo en la producción de B lo que implica que si se opera el reactor durante un
tiempo mayor la producción del producto deseado es menor y, en cambio, la producción del
subproducto no deseado C aumenta llegándose a tiempos en los cuales todo el B producido se
halla descompuesto y se halla transformado en C.

- Para realizar nuestra simulación, primero, inicializamos SIMULINK en el programa

MATLAB y elegimos New model de la barra de herramientas.
- Se sabe que se va a realizar un cambio paso al primer reactor, entonces de la librería
sources de Simulink escogemos step y arrastramos hasta el modelo nuevo.
- Le damos doble clic al step instalado y en la ventana Source block parameters: Step
establecemos las condiciones iniciales y finales de nuestro sistema como se muestra
en la figura 24
 Figura 24: Parámetros de la función step.

- Luego, de la librería Continuous de simulink escogemos el bloque State-Space y lo

llevamos a nuestra ventana de simulación.
- Conectamos el step con el State-Space, luego, hacemos doble click al State-Space
instalado y en la ventana Function Block Parameters: State-Space, colocamos en cada
una de las casillas los elementos correspondientes a cada una de las matrices.
Recordemos que debemos escribirlos entre corchetes y por filas separando cada
elemento con un espacio, y para establecer que los elementos de cada fila han
terminado los separamos con un punto y coma (;) para empezar a escribir los
elementos correspondientes a la siguiente fila; y escribimos las condiciones iniciales
como se muestra en la figura 25. Como la matriz C, es una matriz indentidad de 3x3,
podemos representarle de forma abreviada por la palabra eye y entre paréntesis () se
escribe el tamaño de la matriz. La matriz B y D son matrices de ceros de 3x3.
 -
Figura 25: Parámetros de la función state space

- Para obtener el perfil grafico de la función, instalamos un bloque scope, que se

encuentra en la librería commonly. La función state – space se conecta al bloque
scope. El cual nos da el siguiente perfil grafico.
 Figura 26: Reactor por lotes.

Figura 27: : Esquema de la simulación de reactores en lote