1 Gauss 核と近似単位
命題 1.1.n≥1を整数とし、t>0に対してRn上の Gauss 核を
Gt(x)=(4πt)−n/2exp(−4t∣x∣2)と定める。このときGt>0、∫RnGt(x)dx=1であり、任意のs,t>0について
Gt∗Gs=Gt+sが成り立つ。1≤p<∞、f∈Lp(Rn;C)ならば
t↓0lim∥Gt∗f−f∥p=0.また、f∈L1(Rn;C)ならば
∫Rn(Gt∗f)(x)dx=∫Rnf(x)dx.
証明. 正値性は式から従う。位相因子e−2πix⋅ξによる Fourier 変換について、§E15.26 定理 4.1により
Gt(ξ)=e−4π2t∣ξ∣2,Gt∗h=etΔh(h∈L2(Rn))が成り立つ。ξ=0とすれば∫Gt=1である。Gs∈L2に対するetΔの Fourier 乗数はe−4π2t∣ξ∣2なので、Gt∗Gsの Fourier 変換はe−4π2(t+s)∣ξ∣2=Gt+sである。したがってGt∗Gs=Gt+sをL2で得る。両辺は連続関数であるから、各点でも等しい。
1≤p<∞とする。§E14.7 定理 2.2の Bochner 積分表示とノルム評価により
∥Gt∗f−f∥p≤∫RnGt(y)∥f(⋅−y)−f∥pdy.§E14.7 定理 1.2により、任意のε>0に対して、∣y∣<δならば∥f(⋅−y)−f∥p<εとなるδ>0が存在する。したがって
∥Gt∗f−f∥p≤ε+2∥f∥p∫∣y∣≥δGt(y)dy.y=tzと変数変換すると、最後の積分は∫∣z∣≥δ/tG1(z)dzであり、t↓0で零へ収束する。εは任意なので、Lp収束を得る。
f∈L1のとき、∬Gt(y)∣f(x−y)∣dxdy=∥f∥1である。Fubini の定理により積分順序を交換すると、質量保存の等式を得る。▨
定理 1.2.n≥1、1≤p≤q≤∞、f∈Lp(Rn;C)とし、αを多重指数とする。このとき、u(t,x)=(Gt∗f)(x)は(0,∞)×Rn上で滑らかであり、ut=Δuを満たす。また、tとfに依存しない定数Cn,p,q,αに対して
∥Dxαu(t)∥q≤Cn,p,q,αt−∣α∣/2−2n(1/p−1/q)∥f∥p(t>0)が成り立つ。ただし1/∞=0とする。
証明. 核の拡大縮小から
DαGt(x)=t−(n+∣α∣)/2(DαG1)(x/t)である。DαG1は多項式と Gauss 関数との積なので、任意の1≤r≤∞に対して
∥DαGt∥r=t−∣α∣/2−2n(1−1/r)∥DαG1∥r.p<∞のとき、§E14.7 定理 2.2により
∥(DαGt)∗f∥p≤∥DαGt∥1∥f∥p.1/p+1/p′=1とおく。§E10.22 定理 4.2を各xで適用すると、端点を含む任意のpについて
∣((DαGt)∗f)(x)∣≤∥DαGt∥p′∥f∥p.p<q<∞のとき、h=(DαGt)∗fに対して
∥h∥qq=∫∣h∣q−p∣h∣p≤∥h∥∞q−p∥h∥ppである。二つの端の評価を代入すると、主張のtの指数を得る。p=q=∞のときは
∥(DαGt)∗f∥∞≤∥DαGt∥1∥f∥∞を積分から直接得る。
核の任意の時間・空間微分は、正時間では多項式と Gauss 関数との積である。tが正のコンパクト区間を動き、xがコンパクト集合を動くとき、y↦∂tkDxαGt(x−y)とその差商は、Lp′に属する共通の減衰関数で支配される。p′=∞の場合にも同じ減衰関数による一様評価が成り立つ。Hölder の不等式と優収束により、畳込みの下で任意階の微分を行うことができ、各微分は連続である。直接微分による∂tGt=ΔGtからut=Δuを得る。▨
命題 1.3.0≤f∈L1(Rn;R)が零関数でないならば、任意のt>0、x∈Rnについて(Gt∗f)(x)>0である。特に、零でない非負のコンパクトな台をもつ初期値から出発しても、正時間の熱解の台はRn全体である。
証明.f>0となる集合は正の測度をもち、Gt(x−y)>0はすべてのyで成り立つ。したがって非負関数Gt(x−y)f(y)の積分は正である。積分の有限性と連続性は定理 1.2による。▨
命題 1.5.n≥1とし、Ktを§E15.32 定理 3.1のトーラスTn=(R/Z)nの熱核、Gtを命題 1.1の全空間の Gauss 核とする。Q=[−21,21]nの点をTnの点の代表と見なす。このとき、次の二つが成り立つ。
- nだけに依存する定数Cnが存在し、0<t≤1とx∈Qに対して
0<Kt(x)−Gt(x)≤Cn(4πt)−n/2e−1/(32t)
が成り立つ。特に、supx∈Q∣Kt(x)−Gt(x)∣はt↓0で零へ収束する。
- nだけに依存する定数Cn′が存在し、t≥1とx∈Qに対して
∣Kt(x)−1∣≤Cn′e−4π2t,0<Gt(x)≤(4πt)−n/2
が成り立つ。特に、t→∞でKtは定数1へ一様収束し、Gtは零へ一様収束する。また、u0∈L1(Tn)に対して(Kt∗u0)(x)=∫TnKt(x−y)u0(y)dyとおくと、∣(Kt∗u0)(x)−∫Tnu0dy∣≤Cn′e−4π2t∥u0∥1が成り立つ。
証明.主張 (1)を示す。§E15.32 定理 3.1によりx∈Qに対してKt(x)=∑m∈ZnGt(x+m)であるから、Kt(x)−Gt(x)=∑m=0Gt(x+m)>0である。m=0の最大ノルムをj=∣m∣∞≥1とすると、∣x+m∣≥∣x+m∣∞≥j−21≥j/2であるからGt(x+m)≤(4πt)−n/2e−j2/(16t)である。0<t≤1ではj2/(16t)≥1/(32t)+j2/32が成り立つ。∣m∣∞=jを満たすmの個数は(2j+1)n−(2j−1)n≤2n(2j+1)n−1であるから
Kt(x)−Gt(x)≤(4πt)−n/2e−1/(32t)j≥1∑2n(2j+1)n−1e−j2/32である。右辺の級数は収束するので、その値をCnとすればよい。t−n/2e−1/(32t)→0により後半を得る。
主張 (2)を示す。§E15.32 定理 3.1によりKt(x)=1+∑k=0e−4π2t∣k∣2e2πik⋅xである。k∈Zn∖{0}とt≥1に対してt∣k∣2≥(t−1)+∣k∣2であるから
∣Kt(x)−1∣≤e−4π2(t−1)k=0∑e−4π2∣k∣2となり、Cn′=e4π2∑k=0e−4π2∣k∣2とおけばよい。Gt(x)≤Gt(0)=(4πt)−n/2は式から従う。最後の評価は、∫TnKt(x−y)dy=1により
(Kt∗u0)(x)−∫Tnu0dy=∫Tn(Kt(x−y)−1)u0(y)dyと書き、被積分関数の核の部分を上の評価で抑えれば従う。▨
例 1.6.s>0を固定し、初期値をf=Gsとする。このとき、命題 1.1により、熱解はu(t)=Gt+sである。a=t+sとおくと、Gauss 積分から
∥u(t)∥22=(8πa)−n/2.∇Ga=−xGa/(2a)であり、各座標について部分積分すると
∫Rn∣x∣2e−∣x∣2/(2a)dx=na∫Rne−∣x∣2/(2a)dx.したがって
∥∇u(t)∥22=4an(8πa)−n/2,dtd∥u(t)∥22=−2an(8πa)−n/2=−2∥∇u(t)∥22.これはA=−Δに対する熱のエネルギー等式§E15.26 定理 1.1と一致する。
例 1.7. 一次元でf=1(0,∞)とする。このとき、Φ(a)=∫−∞aG1(z)dzとおくと
(Gt∗f)(x)=Φ(x/t).G1は正で偶関数なので、Φは狭義単調増加であり、Φ(0)=1/2である。x<0では誤差はΦ(x/t)<1/2、x>0では誤差は1−Φ(x/t)<1/2である。両側からx→0とすると誤差は1/2に収束するので
∥Gt∗f−f∥∞=21(t>0).したがって、近似単位のLpノルム初期収束から有限指数という仮定を除くことはできない。
2 一次元の波動解
証明.r=x+t、s=x−tと変数変換すると、utt−uxx=−4ursである。urs=0を積分すると、各古典解はu(t,x)=a(x+t)+b(x−t)の形をもつ。初期条件からa+b=f、a′−b′=gである。よってa′=(f′+g)/2、b′=(f′−g)/2を積分して主張の式を得る。逆に、その式を微分すると方程式と両初期条件が成り立つ。
t>0のときS1(t)の畳込み核は211[−t,t]である。§E15.6 命題 4.2により、この核は因果的波動基本解の時刻tにおける空間核と一致する。積分端を微分すると∂tS1(t)f=(f(⋅+t)+f(⋅−t))/2を得る。▨
命題 2.2.f∈H1(R;C)、g∈L2(R;C)、T>0とする。このとき、定理 2.1の式によって定めたuは、
C([0,T];H1(R))∩C1([0,T];L2(R))に属する唯一のエネルギー解である。任意のT>0とF∈C1([0,T]×R;C)に対して、零初期値の古典解の外力項は
w(t,x)=21∫0t∫x−(t−s)x+(t−s)F(s,y)dydsであり、wtt−wxx=Fを満たす。さらにF∈L1(0,T;L2(R))の場合も、この式の積分をL2値の積分として定めると、得られる関数wは零初期値のエネルギー解である。
証明.τah(x)=h(x+a)とおく。L2における平行移動の等長性と連続性§E14.7 定理 1.2、および積分表示から
∥S1(t)h∥2≤t∥h∥2,∂xS1(t)h=2τth−τ−th,∂tS1(t)h=2τth+τ−thが成り立つ。後二式の空間微分は弱微分、時間微分はL2ノルム微分である。空間微分の式は試験関数に対する積分端の計算から従い、時間微分の式は平行移動の連続性と積分の基本定理から従う。
H1とL2でそれぞれf,gを滑らかなコンパクトな台をもつ関数で近似する。S1(t)とその時間・空間微分のL2評価を近似解の差に適用すると、近似解はC([0,T];H1)∩C1([0,T];L2)で主張の式による関数へ収束する。近似解の方程式を試験関数と組み合わせて極限を取ると、極限も弱方程式と初期条件を満たす。§E15.27 定理 2.2をA=−∂x2に適用すれば一意性を得る。
F∈C1の場合、積分端と時間を微分すると
wt(t,x)wtt(t,x)wxx(t,x)=21∫0t(F(s,x+t−s)+F(s,x−t+s))ds,=F(t,x)+21∫0t(Fx(s,x+t−s)−Fx(s,x−t+s))ds,=21∫0t(Fx(s,x+t−s)−Fx(s,x−t+s))ds.したがって方程式と零初期条件が成り立つ。F∈Lt1Lx2のときはw(t)=∫0tS1(t−s)F(s)dsがH1値の Bochner 積分となる。t>0のときS1(t)は核211[−t,t]との畳込みであり、この核は§E15.6 命題 4.2により因果的波動基本解の時刻tにおける空間核である。その Fourier 変換は
21∫−tte−2πixξdx=2πξsin(2πtξ)であり、これはA=−∂x2の Fourier 乗数λ=4π2ξ2に対するSA(t)の乗数sin(tλ)/λに等しい。有界な乗数が等しいからS1(t)=SA(t)である。したがってw(t)=∫0tSA(t−s)F(s)dsは§E15.27 定理 2.2の Duhamel 公式の外力項であり、この積分はエネルギー解である。▨
例 2.3.f∈Cc∞(R;R)、g=−f′とする。このとき、積分の基本定理と定理 2.1により
u(t,x)=2f(x+t)+f(x−t)−2f(x+t)−f(x−t)=f(x−t).したがって
21∫R(∣ut(t,x)∣2+∣ux(t,x)∣2)dx=∫R∣f′(x−t)∣2dx=∥f′∥22.
3 球面平均と Kirchhoff の公式
補題 3.1.n≥1を整数、h∈C2(Rn;C)とする。σn−1=∣Sn−1∣とし、球面平均を
Mrh(x)=σn−11∫Sn−1h(x+rω)dσ(ω)(r≥0)と定める。n=1ではS0={−1,1}上の数え上げ測度を用いる。このときr>0において
(∂r2+rn−1∂r)Mrh=MrΔh=ΔxMrh.またM0h=h、∂rMrh∣r=0=0である。
証明.n≥2のとき、§E15.6 補題 2.1をU=Rn、f=h、a=xとして適用すると、0<rに対して
σn−1rn−1∂rMrh(x)=∫B(x,r)Δh(y)dyが成り立つ。右辺を極座標で表すとσn−1∫0rsn−1MsΔh(x)dsであり、被積分関数はsについて連続であるから、rで微分して
∂r(rn−1∂rMrh(x))=rn−1MrΔh(x)を得る。rn−1で割れば最初の等式を得る。n=1ではMrh(x)=(h(x+r)+h(x−r))/2の二回微分から同じ式を得る。xに関する微分と球面積分の交換からΔxMrh=MrΔhである。r=0での微分は∇h(x)⋅σn−1−1∫Sn−1ωdσ=0である。▨
補題 3.2.n≥1とし、v∈C2([0,T]×Rn;C)がvtt−Δv=0を満たすとする。このとき、あるx0∈RnとR>0に対し、v(0)=vt(0)=0がB(x0,R)上で成り立つならば
v(t,x)=0(0≤t≤T, ∣x−x0∣<R−t)である。特に、初期値が同じ二つの古典解は一致する。
証明.e=(∣vt∣2+∣∇v∣2)/2、J=ℜ(vt∇v)とおくと、方程式から∂te=divJ、∣J∣≤eである。0≤t<min{T,R}に対し、球面での積分と Gauss の公式から
dtd∫B(x0,R−t)e(t,x)dx=∫∂B(x0,R−t)(J⋅ν−e)dσ≤0.初期エネルギーは零なので、vt=0は縮小する球の中で成り立つ。∣x−x0∣<R−tならば、すべての0≤s≤tでx∈B(x0,R−s)である。よってv(t,x)=v(0,x)+∫0tvs(s,x)ds=0となる。端の時刻は連続性で扱う。同じ初期値をもつ解の差に適用し、各点についてR>tを選べば一意性を得る。▨
証明.h∈C2(R3)に対しv(t,x)=tMth(x)とおく。補題 3.1により、t>0では
vtt=2∂tMth+t∂t2Mth=tΔxMth=Δxv.同じ補題によりv(0)=0、vt(0)=h、vtt(0)=0である。h=f∈C3の場合、vの方程式を時間微分して∂tvも波動方程式を満たす。h=g∈C2の場合のvと加えると、主張のuはC2級であり、方程式と初期条件を満たす。一意性は補題 3.2による。
t>0のとき、tMtの畳込み核は球面∣y∣=t上の測度(4πt)−1dσ(y)である。§E15.6 命題 4.2により、この球面測度は因果的波動基本解の時刻tにおける空間核と一致する。▨
4 Hadamard の降下法と伝播
証明.hを第三座標に依存しない関数H(x1,x2,x3)=h(x1,x2)に延長する。半径tの球面の上下半球をz=(z1,z2)の円板へ射影すると、それぞれのグラフz3=±t2−∣z∣2の面積要素は
dσ=t2−∣z∣2tdzである。したがって、二つの半球を加えることにより
S3(t)H(x,x3)=2π1∫∣z∣<tt2−∣z∣2h(x+z)dz=S2(t)h(x).特異因子は可積分である。実際、z=twと変数変換すれば
S2(t)h(x)=2πt∫∣w∣<11−∣w∣2h(x+tw)dw,∫∣w∣<11−∣w∣2dw=2π.固定した円板上の表示では、データの許す次数まで積分の下で微分することができる。
延長した初期値に定理 3.3を適用する。得られる解は第三座標に依存しないので、三次元の Laplace 作用素は二次元の Laplace 作用素に等しい。よって主張の式が二次元の方程式と初期条件を満たす。一意性は補題 3.2による。第三座標に一定な延長を全空間のL2データと見なす必要はない。有限な依存領域で延長を切断した場合にも、§E15.27 定理 5.1により同じ局所解を得る。▨
命題 4.2. 三次元の斉次波動方程式の初期値をf∈C3(R3)、g∈C2(R3)とする。このとき、t>0、x∈R3に対し、fとgが球面{y:∣y−x∣=t}のある開近傍で零ならば、解はu(t,x)=0を満たす。
証明. 球面はコンパクトであるから、あるε>0が存在して、∣r−t∣<εならば球面∣y−x∣=r上でf=g=0となる。したがってMrf(x)=Mrg(x)=0がそのrの範囲で成り立つ。定理 3.3の二項は、ともに零である。▨
例 4.3. 各次元n∈{1,2,3}でf=0とし、0≤g∈Cc∞(Rn;R)、g=0を取る。suppg⊂B(0,R)とし、xを固定する。このとき、t>∣x∣+Rならば、三次元では球面∣y−x∣=tがsuppgから離れるので
u(t,x)=4πt1∫∣y−x∣=tg(y)dσ(y)=0.二次元では、同じ時刻にsuppg⊂B(x,t)であり、
u(t,x)=2π1∫suppgt2−∣x−y∣2g(y)dy>0.しかもt→∞において、台上で一様にt/t2−∣x−y∣2→1となるので
tu(t,x)⟶2π1∫R2g(y)dy>0.一次元では、区間[x−t,x+t]が台全体を含むため
u(t,x)=21∫Rg(y)dy>0.いずれの次元でも、解に影響する初期値は閉球B(x,t)の中にある。三次元ではさらに球面の近傍だけで十分であるという強い Huygens の原理が成り立つが、上の一次元と二次元の初期値は球面近傍で零であっても解を零にしない。