弱形式




数学において弱形式(じゃくけいしき、英: weak formulation)は、線型代数学の概念を、例えば偏微分方程式などの他の分野において問題を解くために用いることを可能にする、重要な解析上の道具です。弱形式において、方程式の絶対性はもはや要求されず(適切である必要すらない)、代わりにあるテストベクトルあるいはテスト函数に関する弱解が存在します。これは超函数の意味で解を要求する問題を構成することと同値です。


ここでは弱形式に関するいくつかの例を紹介し、その解に対する主要な定理であるラックス=ミルグラムの定理(Lax-Milgram theorem)を述べます。




目次






  • 1 一般の概念


  • 2 例1:線型連立方程式


  • 3 例2 ポアソン方程式


  • 4 ラックス=ミルグラムの定理


    • 4.1 例1への応用


    • 4.2 例2への応用




  • 5 関連項目


  • 6 参考文献


  • 7 外部リンク





一般の概念


V{displaystyle V} をあるバナッハ空間とする。次の方程式の解 u∈V{displaystyle uin V} を見つけたいです。



Au=f{displaystyle Au=f}

但し A:V→V′{displaystyle A:Vto V'} および f∈V′{displaystyle fin V'} であり、V′{displaystyle V'}V{displaystyle V} の双対です。


変分法によると、この問題はすべての v∈V{displaystyle vin V} に対して次を満たす u∈V{displaystyle uin V} を見つけることと同値です。



[Au](v)=f(v){displaystyle [Au](v)=f(v)}.。

ここで v{displaystyle v} をテストベクトルあるいはテスト函数と呼びます。


これを弱形式による一般的な形に書き換えます。すなわち、次を満たす u∈V{displaystyle uin V} を見つけます。



a(u,v)=f(v)∀v∈V{displaystyle a(u,v)=f(v)quad forall vin V}

但し、双線型形式


a(u,v):=[Au](v){displaystyle a(u,v):=[Au](v)}

が定義されます。以上の説明は非常に抽象的であるため、以下ではいくつかの例を見ます。



例1:線型連立方程式


V=Rn{displaystyle V=mathbb {R} ^{n}}A:V→V{displaystyle A:Vto V} を線型写像とします。このとき、方程式


Au=f{displaystyle Au=f}

の弱形式は、すべての v∈V{displaystyle vin V} に対して次の方程式を満たす u∈V{displaystyle uin V} を見つけることとなります。



Au,v⟩=⟨f,v⟩{displaystyle langle Au,vrangle =langle f,vrangle ,}

ここで ,⋅{displaystyle langle cdot ,cdot rangle } は内積を表します。


A{displaystyle A} は線型写像なので、基底ベクトルに対して調べれば十分です。すると


Au,ei⟩=⟨f,ei⟩i=1,…,n{displaystyle langle Au,e_{i}rangle =langle f,e_{i}rangle quad i=1,ldots ,n,}

が得られます。実際、u=∑j=1nujej{displaystyle u=sum _{j=1}^{n}u_{j}e_{j}} と展開することで、次の行列の形式での方程式が得られます。



Au=f.{displaystyle mathbf {A} mathbf {u} =mathbf {f} .}

ここで aij=⟨Aej,ei⟩{displaystyle a_{ij}=langle Ae_{j},e_{i}rangle } および fi=⟨f,ei⟩{displaystyle f_{i}=langle f,e_{i}rangle } です。


この弱形式に関連する双線型形式は、次で与えられます。



a(u,v)=vTAu.{displaystyle a(u,v)=mathbf {v} ^{T}mathbf {A} mathbf {u} .}


例2 ポアソン方程式


ここでの目標は、ある領域 ΩRd{displaystyle Omega subset mathbb {R} ^{d}} 上の次のポアソン方程式


2u=f{displaystyle -nabla ^{2}u=f,}

の解で、境界で u=0{displaystyle u=0} となるようなものを見つけることです。また解空間 V{displaystyle V} は後述の議論で決定します。弱形式の導出のために、次の L2{displaystyle L^{2}}-スカラー内積を用います。



u,v⟩=∫Ωuvdx.{displaystyle langle u,vrangle =int _{Omega }uv,dx.}

微分可能な函数 v{displaystyle v} をテスト函数として用いることで、次が得られます。



Ω(∇2u)vdx=∫Ωfvdx.{displaystyle -int _{Omega }(nabla ^{2}u)v,dx=int _{Omega }fv,dx.}。この方程式の左辺は、グリーンの恒等式を用いた部分積分により、より対称的な次の形式で記述できる。

Ωu⋅vdx=∫Ωfvdx.{displaystyle int _{Omega }nabla ucdot nabla v,dx=int _{Omega }fv,dx.}

これは正しくポアソン方程式の弱形式と通常呼ばれるものである。ここで空間 V{displaystyle V} を定義する必要がある。この空間は、この方程式を導けるものでなければならない。したがってこの空間における導函数は二乗可積分である必要がある。実際、ゼロ境界条件で、弱微分が L2(Ω){displaystyle L^{2}(Omega )} に属す函数からなるソボレフ空間 H01(Ω){displaystyle H_{0}^{1}(Omega )} を考えれば、目的は満たされる。


次のように記号を定めることで、一般的な形を得ることが出来る:


a(u,v)=∫Ωu⋅vdx{displaystyle a(u,v)=int _{Omega }nabla ucdot nabla v,dx}

および


f(v)=∫Ωfvdx.{displaystyle f(v)=int _{Omega }fv,dx.}


ラックス=ミルグラムの定理


これは双線型形式の対称部分の性質に依存するラックス=ミルグラムの定理(Lax-Milgram theorem)の構成である。最も一般的な形という訳ではない。


V{displaystyle V} をヒルベルト空間とし、a(⋅,⋅){displaystyle a(cdot ,cdot )}V{displaystyle V} 上の双線型形式で、次を満たすものとする:




  1. 有界:


  2. 強圧的: a(u,u)≥c‖u‖2.{displaystyle a(u,u)geq c|u|^{2}.}


このとき、任意の f∈V′{displaystyle fin V'} に対して、次の方程式には唯一つの解 u∈V{displaystyle uin V} が存在する。


a(u,v)=f(v).{displaystyle a(u,v)=f(v).}

また次が成立する。


u‖1c‖f‖V′.{displaystyle |u|leq {frac {1}{c}}|f|_{V'}.}


例1への応用


この場合、ラックス=ミルグラムの定理を適用することは明らかに十分すぎるものであるが、他の場合と同様の形にするためにこの定理を使用する。



  • 有界性: Rn{displaystyle mathbb {R} ^{n}} 上のすべての双線型形式は有界である。特に、次が成り立つ。
    |a(u,v)|≤A‖u‖v‖{displaystyle |a(u,v)|leq |A|,|u|,|v|,}


  • 強圧性: これは実際、A{displaystyle A} の固有値の実部が c{displaystyle c} よりも小さくないことを意味する。これは特に、ゼロ固有値が存在しないことを意味するので、系は可解である。


さらに次の評価が得られる。


u‖1c‖f‖,{displaystyle |u|leq {frac {1}{c}}|f|,,}

ここで c{displaystyle c}A{displaystyle A} の固有値の最小実部である。



例2への応用


上述のように、V=H01(Ω){displaystyle V=H_{0}^{1}(Omega )} とし、ノルムは次で定める。


v‖V:=‖v‖{displaystyle |v|_{V}:=|nabla v|}

ここで右辺のノルムは Ω{displaystyle Omega } 上での L2{displaystyle L^{2}}-ノルムである(ポアンカレ不等式により、これは正しく V{displaystyle V} 上のノルムを与える)。しかし、|a(u,u)|=‖u‖2{displaystyle |a(u,u)|=|nabla u|^{2}} であり、コーシー=シュワルツの不等式より次


したがって、任意の f∈[H01(Ω)]′{displaystyle fin [H_{0}^{1}(Omega )]'} に対して、ポアソン方程式の唯一つの解 u∈V{displaystyle uin V} が存在し、次の評価が得られる。


u‖f‖[H01(Ω)]′.{displaystyle |nabla u|leq |f|_{[H_{0}^{1}(Omega )]'}.}


関連項目



  • バフスカ=ラックス=ミルグラムの定理

  • リオンス=ラックス=ミルグラムの定理(英語版)



参考文献



  • Lax, Peter D.; Milgram, Arthur N. (1954). “Parabolic equations”. Contributions to the theory of partial differential equations. Annals of Mathematics Studies, no. 33. Princeton, N. J.: Princeton University Press. pp. 167–190.  MR.mw-parser-output cite.citation{font-style:inherit}.mw-parser-output .citation q{quotes:"""""""'""'"}.mw-parser-output .citation .cs1-lock-free a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/6/65/Lock-green.svg/9px-Lock-green.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .citation .cs1-lock-limited a,.mw-parser-output .citation .cs1-lock-registration a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/d/d6/Lock-gray-alt-2.svg/9px-Lock-gray-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .citation .cs1-lock-subscription a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/a/aa/Lock-red-alt-2.svg/9px-Lock-red-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration{color:#555}.mw-parser-output .cs1-subscription span,.mw-parser-output .cs1-registration span{border-bottom:1px dotted;cursor:help}.mw-parser-output .cs1-ws-icon a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/4/4c/Wikisource-logo.svg/12px-Wikisource-logo.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output code.cs1-code{color:inherit;background:inherit;border:inherit;padding:inherit}.mw-parser-output .cs1-hidden-error{display:none;font-size:100%}.mw-parser-output .cs1-visible-error{font-size:100%}.mw-parser-output .cs1-maint{display:none;color:#33aa33;margin-left:0.3em}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration,.mw-parser-output .cs1-format{font-size:95%}.mw-parser-output .cs1-kern-left,.mw-parser-output .cs1-kern-wl-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right,.mw-parser-output .cs1-kern-wl-right{padding-right:0.2em}
    0067317


外部リンク


  • MathWorld page on Lax–Milgram theorem



Popular posts from this blog

Human spaceflight

Can not write log (Is /dev/pts mounted?) - openpty in Ubuntu-on-Windows?

File:DeusFollowingSea.jpg