§E15.31定常問題の Fourier 級数解

最終更新

矩形や円板の定常境界値問題では、変数分離によって特定の周波数をもつ解が得られる。一般の外力や境界値に対する解をこれらの解の重ね合わせとして表すには、固有関数系の完全性と級数の収束を確かめる必要がある。Fourier 級数は、この重ね合わせを関数や超関数の空間で扱い、微分方程式を係数の方程式として表すための基本的な方法である。例えば、トーラスの Poisson 方程式では、定数関数に対応する係数の方程式から、外力の平均が零であるという可解条件が得られる。このように、定常問題の可解性や解の正則性は、周波数ごとの方程式と係数の減衰を通して調べることができる。本記事では、Fourier 級数による代表的な定常境界値問題の解法を解説する。

1 周期関数と周期超関数

定義 1.1 (周期 Fourier 係数).n∈N≥1n\in\NNとし、Tn=(R/Z)n\mathbb T^n=(\R/\Z)^nに全測度11の Lebesgue 測度を入れる。k∈Znk\in\Z^nに対してek(x)=e2πik⋅xe_k(x)=e^{2\pi i k\cdot x}とおく。f∈L1(Tn)f\in L^1(\mathbb T^n)の Fourier 係数 (Fourier coefficient) を

f^(k)=∫Tnf(x)e−2πik⋅x dx\widehat f(k)=\int_{\mathbb T^n}f(x)e^{-2\pi i k\cdot x}\,dx

と定める。L2L^2の内積は(f,g)2=∫fg‾(f,g)_2=\int f\overline gとし、第一変数について線形とする。

C∞(Tn)C^\infty(\mathbb T^n)には半ノルムpm(φ)=max⁡∣α∣≤m∥∂αφ∥∞p_m(\varphi)=\max_{|\alpha|\leq m}\|\partial^\alpha\varphi\|_\infty、m∈N≥0m\in\Nが定める位相を入れる。その連続な複素線形双対をD′(Tn)\mathcal D'(\mathbb T^n)と書き、その元を 周期超関数 (periodic distribution) という。T∈D′(Tn)T\in\mathcal D'(\mathbb T^n)に対して

T^(k)=⟨T,e−k⟩,⟨∂αT,φ⟩=(−1)∣α∣⟨T,∂αφ⟩\widehat T(k)=\langle T,e_{-k}\rangle,\qquad \langle\partial^\alpha T,\varphi\rangle =(-1)^{|\alpha|}\langle T,\partial^\alpha\varphi\rangle

と定める。双対の作用には共役を入れない。関数f∈L1(Tn)f\in L^1(\mathbb T^n)は⟨f,φ⟩=∫fφ\langle f,\varphi\rangle=\int f\varphiにより周期超関数とみなす。

命題 1.2.n∈N≥1n\in\NNとする。このとき、(ek)k∈Zn(e_k)_{k\in\Z^n}はL2(Tn)L^2(\mathbb T^n)の完全正規直交系であり、すべてのf∈L2(Tn)f\in L^2(\mathbb T^n)について

f=∑k∈Znf^(k)ekL2 で,∥f∥22=∑k∈Zn∣f^(k)∣2f=\sum_{k\in\Z^n}\widehat f(k)e_k\quad\text{$L^2$ で}, \qquad \|f\|_2^2=\sum_{k\in\Z^n}|\widehat f(k)|^2

が成り立つ。

φ∈C∞(Tn)\varphi\in C^\infty(\mathbb T^n)ならば、任意のN∈N≥0N\in\Nに対してsup⁡k(1+∣k∣)N∣φ^(k)∣<∞\sup_k(1+|k|)^N|\widehat\varphi(k)|<\inftyである。逆に、この評価をすべてのNNについて満たす列(ak)(a_k)の級数∑kakek\sum_k a_ke_kは、任意の階数の導関数まで絶対一様収束し、Fourier 係数が(ak)(a_k)である滑らかな関数を定める。滑らかな関数の Fourier 級数はC∞(Tn)C^\infty(\mathbb T^n)で元の関数へ収束する。

証明. 円周の完全系§E14.17 命題 1.2に、積測度のL2L^2との同定§E12.9 定理 8.1とテンソル積の完全系§E12.9 命題 4.1を有限回適用する。各因子の測度は有限であるから、積測度の同定のσ-有限性の仮定を満たす。得られる積関数はeke_kであり、完全性と Parseval の等式が従う。

周期性によって部分積分の境界項が消えるので、任意のM∈N≥0M\in\Nについて

(1+4π2∣k∣2)Mφ^(k)=(1−Δ)Mφ^(k),∣φ^(k)∣≤CM,n(1+∣k∣)−2Mp2M(φ)(1+4\pi^2|k|^2)^M\widehat\varphi(k) =\widehat{(1-\Delta)^M\varphi}(k),\qquad |\widehat\varphi(k)|\leq C_{M,n}(1+|k|)^{-2M}p_{2M}(\varphi)

である。したがって Fourier 係数は任意の次数で減衰する。

2j≤1+∣k∣<2j+12^j\leq1+|k|<2^{j+1}を満たす格子点の個数はCn2jnC_n2^{jn}以下であるから、s>ns>nならば∑k(1+∣k∣)−s<∞\sum_k(1+|k|)^{-s}<\inftyである。急減少する列(ak)(a_k)と任意の多重指数α\alphaに対して、N>∣α∣+nN>|\alpha|+nを選ぶと

∑k∣(2πik)αak∣<∞\sum_k|(2\pi i k)^\alpha a_k|<\infty

を得る。各階の導関数の級数が絶対一様収束するため、その和は滑らかであり、積分と和を交換すると Fourier 係数はaka_kとなる。ak=φ^(k)a_k=\widehat\varphi(k)の場合には、完全性により和とφ\varphiがL2L^2で一致する。両者は連続であるから各点でも一致する。▨

定理 1.3.n∈N≥1n\in\NNとする。このとき、T∈D′(Tn)T\in\mathcal D'(\mathbb T^n)の Fourier 係数は多項式増大する。すなわち、あるC≥0C\geq0とm∈N≥0m\in\Nに対して

∣T^(k)∣≤C(1+∣k∣)m(k∈Zn)|\widehat T(k)|\leq C(1+|k|)^m\qquad(k\in\Z^n)

が成り立つ。逆に、この評価を満たす任意の列(ck)k∈Zn(c_k)_{k\in\Z^n}は

⟨T,φ⟩=∑k∈Znckφ^(−k)(φ∈C∞(Tn))\langle T,\varphi\rangle=\sum_{k\in\Z^n}c_k\widehat\varphi(-k) \qquad(\varphi\in C^\infty(\mathbb T^n))

により Fourier 係数がckc_kである周期超関数を一意に定める。この和は絶対収束し、∑∣k∣≤Rckek\sum_{|k|\leq R}c_ke_kはR→∞R\to\inftyのとき各試験関数に対してTTへ収束する。また、

∂αT^(k)=(2πik)αT^(k)\widehat{\partial^\alpha T}(k)=(2\pi i k)^\alpha\widehat T(k)

が成り立つ。

周期化

Pψ(x)=∑ℓ∈Znψ(x+ℓ)P\psi(x)=\sum_{\ell\in\Z^n}\psi(x+\ell)

を用いた⟨T~,ψ⟩=⟨T,Pψ⟩\langle\widetilde T,\psi\rangle=\langle T,P\psi\rangleは、S(Rn)\mathcal S(\R^n)上の緩増加超関数を定める。T↦T~T\mapsto\widetilde Tは、D′(Tn)\mathcal D'(\mathbb T^n)とZn\Z^nの平行移動で不変なRn\R^n上の超関数との一対一対応である。

証明.TTの零における連続性により、あるm∈N≥0m\in\NとC≥0C\geq0に対して∣⟨T,φ⟩∣≤Cpm(φ)|\langle T,\varphi\rangle|\leq Cp_m(\varphi)が成り立つ。pm(e−k)≤Cm,n(1+∣k∣)mp_m(e_{-k})\leq C_{m,n}(1+|k|)^mを代入すれば係数の多項式増大を得る。

∣ck∣≤C(1+∣k∣)m|c_k|\leq C(1+|k|)^mとし、2M>m+n2M>m+nを選ぶ。命題 1.2の証明中の評価により

∑k∣ckφ^(−k)∣≤C′p2M(φ)∑k(1+∣k∣)m−2M<∞\sum_k|c_k\widehat\varphi(-k)| \leq C' p_{2M}(\varphi)\sum_k(1+|k|)^{m-2M}<\infty

である。したがって定義式は連続な複素線形汎関数を定める。φ=e−j\varphi=e_{-j}を代入するとT^(j)=cj\widehat T(j)=c_jとなる。滑らかな試験関数の Fourier 級数は命題 1.2によりC∞C^\inftyで収束するから、任意の周期超関数SSについて

⟨S,φ⟩=∑jφ^(j)⟨S,ej⟩=∑jS^(−j)φ^(j)\langle S,\varphi\rangle =\sum_j\widehat\varphi(j)\langle S,e_j\rangle =\sum_j\widehat S(-j)\widehat\varphi(j)

である。係数はSSを一意に定める。上の絶対収束する和の末尾を評価すれば部分和の超関数収束も従う。微分の定義へe−ke_{-k}を代入すると微分の係数公式を得る。

ψ∈S(Rn)\psi\in\mathcal S(\R^n)とし、整数N>nN>nを選ぶ。x∈[0,1]nx\in[0,1]^nについて1+∣x+ℓ∣1+|x+\ell|と1+∣ℓ∣1+|\ell|の比はℓ\ellによらず上下から正の定数で抑えられるので、任意のr∈N≥0r\in\Nに対して

pr(Pψ)≤Cr,N,nmax⁡∣α∣≤rsup⁡y∈Rn(1+∣y∣)N∣∂αψ(y)∣.p_r(P\psi)\leq C_{r,N,n} \max_{|\alpha|\leq r}\sup_{y\in\R^n}(1+|y|)^N|\partial^\alpha\psi(y)|.

各階の導関数の級数は基本領域で一様収束するので、PψP\psiは滑らかな周期関数となる。r=mr=mとしてTTの有限階評価を適用し、§E15.4 命題 1.2を用いるとT~∈S′(Rn)\widetilde T\in\mathcal S'(\R^n)を得る。P(ψ(⋅+ℓ))=PψP(\psi(\cdot+\ell))=P\psiであるからT~\widetilde Tは整数平行移動で不変である。

非負関数ρ∈Cc∞(Rn)\rho\in C_c^\infty(\R^n)を[0,1]n[0,1]^n上で正となるように選ぶ。h(x)=∑ℓρ(x+ℓ)h(x)=\sum_\ell\rho(x+\ell)は正の滑らかな周期関数であり、η=ρ/h\eta=\rho/hとおくとη∈Cc∞(Rn)\eta\in C_c^\infty(\R^n)かつPη=1P\eta=1となる。整数平行移動で不変なS∈D′(Rn)S\in\mathcal D'(\R^n)に対して

⟨TS,φ⟩=⟨S,ηφ⟩(φ∈C∞(Tn))\langle T_S,\varphi\rangle=\langle S,\eta\varphi\rangle \qquad(\varphi\in C^\infty(\mathbb T^n))

と定める。右辺ではφ\varphiを周期関数として持ち上げる。固定された台supp⁡η\supp\etaに§E15.1 命題 1.2を適用し、積の微分を評価すると∣⟨TS,φ⟩∣≤Cpr(φ)|\langle T_S,\varphi\rangle|\leq C p_r(\varphi)を得るので、TST_Sは周期超関数である。

ψ∈Cc∞(Rn)\psi\in C_c^\infty(\R^n)に対して、supp⁡η\supp\eta上で周期化の和に寄与する項は有限個である。SSの平行移動不変性により

⟨S,ηPψ⟩=∑ℓ⟨S,η(x)ψ(x+ℓ)⟩=∑ℓ⟨S,η(x−ℓ)ψ(x)⟩=⟨S,ψ⟩.\langle S,\eta P\psi\rangle =\sum_\ell\langle S,\eta(x)\psi(x+\ell)\rangle =\sum_\ell\langle S,\eta(x-\ell)\psi(x)\rangle =\langle S,\psi\rangle.

したがってTS~=S\widetilde{T_S}=SがD′(Rn)\mathcal D'(\R^n)で成り立つ。他方、周期的なφ\varphiに対してP(ηφ)=φP(\eta\varphi)=\varphiであるから、TT~=TT_{\widetilde T}=Tである。したがって二つの構成は互いに逆である。逆写像は一意であるから、TST_Sはη\etaの選択にもよらない。▨

例 1.4.n∈N≥1n\in\NNとする。このとき、Tn\mathbb T^nの原点での Dirac 超関数δ0\delta_0はδ0^(k)=1\widehat{\delta_0}(k)=1を満たす。従って定理 1.3により

δ0=∑k∈Znek,∂αδ0=∑k∈Zn(2πik)αek\delta_0=\sum_{k\in\Z^n}e_k,\qquad \partial^\alpha\delta_0=\sum_{k\in\Z^n}(2\pi i k)^\alpha e_k

が超関数として成り立つ。全空間への延長は

⟨δ0~,ψ⟩=Pψ(0)=∑ℓ∈Znψ(ℓ)\langle\widetilde{\delta_0},\psi\rangle =P\psi(0)=\sum_{\ell\in\Z^n}\psi(\ell)

であるから、δ0~=∑ℓ∈Znδℓ\widetilde{\delta_0}=\sum_{\ell\in\Z^n}\delta_\ellは Dirac comb である。

一方、ck=e∣k∣c_k=e^{|k|}は周期超関数の Fourier 係数ではない。実際、k=(j,0,…,0)k=(j,0,\ldots,0)とすれば、任意のm∈N≥0m\in\Nに対してej/(1+j)m→∞e^j/(1+j)^m\to\inftyとなり、必要な多項式増大の評価に反する。

2 周期 Sobolev 空間と Laplace 作用素

定義 2.1 (周期 Sobolev 空間).n∈N≥1n\in\NN、m∈N≥0m\in\Nとする。Hm(Tn)H^m(\mathbb T^n)を、∣α∣≤m|\alpha|\leq mを満たすすべての多重指数α\alphaについて∂αT\partial^\alpha TがL2(Tn)L^2(\mathbb T^n)の関数で表される周期超関数TTの空間と定め、この空間を 周期 Sobolev 空間 (periodic Sobolev space) という。ノルムは

∥T∥Hm2=∑∣α∣≤m∥∂αT∥22\|T\|_{H^m}^2=\sum_{|\alpha|\leq m}\|\partial^\alpha T\|_2^2

とする。H0(Tn)=L2(Tn)H^0(\mathbb T^n)=L^2(\mathbb T^n)である。

定理 2.2.n∈N≥1n\in\NNとし、ek(x)=e2πik⋅xe_k(x)=e^{2\pi i k\cdot x}、λk=4π2∣k∣2\lambda_k=4\pi^2|k|^2とおく。また、m∈N≥0m\in\Nと周期超関数TTをとる。このとき、

∑k∈Zn(1+λk)m∣T^(k)∣2<∞\sum_{k\in\Z^n}(1+\lambda_k)^m|\widehat T(k)|^2<\infty

であることは、T∈Hm(Tn)T\in H^m(\mathbb T^n)であるための必要十分条件である。m,nm,nのみに依存する正の定数c,Cc,Cが存在し、c∥T∥Hm≤(∑k(1+λk)m∣T^(k)∣2)1/2≤C∥T∥Hmc\|T\|_{H^m}\leq(\sum_k(1+\lambda_k)^m|\widehat T(k)|^2)^{1/2}\leq C\|T\|_{H^m}が成り立つ。

L2(Tn)L^2(\mathbb T^n)上の作用素

Au=−Δu,D(A)=H2(Tn)Au=-\Delta u,\qquad D(A)=H^2(\mathbb T^n)

は非負自己共役作用素であり、(ek)k∈Zn(e_k)_{k\in\Z^n}はその完全正規直交固有系、固有値はλk\lambda_kである。また、

D(A)={u∈L2:∑kλk2∣u^(k)∣2<∞},Au=∑kλku^(k)ekD(A)=\left\{u\in L^2:\sum_k\lambda_k^2|\widehat u(k)|^2<\infty\right\}, \qquad Au=\sum_k\lambda_k\widehat u(k)e_k

であり、作用素の級数はL2L^2で収束する。形式域はD(A1/2)=H1(Tn)D(A^{1/2})=H^1(\mathbb T^n)であり、第一変数について線形な形式は

q(u,v)=∫Tn∇u⋅∇v‾=∑kλku^(k)v^(k)‾(u,v∈H1(Tn))q(u,v)=\int_{\mathbb T^n}\nabla u\cdot\overline{\nabla v} =\sum_k\lambda_k\widehat u(k)\overline{\widehat v(k)} \qquad(u,v\in H^1(\mathbb T^n))

で与えられる。HmH^mの Fourier 部分和はHmH^mで収束する。

証明.T∈HmT\in H^mならば、定理 1.3の微分の係数公式と命題 1.2の Parseval の等式から

∥T∥Hm2=∑kWm(k)∣T^(k)∣2,Wm(k)=∑∣α∣≤m∣(2πk)α∣2\|T\|_{H^m}^2=\sum_k W_m(k)|\widehat T(k)|^2, \qquad W_m(k)=\sum_{|\alpha|\leq m}|(2\pi k)^\alpha|^2

を得る。Wm(k)W_m(k)の各項はC(1+4π2∣k∣2)mC(1+4\pi^2|k|^2)^m以下である。逆向きには、11と各(2πkj)2m(2\pi k_j)^{2m}がWm(k)W_m(k)に含まれることから

cm,n(1+4π2∣k∣2)m≤Wm(k)≤Cm,n(1+4π2∣k∣2)mc_{m,n}(1+4\pi^2|k|^2)^m\leq W_m(k)\leq C_{m,n}(1+4\pi^2|k|^2)^m

となる。m=0m=0の場合には両辺の重みは11である。

重み付き平方和が有限ならば、各列((2πik)αT^(k))k((2\pi i k)^\alpha\widehat T(k))_kは∣α∣≤m|\alpha|\leq mに対して平方可算である。完全正規直交系により、この列を係数にもつL2L^2関数が存在する。定理 1.3の一意性によって、その関数は∂αT\partial^\alpha Tに一致する。従ってT∈HmT\in H^mである。ノルムの等式の末尾を評価すればHmH^m収束も従う。

Fourier 係数への写像はL2(Tn)L^2(\mathbb T^n)からℓ2(Zn)\ell^2(\Z^n)へのユニタリ写像である。この写像の下で−Δ-\Deltaはak↦λkaka_k\mapsto\lambda_ka_kとなり、その最大定義域は∑kλk2∣ak∣2<∞\sum_k\lambda_k^2|a_k|^2<\inftyである。a∈ℓ2a\in\ell^2の下では、この条件は∑k(1+λk)2∣ak∣2<∞\sum_k(1+\lambda_k)^2|a_k|^2<\inftyと同値であるから、作用素域はH2H^2となる。

対角作用素をMMと書く。b∈D(M∗)b\in D(M^*)ならば、各座標の単位ベクトルを随伴の定義へ代入することによって(M∗b)k=λkbk(M^*b)_k=\lambda_kb_kを得る。従って(λkbk)k∈ℓ2(\lambda_kb_k)_k\in\ell^2であり、b∈D(M)b\in D(M)となる。逆にb∈D(M)b\in D(M)ならば Cauchy–Schwarz の不等式により(Ma,b)=(a,Mb)(Ma,b)=(a,Mb)がすべてのa∈D(M)a\in D(M)について成り立つので、b∈D(M∗)b\in D(M^*)である。従ってM=M∗M=M^*であり、(Ma,a)=∑kλk∣ak∣2≥0(Ma,a)=\sum_k\lambda_k|a_k|^2\geq0である。

非負の対角作用素の平方根は対角成分λk\sqrt{\lambda_k}をもち、その定義域は∑kλk∣ak∣2<∞\sum_k\lambda_k|a_k|^2<\inftyである。a∈ℓ2a\in\ell^2の下では、この条件はH1H^1の重み付き条件と同値である。各一次微分への Parseval の等式からqqの積分表示と級数表示が一致する。形式の級数の絶対収束は、∑kλk∣u^(k)v^(k)∣≤q(u,u)1/2q(v,v)1/2\sum_k\lambda_k|\widehat u(k)\widehat v(k)|\leq q(u,u)^{1/2}q(v,v)^{1/2}による。▨

3 矩形とトーラスの Poisson 問題

定理 3.1.n∈N≥1n\in\NN、L1,…,Ln>0L_1,\ldots,L_n>0、Q=∏i=1n(0,Li)Q=\prod_{i=1}^n(0,L_i)とする。AD,ANA_D,A_NをQQの Dirichlet・Neumann の非負 Laplace 作用素とする。このとき、それぞれの正規直交固有系は

ekD(x)=∏i2Lisin⁡πkixiLi(k∈N≥1n),ekN(x)=∏ici,ki(xi)(k∈N≥0n),e_k^D(x)=\prod_i\sqrt{\frac2{L_i}}\sin\frac{\pi k_ix_i}{L_i} \quad(k\in\NN^n), \qquad e_k^N(x)=\prod_i c_{i,k_i}(x_i)\quad(k\in\N^n),ci,0=Li−1/2,ci,j(xi)=2Licos⁡πjxiLi(j≥1),λk=π2∑iki2Li2c_{i,0}=L_i^{-1/2},\qquad c_{i,j}(x_i)=\sqrt{\frac2{L_i}}\cos\frac{\pi jx_i}{L_i}\quad(j\geq1), \qquad \lambda_k=\pi^2\sum_i\frac{k_i^2}{L_i^2}

である。B∈{D,N}B\in\{D,N\}とu∈L2(Q)u\in L^2(Q)に対してuk=(u,ekB)2u_k=(u,e_k^B)_2とおくと、

u∈D(AB)⟺∑kλk2∣uk∣2<∞,u∈D(AB1/2)⟺∑kλk∣uk∣2<∞.u\in D(A_B)\Longleftrightarrow\sum_k\lambda_k^2|u_k|^2<\infty, \qquad u\in D(A_B^{1/2})\Longleftrightarrow\sum_k\lambda_k|u_k|^2<\infty.

形式域はB=DB=DならばH01(Q)H_0^1(Q)、B=NB=NならばH1(Q)H^1(Q)である。

任意のf∈L2(Q)f\in L^2(Q)に対してADu=fA_Du=fの解はただ一つ存在し、

u=∑k∈N≥1nfkλkekDu=\sum_{k\in\NN^n}\frac{f_k}{\lambda_k}e_k^D

で与えられる。ANu=fA_Nu=fが解をもつための必要十分条件は∫Qf=0\int_Q f=0である。この条件の下で平均零の解はただ一つ存在し、

u=∑k∈N≥0n∖{0}fkλkekNu=\sum_{k\in\N^n\setminus\{0\}}\frac{f_k}{\lambda_k}e_k^N

で与えられる。どちらの解の級数も作用素域のグラフノルム∥u∥2+∥ABu∥2\|u\|_2+\|A_Bu\|_2で収束する。Neumann 問題のすべての解は平均零の解に定数を加えたものである。

証明. 固有系と固有値は§E15.19 定理 3.1による。作用素域と形式域の係数表示は§E15.19 定理 1.3による。u∈L2u\in L^2の下では∑k(1+λk)∣uk∣2<∞\sum_k(1+\lambda_k)|u_k|^2<\inftyと∑kλk∣uk∣2<∞\sum_k\lambda_k|u_k|^2<\inftyは同値である。

ABu=fA_Bu=fの各固有関数に沿う係数はλkuk=fk\lambda_ku_k=f_kである。Dirichlet の場合はλk≥π2∑iLi−2>0\lambda_k\geq\pi^2\sum_iL_i^{-2}>0なのでuk=fk/λku_k=f_k/\lambda_kと一意に決まる。Parseval の等式により

∑kλk2∣uk∣2=∑k∣fk∣2=∥f∥22,∑k∣uk∣2≤(π2∑iLi−2)−2∥f∥22.\sum_k\lambda_k^2|u_k|^2=\sum_k|f_k|^2=\|f\|_2^2, \qquad \sum_k|u_k|^2\leq (\pi^2\sum_iL_i^{-2})^{-2}\|f\|_2^2.

従って級数は作用素域の元を定め、ADu=fA_Du=fを満たす。

Neumann の零固有値はk=0k=0だけであり、e0N=∣Q∣−1/2e_0^N=|Q|^{-1/2}である。係数方程式がk=0k=0で成立する条件はf0=∣Q∣−1/2∫Qf=0f_0=|Q|^{-1/2}\int_Q f=0である。k≠0k\ne0ならばλk≥π2min⁡iLi−2>0\lambda_k\geq\pi^2\min_i L_i^{-2}>0であるから、u0=0u_0=0、uk=fk/λku_k=f_k/\lambda_kと定めると、Dirichlet の場合と同じ二つの平方和評価を得る。従って平均零の解が存在する。二つの解の差の係数は零モード以外で零となるので、その差は定数である。各平方和の末尾は零へ収束するため、両級数はグラフノルムで収束する。▨

命題 3.2.n∈N≥1n\in\NN、f∈L2(Tn)f\in L^2(\mathbb T^n)とする。このとき、−Δu=f-\Delta u=fがH2(Tn)H^2(\mathbb T^n)の解をもつための必要十分条件はf^(0)=∫Tnf=0\widehat f(0)=\int_{\mathbb T^n}f=0である。この条件の下で平均零の解は

u=∑k≠0f^(k)4π2∣k∣2eku=\sum_{k\ne0}\frac{\widehat f(k)}{4\pi^2|k|^2}e_k

と一意に与えられる。級数はH2H^2で収束し、すべての解はこの解に定数を加えたものである。

証明.定理 2.2により係数方程式は4π2∣k∣2u^(k)=f^(k)4\pi^2|k|^2\widehat u(k)=\widehat f(k)である。k=0k=0の方程式が必要条件を与える。k≠0k\ne0に対する4π2∣k∣2≥4π24\pi^2|k|^2\geq4\pi^2と Parseval の等式から

∑k≠0(1+4π2∣k∣2)2∣f^(k)4π2∣k∣2∣2≤(1+(4π2)−1)2∥f∥22\sum_{k\ne0}(1+4\pi^2|k|^2)^2 \left|\frac{\widehat f(k)}{4\pi^2|k|^2}\right|^2 \leq (1+(4\pi^2)^{-1})^2\|f\|_2^2

を得る。従って与えた級数はH2H^2の解を定める。係数方程式は零モード以外のすべての係数を一意に定めるから、解の差は定数である。▨

定理 3.3.n∈N≥1n\in\NN、f∈C∞(Tn)f\in C^\infty(\mathbb T^n)とする。このとき、Δu=f\Delta u=fがC∞(Tn)C^\infty(\mathbb T^n)の解をもつための必要十分条件は∫Tnf=0\int_{\mathbb T^n}f=0である。この条件の下で平均零の解は

u^(0)=0,u^(k)=−f^(k)4π2∣k∣2(k≠0)\widehat u(0)=0,\qquad \widehat u(k)=-\frac{\widehat f(k)}{4\pi^2|k|^2}\quad(k\ne0)

と一意に与えられ、その Fourier 級数は各階の導関数まで一様収束する。すべての滑らかな解はこの解に定数を加えたものである。

証明.命題 3.2を外力−f-fに適用すると、可解条件、平均零のH2H^2解の係数表示、一意性を得る。命題 1.2によりf^(k)\widehat f(k)は急減少するので、与えたu^(k)\widehat u(k)も急減少する。同じ命題により、その Fourier 級数は各階の導関数まで一様収束し、滑らかな解を与える。滑らかな解はH2H^2に属するから、一意性も既に得た結論から従う。▨

定理 3.4.n∈N≥1n\in\NN、μ>0\mu>0、f∈C∞(Tn)f\in C^\infty(\mathbb T^n)とする。このとき(Δ−μ)u=f(\Delta-\mu)u=fの滑らかな周期解はただ一つ存在し、

u^(k)=−f^(k)4π2∣k∣2+μ(k∈Zn)\widehat u(k)=-\frac{\widehat f(k)}{4\pi^2|k|^2+\mu}\qquad(k\in\Z^n)

で与えられる。その Fourier 級数は各階の導関数まで一様収束する。特にu^(0)=−f^(0)/μ\widehat u(0)=-\widehat f(0)/\muである。

証明. 分母はすべてのkkについてμ\mu以上であるから、命題 1.2により与えた係数列は急減少する。その級数は滑らかな関数uuを定め、二階までの項別微分により(Δ−μ)u(\Delta-\mu)uの係数はf^(k)\widehat f(k)となる。Fourier 係数の一意性から方程式が成立する。任意の滑らかな解の係数は−(4π2∣k∣2+μ)u^(k)=f^(k)-(4\pi^2|k|^2+\mu)\widehat u(k)=\widehat f(k)を満たすので、解は一意である。▨

例 3.5.n∈N≥1n\in\NNとする。このとき、Δu=1\Delta u=1は周期超関数の解さえもたない。実際、微分の係数公式からΔu^(0)=0\widehat{\Delta u}(0)=0であるが、定数関数11の零次係数は11である。μ>0\mu>0ならば(Δ−μ)u=1(\Delta-\mu)u=1の唯一の滑らかな周期解は定数関数u=−1/μu=-1/\muである。

k∈Zn∖{0}k\in\Z^n\setminus\{0\}とし、f(x)=cos⁡(2πk⋅x)f(x)=\cos(2\pi k\cdot x)とおく。このとき、f^(k)=f^(−k)=1/2\widehat f(k)=\widehat f(-k)=1/2であり、他の係数は零である。従って−Δu=f-\Delta u=fの平均零の解はu=f/(4π2∣k∣2)u=f/(4\pi^2|k|^2)であり、Δu=f\Delta u=fの平均零の解はu=−f/(4π2∣k∣2)u=-f/(4\pi^2|k|^2)である。

例 3.6.b>0b>0、m∈N≥1m\in\NNとする。このとき、矩形Q=(0,π)×(0,b)Q=(0,\pi)\times(0,b)において、上辺の境界値がsin⁡(mx)\sin(mx)、他の三辺の境界値が零である調和関数は

u(x,y)=sin⁡(mx)sinh⁡(my)sinh⁡(mb)u(x,y)=\sin(mx)\frac{\sinh(my)}{\sinh(mb)}

と一意に与えられる。実際、uxx=−m2uu_{xx}=-m^2u、uyy=m2uu_{yy}=m^2uであるからΔu=0\Delta u=0となる。sin⁡(0)=sin⁡(mπ)=0\sin(0)=\sin(m\pi)=0とsinh⁡(0)=0\sinh(0)=0により三辺でu=0u=0であり、y=by=bではu=sin⁡(mx)u=\sin(mx)である。uuは閉矩形上で滑らかである。同じ境界値をもつ二つのC2(Q)∩C(Q‾)C^2(Q)\cap C(\overline Q)の実調和関数に、§E15.18 系 1.2をL=−ΔL=-\Deltaとして両方向に適用すると両者は一致する。複素数値の解の場合は差の実部と虚部に同じ議論を適用する。

4 円板の Poisson 級数

定理 4.1.D={z∈C:∣z∣<1}D=\{z\in\C:|z|<1\}とし、g∈C(T;C)g\in C(\mathbb T;\C)を円周上の境界値g(e2πiθ)=g(θ)g(e^{2\pi i\theta})=g(\theta)とみなす。このとき、唯一のC∞(D)∩C(D‾)C^\infty(D)\cap C(\overline D)の調和拡張は

u(re2πiθ)=∑k∈Zg^(k)r∣k∣e2πikθ(0≤r<1)u(re^{2\pi i\theta})=\sum_{k\in\Z}\widehat g(k)r^{|k|}e^{2\pi i k\theta} \qquad(0\leq r<1)

である。r=0r=0でのk=0k=0の項はg^(0)\widehat g(0)とする。任意の0≤ρ<10\leq\rho<1に対して、この級数は∣z∣≤ρ|z|\leq\rho上で実変数の各階の導関数まで絶対一様収束する。さらに

u(re2πiθ)=∫011−r21−2rcos⁡(2π(θ−t))+r2 g(t) dt,sup⁡θ∈T∣u(re2πiθ)−g(θ)∣⟶0(r↑1)u(re^{2\pi i\theta})=\int_0^1 \frac{1-r^2}{1-2r\cos(2\pi(\theta-t))+r^2}\,g(t)\,dt, \qquad \sup_{\theta\in\mathbb T}|u(re^{2\pi i\theta})-g(\theta)|\longrightarrow0 \quad(r\uparrow1)

が成り立つ。

証明.∣g^(k)∣≤∥g∥∞|\widehat g(k)|\leq\|g\|_\inftyである。z=re2πiθz=re^{2\pi i\theta}とおくと、k≥1k\geq1の項のrke2πikθr^ke^{2\pi ik\theta}はzkz^k、k≤−1k\leq-1の項はz‾∣k∣\overline z^{|k|}である。実変数のqq階の任意の導関数について、k≥qk\geq qの場合にzkz^kとz‾k\overline z^kの導関数の絶対値はkqρk−qk^q\rho^{k-q}以下となる。0<ρ<10<\rho<1ならば∑k≥qkqρk−q<∞\sum_{k\geq q}k^q\rho^{k-q}<\inftyであるから、各階の導関数の級数は閉じた部分円板上で絶対一様収束する。ρ=0\rho=0の場合は正のρ\rhoの場合に含まれる。§E15.14 例 3.2により各項は調和多項式であるから、級数の和はDD上で滑らかな調和関数となる。

0≤r<10\leq r<1とs∈Ts\in\mathbb Tに対して、二つの幾何級数を足すと

∑k∈Zr∣k∣e2πiks=1+re2πis1−re2πis+re−2πis1−re−2πis=1−r21−2rcos⁡(2πs)+r2\sum_{k\in\Z}r^{|k|}e^{2\pi iks} =1+\frac{re^{2\pi is}}{1-re^{2\pi is}} +\frac{re^{-2\pi is}}{1-re^{-2\pi is}} =\frac{1-r^2}{1-2r\cos(2\pi s)+r^2}

となる。固定したr<1r<1での絶対一様収束により、s=θ−ts=\theta-tを代入し、g(t)g(t)を掛けて積分することができる。その積分は主張の級数に一致する。

§E15.29 定理 5.1のn=2n=2、R=1R=1において、z=re2πiθz=re^{2\pi i\theta}、ζ=e2πit\zeta=e^{2\pi it}とおく。σ1=2π\sigma_1=2\pi、dσ(ζ)=2π dtd\sigma(\zeta)=2\pi\,dtであるから

P1(z,ζ) dσ(ζ)=1−r22π∣z−ζ∣2 2π dt=1−r21−2rcos⁡(2π(θ−t))+r2 dt.P_1(z,\zeta)\,d\sigma(\zeta) =\frac{1-r^2}{2\pi|z-\zeta|^2}\,2\pi\,dt =\frac{1-r^2}{1-2r\cos(2\pi(\theta-t))+r^2}\,dt.

従って級数の和は既出の Poisson 積分と一致し、閉円板へ境界値ggをもって連続に延長され、調和拡張は一意である。閉円板上の連続関数uuは一様連続であるため、∣re2πiθ−e2πiθ∣=1−r|re^{2\pi i\theta}-e^{2\pi i\theta}|=1-rから主張の境界への一様収束が従う。▨

注意 4.2.m∈N≥1m\in\NNとし、θ∈T\theta\in\mathbb Tに対してg(θ)=cos⁡(2πmθ)g(\theta)=\cos(2\pi m\theta)とする。このとき、g^(m)=g^(−m)=1/2\widehat g(m)=\widehat g(-m)=1/2であり、他の係数は零である。z=re2πiθz=re^{2\pi i\theta}(0≤r<1)(0\leq r<1)とおくと、ggの Poisson 級数は、§E15.14 例 3.2の角度変数を2πθ2\pi\thetaとした正負mm次の持ち上げの平均rmcos⁡(2πmθ)=ℜ(zm)r^m\cos(2\pi m\theta)=\Re(z^m)に一致する。定数境界値11の Fourier 係数は零次だけが11で他は零であるから、その Poisson 級数は同じ例の零次の持ち上げu=1u=1に一致する。

例 4.3.0<a<10<a<1とし、円周上の関数ggを

g(θ)=1−a21−2acos⁡(2πθ)+a2g(\theta)=\frac{1-a^2}{1-2a\cos(2\pi\theta)+a^2}

と定める。このとき、二つの幾何級数を和にすると

g(θ)=1+∑k=1∞ak(e2πikθ+e−2πikθ)g(\theta)=1+\sum_{k=1}^\infty a^k \bigl(e^{2\pi ik\theta}+e^{-2\pi ik\theta}\bigr)

であり、級数は絶対一様収束する。項別積分によりg^(k)=a∣k∣\widehat g(k)=a^{|k|}(k∈Z)(k\in\Z)であるから、ggは無限個の非零 Fourier 係数をもつ。定理 4.1により調和拡張は、0≤r<10\leq r<1で

u(re2πiθ)=1+∑k=1∞(ar)k(e2πikθ+e−2πikθ)=1−a2r21−2arcos⁡(2πθ)+a2r2u(re^{2\pi i\theta}) =1+\sum_{k=1}^\infty (ar)^k \bigl(e^{2\pi ik\theta}+e^{-2\pi ik\theta}\bigr) =\frac{1-a^2r^2}{1-2ar\cos(2\pi\theta)+a^2r^2}

となる。最後の等式は二つの幾何級数を和にしたものである。原点ではu(0)=1u(0)=1である。分母は0≤r≤10\leq r\leq1において(1−a)2(1-a)^2以上なので、u(re2πiθ)u(re^{2\pi i\theta})はr↑1r\uparrow1のときg(θ)g(\theta)にθ\thetaについて一様収束する。

前提記事

9 本の記事・単元を表示