微分方程式 と幾何学解析 といった数学の分野において、最大値原理 は最も有用かつよく知られた研究ツールの一つである。領域Dにおける 偏微分方程式(あるいはより一般的には微分不等式)の解は、 D の境界で最大値に達するとき、最大値原理を 満たすと言われる。調和関数、そしてより一般的には 楕円型偏微分方程式 の解は、最大値原理を満たす。
最大原理を満たす 2 次元の典型的な関数のグラフ (赤): 最大値 (および最小値) は、ドメイン (青) の境界上で発生します。 最大原理は、微分方程式の解そのものを明示的に知ることなしに、その解に関する情報を得ることを可能にする。特に、最大原理は常微分方程式および偏微分方程式の解の数値近似、およびそのような近似における誤差の境界値を決定する際に有用なツールである。[ 1 ]
単純な2次元の場合、2つの変数u ( x , y ) の関数を考えます。
∂ 2 あなた ∂ × 2 + ∂ 2 あなた ∂ y 2 = 0。 {\displaystyle {\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}=0.} この設定における弱い最大値原理は、 u の定義域の任意の有界開部分集合Mに対して、 M の閉包上のu の最大値はM の境界上で達成される、というものである。強い最大値原理は 、u が 定数関数でない限り、M自身の上ではどこでも最大値は達成されない、というものである。なお、これらの記述は u の最小値についても成り立つ。なぜなら、 - u は 同じ微分方程式を解くからである。
凸最適化 の分野では、コンパクト 凸集合 上の凸関数の最大値は 境界 上で達成されると主張する類似の主張がある。[ 2 ]
直感
ここでは最も単純なケースを考えますが、同じ考え方をより一般的なシナリオに拡張することができます。Mを ユークリッド空間 の開部分集合とし、uを M上の C 2 関数で、
∑ 私 = 1 n ∑ j = 1 n 1つの 私 j ∂ 2 あなた ∂ × 私 ∂ × j = 0 {\displaystyle \sum _{i=1}^{n}\sum _{j=1}^{n}a_{ij}{\frac {\partial ^{2}u}{\partial x^{i}\,\partial x^{j}}}=0} ここで、 1からn までの各i とj に対して、a ij はM 上の関数であり、 a ij = a ji となります。
M におけるx の適当な値を固定する。線型代数 のスペクトル定理 によれば、行列[ a ij ( x )] のすべての固有値は実数であり、ℝ n には固有ベクトルからなる直交基底 が存在する。iが 1 からnまでで ある とき、固有値をλ i 、対応する固有ベクトルをv i で表す。すると、点x における微分方程式は次のように書き直される 。
∑ 私 = 1 n λ 私 d 2 d t 2 | t = 0 ( あなた ( × + t v 私 ) ) = 0。 {\displaystyle \sum _{i=1}^{n}\lambda _{i}\left.{\frac {d^{2}}{dt^{2}}}\right|_{t=0}{\big (}u(x+tv_{i}){\big )}=0.} 最大値原理の本質は、各固有値が正の場合(これは微分方程式の「楕円性」の特定の定式化に相当します)、上記の式は解の方向2次微分に一定の釣り合いを課すという単純な観察にあります。特に、方向2次微分のうちの1つが負の場合、もう1つは必ず正になります。uが 最大となる仮想的な点では、すべての方向2次微分は自動的に非正となり、上記の式で表される「釣り合い」には、すべての方向2次微分が0となることが必要です。
この基本的な推論は、いくつかの追加の仮定 ( a の連続性など) の下で、u が最大化される M の点がある場合、 u は 一定でなければならないことを述べる、強い最大値原理の無限小定式化を表すと主張できます。
上記の推論は、より一般的な偏微分方程式を考えた場合にも影響されないことに注意されたい。
∑ 私 = 1 n ∑ j = 1 n 1つの 私 j ∂ 2 あなた ∂ × 私 ∂ × j + ∑ 私 = 1 n b 私 ∂ あなた ∂ × 私 = 0 、 {\displaystyle \sum _{i=1}^{n}\sum _{j=1}^{n}a_{ij}{\frac {\partial ^{2}u}{\partial x^{i}\,\partial x^{j}}}+\sum _{i=1}^{n}b_{i}{\frac {\partial u}{\partial x^{i}}}=0,} 追加された項は、任意の仮想的な最大点において自動的にゼロとなるためである。より一般的な条件を考慮しても、この推論は影響を受けない。
∑ 私 = 1 n ∑ j = 1 n 1つの 私 j ∂ 2 あなた ∂ × 私 ∂ × j + ∑ 私 = 1 n b 私 ∂ あなた ∂ × 私 ≥ 0 、 {\displaystyle \sum _{i=1}^{n}\sum _{j=1}^{n}a_{ij}{\frac {\partial ^{2}u}{\partial x^{i}\,\partial x^{j}}}+\sum _{i=1}^{n}b_{i}{\frac {\partial u}{\partial x^{i}}}\geq 0,} この条件において、仮想的な最大点において厳密な不等式( ≥ ではなく> )が存在する場合、完全な矛盾が生じるという追加の現象さえも指摘できる。この現象は、古典的な弱最大原理の形式的証明において重要である。
強い最大値原理の適用不可 しかし、上記の推論は、条件を考慮すると当てはまらなくなる。
∑ 私 = 1 n ∑ j = 1 n 1つの 私 j ∂ 2 あなた ∂ × 私 ∂ × j + ∑ 私 = 1 n b 私 ∂ あなた ∂ × 私 ≤ 0 、 {\displaystyle \sum _{i=1}^{n}\sum _{j=1}^{n}a_{ij}{\frac {\partial ^{2}u}{\partial x^{i}\,\partial x^{j}}}+\sum _{i=1}^{n}b_{i}{\frac {\partial u}{\partial x^{i}}}\leq 0,} なぜなら、 u の仮想的な最大点で評価される「釣り合い」条件は、明らかに非正の量の加重平均が非正であるということを示しているに過ぎないからである。これは自明に真であり、したがって、そこから非自明な結論を導くことはできない。これは、例えば以下のような、数多くの具体的な例に反映されている。
∂ 2 ∂ × 2 ( − × 2 − y 2 ) + ∂ 2 ∂ y 2 ( − × 2 − y 2 ) ≤ 0 、 {\displaystyle {\frac {\partial ^{2}}{\partial x^{2}}}{\big (}{-x}^{2}-y^{2}{\big )}+{\frac {\partial ^{2}}{\partial y^{2}}}{\big (}{-x}^{2}-y^{2}{\big )}\leq 0,} そして原点を含む任意の開領域上では関数− x 2 − y 2 は必ず最大値を持ちます。
線形楕円偏微分方程式の古典的な弱最大原理
本質的な考え方 M をユークリッド空間の開集合とする。滑らかな関数が点p で最大化されるとき、以下の関係が成り立つ。 あなた : M → R {\displaystyle u:M\to \mathbb {R} }
( d あなた ) ( p ) = 0 {\displaystyle (du)(p)=0} ( ∇ 2 あなた ) ( p ) ≤ 0 、 {\displaystyle (\nabla ^{2}u)(p)\leq 0,} 行列不等式として。偏微分方程式は、 関数の様々な導関数の間に代数的関係を課すものと見ることができます。したがって、 u が 偏微分方程式の解である場合、u の1次および2次導関数に対する上記の条件は、この代数的関係に矛盾を生じる可能性があります。これが最大値原理の本質です。明らかに、この考え方の適用可能性は、問題となる特定の偏微分方程式に大きく依存します。
例えば、uが 微分方程式を解くとする と
Δ あなた = | d あなた | 2 + 2 、 {\displaystyle \Delta u=|du|^{2}+2,} と が定義域のどの点においても成り立つことは明らかに不可能です。したがって、上記の観察に従えば、u が 最大値を取ることは不可能です。もしu が 微分方程式を解けば、このような矛盾は生じず、これまでの分析からは何ら興味深い結果は得られません。もしu が 微分方程式を解けば、同じ分析からu が 最小値を取らないことが示されます。 Δ あなた ≤ 0 {\displaystyle \Delta u\leq 0} d あなた = 0 {\displaystyle du=0} Δ あなた = | d あなた | 2 {\displaystyle \Delta u=|du|^{2}} Δ あなた = | d あなた | 2 − 2 、 {\displaystyle \Delta u=|du|^{2}-2,}
このような解析の可能性は偏微分方程式に限らない。例えば、 あなた : M → R {\displaystyle u:M\to \mathbb {R} }
Δ あなた − | d あなた | 4 = ∫ M e あなた ( × ) d × 、 {\displaystyle \Delta u-|du|^{4}=\int _{M}e^{u(x)}\,dx,} これは一種の「非局所的」な微分方程式であるため、右辺の自動的な厳密な正値は、上記と同じ分析により、u が 最大値に到達できないことを示します。
この種の解析の適用範囲を様々な方法で拡張する方法は数多く存在する。例えば、u が 調和関数である場合、点p の存在はどこでも要求と矛盾しないため、上述のような矛盾は直接的には生じない。しかし、任意の実数s に対して、次のように定義される 関数u s を考えることもできる。Δ あなた ( p ) ≤ 0 {\displaystyle \Delta u(p)\leq 0} Δ あなた = 0 {\displaystyle \Delta u=0}
あなた s ( × ) = あなた ( × ) + s e × 1 。 {\displaystyle u_{s}(x)=u(x)+se^{x_{1}}.} それは一目瞭然だ
Δ あなた s = s e × 1 。 {\displaystyle \Delta u_{s}=se^{x_{1}}.} 上記の分析によれば、 ならばu s は 最大値に到達できません。u も最大値に到達できないと結論付けるために、 s から 0 への極限を考えることもでき ます。しかし、最大値を持たない関数列の点ごとの極限が最大値を持つ可能性はあります。と は いえ、M が境界を持ち、その境界と共にM がコンパクトである場合、 u を その境界まで連続的に拡張できると仮定すると、u とu s の 両方がで最大値に到達することが直ちに示されます。M上 の関数としてのu s は 最大値を持たないことを示したので、任意のsに対して u s の最大値は 上にあることが示されます。の逐次コンパクト性から、 u の最大値は で達成されます。これは調和関数の弱い最大値原理です。これだけでは、 uの最大値も M 上のどこかで達成される可能性を排除するものではありません。これが「強い最大値原理」の内容であり、さらなる分析が必要です。 s > 0 {\displaystyle s>0} M ∪ ∂ M 。 {\displaystyle M\cup \partial M.} ∂ M 。 {\displaystyle \partial M.} ∂ M 、 {\displaystyle \partial M,} ∂ M 。 {\displaystyle \partial M.}
上記の特定の関数の使用は、本質的にはそれほど重要ではありません。重要なのは、境界まで連続的に拡張され、ラプラシアンが厳密に正となる関数を持つことです。したがって、例えば、 e × 1 {\displaystyle e^{x_{1}}}
あなた s ( × ) = あなた ( × ) + s | × | 2 {\displaystyle u_{s}(x)=u(x)+s|x|^{2}} 同じ効果があります。
線形楕円偏微分方程式の古典的な強最大原理
証明の要約 M をユークリッド空間の開集合とする。 を二回微分可能な関数とし、その最大値C を得るものとする。 あなた : M → R {\displaystyle u:M\to \mathbb {R} }
1つの 私 j ∂ 2 あなた ∂ × 私 ∂ × j + b 私 ∂ あなた ∂ × 私 ≥ 0。 {\displaystyle a_{ij}{\frac {\partial ^{2}u}{\partial x^{i}\,\partial x^{j}}}+b_{i}{\frac {\partial u}{\partial x^{i}}}\geq 0.} 以下のものを発見(または存在を証明)できると仮定します。
M のコンパクト部分集合Ω で、内部が空でなく、Ω の内部にあるすべてのxに対して u ( x ) < C が成り立ち、 Ω の境界上にu ( x 0 ) = C となるx 0 が存在する。Ω の内部で2回微分可能で、h : Ω → R {\displaystyle h:\Omega\to\mathbb{R}} 1つの 私 j ∂ 2 h ∂ × 私 ∂ × j + b 私 ∂ h ∂ × 私 ≥ 0 、 {\displaystyle a_{ij}{\frac {\partial ^{2}h}{\partial x^{i}\,\partial x^{j}}}+b_{i}{\frac {\partial h}{\partial x^{i}}}\geq 0,} そして、 Ω の境界上でu + h ≤ Cが 成り立ち、h ( x 0 ) = 0となる。 すると、Ω の境界上でu + h − C ≤ 0 となる Ω 上でL ( u + h − C ) ≥ 0 となる。弱最大値原理によれば、 Ω 上でu + h − C ≤ 0 となる。これは次のように書き直すことができる。
− あなた ( × ) − あなた ( × 0 ) | × − × 0 | ≥ h ( × ) − h ( × 0 ) | × − × 0 | {\displaystyle -{\frac {u(x)-u(x_{0})}{|x-x_{0}|}}\geq {\frac {h(x)-h(x_{0})}{|x-x_{0}|}}} Ω のすべてのx について。右辺が明らかに正の性質を持つようにh を選択できる場合、 x 0 がM 上のu の最大値であるという事実と矛盾し、その勾配は必ずゼロになる。
証拠 上記の「プログラム」は実行可能である。Ωを 球面環状体とし、その中心x c を 閉集合u −1 ( C ) に近い点、つまり閉集合∂ M に近い点とする。そして、外半径Rをこの中心から u −1 ( C ) までの距離とする。x 0 を 、この距離を実現する後者の集合上の点とする。内半径ρ は任意である。定義する。
h ( × ) = ε ( e − α | × − × c | 2 − e − α R 2 ) 。 {\displaystyle h(x)=\varepsilon {\Big (}e^{-\alpha |x-x_{\text{c}}|^{2}}-e^{-\alpha R^{2}}{\Big )}.} ここで、 Ω の境界は2 つの球面から構成されます。外側の球面では、h = 0です。 R の選択により、この球面ではu ≤ C となり、したがって境界のこの部分ではu + h − C ≤ 0が成立し、さらにh ( x 0 ) = 0 という要件も満たされます。内側の球面では、u < C です。 u の連続性と内側の球面のコンパクトさにより、 u + δ < C となるようにδ > 0 を選択できます。この内側の球面ではh は定数であるため、内側の球面、したがってΩ の境界全体でu + h ≤ C となるようにε > 0 を選択できます。
直接計算すると
∑ 私 = 1 n ∑ j = 1 n 1つの 私 j ∂ 2 h ∂ × 私 ∂ × j + ∑ 私 = 1 n b 私 ∂ h ∂ × 私 = ε α e − α | × − × c | 2 ( 4 α ∑ 私 = 1 n ∑ j = 1 n 1つの 私 j ( × ) ( × 私 − × c 私 ) ( × j − × c j ) − 2 ∑ 私 = 1 n 1つの 私 私 − 2 ∑ 私 = 1 n b 私 ( × 私 − × c 私 ) ) 。 {\displaystyle \sum _{i=1}^{n}\sum _{j=1}^{n}a_{ij}{\frac {\partial ^{2}h}{\partial x^{i}\,\partial x^{j}}}+\sum _{i=1}^{n}b_{i}{\frac {\partial h}{\partial x^{i}}}=\varepsilon \alpha e^{-\alpha |x-x_{\text{c}}|^{2}}\left(4\alpha \sum _{i=1}^{n}\sum _{j=1}^{n}a_{ij}(x){\big (}x^{i}-x_{\text{c}}^{i}{\big )}{\big (}x^{j}-x_{\text{c}}^{j}{\big )}-2\sum _{i=1}^{n}a_{ii}-2\sum _{i=1}^{n}b_{i}{\big (}x^{i}-x_{\text{c}}^{i}{\big )}\right).} 右辺が非負であることが保証される条件はいくつかあります。以下の定理の説明を参照してください。
最後に、円環の内向きの放射状線に沿ったx 0 におけるh の方向微分は正であることに注意する。上記の要約で述べたように、これは x 0 におけるu の方向微分が非ゼロとなることを保証するが、これはx 0 が開集合M上の u の極大点であることと矛盾する。
定理の記述 以下は、ホップ (1927) の元の記述に従って、モーリーとスモラーの著書に記載されている定理の記述です。
M をユークリッド空間 ℝ n の開部分集合とする。1からn までの各i とj に対し、a ij とb i をM 上の連続関数とし、 a ij = a ji とする。Mの任意の x に対して、対称行列[ a ij ] は正定値であるとする。uが M 上の非定数C 2 関数で、
∑ 私 = 1 n ∑ j = 1 n 1つの 私 j ∂ 2 あなた ∂ × 私 ∂ × j + ∑ 私 = 1 n b 私 ∂ あなた ∂ × 私 ≥ 0 {\displaystyle \sum _{i=1}^{n}\sum _{j=1}^{n}a_{ij}{\frac {\partial ^{2}u}{\partial x^{i}\,\partial x^{j}}}+\sum _{i=1}^{n}b_{i}{\frac {\partial u}{\partial x^{i}}}\geq 0} M 上では、u は M 上で最大値に達しません。
連続性仮定の要点は、連続関数がコンパクト集合上で有界であるという点であり、ここでの関連するコンパクト集合とは、証明に現れる球面環状体である。さらに、同じ原理により、環状体内のすべてのx に対して、行列[ a ij ( x )] のすべての固有値が λ 以上となるような数 λ が存在する。 この場合、証明に現れるαは、これらの境界に対して大きいと仮定する。Evans の著書では、やや弱い定式化が用いられており、 M 内のすべてのx に対して[ a ij ] の固有値の下限となる正の数λ が存在すると仮定している。
これらの連続性の仮定は、証明が成立するために考えられる最も一般的な仮定ではないことは明らかです。例えば、以下はギルバーグとトゥルーディンガーによる定理の記述であり、同じ証明に続いています。
M を ユークリッド空間 ℝ n の開集合とする。1からn までの各i とj に対し、a ij とb i をa ij = a ji を満たすM 上の関数とする。Mのすべてのx に対し、対称行列[ a ij ] は正定値であるとし、その最小の固有値をλ(x) で表す。1からn までの各i に対し、 と はM 上の有界関数であるとする。uが M 上の非定数C 2 関数で 、1つの 私 私 λ {\displaystyle \textstyle {\frac {a_{ii}}{\lambda }}} | b 私 | λ {\displaystyle \textstyle {\frac {|b_{i}|}{\lambda }}}
∑ i = 1 n ∑ j = 1 n a i j ∂ 2 u ∂ x i ∂ x j + ∑ i = 1 n b i ∂ u ∂ x i ≥ 0 {\displaystyle \sum _{i=1}^{n}\sum _{j=1}^{n}a_{ij}{\frac {\partial ^{2}u}{\partial x^{i}\,\partial x^{j}}}+\sum _{i=1}^{n}b_{i}{\frac {\partial u}{\partial x^{i}}}\geq 0} M 上では、u は M 上で最大値に達しません。
これらの記述を、1次元の場合で既に見られたように、一般的な2階線形楕円方程式に単純に拡張することはできません。例えば、常微分方程式y ″ + 2 y = 0 は正弦波状の解を持ち、それらは必ず内部最大値を持ちます。これは高次元の場合にも当てはまり、高次元の場合、「固有関数」方程式Δ u + cu = 0の解も内部最大値を持つことがよくあります 。c の符号は、1次元の場合にも見られるように、重要です。例えば、 y ″ - 2 y = 0 の解は指数関数であり、そのような関数の最大値の性質は正弦関数のそれとは全く異なります。
参照
注記
参考文献
研究論文 Calabi, E. E. Hopf の最大原理の拡張とリーマン幾何学への応用 Duke Math. J. 25 (1958), 45–56.Cheng, SY; Yau, ST. リーマン多様体上の微分方程式とその幾何学的応用. Comm. Pure Appl. Math. 28 (1975), no. 3, 333–354. Gidas, B.; Ni, Wei Ming; Nirenberg, L. 最大原理による対称性と関連特性. Comm. Math. Phys. 68 (1979), no. 3, 209–243. Gidas, B.; Ni, Wei Ming; Nirenberg, L. R n における非線形楕円方程式の正値解の対称性. 数学的解析と応用, パートA, pp. 369–402, Adv. in Math. Suppl. Stud., 7a, Academic Press, New York-London, 1981.ハミルトン、リチャードS. 正曲率作用素を持つ4次元多様体. J. 微分幾何学24 (1986), no. 2, 153–179. ホップ、エバーハルト 。 Elementare Bemerkungen Über die Lösungen Partieller Differentialgleichungen zweiter Ordnung vom elliptischen Typus。シトバー。プロイス。アカド。ウィス。ベルリン 19 (1927)、147–152。ホップ、エバーハルト. 2階線形楕円微分方程式に関する考察. アメリカ数学会誌3 (1952), 791–793. ニーレンバーグ、ルイス「放物型方程式の強最大原理」純粋応用数学通信誌6(1953年)、167-177頁。 大森秀樹. リーマン多様体の等長浸漬. 日本数学会誌 19 (1967), 205–214. ヤウ、シン・トン. 完全リーマン多様体上の調和関数. 純粋応用数学通信誌 28 (1975), 201–228. クレイバーグ、HJA 経済プロセスにおける最適制御の最大原理について、1969 年 (ノースカロライナ州トロンハイム、ソシアロコノミスク研究所https://www.worldcat.org/title/on-the-maximum-principle-of-optimal-control-in-economic-processes/oclc/23714026 )
教科書 カファレッリ、ルイス・A. 、ザビエル・カブレ(1995年)『完全非線形楕円方程式 』プロビデンス、ロードアイランド:アメリカ数学会、pp. 31– 41. ISBN 0-8218-0437-5 。エヴァンス、ローレンス・C. 偏微分方程式. 第2版. 大学院数学研究科, 19. アメリカ数学会, プロビデンス, ロードアイランド, 2010. xxii+749 pp. ISBN 978-0-8218-4974-3 フリードマン、アヴナー著『放物型偏微分方程式』Prentice-Hall, Inc., Englewood Cliffs, NJ 1964 xiv+347 pp. ギルバーグ, デイヴィッド; トゥルディンガー, ニール・S. 2階楕円偏微分方程式. 1998年版の再版. 数学の古典. シュプリンガー・フェアラーク, ベルリン, 2001年. xiv+517 pp. ISBN 3-540-41160-7 Ladyženskaja, OA ; Solonnikov, VA; Uralʹceva, NN 「放物型線形方程式および準線形方程式」ロシア語からの翻訳。S. Smith著。Translations of Mathematical Monographs, Vol. 23 American Mathematical Society, Providence, RI 1968 xi+648 pp.Ladyzhenskaya, Olga A. ; Ural'tseva, Nina N. 線形および準線形楕円方程式。ロシア語からの翻訳はScripta Technica, Inc.。翻訳編集者:Leon Ehrenpreis。Academic Press, New York-London 1968 xviii+495 pp.リーバーマン、ゲイリー・M.『2階放物型微分方程式』World Scientific Publishing Co., Inc., River Edge, NJ, 1996. xii+439 pp. ISBN 981-02-2883-X モーリー、チャールズ・B・ジュニア著『変分法における多重積分』1966年版復刻版。数学の古典。シュプリンガー・フェアラーク、ベルリン、2008年。506頁。ISBN 978-3-540-69915-6 プロッター、マレー・H.; ワインバーガー、ハンス・F. 『微分方程式における最大原理』 1967年初版の訂正再版。シュプリンガー・フェアラーク、ニューヨーク、1984年。261頁。ISBN 0-387-96068-6 ロッカフェラー, RT (1970).凸解析 . プリンストン: プリンストン大学出版局.スモーラー、ジョエル。衝撃波と反応拡散方程式。第 2 版。 Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]、258。Springer-Verlag、ニューヨーク、1994。xxiv+632 pp. ISBN 0-387-94259-9