Disenodeunsensor PDF
Disenodeunsensor PDF
Disenodeunsensor PDF
TESIS
Presentada para obtener el grado de
por
Vicente Pea Caballero
Ingeniero Biotecnlogo
Dirigida por
Dr. Agustn Ignacio Cabrera Llanos
Agradezco a:
Mi familia
Mi mama
Mara Caballero Garca por darme la oportunidad de estar en este camino
tan maravilloso que es la vida y por ayudarme entender que la vida podra ser
lo que so, pero sobre todo por tu ejemplo de lucha da a da en tu forma de
ver y entender la realidad y por darme amor con tan slo una mirada.
A mi asesor
Dr. Agustn Ignacio Cabrera Llanos
Un especial agradecimiento por su ayuda y direccin en este trabajo, por
su ejemplo de disciplina y constancia en el trabajo y por los conocimientos
adquiridos pero sobre todo por sus consejos y amistad.
A mis evaluadores
Dr. Claudio Garibay Origel, Dr. Edgar Salgado Manjarrez, Dr. Juan Sil-
vestre Aranda Barradas y Dr. Sergio Garca Salas gracias por su paciencia y
tiempo dedicado a este trabajo. En especial a los Drs. Edgar y Juan por ser un
ejemplo de pasin al trabajo en las aulas y en la ciencia y por su entrega por la
biotecnologa en UPIBI y por ser pilares de lo que fue, lo que es y lo que ser
nuestra casa de estudios.
A mis amigos
Sergio Medina Llanos, un especial agradecimiento por tu ayuda en aque -
llos momentos tan difciles gracias, siempre te estar agradecido y Naur Flores
Uscanga.
ndice general
3
4 NDICE GENERAL
2. Sistemas no Lineales 19
2.1. Introduccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.2. Sistema no lineal . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.3. Redes neuronales artificiales . . . . . . . . . . . . . . . . . . . . . 21
2.3.1. Modelo de la neurona artificial . . . . . . . . . . . . . . . 23
2.3.2. La funcin de activacin () . . . . . . . . . . . . . . . . 24
2.3.3. Estructuras de las redes neuronales artificiales . . . . . . 25
3. Sensores Virtuales 29
3.1. Introduccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.1.1. Monitoreo de procesos de fermentacin . . . . . . . . . . . 30
3.1.2. Estimadores de estados y parmetros en bioprocesos . . . 31
3.2. Estimacin a travs de elementos de ba-lance . . . . . . . . . . . 34
3.2.1. Estimacin a travs de observadores adaptables . . . . . . 34
3.2.2. Tcnicas de filtrado . . . . . . . . . . . . . . . . . . . . . 36
3.2.3. Redes neuronales artificiales . . . . . . . . . . . . . . . . . 36
3.2.4. Redes neuronales diferenciales . . . . . . . . . . . . . . . . 36
5. Resultados y Discusin 69
5.1. Introduccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
5.2. Datos experimentales: Base de datos . . . . . . . . . . . . . . . . 69
5.2.1. Coeficientes de rendimiento de masa . . . . . . . . . . . . 71
5.3. Aplicacin del software sensor a datos experimentales . . . . . . 71
5.4. Diseo de la red neuronal diferencial . . . . . . . . . . . . . . . . 72
5.5. Diseo de los parmetros del neuro-observador diferencial . . . . 73
5.6. Corrida experimental CDA: Reactor lote alimentado (S. cerevisiae) 74
5.6.1. Estados experimentales . . . . . . . . . . . . . . . . . . . 74
5.6.2. Dinmica de alimentacin al biorreactor . . . . . . . . . . 75
5.6.3. Estimacin de los parmetros experimentales o tericos
bExp
t y kL aExp
t . . . . . . . . . . . . . . . . . . . . . . . . 76
5.6.4. Diseo A de los parmetros C1 y C2 . . . . . . . . . . . 76
5.6.5. Observacin del error XtExp X btExp . . . . . . . . . . . . 77
5.6.6. Norma del error de X bt . . . . . . . . . . . . . . . . . . .
Exp
77
5.6.7. Evolucin de la velocidad especfica de crecimiento esti-
mada, bExp
t .. . . . . . . . . . . . . . . . . . . . . . . . . . 78
5.6.8. Coeficiente volumtrico de transferencia de masa teorco
kL aExp
t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
5.6.9. Diseo B los parmetros C1 y C2 . . . . . . . . . . . . . 79
5.6.10. Observacin del error XtExp X btExp . . . . . . . . . . . . 79
b Exp
5.6.11. Norma del error de Xt . . . . . . . . . . . . . . . . . . . 80
5.6.12. Evolucin de la velocidad especfica de crecimiento esti-
mada, bExp
t .. . . . . . . . . . . . . . . . . . . . . . . . . . 80
5.6.13. Coeficiente volumtrico de transferencia de masa terico
kL aExp
t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
5.6.14. Seleccin final de los parmetros C1 y C2 . . . . . . . . . 81
5.6.15. Observacin del error XtExp X btExp . . . . . . . . . . . . 82
5.6.16. Norma del error de X bt . . . . . . . . . . . . . . . . . . .
Exp
82
5.6.17. Evolucin de la velocidad especfica de crecimiento esti-
mada, bExp
t .. . . . . . . . . . . . . . . . . . . . . . . . . . 83
5.6.18. Coeficiente volumtrico de transferencia de masa terico
kL aExp
t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
5.6.19. Estimacin de estados experimentales . . . . . . . . . . . 85
5.6.20. Identificacin de estados experimentales . . . . . . . . . . 85
5.6.21. Dinmica de la ley de aprendizaje . . . . . . . . . . . . . 86
5.6.22. Observacin de estados experimentales . . . . . . . . . . . 87
5.6.23. Biomasa observada, XtObs . . . . . . . . . . . . . . . . . . . 88
5.6.24. Substrato observado, StObs . . . . . . . . . . . . . . . . . . 89
5.6.25. Oxgeno observado, CtObs . . . . . . . . . . . . . . . . . . 90
5.6.26. Volumen observado, VtObs . . . . . . . . . . . . . . . . . . 90
5.6.27. Estimacin del coeficiente volumtrico de transferencia de
obs
masa kd L at . . . . . . . . . . . . . . . . . . . . . . . . . . 91
5.6.28. Velocidad especfica de crecimiento, bObs
t . . . . . . . . . . 91
5.6.29. Diseo de los parmetros C1 y C2 : Diseo A . . . . . . . 91
6 NDICE GENERAL
Apndices 143
.1. Apndice 1: Medios de cultivo: S. cerevisiae . . . . . . . . . . . . 143
.1.1. Medio de conservacin de la levadura . . . . . . . . . . . 143
.1.2. Medio de cultivo sinttico . . . . . . . . . . . . . . . . . . 143
.1.3. Solucin de minerales . . . . . . . . . . . . . . . . . . . . 143
.1.4. La solucin de vitaminas . . . . . . . . . . . . . . . . . . 144
.1.5. Preparacin del inculo . . . . . . . . . . . . . . . . . . . 144
.2. Apndice 2: Estimador para t . . . . . . . . . . . . . . . . . . . 144
.3. Apndice 3: Estimador para kL at . . . . . . . . . . . . . . . . . . 146
.4. Apndice 4: Eigenvalores . . . . . . . . . . . . . . . . . . . . . . . 148
.5. Apndice 5: Sistema lineal de segundo orden . . . . . . . . . . . . 149
.6. Apndice 6: Corolario A . . . . . . . . . . . . . . . . . . . . . . . 151
.7. Apndice 7: Solucin de la ecuacin de Lyapunov . . . . . . . . . 151
.8. Apndice 8: Prueba de la funcin de Lyapunov para la RND . . . 152
.9. Apndice 9: Teora de Estabilidad de Lyapunov . . . . . . . . . 157
0.4. Resumen
En la presente tesis se introduce el uso de sistemas no lineales (observadores)
para estimar el coeficiente volumtrico de transferencia de oxgeno, kL at , de un
biorreactor.
El modelo utilizado para estimar los estados que no pueden ser medidos
en el sistema de fermentacin (biomasa y substrato principalmente) son las
Redes Neuronales Diferenciales (RND). Con las RND es posible realizar la
estimacin de los estados aun cuando existe incertidumbre en la estructura o
en los parmetros, de hecho, es posible realizar la estimacin en ausencia de un
modelo matemtico del sistema y en presencia de perturbaciones acotadas en el
estado y la salida. La velocidad especfica de crecimiento experimetal estimada,
bExp
t , se estima observando el error entre la biomasa experimental y la biomasa
estimada, XtExp X btExp , utilizando el algoritmo de estimacin propuesto por
Dochain (estimador de t ).
El algoritmo matemtico propuesto en este trabajo para estimar el parmetro
kL at , esta basado en la observacin del error entre la dinmica de evolucin del
oxgeno disponible en lnea (mediciones) y la dinmica de estimacin de oxgeno
estimado, CtExp C btExp . El algoritmo de estimacin contiene en su estructura
informacin que debe de estar disponible en lnea (informacin apriori ) del
proceso de biotransformacin, a saber, biomasa, substrato y oxgeno. Razn por
la cual, el algoritmo diseado se acopla con el modelo de la RND y el estimador
de Dochain (para estimar t ) con lo cual se logra un sistema de estimacin com-
pletamente determinado, es decir, el sensor virtual. Para probar el desempeo de
el observador propuesto se emplean en una base de datos experimentales, para
dos tipos de cultivos (S.cerevisiae y A. vinelandii ), reactor lote alimentado y
reactor lote, respectivamente. Los resultados obtenidos muestran un aceptable
desempeo del estimador para kL at en toda la evolucin de los cultivos. Para A.
vinelandii el porcentaje de error del proceso de estimacin fue de 1.45 % para
el cultivo de A. vinelandii y 7.1 % para S.cerevisiae.
Abstract
Mass transfer from the gas to liquid phase has a decisive importance for the
description of systems involving absorption, chemical reactions and
fermentation. Usually the mass transfer rate is described as proportional to the
concentration gradient, where the proportionality is given by the volumetric
mass transfer coefficient (kLa). The oxygen ( C (t ) ) transfer represents the most
important parameter implied on the design and operation of mixing-sparging
equipment of the bioreactors. It can be described and analyzed by means of
kLa.
In this work we use the Differential Neural Networks (DNN) structure and a non
linear observer, to estimate kLa in the dynamic evolution for the fed-batch model
reactor.
The estimation of the kLa is driven by observation error C (t ) C (t ) . To identify a
Saccharomyces cerevisiae fed-batch fermentation process, it is suggested that
a DNN (dynamic in continuous time) be implemented and then the states are
using to estimate the kLa parameter with the application of a non linear
observer. We assume the oxygens evolution in the aerobic system
fermentation is the only state measure on line. This evolution is used in the non
linear observer in order to estimate the kLa parameter and the results obtained
proved the applicability and efficiency of this methodology because the
corresponding experimental results demonstrate the good performance: the
identification (and observation) error has a small value.
16 NDICE GENERAL
0.5. Justificacin
La transferencia de masa en los sistemas de fermentacin desempea un
papel importante en la economa del proceso. Particularmente en los cultivos
aerobios, donde la transferencia de oxgeno es un factor clave a considerar para
el xito del proceso de biotransformacin a una escala planta piloto o indus-
trial. En un biorreactor tipo tanque agitado la capacidad de transferencia de
oxgeno es caracterizada por el coeficiente volumtrico de transferencia de masa
,kL at . Numerosos estudios se han realizado para determinar el parmetro kL at
en situaciones particulares, utilizando mtodos como oxidacin del sulfito o el
dinmico, entre otros. Tambin, se hace uso de correlaciones empricas. Recien-
temente modelos que utilizan conceptos fundamentales o la aplicacin de cien-
cia artificial, sin embargo, la metodologa existente para estimar el kL at esta
relacionada con grandes esfuerzos de trabajo experimental o altos costos para
instrumentacin, los resultados son particulares y depende de las caractersticas
de operacin para cada caso. Por otro lado el valor estimado es puntal, es decir,
en los casos donde si hay una evolucin completa del parmetro kL at de manera
dinmica a lo largo de todo el proceso implica una fuerte imversin de recursos
en equipo. Por ello el contar con una nueva metodologa que permita estimar el
parmetro kL at para toda la evolucin de la dinmica de biotransformacin de
los substratos a clulas y productos en un biorreactor de manera ms prctica
y econmica , es importante y necesario para el diseo y control de los procesos
a diferentes escalas de operacin.
En los procesos de biotransformacin se tiene conocimiento previo que los
estados pueden tomar un conjunto restringido de valores, de acuerdo a la natu-
raleza del proceso que se estudia; por ejemplo, en el caso de sistemas bilogicos
la velocidad especifca de crecimiento no puede tomar valores negativos; as mis-
mo, ninguna de las concentraciones de los componentes nunca pueden susperar
el valor correspondiente al balance de materia. Por lo tanto, se puede aprovechar
el conocimiento previo para disear observadores que permitan estimar estados
y/o parmetros que no pueden ser medidos en lnea.
0.6. OBJETIVO GENERAL 17
1.1. Introduccin
Una gran variedad de tcnicas de cultivo de microorganismos han sido de-
sarrolladas para cultivos sumergidos, por ejemplo, los cultivos en lote y lote
alimentado, los cuales han sido parte del estado del arte en el campo de las bio-
transformaciones antiguamente y que son usados en la actualidad por muchas
industrias de fermentacin[1]. Las tcnicas de cultivos en continuo son relati-
vamente nuevas y no han encontrado fuertes aplicaciones a nivel industrial, sin
embargo, son una herramienta invaluable en los estudios de laboratorio. En este
captulo se presenta de manera breve los diferentes modelos matemticos (balan-
ces de materia) para cada tipo de biorreactor: reactor lote (batch), lote alimenta-
do (feed-batch o fed-batch) y continuo (continuous) que son el punto de partida
en el anlisis del sistema (biorreactor). Tambin, se describe la importancia del
proceso de transferencia de oxgeno en un sistema de fermentacin aerobio y
sus efectos sobre el desempeo de los cultivos y finalmente las metodologas exis-
tentes en la literatura para estimar el parmetro llamado coeficiente volumtri-
co de transferencia de masa, kL a, siendo este el parmetro ms utilizado para
describir la transferencia de oxgeno en el biorreactor.
1.2. Consideraciones
Dado que el conocimiento de la dinmica del fluido (volumen de la mezcla
de reaccin) en el biorreactor es limitado, no se puede establecer una teora ver-
daderamente satisfactoria sobre este tema. Esta dificultad se debe, al conocimien-
to insuficiente de la estructura del fluido y de las fuerzas intermoleculares. Para
superar esta dificultad, ha sido frecuente enfocar la atencin al criterio de un
mezclado perfecto (estado ideal), sobre todo cuando se consideran pequeos
1
2 CAPTULO 1. ANLISIS DEL SISTEMA: BIORREACTOR
X t =0 X t =0 Ft 2 X t1 Xt2
Ft Et Ft1 Et
1 2
a b
XtE
Et
Ft X t =0 Et 1
Xt 1
X t =0 Ft 2 X t1 XtE
Ft1 2 Et
1
XtR R
c d
Figura 1-1: Arreglos para operacin del biorreactor tipo tanque agitado.
Vt = Vt L + Vt G Vt L
E
kg C
h Et G E G
G ,t Mezcla de Vt G
CO E G ,t reaccin
Ct
kg
Lh COt
Xt
d [ Masa ] St
= acumulacin
dt Pt
F
Ct E
Ct
kg COt E kg
h
Ft F
COt F E
Et E h
S o ,t
S F o ,t E
Xt
Qt Pt E
C F G ,t kg
F
Ft G G F h
CO G ,t
Xt = XtE (1.1)
St = StE (1.2)
Pt = PtE (1.3)
1.4.2. Biomasa
La velocidad de cambio de la concentracin celular en el biorreactor se ex-
presa por el siguiente balance general de masa,
d Ft Et d Xt d
Xt = Xt,0 Xt + Xt Vt (1.4)
dt Vt Vt dt crecimiento Vt dt
d
Donde dt Xt es trmino de acumulacin, es decir, la velocidad de cambio neta
de concentracin
d de clulas en la mezcla de reaccin debido a todas las contribu-
ciones. dt Xt crecimiento contribucin de la multiplicacin celular al cambio de
6 CAPTULO 1. ANLISIS DEL SISTEMA: BIORREACTOR
d
Vt = Ft Et (1.6)
dt
1.4.3. Substrato
El balance de masa general para el substrato,
d Ft Et d St d
St = St,0 St + St Vt (1.7)
dt Vt Vt dt reacci
on Vt dt
1.4.4. Producto
Velocidad de cambio de la concentracin de producto es,
d Ft Et d Pt d
Pt = Pt,0 Pt + Pt Vt (1.8)
dt Vt Vt dt reacci
on Vt dt
Los balances (1.4 -1.8) son las ecuaciones generales de balance . De estos balances
se definen los diferentes modos de operacin del biorreactor: reactor lote, lote
alimentado, y continuo.
d
St = qSt Xt (1.9)
dt
Donde qS = la velocidad especfica de consumo de substrato. Para un cultivo
con un apreciable factor de mantenimiento, entonces
St
qSt = Y S m +m (1.10)
X Ks + St
d
Pt = qPt Xt (1.11)
dt
Donde qPt velocidad especfica de aparicin del producto.
En el caso donde hay sntesis de producto no asociado al crecimiento,
St
qPt = Y P m +B (1.12)
X KS + St
Volumen
d
Vt = Ft = Dt Vt (1.14)
dt
Ft
Donde Dt = Vt , es el trmino de dilucin del cultivo.
Biomasa
d Ft d Xt d
Xt = Xt |t=0 + Xt Vt (1.15)
dt Vt dt crecimiento Vt dt
Considerando la ecuacin (1.14), la ecuacin (1.15) se hace
d Ft d Ft Xt
Xt = Xt |t=0 + Xt (1.16)
dt Vt dt crecimiento Vt
Substrato
d Ft d St d
St = St |t=0 + St Vt (1.17)
dt Vt dt reaccion Vt dt
De manera similar considerando (1.14) la dilucin la ecuacin (1.18) se hace
d Ft d Ft St
St = St |t=0 + St (1.18)
dt Vt dt reaccion Vt
Producto
d Ft d Ft Pt
Pt = Pt |t=0 + Pt (1.19)
dt Vt dt reaccion Vt
d d
En estado estable, ambas dt St y dt Xt son cero. De esta manera de la ecuacin
(1.25),
Dt Xt = t Xt (1.26)
Dt = t
10 CAPTULO 1. ANLISIS DEL SISTEMA: BIORREACTOR
y de la ecuacin (1.26),
Dt ( St |t=0 St ) = Y S t Xt (1.27)
X
Pero Dt = t , entonces
( St |t=0 St ) = Y S Xt (1.28)
X
m Y S Xt + KS Dt Dt Y S Xt
X X
St |t=0 =
m Dt
donde St |t=0 =concentracin de substrato requerida en la corriente de alimentacin,
Ft , para producir una concentracin de clulas, Xt , en la corriente de producto,
Et .
En este trabajo se utilizarn los modelos para un rector lote y lote alimen-
tado.
En los cultivos aerobios se suministra aire estril al biorreactor para cubrir las
necesidades de oxgeno requeridas por los microorganismos. El oxgeno contenido
en las burbujas de aire es transportado de la fase lquida (medio de cultivo)
verificndose el fenmeno de transferencia de oxgeno. La cantidad de oxgeno
que se transfiere a la fase liquida en el tiempo se expresa como la velocidad
volumtrica de transferencia de oxgeno. La transferencia de oxgeno desde una
burbuja hacia el seno del lquido, puede describirse de acuerdo con la teora de
la doble pelcula, desarrollada por Lewis y Withman en 1924.
La teora de la doble pelcula se basa en las siguientes consideraciones: i)
se forma una pelcula en cada lado de la interfase gas-lquido (una pelcula de
gas en el lado del gas y una pelcula de lquido en el lado del lquido) y la
12 CAPTULO 1. ANLISIS DEL SISTEMA: BIORREACTOR
mol
JAz = flujo molar de A en la direccin z m2 h h i
m2
DAB = difusividad (o el coeficiente de difusin) de A en el medio B, h
CA = concentracin de A mol m3
z= distancia en la direccin z [m]
En reactores biolgicos, B es el medio de crecimiento. Una primera aproxi-
macin para un lquido acuoso, puede ser asumido a ser agua con una correlacin
aplicada para la viscosidad del medio actual. Una ecuacin ampliamente usada
para aproximacin de DAB para una solucin diluida de A (no electrolito) en
un solvente B es la ecuacin de Wilke-Chang:
1
7,4 108 (MB ) 2 T
DAB = (1.31)
B VA0,6
donde
= factor de asociacin del solvente
B adimensional
g
MB = peso molecular del solvente B mol
T = temperatura [ K] g
B = viscosidad del solvente B Ls h i
cm3
VA = volumen molar del soluto A mol
1.5. TRANSFERENCIA DE MASA GAS-LQUIDO 13
Pelcula Pelcula
de lquido de gas
Ct
Interfase
Po
Ct
Lquido + oxgeno
disuelto C t Oxgeno
Ct En equilibrio
C t HP o
0
Fase liquida Fase gaseosa
Figura 1-3: Modelo fsico para la transferencia de oxgeno a travs de la pelcula de lquido.
O2
d Et Ct
Ct = kL at [Ct Ct ] OU Rt + (1.34)
dt Vt
En el balance (1.34) la dinmica del oxgeno es mucho ms rpida que la dinmi-
ca de otros procesos relevantes (biomasa, substrato) y la contribucin del tr-
Et Ct
mino de dilucin, , es pequeo comparado con la velocidad de cambio de
Vt
oxgeno disuelto. Como consecuencia en (1.34) se puede expresar como:
d
Ct = kL at [Ct Ct ] OU Rt (1.35)
dt
d
VtG F
CG,t = CG,t FtG EtG CG,t
E
kL at [Ct Ct ] VtG (1.36)
dt
presenta este mtodo con respeto a los de balance de oxgeno en estado estable,
es que no es dependiente de la exactitud de las mediciones de la concentracin
de oxgeno en las corrientes de alimentacin y salida del gas [7],[18],
Adems, de estos mtodos se han reportado una cantidad importantes de
correlaciones para estimar el parmetro kL at en situaciones particulares [19].
En Meenal y col. [20] se estudio el efecto sobre el kL at para diferentes medios y
diferentes condiciones de operacin utilizando correlaciones para medir el kL at
. Estas correlaciones estn relacionadas con:
1. Diferentes geometras del biorreactor.
2. Diferentes propiedades fsicas del medio.
3. Diferentes condiciones de operacin del biorreactor.
Mtodo dinmico. Otras investigaciones se han realizado para desarrollar
nuevas correlaciones para medir el kL at en situaciones muy particulares, vase
[21], [20].
Ct
F t constante Ft 0 F t constante
Pendiente UORt
Ct 1 d C UOR C
d C t UORt
t t
k Lat dt
dt
C critica
t
0
t
Figura 1-5: Mtodo dinmico.
Ct
Ft 0 F t constante
Pendiente 1
k Lat
C critica
t
0
d C UOR
t t
dt
Sistemas no Lineales
2.1. Introduccin
En este captulo se presentan los conceptos bsicos sobre las Redes Neu-
ronales Artificiales (RNA), particularmente para los dos tipos de RNA, Redes
Neuronales Estticas (RNE), y las Redes Neuronales Dinmicas (RNDi). De
este ltimo tipo surgen las Redes Neuronales Diferenciales (RND) que sern
empleadas en este trabajo para estimar los estados del sistema que no son sus-
ceptibles a medicin fsica. La descripcin matemtica de la estructura de la
RND, tanto para el proceso de identificacin de estados como para el proceso de
observacin de estados (RND como observador no lineal) se presentaran en el
captulo 4.
T
xt : = x1,t , ..., xn,t (2.1)
xt Rn ; n = 4 (2.2)
xt : = (Xt , St, Ct , Vt )T (2.3)
2. Acciones de control:
19
20 CAPTULO 2. SISTEMAS NO LINEALES
T
ut : = u1,t , ..., uk,t (2.4)
T
yt : = y1,t , ..., ym,t (2.5)
Factores no controlables
(Ruido y perturbaciones)
Sistema
xt = f (xt , t) + t , xt |t=0 conocida (2.7)
yt = Cxt + t
con
xt Rn , yt Rm t 0
C Rmn , t Rn1 y t Rm1
kCxt k > kt k
2.3. REDES NEURONALES ARTIFICIALES 21
d
bt = F (b
x bt |t=0 conocida
xt , yt , t) , x (2.8)
dt
estructura de una neurona se puede describir por tres partes, un cuerpo celular
o soma, que se caracteriza por contener al ncleo de la clula. El axn que es la
conexin cilndrica sencilla y larga, el cual se extiende del cuerpo de la clula y
lleva los impulsos elctricos a otras neuronas, cada final de un axn se conecta a
una ramificacin estndar muy fina. Cada terminacin tiene un bulbo llamado
sinpsis, donde la red neuronal induce a las neuronas vecinas.
Las redes neuronales artificiales (RNA) son sistemas de informacin inspira-
dos en las redes neuronales biolgicas imitando estas dos caractersticas princi-
pales [26]
Una caracterstica importante es que una RNA puede ser considerada como
una caja negra (o dispositivo de caja negra) que acepta entradas y produce
salidas [12].
Las operaciones que pueden realizarse con una RNA se tienen las siguientes:
Wk1
X1
Funcin de
activacin
Wk 2 uk yk
Seales de
X2 f () Salida
entrada .
.
Wk n Umbral
.
Xn
Pesos
sinpticos
n
X
i = i,j uj (2.10)
j=1
yi = ( i i )
24 CAPTULO 2. SISTEMAS NO LINEALES
Donde:
i j-sima componente de la entrada
i,j uj es el peso de la conexin de la j-sima componente de la entrada con
la neurona ies presentada
i salida del sumador
i umbral de activacin
() funcin de activacin no lineal
yi salida de la neurona i
0.75
0.5
0.25
-5 -2.5 0 2.5 5
x
Figura 2-3: Funcin sigmoide
Wk 2 Wk 2
X2 X2
Wk3 Wk 3
X3 . X3
.
. .
Wk n Wk n
Xn
. . Capa de
Xn salida
Capa de Capa de Capa de Primera Primera
entrada salida entrada Capa oculta Capa oculta
Wk 1
a) G b)
. W1
X1
.
X2 W2 Capa de
G
. P (x) entrada
X3 .
. Capa de
Wp
X p salida
G.
Capa oculta
de funciones de base radial
c)
Figura 2-4: Estructuras de las RNA
La cualidad ms importante de la RNA es su capacidad de aprendizaje,
en este contexto se entiende el aprendizaje como la forma de actualizacin del
valor de los pesos sinpticos. El problema de aprendizaje en las RNA pude
formularse como un problema de aprendizaje de libre modelo, donde el objetivo
es entrenar a una hipotsis RNA definida como hm , donde m denota el nmero
2.3. REDES NEURONALES ARTIFICIALES 27
Sensores Virtuales
3.1. Introduccin
Los sistemas de control son usados para regular una enorme cantidad de
equipos, productos y procesos. Las variables a controlar pueden ser la tempe-
ratura, el movimiento, flujo de vapor, de fluidos, presin de fluidos, tensin,
voltaje y corriente. En el campo de los bioprocesos, uno de los equipos ms
importantes es el biorreactor, que en muchas aplicaciones prcticas, solamente
algunas variables del proceso estn disponibles en lnea, por ejemplo, concen-
tracin de oxgeno disuelto en el biorreactor, temperatura, y la velocidad de
flujo de los gases mientras que los valores de las concentraciones de productos,
reactantes y biomasa estn disponibles a menudo va anlisis. Una alternativa
interesante que evita el uso de un modelo y explota un conjunto limitado de
mediciones es el uso de observadores de estados. El diseo y la aplicacin de
los observadores de estados en bioprocesos han sido una rea activa las pasadas
dcadas.
Muchos conceptos en la teora de control estn basados en el hecho de tener
disponibles los sensores para medir la variable a controlar. De hecho la teora
de control, est a menudo basada en que estn disponibles las seales de retro-
alimentacin. Desafortunadamente, tal suposicin no es vlida en su totalidad.
Los sensores tienen deficiencia que puede degradar un sistema de control.
Hay al menos cuatro problemas relacionados por los sensores. Primero, los
sensores son caros. El costo de un sensor puede sustancialmente elevar el costo
del sistema de control. En muchos casos, el sensor y el cableado son los com-
ponentes ms caros en el sistema de control. Segundo, los sensores y su cable
asociado reducen la confiabilidad del sistema de control. Tercero, algunas seales
son imprcticas para medirlas. Cuarto, los sensores usualmente inducen errores
significativas tales como ruido estocstico, errores cclicos, y respuesta limita-
da. Por tanto, las tcnicas de estimacin de estados juegan un rol importante en
muchas reas de la ciencia, ingeniera y particularmente en biotecnologa, donde
el hardware (sensores) son extermadamente costoso y son usados en condiciones
29
30 CAPTULO 3. SENSORES VIRTUALES
Medicin y Medicin
control
Alimentacin pH
Temperatura pO2
Aireacin Presin
Consumo de Composicin
potencia de exhaustos
y = h(x) (3.2)
b
e=xx (3.5)
d
e = f (b
x + e, u) f (b
x, u) K (b
x) [h(b
x + e) h(b
x)] (3.6)
dt
32 CAPTULO 3. SENSORES VIRTUALES
d
e = [A(b
x) K(b
x)L(b
x)] e (3.7)
dt
Donde A(b
x) y L(b
x) es respectivamente:
f h
A(bx) = , L(b
x) = (3.8)
x x=x x x=x
1. La matriz A(b
x) K(b
x)L(b
x) y su derivada en el tiempo esta acotada:
kA(b
x) K(b
x)L(b
x)k C1 , b
x (3.9)
d
k [A(b
x) K(b x)] k C2 , b
x)L(b x (3.10)
dt
2. Los eigenvalores de A(bx) K(b
x)L(b
x) tienen estrictamente partes reales y
negativas:
Re (i [A(b
x) K(b
x)L(bx)]) C3 < 0, bx
y i = 1 a n, donde n es el orden del sistema descrito por (3.1)
f (x1 , ..., xn ) = 0
f f
f (x1 , ..., xn )[f (x1 , ..., xn )]0 = f (x1 , ..., x10 )+...+ f (xn , ..., xn0 )
x1 0 xn 0
ser medidas, parmetros del modelo o para vencer retardos en las mediciones.
Muchas tcnicas de estimacin han sido propuestas en la literatura. Entre estas
tcnicas, cuatro de ellas han sido reconocidas para tener fuerte potencial en la
estimacin en lnea en los bioprocesos, a saber, (1) estimacin a travs de ele-
mentos de balance; (2) observadores adaptables; (3) tcnicas de filtrado (Filtro
de Kalman, FK, filtro extendido de Kalman, EFK (por sus siglas en ingls));
(4) redes neuronales artificiales [41]. Una quinta tcnica que se propone en este
trabajo son las Redes Neuronales Diferenciales o Dinmicas (RND). Adems,
pueden ser usados otros algoritmos de estimacin [42], en Docahain [40] describe
ocho diferentes algoritmos de estimacin para estados o parmetros en proce-
sos bioqumicos: Observador de Luenberger Extendido (ELO), por sus siglas
en ingls), Observador de Kalman Extendido (EKO, por sus siglas en ingls),
observadores adaptables no lineales, observador asinttico, observador por in-
tervalos [43], EKO para estimar estados y parmetros combinados, observador
con parmetros desconocidos como parmetro de diseo, observador adaptable.
En Borgaerts y col. [44] se utilizan estimadores basados en el concepto de la ob-
servabilidad cannica para sistemas no lineales. El diseo y la aplicacin de un
software sensor se han aplicado en sistemas biolgicos en [45].
d
xe = g(xe , u, ke ) + k1 [(M f (xe , ke )]
dt
d
ke = k2 [(M f (xe , ke )]
dt
Donde k1 y k2 son los parmetros de diseo del observador, los cuales tienen que
ser definidos para asegurar la convergencia de estimacin y el desempeo. En los
sistemas de control el observador de Luenberger se ha aplicado exitosamente,
en la Figura 3-2 se muestra un diagrama a bloques de la estructura de este
observador.
Biorreactor Estados
actuales Sensor
Sistema fsico
+
Error observado
Compensador -
observador
+ Estimacin:
Modelo de estados Sensor
estimacin parmetros
+
Modelo del sistema
4.1. Introduccin
En este captulo se presentan primeramente la metodologa de la experi-
mentacin de los diferentes cultivos que integraron la base de datos experimen-
tales, la cual se utiliz para validar el software sensor diseado y finalmente se
presentan las diferentes tcnicas que fueron utilizadas para disear el estimador
para el coeficiente volumtrico de transferencia de masa, kL at utilizada para el
diseo del software sensor. La metodologa utilizada para disear el estimador
se realiz llevando a cabo dos etapas de trabajo, siendo la primera el uso de las
tcnicas de Redes Neuronales Diferenciales (RND), esta metodologa se utiliz
para identificar (entrenamiento) y observar (estimar) los estados del modelo de
un reactor lote y lote alimentado y a los estados experimentales obtenidos en
procesos de biotransformacin en dos diferentes modos de operacin: reactor
lote alimentado y en lote para los sistemas biolgicos Saccharomyces cerevisiae
y Azotobacter vinelandii respectivamente. En la segunda etapa se detallan los
algoritmos utilizados y diseados para estimar la evolucin del parmetro kL at
en cada tipo de cultivo. Esta ltima etapa se acopla con las tcnicas de RND
y en conjunto se desarrollo el software sensor para estimar el parmetro kL at .
37
38CAPTULO 4. DISEO DEL ESTIMADOR PARA EL PRAMETRO KL AT
Medios de cultivo Los medios de cultivo utilizados para las corridas expe-
rimentales se reportan en el anexo apndice 1.
Biorreactor o
Fermentador de
produccin Sensores:
analizadores y
aparatos de
medicin
Substrato
Sistema
a)
b)
RND
T
Observador x t0 X S C V Identificador
T T
RND RND
yt Xt St Ct Vt yt Xt S t Ct Vt
c) Estimadores
ut
Velocidad
Velocidad Especfica
Especfica de
de Crecimiento
Crecimiento
k Lat
Coeficiente Volumtrico de
Coeficiente Volumtrico de
Transferencia de Oxgeno
Transferencia de Oxgeno
Factores no controlables
(Ruido y perturbaciones)
Factores Informacin
controlables disponible
(Acciones de control) (Salidas)
Biorreactor
d
x = F1 (x) + F2 (x) (4.4)
dt
donde x <n es el vector estados, Rp es el vector de parmetros descono-
cidos y F (x)1 y F (x)2 son (matrices) generalmente no lineales, de funciones del
4.5. REDES NEURONALES DIFERENCIALES 43
d
x = KG (x) + F Q Dx (4.5)
dt
Donde xt , K, G(x), F, Q y D son el vector de estados del proceso gL1 (ms
la temperatura K cuando se considera un balance de energa), matriz de
coeficientes estequiomtricos, el vector
de la velocidad
de reaccin gL1 h1 ,
el vector de velocidad
1alimentacin
gL1 h1 , el vector de velocidad para los
1
1
gases exhaustos gL h y la velocidad de dilucin h respectivemente.
Los balances de masa y energa dados en el modelo descrito por (4.5) esta
compuesto de dos trminos
d
xt bt = f (xt , ut , t) + 1,t , x0 conocida
= x (4.6)
dt
yt = Cxt + 2,t
xt Rn ,yt Rp t 0, ut Rm
C Rpn , t,1 Rn1 y 2,t Rp1
kCxt kik 2,t k
conocimiento a priori , 1,t y 2,t son los ruidos que no pueden ser medidos
(pero si acotados) en la dinmica del estados y la de la salida respectivamente:
xbt = f0 (xt , ut , t | ) + fet + 1,t (4.8)
fet : = f (xt , ut , t) f0 (xt , ut , t | )
: = W10 , W20 , A0 <nn , W10 , W20 <nn , (x) <l1 , (t) <nm
(4.10)
Las funciones de activacin i () y () son seleccionadas como funciones sig-
moides
1
Xn
j (t) : = aj 1 + bj exp cj xj (4.11)
j=1
1
Xn
kl (t) : = aj 1 + bkl exp ckl xl (4.12)
j=1
j = 1, n, k = 1, n, l = 1, m
AT P + P A + P RP + Q = 0 (4.15)
Teorema 1. Consideramos el sistema no lineal (4.6) el modelo del identificador
de la red neuronal cuyos pesos son ajustados por las siguiente ley de aprendizaje
[26]:
d T
W1,t = K1 P t (b
xt ) (4.16)
dt
d T T
W2,t = K2 P t (b
ut ) (b
xt ) (4.17)
dt
t : bt
= xt x (4.18)
bt , yt ) = f (b
F (t, x bt ) ,
xt , t) + K (yt C x K Rnm (4.25)
d
W2,t = k(ut )k2 (1,t + 2 ,t ) [W2,t W2 ] (b xt )T M2,t
xt ) (b (4.29)
dt
con
+T 1
1,t = P N T C0+ (1 1
C + 2 ) C0 N P (4.31)
2 ,t = P N T N 1 P (4.32)
xt ) yet C0+
M1,t = 2P N T (b (4.33)
1 hfT f i 1 hfT f i
Vt := Tt P t + tr W1,t W1,t + tr W2,t W2,t (4.36)
2 2
d d d
= x bt xt (4.37)
dt dt dt
ut Planta xt
Biorreactor +
Yt
RND xt - C
IDENTIFICADOR +
xt xt Yt Yt
-
RND Yt
C
OBSERVADOR
K
Figura 4-4: Diagrama a bloque de la RND. Identificador y observador
a) Entrenamiento de la RND
ut Planta xt
Biorreactor
xt = [ X t Ct Vt ]
T
St +
T
xt xt
xt = X t St C t Vt
-
xt
RND
IDENTIFICADOR
ut Planta Yt = [ X t St Ct Vt ]
Biorreactor
C = [ Xt Vt ] C
T
St Ct
+
C = [ 0 0 1 0]
T
Yt Yt
RND -
C
OBSERVADOR
Yt
K
Yt = [ 0 0 Ct 0]
Figura 4-6: Diagrama a bloques de la estructura del observado
1. Reactor lote
d
Xt = t Xt (4.41)
dt
d
St = k1 t Xt (4.42)
dt
d
Vt = Constante (4.43)
dt
OT Rt ' OU Rt (4.49)
kL at Ct DOtSensor = OU Rt (4.50)
Resolviendo para kL at en la ecuacin (4.50)
1
kL at = OU Rt (4.51)
Ct DOtSensor
Y cuando DOtSensor 0 De la igualdad (4.51) la velocidad mxima de trans-
ferencia de oxgeno esta dada por
C
kL at C = OU Rt (4.52)
C DOtSensor
y cuando DOtSensor tiende a cero en (4.52), se tiene que kL at C = OU Rt . Las
lecturas del sensor, DOt , tiene un tiempo de respuesta Sensor [s] y es modelado
con un sistema de primer orden con
d 1 1
DOtSensor = Sensor DOtSensor + Sensor DOt (4.53)
dt
El parmetro kL at se calculo utilizando los datos experimentales a partir de la
expresin (corrida CDA-IBT)
1 1 1
kL at = OU Rt = t (4.54)
C DOtSensor YO2 C DOtSensor
4.7. ESTADO ESTABLE O CUASI-ESTABLE 53
CDA CDA-1
10 12
Reactor: lote alimentado Reactor: lote alimentado
10
8
8
S(t) [g*L-1]
S(t) [g*L-1]
6
6
4
4
2
2
0 0
0 5 10 15 20 25 30 0 5 10 15 20 25 30
Tiempo [h] Tiempo [h]
S(t)
CDA-2 CDA-IBT
12 25
Reactor: lote alimentado Reactor: lote alimentado
10
20
8
S(t) [g*L-1]
S(t) [g*L-1]
15
6
10
4
5
2
0 0
0 5 10 15 20 25 30 0 20 40 60 80
Tiempo [h] Tiempo [h]
Entonces de (4.45)
1
0' Xt + Dt Sin (4.56)
YXS t
1
Dt Sin = Xt (4.57)
YXS t
54CAPTULO 4. DISEO DEL ESTIMADOR PARA EL PRAMETRO KL AT
Teorema 1 [57]
Para la dinmica bsica del modelo (reactor lote alimentado) descrito por
(4.44-46), bajo las suposicin A1 a A7, la concentracin de biomasa Xt y la
concentracin del substrato, St , estn acotadas como sigue para todo t > 0:
Stmx
0 < St < Stmx 0 < Xt < Xtmx = (4.62)
k1
Prueba
d
1. Cuando Xt = 0, tenemos que dt Xt = 0. Entonces, la conclusin Xt > 0,
para todo t siguiendo de A2.
2. Cuando St = 0, temos que usando A3, A4, A6:
d
St = Dt Sin > 0 (4.63)
dt
3. Cuando St = Stmx , tenemos, usando A1, A3, A4, A5 en (1):
d
St = k1t Xt + Dt [Sin Stmx ] < 0 (4.64)
dt
Entonces, la conclusin St < Stmx para todo t seguida de A7.
4. Se define Zt = k1 + St . cuando Zt = Stmx y usando (4.6), (4.7) y A4:
d
Zt = Dt [Sin Stmx ] < 0 (4.65)
dt
La conclusin Zt < Stmx para todo t seguido de A7. puesto que St > 0 por
(2), claramente se tiene para todo t:
Stmx
Xt < = Xtmx (4.66)
k1
X bt
et = Xt X (4.69)
bt
et = t (4.70)
Entonces, la dinmica del error puede ser deducidad de
d e d d b
Xt = Xt X t (4.71)
dt dt dt
d d d
e =
b (4.72)
dt t dt t dt t
Desarrollando (4.71) con (4.44) y (4.67)
d e et
Xt = t Xt Dt Xt bt Xt + Dt Xt C1 X (4.73)
dt
d e et
et C1 X
Xt = Xt (4.74)
dt
y tambin desarrolla (4.72) con (4.68)
d d et
e = C2 X
(4.75)
dt t dt t
El error del sistema (4.74) y (4.75) se rescribe como
" #
d X et et
et C1 X
Xt
= d et (4.76)
dt et dtt C2 X
4.10. ESTIMADOR EN LINEA PARA LA VELOCIDAD ESPECFICA DE CRECIMIENTO, T 57
T
B(t) = 0 1 (4.80)
d
u =
dt t
El diseo de los parmetros C1 y C2 se visualiza al obtener la ecuacin carac-
terstica de la matriz A, det (I A),
( + C1 ) Xt
det (I A) = det (4.81)
C2
el determinante de la matriz en (4.81) es la ecuacin
2 + C1 + C2 Xt = 0 (4.82)
Esta es una ecuacin de segundo orden que se resuelve de forma general como
sigue p
C1 2 C12 4C2 Xt
i = (4.83)
2
una importante guia es seleccionar los polos reales, esto se logra estableciendo
la siguiente desiguladad (4.84) (vase apndice 4 y 5).
C12
C2 = (4.85)
4Xt
con
d h i
bt
bt = C2 Xt Xt X (4.88)
dt
por tanto, con el sistema (4.71-72) se deduce el error del sistema de manera
semejante al caso anterior y sudinmica de error es
d
x (t) = Xt Ax (t) + B (t) u (t) (4.89)
dt
Con parmetros C1 Lh1 g 1 y C2 L2 h2 g 2 son estrictamente positivos, es
decir, C1 > 0 y C2 > 0.
h iT
x (t) = X et
et (4.90)
C1 1
A= (4.91)
C2 0
T
B(t) = 0 1 (4.92)
y utilizando A de (4.91), el diseo de los polos se obtiene al resolver la ecuacin
caracterstica de la matriz A, det (I A), por lo tanto, C2 ya no es una funcin
de la concentracin de biomasa, como se muestra en la siguiente relacin:
C12
C2 = (4.93)
4
Las propiedades de estabi-lidad y la convergencia del algoritmo de estimacin
(4.87-88) son analizadas bajo las siguientes suposiciones:
d
A8 La derivada de t en el tiempo esta acotada para todo t k dt t k < M1
A9 Xt es estrictamente positiva para todo t: Xt > 0
Estabilidad
Teorema 2.1 (Bastin y Dochain, 1986)
Bajo las suposiciones A1 a A9, existe una constante positiva k0 y k1 , tal que
el vector de x (t) esta acotado como sigue, para todo t:
Prueba de estabilidad
Se considera el sistema autnomo descrito por
d
= Xt A t (4.95)
dt t
4.10. ESTIMADOR EN LINEA PARA LA VELOCIDAD ESPECFICA DE CRECIMIENTO, T 59
Convergencia
El procedimiento del teorema muestra que la estimacin del error et = t b
t
esta acotado con tal que C1 y C2 tenga valores positivos. Pero esto no nos da
informacin sobre la precisin de la estimacin de t . En este prrafo para un
seleccin ligeramente restringida de C1 y C2 , se calcula explcitamente una cota
et .
superior de la precisin asinttica de
Prueba:
Puesto que Xt (s) > > 0 para toda s (suposicin A9), el siguiente cambio
de escala de tiempo:
dt = Xs ds (4.103)
Permite rescribir el sistema (4.89), (4.90), (4.91), (4.92), (4.93) como sigue:
d
x(t) = Ax(t) + B o (t) u(t) (4.104)
dt
con
1 1 0
B o (t) = B (t) = (4.105)
x(t) x(t) 1
d b h i
Xt = bt
bt Xm,t Dt Xm,t + C1 Xm,t Xm,t + X (4.107)
dt
d h i
b
= C2 Xm,t Xm,t + X bt (4.108)
dt t
y el error se hace
d
x(t) = Xm,t Ax(t) + B(t)u(t) (4.109)
dt
con x(t) y Xt definidas de igual forma que el caso anterior (libre de perturba-
ciones):
h iT
x(t) = Xet
et (4.110)
C1 1
A = (4.111)
C2 0
A11 La medicin del ruido y,t esta acotada, pata todo t: k y,t k < M2
A12 Xm,t es estrictamente positiva, para todo t: Xm,t > > 0
A13 La velocidad de dilucin Dt esta acotada, para todo t: 0 < Dt < Dtmx
Estabilidad
Teorema 2.3 (Bastin y Dochain, 1986)
Bajo las suposiciones A1 a A8 , A10 a A13, existen constantes positivas k0 ,
k1 , k2 tales que el vector de error x(t) es acotado como sigue, para todo t:
Convergencia
Teorema 2.4 (B y D,19986)
Bajo las suposiciones A1 a A8, A10 a A13 el error X et y
et estan asimpto-
ticamente acotados como se muestra a continuacin:
e 1 M1
lm sup | Xt | + ( + 2 C1 ) + 3 + 1 + 4 C1 M2 M2
t C2 1 2 1 1 2
(4.115)
1 C1
et |
lm sup | M1 + ( 5 C1 M2 + 6 + 7 C1 ) M2 (4.116)
t C2
con
1 2
2 1 2 2 1 2
, (4.117)
1 1
y 1 (i = 1 a 7) son constantes positivas independientes del diseo de los
parmetros C1 y C2 .
Prueba
De manera similar que en el teorema 2.2 y si se considera el siguientes cambio
de escala:
dt = Xm,s ds (4.118)
y denotando aqui que:
M2
B1 = max {Dtmx + mx
t + C1 (Xm,t + M2 )} (4.119)
62CAPTULO 4. DISEO DEL ESTIMADOR PARA EL PRAMETRO KL AT
1
B2 = {M1 + C2 (Xm,t + M2 )} (4.120)
1 (i = 1 a 7) pueden ser facilmente calculados de (B.8) y B.9) del apndice 6
kg
L at = kL at kdL at (4.126)
dg d d
kL at = kL at kd L at
dt dt dt
Las ecuaciones (4.125) y (4.126) se expresan de una forma muy parecida a la
ecuacin diferencial de estado, ya que el estado del sistema es de primer orden
escritas en funcin de las variables de estado, es decir, la forma matricial de la
ecuacin.
d
x (t) = Ax (t) + B (t) u (t) (4.127)
dt
Donde las matrices x (t), A, B (t) se definen como:
h iT
x (t) = X et kg
L at (4.128)
C1 [C Ct ]
A= (4.129)
C2 0
T
B(t) = 0 1 (4.130)
el diseo de los parmetros C1 y C2 se visualiza al obtener la ecuacin carac-
terstica de la matriz A, det (I A) .
0 C1 [C Ct ]
det (I A) = det (4.131)
0 C2 0
( + C1 ) [Ct C ]
det (I A) = det (4.132)
C2
Y la ecuacin caracterstica se obtiene al resolver la ecuacin (4.132):
2 + C1 + C2 [Ct C ] = 0 (4.133)
Y por lo tanto
C12
C2 = (4.134)
4 [Ct C ]
En (4.134) se tiene que C2 = f (Ct ). Es necesario que la matriz A tenga la forma:
C1 1
A= (4.135)
C2 0
con
AT P + P A = Q
Nota: la sulucin de AT P + P A = Q se presenta en el apndice 7
con positiva ( > 0) ; = [C Ct ] entonces
T
V [x] = [C Ct ] x (t) Qx (t) (4.153)
Por tanto la funcin de Lyapunov candidata descrita por (4.151) cumple con las
propiedades:
Resultados y Discusin
5.1. Introduccin
En este captulo se presentan los resultados y la discusin obtenidos des-
pus de haber aplicado el software sensor desarrollado a datos experimentales.
El diseo del software sensor se realiz por medio de la combinacin de al-
goritmos matemticos (sistemas no lineales): Redes Neuronales Diferenciales
(RND) identificador y observador de estados basados en RND, estimadores no
lineales para la velocidad especfica de crecimiento, t , y para el coeficiente
volumtrico de transferencia de masa, kL at . Los datos experimentales utiliza-
dos, fueron obtenidos de cultivos aerobios para Saccharomices cerevisiae y Azo-
tovacter vinelandii en sistemas de cultivos reactor lote alimentado y lote respec-
tivamente.
El captulo esta organizado como sigue: los algoritmos utilizados para estimar
el parmetro kL at fueron validados primero para los cultivos en reactor lote
alimentado (corridas experimentales identificadas como CDA, CDA-1 y CDA-
2), posteriormente la validacin se realiz para el cultivo en reactor lote (corrida
experimental identificada como CDA-IBT).
Los algoritmos (estimadores y la RND) se disearon en su estrctura para
su mejor desempeo, a travs de un proceso de seleccin de valores en sus
parmetros, es decir, para cada diseo en los parmetros es una nueva seleccin
de valores de estos y no se refiere a un cambio en la estructura de los algoritmos.
69
70 CAPTULO 5. RESULTADOS Y DISCUSIN
bObs
5. Estimacin de la velocidad especfica de crecimiento observada, t : Esta
velocidad se estim al aplicar algoritmo de estimacin descrito por (4.87)
y (4.88) y utilizando la biomasa observada, XtObs .
obs
6. Estimacin del coeficiente volumtrico de transferencia de masa, kd L at :
obs
El parmetro kd
La se estim utilizando el algoritmo descrito por (4.123)
t
bObs
y (4.124) utilizado t y XtObs .
10,2 0 0 0
0 82 0 0
A=
0
o A = diag 10,2 82 168 92,3
0 168 0
0 0 0 92,3
Dado que los pesos de la red neuronal, son ajustados mediante dos ecuaciones
diferenciales, es muy importante establecer de manera adecuada sus compo-
nentes. Como se ha mencinado previamente, la adaptacin de la red se ve al-
tamente influenciada por la seleccin correcta de la solucin de la ecuacin de
Reccati (4.15)
Y con los valores anteriores se resolvi la ecuacin de Ricatti (4.15)
P = diag 24 32 65 53
4,54 5,34 5,40 2,50 4,54 5,34 1,40 2,50
5,60 2,34 4,34 3,50 5,60 2,34 4,34 3,50
W1,t=0 =
3,40
W =
5,40 4,50 4,60 2,t=0 3,40 5,40 4,50 4,60
2,30 4,60 5,50 2,80 2,30 4,60 1,50 2,80
20
Reactor: Lote Reactor: Lote alimentado
18 CDA
16
14
X, S [g*L-1], V [L]
12
10
8 X-Exp
6 S-Exp
4 V-Exp
0
0 8 10 15 20 25 30
Tiempo [h]
1.4
Reactor: Lote Reactor: Lote alimentado
CDA
1.2
0.8
F(t) [Lh-1]
F(t)
0.6
0.4
0.2
0 5 10 15 20 25 30
Tiempo [h]
5.6.5. btExp
Observacin del error XtExp X
12
10
x1, x2 [gL-1]
0
0 5 10 15 20 25 30
Tiempo [h]
5.6.6. btExp .
Norma del error de X
btExp ,
La norma del error para el estado de biomasa experimental estimada, X
se presenta en la Figura 5-4
78 CAPTULO 5. RESULTADOS Y DISCUSIN
2
Reactor: Lote Reactor: Lote alimentado
1.8 CDA
1.6
1.2
0.8
0.6
0.4
0.2
0
0 5 10 15 20 25 30
Tiempo [h]
0.6
y1 [h-1]
0.4
0.2
-0.2
-0.4
-0.6
0 5 10 15 20 25 30
Tiempo [h]
1500
1000
y2 [h-1]
500
-500
-1000
-1500
0 5 10 15 20 25 30
Tiempo [h]
5.6.10. btExp
Observacin del error XtExp X
16
Reactor: Lote Reactor: Lote alimentado
CDA
14
X- (t) - Biomasa experimental
0
0 5 10 15 20 25 30
Tiempo [h]
bt X
Figura 5-7: CDA. Norma del error de X btExp con C1 = 0,07 Lh1 g 1 y C2 = 0,1225 L2 h2 g 2
5.6.11. btExp .
Norma del error de X
La norma del error para el estado de biomasa estimada se presenta en la
Figura 5-8, ntese que se conservan las misma escalas que tiene la Figura 5-4
del Diseo A.
2
Reactor: Lote Reactor: Lote alimentado
1.8 CDA
1.6
Norma del error de x1 [gL-1]
1.4
1.2
0.8
0.6
0.4
0.2
0
0 5 10 15 20 25 30
Tiempo [h]
0.6
0.4
0.2
-0.2
0 5 10 15 20 25 30
Tiempo [h]
bExp
Figura 5-9: CDA. Evolucin de t con C1 = 0,07 Lh1 g 1 y C2 = 0,1225 L2 h2 g 2
800
Reactor: Lote Reactor: Lote alimentado
700 CDA
y2 Coeficiente volumtrico de transferencia de masa estimado
600
500
400
y2 [h-1]
300
200
100
-100
0 5 10 15 20 25 30
Tiempo [h]
La evolucin de kL aExp
t se presenta en la Figura 5-10.
5.6.15. btExp
Observacin del error XtExp X
Como se observa en la Figura 5-11 prcticamente la evolucin de X btExp es
btExp ' 0, se utiliz un marcador para poder
casi igual a XtExp , es decir, XtExp X
distinguir la trayectoria de X bt .
Exp
16
Reactor: Lote Reactor: Lote alimentado
CDA
14 x1 Biomasa estimada
x2 Biomasa experimental
12
10
x1, x2 [gL-1]
0
0 5 10 15 20 25 30
Tiempo [h]
5.6.16. btExp .
Norma del error de X
2
Reactor: Lote Reactor: Lote alimentado
1.8 CDA
1.6
Norma del error de x1 [gL-1]
1.4
1.2
0.8
0.6
0.4
0.2
0
0 5 10 15 20 25 30
Tiempo [h]
La grfica de la norma del error (Figura 5-12) muestra que el error del
estimador para la velocidad especfica de crecimiento experimental ,b
(t)Exp ,
converge rpidamente a cero durante las primeras horas, es decir, X(t)Exp
5.6. CORRIDA EXPERIMENTAL CDA: REACTOR LOTE ALIMENTADO (S. CEREVISIAE)83
La evolucin de bExp
t se muestra en la Figura 5-13, este comportamiento de
bExp
t se tom como referencia (set point), en la seccin donde se muestran los
resultados del desempeo del estimador para kL at . Adems, esta dinmica fue
una funcin del flujo de alimentacin, Ft .
1
Reactor: Lote Reactor: Lote alimentado
0.9 CDA
0.8
0.6
y1 [h-1]
0.5
0.4
0.3
0.2
0.1
0 5 10 15 20 25 30
Tiempo [h]
bExp
Figura 5-13: CDA. Evolucin de t con C1 = 7,7 Lh1 g 1 y C2 = 14,8225 L2 h2 g 2
500
Reactor: Lote Reactor: Lote alimentado
450
CDA y2 Coeficiente volumtrico de transferencia de masa estimado
400
350
300
y2 [h-1]
250
200
150
100
50
0
0 5 10 15 20 25 30
Tiempo [h]
0.1
Norma del error de X(t)-Est
0.08
0.06
0.04
0.02
0
0 5 10 15 20 25 30
Tiempo [h]
W 2 . La dinmica de los estados fue obtenida usando CDA con las condiciones
iniciales del cultivo, cuando se obtuvo un mal desempeo en la convergencia
del error, un nuevo arreglo de las matrices propuestas empricamente fue nece-
sario en la prxima vez y as sucesivamente. Entonces cuando se logr la mejor
convergencia del error, entonces la RND emula la dinmica del cultivo.
16
Reactor: Lote Reactor: Lote alimentado
14 CDA
12
X(t)-Exp
Estados Experimentales
10 S(t)-Exp
O(t)-Exp
8 V(t)-Exp
X(t)-Red
6 S(t)-Red
O(t)-Red
4 V(t)-Red
0 5 10 15 20 25 30
Tiempo [h]
450
Reactor: Lote Reactor: Lote alimentado
400 CDA
350
300
Magnitude de los pesos
250
200
150
100
50
-50
0 5 10 15 20 25 30
Tiempo [h]
6.5
Reactor: Lote Reactor: Lote alimentado
6 CDA
5.5
5
Magnitud de los pesos
4.5
3.5
2.5
1.5
1
0 5 10 15 20 25 30
Tiempo [h]
La dinmica de los pesos W1,t y W2,t convergen a sus mejores valores, este
hecho se observa en la Figura 5-18 para W2,t , es decir, la red se ajusto a sus
5.6. CORRIDA EXPERIMENTAL CDA: REACTOR LOTE ALIMENTADO (S. CEREVISIAE)87
mejores valores. Por otro lado la evolucin de W1,t se ajusta a sus mejores
valores tambin, slo que la dinmica de W1,t se reajusta en el tiempo, como se
puede ver en la Figura 5-17, para el rector lote W1,t se ajusta encontrando el
mejor valor, pero cuando el cultivo pasa de un rector lote a lote alimentado, el
sistema es perturbado y la dinmica misma del cultivo obliga a que W1,t se este
reajustando durante todo el cultivo.
Despus de entonar la RND descrita por (4.14), la dinmica final de la ley
de aprendizaje descrita por (4.16-17) entrego los siguientes resultados para W1,t
y W2,t :
Las matrices de pesos finales (5.1) y (5.2) son para la mejor convergencia del
error en el proceso de entrenamiento de la RND, fueron las utilizadas en la ley
de aprendizaje del NOD descrita por (4.28-29) como W 1 y W 2 ,lo anterior se
define como
W 1 := W1,t (5.3)
W 2 := W2,t (5.4)
Los valores de las matrices W 1 y W 2 se utilizaron en la estructura de la ley de
aprendizaje
450
Reactor: Lote Reactor: Lote alimentado
400 CDA
350
300
200
150
100
50
-50
0 5 10 15 20 25 30
Tiempo [h]
CDA
8
6
Magnitud de los pesos
-2
-4
0 5 10 15 20 25 30
Tiempo [h]
16
Reactor: Lote Reactor: Lote alimentado
CDA
14
12 X(t)- Experimental
X(t)- Neuro-observador
10
Biomasa [g*L-1]
0 5 10 15 20 25 30
Tiempo [h]
16
Reactor: Lote Reactor: Lote alimentado
CDA
14
12 S(t) - Experimental
S(t) - Neuro-observador
10
Substrato [g*L-1]
0 5 10 15 20 25 30
Tiempo [h]
1.6
C(t)- Referencia (set point)
1.4
C(t)- Neuro-observador
Oxgeno [g*L-1]
1.2
0.8
0.6
0.4
0.2
0
0 5 10 15 20 25 30
Tiempo [h]
16
Reactor: Lote Reactor: Lote alimentado
CDA
14
V(t)- Experimental
12 V(t)- Neuro-observador
Volumen [L]
10
4
0 5 10 15 20 25 30
Tiempo [h]
5.6.28. bObs
Velocidad especfica de crecimiento, t .
0.7
Reactor: Lote Reactor: Lote alimentado
0.6
y1 Velocidad especifca de crecimiento Estimada
0.4
y1, y2 [h-1]
0.3
0.2
0.1
0
0 5 10 15 20 25 30
Tiempo [h]
bObs
Figura 5-25: CDA. Evolucin de t con C1 = 0,7 Lh1 g 1 y C2 = 0,1225 L2 h2 g 2
16
Reactor: Lote Reactor: Lote alimentado
14
x1 Biomasa estimada
12 x2 Biomasa estimada RND
10
x1, x2 [gL-1]
0
0 5 10 15 20 25 30
Tiempo [h]
0.25
0.15
0.1
0.05
0
0 5 10 15 20 25 30
Tiempo [h]
Diseo C1 Lh1 g 1 C2 L2 h2 g 2 Figura
B 1.4 0.49 5-28
C 2.1 1.1025 5-29
D 2.8 1.96 5-30
E 3.5 3.0625 5-31
F 4.2 4.41 5-32
bObs
Los resultados obtenidos para el estimador de t con respecto a los valores
de C1 y C2 para cada diseo en la Tabla 5.1, se presentan a continuacin.
94 CAPTULO 5. RESULTADOS Y DISCUSIN
0.7
Reactor: Lote Reactor: Lote alimentado
0.6
y1 Velocidad especifca de crecimiento estimada
y2 Velocidad especifca de crecimiento experimental
0.5
0.4
y1, y2 [h-1]
0.3
0.2
0.1
0
0 5 10 15 20 25 30
Tiempo [h]
bObs
Figura 5-28: CDA. Evolucin de t con C1 = 1,4 Lh1 g 1 y C2 = 0,49 L2 h2 g 2
0.7
Reactor: Lote Reactor: Lote alimentado
0.6
y1 Velocidad especifca de crecimiento estimada
0.4
y1, y2 [h-1]
0.3
0.2
0.1
0
0 5 10 15 20 25 30
Tiempo [h]
bObs
Figura 5-29: CDA. Evolucin de t con C1 = 2,1 h1 y C2 = 1,1025 g 2 L2 h1
5.6. CORRIDA EXPERIMENTAL CDA: REACTOR LOTE ALIMENTADO (S. CEREVISIAE)95
0.7
Reactor: Lote Reactor: Lote alimentado
0.6
y1 Velocidad especifca de crecimiento estimada
0.4
y1, y2 [h-1]
0.3
0.2
0.1
0
0 5 10 15 20 25 30
Tiempo [h]
bObs
Figura 5-30: CDA. Evolucin de t con C1 = 2,8 Lh1 g 1 y C2 = 1,96 L2 h2 g 2
0.7
Reactor: Lote Reactor: Lote alimentado
0.6
y1 Velocidad especifca de crecimiento estimada
y2 Velocidad especifca de crecimiento experiemental
0.5
0.4
y1, y2 [h-1]
0.3
0.2
0.1
0
0 5 10 15 20 25 30
Tiempo [h]
bObs
Figura 5-31: CDA. Evolucin de t con C1 = 3,5 Lh1 g 1 y C2 = 3,0625 L2 h2 g 2
96 CAPTULO 5. RESULTADOS Y DISCUSIN
0.7
Reactor: Lote Reactor: Lote alimentado
0.6
y1 Velocidad especifca de crecimiento estimada
0.3
0.2
0.1
0
0 5 10 15 20 25 30
Tiempo [h]
bObs
Figura 5-32: CDA. Evolucin de t con Diseo F: C1 = 4,2 Lh1 g 1 y C2 = 4,41 L2 h2 g 2
-3
x 10
6
5.5
5
Error
4.5
3.5
3
0.5 1 1.5 2.1 2.5 3 3.5 4 4.5
C1 [h-1]
obs
Para estimar el coeficiente volumtrico de transferencia de masa, kdL at ,con
el algoritmo diseado y descrito por (4.136-37). La velocidad de convergencia del
estimador se logr mediante el diseo apropiado de los parmetros C1 y C2 , es
decir, el mejor desempeo del estimador se obtiene a travs del procedimiento de
entonamiento del estimador. C1 y C2 . Es importante mencionar que la dinmica
obs
de bObs que se utiliz para estimar el kd
t La , corresponde al mejor desempeo t
bObs
en el proceso de estimacin de t .
600
Reactor: Lote Reactor: Lote alimentado
400
y3, y4 [h-1]
300
200
100
0
0 5 10 15 20 25 30
Tiempo [h]
obs
Figura 5-34: CDA. Evolucin de kd
L at para C1 = 100 Lh1 g 1 y C2 = 250 103 L2 h2 g 2
600
Reactor: Lote Reactor: Lote alimentado
400
y3, y4 [h-1]
300
200
100
0
0 5 10 15 20 25 30
Tiempo [h]
obs
Figura 5-35: CDA. Evolucin de kd
L at para C1 = 3000 Lh1 g 1 y C2 = 2250 103 L2 h2 g 2
98 CAPTULO 5. RESULTADOS Y DISCUSIN
600
Reactor: Lote Reactor: Lote alimentado
400
y3, y4 [h-1]
300
200
100
0
0 5 10 15 20 25 30
Tiempo [h]
obs
Figura 5-36: CDA. Evolucin de kd
L at para C1 = 7000 Lh1 g 1 y C2 = 12250 103 L2 h2 g 2
600
Reactor: Lote Reactor: Lote alimentado
400
y3, y4 [h-1]
300
200
100
0
0 5 10 15 20 25 30
Tiempo [h]
obs
Figura 5-37: CDA. Evolucin de kd
L at para C1 = 10000 Lh1 g 1 y C2 = 25000 103 L2 h2 g 2
5.7. CORRIDA EXPERIMENTAL CDA-1: REACTOR LOTE ALIMENTADO (S. CEREVICIAE)99
-3
x 10
9
8
Error Promedio Relativo Absoluto
3
1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
C1 [Lg-1*h-1]
Los resultados obtenidos al comparar los valores para el error relativo ab-
soluto para los diferentes diseos de los parmetros C1 y C2 se muestran en la
Figura 5.38 se puede ver que el mejor diseo fue para C1 = 7000 Lh1 g 1 y
C2 = 12250 103 L2 h2 g 2 .
16
Reactor: lote Reactor: lote alimentado
14 CDA-1
12
10
X, S [g*L-1], V [L]
8
X-Exp
6 S-Exp
V-Exp
4
0
0 5 10 15 20 25
Tiempo [h]
16
Reactor: lote Reactor: lote alimentado
14 CDA-1
Estados Experimentales e Idetificados por la Red
12
10
X(t)- Exp
8
S(t)- Exp
V(t)- Exp
6
C(t)- Setpoint
X(t)- Red
4
S(t)- Red
V(t) -Red
2
C(t)- Red
0 5 10 15 20 25
Tiempo [h]
500
Reactor: lote Reactor: lote alimentado
CDA-1
400
300
Magnitud de los pesos
200
100
-100
0 5 10 15 20 25
Tiempo [h]
CDA-1
6
5
Magnitud de los pesos
1
0 5 10 15 20 25
Tiempo [h]
CDA-1
0.5
F(t)
0.4
0.3
F(t) [Lh-1]
0.2
0.1
-0.1
0 5 10 15 20 25
Tiempo [h]
0.5
Reactor: lote Reactor: lote alimentado
0.4
0.3
0.25
0.2
0.15
0.1
0.05
0
0 5 10 15 20 25
Tiempo [h]
0.5
Reactor: lote Reactor: lote alimentado
0.45 CDA-1
y1 Velocidad especifca de crecimiento estimada
0.4
0.5
a)
0.35 0.4
0.3
0.3
y1 [h-1]
0.2
0.25
0.1
0.2 0
0 5 10 15 20 25
0.15
0.1
0.05
0 5 10 15 20 25
Tiempo [h]
bExp
Figura 5-45:CDA-1. Evolucin de t con C1 = 7,7 Lh1 g 1 y C2 = 14,8225 L2 h2 g 2
5.7.8. btExp
Observacin del error XtExp X
bExp
La evolucin de t se logr minimizando el error de XtExp X btExp por
medio del diseo de los parmetros C1 y C2 en la estructura del estimador
btExp se muestra en la Figura 5.46.
(4.87-88). La observacin del error XtExp X
16
Reactor: lote Reactor: lote alimentado
14 CDA-1
15 a)
12
10
10
x1, x2 [gL-1] 5
8
0
0 5 10 15 20 25
4
x1 Biomasa estimada
x2 Biomasa experimental
2
0
0 5 10 15 20 25
Tiempo [h]
5.7.9. btExp .
Norma del error de X
CDA-1
1.5
a)
Norma del error de X(t) [gL-1]
1
0.1
0.5
0
0 5 10 15 20 25
0.05
0
0 5 10 15 20 25
Tiempo [h]
La evolucin de kL aExp
t calculado se muestra en la Figura 5.48. La evolucin
de kL aExp
t en la etapa de reactor lote este incrementa con la edad del cultivo, y
por tanto por el aumento de la poblacin celular en el biorreactor. Durante el
cultivo en lote alimentado la dinmica de kL aExp
t esta relacionada directamente
con la dinmica del flujo de alimentacin al biorreactor. Finalmente la evolucin
de kL aExp
t tambin se consider como la dinmica de referencia para validar el
estimador para el coeficiente volumtrico de transferencia de masa.
250
y2 [h-1]
200
150
100
50
0
0 5 10 15 20 25
Tiempo [h]
Se utiliz el modelo del NOD para estimar los estados del proceso de fer-
mentacin. La evolucin de biomasa, substrato, oxgeno, y el volumen observa-
dos se muestran en las Figuras: 5-49, 5-50, 5-51 y 5- 52 respectivamente.
106 CAPTULO 5. RESULTADOS Y DISCUSIN
16
Reactor: lote Reactor: lote alimentado
14 CDA-1
12
X(t) Experimental
x(t)-Neuro-observador
10
x1, x2 [gL-1]
8
0
0 5 10 15 20 25
Tiempo [h]
16
Reactor: lote Reactor: lote alimentado
14 CDA-1
S(t)- Experimental
12
S(t)- Neuro-observador
10
Substrato [g*L-1]
0 5 10 15 20 25
Tiempo [h]
16
Reactor: lote Reactor: lote alimentado
15
CDA-1
14
V(t)- Experimental
13 V(t)- Neuro-observador
12
Volumen [L]
11
10
5
0 5 10 15 20 25
Tiempo [h]
-3
x 10
5
Reactor: lote Reactor: lote alimentado
4.5 CDA-1
C(t)- Setpoint
4 C(t)- Neuro-observador
3.5
Oxgeno [g*L-1]
2.5
1.5
0.5
0
0 5 10 15 20 25
Tiempo [h]
500
Reactor: lote Reactor: lote alimentado
CDA-1
400
300
100
-100
0 5 10 15 20 25
Tiempo [h]
12
Reactor: lote Reactor: lote alimentado
10 CDA-1
8
Magnitud de los pesos
-2
-4
0 5 10 15 20 25
Tiempo [h]
5.7.12. bObs
Velocidad especifica de crecimiento observada, t .
0.45 CD-1
0.4
0.35
0.2
0.15
0.1
0.05
0
0 5 10 15 20 25
Tiempo [h]
bObs
Figura 5-55: CDA-1.Evolucin de t con C1 = 7,7 Lh1 g 1 y C2 = 14,8225 L2 h2 g 2
350
Reactor: lote Reactor:lote alimentado
CD-1
300
y3 Coeficiente volumtrico de transferencia de oxgeno estimado
y4 Coeficiente volumtrico de transferencia doxgeno experimental
250
y3, y4 [h-1]
200
150
100
50
5 6 10 15 20 25
Tiempo [h]
obs
Figura 5-56: CDA-1. Evolucin de kd
L at para C1 = 10 103 Lh1 g 1 y C2 = 25 108 . L2 h2 g 2
110 CAPTULO 5. RESULTADOS Y DISCUSIN
16
Reactor: lote Reactor: lote alimentado
CDA-2
14
12
10
F(t) [L*h-1]
6
X- Exp
S- Exp
4
V- Exp
0
0 5 11 15 20 25 30
Tiempo [h]
0.7
Reactor: lote Reactor: lote alimentado
CDA-2
0.6
F(t)
0.5
0.4
F(t) [L*h-1]
0.3
0.2
0.1
-0.1
0 5 11 15 20 25 30
Tiempo [h]
0.5
Reactor: lote Reactor: lote alimentado
Velocidad especfica de crecimiento experimental [h-1]
0.45 CDA-2
0.4
Velocidad especfica de crecimiento experimental
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0
0 5 11 15 20 25 30
Tiempo [h]
0.5
Reactor: lote Reactor: lote alimentado
Velocidad especfica de crecimiento experimental [h-1]
0.45 CDA-2
a)
0.4 0.25
0.2
0.35
0.15
0.3 0.1
0.05
0.25
0
0 5 11 15 20 25 30
0.2
0.15
0.1
0
0 5 11 15 20 25 30
Tiempo [h]
bExp
Figura 5-60: CDA-2. Evolucin de t con C1 = 5,6 Lh1 g 1 y C2 = 7,84 .L2 h2 g 2
5.8.5. btExp
Observacin del error XtExp X
16
Reactor: lote Reactor: lote alimentado
12 CDA-2
14 a)
10
12 8
6
10
4
x1, x2 [g*L-1]
2
8
0
0 5 11 15 20 25 30
x1 Biomasa estimada
2 x2 Biomasa experimental
0
0 5 10 15 20 25 30
Tiempo [h]
CDA-2
1.5
a)
Norma del error de X(t) [g*L-1]
0.1 1
0.5
0
0 5 10 15 20 25 30
0.05
0
0 5 10 15 20 25 30
Tiempo [h]
Figura 5-62: CDA-2. Norma del error de XtExp con C1 = 5,6 Lh1 g 1 y C2 = 7,84. L2 h2 g 2
400
Reactor: lote Reactor: lote alimentado
CDA-2
350
y2 Coeficiente volumtrico de transferencia de oxgeno
300
250
y2 [h-1]
200
150
100
50
0
0 5 11 15 20 25 30
Tiempo [h]
T
bt=0 = 0,4482 10,3 0,0007 6
x
16
Reactor: lote Reactor: lote alimentado
14 CDA-2
X(t)- Exp
S(t)- Exp
12 V(t)- Exp
Estados Identificados por la Red
C(t)- Exp
10 X(t)- Red
S(t)- Red
8 V(t)- Red
C(t)- Red
6
0 5 11 15 20 25 30
Tiempo [h]
400
Reactor: lote Reactor: lote alimentado
350 CDA-2
300
250
Magnitud de los pesos
200
150
100
50
-50
0 5 11 15 20 25 30
Tiempo [h]
7
Reactor: lote Reactor: lote alimentado
CDA-2
6
1
0 5 11 15 20 25 30
Tiempo [h]
16
Reactor: lote Reactor: lote alimentado
CDA-2
14
X(t)- Experimental
12
X(t)- Neuro-observador
10
Biomasa [g*L-1]
0
0 5 11 15 20 25 30
Tiempo [h]
16
Reactor: lote Reactor: lote alimentado
14 CDA-2
12
10
Substrato [g*l-1]
S(t)- Experimental
8 S(t)- Neuro-observador
0 5 11 15 20 25 30
Tiempo [h]
16
Reactor: lote Reactor: lote alimentado
15 CDA-2
V(t)- Experimental
14
V(t)- Neuro-observador
13
12
Volumen [L]
11
10
5
0 5 11 15 20 25 30
Tiempo [h]
-3
x 10
5
Reactor: lote Reactor: lote alimentado
4.5 CDA-2
4
C(t)-Setpoint
3.5
C(t)- Neuro-observador
Oxgeno [g*l-1]
3
2.5
1.5
0.5
0
0 5 10 15 20 25 30
Tiempo [h]
La evolucin de los pesos W1,t y W2,t se muestran en las Figuras 571 y 5-72.
450
Reactor: lote Reactor: lote alimentado
400 CDA-2
350
300
Magnitud de los pesos
250
200
150
100
50
-50
0 5 11 15 20 25 30
Tiempo [h]
10
Reactor: lote Reactor: lote alimentado
CDA-2
8
6
Magnitiud de los pesos
-2
-4
0 5 10 15 20 25 30
Tiempo [h]
1,3 3,4 4,5 4,5
3,5 4,5 2,8 5,6
W1,t=0 = W2,t=0 =
2,8
3,7 1,5 5,6
2,8 3,7 1,5 5,6
5.8.11. bObs
Velocidad especifica de crecimiento observada, t
bObs
Los resultados para t se desarrollaron para cuatro diseos diferentes de
los parmetros C1 y C2 , estos se muestran en las figuras a), b), c) y d) de
la Figura 5-73. El mejor desempeo de para el estimador fue para C1 = 3,5
Lh1 g 1 y C2 = 3,0625 L2 h2 g 2 , este resultado corresponde a la figura c).
Para cada uno de los cuatro diseos se determino la norma del error.
120 CAPTULO 5. RESULTADOS Y DISCUSIN
a) b)
0.4 0.4
Reactor: lote Reactor: lote alimentado Reactor: lote Reactor: lote alimentado
0.3
0.3
y1, y2 [h-1]
y1, y2 [h-1]
0.2
0.2
0.1
0.1
0
-0.1 0
0 5 11 15 20 25 30 0 5 11 15 20 25 30
Tiempo [h] y1 Tiempo [h]
c) y2 d)
0.4 0.4
Reactor: lote Reactor: lote alimentado Reactor: lote Reactor: lote alimentado
0.3 0.3
y1, y2 [h-1]
y1, y2 [h-1]
0.2 0.2
0.1 0.1
0 0
0 5 11 15 20 25 30 0 5 11 15 20 25 30
Tiempo [h] Tiempo [h]
bObs
Figura 5-73: CDA-2. Evolucin de t con C1 = 3,5 Lh1 g 1 y C2 = 3,0625 L2 h2 g 2 para c)
a) b)
1 1
Reactor: lote Reactor: lote alimentado Reactor: lote Reactor: lote alimentado
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 5 10 15 20 25 30 0 5 10 15 20 25 30
Tiempo [h] Tiempo [h]
c) d)
1 1
Reactor: lote Reactor: lote alimentado Reactor: lote Reactor: lote alimentado
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 5 10 15 20 25 30 0 5 11 15 20 25 30
Tiempo [h] Tiempo [h]
bObs
Figura 5-74: CDA-2. Norma del error de t con C1 = 3,5 Lh1 g 1 y C2 = 3,0625 L2 h2 g 2 para c)
a) b)
1000
Reactor: lote Reactor: lote alimentado 400 Reactor: lote Reactor: lote alimentado
800
300
y3, y4 [h-1]
y3, y4 [h-1]
600
200
400
200 100
0 0
0 5 11 15 20 25 30 0 5 11 15 20 25 30
Tiempo [h] Tiempo [h]
y3
c) y4 d)
Reactor: lote Reactor: lote alimentado Reactor: lote Reactor: lote alimentado
400 400
300 300
y3, y4 [h-1]
y3, y4[h-1]
200 200
100 100
0 0
0 5 11 15 20 25 30 0 5 11 15 20 25 30
Tiempo [h] Tiempo [h]
obs
Figura 5-75: CDA-2. Evolucin de kd
L at para C1 = 11 104 Lh1 g 1 y C2 = 30,25 108 L2 h2 g 2 para c)
122 CAPTULO 5. RESULTADOS Y DISCUSIN
Biomasa Substrato
5 25
Reactor: lote Reactor: lote
4 20
X(t)- Exp [g*L-1]
3 S(t)-Exp [g*L-1] 15
X(t)- Exp S(t)- Exp
2 10
1 5
0 0
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
-3
x 10 Oxgeno Volumen
8 4
Reactor: lote Reactor: lote
6 3.5
C(t)- Exp [g*L-1]
V(t)-Exp
4 C(t)-Exp 3
2 2.5
0 2
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
0.35
Reactor: lote
0.1
y1- F)
0.05
0.15 y1- G)
0
0 20 40 60 80 y1- H)
y1- I)
0.1
y1- J)
y1- K)
0.05
y1- L)
0
0 10 20 30 40 50 60 70 80
Tiempo [h]
bExp
Figura 5-77: CDA-IBT. Evolucin de t con diferentes valores de C1 y C2
bExp
En la Figura 5-78 se muestra la evolucin final obtenida para t (Diseo
Exp
J) y la velocidad especfica de crecimiento experimenta, t .
124 CAPTULO 5. RESULTADOS Y DISCUSIN
0.35
Reactor lote
CDA-IBT
0.3
0.15
0.1
0.05
0
0 10 20 30 40 50 60 70 80
Tiempo [h]
bExp
Figura 5-78: CDA-IBT. Evolucin de t con diferentes valores de C1 y C2
22
Reactor: lote
20
CDA-IBT
Estados experimentales identificados por la red
18
16
14
X(t)- Exp
12
S(t)- Exp
10 V(t)- Exp
X(t)- Red
8
S(t)- Red
6 V(t)- Red
0 10 20 30 40 50 60 70 80
Tiempo [h]
Figura 5-79: CDA-IBT.Estados experimentales ( XtExp , StExp , VtExp ) identificados con la RND
5.9. CORRIDA EXPERIMENTAL CDA-IBT: REACTOR LOTE ( A. VINELANDII)125
-3
x 10
Reactor: lote
8
C(t)- Exp
5
C(t)- Red
4
-1
0 10 20 30 40 50 60 70 80
Tiempo [h]
Reactor: lote
500
400
Magnitud de los pesos
300
200
100
-100
0 10 20 30 40 50 60 70 80
Tiempo [h]
7
Reactor: lote
0
0 10 20 30 40 50 60 70 80
Tiempo [h]
7,84 42,22 0,5268 37,45
10,2 141,1 0,2767 72,97
W1,t =
6,686 70,73
2,377 32,13
7,982 108,2 0,6687 63,62
4,543 5,556 3,397 2,332
5,344 2,468 5,395 4,663
W2,t =
1,401
4,401 4,498 1,524
2,505 3,596 4,595 2,8556
Reactor: lote
5
CDA-IBT
4
Biomasa [g*L-1]
X(t)- Experimental
2 X(t)- Neuro-observador
0 10 20 30 40 50 60 70 80
Tiempo [h]
5.9.6. Substrato
Reactor: lote
CDA-IBT
20
15
Substrato [g*L-1]
S(t)- Experimental
S(t)- Neuro- observador
10
0 10 20 30 40 50 60 70 80
Tiempo [h]
5.9.7. Oxgeno
128 CAPTULO 5. RESULTADOS Y DISCUSIN
-3
x 10
Reactor: lote
8 CDA-IBT
Oxgeno [g*L-1]
5
4
C(t)- Experimental
-1
0 10 20 30 40 50 60 70 80
Tiempo [h]
Volumen
5
Reactor: lote
4.5 CDA-IBT
4 data1
V(t)- Neuro- observador
3.5
Volumen [L]
2.5
1.5
0.5
0 10 20 30 40 50 60 70 80
Tiempo [h]
Reactor: lote
500
CDA-IBT
400
Magnitud de los pesos
300
200
100
-100
0 10 20 30 40 50 60 70 80
Tiempo [h]
Reactor: lote
10
CDA-IBT
Magnitud de los pesos
-5
0 10 20 30 40 50 60 70 80
Tiempo [h]
5.9.9. bObs
Velocidad especfica de crecimiento observada, t
bObs
En el proceso de estimacin de t se realiz diseando los parmetros C1
y C2 en la estructura del algoritmo de estimacin descrito por (4.123-124). Los
resultados obtenidos se muestran en la Figura 5-89.
Los resultados que fueron obtenidos en el proceso de estimacin para bObs
t
mostrados en la Figura 5.87 corresponden al diseo de los parmetros C1 y C2
que se muestra en la Tabla 5.6, donde el mejor desempeo fue para el Diseo
bObs
D, este hecho fue confirmado al calcular la norma del error de t con respecto
Exp
bt .
a la seal de referencia
bObs
Tabla 5.6. Diseo de C1 y C2 , estimacin de t en CDA-IBT
130 CAPTULO 5. RESULTADOS Y DISCUSIN
Diseo Figura C1 Lh1 g 1 C2 L2 h2 g 2
A a) 1,4 0,49
B b) 4,2 4,41
C c) 6,3 9,92
D d) 9,8 24,41
a) CDA-IBT b)
0.4
Reactor: lote Reactor: lote
0.3
0.3
0.2
0.2
y1, y2
y1, y2
0.1
0.1
0 0
-0.1 -0.1
0 20 40 60 80 0 20 40 60 80
Tiempo [h] y1 Tiempo [h]
c) y2 d)
Reactor: lote 0.8 Reactor: lote
0.5
0.4 0.6
0.3
y1, y2
y1, y2
0.4
0.2
0.1 0.2
0
0
-0.1
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
bObs
Figura 5-89: CDA-IBT. Evolucin de t para diferentes diseos de C1 y C2
5.9.10. bExp
Norma del error t bObs
t
a) CDA-IBT b)
0.1 0.08
Norma del error de y2-y1 [h-1] Reactor: lote Reactor: lote
0.06
0.04
0.04
0.02
0.02
0 0
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
c) d)
0.25 0.8
Norma del erros de y2-y1 [h-1]
0.15
0.4
0.1
0.2
0.05
0 0
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
bObs
Figura 5-90: CDA-IBT. Norma del error de t
5.9.11. btObs
Observacin del error de XtObs X
btObs a cero, se muestran en la Figura 5-91.
La convergencia del error XtObs X
b
La evolucin de Xt Obs
para los diferentes diseos de los parmetros C1 y C2 de
la Tabla 5.6.
a) CDA-IBT b)
5 5
Reactor: lote Reactor: lote
4 4
Biomasa [g*L-1]
Biomasa [g*L-1]
3 3
2 2
1 1
0 0
-1 -1
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
X(t)- Neuro- observador
X(t)- Estimada
c) d)
5 5
Reactor: lote Reactor: lote
4 4
Biomasa [g*L-1]
Biomasa [g*L-1]
3 3
2 2
1 1
0 0
-1 -1
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
5.9.12. btObs
Norma del error de X
a) CDA-IBT b)
Reactor lote Reactor lote
0.5 0.5
Norma del error [gL-1]
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
c) d)
Reactor lote Reactor lote
0.5 0.5
Norma del error [gL-1]
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
btObs
Figura 5-92: CDA-IBT. Norma del error de X
5.9. CORRIDA EXPERIMENTAL CDA-IBT: REACTOR LOTE ( A. VINELANDII)133
a) CDA-IBT b)
1 1
Biomasa [g*L-1] Reactor: lote Reactor: lote
Biomasa [g*L-1]
0.5 0.5
0 0
-0.5 -0.5
0 2 4 6 8 10 0 2 4 6 8 10
Tiempo [h] X(t)- Neuro- observador Tiempo [h]
c) X(t)- Estimada d)
1 1
Reactor: lote Reactor: lote
Biomasa [g*L-1]
Biomasa [g*L-1]
0.5 0.5
0 0
-0.5 -0.5
0 2 4 6 8 10 0 2 4 6 8 10
Tiempo [h] Tiempo [h]
btObs acercamiento en 10 h t 0 h
Figura 5-93: CDA-IBT.Evolucin X
obs
Tabla 5.7. Diseo de C1 y C2 , estimacin de kd
L at en CDA-IBT
1 1 2 2 2
Diseo Figura C1 Lh g C2 L h g
A a) 1000 250 103
B b) 5000 625 104
3
C c) 10 10 25 106
3
D d) 11 10 3025 105
134 CAPTULO 5. RESULTADOS Y DISCUSIN
CDA-IBT
a) b)
40
Reactor: lote 40 Reactor: lote
30
30
y3, y4
y3, y4
20 20
10 10
0 0
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
y3
y4
c) d)
40 Reactor: lote 40 Reactor: lote
30 30
y3, y4
y3, y4
20 20
10 10
0 0
0 20 40 60 80 0 20 40 60 80
Tiempo [h] Tiempo [h]
obs
Figura 5.94: CDA-IBT. Estimacin de kd
L at
CDA-IBT
a) -3
b)
x 10
0.01 10
Reactor: lote Reactor: lote
8
0.005
6
y3, y4
0
4
y3, y4
-0.005
2
-0.01
0
-0.015 -2
0 20 40 60 80 0 20 40 60 80
Tiempo [h] y3 Tiempo [h]
y4
-3 c) -3
d)
x 10 x 10
10 10
Reactor: lote Reactor: lote
8 8
6 6
y3, y4
y3, y4
4 4
2 2
0 0
-2 -2
0 20 40 60 80 0 20 40 60 80
Tiempo [h]
Para una situacin prctica slo se podr disponer de una referencia prome-
dio para la evolucin del parmetro kL at de acuerdo a una serie de corridas
5.9. CORRIDA EXPERIMENTAL CDA-IBT: REACTOR LOTE ( A. VINELANDII)135
del proceso como tal. Sin embargo, la dinmica de un cultivo entre lote y lote
presenta desviaciones y entonces el comportamiento del parmetro kL at segura-
mente tambin las presentar. Por tanto, el desempeo del estimador propuesto
esta basado en la observacin del error entre el estado disponible en lnea, Ct ,
y el correspondiente estado estimado, CbtObs .
Volumen:
N h
X i2
EV Exp = 100 VbtExp VtExp (5.9)
t
k=1
EV Exp = 1,186 107 L2
t
N h
X i2
%EObs
t
= 100 bObs
t Exp
t (5.14)
k=1
1
%EExp = 0,0001263
t h2
con Exp
t ybObs
t el valor de la velocidad especifica de crecimiento experimental
medida y estimada. Para bObs
t esta dinmica fue estimada asumiendo que Obs
t se
temia en lnea, la que se obtuvo utilizando la dinmica de la biomasa observada,
XtObs .
N
X 2
obs
%Ek obs = 100 kd
L at kL aExp
t (5.15)
L at
k=1
1
%Ek obs = 0,1571
L at h2
obs
con kL aExp
t y kd L at el valor del coeficiente volumtrico de transferencia
de masa medido y su estimado. La desviacin en el proceso de estimacin del
obs
parmetro estimado, kd L at , para cada diseo en la Tabla 5.7 aplicando la
definicin (5.15) los resultados se resume en la Tabla 5.8.
obs
Tabla 5.8 Desviacin de kL aExp
t y kd
L at
%Ek aobs
2 L t
Diseo h
A 2.426
B 0.1571
C 0.3959
D 0.4232
-5
x 10
8
2
Magnitud de la seal
-2
-4
-6
-8
-10
0 10 20 30 40 50 60 70 80
Tiempo [h]
La seal del ruido mostrada en la Figura 5-96 interfiere con la seal que
entrega el sensor para oxigeno disuelto, ambas seales la seal del estado Ct
cuando se considera que no hay perturbaciones y la seal afectada por la per-
turbacin, Cm,t , se muestra en la Figura 5-97.
5.9. CORRIDA EXPERIMENTAL CDA-IBT: REACTOR LOTE ( A. VINELANDII)139
-3
x 10
10
8
C(t) + perturbacin
C(t)- Experimental
6
Oxgeno [g*L-1]
-2
0 10 20 30 40 50 60 70 80
Tiempo [h]
btObs
Figura 5-97 Onservacin de la desviacin entre Cm,t y C
40
35
30
y3 - perturbacin
25
y4
y3, y4
20
15
10
0
0 10 20 30 40 50 60 70 80
Tiempo [h]
obs
Figura 5-98 Evolucin de kd
L at y kL aExp
t en presencia de perturbaciones
obs
En la Figura 5-98 se muestra la evolucin simultanea de kd L at y kL aExp
t ,
para esta dinmica las desvicines entre los valores estimados y los valores ex-
perimentales fue de
1
%Ek aobs = 0,4248 2
L t h
para el resultado anterior el diseo de los parmetros C1 y C1 fueron los mismos
que para el diseo B de la Tabla 5.8, sin embargo, el porciento de la desviacin
fue mayor en el caso de presencia de perturbaciones en la seal de medida para el
140 CAPTULO 5. RESULTADOS Y DISCUSIN
Conclusiones Generales
141
142 CAPTULO 6. CONCLUSIONES GENERALES
143
144 APNDICES
d e h h ii
Xt = t Xt Dt Xt bt
bt Xm,t Dt Xm,t + C1 Xm,t Xm,t + X (5)
dt
d e h i
Xt bt Xm,t Dt Xt + Dt Xm,t C1 Xm,t Xm,t + X
= t Xt bt (6)
dt
.2. APNDICE 2: ESTIMADOR PARA T 145
Desarrollando (7)
d e
Xt bt Xm,t
= t Xt (10)
dt
bt Xm,t
= t (Xm,t y,t ) (11)
et Xm,t t y,t
= (12)
Dt Xt + Dt Xm,t (13)
+Dt y,t (14)
y para (9)
h i
bt
C1 Xm,t Xm,t + X (15)
C1 Xm,t y,t
et
C1 Xm,t X (16)
Con los resultados de (7), (8) y (9), Entonces (6) tiene la forma siguiente:
d e et
et Xm,t t
Xt = y,t + Dt y,t C1 Xm,t y,t C1 Xm,t X (17)
dt
Finalmente se desarrolla (4):
d d h i
e
= t C2 Xm,t Xm,t + X bt (18)
dt t dt
d et
= C2 Xm,t y,t C2 Xm,t X (19)
dt t
Con el desarrollo de (17) y (19) la dinmica de error para el estimador es
" #
d X et et
et Xm,t t y,t + Dt y,t C1 Xm,t y,t C1 Xm,t X
= d e (20)
dt et dt t C2 Xm,t y,t C2 Xm,t Xt
kg d
L at = kL at kL at (28)
ft :
Derivando el error C
df d d b
Ct = Ct C (29)
dt dt dt
Desarrollando cada trmino de (29)
df
Ct = kL at [C Ct ] OU Rt (30)
dt
h h ii
kda
L t [C
Cm,t ] OU Rt + C1 [C
Cm,t ] Cm,t b
Ct (31)
df
Ct = kL at [C Ct ] kd
L at [C Cm,t ] (32)
dt h i
C1 [C Cm,t ] Cm,t C bt (33)
.3. APNDICE 3: ESTIMADOR PARA KL AT 147
df
Ct = kL at [C Cm,t ] + kL at y,t kd
L at [C Cm,t ] (36)
dt h i
C1 [C Cm,t ] Cm,t C bt (37)
df
Ct = kg
L at [C Cm,t ] + kL at y,t (38)
dt h i
C1 [C Cm,t ] Cm,t C bt (39)
df
Ct = kg
L at [C Cm,t ] + kL at y,t (42)
dt
ft C1 [C Cm,t ]
C1 [C Cm,t ] C y,t (43)
d
x(t) = [C Ct ] Ax(t) + B(t)u(t) (49)
dt
Entonces A, B(t) y u(t) tienen la forma siguiente:
C1 [C Cm,t ] [C Cm,t ]
A = (50)
C2 [C Cm,t ] 0
C1 1
= [C Cm,t ] A = [C Cm,t ] (51)
C2 0
k a C1 [C Cm,t ] 0
B(t) = L t (52)
C2 [C Cm,t ] 1
u(t) = d y,t (53)
dt kL at
Por tanto el sistema (49) es el mismo sistema en (44) con lo cual queda demostra-
do.
Z
= | 1 exp (1 s) 2 exp (2 s) | ds (59)
0
Z t
+ | 1 exp (1 s) 2 exp (2 s) | ds (60)
.5. APNDICE 5: SISTEMA LINEAL DE SEGUNDO ORDEN 149
d
x(t) = Ax(t) + b(t)u(t) (61)
dt
Con
x1 (t)
x(t) = (62)
x2 (t)
C1 1
A= (63)
C2 0
b1 (t)
B(t) = (64)
b2 (t)
Asumiendo que b1 (t) y b2 (t) estn acotadas:
y que la matriz A tiene eigenvalores distintos y reales 1 < 2 < 0 (por ejemplo,
C 2
C2 = 4 1 ).
Entonces,
1 1
2B1 1 1 2 1 1 2 B2
lm sup | x1 (t) | + B2 (66)
t 1 2 2 2 1 2
1 + 2
lm sup | x2 (t) | B1 B2 (67)
t 1 2
Prueba
La solucin de x(t) en (61) puede escribirse como sigue:
Z t
x(t) = P J(t)P 1 x(0) + P J (t ) P 1 b ( ) d (68)
0
con :
1 1
P = (69)
1 2
El eigenvector de la matriz de A
1 1 1
P 1 = (70)
1 2 2 1
exp (1 t) 0
J(t) = (71)
0 exp (2 t)
150 APNDICES
C12
C2 = , 0<<1 (89)
4
Entonces las expresiones
2B1
lm sup | x1 (t) | [1 + 2 ] + 3 (90)
C1
2
1
con
1 21
1 21 221
1 = (91)
1+ 21
1+ 21
1 21 221
2 = (92)
1+ 21
4
3 = B2 (93)
C1
Prueba
C1
1 = 1 21 (94)
2
C1
2 = 1+ 21 (95)
2
AT P + P A + Q = 0 (96)
Supngase que tenemos una matriz Q > 0. Para encontrar una matriz P apropi-
ada, se utiliza Z
T
P = eA QeA d (97)
0
Notamos que
d A
e = AeA = eA A (98)
dt
152 APNDICES
d At T T T
e QeAt = eAt A QeAt + eAt Q eAt A (99)
dt
y relacionando (97) con (96) se tiene que
Z Z
T T
AT P + P A = AT eA QeA d + eA QeA d A (100)
0
Z T Z 0
A A T
= e A Qe d + eA Q eA A d (101)
Z0 T
0
T
A A A
= e A Qe + e Q eA A d (102)
0
Z
d T
= eA QeA d (103)
0 dt
El sistema es estable si y solo si eA = 0. De esta manera tenemos
0 T
AT P + P A = eA QeA(0) = (I) Q(I) = Q (106)
Derivando el error
bt xt
t = x (109)
sustituyendo en (108 y 109) en (110)
t = Ab
xt + W1,t (b
xt ) + W2,t (b
xt ) (u (b
xt ))
Axt + W1,t (xt ) + W2,t (xt ) (u (xt )) + f (110)
t = At + W1,t (bxt ) + W2,t (bxt ) (u (b
xt ))
W1,t (xt ) W2,t (xt ) (u (xt )) f
Definimos
e : = (b
xt ) (xt ) (111)
e xt ) (xt )
: = (b (112)
t = At +W1,t (b
xt )W1,t (xt )+W2,t (b
xt ) (u (b
xt ))W2,t (xt ) (u (xt ))f
1. V [x] = 0
2. V [x] > 0, para toda x 6= x
3. V [x] < 0,a lo largo de la trayectoria x = f (x), en el apndice 7 se presenta
la teora de Estabilidad de Lyapunov.
Se parte del hecho que un sistema no lineal x = f (xt , ut , t) siempre se puede
expresar alternativamente por un trmino conocido ms un trmino desconocido
en la siguiente forma:
Definimos un sistema no lineal de la forma
x = f (x, u) (113)
ahora supones que existe una red que se ajusta al sistema con la siguiente forma
x = f (x, u) = Ax + W1 (x) + W2 (x) (ut ) + f (115)
eT e
T Tt DT t (117)
e = (b
x) (x) (118)
T eT e (b
(ut ) (b
xt ) (ut ) xt ) Tt D t kk (u) k2 u (119)
e = (b
x) (x) (120)
e = (b
x) (x)
W1 W1 W1
W2 W2 W2
La ecuacin de Riccati
P A + AT P + P RP + Q0 = 0
t = Ab
x + W1,t (b x) (ut ) [Ax + W1 (x) + W2 (x) (ut ) + f ]
x) + W2,t (b
agrupando trminos
t = A (b
x x) + W1,t (b x) (ut ) W1 (x) W2 (x) (ut ) f
x) + W2,t (b
se derivando
" T # " T #
T
d T
f1 K W
1 f1 K W
f1 + W T 1 f1 +tr W
f2 K W
1 f2 K W
f2 + W T 1 f2
Vt = t P t +t P t +tr W 1 1 2 2
dt
" T # " T #
d T
f 1 f f 1 f
Vt = 2t P t + 2tr W 1 K1 W1 + 2tr W 2 K2 W2
dt
trabajando con 2Tt P t y al sustiruir t se tiene
h i
2Tt P t = 2Tt P At + W1,t
e+W e (ut ) + W
f1,t (xt ) + W2,t f2,t (xt ) (ut ) f
desarrollamos b)
2Tt P At
como T T
t P At = Tt P At
utilizando la siguiente desigualdad matricial
T
XT + XT Y X T 1 X + Y T Y
Y e
=
XT = Tt P W1,t
entonces
2Tt P W1,t eT e
Tt P W1,t 1 +
T 1 T
eT e
t P W1,t W1,t P t +
T
T
t P W1 P t + t Dt t
156 APNDICES
e (ut )
2Tt P W2,t
X = Tt P W2,t
Y = e (ut )
T T
e (ut )
2Tt P W2,t Tt P W2,t + 1 Tt P W2,t + e (ut )
e (ut )
e
Tt P W2 P t + k (ut ) k2
Tt P W2 P t + t D t u
c)
f1,t (xt )
2Tt P W
h i
f1,t (xt ) = 2tr P t (b
2Tt P W f1,t
xt ) W
d)
e (ut )
2Tt P W2,t
h i
e (ut ) = 2tr P t (ut )T (b
2Tt P W2,t f2,t
xt )T W
" T #
f2,t K 1 W
+2tr W f2,t + P t (ut )T (b f2,t
xt ) W
2
factorizando t y Tt
V Tt AT P + P A + P W 1 + W 2 P + D + D + Q0 Tt Tt Q0 t
"" # #
T
+2tr W f1,t K + P t (b
1
xt ) Wf1,t
1
"" # #
T
+2tr f 1 T
W2,t K2 + P t (ut ) (b f
xt ) W2,t
"" T # #
2tr f1,t K 1
W + P t (b f1,t = 0
xt ) W
1
T
f1,t K 1 + P t (b
W xt ) = 0
1
.9. APNDICE 9: TEORA DE ESTABILIDAD DE LYAPUNOV 157
T
f1,t = K1 [P t (b
W xt )]T
"" T # #
+2tr f2,t K 1
W
T
+ P t (ut ) (b f2,t = 0
xt ) W
2
T
f2,t K 1 + P t (ut )T (b
W xt ) = 0
2
h iT
f2,t = K2 P t (ut )T (b
W xt )
T
punov
El concepto de estabilidad de acuerdo a Lyapunov ha tenido muchas apli-
caciones en sistemas de control; de hecho toda la teora de sistemas dinmicos
esta basada en los mtodos de Lyapunov.
Con
AT P + P A + Q = 0 (133)
Es la ecuacin matricial de Lyapunov, si Q es una matriz definida positiva,
entonces el sistema es asimptoticamente estable.
El criterio para el anlisis de estabilidad para una funcin V [x] es el criterio
de estabilidad de Lyaponov.
Una matriz A es una matriz estable para todos los eigenvalores de A, si
y solo si para cualquier matriz (simtrica) semidefinida positiva Q existe una
matriz (simtrica) definida positiva P que satisface
AT P + P A = Q (134)
162 APNDICES
Bibliografa
163
164 BIBLIOGRAFA
[17] C. P. A.-P. Zeng, T.-G-Byun and W.-D. Deckwer, Use of respiratory quo-
tient as a control parameter for optimum oxygen supply and scale-up of
2,3-butanediol production under microaerobic conditions, Biotechnology
and Bioengineering, vol. 44, pp. 11071114, 1994.
[26] E. S. Poznyak, A.S. and W. Yu, Dynamic Neural Networks for Nonlinear
Control: Identification, States Estimation and Trajectory Tracking, N. Jer-
sey, Ed. World Scientific, 2001.