La PizarraMatemática.

Máster y doctorado

Control de ecuaciones en derivadas parciales

Para calentar una barra no basta con indicar una temperatura: hay una en cada punto. Un control puede aportar calor en una zona o actuar en un extremo. Empezaremos por dibujar ese estado, después aislaremos una oscilación y finalmente veremos por qué controlar infinitas oscilaciones a la vez exige algo más.

1. El estado es un perfil y el actuador dice dónde podemos intervenir

En una barra 0<x<π, y(t,x) representa temperatura respecto a una referencia. En el modelo y_{t}=y_{xx}+u, con extremos fijados a cero y difusividad uno, y_{xx} describe la difusión y u añade una fuente. Fijar t deja una curva de temperaturas; fijar x deja una evolución temporal.

El estado es una temperatura en cada punto de la barra
yₜ=yₓₓ; extremos a cero0πt=0: sin(x)t=1: e⁻¹sin(x)Cada curva = un estado completo.
Perfiles de la solución libre y(t,x)=exp(−t)sin(x) en una barra de longitud pi con temperatura cero en los extremos y difusividad uno. Azul para tiempo cero y verde para tiempo uno, con los mismos ejes espaciales y de temperatura. Cada curva es un estado completo en un instante: sus puntos pertenecen a posiciones distintas de la barra. La amplitud baja de uno a exp menos uno. No son dos trayectorias temporales de una partícula ni se presupone que cualquier perfil sea un solo modo.

El espacio H=L²(0,π) mide el tamaño de un perfil por su norma cuadrática. Con control distribuido en toda la barra, u(t) también puede ser un perfil en H. Si solo tenemos una resistencia con forma b(x), la fuente es b(x)v(t): una señal escalar, con muchas menos decisiones.

El control de frontera requiere otro marco de operadores y trazas.

Verificar un perfil libre

y(t,x)=exp(−t)sin(x), sin fuente.

  1. Tiempo

    y_t=-e^{-t}\sin x.
  2. Espacio

    y_{xx}=-e^{-t}\sin x.
  3. Datos

    Coinciden las derivadas, los extremos son cero y al inicio queda sin(x).

Resultado

La curva completa decae sin cambiar de forma.

Lo comprobamos

No se está siguiendo una partícula de la barra.

Comparar dos actuadores

Fuente u(t,x) libre frente a sin(x)v(t).

  1. Libre

    Permite elegir un perfil espacial para cada instante.

  2. Fija

    Solo permite multiplicar sin(x) por una amplitud temporal.

  3. Consecuencia

    Una oscilación ortogonal a sin(x) no recibe directamente esa segunda fuente.

Resultado

La ecuación y el lugar de actuación forman juntos el sistema de control.

Lo comprobamos

No se transfieren resultados entre actuadores distintos.

2. La solución final suma lo que ocurre solo y lo que añade el control

Para y′=Ay+Bu, A genera el semigrupo S(t), que lleva un estado inicial a su evolución libre. En Hilbert, con B acotado de U a H y u en L²(0,T;U), la solución suave en sentido de semigrupos se escribe y(T)=S(T)y_0+R_Tu,\quad R_Tu=\int_0^T S(T-s)Bu(s)\,ds. No exige que el dato inicial esté en el dominio de A.

El control aplicado en s se propaga durante T−s antes de medir el resultado.

Para llegar a z hay que resolver R_Tu=z-S(T)y_0.

B acotado cubre aquí controles distribuidos; la misma escritura con controles de frontera necesita justificar admisibilidad, no sustituir B por un símbolo sin dominio.

Un modo con fuente constante

a′=−a+2, a(0)=0.

  1. Propagar

    El aporte de tiempo s llega multiplicado por exp(−(T−s)).

  2. Integrar

    a(T)=2\int_0^T e^{-(T-s)}ds=2(1-e^{-T}).
  3. Comprobar

    La función a(t)=2(1−exp(−t)) cumple a′=2exp(−t)=−a+2.
Resultado

Los aportes antiguos se han disipado más.

Lo comprobamos

Sumar simplemente 2T ignoraría la difusión.

Calcular el cambio necesario

a′=−2a+v, a(0)=1 y objetivo a(T)=0.

  1. Sin control

    La amplitud final sería exp(−2T).

  2. Diferencia

    El control debe aportar d=−exp(−2T).
  3. Ecuación

    \int_0^T e^{-2(T-s)}v(s)ds=-e^{-2T}.
Resultado

Cancelar el estado libre es un problema de alcanzabilidad.

Lo comprobamos

El objetivo cero no implica que el control sea cero.

3. Los modos separan ecuaciones, pero pueden compartir la misma señal

Los modos ortonormales del calor en esta barra son e_j(x)=\sqrt{2/\pi}\sin(jx), con decaimiento λⱼ=j².

Si y=\sum_j a_je_j y la fuente es b(x)v(t), proyectar da a_{j}′=−j²a_{j}+b_jv.

El coeficiente b_{j} mide cuánto toca el actuador a ese modo.

b_{j}=0 deja ese modo libre.

Si b_{j}≠0 la ecuación aislada puede llevarse a cualquier amplitud en tiempo positivo. En el conjunto, la señal v es la misma para todos: no podemos escoger una v distinta para cada ecuación y luego llamarlas un único control.

Un modo que no recibe señal

b(x)=e₁(x) y estado inicial e₂(x).

  1. Proyectar

    b₁=1 y b₂=0 por ortogonalidad.
  2. Segundo modo

    a₂′=−4a₂, con a₂(0)=1.
  3. Resultado

    a_2(T)=e^{-4T}>0 para T finito, con cualquier v.
Resultado

Ese actuador no lleva este estado exactamente a cero en tiempo finito.

Lo comprobamos

Decaer hacia cero cuando T crece es otra afirmación.

Dos ecuaciones pueden ocultar una diferencia incontrolable

a₁′=−a₁+v y a₂′=−a₂+v.

  1. Restar

    d=a₁−a₂ satisface d′=−d.
  2. Cancelar

    La señal común desaparece de la diferencia.

  3. Llegada

    Desde cero siempre a₁(T)=a₂(T); no se llega a (1,0).
Resultado

Tener dos coeficientes de actuación no nulos no basta en un sistema acoplado por una sola entrada.

Lo comprobamos

La multiplicidad de tasas importa en el análisis de rango.

4. El Gramiano de un modo traduce un objetivo en energía mínima

En a′=−λa+bv, la diferencia final d es el producto interior temporal entre v y h(s)=b exp(−λ(T−s)).

Su norma al cuadrado es W_T=\int_0^T b^2e^{-2\lambda s}ds.
Si λ>0, vale b²(1−exp(−2λT))/(2λ); si λ=0, vale b²T.
Cauchy–Schwarz da d²≤W_{T}∫v².
Si W_{T}>0 se alcanza la igualdad con v(s)=d h(s)/W_{T}.

Por tanto el coste mínimo ∫v² es d²/W_{T}. Aquí no hay factor 1/2 en el coste.

W_{T}=0 permite únicamente d=0.

Construir un control, no solo su coste

λ=1, b=2, T=1 y diferencia d=1.

  1. Gramiano

    W=2(1−exp(−2)).
  2. Control

    v(s)=2e^{-(1-s)}/W.
  3. Verificar

    El aporte final es (1/W)∫4exp(−2(1−s))ds=1; su energía es 1/W.
Resultado

La fórmula produce exactamente la diferencia requerida.

Lo comprobamos

Es mínima entre todas las señales L² del modo.

El caso sin disipación

λ=0, b=2, T=3 y d=6.

  1. Gramiano

    W=4·3=12.
  2. Control

    v=d·b/W=1, constante.
  3. Energía

    ∫₀³1²ds=3=d²/W.
Resultado

La fórmula integral también funciona cuando λ es cero.

Lo comprobamos

No se usa una expresión con división por λ=0.

5. Exacto, aproximado y nulo piden cosas distintas

Desde cero, control exacto significa R_{T} sobreyectivo sobre H. Aproximado significa que su rango es denso: cada objetivo se aproxima con la precisión pedida, sin una cota uniforme del coste. Control nulo exige poder cancelar S(T)y₀ para todo y₀, es decir \operatorname{Ran}S(T)\subseteq\operatorname{Ran}R_T.

Un ejemplo sencillo de rango denso no cerrado es D:ℓ²→ℓ², (Dx)_j=x_{j}/j.

Todo vector con un número finito de componentes está en su rango. Pero y_{j}=1/j pertenece a ℓ² y su única preimagen formal es x_{j}=1, que no pertenece a ℓ². Aproximar no es alcanzar exactamente.

Una sucesión de aproximaciones

Objetivo y_{j}=1/j para el operador D anterior.

  1. Truncar

    Tomamos x⁽ᴺ⁾ con sus N primeras componentes iguales a 1 y las otras cero.

  2. Resultado

    Dx⁽ᴺ⁾ coincide con y en esas N componentes; el error cuadrado es \sum_{j>N}1/j^2 y tiende a cero.

  3. Coste

    La norma de x⁽ᴺ⁾ es √N y crece sin límite.

Resultado

El objetivo se aproxima, aunque no hay una solución en ℓ².

Lo comprobamos

La convergencia del resultado no obliga a controles acotados.

Comparar dos peticiones

Desde y₀, pedimos terminar en cero o en cualquier z∈H.

  1. Cero

    Solo hay que producir −S(T)y₀.

  2. Cualquier z

    Hay que producir todas las diferencias z−S(T)y₀.

  3. Distinción

    La evolución libre puede haber suavizado mucho el dato que hay que cancelar.

Resultado

El control nulo puede existir sin control exacto sobre todo H.

Lo comprobamos

El calor proporciona el ejemplo explícito de la sección final.

6. El adjunto detecta direcciones y la observabilidad exige una cota uniforme

Para B acotado, (R_T^*z)(s)=B^*S(T-s)^*z.

Es la señal observada al propagar un dato terminal por el sistema adjunto.

El rango de R_{T} es denso si y solo si R_T^*z=0 obliga a z=0.

Eso expresa detección de toda dirección no nula.

El control exacto exige más: una constante C válida para todos los z con \|z\|^2\le C\int_0^T\|B^*S(T-s)^*z\|^2ds. Para control nulo se reemplaza el lado izquierdo por ∥S(T)*z∥². Detectar una señal pequeñísima no garantiza reconstruirla o controlarla con coste uniforme.

Detectable pero cada vez más débil

D de la sección anterior es autoadjunto y De_j=e_{j}/j.

  1. Núcleo

    Dz=0 implica z=0.
  2. Normas

    ∥e_{j}∥²=1 y ∥De_j∥²=1/j².
  3. Uniformidad

    La desigualdad 1≤C/j² no puede valer para todos los j con un mismo C.

Resultado

Hay densidad del rango, sin sobreyectividad.

Lo comprobamos

No basta con que cada cantidad de la derecha sea positiva.

Un modo escalar observado

a′=−λa+bv, dato adjunto terminal z.

  1. Observación

    Es b exp(−λ(T−s))z.

  2. Integrar

    Su energía es W_{Tz}².

  3. Cota

    Con W_{T}>0, z²≤(1/W_{T})W_{Tz}²; la constante es 1/W_{T}.
Resultado

El inverso del Gramiano mide la dificultad de esa dirección.

Lo comprobamos

Pasar a infinitos modos exige una constante independiente del índice.

7. En una EDP la geometría y el tiempo también forman parte del problema

En calor sobre un dominio acotado, conectado y suficientemente regular, con Dirichlet homogéneo y control interior en un abierto no vacío, hay teoremas de control nulo para todo T>0. La prueba cuantitativa requiere estimaciones globales; el coste puede crecer mucho al reducir T.

Las ondas se propagan a velocidad finita. Una observación que no recoge ciertas trayectorias de propagación puede fallar. Las condiciones geométricas de control y la duración dependen del modelo y sus fronteras. Una truncación con constantes que empeoran al aumentar modos no demuestra una propiedad uniforme de la EDP.

Una distancia impone una espera

Una señal viaja a velocidad c=2 y un punto está a distancia 3 de la zona actuada.

  1. Recorrido

    El tiempo para esa llegada directa es distancia/velocidad.

  2. Calcular

    3/2 unidades de tiempo.

  3. Alcance

    Antes de ese tiempo, esa influencia no ha llegado por ese recorrido.

Resultado

La propagación aporta una restricción temporal.

Lo comprobamos

Este cálculo no es por sí solo un teorema de control exacto de una onda con reflexiones.

Detectar una falsa conclusión por truncación

La estimación de N modos tiene constante C_{N}=N².

  1. Cada N

    Es una cota finita para ese problema truncado.

  2. Límite

    No hay una misma constante que cubra todos los N.

  3. Conclusión

    Falta una estimación uniforme para pasar a la EDP completa.

Resultado

Más modos pueden revelar una dificultad que un modelo pequeño ocultaba.

Lo comprobamos

No se reemplaza una prueba de operadores por estabilidad de una matriz aislada.

8. Estabilizar busca decaimiento, alcanzar fija un estado a una hora

En el calor de toda la barra con B=I, elegir u=−ky con k≥0 da a_{j}′=−(j²+k)a_{j}.

Por Parseval, ∥y(t)∥≤exp(−(1+k)t)∥y₀∥.

La cota vale para todos los modos y todos los datos en L², no solo para un dibujo.

Aunque la norma tienda a cero, para un modo no nulo queda una amplitud no nula a cualquier tiempo finito. Esa realimentación estabiliza y no ejecuta por sí sola un control nulo a tiempo fijado. Con actuación restringida, ni siquiera podemos aplicar arbitrariamente u=−ky como fuente en toda la barra.

Una garantía uniforme

B=I, k=2 y estado inicial de norma 4.

  1. Menor tasa

    j²+2≥3 para todo j≥1.
  2. Cota

    ∥y(t)∥≤4exp(−3t).
  3. Tiempo uno

    La norma es como máximo 4exp(−3).

Resultado

El modo más lento permite acotar simultáneamente los demás.

Lo comprobamos

La desigualdad no dice que todos los perfiles alcancen la igualdad.

Un límite no es una llegada finita

Estado inicial e₁ y k=2.

  1. Amplitud

    a₁(t)=exp(−3t).
  2. Finito

    Es estrictamente positiva para todo t finito.

  3. Límite

    Tiende a cero cuando t→∞.
Resultado

Estabilidad asintótica y control nulo son propiedades diferentes.

Lo comprobamos

Para cancelar exactamente en T se escoge otro control.

9. El calor con control completo muestra la diferencia de forma calculable

Volvamos a H=L²(0,π), B=I y controles L² en tiempo y espacio.

Ahora cada modo tiene su señal independiente.

Su Gramiano es W_j(T)=\frac{1-e^{-2j^2T}}{2j^2}.

Un objetivo d desde cero es alcanzable precisamente cuando Σd_{j}²/W_{j}(T)<∞; la suma es el coste mínimo, por la construcción escalar y Parseval.

Un modo alto puede alcanzarse, pero cuesta cada vez más
Calor; control en toda la barraT=1; coste mínimo ∫‖u‖²j=12.31j=28.00j=318.00Wⱼ=(1−e⁻²ʲ²)/(2j²)Coste de llegar a eⱼ: 1/WⱼMás oscilaciones, más esfuerzo.
Modelo del calor en cero pi con controles espaciales independientes en todos los modos ortonormales. Desde cero se pide un único modo final e j de norma uno, con horizonte uno. El coste mínimo es uno sobre W j, donde W j es uno menos exp menos dos j cuadrado, dividido por dos j cuadrado. Las barras para j uno, dos y tres están en una misma escala vertical y valen aproximadamente 2,31, 8,00 y 18,00. Cada objetivo individual es alcanzable; no existe una cota uniforme del coste por la norma del objetivo. La falta de control exacto de todo L2 se justifica después con un objetivo de infinitos modos, no solo por estas tres barras.
Para cualquier T>0, W_{j}→0.

El objetivo d_{j}=1/j está en ℓ², pero su coste exigiría sumar 2/(1−exp(−2j²T)), que diverge. En cambio, cancelar un estado inicial a requiere d_{j}=−exp(−j²T)a_{j}. El cociente exp(−2j²T)/W_{j} es acotado en j: su coste queda controlado por una constante dependiente de T por ∥a∥². Así hay control nulo y aproximado, pero no exacto sobre todo L².

Objetivos de norma uno con costes distintos

T=1 y objetivo e_{j} desde cero.

  1. Amplitud

    Solo la componente j vale uno.

  2. Coste

    Es 2j²/(1−exp(−2j²)).

  3. Comparar

    Para j=1 es aproximadamente 2,31; para j=2, aproximadamente 8,00.
Resultado

Cada modo individual se alcanza, con coste creciente.

Lo comprobamos

La imposibilidad para todos los objetivos se completa con d_{j}=1/j, no con solo dos valores.

Construir un control nulo de todo el perfil

Dato a∈ℓ² y T>0.

  1. Diferencias

    d_{j}=−exp(−j²T)a_{j}.
  2. Señales

    Usamos u_{j}(s)=d_{j} exp(−j²(T−s))/W_{j}(T).
  3. Suma

    Σ∫u_{j}²=Σ[2j²/(exp(2j²T)−1)]a_{j}²≤∥a∥²/T, porque exp(r)−1≥r.
Resultado

El control pertenece a L² y deja todos los modos a cero a la vez.

Lo comprobamos

La cota 1/T es suficiente en este modelo con actuación completa; no es el coste general con actuación localizada.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Modo invisible

Para a′=−2a+0·v, a(0)=1, halla a(t).

  1. Leer el dato

    La ecuación es a′=−2a+0·v y a(0)=1.
  2. Identificar actuación

    El término de control vale cero para cualquier señal v.

  3. Resolver

    La solución de a′=−2a es a(t)=Cexp(−2t).
  4. Usar el inicio

    a(0)=C=1.
  5. Resultado

    a(t)=exp(−2t), independientemente de v.
  6. Interpretar

    El actuador no elige la evolución de este modo; decae sin llegar a cero en tiempo finito.

Resultado

No se puede elegir su evolución mediante ese actuador.

Lo comprobamos

Cambiar v no modifica la ecuación.

EJEMPLO RESUELTO 2

Gramiano a tiempo infinito

Con λ=2,b=2, calcula lim_(T→∞)W_T.

  1. Fijar parámetros

    λ=2 y b=2.
  2. Escribir la integral

    W_T=\int_0^T4e^{-4s}ds.
  3. Integrar

    W_{T}=1−exp(−4T).
  4. Tomar el límite

    exp(−4T)→0 al crecer T.
  5. Resultado

    W_∞=1.
  6. Comprobar

    Es positivo porque b≠0 y finito porque la tasa de disipación es positiva.

Resultado

El Gramiano límite es 1.

Lo comprobamos

Es positivo: el modo está actuado.

Los símbolos que hemos usado

y′=Ay+Bv
Evolución abstracta en un Hilbert de estados con control v.
S(t)
Semigrupo que propaga el estado sin control.
R_T
Operador que lleva un control a su contribución al estado final.
W_T=R_TR_T*
Gramiano de control.
Control exacto
Alcanzar cualquier estado terminal del espacio indicado.
Control nulo
Llevar cualquier estado inicial permitido al estado cero.
Control aproximado
Aproximar estados terminales con la precisión del espacio escogido.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar