Óptimo no significa el mejor control para cualquier situación: significa el que minimiza un coste y cumple unas condiciones concretas. Robusto significa que una garantía se mantiene para toda una familia declarada de modelos. Vamos a construir primero un coste, deducir una ley LQR pequeña y después separar lo que prueba para el nominal de lo que exige una incertidumbre.
Si necesitas repasar:
1. Antes de comparar controles fijamos qué problema resuelven
Un problema de control incluye dinámica, dato inicial, horizonte, coste y restricciones. Minimizar tiempo, esfuerzo o error puede producir decisiones diferentes. Un horizonte finito termina en T y puede penalizar o fijar el estado terminal. En horizonte infinito hay que revisar convergencia del coste y la clase de políticas admisibles.
Esta introducción deducirá LQR de regulación al origen con modelo lineal conocido, estado disponible y sin saturación explícita. Los primeros ejemplos finitos muestran por qué cambiar un dato del problema cambia el significado de óptimo. No se comparan números de costes diferentes como si midieran lo mismo.
Llegar al mismo lugar no implica gastar el mismo esfuerzo
x′=u, x(0)=2, x(1)=0; coste ∫₀¹u²dt.
Control uniforme
u=−2 durante una unidad de tiempo llega a cero y cuesta cuatro.
Control concentrado
u=−4 durante la primera mitad y cero después también llega, pero cuesta 16·1/2=8.
Certificar
∫u=−2.Por Cauchy–Schwarz, (∫u)²≤∫u²·1; todo candidato cuesta al menos cuatro, que alcanza el uniforme.
El objetivo y el límite terminal permiten probar cuál es mejor.
El área de u es la misma en ambos, pero la integral de su cuadrado es distinta.
Un coste terminal blando produce otra decisión
x′=u, x(0)=1, horizonte uno; J=∫₀¹u²dt+x(1)².
Resumir el efecto
Sea v=∫₀¹u.Entonces x(1)=1+v y ∫u²≥v².Minimizar
J≥v²+(1+v)²=2(v+1/2)²+1/2.Alcanzar
u=−1/2 constante da v=−1/2, x(1)=1/2 y J=1/2.No obliga a terminar exactamente en cero.
Penalizar el estado final no equivale a imponerlo como igualdad.
No se usa una solución estacionaria de horizonte infinito para sustituir este problema finito.
2. El coste cuadrático separa tamaño del estado y esfuerzo de entrada
Q es simétrica semidefinida positiva y R simétrica definida positiva. El primer término puede dejar direcciones sin penalizar; el segundo hace costosa cualquier entrada no nula y permite invertir R. Las variables y pesos deben conservar unidades o una normalización indicada.
Cambiar pesos cambia el problema. Multiplicar Q y R por la misma constante positiva multiplica todo J y no cambia su minimizador. Penalizar más el control puede favorecer otra política; el efecto concreto se comprueba en el modelo, no se interpreta el peso como un límite duro de fuerza.
Leer una forma cuadrática con valores concretos
x=(1,2)ᵀ, Q=diag(4,1), u=3 y R=2.
Estado
xᵀQx=4·1²+1·2²=8.Entrada
uᵀRu=2·3²=18.Instante
El coste instantáneo es veintiséis. Integrarlo requiere la trayectoria de x y u, no solo este valor aislado.
Los dos términos expresan penalizaciones elegidas y compatibles.
Con Q=diag(1,0), el estado (0,5) tendría coste de estado cero: semidefinida no significa estrictamente positiva.
Comparar políticas dentro de un mismo coste
Integrador x′=u, x(0)=1; política u=−kx con k>0, Q=4.
Trayectoria
x=exp(−kt), así J(k)=(4+Rk²)/(2k).Con R=1
k=1 cuesta 5/2; k=2 cuesta dos; k=3 cuesta 13/6.La comparación favorece dos entre esos candidatos.
Con R=4
k=1 cuesta cuatro y k=2 cuesta cinco.Al cambiar R hemos cambiado el problema, no demostrado que un mismo coste tenga dos óptimos diferentes.
La selección depende de la penalización relativa, además de estabilizar.
Multiplicar Q=1, R=1 por cuatro da Q=4, R=4 y conserva k=1 como óptimo, escalando su valor.
3. La solución estabilizante necesita que los modos peligrosos se puedan tratar
Las hipótesis estándar son estabilizabilidad de (A, B) y detectabilidad de (Q¹ᐟ², A), con Q=(Q¹ᐟ²)ᵀQ¹ᐟ². La primera exige que los modos incontrolables ya sean estables. La segunda exige que los modos no penalizados por el coste de estado sean estables. R sigue siendo definida positiva.
Bajo esas condiciones existe una solución simétrica semidefinida positiva estabilizante de Riccati. No se afirma que P sea siempre definida positiva: una dirección estable sin penalización puede tener valor óptimo cero. Resolver una ecuación sin revisar estos modos podría confundir coste bajo con regulación al origen.
Un modo creciente sin actuador impide regular
x′=x, sin acción de control, Q=1 y x(0)≠0.
Trayectoria
x=x₀exp(t).Coste
∫₀∞x²dt diverge y ninguna entrada puede modificarlo.
Hipótesis
El modo +1 es incontrolable e inestable; la estabilizabilidad falla antes de intentar resolver Riccati.
Una matriz de pesos no crea una acción de control que el modelo no tiene.
El caso x₀=0 trivial no garantiza regular otros estados iniciales.
Una dirección estable no penalizada puede tener valor cero
x′=−x+u, Q=0 y R=1.
Política
u=0 deja x=x₀exp(−t), que ya converge al origen.Coste
J=∫u²=0, el menor valor posible.P=0 describe esa función de valor.Condición
El modo que Q no ve es −1, estable. Si A fuese +1 con el mismo Q=0, u=0 también costaría cero, pero no regularía: la detectabilidad habría fallado.
Semidefinido positivo admite cero sin invalidar el problema bajo sus hipótesis.
El segundo caso explica por qué un coste que ignora un modo creciente necesita especial atención.
4. Riccati selecciona una función de valor y una ganancia estabilizante
Su derivada es xᵀ(AᵀP+PA)x+2xᵀPBu. Minimizar coste instantáneo más esa derivada respecto a u da 2Ru+2BᵀPx=0, por tanto u=−Kx con K=R⁻¹BᵀP. R positiva hace único ese mínimo en u.
Se llama ecuación algebraica de Riccati. Entre posibles raíces se escoge la que hace A−BK Hurwitz y satisface la estructura de valor. Una raíz algebraica no estabilizante no se entrega como control óptimo de regulación.
Dos raíces escalares no representan dos reguladores óptimos
A=1, B=1, Q=3, R=1.
Ecuación
2P−P²+3=0, equivalente a (P−3)(P+1)=0.Raíz correcta
P=3 da K=3 y polo 1−3=−2, estable.Otra raíz
P=−1 daría u=+x y polo +2; su trayectoria y coste de estado crecen. No es una función de valor no negativa ni la solución estabilizante.
La selección usa estabilidad y significado del coste, no solo resolver el polinomio.
La ganancia es P aquí porque B=R=1; en otro problema se conserva la fórmula completa.
El peso del control también aparece en el divisor de K
Integrador A=0, B=1, Q=4 y R=4.
Riccati
−P²/4+4=0, con raíces ±4.Estabilizante
P=4 produce K=P/R=1 y polo −1.Comparar otro coste
Con R=1 y el mismo Q, P=2 y K=2.Son las dos curvas del dibujo; aumentar R cambió la política de este ejemplo.
P y K no son siempre el mismo número.
Sustituir P=4 en −P²/4+4 da cero, un control directo del factor R.
5. Completar cuadrados prueba la optimalidad conservando el término final
Es una identidad finita y exacta; el término final no se borra antes de revisar la trayectoria.
Para la política estabilizante, x(T)→0 y el cuadrado añadido es cero, así J=V(x₀). Para comparar con otras políticas admisibles se justifica también su condición terminal, por ejemplo en la clase de trayectorias estabilizantes con coste finito. Bajo las hipótesis estándar de LQR, los argumentos de detectabilidad permiten extender la caracterización a la clase de coste finito pertinente.
No se usa una trayectoria divergente para descartar sin explicación V final.
El valor predicho coincide con la integral del control óptimo
A=1, B=1, Q=3, R=1, P=3; x₀=2.
Trayectoria
u=−3x, x=2exp(−2t) y u=−6exp(−2t).Integral
3x²+u²=48exp(−4t), cuya integral es doce.Valor
V(x₀)=3·2²=12.Además 3x²+u²+d(3x²)/dt=(u+3x)², que se anula con esa política.
La ecuación de Riccati y la cuenta temporal dan el mismo valor.
El término terminal 3x(T)² tiende a cero por la trayectoria estable.
Una política estable diferente paga el cuadrado adicional
Mismo problema e inicio, pero usamos u=−4x.
Trayectoria
x=2exp(−3t), u=−8exp(−3t).Coste directo
∫(3x²+u²)=∫76exp(−6t)=38/3.Brecha
u+3x=−x, así la penalización extra integra 4/6=2/3.Doce más 2/3 coincide con 38/3; el término terminal vuelve a desaparecer.
Estabilizar no basta para minimizar el coste elegido.
Ambas políticas se comparan con las mismas Q, R, dinámica e inicio.
6. La política nominal se revisa con estados disponibles y límites reales
Después de resolver P se forma K, se verifican polos, trayectorias y esfuerzo. Los pesos no imponen por sí solos un máximo de fuerza, corriente o desplazamiento. Si se satura la orden, la ley aplicada cambia y se analiza esa ejecución. Cambiar los pesos crea un nuevo problema óptimo, aunque pueda ayudar a conseguir controles más suaves.
LQR de estado completo supone conocer x. Sustituirlo por x̂ requiere un observador y analizar error, ruido y coste de la combinación. La separación puede garantizar estabilidad ideal bajo condiciones, pero no convierte automáticamente esa combinación en el mismo óptimo de información completa. Diseños estadísticos necesitan un modelo de ruido adicional.
Un coste puede suavizar una orden, pero no es una restricción dura
Integrador, Q=4, x₀=2 y actuador limitado a |u|≤1.
Con R=1
K=2 y u(0)=−4: la orden nominal no cabe.Otro problema
Con R=16, P=8 y K=1/2.Para este dato concreto, u=−exp(−t/2), que sí cumple la cota.
Alcance
No se ha impuesto una desigualdad de actuador dentro de Riccati. Un dato inicial mayor podría volver a romperla; se revisan los datos y la política realmente aplicada.
Una elección de pesos puede ser útil sin sustituir el análisis de restricciones.
Los valores de coste de R=1 y R=16 no se comparan como si fueran el mismo objetivo.
Usar un estado estimado cambia el problema de información
A=1, B=1, Q=3, R=1; K=3, estado real x₀=0 y estimación x̂₀=1.
Estado completo
Con x₀=0, la política óptima conocida aplicaría u=0 y mantendría coste cero.
Estimación
Usar u=−3x̂ ordena inicialmente −3 y mueve la planta aunque su estado real empezara en el origen.
Conclusión
Que un observador converja después no elimina ese transitorio ni demuestra el mismo coste óptimo de estado completo.
Estabilidad de la combinación y optimalidad con otra información son propiedades diferentes.
No se atribuye una garantía estadística sin especificar probabilidades ni características del ruido.
7. Robustez cuantifica una propiedad para una familia declarada
Un nominal es un modelo de referencia. Podemos permitir un parámetro en un intervalo, un error aditivo o un error multiplicativo con forma y tamaño definidos. Estabilidad robusta pide estabilidad para todos los miembros; prestaciones robustas añaden una cota de error, esfuerzo u otra salida. Una familia de ganancias no incluye automáticamente retardos o dinámicas que no declaramos.
La garantía puede ser uniforme: las mismas constantes sirven para toda la familia. Eso es más informativo que simular unos pocos valores. Para familias matriciales generales, comprobar solo extremos no basta sin una propiedad que lo justifique; los próximos ejemplos escalares sí permiten una desigualdad para todos sus parámetros.
Certificar un intervalo completo con una desigualdad
x′=ax+u, a∈[1/2,3/2]; elegir u=−3x.
Polos
a−3 está entre −5/2 y −3/2 para todo a permitido.
Certificado
Con V=x², V′=2(a−3)x²≤−3V.Cota uniforme
|x(t)|≤exp(−3t/2)|x₀| para toda la familia nominal declarada, sin saturación.
La prueba abarca cada parámetro del intervalo, no solo tres simulaciones.
No extiende la familia por sí sola a un control retrasado o a otras ecuaciones.
Una incertidumbre en el actuador tiene otra forma
x′=x+bu, b∈[4/5,6/5], con la misma ley u=−3x.
Dinámica
El polo es 1−3b.
Intervalo
Varía entre −13/5 y −7/5, siempre negativo.
Interpretar
Hay estabilidad para esa familia de ganancia del actuador. No es exactamente la familia de a del ejemplo anterior; tiene otra ecuación y otro margen uniforme.
Dos garantías distintas deben indicar qué parámetro cambia y dónde aparece.
El extremo más lento aquí es −7/5, no −3/2.
8. La ganancia inducida especifica cómo medimos las señales
Una cota inducida ||salida||≤γ||entrada|| necesita una norma declarada y estado inicial cero. Podemos medir picos o energía; no son intercambiables en general. Para un LTI estable, la norma H infinito es su ganancia inducida L²: en escalar, sup_ω|H(iω)|. Para varias señales se usa la mayor singular de la matriz frecuencial.
Plancherel explica la cota energética: multiplicar cada frecuencia por una matriz de norma como máximo γ limita la energía total por γ². El supremo no se demuestra mirando una sola frecuencia arbitraria. Una senoidal permanente o una constante no tiene energía L² finita en toda la recta; pueden usarse señales largas de energía finita para aproximar algunas cotas, distinguiendo límite de máximo alcanzado.
Encontrar un supremo frecuencial completo
H=2/(s+2), estable.
Módulo
|H(iω)|=2/√(4+ω²)≤1 para todo ω real.Supremo
En ω=0 vale uno, así ||H||∞=1.Señales
La ganancia L² es uno. No se prueba usando como entrada una constante de energía infinita; señales de baja frecuencia y energía finita pueden aproximar ese cociente.
Se ha acotado toda la curva y localizado su supremo.
Que otro punto tenga módulo 0,4 no rebaja por sí solo esta norma global a 0,4.
La ganancia MIMO no es el mayor elemento de una matriz
H=[[1,1],[1,1]], paso directo estable.
Entrada de dirección fija
Para una señal escalar f de energía finita, ponemos w=f(1,1)ᵀ/√2; su norma L² es ||f||₂.
Salida
Hw=f(√2,√2)ᵀ, de norma 2||f||₂.Norma matricial
La mayor singular de H es dos, aunque cada elemento solo valga uno. Esa es la ganancia inducida energética.
Varias componentes pueden combinarse y amplificar una dirección más que un elemento aislado.
La dirección (1,−1) se anula, pero eso no reduce el supremo sobre todas las entradas a cero.
9. Pequeña ganancia controla las vueltas del lazo con una serie convergente
Consideramos z=M(w+Δz), con bloques causales estables, interconexión bien definida y ganancias inducidas compatibles γ₁ y γ₂.
En horizontes truncados y con causalidad se formula la garantía sin presuponer ya acotada la solución global.
Cada vuelta produce una cota geométrica de razón γ₁γ₂. La condición es estricta y suficiente; no confía en cancelaciones favorables de fase. Deben cumplirse estabilidad de los bloques y buena definición del lazo: una manipulación de normas aislada no construye por sí sola una interconexión válida.
Una familia de perturbaciones dinámicas queda cubierta por una cota
M=2/(s+2), ||M||∞=1, y Δ estable con ||Δ||∞≤2/5.
Producto
γ₁γ₂≤2/5<1.Ganancia cerrada
||z||₂≤[1/(1−2/5)]||w||₂=(5/3)||w||₂.Caso que alcanza el margen de amplitud baja
Si Δ=2/5 constante, la transferencia cerrada es 2/(s+6/5), con norma H infinito 5/3. M es estrictamente propia, evitando aquí un lazo algebraico directo.
La prueba cubre toda la familia estable de norma declarada, no solo esa constante.
La norma elegida es energética; no se reemplaza sin prueba por cualquier otra medida de señales.
Calcular producto y cota son dos operaciones distintas
γ₁=1/2, γ₂=6/5 y ||w||₂=2.
Producto
(1/2)(6/5)=3/5<1.Factor cerrado
γ₁/(1−γ₁γ₂)=(1/2)/(2/5)=5/4.Salida
||z||₂≤(5/4)·2=5/2.Por ejemplo M=1/(s+2) y Δ=6/5 constante en el segundo bloque dan norma cerrada 5/4.
El factor de la cota no es simplemente el producto de las dos ganancias.
El denominador positivo conserva el efecto de todas las vueltas del lazo.
10. Fallar una condición suficiente no demuestra inestabilidad de cada caso
Si el producto vale uno o más, pequeña ganancia no certifica con esa cota. Una interconexión concreta puede seguir siendo estable por fase o estructura. A la vez, no se declara robusta una familia grande solo porque uno de sus miembros pase una comprobación favorable.
Diseñar para reducir ||M||∞ permite obtener garantías frente a perturbaciones Δ incluidas en la familia. Estabilidad y prestaciones son niveles distintos: todos los modelos pueden ser estables y alguno incumplir una cota de error. Esta introducción no sustituye algoritmos de síntesis H infinito ni análisis de incertidumbre estructurada; deja claros los datos y la lógica de la garantía.
El mismo producto grande permite dos conclusiones temporales distintas
M=2/(s+2); comparar Δ=−2 con Δ=+2 en z=M(w+Δz).
Cota
Ambos tienen producto de normas dos, que no cumple el criterio estricto.
Signo negativo
Con Δ=−2, la transferencia cerrada es 2/(s+6), estable.Signo positivo
Con Δ=+2, es 2/(s−2), inestable.La familia que permite ambos no tiene garantía robusta, aunque uno de ellos funcione.
No certificación y demostración de inestabilidad no son sinónimos.
Con producto exactamente uno también habría que analizar el caso: el criterio seguiría sin certificar.
Una familia estable puede fallar una prestación más exigente
x′=(a−3)x+w, y=x, a∈[1/2,3/2].
Estabilidad
Todos los polos son negativos, como ya certificamos.
Ganancia de perturbación
H_a=1/(s+3−a), con ||H_a||∞=1/(3−a).El máximo de la familia es 2/3.
Prestación
Una cota energética 7/10 está garantizada para toda la familia; una cota 1/2 no, porque el miembro a=3/2 tiene ganancia 2/3.
Decir robustamente estable no incluye automáticamente cualquier objetivo de precisión.
La propiedad y el tamaño de la familia se especifican junto con la norma.
Ahora lo resolvemos de principio a fin.
EJEMPLO RESUELTO 1
LQR escalar
x′=u y coste de horizonte infinito ∫₀∞(x²+u²)dt. Calcula la solución estabilizante de Riccati.
Datos
Integrador A=0, B=1 y pesos Q=R=1, horizonte infinito.Riccati
−P²+1=0, con raíces uno y menos uno.Seleccionar
P=1 da K=1 y dinámica x′=−x, estable.P=−1 daría crecimiento.Control
u=−x.Certificado
x²+u²+d(x²)/dt=(u+x)².Valor
Con ese control el cuadrado adicional es cero y x² tiende a cero. El coste óptimo es x₀².
P=1 y u=−x.
La identidad x²+u²+d(x²)/dt=(u+x)² confirma que u=−x anula el término no negativo y la trayectoria estabilizante hace x² tender a cero.
EJEMPLO RESUELTO 2
Comprobar pequeña ganancia
Dos bloques estables tienen ganancias inducidas acotadas por 0,4 y 2. ¿Se cumple el criterio estricto?
Hipótesis
Se parte de bloques estables, normas inducidas compatibles e interconexión bien definida.
Producto
0,4·2=0,8.Criterio
0,8<1, así se cumple la condición estricta.Razón geométrica
Las cotas de vueltas sucesivas tienen razón 0,8.
Suma
1+0,8+0,8²+…=1/(1−0,8)=5.Conclusión
El criterio sí certifica bajo las hipótesis indicadas; el factor de una señal cerrada concreta conserva además dónde se aplica cada bloque.
Sí, el producto satisface el criterio.
La serie 1+0,8+0,8²+… converge a 1/(1−0,8)=5, coherente con la cota de realimentación del teorema.
Los símbolos que hemos usado
- x, u, A, B
- Estado, entrada y matrices constantes de la dinámica x′=Ax+Bu.
- J
- Coste acumulado que se intenta minimizar.
- Q, R
- Matrices simétricas de peso del estado y del control: Q semidefinida positiva y R definida positiva.
- P
- Matriz simétrica de la función de valor V(x)=xᵀPx, elegida estabilizante.
- K
- Ganancia de control óptimo, R⁻¹BᵀP, usada en u=−Kx.
- Δ
- Bloque desconocido que representa una incertidumbre dentro de una clase declarada.
- γ
- Cota de ganancia inducida: cuánto puede amplificar un bloque el tamaño de una señal según una norma fijada.
- ‖H‖∞
- Norma H infinito de una transferencia estable; es el máximo módulo frecuencial en el caso escalar y la mayor singular máxima en el matricial.
Ahora te toca a ti.
Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.
Practicar este tema →