Weak formulations are important tools for the analysis of mathematical equations that permit the transfer of concepts of linear algebra to solve problems in other fields such as partial differential equations. In a weak formulation, equations or conditions are no longer required to hold absolutely (and this is not even well defined) and has instead weak solutions only with respect to certain "test vectors" or "test functions". In a strong formulation, the solution space is constructed such that these equations or conditions are already fulfilled.
The Lax–Milgram theorem, named after Peter Lax and Arthur Milgram who proved it in 1954, provides weak formulations for certain systems on Hilbert spaces.
Let V {\displaystyle V} be a Banach space, let V ′ {\displaystyle V'} be the dual space of V {\displaystyle V} , let A : : --> V → → --> V ′ {\displaystyle A\colon V\to V'} ,[clarification needed] and let f ∈ ∈ --> V ′ {\displaystyle f\in V'} . A vector u ∈ ∈ --> V {\displaystyle u\in V} is a solution of the equation
A u = f {\displaystyle Au=f}
if and only if for all v ∈ ∈ --> V {\displaystyle v\in V} ,
( A u ) ( v ) = f ( v ) . {\displaystyle (Au)(v)=f(v).}
A particular choice of v {\displaystyle v} is called a test vector (in general) or a test function (if V {\displaystyle V} is a function space).
To bring this into the generic form of a weak formulation, find u ∈ ∈ --> V {\displaystyle u\in V} such that
a ( u , v ) = f ( v ) ∀ ∀ --> v ∈ ∈ --> V , {\displaystyle a(u,v)=f(v)\quad \forall v\in V,}
by defining the bilinear form
a ( u , v ) := ( A u ) ( v ) . {\displaystyle a(u,v):=(Au)(v).}
Now, let V = R n {\displaystyle V=\mathbb {R} ^{n}} and A : V → → --> V {\displaystyle A:V\to V} be a linear mapping. Then, the weak formulation of the equation
involves finding u ∈ ∈ --> V {\displaystyle u\in V} such that for all v ∈ ∈ --> V {\displaystyle v\in V} the following equation holds:
⟨ ⟨ --> A u , v ⟩ ⟩ --> = ⟨ ⟨ --> f , v ⟩ ⟩ --> , {\displaystyle \langle Au,v\rangle =\langle f,v\rangle ,}
where ⟨ ⟨ --> ⋅ ⋅ --> , ⋅ ⋅ --> ⟩ ⟩ --> {\displaystyle \langle \cdot ,\cdot \rangle } denotes an inner product.
Since A {\displaystyle A} is a linear mapping, it is sufficient to test with basis vectors, and we get
⟨ ⟨ --> A u , e i ⟩ ⟩ --> = ⟨ ⟨ --> f , e i ⟩ ⟩ --> , i = 1 , … … --> , n . {\displaystyle \langle Au,e_{i}\rangle =\langle f,e_{i}\rangle ,\quad i=1,\ldots ,n.}
Actually, expanding u = ∑ ∑ --> j = 1 n u j e j {\displaystyle u=\sum _{j=1}^{n}u_{j}e_{j}} , we obtain the matrix form of the equation
A u = f , {\displaystyle \mathbf {A} \mathbf {u} =\mathbf {f} ,}
where a i j = ⟨ ⟨ --> A e j , e i ⟩ ⟩ --> {\displaystyle a_{ij}=\langle Ae_{j},e_{i}\rangle } and f i = ⟨ ⟨ --> f , e i ⟩ ⟩ --> {\displaystyle f_{i}=\langle f,e_{i}\rangle } .
The bilinear form associated to this weak formulation is
a ( u , v ) = v T A u . {\displaystyle a(u,v)=\mathbf {v} ^{T}\mathbf {A} \mathbf {u} .}
To solve Poisson's equation
− − --> ∇ ∇ --> 2 u = f , {\displaystyle -\nabla ^{2}u=f,}
on a domain Ω Ω --> ⊂ ⊂ --> R d {\displaystyle \Omega \subset \mathbb {R} ^{d}} with u = 0 {\displaystyle u=0} on its boundary, and to specify the solution space V {\displaystyle V} later, one can use the L 2 {\displaystyle L^{2}} -scalar product
⟨ ⟨ --> u , v ⟩ ⟩ --> = ∫ ∫ --> Ω Ω --> u v d x {\displaystyle \langle u,v\rangle =\int _{\Omega }uv\,dx}
to derive the weak formulation. Then, testing with differentiable functions v {\displaystyle v} yields
− − --> ∫ ∫ --> Ω Ω --> ( ∇ ∇ --> 2 u ) v d x = ∫ ∫ --> Ω Ω --> f v d x . {\displaystyle -\int _{\Omega }(\nabla ^{2}u)v\,dx=\int _{\Omega }fv\,dx.}
The left side of this equation can be made more symmetric by integration by parts using Green's identity and assuming that v = 0 {\displaystyle v=0} on ∂ ∂ --> Ω Ω --> {\displaystyle \partial \Omega } :
∫ ∫ --> Ω Ω --> ∇ ∇ --> u ⋅ ⋅ --> ∇ ∇ --> v d x = ∫ ∫ --> Ω Ω --> f v d x . {\displaystyle \int _{\Omega }\nabla u\cdot \nabla v\,dx=\int _{\Omega }fv\,dx.}
This is what is usually called the weak formulation of Poisson's equation. Functions in the solution space V {\displaystyle V} must be zero on the boundary, and have square-integrable derivatives. The appropriate space to satisfy these requirements is the Sobolev space H 0 1 ( Ω Ω --> ) {\displaystyle H_{0}^{1}(\Omega )} of functions with weak derivatives in L 2 ( Ω Ω --> ) {\displaystyle L^{2}(\Omega )} and with zero boundary conditions, so V = H 0 1 ( Ω Ω --> ) {\displaystyle V=H_{0}^{1}(\Omega )} .
The generic form is obtained by assigning
a ( u , v ) = ∫ ∫ --> Ω Ω --> ∇ ∇ --> u ⋅ ⋅ --> ∇ ∇ --> v d x {\displaystyle a(u,v)=\int _{\Omega }\nabla u\cdot \nabla v\,dx}
and
f ( v ) = ∫ ∫ --> Ω Ω --> f v d x . {\displaystyle f(v)=\int _{\Omega }fv\,dx.}
This is a formulation of the Lax–Milgram theorem which relies on properties of the symmetric part of the bilinear form. It is not the most general form.
Let V {\displaystyle V} be a Hilbert space and a ( ⋅ ⋅ --> , ⋅ ⋅ --> ) {\displaystyle a(\cdot ,\cdot )} a bilinear form on V {\displaystyle V} , which is
Then, for any bounded f ∈ ∈ --> V ′ {\displaystyle f\in V'} , there is a unique solution u ∈ ∈ --> V {\displaystyle u\in V} to the equation
a ( u , v ) = f ( v ) ∀ ∀ --> v ∈ ∈ --> V {\displaystyle a(u,v)=f(v)\quad \forall v\in V}
and it holds
‖ ‖ --> u ‖ ‖ --> ≤ ≤ --> 1 c ‖ ‖ --> f ‖ ‖ --> V ′ . {\displaystyle \|u\|\leq {\frac {1}{c}}\|f\|_{V'}\,.}
Here, application of the Lax–Milgram theorem is a stronger result than is needed.
Additionally, this yields the estimate ‖ ‖ --> u ‖ ‖ --> ≤ ≤ --> 1 c ‖ ‖ --> f ‖ ‖ --> , {\displaystyle \|u\|\leq {\frac {1}{c}}\|f\|,} where c {\displaystyle c} is the minimal real part of an eigenvalue of A {\displaystyle A} .
Here, choose V = H 0 1 ( Ω Ω --> ) {\displaystyle V=H_{0}^{1}(\Omega )} with the norm ‖ ‖ --> v ‖ ‖ --> V := ‖ ‖ --> ∇ ∇ --> v ‖ ‖ --> , {\displaystyle \|v\|_{V}:=\|\nabla v\|,}
where the norm on the right is the L 2 {\displaystyle L^{2}} -norm on Ω Ω --> {\displaystyle \Omega } (this provides a true norm on V {\displaystyle V} by the Poincaré inequality). But, we see that | a ( u , u ) | = ‖ ‖ --> ∇ ∇ --> u ‖ ‖ --> 2 {\displaystyle |a(u,u)|=\|\nabla u\|^{2}} and by the Cauchy–Schwarz inequality, | a ( u , v ) | ≤ ≤ --> ‖ ‖ --> ∇ ∇ --> u ‖ ‖ --> ‖ ‖ --> ∇ ∇ --> v ‖ ‖ --> {\displaystyle |a(u,v)|\leq \|\nabla u\|\,\|\nabla v\|} .
Therefore, for any f ∈ ∈ --> [ H 0 1 ( Ω Ω --> ) ] ′ {\displaystyle f\in [H_{0}^{1}(\Omega )]'} , there is a unique solution u ∈ ∈ --> V {\displaystyle u\in V} of Poisson's equation and we have the estimate
‖ ‖ --> ∇ ∇ --> u ‖ ‖ --> ≤ ≤ --> ‖ ‖ --> f ‖ ‖ --> [ H 0 1 ( Ω Ω --> ) ] ′ . {\displaystyle \|\nabla u\|\leq \|f\|_{[H_{0}^{1}(\Omega )]'}.}
Lokasi Pengunjung: 3.145.52.10