Contro I - Cap. II-A
-
Upload
willfredo-william-walas -
Category
Documents
-
view
17 -
download
0
Transcript of Contro I - Cap. II-A
-
5/22/2018 Contro I - Cap. II-A
1/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
CAPITULO II
II. MODELOS DE SISTEMAS Y APLICACIONES
2.0 MODELOS DE SISTEMAS
Todo sistema es de naturaleza dinmica que proviene de leyes fsicas quegobiernan su comportamiento.El estudio de un sistema de control se realiza a travs de modelos los cualesdeben ser linealizadas para analizar y controlar las relaciones que existenentre las variables de un sistema.El anlisis del modelo se realiza con la ayuda de herramientas matemticasque permita encontrar soluciones que describe la operacin y comportamientodel sistema.Estos modelos utilizan los bloques funcionales en sus anlisis. Estos bloquesfuncionales representan el comportamiento de las relaciones ausa!Efecto" y
de las relaciones de Entrada!#alida del sistema" por lo que estos sistemas queutiliza los bloques funcionales se les conocen como sistemas de parmetrosconcentrados" debido a que cada parmetro se considera en formaindependiente.Todos los modelos tienen similitudes en el comportamiento de los bloquesfuncionales como son en los sistemas mecnicos" elctricos" trmicos y defluidos ya que en base a estos bloques funcionales bsicos y suscombinaciones se producen los modelos matemticos para sistemas fsicosreales.2.1 CLASIFICACION DE LOS MODELOS
$os modelos se clasifican%
a. &odelo Terico.b. &odelo estrictamente Emprico.c. &odelos #emi Emprico.
a. MODELO TEORICO:
'asados en las leyes de la fsica y qumica que emplean en su desarrollobalance de masa" energa" termodinmica" equilibrio qumico" leyes elctricas"
hidrulica" etc.
b. MODELO ESTRICTAMENTE EMPIRICOS
onocidos como modelos de ca(a negra" emplean tcnicas estadsticas querelacionan los datos de Entrada!#alida para cualquier analizas fsico!qumicodel proceso.
c. MODELO SEMI EMPIRICO
onocido como modelo de ca(a gris" son modelos hbridos de los modelostericos y de los modelos estrictamente empricos.
1
-
5/22/2018 Contro I - Cap. II-A
2/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
2.2 MODELOS DINAMICOS Y MODELOS EN ESTADO ESTACIONARIO
Todos los modelos consideran la relacin Entrada!#alida que representaganancia de un sistema y que se le conoce como funcin de transferencia
(F.T),esta dado por la relacin de las salidas y entradas del sistema.
$a funcin de transferencia describe el comportamiento transitorio y estable deun sistema de control. ) esta dado por ecuaciones diferenciales que involucranderivadas de primer orden" segundo orden" tercer orden etc.En este sentido los modelos matemticos de un sistema de control involucranlos modelos dinmicos y modelos de estado estacionario.
2.2.1 MODELOS DE ESTADO ESTACIONARIO
#on aquellos modelos que no varan con el tiempo llamadas tambin estticos"las variables dependientes permanecen constantes en el tiempo.
2.2.2 MODELOS DINAMICOS
#on aquellos en los cuales sus variables dependientes cambian con el tiempo.&ediante tcnicas matemticas y evaluando los parmetros de funcionamientode un modelo dinmico se puede dise*ar el controlador adecuado.
2.3 LINEALIZACION DE LOS MODELOS
omo se ha visto todos los modelos son de naturaleza dinmica que varan enel tiempo" esas variaciones son tan peque*as que pueden despreciarse yobviarse para los clculos de ingeniera" en otros casos no es posibledespreciarse ya que influyen sobre los resultados.El ob(etivo de la linealizacin del modelo es suponer que las variables sedesvan poco dentro de las condiciones de operacin.Este comportamiento no lineal del sistema puede linealizarse aplicando
algunas condiciones y herramienta matemtica como son%
a. $inealizacion utilizando las matrices +acobianas.b. $inealizacion utilizando el desarrollo de la serie de Taylor.
2.4 LINEALIZACION UTILIZANDO MATRICES JACOBIANAS
Este mtodo" supone que las variables se desvan poco y adems la dinmicade un proceso ya sea MIMO ,&-ltiple nput y &ultiplo /upout0" SISO ,/nenput y /ne /upout0" en el tiempo continuo puede ser representado en elespacio de estado mediante dos con(untos de ecuaciones diferenciales
ordinaria de primer orden" que se les denomina ecuaciones de estado y desalida tales como%
2
Entrada
SalidaTF =.
-
5/22/2018 Contro I - Cap. II-A
3/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
1111.. (1
2onde
F, h % 3unciones vectoriales de orden n y r.X % 4ector de estado de orden n.u % 4ector de control de orden m.y % 4ector de salida de orden r.v % 4ector de disturbio en los estados de orden n.w % 4ector de disturbio en los estados de orden r.
t %Tiempo.
5ara linealizar al sistema en los puntos de equilibrio en el espacio de estadodel proceso descrito en la ecuacin (1),se determinara cuando el estado del
proceso no cambia debido a la accin de la fuerza de control U, porconsiguiente los puntos de equilibrio X se determina de%
on este criterio y considerando que la ecuacin (1! es invariante en eltiempo" podemos obtener su representacin linealizada de la forma.
2onde%
A : &atriz de estado ,n x n0B : &atriz de control ,n x m0C : &atriz de salida de los estados ,r x n0D : &atriz de salida de las entradas ,r x m0E : &atriz de disturbio de los estados ,n x n0
F : &atriz de disturbio de las salidas ,r x r0
#i consideramos nula la presencia de las perturbaciones es decir " # $ # O) adems operamos alrededor del estado de equilibrio (% ! U" se puedenobtener las matrices A" B" Cy Dal evaluar las matrices +acobianas.
3
),,,,(
),,,(
twuxhy
tvuxfXo
==
( ) 0, == UXfx
FwDuCxy
EvBuAxxo
++=++=
-
5/22/2018 Contro I - Cap. II-A
4/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
6 7
1
1
2
1
1
X
f
X
f
X
f
n
),(2
2
2
2
1
2
1
................
................
................
UXn
nn
n
n
X
f
X
f
Xf
Xf
X
f
X
f
' 7
1
1
2
1
1
U
f
U
f
U
f
n
( )UXn
nn
n
n
U
f
U
f
U
f
U
f
aU
f
U
f
,2
2
2
2
1
2
1
................
................
................
4
-
5/22/2018 Contro I - Cap. II-A
5/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
7
1
1
2
1
1
x
h
x
h
x
h
n
( )UXn
nn
n
n
x
h
x
h
xh
xh
x
h
x
h
,2
2
2
2
1
2
1
................
...............
...............
2 7
1
1
2
1
1
U
h
U
h
U
h
n
( )UXn
nn
n
n
U
h
U
h
U
h
U
h
U
h
U
h
;2
2
2
2
1
2
1
....................
...................
..................
PROBLEMA
2ado el siguiente sistema de ecuaciones.
5
-
5/22/2018 Contro I - Cap. II-A
6/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
121
0
1 82 uXSenXX ++= " 121 ==uu
3
2
11
0
2 )32(3 XXXX +=
23
0
3 2 uXX +=
#e pide%a. $os puntos de equilibrio , 321 ,, xxx 0.b. $a matriz 6" '" " 2" del sistema linealizado" considerando los puntos de
equilibrio determinado en el caso ,a0.
SOLUCI&N.
121
0
1 82 uXSenXX ++= .................. , f8 0
3211
0
2 )32(3 XXXX += ................... , f9 0 " 121 ==uu
23
0
3 2 uXX += .................... , f: 0
a. 'aa)*+ +, -)/+, * b+.
082 121 =++ uXSenX .................... ,80
0)32(3 32
11 =+ XXX .................... ,90
02 23 =+ uX .................... ,:0
D a cac) ( 3
5.02
1
2 23 === u
X
D a cac) ( 2
0323 32131 =+ XXXX
0233 332
11 =+ XXXX
015.132
11 =+ XX
a
acbbX
XX
2
4
0135.1
2
1
1
2
1
=
=+
6
-
5/22/2018 Contro I - Cap. II-A
7/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
( ) ( ) ( )( )
29.23
1533
693
29.03
153
3
693
5.12
15.143)3(
1
1
2
1
==+=
=+
=++
=
+=
X
X
X
D a cac) ( 1 .
8
2 112
uXXSen
+=
8
129.022
+=XSen
1975.08
58.12 ==XSen
( ) == 39.111975.0.2 SenArcX
1988.0360
239.112 ==
X
( )4475.0
8
129.22
8
2 112 =
+=
+=
uXXSen
( ) 464.0360
258.2658.264475.0.
22 ==== XSenArcX
$os puntos de equilibrio son%
, ;.9
-
5/22/2018 Contro I - Cap. II-A
8/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
( ) ( ) ( )
( )[ ] ( )[ ] ( )[ ]
( ) ( ) ( )
+
+
+
++
+
++
++
++
=
=
=
232323
3
2
113
2
113
2
11
121121121
3
3
2
3
1
3
3
2
2
2
1
2
3
1
2
1
1
1
222
323323323
828282
321
321
321
uxuXux
XxXXxxXxx
usenXxusenXxusenXX
X
f
X
f
X
f
X
f
X
f
X
f
X
f
X
f
X
f
A
XXX
XXX
XXX
=
=
== 200
75.108.3
084.72
200
32063
0cos82
1988.0
29.0X
2
131
2
2
1
X
XXX
X
A
8
-
5/22/2018 Contro I - Cap. II-A
9/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
=
=10
00
01
2
3
1
3
2
2
1
2
2
1
1
1
uf
uf
u
f
u
fu
f
u
f
B
PROBLEMA
2ado el siguiente sistema no lineal.
13
2
21
0
3 uXXX += " 121 ==uu
212
0
2 XSenXX =
223
0
34 uXX =2
1Xy=
a. $os puntos de equilibrio ( 321 x,x,xb. $as matrices 6" '" " 2 del sistema linealizado en los puntos de
equilibrio" determinados en , a 0
SOLUCI&N.
$as ecuaciones de estado%
13
2
21
0
3 uXXX += 111 , f 8 0
212
0
2 XSenXX = 111. , f 90 " 121 ==uu
223
0
34 uXX = 111. , f :02
1Xy= 111. , h 0
a. 'aa)*+ +, -)/+, * b+.
9
-
5/22/2018 Contro I - Cap. II-A
10/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
03X 132
2 =+ uX 111. ,80
02 21 = XSenX 111.. ,90 034 22 = uX 111.. ,:0
D a cac) (3
22 4
3uX = 12 =u
75.02 =X
D a cac) (2
21 X2
1SenxX =
)9718.42(2
1
2
360.75.0
2
1 01 Sen
RadxRadSenX =
=
3408.06816.02
11 == xX
D a cac) ( 1 .
22
213 )75.0()1(33 +=+= XuX
5625.33=X
L+, -)/+, * b+ ,+):
( ) )5625.3,75.0,3408.0(,, 321 =xxx
b. 'aa)*+ a, 9a/c, * ,,/9a )a a*+ a**+* +, -)/+, * b+.
10
DuCXy
BuAXX
+=+=
-
5/22/2018 Contro I - Cap. II-A
11/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
( ) ( ) ( )
[ ] [ ] [ ]
( ) ( ) ( )
++
+
=
=
222222
212121
122213
2213
22
3
3
2
3
1
3
3
2
2
2
1
2
3
1
2
1
1
1
343434
222
33
321
321
321
uXuXuX
SenXXSenXXSenXX
uXXuXXuXX
Xf
Xf
Xf
X
f
X
f
X
fX
f
X
f
X
f
A
XXX
XXX
XXX
=
=
040
07317.02
15.10
040
02
120
2
2
CosX
X
A
11
-
5/22/2018 Contro I - Cap. II-A
12/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
( )
( ) ( )
( ) ( )
( ) ( )
=
++
=
=30
00
03
3434
22
33
22
2
22
1
21
2
21
1
132
2
132
1
,,2
3
1
3
2
2
1
2
2
1
1
1
321
XuXu
SenXXuSenXXu
XXXXu
u
f
u
f
u
f
u
f
u
f
u
f
B
XXX
( ) ( )321 13
2
2
2
1
2
321
002)()()(
XXXX
X
X
X
X
X
X
X
h
X
h
X
hC =
=
=
( )006816.0=C
( )00)()(
2
2
1
1
2
1
21 =
=
=u
X
u
Xu
h
u
h
D
12
-
5/22/2018 Contro I - Cap. II-A
13/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
+
=
2
1
3
2
10
3-0
00
03
040
00.7315-2
11.5-0
u
u
x
x
x
X
13
-
5/22/2018 Contro I - Cap. II-A
14/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
[ ] [ ]
+
=2
1
3
2
1
00006816.0u
u
x
x
x
y
PROBLEMA
2ado un proceso representado por el siguiente sistema no lineal.
1221
0
1 3,2 += xxxx2
21
0
2 3 uxx =
y 7 212
1 32 x
'aa.
a. $os puntos de equilibrio.b. $as matrices y (acobiano 6" '" " 2 del sistema linealizado alrededor de los
puntos de equilibrio.,onsiderar el valor ms peque*o de alguna variable en el punto de equilibrio"en caso de que tuvieran dos races..
SOLUCI&N.
1221
0
1 3,2 += xxxx f8
14
-
5/22/2018 Contro I - Cap. II-A
15/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
2
21
0
2 3 uxx = f9
y 7 212
1 32 x f: 121 ==uu
a. 'aa)*+ +, -)/+, * b+.
(1)03.2 1221 =+ XXX
)(2032
21 =+ X
D+)*:
121 ==uu
2e la ecuacin ,90
33.03
1133 11
2
21 ==== XXX
R9-aa)*+ ) a cac) (1:
( )1366.0
0333.02
22
122
=
=+
XX
XX
E;a)*+ a ca*a*+.
( )222 1344.0 = XX
0144.69
16944.0
2
2
2
2
2
22
=+
+=
XX
XXX
D,a+a)*+ ,a*a.
;.@>B>
2X
;.99B>
D+)* +, -)/+, * +-ac) , D ,a)*+ (0.33! 0.22=7.
15
-
5/22/2018 Contro I - Cap. II-A
16/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
( )
=
=
=03
3086.29546.0
03
32
2278.0,33.0
2
12
,,2
2
1
2
2
1
1
1
21
X
XX
FF
FF
A
XX
=
=
=20
01
20
01
22
2
1
2
2
1
1
1
FF
FF
B
( ) ( )032.104 12
1
1
1 ==
= XXX
C
( )312
1
1
1 =
=
C
2.4.1 MATRIZ DE TRANSFERENCIA
#e aplica para procesos multivariados" en el espacio de estado MIMO.omo se ha visto el espacio de estado MIMOesta formado por las ecuacionesde estado y de salida%
6plicando la transformada de laplace y considerando las condiciones inicialesigual a cero" se obtiene.
de (1!factorizando y ordenando.
( ) )()()(. ssss UBXAIS =
#e agrega la matriz identidad para poder operar.
16
SalidadeEcuacDuCxy
EstadodeEcuacBuAxxo
....
....
+=
+=
)2(.......
)1(......
)()()()()(
)()()()()(
sssss
sssss
UDXC
UBXASX
+=
+=
-
5/22/2018 Contro I - Cap. II-A
17/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
2espe(ando C ( )s :
Demplazando (3) en(2)
Delacionando la salida con la entrada se obtiene la matriz de transferencia.2onde
2.4.2 MATRIZ IN"ERSA
2onde%
ASI : Ecuc!"# c$ct%$&'t!c
2.4.3 ECUACION CARACTERISTICA
#e obtiene a partir de la matriz de transferencia y se obtiene al igualar ladeterminante del denominador de la matriz de la transferencia ,&atriz inversa0igual a cero" obtenindose las races de la &atriz de transferencia conocidocomo polos.
17
=
1..........00
.
.
.
0..........10
0..........01
I)3(.....)( )()(
1
)()( ssss UBASIX =
)()()()(
1
)()()( ])[( sssssss UDUBASIC += [ ] )()()(1)()()( )( ssssss UDBASIC +=
ciaTransferende!atri"DBASICU
# ssss
s
s.....)( )()(
1
)()(
)(
)( +==
( )ASI
ASId$ASI
T
= 1)(
( ) 0. == AsIASIDet
-
5/22/2018 Contro I - Cap. II-A
18/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
E?9-+
2eterminar la matriz de transferencia del sgt.#istema expresado en el espacio de estado"mediante las ecuaciones de estado y de salidadefinida por%
S+c):
S,/9aMIMO: 2)/a*a > 2,a*a,.
18
+=
2
1
2
1o
01
10
X
53
01
X
u
u
X
A
[ ]0,10
01
2
1
)( =
= D
X
Xs
-
5/22/2018 Contro I - Cap. II-A
19/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
( )ASI = S
10
01 -
53
01=
+
+53
01
S
S
( ) ( )( ) ( )51
1S3
05
1
++
+
+
==
SS
S
ASIASIad$ASI
T
Ma/ * /a),@)ca
PROBLEMA
allar C8, t 0 y C9, t 0 si
on las condiciones iniciales.
19( ) ( )
[ ]051
01
10
13
05
10
01
)(
)( +
++
+
+
=SS
S
S
U
s
s [ ]0)5(1)(S
01
10
1S3
05
)(
)( +++
+
+
=S
S
U
s
s
( ) ( )51
31
50
)(
)(
++
+ +=SS
S
S
U
s
s
=
11
)0(2
)0(1
XX
=
2
1
2
01
0
23
10
X
X
X
X
-
5/22/2018 Contro I - Cap. II-A
20/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
SOLUCIONES .
ORDENANDO Y ARRELANDO
20
=23
10A
( ) 3223
1
23
10
10
01 ++=
+
=
= SS
S
SSASI
( ) ( ) ( )22 2132 ++=++= SSSASI
( )[ ] ( )
== ASI
ASIad$lASIle
T
tA 111
( )
+
++=
3
12
32
11
S
S
SSle tA
( )
( ) ( ) ( ) ( )
( ) ( )( )
( ) ( ) 2222
22221
21
11
21
321
1
21
11
++
+
++
++++
++
=
S
S
S
SS
S
le tA
2e2cos22
3
2212
212cos
t- tsentetsene
tsenetsenee
ett
ttt
At
+=
1
1)0()(
==
AtAt
t eXeX
( ) ( ) ( ) ( ) ( ) ( )
( ) ( )( )
( ) ( ) ( ) ( )
21
2
2
1
21
1
21
2
2
3
21
2
2
1
21
2
2
1
21
1
222222
2222221
++
++
+
++
++
+++
++
+
=
Sx
S
S
S
SSS
S
le tA
1
1
22
12cos2
2
3
22
12
2
12cos
)(
+=
tsenetetsene
tsenetsenee
Xttt
ttt
t
=
tsenete
teX
tT
t
t222cos
2cos)(
-
5/22/2018 Contro I - Cap. II-A
21/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
PROBLEMA
2ado un sistema &&/ con las siguientes ecuaciones de entrada y salida.
SE PIDE:
la matriz de transferencia.
SOLUCION
2onde%
( )
+
=
=
32
1
32
10
10
01
S
SSASI
( ) ( )( ) ( )( )12
2
13
23
2
13
23
2
13
2 ++
+
=++
+
=++
+
=
=SS
S
S
SS
S
S
SS
S
S
ASI
ASIad$ASI
T
( ) ( ) ( )0
10
01
12
2
13
10
01
)(
)( +
++
+
==
SS
S
S
%
#
s
s
S
( ) ( )( )( )( ) ( )( )
( )( ) ( ) ( )
++++
+++++
=++
+
=
1212
232
1
12
3
12
2
13
SS
S
SS
SSSS
S
SS
S
S
# s
PROBLEMA
21
+
=
2
1
2
1
2
01
0
10
01
32
10
U
U
X
X
X
X
=
2
1
2
1
10
01
X
X
( ) ( ) DBASIC%
#
s
s
s +== 1
)(
)(
-
5/22/2018 Contro I - Cap. II-A
22/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
Da*+ ,,/9a:
#e pide" la matriz de transferencia.
S+c) .
2el sistema la salida )8,#0
[ ] (1)21 )(2)(1)(1)(2)(1 ssssS US
++=
La sa!"a #2(5).
[ ] (2)1 )(2)(2)(2 sss US
=
$e (1).
)(1)(1)(2)(1 ssss US +=
( ) (3)1)(2)(1)(1 sss
US
=+
$e (2).
)(2)(2)(2 sss US =
( ) )(2)(21 ss US =+
( )(4).....
1
)(2
)(2 +=
S
U
s
s
22
-
5/22/2018 Contro I - Cap. II-A
23/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
(4) e% (3).
( ) ( )11)(2
)(1)(1 +=+ SU
US sss
( ) ( )(5).....
112
)(21(s)
)(1 +
+=
S
U
S
U
s
s
/rdenando la ecuacin ,@0
( )(6).....
1
10 )(2)(1)(2 sSs US
U+
+=
2e las ecuaciones ,?0 y ,A0
( )
1
1
0
1
1
1
12
)(2
)(1
+
+
+=
S
SS
s
s
2onde la matriz de transferencia ser%
2.4.4 REPRESENTACION EN EL ESPACIO DE ESTADO EN FORMACANONICA DE LOS SISTEMAS BASADOS EN LA FUNCION DE
TRANSFERENCIA
23
( )
+
+
+
1
10
1
1
1
12
S
SS
-
5/22/2018 Contro I - Cap. II-A
24/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
#e presentan dos casos%
a. CUANDO LA FUNCI&N E%CITADORA NO INCLUYE TRMINOSDERI"ATI"OS:
#i se tiene un sistema de orden )donde%
Esta ecuacin se puede convertir en )ecuaciones diferenciales de primerorden" para ello se tiene que elegir )variables con la siguiente asignacin&
$a
ecuacin ,80 ser%
24
( ))1(...................... )()(
)(
11
)(
1
1 ttn
t
nn
t
n
u
t
n
uyadt
dya
dt
yda
dt
yd=++++
+
(2).......
..........
.
.
.
.
)()(1)(21)(1)(1
)(
1
)(
)(4)(3(t))(3
)(3)(2(t))(2
)(2)(1)()(1
+==
==
==
==
ttntntntn
o
n
t
t
n
tn
Tt
ooo
t
tt
oo
t
tt
o
tt
uXaXaXaX
d
ydX
XXyX
XXyX
XXyX
)()(1)(21)(1)( .......... ttntntntno
uXaXaXaX =++++
-
5/22/2018 Contro I - Cap. II-A
25/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
$as ecuaciones ,90 #e pueden expresar en forma matricial.
#i consideramos que la
salida del sistema > esla variable %1 " entonces
dicha salida se puedeescribir.
En su forma compactaseria. ,Ecuacin de salidadel sistema0
25
u
X
X
X
aaaaX
X
X
nnnntn
o
to
t
o
1
.
.
.
0
0
.
.
.
...........-
.
.
.
0...........100
0...........010
...
2
1
121)(
)(2
)(1
+
=
[ ] [ ] u
X
X
X
y
n
0
.
.
.0..........01
2
1
+
=
-
5/22/2018 Contro I - Cap. II-A
26/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
$o%"e
2iagrama del bloque de ecuaciones de estado y salida son%
Tambin se puede mostrar el diagrama de bloque completo del sistema de
orden )
2onde%
26
= 0.........01C
=
1
0
0
B
[ ]0=D
=
121 .........
.
.
.
0..........100
0..........010
aaaa
A
nnn
= 1..............00TB
-
5/22/2018 Contro I - Cap. II-A
27/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
#iendo%
( )
o
tnX= -
naX ( )t1 - 1na
X ( )t2 - 2na X ( )t
'''.- 2a 1nX ( )t - ( ) ( )ttnXa +1
/bservacin%
Fn modelo de estado de un sistema lineal con cuatro matrices definido de
5rimer orden de una ecuacin. 2iferencial de grado )esta definido por%
b. CUANDO SU FUNCIONE%CITADORA INCLUYE TERMINOS DERI"ATI"OS
#i se tiene un sistema en orden. )
( ) ( )
( )( )
( ) ( ) ( )
( )tn
t
nn
t
n
n
t
n
otn
t
t
nn
n
n
t
n
ubdt
udb
dt
udb
dt
udbya
d
dya
dt
yda
dt
yd++++=++++
11
1
111
1
1 ................
En laplace esta ecuacin puede escribirse%
C+) ,a*a:
27
[ ] )(1)(1221 tt yXyyXXX ===== )(
1
2)(111
332 tt yXyyXXX =====
[ ] )()()()( tntntntnn Xdt
dXXXX ====
( )
( )0)(
)()()()(
)()()()(
0XX
uDXCy
uBXAX
t
ttttt
tttt
o
t =
+=
+=
ububububyayayay n
o
n
nn
on
o
n
nn
++++=++++
+
1
1
1
1
1 .........,.....
SalidalaEsy
EntradalaEsu
Donde
&
&
&
( )
( ))(........
........
.....
1
1
1
1
1
1 nn
nn
nn
nn
o
s
s
aSaSaS
bSbSbSb
U
++++++++
=
[ ] 2)(1)(
21
)(1
===== ntntnntnnn yXyyXXX
-
5/22/2018 Contro I - Cap. II-A
28/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
2onde las variables de estado debido a los trminos derivados del segundo
miembro" deben ser tal que eliminan la derivada de en la ecuacin deestado%Fna forma de obtener una ecuacin de estado y una ecuacin de salida es
definir las siguientes )variables como un con(unto de ) variables deestado.
D+)*:
P+ + /a)/+:
(1 en (2
(3en (4.
A*9
-
5/22/2018 Contro I - Cap. II-A
29/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
$a representacin en el espacio estado para la ecuacin de entrada y desalida es%
29
)(.....
.
.
.
.
.
.
..........
1..........0000000
.
.
.0..........0000100
0..........0000010
.
.
.
1
1
0
1
2
1
1654321
2
1
Iu
X
X
X
X
aaaaaaaaX
X
X
n
n
n
n
nnnnnnnn
o
o
o
+
=
[ ]
.
.
.00........0001
2
1
=
nX
X
X
y
ticaCaracterisEcuacAsI .0=
-
5/22/2018 Contro I - Cap. II-A
30/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
PROBLEMA
2ado el sistema encontrar una representacin en variable de estado.
6dems se cumple%
SOLUCI&N:
62E&G#%
2el dato
(3)
30
UXX = 133
21 1X
SX =
21 XXS =
(1)X00 32110
21
0
++== XXXXX
( ) (2)3
2132 XX
SX
I
+=
UXX +=3
1
3UXX = 133
-
5/22/2018 Contro I - Cap. II-A
31/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
: en 9
2el dato%
UXX = 133
2el 2iagrama de 'loque
/rdenando%
$a salida%
31
( )132
3
2XUX
SX +
+=
( ) 132 2223 XUXSX +=+ 1322 2223 XUXXSX +=+ UXXXSX 2232 3212 ++=
)(42232 32120
UXXXX ++=
UUS
SX 5
23 +
+=
USUSSX )5()2()5(3 ++=+ USUUSUXSX 525 33 +=+
UXSX 35 33 =
(5)3500 321
0
3 UXXXX +=
U
X
X
X
X
X
X
3
2
0
5-00
23-2-
010
3
2
1
2
0
2
0
1
0
+=
[ ]
3
2
1
001
X
X
X
y=1X=
-
5/22/2018 Contro I - Cap. II-A
32/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
PROBLEMA
2ada la 3.T.
)545(
43022
2
)(
)(
SS
S
%
S
S
=
allar%
a. )()()( ttt BuAxX +=
)()( tt Cx =
N+/a.6para la solucin use el modelo de grafica de flu(os de se*ales conalimentacin directa de las variables de estado a la se*al de salida.
SOLUCION
2onde%
24
2
22
2
)(
)(
455
430
)545(
430
SS
S
SS
S
%
S
S
+
=
=
32
-
5/22/2018 Contro I - Cap. II-A
33/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
004505
300400234
234
)(
)(
+++++++=SSSS
SSSS
%
S
S
300400004505 234234 +++=++++ uuuuyyyy
1a 2a 3a 4a 0b 1b 2b 3b 4b
6dems% Deemplazando datos%
210)4(4530
0)4(00
4
0
0
4
3
22
1
0
====
====
b
6dems%
00
0
1
0
uyX =
2onde% uXX 110
2 =
uXX 1210
+=
#imilar para% uXX 2320
+=
uXX 3430
+=
uXaXaXaXaX 44132231440
+++=Deemplazando datos%
33
0413223144
03122133
021122
0111
00
aaaab
aaab
aab
ab
b
==
==
=
uuyX
uyX
1
0
0
0
2
01
=
=
uXXXXX
uXXXXuXX
uXXXXuXX
uXXXXuXX
21004500
00000
40004
00000
43214
0
432143
0
432132
0
432121
0
+=
++++=+=
+++==
++++=+=
-
5/22/2018 Contro I - Cap. II-A
34/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
$a salida y 7HH.2.8 LINEALIZACI&N USANDO
EL DESARROLLO DE LASERIE DE TAYLOR
5ara aplicar este mtodo tambinse debe suponer que las variablesse desvan poco dentro de unacondicin de operacin. #i setiene un sistema no lineal cuya
entrada es % ( )t y cuya salida >( )t " donde su relacin esta dado
por.
34
u
X
X
X
X
X
X
X
X
210
0
4
0
04500
1000
0100
0010
4
3
2
1
4
0
3
0
2
0
1
0
+
=
-
5/22/2018 Contro I - Cap. II-A
35/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
En la linealizacin se considera solo una parte de la curva como una lnearecta" que es pendiente a esa curva en el puntoX s ,que va a ser la condicinnormal de la operacin en el punto de equilibrioX s , Y s .X s : Es el valor de la variable en el estado estacionario.
6plicando el desarrollo de la serie de Taylor y suponiendo que las variables sedesvan poco dentro de una ecuacin de operacin.
( )xfy=
$a condicin normal de operacin se da el punto de equilibrio" SS X ,
2onde%
#i la evaluacin ( )sXX es peque*a se podra despreciar los trminos deorden superior en ( )SXX
#iendo
35
( ) ( )S
xxXXenevaluadosSon
dx
Fd
dx
dF=,
2
2
( ) ( )sX
S XXdx
fdyy +=
( ) ( ) (2).........sx
s XXdx
dfyy =
( ) ( )( ) ( ) ( ) ( ) ......
12
1 22
2
+
+
+== S
X
x
o
S
X
x
Xx XXdX
fdXX
dX
dfffy
SS
S
( )sS xfy =
( )xfy=
-
5/22/2018 Contro I - Cap. II-A
36/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
4emos que ( )syy es proporcional a ( )sXX .5or lo tanto la ecuacin (2 da un modelo matemtico lineal de un sistema nolineal dado.
LA SERIE DE TAYLOR APLICADOS A SISTEMA NO LINEAL CUYA SALIDAES UNA FUNCI&N DE DOS ENTRADAS 21 XyX
#u aplicacin lineal ser alrededor del punto normal de operacin6plicando Taylor%
2onde%
$as derivadas parciales se calculan%
sS XXXX 2211 , == " adems cerca del punto normal de operacin se puededespreciar los trminos de orden superior.El modelo matemtico lineal aproximado ser%
6l eliminarlostrminosde orden
superior se tiene.
36
( )21,XXfy=
( )SS
XX 21 ...
( ) ( ) +
+
+= SXX
S
XX
xx XXx
fXX
x
ffy
SssS
ss 22
2
11
1
),(
2,12,1
21
( ) ( )( ) ( ) ] ...22
1 2222
2
2
2211
21
22
112
1
2
212,12,1
+
+
+
+ SXX
ss
XX
s
XX
XXx
fXXXX
xx
fXX
x
f
SSSSSS
( ) ( )
( ) ( )( ) ( )
+
+
+
+
+
=
2
222
2
2
2211
21
22
112
1
2
22
2
11
2
2
1sssS
sss
XXx
fXXXX
xx
fXX
x
f
XXx
fXX
X
fyy
( ) ( ) )3(...........222
11
1
sss XXxfXX
xfyy +=
-
5/22/2018 Contro I - Cap. II-A
37/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
$a ecuacin (3 da un modelo matemtico lineal de un sistema no lineal dado.
2. "ARIABLE DE DES"IACION
onocida tambin como variable de perturbacin el cual indica cuanto sedesvan las variables reales del punto en estado estacionario.5ara analizar los modelamientos de un sistema de control es importanterealizar el cambio de variables reales a variables desviacin donde%
1X % 4ariable de
desviacin
( )tX % 4ariable real en el instante ItJ
C s % variable de estado estacionario
EJEMPLO
Fna reaccin qumica que se inicia a Co45 en un instante t o " que es la
temperatura sT en el estado estacionario.2urante la reaccin qumica esta temperatura no se mantiene constante sinova a cambiar debido a las perturbaciones as%
! En 1t se encontrara a Co44 ! En 2t se encontrara a Co46
4emos que las variables reales de la temperatura en los tiempos 210 ,, ttt sonCCC ooo 46,44,45 respectivamente.
En trminos de variable de desviacin con respecto al valor deseado que es@? Co para%
! En ot no hay desviacin
! En 1t hay una desviacin de ( )( )CCXXXC ooSto 45441 1 ==! En 2t hay una desviacin de ( )( )CCXXXC ooSto 45461 1 ==
37
( ) st XXX =1
-
5/22/2018 Contro I - Cap. II-A
38/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
PROBLEMA
Fn dispositivo no lineal se representa con la funcin%
21)( X'
X ==
#e pide linealizar el dispositivo" cuando el punto de operacin para la entrada
X es 21
0=X .
6sumir que la entrada varia en un peque*o intervalo alrededor de 0X .
SOLUCI&N:
S!
21
)( X'
X ==
El dispositivoesta dado%
38
-
5/22/2018 Contro I - Cap. II-A
39/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
5ara la linealizacin del dispositivo no lineal se considera solo una parte de lacurva como una lnea recta" pendiente a esa curva en el punto 0X " que es lacondicin normal de operacin en el punto de equilibrio ),( SS X .6plicando la serie de Taylor y suponiendo que las variables se desvan un pocodentro de una ecuacin de operacin.
)(X'=
2onde los trminos derivativos2
)(2
)(,
dx
'd
dx
d' XX " etc. son evaluados en 0XX= .
omo la evaluacin , 0XX 0 es peque*a se puede despreciar lostrminos de orden superior en , 0XX 0.
Kuedara% Xy=
)X-X( 0)(
0)(dx
d'y'y X
X +==
2onde%2
10)(0 === X'y X
).( 0)(
0
0
XXdx
d'yy
XX
X +==
39
-
5/22/2018 Contro I - Cap. II-A
40/40
CURSO : CONTROL I DOCENTE : ING. ROBERTO JAIME QUIROZ SOSA
)(1
2
1
2
2( 0
21
)00
00
XXX
xXXdx
Xdyy
XXXX
+=+====
)2
1(
21
1
2
1
2
1+= XxxX
2
2)
21(
2
2)( +== Xy' X
)2
11(2
2)( += X' X
2ispositivo $inealizado
40
4
2
2
2)( += X' X