Las formulaciones débiles son herramientas importantes para el análisis de ecuaciones matemáticas que permiten la transferencia de conceptos del álgebra lineal para resolver problemas en otros campos como las ecuaciones diferenciales parciales . En una formulación débil, ya no es necesario que las ecuaciones o condiciones se cumplan absolutamente (y esto ni siquiera está bien definido) y, en cambio, tiene soluciones débiles sólo con respecto a ciertos "vectores de prueba" o " funciones de prueba ". En una formulación fuerte , el espacio de soluciones se construye de manera que estas ecuaciones o condiciones ya se cumplan.
El teorema de Lax-Milgram , que lleva el nombre de Peter Lax y Arthur Milgram , quienes lo demostraron en 1954, proporciona formulaciones débiles para ciertos sistemas en espacios de Hilbert .
Concepto general
Sea un espacio de Banach , sea el espacio dual de , let , [ se necesita aclaración ] y let . Un vector es una solución de la ecuación.![{\displaystyle V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle V'}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle A\dos puntos V\a V'}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle f\en V'}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle u\en V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle Au=f}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
si y sólo si para todos ,![{\displaystyle v\en V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle [Au](v)=f(v).}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Aquí, se denomina vector de prueba (en general) o función de prueba (si es un espacio funcional).![{\displaystyle v}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Para llevar esto a la forma genérica de una formulación débil, encuentre tal que![{\displaystyle u\en V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle a(u,v)=f(v)\quad \forall v\in V,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
definiendo la forma bilineal
![{\displaystyle a(u,v):=[Au](v).}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ejemplo 1: sistema lineal de ecuaciones
Ahora, sea y sea un mapeo lineal . Entonces, la formulación débil de la ecuación![{\displaystyle V=\mathbb {R} ^{n}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle A:V\a V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle Au=f}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
implica encontrar tal que para todos se cumpla la siguiente ecuación:![{\displaystyle u\en V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle v\en V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \langle Au,v\rangle =\langle f,v\rangle ,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde denota un producto interno .![{\displaystyle \langle \cdot ,\cdot \rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Como es un mapeo lineal, es suficiente probar con vectores base y obtenemos![{\displaystyle A}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \langle Au,e_{i}\rangle =\langle f,e_{i}\rangle ,\quad i=1,\ldots ,n.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
En realidad, ampliando , obtenemos la forma matricial de la ecuación![{\displaystyle u=\sum _ {j=1}^{n}u_ {j}e_ {j}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbf {A} \mathbf {u} =\mathbf {f},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
dónde y .![{\displaystyle a_{ij}=\langle Ae_{j},e_{i}\rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle f_{i}=\langle f,e_{i}\rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La forma bilineal asociada a esta formulación débil es
![{\displaystyle a(u,v)=\mathbf {v} ^{T}\mathbf {A} \mathbf {u} .}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ejemplo 2: ecuación de Poisson
Para resolver la ecuación de Poisson
![{\displaystyle -\nabla ^{2}u=f,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
en un dominio con su límite , y para especificar el espacio de solución más adelante, se puede usar el producto escalar![{\displaystyle \Omega \subset \mathbb {R} ^{d}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle u=0}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle L^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \langle u,v\rangle =\int _{\Omega }uv\,dx}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
para derivar la formulación débil. Luego, al probar con funciones diferenciables se obtiene![{\displaystyle v}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle -\int _{\Omega }(\nabla ^{2}u)v\,dx=\int _{\Omega }fv\,dx.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
El lado izquierdo de esta ecuación se puede hacer más simétrico integrando por partes usando la identidad de Green y asumiendo que en :![{\displaystyle v=0}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \partial \Omega }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \int _{\Omega }\nabla u\cdot \nabla v\,dx=\int _{\Omega }fv\,dx.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Esto es lo que se suele llamar formulación débil de la ecuación de Poisson . Las funciones en el espacio de solución deben ser cero en el límite y tener derivadas integrables al cuadrado . El espacio apropiado para satisfacer estos requisitos es el espacio de Sobolev de funciones con derivadas débiles en y con condiciones de frontera cero, entonces .
![{\displaystyle H_{0}^{1}(\Omega)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle L^{2}(\Omega)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle V=H_{0}^{1}(\Omega)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La forma genérica se obtiene asignando
![{\displaystyle a(u,v)=\int _{\Omega }\nabla u\cdot \nabla v\,dx}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y
![{\displaystyle f(v)=\int _{\Omega }fv\,dx.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
El teorema de Lax-Milgram
Ésta es una formulación del teorema de Lax-Milgram que se basa en las propiedades de la parte simétrica de la forma bilineal . No es la forma más general.
Sea un espacio de Hilbert y una forma bilineal en , que es![{\displaystyle V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle a(\cdot,\cdot)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
- acotado : y
![{\displaystyle |a(u,v)|\leq C\|u\|\|v\|\,;}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
- coercitivo :
![{\displaystyle a(u,u)\geq c\|u\|^{2}\,.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Entonces, para cualquier , existe una única solución a la ecuación![{\displaystyle f\en V'}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle u\en V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle a(u,v)=f(v)\quad \forall v\in V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y se sostiene
![{\displaystyle \|u\|\leq {\frac {1}{c}}\|f\|_{V'}\,.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Aplicación al ejemplo 1
En este caso, la aplicación del teorema de Lax-Milgram es un resultado más sólido de lo necesario.
- Acotación: todas las formas bilineales están acotadas. En particular, tenemos
![{\displaystyle \mathbb {R} ^{n}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle |a(u,v)|\leq \|A\|\,\|u\|\,\|v\|}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
- Coercitividad: esto en realidad significa que las partes reales de los valores propios de no son menores que . Dado que esto implica en particular que ningún valor propio es cero, el sistema tiene solución.
![{\displaystyle A}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle c}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Además, esto produce la estimación
![{\displaystyle \|u\|\leq {\frac {1}{c}}\|f\|,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
.![{\displaystyle c}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle A}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Aplicación al ejemplo 2
Aquí, elige con la norma.![{\displaystyle V=H_{0}^{1}(\Omega)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \|v\|_{V}:=\|\nabla v\|,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde la norma de la derecha es la norma - (esto proporciona una norma verdadera según la desigualdad de Poincaré ). Pero vemos que y por la desigualdad de Cauchy-Schwarz ,.![{\displaystyle L^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Omega}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle |a(u,u)|=\|\nabla u\|^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle |a(u,v)|\leq \|\nabla u\|\,\|\nabla v\|}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Por lo tanto, para cualquier , existe una solución única de la ecuación de Poisson y tenemos la estimación![{\displaystyle f\in [H_{0}^{1}(\Omega)]'}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle u\en V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \|\nabla u\|\leq \|f\|_{[H_{0}^{1}(\Omega )]'}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ver también
Referencias
- Lax, Peter D .; Milgram, Arthur N. (1954), "Ecuaciones parabólicas", Contribuciones a la teoría de ecuaciones diferenciales parciales , Annals of Mathematics Studies, vol. 33, Princeton, Nueva Jersey : Princeton University Press , págs. 167–190, doi :10.1515/9781400882182-010, ISBN 9781400882182, SEÑOR 0067317, Zbl 0058.08703
enlaces externos
- Página de MathWorld sobre el teorema de Lax-Milgram