偏微分方程式 (へんびぶんほうていしき、英 :partial differential equation,PDE )は、未知関数の偏導関数 を含む微分方程式 である。偏微分方程式の厳密解を明示的な公式として得ることは、限られた場合を除いて困難である。そのため、解の存在性や一意性、正則性といった定性的な性質を研究する理論的アプローチ、近似的な解析解を求める摂動論や漸近解析、そしてコンピュータを使用して解を数値的に近似する 手法が、それぞれ重要な役割を果たしている。偏微分方程式は純粋数学の研究 においても大きな分野を占めており、そこでは存在性、一意性、正則性、安定性といった、様々な偏微分方程式の解の定性的な特徴に焦点が当てられている[ 1] 。数多くの未解決問題の中には、2000年にミレニアム懸賞問題 の一つとして挙げられた、ナビエ–ストークス方程式 の解の存在と滑らかさ がある。
偏微分方程式は、物理学 や工学 といった数学指向の科学分野で広く現れる。例えば、音 、熱 、拡散 、静電気学 、電磁気学 、熱力学 、流体力学 、弾性 、一般相対性理論 、量子力学 (シュレーディンガー方程式 、パウリ方程式 など)の現代的な科学的理解において、偏微分方程式は基礎となっている。また、微分幾何学 や変分法 などの純粋な数学的考察からも発生する。特筆すべき応用として、幾何学的トポロジー におけるポアンカレ予想 の証明において基本的な道具となったことが挙げられる。現代では、画像処理 、コンピュータビジョン 、機械学習 、金融工学 などの分野でも偏微分方程式が重要な役割を果たしている。
このような発生源の多様性もあり、偏微分方程式の種類は広範囲にわたる。個々の方程式を扱うために、多くの異なる手法が開発されてきた。そのため、偏微分方程式の「万能な理論」は存在せず、専門知識はいくつかの異なる下位分野に分かれている[ 2] 。
常微分方程式 は、一変数関数 に対応する偏微分方程式のサブクラスと見なすことができる。確率偏微分方程式 (英語版 ) や非局所方程式 は、「PDE」の概念を拡張したものであり、広く研究されている。より古典的なトピックとして、現在も活発な研究が行われているものには、楕円型偏微分方程式 、放物型偏微分方程式 、流体力学 におけるナビエ–ストークス方程式 、ボルツマン方程式 、分散型偏微分方程式 (英語版 ) などがある[ 3] 。
偏微分方程式 (PDE) の歴史は、物理現象を数学的に記述しようとする試みと、それに伴う解析学の基礎概念(関数、微分、積分)の確立と密接に関わって発展してきた。
偏微分方程式の研究は、1740年代にニュートン力学 を連続体へと拡張する過程で本格化した。中心となったのは「弦の振動」の問題である。ジャン・ル・ロン・ダランベール は1747年に一次元波動方程式 を見出し、その一般解(ダランベールの解)を示した。直後にレオンハルト・オイラー やダニエル・ベルヌーイ もこの問題に取り組み、解の表現(進行波による表現か、三角関数の重ね合わせか)を巡って論争が起きた。この論争は、当時の「関数」の定義(滑らかな曲線のみを指すのか、角を持つ折れ線も許容するのか)を再考させる契機となった。
また、流体力学 や重力ポテンシャル の研究において、ピエール=シモン・ラプラス はラプラス方程式 を定式化し、これが楕円型偏微分方程式 の原型となった。ジョゼフ=ルイ・ラグランジュ らは一階偏微分方程式の一般論を整備した。
19世紀に入ると、ジョゼフ・フーリエ による熱伝導方程式 の研究(1822年)が大きな転換点をもたらした。フーリエは、任意の関数を三角関数の級数(フーリエ級数 )で表現できると主張し、これにより偏微分方程式を変数分離法で解く道が開かれた。これは同時に、不連続な関数を含む広範な対象を扱うための厳密な解析学(イプシロン-デルタ論法 など)の整備を、ベルンハルト・リーマン やカール・ワイエルシュトラス らに促すこととなった。
19世紀後半には、オーギュスタン=ルイ・コーシー とソフィア・コワレフスカヤ によって解の局所的な存在定理(コーシー=コワレフスカヤの定理 )が証明されたほか、ジョージ・グリーン 、ガウス 、ディリクレ らによって、境界値問題 (ポテンシャル論)の手法が飛躍的に発展した。
20世紀初頭、ダフィット・ヒルベルト はヒルベルトの23の問題 において、偏微分方程式の解の滑らかさや存在条件を主要な未解決問題として挙げた。これに応える形で、偏微分方程式論は「解を見つける」ことから「関数空間の構造を調べる」方向へと大きくシフトした。
古典的な微分が定義できないような現象(衝撃波など)や、不連続な境界値を持つ問題を扱うために、解の概念そのものを拡張する必要が生じた。1930年代以降、セルゲイ・ソボレフ によるソボレフ空間 の導入や、ローラン・シュヴァルツ による超関数 (ディストリビューション)の理論により、「弱解 」という概念が確立された。これにより、偏微分方程式は関数解析学 の枠組みで統一的に扱われるようになり、現代の線形 および非線形偏微分方程式 (英語版 ) の理論的基礎が完成した。
偏微分方程式とは、n ≥ 2 {\displaystyle n\geq 2} 個の変数の未知関数とその偏微分を含む方程式である[ 4] 。方程式に含まれる導関数のうち、最高の階数をその偏微分方程式の階数 (order) と呼ぶ。例えば、二階導関数を含み、三階以上の導関数を含まないものは二階偏微分方程式となる。
R n {\displaystyle \mathbb {R} ^{n}} の開集合U {\displaystyle U} に属する変数x = ( x 1 , … , x n ) {\displaystyle x=(x_{1},\dots ,x_{n})} の未知関数u : U → R , {\displaystyle u:U\rightarrow \mathbb {R} ,} に対して、一般にk {\displaystyle k} 階の偏微分方程式は以下のように記述される。F ( x , u , D u , D 2 u , … , D k u ) = 0 , {\displaystyle F(x,u,Du,D^{2}u,\ldots ,D^{k}u)=0,} ここで、F {\displaystyle F} は適当な領域からR {\displaystyle \mathbb {R} } への関数であり、u {\displaystyle u} のk {\displaystyle k} 階までの全ての偏微分、u {\displaystyle u} 自身、および独立変数x {\displaystyle x} を変数とする。D j u {\displaystyle D^{j}u} はj {\displaystyle j} 階の偏微分の全体を表す。
未知関数がベクトル値をとる場合、すなわちu : U → R m {\displaystyle u:U\rightarrow \mathbb {R} ^{m}} の場合には、偏微分方程式は方程式系となり、F {\displaystyle F} はR m {\displaystyle \mathbb {R} ^{m}} 値関数となる。このような方程式系は流体力学 や電磁気学 などで自然に現れる。
PDEを書く際、偏微分を下付き文字で表すことが一般的である。例えば:u x = ∂ u ∂ x , u x x = ∂ 2 u ∂ x 2 , u x y = ∂ 2 u ∂ y ∂ x = ∂ ∂ y ( ∂ u ∂ x ) . {\displaystyle u_{x}={\frac {\partial u}{\partial x}},\quad u_{xx}={\frac {\partial ^{2}u}{\partial x^{2}}},\quad u_{xy}={\frac {\partial ^{2}u}{\partial y\partial x}}={\frac {\partial }{\partial y}}\left({\frac {\partial u}{\partial x}}\right).} u がn 変数の関数である一般的な状況では、u i またはu xi はi 番目の独立変数x i に関する一階偏微分を表し、u ij は二階偏微分を表す。なお、未知関数がベクトル値をとる場合の成分表示と混同を避けるため、添字の前にカンマを打つ、あるいは変数を明示する記法も用いられる。
ギリシャ文字Δ はラプラス作用素 を表す。もしu がn 変数の関数であれば、以下のようになる。Δ u = u 11 + u 22 + ⋯ + u n n . {\displaystyle \Delta u=u_{11}+u_{22}+\cdots +u_{nn}.} 物理学の文献では、ラプラス作用素はしばしば∇ 2 {\displaystyle \nabla ^{2}} で表される。数学の文献でもこの記法は用いられるが、多変数関数の二階導関数全体を指すヘッセ行列 (D 2 u {\displaystyle D^{2}u} またはH ( u ) {\displaystyle \mathbf {H} (u)} と記される)との混同を避けるため、Δ {\displaystyle \Delta } が好まれる傾向にある。
PDEは、未知数とその微分に関して線形である場合、線形 と呼ばれる。例えば、x とy の関数u に対して、二階線形PDEは以下の形式をとる。a 1 ( x , y ) u x x + a 2 ( x , y ) u x y + a 3 ( x , y ) u y x + a 4 ( x , y ) u y y + a 5 ( x , y ) u x + a 6 ( x , y ) u y + a 7 ( x , y ) u = f ( x , y ) {\displaystyle a_{1}(x,y)u_{xx}+a_{2}(x,y)u_{xy}+a_{3}(x,y)u_{yx}+a_{4}(x,y)u_{yy}+a_{5}(x,y)u_{x}+a_{6}(x,y)u_{y}+a_{7}(x,y)u=f(x,y)} ここでai とf は独立変数x とy のみの関数である。(多くの場合、混合偏微分uxy とuyx は等置されるが、これは線形性の議論には必須ではない。)もしai が定数(x とy に依存しない)であれば、そのPDEは定数係数線形 と呼ばれる。もしf がどこでもゼロであれば、その線形PDEは斉次 であり、そうでなければ非斉次 である。(これは、PDEの解に対する高周波振動の影響を研究する漸近的均質化 (英語版 ) とは別のものである。)
線形PDEに最も近いのが半線形 (semi-linear) PDEであり、最高階微分のみが線形項として現れ、その係数は独立変数の関数である。低階の微分と未知関数は任意に現れてよい。例えば、2変数の一般的な二階半線形PDEは以下のようになる。a 1 ( x , y ) u x x + a 2 ( x , y ) u x y + a 3 ( x , y ) u y x + a 4 ( x , y ) u y y + f ( u x , u y , u , x , y ) = 0 {\displaystyle a_{1}(x,y)u_{xx}+a_{2}(x,y)u_{xy}+a_{3}(x,y)u_{yx}+a_{4}(x,y)u_{yy}+f(u_{x},u_{y},u,x,y)=0}
準線形 (quasilinear) PDEでは、同様に最高階微分は線形項としてのみ現れるが、係数は未知数および低階微分の関数であってもよい。a 1 ( u x , u y , u , x , y ) u x x + a 2 ( u x , u y , u , x , y ) u x y + a 3 ( u x , u y , u , x , y ) u y x + a 4 ( u x , u y , u , x , y ) u y y + f ( u x , u y , u , x , y ) = 0 {\displaystyle a_{1}(u_{x},u_{y},u,x,y)u_{xx}+a_{2}(u_{x},u_{y},u,x,y)u_{xy}+a_{3}(u_{x},u_{y},u,x,y)u_{yx}+a_{4}(u_{x},u_{y},u,x,y)u_{yy}+f(u_{x},u_{y},u,x,y)=0} 物理学における基本的なPDEの多くは準線形である。例えば、一般相対性理論 のアインシュタイン方程式 や、流体運動を記述するナビエ–ストークス方程式 などが挙げられる。
線形性の性質を持たないPDEは完全非線形 (fully nonlinear) と呼ばれ、最高階微分の1つ以上に非線形性を持つ。例として、微分幾何学 で現れるモンジュ=アンペール方程式 (英語版 ) がある[ 5] 。
一般に非線形な偏微分方程式として重要なものに、
などがある。
uxy =uyx を仮定すると、2つの独立変数を持つ一般的な線形二階偏微分方程式は以下の形式を持つ。
A u x x + 2 B u x y + C u y y + ⋯ = 0 , {\displaystyle Au_{xx}+2Bu_{xy}+Cu_{yy}+\cdots =0,} ここで係数A ,B ,C ... はx とy に依存してもよい。xy 平面のある領域上でA 2 +B 2 +C 2 > 0 であれば、PDEはその領域で二階である。この形式は円錐曲線の式に類似している。A x 2 + 2 B x y + C y 2 + ⋯ = 0. {\displaystyle Ax^{2}+2Bxy+Cy^{2}+\cdots =0.}
より正確には、∂x をX に置き換え、他の変数についても同様にすることで(形式的にはフーリエ変換 によって行われる)、定数係数PDEを同次数の多項式に変換でき、最高次数の項(斉次多項式 、ここでは二次形式 )が分類にとって最も重要となる。
円錐曲線 と二次形式を判別式 B 2 − 4AC に基づいて放物型、双曲型、楕円型に分類するのと同様に、与えられた点における二階PDEについても同じことができる。ただし、PDEの判別式 は、xy の項がB ではなく2B という慣例のため、B 2 −AC で与えられる。形式的には、(関連する二次形式の)判別式は(2B )2 − 4AC = 4(B 2 −AC ) であるが、単純化のために係数4は落とされる。
B 2 −AC < 0 (楕円型偏微分方程式 ): 楕円型PDEの解は、方程式と解が定義されている領域の内部において、係数が許す限り滑らかである。代表的な例としてラプラス方程式 がある。例えば、ラプラス方程式 の解は定義されている領域内で解析的であるが、境界値は滑らかでない場合がある。物理的な例として、定常な亜音速流の支配方程式は楕円型となる。例えば、オイラー=トリコミ方程式 (英語版 ) (u x x = x u y y {\displaystyle u_{xx}=xu_{yy}} )は、x < 0 の領域において楕円型である。定数係数の楕円型方程式の場合、適切な変数変換により、方程式を次の標準形で表現できる:u x x + u y y + ⋯ = 0 , {\displaystyle u_{xx}+u_{yy}+\cdots =0,} ここで x と y は変換後の変数に対応する。これは、ラプラス方程式 がこのタイプの例であることを裏付けている[ 10] 。B 2 −AC = 0 (放物型偏微分方程式 ): すべての点で放物型である方程式は、独立変数の変換によって熱方程式 に類似した形式に変換できる。解は変換された時間変数の増加とともに滑らかになる。オイラー=トリコミ方程式はx = 0 の直線上で放物型である。定数係数の放物型方程式の場合、適切な変数変換により、方程式を次の標準形で表現できる:u t − u x x + ⋯ = 0 , {\displaystyle u_{t}-u_{xx}+\cdots =0,} ここで t は時間変数、x は空間変数に対応する。これは熱方程式 がこのタイプの例であることを示している[ 10] 。B 2 −AC > 0 (双曲型偏微分方程式 ): 双曲型方程式は、初期データにおける関数や微分の不連続性を保持する。例として波動方程式 がある。超音速での流体の運動は双曲型PDEで近似でき、オイラー=トリコミ方程式はx > 0 において双曲型である。変数変換により、方程式は常に次の形式で表現できる:u x x − u y y + ⋯ = 0 , {\displaystyle u_{xx}-u_{yy}+\cdots =0,} ここで x と y は変換後の変数に対応する。これは、波動方程式 がこのタイプの例であることを正当化する[ 10] 。n 個の独立変数x 1 ,x 2 , …,x n がある場合、一般的な二階線形偏微分方程式は以下の形式を持つ。L u = ∑ i = 1 n ∑ j = 1 n a i , j ∂ 2 u ∂ x i ∂ x j + ∑ i = 1 n b i ∂ u ∂ x i + c u = 0. {\displaystyle Lu=\sum _{i=1}^{n}\sum _{j=1}^{n}a_{i,j}{\frac {\partial ^{2}u}{\partial x_{i}\partial x_{j}}}+\sum _{i=1}^{n}b_{i}{\frac {\partial u}{\partial x_{i}}}+cu=0.}
分類は係数行列a i ,j の固有値 の符号(シグネチャ)に依存する。
楕円型: 固有値がすべて正、またはすべて負。 放物型: 固有値がすべて正、またはすべて負だが、1つだけゼロがある。 双曲型: 負の固有値が1つだけで残りがすべて正、または正の固有値が1つだけで残りがすべて負。 超双曲型 (Ultrahyperbolic): 正の固有値が複数かつ負の固有値が複数あり、ゼロ固有値がない[ 11] 。 楕円型、放物型、双曲型方程式の理論は、主にラプラス方程式 、熱方程式 、波動方程式 という標準的な例を中心に、あるいはそれらに基づいて、何世紀にもわたって研究されてきた。
二階偏微分方程式の楕円型・放物型・双曲型の分類は、PDEの理論と応用において極めて重要である。この分類は、解の定性的な性質、適切な初期条件 および境界条件 の設定、数値解法の選択、そして解の滑らかさ や伝播特性に関する指針を与える。楕円型方程式は通常境界値問題として定式化され、放物型および双曲型方程式は初期値問題または初期境界値問題として扱われる。また、この分類は解の情報伝播の仕方とも密接に関連しており、双曲型では有限速度での伝播、放物型では無限速度での伝播、楕円型では瞬時の影響が特徴的である。
しかし、この分類は二階項の線形性にのみ依存するため、半線形および準線形PDEにも適用可能である。基本的なタイプは、領域によって楕円型から双曲型に変化するオイラー=トリコミ方程式 (英語版 ) のような混合型や、高階PDEにも拡張される。
これらの方程式の型は、適切な問題設定(境界条件や初期条件)と密接に関係している。一般に、楕円型方程式は境界値問題として定式化され、領域の境界上で未知関数や微分に関する条件を与えることで領域内部の解が一意に定まることが多い。一方、放物型および双曲型方程式は、時間発展を記述するため、初期値問題または初期境界値問題として扱われるのが通常である。
偏微分方程式の分類は一階方程式系にも拡張できる。ここで未知数u はm 個の成分を持つベクトル であり、係数行列Aν はν = 1, 2, …,n に対してm ×m 行列である。偏微分方程式は以下の形式をとる。L u = ∑ ν = 1 n A ν ∂ u ∂ x ν + B = 0 , {\displaystyle Lu=\sum _{\nu =1}^{n}A_{\nu }{\frac {\partial u}{\partial x_{\nu }}}+B=0,} ここで係数行列Aν とベクトルB はx とu に依存してもよい。もし超曲面 S が陰関数形式φ ( x 1 , x 2 , … , x n ) = 0 , {\displaystyle \varphi (x_{1},x_{2},\ldots ,x_{n})=0,} で与えられ、φ が非ゼロの勾配を持つとき、特性形式 (characteristic form) が消えるならば、S は与えられた点における作用素 L の特性曲面 (characteristic surface) である。Q ( ∂ φ ∂ x 1 , … , ∂ φ ∂ x n ) = det [ ∑ ν = 1 n A ν ∂ φ ∂ x ν ] = 0. {\displaystyle Q\left({\frac {\partial \varphi }{\partial x_{1}}},\ldots ,{\frac {\partial \varphi }{\partial x_{n}}}\right)=\det \left[\sum _{\nu =1}^{n}A_{\nu }{\frac {\partial \varphi }{\partial x_{\nu }}}\right]=0.}
この条件の幾何学的解釈は以下の通りである。S 上でu のデータが規定されている場合、微分方程式からS 上のu の法線微分を決定できる可能性がある。もしS 上のデータと微分方程式がS 上のu の法線微分を決定する場合、S は非特性的である。もしS 上のデータと微分方程式がS 上のu の法線微分を決定しない 場合、その曲面は特性的 であり、微分方程式はS 上のデータを制限する。すなわち、微分方程式はS に「内在的」である。
一階方程式系Lu = 0 は、L に対して特性的な曲面が存在しない場合、楕円型 である。S 上のu の値と微分方程式は、常にS 上のu の法線微分を決定する。 一階方程式系は、ある点において法線ξ を持つ空間的 (spacelike) 曲面S が存在する場合、その点で双曲型 である。これは、ξ に直交する任意の非自明なベクトルη とスカラー乗数λ に対して、方程式Q (λξ +η ) = 0 がm 個の実根λ 1 ,λ 2 , …,λ m を持つことを意味する。これらの根が常に異なる場合、システムは狭義双曲型 である。この条件は幾何学的には以下のように解釈される。特性形式Q (ζ ) = 0 は、ζ 空間において原点を頂点とする錐体(法線錐)を定義する。双曲型の場合、この法線錐は実数領域に非自明な解を持ち、時間的な方向軸はこの錐体の内側を通る。一方、楕円型の場合、特性形式Q (ζ ) = 0 を満たす実ベクトルζ ≠ 0 は存在せず(法線錐が一点に縮退する、あるいは虚の方向にしか存在しない)、したがって実特性曲面を持たない。 偏微分方程式の一般解は通常、任意関数を含むため、特定の物理現象を一意に記述するには、方程式本体に加えて適切な付加条件が必要となる。これらの方程式と条件の組み合わせが、解の存在性と一意性が保証された良設定問題 を構成する。
時間発展を含む問題(放物型や双曲型)において、解析を開始する時刻(通常は t = 0)における未知関数の状態を指定する条件。
解析領域の境界における振る舞いを規定する条件であり、主に以下の3種類がある。
ディリクレ境界条件 :境界上での未知関数の値を直接指定する(例:境界の温度を固定する)。ノイマン境界条件 :境界上での法線方向の微分値を指定する(例:境界を通る熱流束や流体の流入量を指定する)。ロビン境界条件 :未知関数とその微分の線形結合を指定する(例:境界での熱伝達・対流を考慮する)。純粋数学 において、PDEの理論的研究は、解の公式を見つけることよりも解が存在するための基準や解の性質に焦点を当てられる。
良設定性 (well-posedness) とは、アダマール によって提唱された、偏微分方程式の問題設定が物理的・数学的に適切であることを示す概念である。問題が良設定であるためには、以下の条件を満たす必要がある。
解の存在性 :与えられたデータ(初期条件や境界条件)に対し、少なくとも一つの解が存在する。解の一意性 :解が存在する場合、それはただ一つである。解の連続依存性 (安定性):入力データを連続的に変化させたとき、対応する解も連続的に変化する。偏微分方程式においてこれらの条件を確認することは、常微分方程式 (ODE) よりもはるかに困難かつ重要である。ODEの一般解が有限個の任意定数(パラメータ)を含むのに対し、PDEの一般解は任意関数を含むためである。
例えば、2変数関数v ( x , y ) {\displaystyle v(x,y)} に対する単純な方程式v x y = 0 {\displaystyle v_{xy}=0} の一般解は、任意の1変数関数f , g {\displaystyle f,g} を用いてv ( x , y ) = f ( x ) + g ( y ) {\displaystyle v(x,y)=f(x)+g(y)} と表される。このように無数の解が存在するため、特定の解を定めるには定義域全体を含めた適切な条件設定が不可欠である。
方程式の型によって、解の一意性を保証するために必要な条件は異なる。以下にその著しい対比を示す。
さらに、非線形方程式ではより複雑な現象が起こりうる。例えば、微分幾何学 に由来するベルンシュタインの問題 (英語版 ) (極小曲面方程式)は、∂ ∂ x u x 1 + u x 2 + u y 2 + ∂ ∂ y u y 1 + u x 2 + u y 2 = 0 {\displaystyle {\frac {\partial }{\partial x}}{\frac {u_{x}}{\sqrt {1+u_{x}^{2}+u_{y}^{2}}}}+{\frac {\partial }{\partial y}}{\frac {u_{y}}{\sqrt {1+u_{x}^{2}+u_{y}^{2}}}}=0} という非線形方程式であるが、全平面R 2 {\displaystyle \mathbb {R} ^{2}} で定義された解は、u ( x , y ) = a x + b y + c {\displaystyle u(x,y)=ax+by+c} という1次式に限られる(ベルンシュタインの定理)。これは、任意関数を含む一般解を持つ線形方程式とは対照的な性質である。
弱解 (weak solution) は、PDEを満たすが、通常の古典的な意味(すべての微分が通常の意味で存在し連続である)とは異なる意味で満たす関数のことである。弱解の概念は、古典解が存在しない場合や、不連続性を持つ解を扱う際に特に重要となる。この用語の意味は文脈によって異なる場合があるが、最も一般的に使用される定義の一つは、超関数 (distribution) の概念とソボレフ空間 に基づいている。ソボレフ空間とは、関数とその弱微分が可積分であることを要求する関数空間であり、弱解の理論において中心的な役割を果たす。
弱解の定義の例[ 12] は以下の通りである。
次で与えられる境界値問題を考える。L u = f in U , u = 0 on ∂ U , {\displaystyle {\begin{aligned}Lu&=f\quad {\text{in }}U,\\u&=0\quad {\text{on }}\partial U,\end{aligned}}} ここでL u = − ∑ i , j ∂ j ( a i j ∂ i u ) + ∑ i b i ∂ i u + c u {\displaystyle Lu=-\sum _{i,j}\partial _{j}(a^{ij}\partial _{i}u)+\sum _{i}b^{i}\partial _{i}u+cu} は発散形式 の二階偏微分作用素を表す。
u ∈ H 0 1 ( U ) {\displaystyle u\in H_{0}^{1}(U)} が弱解であるとは、すべてのv ∈ H 0 1 ( U ) {\displaystyle v\in H_{0}^{1}(U)} に対して∫ U [ ∑ i , j a i j ( ∂ i u ) ( ∂ j v ) + ∑ i b i ( ∂ i u ) v + c u v ] d x = ∫ U f v d x {\displaystyle \int _{U}{\bigg [}\sum _{i,j}a^{ij}(\partial _{i}u)(\partial _{j}v)+\sum _{i}b^{i}(\partial _{i}u)v+cuv{\bigg ]}dx=\int _{U}fvdx} が成り立つことである。これは形式的な部分積分によって導出できる。
弱解の例は以下の通りである。ϕ ( x ) = 1 4 π 1 | x | {\displaystyle \phi (x)={\frac {1}{4\pi }}{\frac {1}{|x|}}} は、− ∇ 2 ϕ = δ in R 3 {\displaystyle -\nabla ^{2}\phi =\delta {\text{ in }}\mathbb {R} ^{3}} を超関数の意味で満たす弱解である。
なぜなら、形式的に∫ R 3 ( − ∇ 2 ϕ ( x ) ) ψ ( x ) d x = ∫ R 3 ϕ ( x ) ( − ∇ 2 ψ ( x ) ) d x = ψ ( 0 ) for ψ ∈ C c ∞ ( R 3 ) {\displaystyle \int _{\mathbb {R} ^{3}}(-\nabla ^{2}\phi (x))\psi (x)dx=\int _{\mathbb {R} ^{3}}\phi (x)(-\nabla ^{2}\psi (x))dx=\psi (0){\text{ for }}\psi \in C_{c}^{\infty }(\mathbb {R} ^{3})} となるからである。
ここでC c ∞ ( R 3 ) {\displaystyle C_{c}^{\infty }(R^{3})} は、コンパクトな台を持つ滑らかな関数の空間(テスト関数空間)を表す。
正則性 (regularity) は、弱解の可積分性と微分可能性を指し、多くの場合ソボレフ空間 によって表現できる。
この問題は、古典解を探すのが難しいために生じる。研究者は、最初は弱解を見つけ、その後、それが古典解としての資格を得るのに十分滑らかであるかどうかを確認する傾向がある。
この研究分野では、関数解析学 の結果が頻繁に使用される。
偏微分方程式の解の振る舞いを調べる上で、以下の概念が重要である。
最大値原理 :楕円型や放物型方程式の解が、領域の内部で最大値(または最小値)を取らず、必ず境界上で取るという性質。解の一意性や安定性の証明に用いられる。エネルギー法 :解の「エネルギー」(通常はL2ノルムなど)の物理的な保存則や減衰を数学的に定式化する手法。双曲型方程式の解の評価などに多用される。偏微分方程式の一般解は、常微分方程式とは異なり、任意定数ではなく任意関数を含むため、無数の解を持つ自由度の高いものである。したがって、物理的あるいは数学的に意味のある特定の解を一意に決定するためには、方程式に加えて適切な境界条件 や初期条件 を課して定式化する必要がある(良設定問題 )。
代表的な境界条件には、境界上で未知関数の値を指定するディリクレ境界条件 、境界上での法線方向微分(流束など)を指定するノイマン境界条件 、未知関数とその法線微分の線形結合を指定するロビン境界条件 などがある。
以下の節では、特別な記述がない限り、変数は時間 t と3次元空間 (x ,y ,z ) とする。ただし、これらの方程式は数学的には一般の次元に拡張できる。
楕円型偏微分方程式 の典型的な例として、
ψ x x + ψ y y + ψ z z = 0 {\displaystyle \psi _{xx}+\psi _{yy}+\psi _{zz}=0} で表されるラプラス方程式 がある。これはまた、∇(ナブラ )や Δ(ラプラス作用素 )といった微分作用素を用いて
∇ 2 ψ = 0 , Δ ψ = 0 {\displaystyle {\nabla }^{2}\psi =0,\quad \Delta \psi =0} のようにも書かれる。ラプラス方程式の解は調和関数 と呼ばれ、重力場や静電場 といった物理的なベクトル場 のポテンシャル を与える。
ラプラス方程式は既知の関数f (x ,y ,z ) に関する微分方程式
∇ 2 ψ = ψ x x + ψ y y + ψ z z = f ( x , y , z ) {\displaystyle \nabla ^{2}\psi =\psi _{xx}+\psi _{yy}+\psi _{zz}=f(x,y,z)} に一般化される。この偏微分方程式をポアソン方程式 という[ 13] 。これは質量の存在する重力場や、電荷の存在する静電場など、場に発生源がある場合のポテンシャルを記述する方程式である。
次の方程式のことをいう。電磁波 の放射 、地震学 、音響学 などで用いられる。
( ∇ 2 + k 2 ) ψ = 0 {\displaystyle (\nabla ^{2}+k^{2})\psi =0} 波動方程式 は時間変数t を含む双曲型偏微分方程式
ψ t t = c 2 ∇ 2 ψ = c 2 ( ψ x x + ψ y y + ψ z z ) {\displaystyle \psi _{tt}=c^{2}\nabla ^{2}\psi =c^{2}(\psi _{xx}+\psi _{yy}+\psi _{zz})} のことである[ 14] [ 15] 。この方程式は光波 や音波 といった波 を記述するもので、定数c は波の速さを示している。より身近な現象として、ひも の振動であるとか太鼓 の鼓面の振動などといったものもこの方程式に従う。波動方程式の解は基本的には正弦波 を重ね合わせることによって得られる。
移流 方程式 (en:Advection )は速度場u = (u ,v ,w )のもとでの保存スカラー量ψの輸送を記述するもので、方程式は
ψ t + ∇ ⋅ ( u ψ ) = ψ t + ( u ψ ) x + ( v ψ ) y + ( w ψ ) z = 0 {\displaystyle \psi _{t}+\nabla \cdot (\mathbf {u} \psi )=\psi _{t}+(u\psi )_{x}+(v\psi )_{y}+(w\psi )_{z}=0} であたえられる。もし速度場u が管状ベクトル場、すなわち ∇・u = 0 ならば方程式は
ψ t + u ⋅ ∇ ψ = ψ t + u ψ x + v ψ y + w ψ z = 0 {\displaystyle \psi _{t}+\mathbf {u} \cdot \nabla \psi =\psi _{t}+u\psi _{x}+v\psi _{y}+w\psi _{z}=0} と簡略化される。
一次元線形移流方程式
ψ t + u ψ x = 0 {\displaystyle \psi _{t}+u\psi _{x}=0} (u は定数)は、初期波形が速度u で平行移動する解を持つ基本的な方程式である。
次の方程式は、流体力学における非線形性と粘性の相互作用を記述する。
u t + u u x = ν u x x {\displaystyle u_{t}+uu_{x}=\nu u_{xx}} ここでν {\displaystyle \nu } は粘性係数である。左辺の非線形項(対流項)により衝撃波が形成され、右辺の拡散項により滑らかさが保たれる。非線形波動や乱流の基礎モデルとして重要である。
熱方程式 (または拡散方程式 )は、与えられた領域において時間とともに変化する熱や物質の拡散を記述する放物型偏微分方程式 の代表例である。
u t = k ∇ 2 ψ = k ( ψ x x + ψ y y + ψ z z ) {\displaystyle u_{t}=k\nabla ^{2}\psi =k(\psi _{xx}+\psi _{yy}+\psi _{zz})} によって与えられる。ここでu は温度場や濃度場を表し、定数k は熱拡散率や拡散係数を示す。解は時間の経過とともに滑らかになり、空間的に均一な状態へと近づく性質を持つ。定数k は物質の熱伝導性や拡散係数などを示している。解は時間の増加とともに大体均一に分布するように変化し、t →∞で調和関数に近づく。
量子力学 の基礎方程式であるシュレーディンガー方程式 も重要な偏微分方程式である[ 16] 。時間に依存するシュレーディンガー方程式i ℏ ∂ ψ ∂ t = − ℏ 2 2 m ∇ 2 ψ + V ψ {\displaystyle i\hbar {\frac {\partial \psi }{\partial t}}=-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi +V\psi } は、虚数単位i を含む複素係数の方程式であり、純粋な実係数の二階方程式とは性質が異なる。数学的には分散型偏微分方程式 (英語版 ) の代表的な例に位置づけられ、解は時間とともに波動的に振る舞い、ノルムを保存する(ユニタリ発展 )という性質を持つ。一方、時間に依存しない(定常状態の)シュレーディンガー方程式は楕円型偏微分方程式 である。
電磁気学 の基礎となる4つの方程式系である。真空中の電場E {\displaystyle \mathbf {E} } と磁場B {\displaystyle \mathbf {B} } に対して以下のように記述される。
{ ∇ ⋅ B ( t , r ) = 0 ∇ × E ( t , r ) = − ∂ B ( t , r ) ∂ t ∇ ⋅ D ( t , r ) = ρ ( t , r ) ∇ × H ( t , r ) = j ( t , r ) + ∂ D ( t , r ) ∂ t {\displaystyle {\begin{cases}{\begin{aligned}\nabla \cdot {\boldsymbol {B}}(t,{\boldsymbol {r}})&=0\\\nabla \times {\boldsymbol {E}}(t,{\boldsymbol {r}})&=-{\dfrac {\partial {\boldsymbol {B}}(t,{\boldsymbol {r}})}{\partial t}}\\\nabla \cdot {\boldsymbol {D}}(t,{\boldsymbol {r}})&=\rho (t,{\boldsymbol {r}})\\\nabla \times {\boldsymbol {H}}(t,{\boldsymbol {r}})&={\boldsymbol {j}}(t,{\boldsymbol {r}})+{\dfrac {\partial {\boldsymbol {D}}(t,{\boldsymbol {r}})}{\partial t}}\end{aligned}}\end{cases}}} これは一階の偏微分方程式系であり、電磁波の存在を予言する。
流体の運動を記述する非線形偏微分方程式である。非圧縮性粘性流体の場合は次のように書かれる。
ρ ( ∂ v ∂ t + ( v ⋅ ∇ ) v ) = − ∇ p + μ Δ v + f {\displaystyle \rho \left({\frac {\partial \mathbf {v} }{\partial t}}+(\mathbf {v} \cdot \nabla )\mathbf {v} \right)=-\nabla p+\mu \Delta \mathbf {v} +\mathbf {f} } この方程式は、左辺の非線形項(対流項)のために解析が極めて困難であり、3次元における解の存在と滑らかさはミレニアム懸賞問題 の一つとなっている。
一般相対性理論 において、時空の歪みと物質・エネルギーの分布を関連付ける方程式である。
G μ ν + Λ g μ ν = κ T μ ν {\displaystyle G_{\mu \nu }+\Lambda g_{\mu \nu }=\kappa T_{\mu \nu }} これは計量テンソルg μ ν {\displaystyle g_{\mu \nu }} に関する高度に非線形な二階偏微分方程式系であり、宇宙論 やブラックホールの研究の基礎となる。
金融工学 において、オプション価格の理論的な値を算出するための方程式である。
∂ V ∂ t + 1 2 σ 2 S 2 ∂ 2 V ∂ S 2 + r S ∂ V ∂ S − r V = 0 {\displaystyle {\frac {\partial V}{\partial t}}+{\frac {1}{2}}\sigma ^{2}S^{2}{\frac {\partial ^{2}V}{\partial S^{2}}}+rS{\frac {\partial V}{\partial S}}-rV=0} この方程式は適切な変数変換によって熱方程式 に帰着させることができ、偏微分方程式の理論が社会科学分野へ応用された顕著な例である。
化学物質の反応と拡散を記述する方程式系である。
∂ u ∂ t = D ∇ 2 u + R ( u ) {\displaystyle {\frac {\partial \mathbf {u} }{\partial t}}=D\nabla ^{2}\mathbf {u} +\mathbf {R} (\mathbf {u} )} 1952年にアラン・チューリング は、この方程式によって、生物の縞模様や斑点模様(チューリング・パターン )が自己組織的に形成されることを理論的に示した。
浅い水域を伝わる波などを記述する非線形偏微分方程式である。
u t + 6 u u x + u x x x = 0 {\displaystyle u_{t}+6uu_{x}+u_{xxx}=0} この方程式は、非線形効果と分散効果が釣り合うことで、形を崩さずに一定速度で進む孤立波ソリトン を解として持つ。ソリトン同士は衝突しても互いに形を変えずに通り抜けるという粒子のような性質を持ち、この発見は可積分系 という分野を切り拓いた。現代では光ファイバー 中の通信パルスの挙動(非線形シュレディンガー方程式による記述)などに応用されている。
以下の表に、代表的な偏微分方程式とその主な適用分野を示す。
線型偏微分方程式はその解を基底関数 系で線型展開したもので近似することにより一般的に解かれる事が多い(たとえば正弦波 関数を使ったフーリエ級数 展開[ 17] )。展開した個々の解の線型結合 もまた元の方程式の解になる。また線型偏微分方程式が変数分離可能な構造を持つ場合には変数分離法 により低い次元の微分方程式の問題に帰着して解くことができる場合がある。
非線形偏微分方程式に対しては、線形方程式のような統一的な解法理論は存在しない。しかし、方程式の種類に応じて有効な手法が多数開発されている。たとえば、ホモトピー原理 (英語版 ) は劣決定方程式系を解くための強力な方法であり、ソリトン 方程式に対しては広田の方法 が標準的な手法として用いられる[ 18] 。また、変分法 、リー群 の対称性、特性曲線法 なども、特定のクラスの非線形方程式に対して有効である。さらに、解の存在性や一意性、正則性に関する理論的研究も活発に行われている。
あるいは偏微分方程式を容易に解ける方程式から変化したものであるとみなせる場合にはその解を摂動 的な展開を用いて表すことで近似解が得られる場合がある。またそれ以外に数値的な近似解法として、有限差分スキーム [ 19] [ 20] や有限要素法 [ 21] [ 22] [ 23] などが挙げられる。 具体的な解が必要であるが解析的な手段を用いては解くことのできない多くの科学や工学上の問題は,このような近似手法により離散化された大規模な数値計算の問題に変換され,もっぱらコンピュータを用いた計算により解かれる[ 24] [ 25] [ 26] 。境界条件も含めて近似解を表すための自由度は非常に大きいものになることが普通であるので,計算を行うために扱うデータの量や演算の回数の多さから,計算機の性能向上と共に効率的な近似法や数値解法の開発が進められ発展してきた。
線型偏微分方程式は、解を基底関数 系の線型結合として表す手法がよく用いられる(たとえば正弦波 関数を使ったフーリエ級数 展開など[ 27] )。
線形PDEは、変数分離という重要な手法によって常微分方程式系に帰着させることができる。この手法は微分方程式の解の特徴に基づいている。すなわち、方程式を解き、かつ境界条件を満たす解を何か1つ見つけることができれば、それが「唯一の」解である(これはODEにも当てはまる)。解が各独立変数の関数の積(例えばu ( x , t ) = X ( x ) T ( t ) {\displaystyle u(x,t)=X(x)T(t)} )の形で書けると仮定し、元の偏微分方程式を各変数の常微分方程式へと分解する。[ 28] 。
変数分離法では、多変数の偏微分方程式を、より変数の少ない微分方程式(典型的には常微分方程式)の系に帰着させる。これにより、方程式をより容易に解くことが可能となる。
この手法は、変数分離 可能な形式の方程式(多くは線形)であり、かつ定義域が直方体(区間の直積)や円、球といった対称性を持つ領域である場合に適用できる。変数分離によって現れる常微分方程式の境界値問題は、多くの場合、微分作用素の固有値問題(スツルム=リウヴィル型微分方程式 など)となり、その固有関数系を用いて解を展開・構成することができる。これは線形代数において、行列を対角化して扱う手法と類似した構造を持つ。
変数分離法は、方程式の型によって適用される文脈が異なる。ラプラス方程式 のような楕円型偏微分方程式 では、円や矩形などの対称性を持つ領域上の境界値問題を解くために用いられ、動径方向と角度方向などで変数を分離する。一方、熱方程式 (放物型)や波動方程式 (双曲型)では、時間変数と空間変数を分離することで、空間部分を固有値問題(スツルム=リウヴィル型微分方程式 など)に帰着させ、時間発展を解析する手法として定石となっている。
n =2 次元空間における特性曲面は、特性曲線 と呼ばれる[ 29] 。特別な場合として、一階PDEがODEに帰着する特性曲線を見つけることができる。定義域の座標を変換してこれらの曲線を直線化することで変数分離が可能になり、これは特性曲線法 と呼ばれる。
より一般的に、高次元の一階PDEにこの手法を適用すると、特性曲面を見つけることができる。
積分変換 は、PDEをより単純なもの、特に変数分離可能なPDEに変換することができる。これは作用素の対角化に対応する。
この重要な例がフーリエ解析 であり、正弦波の固有基底 を使用して、定数係数の線形偏微分方程式(熱方程式や波動方程式など)を対角化する。
対象とする領域によって手法は使い分けられる。定義域が有限または周期的であれば、解を三角関数の無限和で表すフーリエ級数 展開が用いられ、これは変数分離法とも密接に関連する。一方、無限領域上の問題では、フーリエ変換 (フーリエ積分)を用いることで、微分演算を代数演算(多項式の掛け算)に変換して解く手法が一般的である。上記の熱方程式の点源に対する解は、フーリエ積分の使用例である。
多くの場合、適切な変数変換 (英語版 ) によって、PDEを既知の解を持つより単純な形式に帰着させることができる。例えば、ブラック–ショールズ方程式 ∂ V ∂ t + 1 2 σ 2 S 2 ∂ 2 V ∂ S 2 + r S ∂ V ∂ S − r V = 0 {\displaystyle {\frac {\partial V}{\partial t}}+{\tfrac {1}{2}}\sigma ^{2}S^{2}{\frac {\partial ^{2}V}{\partial S^{2}}}+rS{\frac {\partial V}{\partial S}}-rV=0} は、次の変数変換によってV ( S , t ) = v ( x , τ ) , x = ln ( S ) , τ = 1 2 σ 2 ( T − t ) , v ( x , τ ) = e − α x − β τ u ( x , τ ) . {\displaystyle {\begin{aligned}V(S,t)&=v(x,\tau ),\\[5px]x&=\ln \left(S\right),\\[5px]\tau &={\tfrac {1}{2}}\sigma ^{2}(T-t),\\[5px]v(x,\tau )&=e^{-\alpha x-\beta \tau }u(x,\tau ).\end{aligned}}} 熱方程式 ∂ u ∂ τ = ∂ 2 u ∂ x 2 {\displaystyle {\frac {\partial u}{\partial \tau }}={\frac {\partial ^{2}u}{\partial x^{2}}}} に帰着可能である[ 30] 。
非斉次方程式は、多くの場合(定数係数PDEの場合は常に)、基本解 (点源に対する解P ( D ) u = δ {\displaystyle P(D)u=\delta } )を見つけ、非斉次項(線源)との畳み込み をとることによって解くことができる。
特定の境界条件を満たす境界値問題を解く場合、基本解の概念はグリーン関数 へと拡張される。グリーン関数は、基本解と同様にデルタ関数に対する応答であるが、領域の境界条件も同時に満たすように構成されたものである。これは信号処理 において、フィルタをそのインパルス応答 によって理解することに類似している。
重ね合わせの原理は、PDEの線形系を含むあらゆる線形系に適用される。この概念の一般的な可視化は、同位相の2つの波が結合してより大きな振幅になる相互作用である(例:sinx + sinx = 2 sinx )。同じ原理は、解が実数または複素数で加法性を持つPDEでも観察される。u 1 とu 2 がある関数空間R における線形PDEの解である場合、任意の定数c 1 およびc 2 に対してu =c 1 u 1 +c 2 u 2 もまた、同じ関数空間におけるそのPDEの解である。
非線形PDEを解くための一般的に適用可能な解析的手法は存在しない。それでも、存在と一意性の結果(コーシー=コワレフスカヤの定理 など)はしばしば可能であり、解の重要な定性的および定量的性質の証明も可能である(これらの結果を得ることは解析学 の主要な部分である)。
それにもかかわらず、いくつかのタイプの方程式には使用できる手法がある。h -原理(英語版 ) は、劣決定系 (英語版 ) 方程式を解くための最も強力な手法である。リキエ=ジャネ理論 は、多くの解析的過剰決定系 (英語版 ) に関する情報を得るための有効な手法である。
特性曲線法 は、非線形偏微分方程式を解くために、いくつかの非常に特殊な場合にのみ使用できる[ 31] 。
場合によっては、解を既知の解を持つ方程式への補正とみなす摂動解析 によってPDEを解くことができる。代替手段として、単純な有限差分法 スキームから、より成熟したマルチグリッド法 や有限要素法 に至る数値解析 技術がある。科学や工学における多くの興味深い問題は、コンピュータ (時には高性能スーパーコンピュータ )を使用してこの方法で解かれている。
PDEを解くための一般的なアプローチは、微分方程式の対称性、つまり解から解への連続的な微小変換 (英語版 ) (リー理論 (英語版 ) )を利用する。19世紀後半にソフス・リー によって開発されたこの手法は、古典的な積分理論を統一的な枠組みで理解することを可能にした。連続群論 、リー代数 、微分幾何学 は、線形および非線形偏微分方程式の構造を理解し、ラックス対 、再帰作用素、ベックルント変換 (英語版 ) を見つけ、最終的にPDEの厳密な解析解を見つけるために使用される。
対称性手法は、数学、物理学、工学、および他の多くの分野で発生する微分方程式を研究するための強力な道具として認識されている。
アドミアン分解法 (英語版 ) [ 32] 、リャプノフ の人工小パラメータ法、および彼のホモトピー摂動法は、すべてより一般的なホモトピー解析法 (英語版 ) の特殊なケースである[ 33] 。これらは級数展開に基づく手法である。リャプノフ法を除き、これらは標準的な摂動論 とは異なり微小な物理パラメータの存在を前提としない。そのため、より柔軟で一般性の高い解を導出できるという利点がある。
偏微分方程式の数値解法 としては、有限要素法 (FEM)、有限体積法 (FVM)、有限差分法 (FDM) の3つが代表的である。各手法は問題の性質に応じて使い分けられる。例えば、有限要素法は複雑な境界形状を持つ構造解析などに強みを持ち、有限体積法は保存則が重要な流体解析(CFD)の分野で広く普及している。近年では、これらの手法の弱点を補うためにメッシュフリー法 などの研究も進められている。
有限要素法 (FEM)(その実用的な応用は有限要素解析 (FEA) としても知られる)は、偏微分方程式 (PDE) や積分方程式の解を、有限個の関数のセットを使用して近似するための数値技術である[ 34] [ 35] 。基本的なアプローチは、微分作用素を離散化することで偏微分方程式を代数方程式系に帰着させる(定常問題)か、時間微分のみを残して常微分方程式系に変換し、それをオイラー法やルンゲ=クッタ法などの標準的な手法で数値積分することに基づいている。
有限差分法は、有限差分 方程式を使用して微分を近似することにより、微分方程式の解を近似する数値手法である。
有限差分法や有限要素法と同様に、メッシュ化されたジオメトリ上の離散的な場所で値が計算される。「有限体積」とは、メッシュ上の各節点を囲む小さな体積を指す。有限体積法では、発散項を含む偏微分方程式の面積分は、発散定理 を使用して体積積分に変換される。これらの項は、各有限体積の表面におけるフラックスとして評価される。与えられた体積に入ってくるフラックスは、隣接する体積から出ていくフラックスと同一であるため、これらの手法は設計上、質量を保存する。
ウィキメディア・コモンズには、
偏微分方程式 に関連するカテゴリがあります。
^ “Regularity and singularities in elliptic PDE's: beyond monotonicity formulas | EllipticPDE Project | Fact Sheet | H2020 ” (英語). CORDIS | European Commission . 2024年2月5日閲覧。 ^ Klainerman, Sergiu (2010). “PDE as a Unified Subject”. In Alon, N.; Bourgain, J.; Connes, A. et al.. Visions in Mathematics . Modern Birkhäuser Classics. Basel: Birkhäuser. pp. 279–315. doi :10.1007/978-3-0346-0422-2_10 . ISBN 978-3-0346-0421-5 ^ Erdoğan, M. Burak; Tzirakis, Nikolaos (2016). Dispersive Partial Differential Equations: Wellposedness and Applications . London Mathematical Society Student Texts. Cambridge: Cambridge University Press. ISBN 978-1-107-14904-5 . https://www.cambridge.org/core/books/dispersive-partial-differential-equations/2DC65286BA080B54EB659E42A553CA88 ^ Evans 1998 , pp. 1–2.^ Klainerman, Sergiu (2008), Partial Differential Equations , Princeton University Press, pp. 455–483 ^ Roger Temam (1984): "Navier–Stokes Equations: Theory and Numerical Analysis", ACM Chelsea Publishing. ^a b Korteweg-de Vries and Nonlinear Schrödinger Equations: Qualitative Theory (2001), Zhidkov, Peter E., Springer. ^ The Nonlinear Schrödinger Equation (1999) -Self-Focusing and Wave Collapse-, Sulem, Catherine, Sulem, Pierre-Louis, Springer. ^ The Nonlinear Schrödinger Equation -Singular Solutions and Optical Collapse- (2015), Gadi Fibich, Springer. ^a b c Levandosky, Julie. “Classification of Second-Order Equations ”. 2026年2月5日閲覧。 ^ Courant and Hilbert (1962), p.182. ^ Evans 1998 , chpt. 6. Second-Order Elliptic Equations.^ R. Courant, D. Hilbert,Methoden Der Mathematischen Physik ,R. クーラン ,D. ヒルベルト (著), 丸山 滋弥, 斎藤 利弥 (翻訳)『数理物理学 の方法』東京図書 ^ 恒藤敏彦『弾性体と流体』岩波書店 〈物理入門コース 8〉、1983年。ISBN 4000076485 。 ^ 際本泰士『振動・波動論講義―物理実験を取り入れて』コロナ社、2005年。ISBN 4339066095 。 ^ シュレーディンガー方程式I, II (朝倉数学大系) by谷島賢二 (著),砂田利一 (編集),増田久弥 (編集), 堀田良之 (編集),朝倉書店 . ^ 大石進一 『フーリエ解析』岩波書店 〈理工系の数学入門コース 6〉、1989年。ISBN 4000077767 。^ 広田良吾 『直接法によるソリトンの数理』(岩波書店 、1992年)^ Smith, G. D. (1985). Numerical solution of partial differential equations: finite difference methods.en:Oxford university press . ^ Strikwerda, J. C. (2004). Finite difference schemes and partial differential equations.SIAM . ^ 森正武 . (1986)有限要素法 とその応用.岩波書店 .^ 菊池文雄. (1999). 有限要素法概説 [新訂版].サイエンス社 . ^ 菊池文雄. (1994). 有限要素法の数理.培風館 . ^ 田端正久; 偏微分方程式の数値解析 , 2010.岩波書店 . ^ 偏微分方程式の数値解法, 編集委員:伊理正夫 ・杉原厚吉 ・速水謙 ・今井浩, 神谷紀生 & 北栄輔著, 工系数学講座 第11巻,978-4-320-01610-1, 1998年03月,共立出版 . ^ 登坂宣好, & 大西和榮. (2003). 偏微分方程式の数値シミュレーション.東京大学出版会 . ^ 大石進一 『フーリエ解析』岩波書店 〈理工系の数学入門コース 6〉、1989年。ISBN 4000077767 。^ Gershenfeld, Neil (2000). The nature of mathematical modeling (Reprinted (with corr.) ed.). Cambridge: Cambridge University Press. p. 27 . ISBN 0521570956 . https://archive.org/details/naturemathematic00gers_334 ^ Zachmanoglou & Thoe 1986 , pp. 115–116.^ Wilmott, Paul; Howison, Sam; Dewynne, Jeff (1995). The Mathematics of Financial Derivatives . Cambridge University Press. pp. 76–81. ISBN 0-521-49789-2 . https://books.google.com/books?id=VYVhnC3fIVEC&pg=PA76 ^ Logan, J. David (1994). “First Order Equations and Characteristics”. An Introduction to Nonlinear Partial Differential Equations . New York: John Wiley & Sons. pp. 51–79. ISBN 0-471-59916-6 ^ Adomian, G. (1994). Solving Frontier problems of Physics: The decomposition method . Kluwer Academic Publishers. ISBN 9789401582896 . https://books.google.com/books?id=UKPqCAAAQBAJ&q=%22partial+differential%22 ^ Liao, S. J. (2003). Beyond Perturbation: Introduction to the Homotopy Analysis Method . Boca Raton: Chapman & Hall/ CRC Press. ISBN 1-58488-407-X ^ Solin, P. (2005). Partial Differential Equations and the Finite Element Method . Hoboken, New Jersey: J. Wiley & Sons. ISBN 0-471-72070-4 ^ Solin, P.; Segeth, K. & Dolezel, I. (2003). Higher-Order Finite Element Methods . Boca Raton: Chapman & Hall/CRC Press. ISBN 1-58488-438-X Lawrence C. Evans:Partial Differential Equations (Graduate Studies in Mathematics)American Mathematical Society , 2010/04/. Dorina Mitrea:Distributions, Partial Differential Equations, and Harmonic Analysis (Universitext) , Springer, 2019/01. Michael Taylor:Partial Differential Equations I-III (Applied Mathematical Sciences), Springer. Egorov, Y. V., & Shubin, M. A. (2013):Foundations of the classical theory of partial differential equations .en:Springer Science & Business Media . Egorov, Komech and Shubin :Elements of the Modern Theory of Partial Differential Equations (1999) Springer. Sommerfeld, A. (1949):Partial differential equations in physics .en:Academic Press . Renardy, M., & Rogers, R. C. (2006):An introduction to partial differential equations .en:Springer Science & Business Media . Logan, J. D. (2008):An introduction to nonlinear partial differential equations . John Wiley & Sons. Olver, P. J.:Introduction to partial differential equations . Berlin: Springer. Fritz Schwarz:Loewy Decomposition of Linear Differential Equations , Springer,ISBN 978-3-7091-1286-1 (2012).