ラプラス方程式の基本解
📂偏微分方程式ラプラス方程式の基本解
ビルドアップ
ラプラス方程式は回転変換に対して不変なので、u(x)の変数を半径に変更して考えます。そうすれば、次のような過程を経て、微分方程式をより簡単な形にすることができます。
u=u(x)をラプラス方程式の解としましょう。
Δu=0
そして、r=∣x∣=(x12+⋯+xn2)1/2と置き、v∈C2で、u(x)=v(∣x∣)=v(r)(x∈Rn∖{0})とします。
v(r)Δv=u(x)=0
今、uのラプラス方程式を表すために、次の微分を計算しましょう。
∂xi∂r=∂xi∂(x12+⋯xn2)1/2=2(x12+⋯xn2)1/212xi=(x12+⋯xn2)1/2xi=rxianduxi(x)=∂xi∂v(r)=drdv(r)∂xi∂r=v′(r)rxi
uxixi=∂xi∂(v′(r)rxi)=∂xi∂v′(r)rxi+v′(r)∂xi∂(rxi)=drdv′(r)∂xi∂rrxi+v′(r)(r1+xi∂xi∂(rxi))=v′′(r)r2xi2+v′(r)(r1+xidrd(r1)∂xi∂r)=v′′(r)r2xi2+v′(r)(r1+xi(−r21)rxi)=v′′(r)r2xi2+v′(r)(r1−r3xi2)
そうすると、ラプラス方程式は次のようになります。
Δu=i∑nuxixi=r2v′′(r)(x12+⋯xn2)+v′(r)(rn−r3x12+⋯xn2)=r2v′′(r)r2+v′(r)(rn−r3r2)=v′′(r)+rn−1v′(r)x=0
したがって、以下の二つの式は等しいです。
Δu=0 in Rn∖{0}⟺v′′(r)+rn−1v′(r)=0(r>0)
u(x)のラプラス方程式を解く問題が、v(r)の2次常微分方程式を解く問題に変わったわけです。
今、v′′+rn−1v′=0 in (0,∞)と、v′=0 in (0,∞)と仮定しましょう。
そうすると、整理して次の式を得ます。
v′v′′=r1−n
w=v′∈C1(0,∞)と置換すると、次を得ます。
ww′=r1−n
左辺を積分すると、以下のようになります。
∫1sww′dr=log∣w(s)∣−log∣w(1)∣=logw(1)w(s)
右辺を積分すると、以下のようになります。
∫1sr1−ndr=(1−n)[logs−log1]=(1−n)logs=logs1−n
したがって、次を得ます。
∣w(1)∣∣w(s)∣=s1−n⟹∣w(s)∣=∣w(1)∣s1−n(s>0)
再びvで表すと、次のようになります。
v′(s)=w(s)=w(1)s1−n=v′(1)s1−n
ここで、v(r)−v(1)に微分積分学の基本定理を適用し、上の式を代入して整理すると、次を得ます。
v(r)−v(1)=∫1rv′(s)ds=v′(1)∫1rs1−nds=⎩⎨⎧v′(1)logrv′(1)2−n1(r2−n−1)n=2n≥3
v(r)に対して整理すると、次のようになります。
v(r)=⎩⎨⎧v′(1)logr+v′(1)v′(1)2−n1rn−21+(v(1)+n−2v′(1))n=2n≥3
定数部分をb,cと表示すると、次のように簡単に整理できます。
v(r)=⎩⎨⎧blogr+crn−2b+cn=2n≥3
これらの結果から、ラプラス方程式の基本解を定義します。
定義
x∈Rnで、x=0に対して、以下の関数Φをラプラス方程式の基本解fundamental solutionと定義します。
Φ(x):=⎩⎨⎧−2π1log∣x∣n(n−2)α(n)1∣x∣n−21n=2n≥3
ここで、α(n)はRnのユニットボールB(0, 1)の体積です。nα(n)はRnのユニットボールの表面積です。