Articulo de referencia

Linealización por retroalimentación

Diagrama de bloques que ilustra la linealización por retroalimentación de un sistema no lineal. La linealización por retroalimentación es una estrategia común empleada en el con...

Diagrama de bloques que ilustra la linealización por retroalimentación de un sistema no lineal.

La linealización por retroalimentación es una estrategia común empleada en el control no lineal para controlar sistemas no lineales . Las técnicas de linealización por retroalimentación se pueden aplicar a sistemas de control no lineales de la forma

dóndeincógnita(t)Rnorte{\displaystyle x(t)\in \mathbb {R} ^{n}}es el estado,1(t),,metro(t)R{\displaystyle u_{1}(t),\ldots ,u_{m}(t)\in \mathbb {R} }son las entradas. El enfoque implica transformar un sistema de control no lineal en un sistema de control lineal equivalente a través de un cambio de variables y una entrada de control adecuada. En particular, se busca un cambio de coordenadasz=Φ(incógnita){\displaystyle z=\Phi (x)}y entrada de control=a(incógnita)+b(incógnita)v,{\displaystyle u=a(x)+b(x)\,v,}para que la dinámica deincógnita(t){\displaystyle x(t)}en las coordenadasz(t){\displaystyle z(t)} adoptan la forma de un sistema de control lineal y controlable,

Posteriormente, se puede aplicar una estrategia de control de bucle externo al sistema de control lineal resultante para lograr el objetivo de control.

Linealización por retroalimentación de sistemas SISO

Aquí, consideremos el caso de linealización por retroalimentación de un sistema de entrada única y salida única (SISO). Resultados similares pueden extenderse a sistemas de entrada múltiple y salida múltiple (MIMO). En este caso,R{\displaystyle u\in \mathbb {R} }yyR{\displaystyle y\in \mathbb {R} }El objetivo es encontrar una transformación de coordenadas.z=T(incógnita){\displaystyle z=T(x)}que transforma el sistema (1) en la llamada forma normal que revelará una ley de retroalimentación de la forma

que generará un mapa lineal de entrada - salida a partir de la nueva entrada.vR{\displaystyle v\in \mathbb {R} }a la saliday{\displaystyle y}Para garantizar que el sistema transformado sea una representación equivalente del sistema original, la transformación debe ser un difeomorfismo . Es decir, la transformación no solo debe ser invertible (es decir, biyectiva), sino que tanto la transformación como su inversa deben ser suaves , de modo que la diferenciabilidad en el sistema de coordenadas original se conserve en el nuevo sistema de coordenadas. En la práctica, la transformación solo puede ser localmente difeomorfa y los resultados de la linealización solo son válidos en esta región más pequeña.

Se requieren varias herramientas para resolver este problema.

Derivada de Lie

El objetivo de la linealización por retroalimentación es producir un sistema transformado cuyos estados sean la salida.y{\displaystyle y}y su primera(norte1){\displaystyle (n-1)}derivadas. Para comprender la estructura de este sistema objetivo, utilizamos la derivada de Lie . Consideremos la derivada temporal de (2), que se puede calcular utilizando la regla de la cadena .

y˙=dh(incógnita)dt=h(incógnita)incógnitaincógnita˙=h(incógnita)incógnitaF(incógnita)+h(incógnita)incógnitagramo(incógnita){\displaystyle {\begin{aligned}{\dot {y}}={\frac {{\mathord {\operatorname {d} }}h(x)}{{\mathord {\operatorname {d} }}t}}&={\frac {\partial h(x)}{\partial x}}{\dot {x}}\\&={\frac {\partial h(x)}{\partial x}}f(x)+{\frac {\partial h(x)}{\partial x}}g(x)u\end{aligned}}}

Ahora podemos definir la derivada de Lie deh(incógnita){\displaystyle h(x)}a lo largo deF(incógnita){\displaystyle f(x)}como,

LFh(incógnita)h(incógnita)incógnitaF(incógnita),{\displaystyle L_{f}h(x)\triangleq {\frac {\partial h(x)}{\partial x}}f(x),}

y de manera similar, la derivada de Lie deh(incógnita){\displaystyle h(x)}a lo largo degramo(incógnita){\displaystyle g(x)}como,

Lgramoh(incógnita)h(incógnita)incógnitagramo(incógnita).{\displaystyle L_{g}h(x)\triangleq {\frac {\partial h(x)}{\partial x}}g(x).}

Con esta nueva notación, podemos expresary˙{\displaystyle {\dot {y}}}como,

y˙=LFh(incógnita)+Lgramoh(incógnita){\displaystyle {\dot {y}}=L_{f}h(x)+L_{g}h(x)u}

Nótese que la notación de derivadas de Lie es conveniente cuando tomamos derivadas múltiples con respecto al mismo campo vectorial o a uno diferente. Por ejemplo,

LF2h(incógnita)=LFLFh(incógnita)=(LFh(incógnita))incógnitaF(incógnita),{\displaystyle L_{f}^{2}h(x)=L_{f}L_{f}h(x)={\frac {\partial (L_{f}h(x))}{\partial x}}f(x),}

y

LgramoLFh(incógnita)=(LFh(incógnita))incógnitagramo(incógnita).{\displaystyle L_{g}L_{f}h(x)={\frac {\partial (L_{f}h(x))}{\partial x}}g(x).}

Grado relativo

En nuestro sistema linealizado de retroalimentación compuesto por un vector de estado de la saliday{\displaystyle y}y su primera(norte1){\displaystyle (n-1)}derivados, debemos entender cómo la entrada{\displaystyle u}entra en el sistema. Para ello, introducimos la noción de grado relativo. Se dice que nuestro sistema dado por (1) y (2) tiene grado relativo.rW{\displaystyle r\in \mathbb {W} }en un puntoincógnita0{\displaystyle x_{0}}si,

LgramoLFkh(incógnita)=0incógnita{\displaystyle L_{g}L_{f}^{k}h(x)=0\qquad \forall x}en un barrio deincógnita0{\displaystyle x_{0}}y todokr2{\displaystyle k\leq r-2}
LgramoLFr1h(incógnita0)0{\displaystyle L_{g}L_{f}^{r-1}h(x_{0})\neq 0}

Considerando esta definición de grado relativo a la luz de la expresión de la derivada temporal de la saliday{\displaystyle y}Podemos considerar que el grado relativo de nuestro sistema (1) y (2) es el número de veces que tenemos que diferenciar la salida.y{\displaystyle y}antes de la entrada{\displaystyle u}aparece explícitamente. En un sistema lineal invariante en el tiempo (LTI) , el grado relativo es la diferencia entre el grado del polinomio denominador de la función de transferencia (es decir, el número de polos ) y el grado de su polinomio numerador (es decir, el número de ceros ).

Linealización por retroalimentación

Para la discusión que sigue, asumiremos que el grado relativo del sistema esnorte{\displaystyle n}En este caso, después de diferenciar la salidanorte{\displaystyle n}tiempos que tenemos,

y=h(incógnita)y˙=LFh(incógnita)y¨=LF2h(incógnita)y(norte1)=LFnorte1h(incógnita)y(norte)=LFnorteh(incógnita)+LgramoLFnorte1h(incógnita){\displaystyle {\begin{aligned}y&=h(x)\\{\dot {y}}&=L_{f}h(x)\\{\ddot {y}}&=L_{f}^{2}h(x)\\&\vdots \\y^{(n-1)}&=L_{f}^{n-1}h(x)\\y^{(n)}&=L_{f}^{n}h(x)+L_{g}L_{f}^{n-1}h(x)u\end{aligned}}}

donde la notacióny(norte){\displaystyle y^{(n)}}indica elnorte{\displaystyle n}derivada dey{\displaystyle y}. Porque asumimos que el grado relativo del sistema esnorte{\displaystyle n}, los derivados de Lie de la formaLgramoLFih(incógnita){\displaystyle L_{g}L_{f}^{i}h(x)}parai=1,,norte2{\displaystyle i=1,\dots ,n-2}son todos cero. Es decir, la entrada{\displaystyle u}no tiene ninguna contribución directa a ninguno de los primeros(norte1){\displaystyle (n-1)}derivadas de th.

La transformación de coordenadasT(incógnita){\displaystyle T(x)}que pone el sistema en forma normal proviene del primero(norte1){\displaystyle (n-1)}derivados. En particular,

z=T(incógnita)=[z1(incógnita)z2(incógnita)znorte(incógnita)]=[yy˙y(norte1)]=[h(incógnita)LFh(incógnita)LFnorte1h(incógnita)]{\displaystyle z=T(x)={\begin{bmatrix}z_{1}(x)\\z_{2}(x)\\\vdots \\z_{n}(x)\end{bmatrix}}={\begin{bmatrix}y\\{\dot {y}}\\\vdots \\y^{(n-1)}\end{bmatrix}}={\begin{bmatrix}h(x)\\L_{f}h(x)\\\vdots \\L_{f}^{n-1}h(x)\end{bmatrix}}}

transforma trayectorias desde la originalincógnita{\displaystyle x}sistema de coordenadas en el nuevoz{\displaystyle z}sistema de coordenadas. Siempre que esta transformación sea un difeomorfismo , las trayectorias suaves en el sistema de coordenadas original tendrán contrapartes únicas en el sistema de coordenadas.z{\displaystyle z}sistema de coordenadas que también son suaves. Esosz{\displaystyle z}Las trayectorias serán descritas por el nuevo sistema,

{z˙1=LFh(incógnita)=z2(incógnita)z˙2=LF2h(incógnita)=z3(incógnita)z˙norte=LFnorteh(incógnita)+LgramoLFnorte1h(incógnita).{\displaystyle {\begin{cases}{\dot {z}}_{1}&=L_{f}h(x)=z_{2}(x)\\{\dot {z}}_{2}&=L_{f}^{2}h(x)=z_{3}(x)\\&\vdots \\{\dot {z}}_{n}&=L_{f}^{n}h(x)+L_{g}L_{f}^{n-1}h(x)u\end{cases}}.}

Por lo tanto, la ley de control de retroalimentación

=1LgramoLFnorte1h(incógnita)(LFnorteh(incógnita)+v){\displaystyle u={\frac {1}{L_{g}L_{f}^{n-1}h(x)}}(-L_{f}^{n}h(x)+v)}

genera un mapa de entrada - salida lineal a partir dev{\displaystyle v}az1=y{\displaystyle z_{1}=y}El sistema linealizado resultante

{z˙1=z2z˙2=z3z˙norte=v{\displaystyle {\begin{cases}{\dot {z}}_{1}&=z_{2}\\{\dot {z}}_{2}&=z_{3}\\&\vdots \\{\dot {z}}_{n}&=v\end{cases}}}

es una cascada denorte{\displaystyle n}integradores y un control de bucle externov{\displaystyle v}puede elegirse utilizando la metodología estándar de sistemas lineales. En particular, una ley de control de realimentación de estado de

v=Kz,{\displaystyle v=-Kz\qquad ,}

donde el vector de estadoz{\displaystyle z}es la saliday{\displaystyle y}y su primera(norte1){\displaystyle (n-1)}derivados, resultados en el sistema LTI

z˙=Az{\displaystyle {\dot {z}}=Az}

con,

A=[010000100001k1k2k3knorte].{\displaystyle A={\begin{bmatrix}0&1&0&\ldots &0\\0&0&1&\ldots &0\\\vdots &\vdots &\vdots &\ddots &\vdots \\0&0&0&\ldots &1\\-k_{1}&-k_{2}&-k_{3}&\ldots &-k_{n}\end{bmatrix}}.}

Así pues, con la elección adecuada deK{\displaystyle K}Podemos colocar arbitrariamente los polos del lazo cerrado del sistema linealizado.

Dinámica cero inestable

La linealización por retroalimentación se puede lograr con sistemas que tienen un grado relativo menor quenorte{\displaystyle n}Sin embargo, la forma normal del sistema incluirá dinámicas nulas (es decir, estados que no son observables a partir de la salida del sistema) que pueden ser inestables. En la práctica, las dinámicas inestables pueden tener efectos perjudiciales en el sistema (por ejemplo, puede ser peligroso que los estados internos del sistema crezcan sin límite). Estos estados no observables pueden ser controlables o, al menos, estables, por lo que se pueden tomar medidas para asegurar que no causen problemas en la práctica. Los sistemas de fase mínima proporcionan información sobre las dinámicas nulas.

Linealización por retroalimentación de sistemas MIMO

Aunque NDI no se limita necesariamente a este tipo de sistema, consideremos un sistema MIMO no lineal que es afín en la entrada.{\displaystyle \mathbf {\mathbf {u} } }, como se muestra a continuación.

Se supone que la cantidad de entradas es la misma que la cantidad de salidas. Digamos que haymetro{\displaystyle m}entradas y salidas. LuegoGRAMO=[gramo1gramo2gramometro]{\displaystyle G=[\mathbf {g} _{1}\,\mathbf {g} _{2}\,\cdots \,\mathbf {g} _{m}]}es unnorte×metro{\displaystyle n\times m}matriz, dondegramoj{\displaystyle \mathbf {g} _{j}}son los vectores que componen sus columnas. Además,Rmetro{\displaystyle \mathbf {u} \in \mathbb {R} ^{m}}yyRmetro{\displaystyle \mathbf {y} \in \mathbb {R} ^{m}}Para utilizar una derivación similar a la de SISO, el sistema de la ecuación 4 se puede dividir aislando cada uno de los componentes.i{\displaystyle i}salidayi{\displaystyle y_{i}}, como se muestra en la ecuación 5.

De manera similar a SISO, se puede demostrar que hasta el(ri1){\displaystyle (r_{i}-1)}'º derivado deyi{\displaystyle y_{i}}, el términoLgramojhi(incógnita)=0{\displaystyle L_{g_{j}}h_{i}(\mathbf {x} )=0}. Aquíri{\displaystyle r_{i}}se refiere al grado relativo de lai{\displaystyle i}'th salida. Análogamente, esto da

Si se resuelve de la misma manera que con SISO, se descubre que definir una entrada virtualvi{\displaystyle v_{i}}de tal manera que

linealiza estoi{\displaystyle i}sistema. Sin embargo, simetro>1{\displaystyle m>1},{\displaystyle \mathbf {u} }Obviamente no se puede resolver dado un valor paravi{\displaystyle v_{i}}Sin embargo, establecer tal ecuación para todosmetro{\displaystyle m}resultados,y1,y2,,ymetro{\displaystyle y_{1},y_{2},\ldots ,y_{m}}, resulta enmetro{\displaystyle m}ecuaciones de la forma mostrada en la ecuación 7. La combinación de estas ecuaciones da como resultado una ecuación matricial, que generalmente permite resolver para la entrada.{\displaystyle \mathbf {u} }, como se muestra a continuación.

Véase también

Lecturas adicionales

  • A. Isidori, Sistemas de control no lineales, tercera edición, Springer Verlag, Londres, 1995.
  • HK Khalil, Sistemas no lineales, tercera edición, Prentice Hall, Upper Saddle River, Nueva Jersey, 2002.
  • M. Vidyasagar, Análisis de sistemas no lineales , segunda edición, Prentice Hall, Englewood Cliffs, Nueva Jersey, 1993.
  • B. Friedland, Diseño avanzado de sistemas de control , edición facsímil, Prentice Hall, Upper Saddle River, Nueva Jersey, 1996.

Referencias

  1. Isidori, Alberto (1995). Sistemas de control no lineales (Tercera  ed.). Springer-Verlag Londres. pag.  5.ISBN 978-1-4471-3909-6.
  2. ^ H. Nijmeijer y A. van der Shaft, Sistemas de control dinámico no lineal, Springer-Verlag, p. 163, 2016.
  3. Isidori, Alberto (1995). Sistemas de control no lineales (Tercera ed.). Springer-Verlag Londres. pag. 147.ISBN   978-1-4471-3909-6.
  • Fabio Celani y Alberto Isidori, eds. (2009). "Linealización por retroalimentación" . Scholarpedia . Recuperado el 31 de diciembre de 2022 .
  • ECE  758: Modelado y control no lineal de un manipulador de articulación flexible de un solo eslabón : ofrece una explicación y una aplicación de la linealización por retroalimentación. 
  • ECE  758: Ejemplo de linealización de bola en tubo : aplicación trivial de la linealización para un sistema que ya está en forma normal (es decir, no es necesaria ninguna transformación de coordenadas). 
  • El lenguaje Wolfram incluye funciones para realizar la linealización por retroalimentación , calcular órdenes relativos y determinar la dinámica cero .