§E12.46変分的発展方程式

最終更新

線形発展方程式では、初期値から時間とともに変化する状態を定める。係数が時間に依存すると、一般には一つの生成作用素の半群による表示を用いることができず、係数に跳びがある場合には解の時間微分も連続とは限らない。変分的発展方程式は、解と時間微分を異なる空間で評価し、このような問題の弱解を扱う基本的な枠組みである。本記事では、時間依存形式に対する変分解の基本的な性質と一意存在を解説する。

1 Hilbert 三つ組と解空間

定義 1.1.K∈{R,C}\mathbb K\in\{\R,\C\}とし、V,HV,HをK\mathbb K上の可分 Hilbert 空間、j:V→Hj:V\to Hを像が稠密な連続単射線形写像とする。内積は第一変数について線形とする。V×V^\timesは実の場合には連続双対、複素の場合には連続反双対とし、双対ノルムを入れる。§E12.43 定義 1.1の写像

I:H⟶V×,(Ih)(v)=⟨h,jv⟩HI:H\longrightarrow V^\times,\qquad (Ih)(v)=\langle h,jv\rangle_H

によってV⊂H⊂V×V\subset H\subset V^\timesと同一視し、この三つ組を Hilbert 三つ組という。cj=∥j∥c_j=\|j\|とおくと、§E12.43 命題 1.2により包含は連続な単射であり、

∥v∥H≤cj∥v∥V,∥h∥V×≤cj∥h∥H\|v\|_H\leq c_j\|v\|_V,\qquad \|h\|_{V^\times}\leq c_j\|h\|_H

である。

T>0T>0に対し、変分的発展方程式の解空間 (variational evolution space) を

W(0,T)={u∈L2(0,T;V):u′∈L2(0,T;V×)},∥u∥W2=∥u∥L2(0,T;V)2+∥u′∥L2(0,T;V×)2W(0,T)=\{u\in L^2(0,T;V):u'\in L^2(0,T;V^\times)\},\qquad \|u\|_W^2=\|u\|_{L^2(0,T;V)}^2+\|u'\|_{L^2(0,T;V^\times)}^2

と定める。ここでu′u'はV×V^\times値の弱微分であり、すべてのφ∈Cc∞(0,T)\varphi\in C_c^\infty(0,T)について

∫0Tu(t)φ′(t) dt=−∫0Tu′(t)φ(t) dtが V× で成り立つ\int_0^T u(t)\varphi'(t)\,dt=-\int_0^T u'(t)\varphi(t)\,dt \quad\text{が }V^\times\text{ で成り立つ}

という意味である。弱微分の定義は§E10.48 定義 3.1と同じであり、W(0,T)⊂H1(0,T;V×)W(0,T)\subset H^1(0,T;V^\times)である。

命題 1.2.V⊂H⊂V×V\subset H\subset V^\timesを実または複素の可分 Hilbert 三つ組、T>0T>0とする。このとき、W(0,T)W(0,T)はノルム∥⋅∥W\|\cdot\|_Wに関して Hilbert 空間である。また、写像

F⟼[u⟼∫0TF(t)(u(t)) dt]F\longmapsto\left[u\longmapsto\int_0^T F(t)(u(t))\,dt\right]

はL2(0,T;V×)L^2(0,T;V^\times)からL2(0,T;V)L^2(0,T;V)の実双対または複素反双対への線形等長同型である。

証明.C:V×→V∗C:V^\times\to V^*を(Cℓ)(v)=ℓ(v)‾(C\ell)(v)=\overline{\ell(v)}と定める。ただし実数体では複素共役を恒等写像とする。CCは全単射等長写像であり、複素数体では共役線形である。§E10.5 系 1.2の Riesz 対応とC−1C^{-1}の合成は

R×:V⟶V×,R×v=⟨v,⋅⟩VR_\times:V\longrightarrow V^\times,\qquad R_\times v=\langle v,\cdot\rangle_V

という線形等長同型である。この写像で内積を移すとV×V^\timesは Hilbert 空間となる。§E10.25 定理 5.1をVVとV×V^\timesに適用すると、二つのL2L^2空間は Hilbert 空間である。同定理の通常双対への積分表示を、入力のFFにはCC、出力の汎関数には複素共役を施して移すと、主張の線形等長同型を得る。二つの共役線形写像の合成なので、複素数体の場合にも線形である。

un→uu_n\to uをL2(0,T;V)L^2(0,T;V)で、un′→gu_n'\to gをL2(0,T;V×)L^2(0,T;V^\times)で満たす列をとる。包含V→V×V\to V^\timesのノルムはcj2c_j^2以下であるから、任意のφ∈Cc∞(0,T)\varphi\in C_c^\infty(0,T)に対して

∥∫0T(un−u)φ′∥V×≤cj2∥un−u∥L2V∥φ′∥2⟶0,∥∫0T(un′−g)φ∥V×≤∥un′−g∥L2V×∥φ∥2⟶0.\left\|\int_0^T(u_n-u)\varphi'\right\|_{V^\times} \leq c_j^2\|u_n-u\|_{L^2V}\|\varphi'\|_2\longrightarrow0, \qquad \left\|\int_0^T(u_n'-g)\varphi\right\|_{V^\times} \leq\|u_n'-g\|_{L^2V^\times}\|\varphi\|_2\longrightarrow0.

したがってu′=gu'=gであり、弱微分のグラフは二つのL2L^2空間の積の閉線形部分空間である。写像u↦(u,u′)u\mapsto(u,u')によってその積内積を引き戻すと、∥⋅∥W\|\cdot\|_Wを誘導する完備な内積を得る。▨

2 連続代表と部分積分

補題 2.1.V⊂H⊂V×V\subset H\subset V^\timesを実または複素の可分 Hilbert 三つ組、T>0T>0とする。このとき、任意のu∈W(0,T)u\in W(0,T)に対して、[0,T][0,T]の近傍で定義された滑らかなVV値関数の列(un)(u_n)が存在して

un⟶uin L2(0,T;V),un′⟶u′in L2(0,T;V×)u_n\longrightarrow u\quad\text{in }L^2(0,T;V),\qquad u_n'\longrightarrow u'\quad\text{in }L^2(0,T;V^\times)

を満たす。

証明.u∈W(0,T)⊂H1(0,T;V×)u\in W(0,T)\subset H^1(0,T;V^\times)に§E10.48 定理 4.2を適用し、V×V^\times値連続代表とその端点値u(0),u(T)u(0),u(T)をとる。−T<t<2T-T<t<2Tにおいて

u^(t)={u(−t),−T<t<0,u(t),0≤t≤T,u(2T−t),T<t<2T,g^(t)={−u′(−t),−T<t<0,u′(t),0<t<T,−u′(2T−t),T<t<2T\widehat u(t)= \begin{cases} u(-t),&-T<t<0,\\ u(t),&0\leq t\leq T,\\ u(2T-t),&T<t<2T, \end{cases} \qquad \widehat g(t)= \begin{cases} -u'(-t),&-T<t<0,\\ u'(t),&0<t<T,\\ -u'(2T-t),&T<t<2T \end{cases}

とおく。V×V^\timesにおける積分表示を各区間で変数変換すると、すべての−T<s<t<2T-T<s<t<2Tについて

u^(t)−u^(s)=∫stg^(r) dr\widehat u(t)-\widehat u(s)=\int_s^t\widehat g(r)\,dr

となる。区間が00またはTTをまたぐ場合も、その点で両側のV×V^\times値が一致するので、積分を足せば同じ等式を得る。また、u^∈L2(−T,2T;V)\widehat u\in L^2(-T,2T;V)、g^∈L2(−T,2T;V×)\widehat g\in L^2(-T,2T;V^\times)である。したがって§E10.48 定理 4.2によりu^′=g^\widehat u'=\widehat gである。

χ∈Cc∞(−T,2T)\chi\in C_c^\infty(-T,2T)を[0,T][0,T]の近傍で11となるようにとり、χu^\chi\widehat uをR\Rに零延長してUUとする。§E10.48 定理 5.1の積の公式と、χ\chiが両端の近傍で零であることから、

U∈L2(R;V),U′=G∈L2(R;V×),U\in L^2(\R;V),\qquad U'=G\in L^2(\R;V^\times),

ただしGGはχ′u^+χg^\chi'\widehat u+\chi\widehat gの零延長である。非負の実数値核ρ∈Cc∞(−1,1)\rho\in C_c^\infty(-1,1)、∫ρ=1\int\rho=1に対し、ρε(t)=ε−1ρ(t/ε)\rho_\varepsilon(t)=\varepsilon^{-1}\rho(t/\varepsilon)、uε=ρε∗Uu_\varepsilon=\rho_\varepsilon*Uとおく。§E10.48 命題 1.1によりuε∈C∞(R;V)u_\varepsilon\in C^\infty(\R;V)であり、uε→Uu_\varepsilon\to UがL2(R;V)L^2(\R;V)で成り立つ。

§E10.18 命題 5.2を包含V→V×V\to V^\timesに適用すると、uεu_\varepsilonはV×V^\timesで計算した畳み込みにも一致する。U′=GU'=Gの弱微分の等式に、各ttについて試験関数s↦ρε(t−s)s\mapsto\rho_\varepsilon(t-s)を代入すると

uε′(t)=∫Rρε′(t−s)U(s) ds=∫Rρε(t−s)G(s) dsu_\varepsilon'(t)=\int_\R \rho_\varepsilon'(t-s)U(s)\,ds =\int_\R\rho_\varepsilon(t-s)G(s)\,ds

がV×V^\timesで成り立つ。§E10.48 命題 1.1をGGにも適用するとuε′→Gu_\varepsilon'\to GがL2(R;V×)L^2(\R;V^\times)で成り立つ。[0,T][0,T]上ではU=uU=u、G=u′G=u'であるから、ε↓0\varepsilon\downarrow0の列をとればよい。▨

定理 2.2.V⊂H⊂V×V\subset H\subset V^\timesを実または複素の可分 Hilbert 三つ組、T>0T>0とする。このとき、W(0,T)W(0,T)の各元はただ一つのHH値連続代表をもち、この代表をとる写像W(0,T)→C([0,T];H)W(0,T)\to C([0,T];H)は連続である。連続代表をuuと書くと、

sup⁡0≤t≤T∥u(t)∥H2≤1T∥u∥L2(0,T;H)2+2∥u′∥L2(0,T;V×)∥u∥L2(0,T;V).\sup_{0\leq t\leq T}\|u(t)\|_H^2 \leq \frac1T\|u\|_{L^2(0,T;H)}^2 +2\|u'\|_{L^2(0,T;V^\times)}\|u\|_{L^2(0,T;V)}.

また、t↦∥u(t)∥H2t\mapsto\|u(t)\|_H^2は絶対連続であり、すべての0≤s≤t≤T0\leq s\leq t\leq Tについて

∥u(t)∥H2−∥u(s)∥H2=2ℜ∫stu′(r)(u(r)) dr.\|u(t)\|_H^2-\|u(s)\|_H^2 =2\Re\int_s^t u'(r)(u(r))\,dr.

特に、ほとんどすべてのttでddt∥u(t)∥H2=2ℜ(u′(t)(u(t)))\frac{d}{dt}\|u(t)\|_H^2=2\Re\bigl(u'(t)(u(t))\bigr)である。

証明.[0,T][0,T]の近傍で滑らかなVV値関数wwをとり、内積を微分すると

ddt∥w(t)∥H2=2ℜ⟨w′(t),w(t)⟩H=2ℜ(w′(t)(w(t)))\frac{d}{dt}\|w(t)\|_H^2 =2\Re\langle w'(t),w(t)\rangle_H =2\Re\bigl(w'(t)(w(t))\bigr)

を得る。s∈[0,T]s\in[0,T]を∥w(s)∥H2≤T−1∥w∥L2H2\|w(s)\|_H^2\leq T^{-1}\|w\|_{L^2H}^2となるように選ぶ。上の等式をssから各ttへ積分し、積分の向きに応じて絶対値をとると、

sup⁡t∥w(t)∥H2≤T−1∥w∥L2H2+2∫0T∥w′(r)∥V×∥w(r)∥V dr≤T−1∥w∥L2H2+2∥w′∥L2V×∥w∥L2V.\sup_t\|w(t)\|_H^2 \leq T^{-1}\|w\|_{L^2H}^2+2\int_0^T\|w'(r)\|_{V^\times}\|w(r)\|_V\,dr \leq T^{-1}\|w\|_{L^2H}^2+2\|w'\|_{L^2V^\times}\|w\|_{L^2V}.

u∈W(0,T)u\in W(0,T)を補題 2.1の列unu_nで近似する。直前の評価をun−umu_n-u_mに適用し、∥un−um∥L2H≤cj∥un−um∥L2V\|u_n-u_m\|_{L^2H}\leq c_j\|u_n-u_m\|_{L^2V}を用いると、(un)(u_n)はC([0,T];H)C([0,T];H)の Cauchy 列である。その一様極限をu~\widetilde uとすると、L2(0,T;H)L^2(0,T;H)の極限の一意性によりu~=u\widetilde u=uがほとんど至る所で成り立つ。二つの連続代表の差が一点で零でなければ、ある相対近傍で差のノルムが正となり、ほとんど至る所での一致に反する。したがって連続代表は一意である。滑らかな関数の評価を極限へ通すと、主張の評価と連続埋め込みを得る。

双対の評価と Cauchy–Schwarz の不等式により

∥un′(un)−u′(u)∥L1≤∥un′−u′∥L2V×∥un∥L2V+∥u′∥L2V×∥un−u∥L2V⟶0.\|u_n'(u_n)-u'(u)\|_{L^1} \leq\|u_n'-u'\|_{L^2V^\times}\|u_n\|_{L^2V} +\|u'\|_{L^2V^\times}\|u_n-u\|_{L^2V}\longrightarrow0.

滑らかな関数の積分等式を極限へ通すと、すべてのs,ts,tについて主張の等式を得る。被積分関数はL1L^1に属するから、得られた積分等式により、ノルムの二乗は絶対連続であり、ほとんど至る所で主張の微分公式が成り立つ。▨

命題 2.3.V⊂H⊂V×V\subset H\subset V^\timesを実または複素の可分 Hilbert 三つ組、T>0T>0とし、u,v∈W(0,T)u,v\in W(0,T)をそれぞれのHH値連続代表と同一視する。このとき、すべての0≤s≤t≤T0\leq s\leq t\leq Tについて

⟨u(t),v(t)⟩H−⟨u(s),v(s)⟩H=∫st(u′(r)(v(r))+v′(r)(u(r))‾) dr.\langle u(t),v(t)\rangle_H-\langle u(s),v(s)\rangle_H =\int_s^t\left(u'(r)(v(r))+\overline{v'(r)(u(r))}\right)\,dr.

実数体の場合、複素共役は恒等写像とする。

証明.u,vu,vを補題 2.1により滑らかなVV値関数un,vnu_n,v_nでWWにおいて近似する。内積の積の微分公式を積分すると、

⟨un(t),vn(t)⟩H−⟨un(s),vn(s)⟩H=∫st(un′(r)(vn(r))+vn′(r)(un(r))‾) dr.\langle u_n(t),v_n(t)\rangle_H-\langle u_n(s),v_n(s)\rangle_H =\int_s^t\left(u_n'(r)(v_n(r))+\overline{v_n'(r)(u_n(r))}\right)\,dr.

定理 2.2により左辺は主張の左辺へ収束する。右辺の第一項については

∥un′(vn)−u′(v)∥1≤∥un′−u′∥L2V×∥vn∥L2V+∥u′∥L2V×∥vn−v∥L2V⟶0.\|u_n'(v_n)-u'(v)\|_1 \leq\|u_n'-u'\|_{L^2V^\times}\|v_n\|_{L^2V} +\|u'\|_{L^2V^\times}\|v_n-v\|_{L^2V}\longrightarrow0.

第二項にもuuとvvを入れ替えた同じ評価を適用すると、積分の極限を得る。▨

3 可測な時間依存形式

定義 3.1.V⊂H⊂V×V\subset H\subset V^\timesを実または複素の可分 Hilbert 三つ組、T>0T>0とする。関数a:(0,T)×V×V→Ka:(0,T)\times V\times V\to\mathbb Kに対し、次を時間依存形式の有界性と強圧性の条件とする。

  1. 各v,w∈Vv,w\in Vについてt↦a(t;v,w)t\mapsto a(t;v,w)は Lebesgue 可測である。
  2. 定数M,α>0M,\alpha>0、λ≥0\lambda\geq0と零集合N⊂(0,T)N\subset(0,T)が存在する。すべてのt∉Nt\notin Nについて、a(t;⋅,⋅)a(t;\cdot,\cdot)は実の場合には双線形、複素の場合には第一変数線形、第二変数共役線形であり、すべてのv,w∈Vv,w\in Vに対して ∣a(t;v,w)∣≤M∥v∥V∥w∥V,ℜa(t;v,v)≥α∥v∥V2−λ∥v∥H2|a(t;v,w)|\leq M\|v\|_V\|w\|_V,\qquad \Re a(t;v,v)\geq\alpha\|v\|_V^2-\lambda\|v\|_H^2 を満たす。後者の評価を Gårding の不等式という。

t∉Nt\notin NにおいてA(t):V→V×\mathcal A(t):V\to V^\timesを(A(t)v)(w)=a(t;v,w)(\mathcal A(t)v)(w)=a(t;v,w)で定め、NN上ではA(t)=0\mathcal A(t)=0とする。

補題 3.2.V⊂H⊂V×V\subset H\subset V^\timesを実または複素の可分 Hilbert 三つ組、T>0T>0とし、aaが定義 3.1の条件を満たすとする。このとき

(Au)(t)=A(t)u(t)(\mathcal Au)(t)=\mathcal A(t)u(t)

はL2(0,T;V)L^2(0,T;V)からL2(0,T;V×)L^2(0,T;V^\times)への有界線形作用素を定め、∥Au∥L2V×≤M∥u∥L2V\|\mathcal Au\|_{L^2V^\times}\leq M\|u\|_{L^2V}を満たす。

証明.§E10.3 定理 3.6によりVVの高々可算な完全正規直交系(ek)(e_k)をとる。有限次元の場合には以下の和をその次元で止め、零空間の場合には空和とする。固定したv∈Vv\in Vについて、Riesz 同型R×v=⟨v,⋅⟩VR_\times v=\langle v,\cdot\rangle_Vを用いると、t∉Nt\notin Nにおいて

A(t)v=lim⁡n→∞∑k=1na(t;v,ek)R×ekin V×.\mathcal A(t)v=\lim_{n\to\infty}\sum_{k=1}^n a(t;v,e_k)R_\times e_k \quad\text{in }V^\times.

各係数は可測であるから有限和は強可測であり、その極限も強可測である。NN上の零への変更はこの結論を保つ。

強可測なu:(0,T)→Vu:(0,T)\to Vに対し、可測単関数sns_nをsn(t)→u(t)s_n(t)\to u(t)がほとんど至る所で成り立つようにとる。各A(⋅)sn(⋅)\mathcal A(\cdot)s_n(\cdot)は前段の有限個の強可測関数を可測集合上で貼り合わせた関数である。また、t∉Nt\notin Nについて

∥A(t)sn(t)−A(t)u(t)∥V×≤M∥sn(t)−u(t)∥V⟶0.\|\mathcal A(t)s_n(t)-\mathcal A(t)u(t)\|_{V^\times} \leq M\|s_n(t)-u(t)\|_V\longrightarrow0.

したがってAu\mathcal Auは強可測であり、∥A(t)u(t)∥V×≤M∥u(t)∥V\|\mathcal A(t)u(t)\|_{V^\times}\leq M\|u(t)\|_Vを二乗積分すると主張の評価を得る。代表元を変更しても出力はほとんど至る所で変わらず、線形性は各A(t)\mathcal A(t)の線形性から従う。▨

命題 3.3.V⊂H⊂V×V\subset H\subset V^\timesを実または複素の可分 Hilbert 三つ組、T>0T>0とし、aaが定義 3.1の条件を定数M,α>0M,\alpha>0、λ≥0\lambda\geq0について満たすとする。u∈W(0,T)u\in W(0,T)がf∈L2(0,T;V×)f\in L^2(0,T;V^\times)に対してu′+Au=fu'+\mathcal Au=fを満たすとし、u0=u(0)∈Hu_0=u(0)\in Hとする。このとき、すべての0≤s≤t≤T0\leq s\leq t\leq Tについて

12∥u(t)∥H2+∫stℜa(r;u(r),u(r)) dr=12∥u(s)∥H2+ℜ∫stf(r)(u(r)) dr.\frac12\|u(t)\|_H^2+\int_s^t\Re a(r;u(r),u(r))\,dr =\frac12\|u(s)\|_H^2+\Re\int_s^t f(r)(u(r))\,dr.

また、T,M,α,λ,cjT,M,\alpha,\lambda,c_jだけに依存する定数CCが存在して

∥u∥C([0,T];H)+∥u∥L2(0,T;V)+∥u′∥L2(0,T;V×)≤C(∥u0∥H+∥f∥L2(0,T;V×))\|u\|_{C([0,T];H)}+\|u\|_{L^2(0,T;V)}+\|u'\|_{L^2(0,T;V^\times)} \leq C\left(\|u_0\|_H+\|f\|_{L^2(0,T;V^\times)}\right)

を満たす。同じ形式に対する二つの解の差にも、右辺を初期値と外力の差のノルムに替えた評価が成り立つ。

証明.定理 2.2の積分等式へu′=f−Auu'=f-\mathcal Auを代入すると、主張のエネルギー等式を得る。y(t)=∥u(t)∥H2y(t)=\|u(t)\|_H^2とおく。Gårding の不等式と

2∣f(t)(u(t))∣≤α∥u(t)∥V2+α−1∥f(t)∥V×22|f(t)(u(t))|\leq\alpha\|u(t)\|_V^2+\alpha^{-1}\|f(t)\|_{V^\times}^2

により、ほとんど至る所で

y′(t)+α∥u(t)∥V2≤2λy(t)+α−1∥f(t)∥V×2.y'(t)+\alpha\|u(t)\|_V^2 \leq2\lambda y(t)+\alpha^{-1}\|f(t)\|_{V^\times}^2.

非負のVVノルムの項を除いて積分し、§E13.7 定理 1.1を適用する。K=∥u0∥H2+α−1∥f∥L2V×2K=\|u_0\|_H^2+\alpha^{-1}\|f\|_{L^2V^\times}^2とおくと

sup⁡ty(t)≤e2λTK,α∥u∥L2V2≤K+2λTe2λTK.\sup_t y(t)\leq e^{2\lambda T}K, \qquad \alpha\|u\|_{L^2V}^2\leq K+2\lambda T e^{2\lambda T}K.

二つ目の評価は微分不等式を00からTTへ積分して得られる。補題 3.2により

∥u′∥L2V×≤∥f∥L2V×+M∥u∥L2V\|u'\|_{L^2V^\times}\leq\|f\|_{L^2V^\times}+M\|u\|_{L^2V}

であるから、三つのノルムの評価が従う。二つの解u,vu,vの差はu−v∈W(0,T)u-v\in W(0,T)と(u−v)′+A(u−v)=fu−fv(u-v)'+\mathcal A(u-v)=f_u-f_vを満たすので、同じ評価を適用する。▨

4 Galerkin 近似と一意存在

補題 4.1.V⊂H⊂V×V\subset H\subset V^\timesを実または複素の可分 Hilbert 三つ組、T>0T>0とし、aaが定義 3.1の条件を満たすとする。u0∈Hu_0\in H、f∈L2(0,T;V×)f\in L^2(0,T;V^\times)とする。このとき、増加する有限次元部分空間En⊂VE_n\subset Vとun∈H1(0,T;En)u_n\in H^1(0,T;E_n)が存在し、⋃nEn\bigcup_nE_nはVVで稠密であって、Pn:H→EnP_n:H\to E_nをHH直交射影として

un(0)=Pnu0,⟨un′(t),w⟩H+a(t;un(t),w)=f(t)(w)(w∈En)u_n(0)=P_nu_0,\qquad \langle u_n'(t),w\rangle_H+a(t;u_n(t),w)=f(t)(w) \quad(w\in E_n)

がほとんど至る所で成り立つ。さらに、nnに依存しない定数CCに対して

sup⁡t∥un(t)∥H+∥un∥L2(0,T;V)≤C(∥u0∥H+∥f∥L2(0,T;V×)).\sup_t\|u_n(t)\|_H+\|u_n\|_{L^2(0,T;V)} \leq C\left(\|u_0\|_H+\|f\|_{L^2(0,T;V^\times)}\right).

証明.VVの可算稠密列をとり、HHの内積による Gram–Schmidt の操作を順に施す。零になる剰余を除くと、HHで正規直交する列(wk)(w_k)を得る。各操作は既に選んだ有限個のベクトルの線形結合を引くものであるから、各段階で元の列と得られたベクトルの有限線形包は一致する。したがってEn=span⁡{w1,…,wn}E_n=\operatorname{span}\{w_1,\ldots,w_n\}の合併はVVで稠密であり、HHでも稠密である。列が有限で止まる場合には、その後のEnE_nをすべてV=HV=Hとする。V={0}V=\{0\}の場合にはH={0}H=\{0\}であり、un=0u_n=0とすればよい。

un(t)=∑k=1dnck(t)wku_n(t)=\sum_{k=1}^{d_n}c_k(t)w_k、dn=dim⁡End_n=\dim E_nとおく。係数方程式は

ci′(t)+∑k=1dna(t;wk,wi)ck(t)=f(t)(wi),ci(0)=⟨u0,wi⟩H(1≤i≤dn)c_i'(t)+\sum_{k=1}^{d_n}a(t;w_k,w_i)c_k(t)=f(t)(w_i), \qquad c_i(0)=\langle u_0,w_i\rangle_H \quad(1\leq i\leq d_n)

である。係数行列は有界可測であり、右辺はL2(0,T;Kdn)L^2(0,T;\mathbb K^{d_n})に属する。§E13.7 系 4.2により[0,T][0,T]全体に絶対連続な解が存在する。係数ベクトルは閉区間上で連続なので有界であり、方程式からci′∈L2(0,T)c_i'\in L^2(0,T)となる。したがってun∈H1(0,T;En)u_n\in H^1(0,T;E_n)である。

EnE_nにVVの内積を入れた空間と、同じEnE_nにHHの内積を入れた空間から有限次元の Hilbert 三つ組を作る。形式の制限は同じM,α,λM,\alpha,\lambdaを満たし、包含のノルムはcjc_j以下である。外力の制限fn(t)=f(t)∣Enf_n(t)=f(t)|_{E_n}は∥fn∥L2En×≤∥f∥L2V×\|f_n\|_{L^2E_n^\times}\leq\|f\|_{L^2V^\times}を満たす。係数方程式に命題 3.3をこの三つ組で適用し、∥Pnu0∥H≤∥u0∥H\|P_nu_0\|_H\leq\|u_0\|_Hを用いると、主張の二つのノルムの一様評価を得る。▨

定理 4.2.V⊂H⊂V×V\subset H\subset V^\timesを実または複素の可分 Hilbert 三つ組、T>0T>0とし、aaが定義 3.1の条件を満たすとする。このとき、任意のu0∈Hu_0\in Hとf∈L2(0,T;V×)f\in L^2(0,T;V^\times)に対して

u(0)=u0,u′+Au=fin L2(0,T;V×)u(0)=u_0,\qquad u'+\mathcal Au=f\quad\text{in }L^2(0,T;V^\times)

を満たすu∈W(0,T)u\in W(0,T)がただ一つ存在する。初期値はHH値連続代表の値である。この解を形式aaに対する変分解 (variational solution) という。

証明.補題 4.1の(En,un)(E_n,u_n)をとる。L2(0,T;V)L^2(0,T;V)は§E10.25 定理 5.1により Hilbert 空間であり、§E12.32 例 2.5により回帰的である。unu_nはこの空間で有界であるから、§E12.35 定理 5.2により、unk⇀uu_{n_k}\rightharpoonup uを満たす部分列(unk)(u_{n_k})とu∈L2(0,T;V)u\in L^2(0,T;V)が存在する。

w∈Emw\in E_mとφ∈Cc∞(0,T)\varphi\in C_c^\infty(0,T)を固定する。nk≥mn_k\geq mのとき、Galerkin 方程式を部分積分すると

−∫0T⟨unk(t),w⟩Hφ′(t) dt+∫0Ta(t;unk(t),w)φ(t) dt=∫0Tf(t)(w)φ(t) dt.-\int_0^T\langle u_{n_k}(t),w\rangle_H\varphi'(t)\,dt +\int_0^T a(t;u_{n_k}(t),w)\varphi(t)\,dt =\int_0^T f(t)(w)\varphi(t)\,dt.

左辺の二項はunk∈L2Vu_{n_k}\in L^2Vの連続線形汎関数である。実際、その絶対値はそれぞれ

cj∥w∥H∥φ′∥2∥unk∥L2V,M∥w∥V∥φ∥2∥unk∥L2Vc_j\|w\|_H\|\varphi'\|_2\|u_{n_k}\|_{L^2V},\qquad M\|w\|_V\|\varphi\|_2\|u_{n_k}\|_{L^2V}

以下である。したがって弱極限をとると同じ等式がuuに対して成り立つ。⋃mEm\bigcup_mE_mはVVで稠密であり、各項はwwについて連続なので、等式はすべてのw∈Vw\in Vに拡張される。

補題 3.2によりg=f−Au∈L2(0,T;V×)g=f-\mathcal Au\in L^2(0,T;V^\times)とおくことができる。直前の等式は、すべてのw∈Vw\in Vについて

(∫0Tu(t)φ′(t) dt+∫0Tg(t)φ(t) dt)(w)=0\left(\int_0^T u(t)\varphi'(t)\,dt+\int_0^Tg(t)\varphi(t)\,dt\right)(w)=0

を意味するので、括弧内のV×V^\timesの元は零である。したがってu′=gu'=g、u∈W(0,T)u\in W(0,T)であり、定理 2.2によりuuはHH値連続代表をもつ。

w∈Emw\in E_mを固定し、今度はφ∈C∞([0,T];R)\varphi\in C^\infty([0,T];\R)、φ(T)=0\varphi(T)=0とする。nk≥mn_k\geq mでは⟨Pnku0,w⟩H=⟨u0,w⟩H\langle P_{n_k}u_0,w\rangle_H=\langle u_0,w\rangle_Hであるから、端点項を残した Galerkin 方程式の部分積分と同じ弱極限により

−∫0T⟨u(t),w⟩Hφ′(t) dt=∫0Tg(t)(w)φ(t) dt+⟨u0,w⟩Hφ(0)-\int_0^T\langle u(t),w\rangle_H\varphi'(t)\,dt =\int_0^T g(t)(w)\varphi(t)\,dt+\langle u_0,w\rangle_H\varphi(0)

となる。一方、v(t)=φ(t)wv(t)=\varphi(t)wに命題 2.3を適用すると、左辺は

∫0Tu′(t)(w)φ(t) dt+⟨u(0),w⟩Hφ(0)\int_0^T u'(t)(w)\varphi(t)\,dt+\langle u(0),w\rangle_H\varphi(0)

に等しい。u′=gu'=gかつφ(0)=1\varphi(0)=1とすることができるから、⟨u(0)−u0,w⟩H=0\langle u(0)-u_0,w\rangle_H=0である。⋃mEm\bigcup_mE_mのHHにおける稠密性からu(0)=u0u(0)=u_0を得る。

同じデータに対する二つの解の差に命題 3.3を適用すると、差のC([0,T];H)C([0,T];H)ノルムは零である。したがって解は一意である。▨

5 自律閉形式と半群

定理 5.1.HHを可分な複素 Hilbert 空間、qqをHH上の稠密に定義された下に有界な閉対称形式とする。q(x,x)≥m∥x∥H2q(x,x)\geq m\|x\|_H^2を満たすm∈Rm\in\Rを固定し、形式域V=D(q)V=D(q)に

⟨x,y⟩V=q(x,y)+(1−m)⟨x,y⟩H\langle x,y\rangle_V=q(x,y)+(1-m)\langle x,y\rangle_H

を入れる。VVはこの内積に関して可分であるとする。AAをqqに対応する自己共役作用素、T>0T>0とする。このとき、任意のu0∈Hu_0\in Hとf∈L2(0,T;H)f\in L^2(0,T;H)に対して

u(t)=e−tAu0+∫0te−(t−s)Af(s) ds(0≤t≤T)u(t)=e^{-tA}u_0+\int_0^t e^{-(t-s)A}f(s)\,ds\qquad(0\leq t\leq T)

はW(0,T)W(0,T)に属し、形式a(t;v,w)=q(v,w)a(t;v,w)=q(v,w)、外力IfIf、初期値u0u_0に対する変分解に一致する。

HHが実 Hilbert 空間でqqが実対称形式である場合にも、同じ仮定の下で同じ結論が成り立つ。この場合、e−tAe^{-tA}は複素化した形式に対応する半群の実部分空間への制限を表す。

証明.qqの形式ノルムに§E12.13 補題 1.2 (1)を適用すると、V→HV\to Hは連続であり、qqはVV上で有界である。また

q(v,v)=∥v∥V2−(1−m)∥v∥H2≥∥v∥V2−max⁡{1−m,0}∥v∥H2.q(v,v)=\|v\|_V^2-(1-m)\|v\|_H^2 \geq\|v\|_V^2-\max\{1-m,0\}\|v\|_H^2.

したがってqqは定義 3.1の条件を満たす。§E12.13 定理 2.1により、AAは自己共役で下に有界であり、D(A)D(A)はVVで稠密である。§E12.25 定理 2.1により、S(t)=e−tAS(t)=e^{-tA}は生成作用素−A-Aをもつ強連続半群であり、∥S(t)∥≤e−mt\|S(t)\|\leq e^{-mt}である。

x∈D(A)x\in D(A)、h∈C1([0,T];H)h\in C^1([0,T];H)をとる。§E12.26 定理 3.1 (2)を生成作用素−A-Aに適用すると、

z(t)=S(t)x+∫0tS(t−s)h(s) dsz(t)=S(t)x+\int_0^tS(t-s)h(s)\,ds

はz′+Az=hz'+Az=hを満たす古典解であり、z,z′,Az=h−z′z,z',Az=h-z'はHH値連続関数である。v∈D(A)v\in D(A)に対して

∥v∥V2=⟨Av,v⟩H+(1−m)∥v∥H2≤∥Av∥H∥v∥H+∣1−m∣∥v∥H2\|v\|_V^2=\langle Av,v\rangle_H+(1-m)\|v\|_H^2 \leq\|Av\|_H\|v\|_H+|1-m|\|v\|_H^2

である。この評価をz(t)−z(s)z(t)-z(s)に適用すると、zzはVV値連続でもある。したがってz∈L2Vz\in L^2V、z′∈L2H⊂L2V×z'\in L^2H\subset L^2V^\timesであり、q(z(t),v)=⟨Az(t),v⟩Hq(z(t),v)=\langle Az(t),v\rangle_Hによりz′+Az=Ihz'+\mathcal Az=Ihである。よってzzは変分解である。

xn∈D(A)x_n\in D(A)をxn→u0x_n\to u_0がHHで成り立つようにとる。ffをR\Rへ零延長し、§E10.48 命題 1.1による畳み込みを[0,T][0,T]に制限すると、hn∈C∞([0,T];H)h_n\in C^\infty([0,T];H)でhn→fh_n\to fがL2HL^2Hで成り立つ列を得る。前段の古典解をznz_nとすると、命題 3.3の差の評価と∥Ih∥V×≤cj∥h∥H\|Ih\|_{V^\times}\leq c_j\|h\|_Hにより(zn)(z_n)はWWの Cauchy 列である。命題 1.2によりzn→z∈Wz_n\to z\in Wであり、定理 2.2と補題 3.2により

z(0)=u0,z′+Az=If.z(0)=u_0,\qquad z'+\mathcal Az=If.

一方、L2(0,T;H)⊂L1(0,T;H)L^2(0,T;H)\subset L^1(0,T;H)であるから、§E12.26 定理 2.1の一様評価により、znz_nは主張の半群表示へC([0,T];H)C([0,T];H)で収束する。また、znz_nはWWでzzに収束するので、C([0,T];H)C([0,T];H)でもzzに収束する。したがって、二つの極限は一致する。定理 4.2の一意性により、表示された関数は変分解である。

H,qH,qが実である場合、HC=H+iHH_{\C}=H+iH、VC=V+iVV_{\C}=V+iVとし、

qC(x+iy,z+iw)=q(x,z)+q(y,w)+i(q(y,z)−q(x,w))q_{\C}(x+iy,z+iw)=q(x,z)+q(y,w)+i\bigl(q(y,z)-q(x,w)\bigr)

で複素化する。qC(x+iy,x+iy)=q(x,x)+q(y,y)q_{\C}(x+iy,x+iy)=q(x,x)+q(y,y)であり、形式ノルムの二乗も∥x∥V2+∥y∥V2\|x\|_V^2+\|y\|_V^2となる。したがってqCq_{\C}は同じ下界をもつ稠密に定義された閉対称形式であり、その形式域は可分である。第一表現定理の定義を実部と虚部に適用すると、対応作用素は

D(AC)=D(A)+iD(A),AC(x+iy)=Ax+iAyD(A_{\C})=D(A)+iD(A),\qquad A_{\C}(x+iy)=Ax+iAy

である。複素共役J(x+iy)=x−iyJ(x+iy)=x-iyはこの定義域を保ち、JACJ=ACJA_{\C}J=A_{\C}を満たす。

SC(t)=e−tACS_{\C}(t)=e^{-tA_{\C}}とおくと、JSC(t)JJS_{\C}(t)Jも複素線形な強連続半群である。その生成作用素は、t↓0t\downarrow0の差商にJJを作用させることにより−JACJ=−AC-JA_{\C}J=-A_{\C}となる。§E12.23 系 5.5によりJSC(t)J=SC(t)JS_{\C}(t)J=S_{\C}(t)であるから、SC(t)S_{\C}(t)はJJの固定部分空間HHを保つ。実データの半群表示もHH値である。複素の場合の結論により、その関数はL2VCL^2V_{\C}に属し、実部をとる有界写像を弱微分の等式へ作用させると、実のW(0,T)W(0,T)に属して実形式の変分方程式を満たす。実の場合の一意性により、実変分解との一致を得る。▨

6 例と強圧性の役割

例 6.1.K∈{R,C}\mathbb K\in\{\R,\C\}、T>0T>0とし、

H=ℓ2(N≥1;K),V={x∈H:∑n=1∞n2∣xn∣2<∞},⟨x,y⟩V=∑n=1∞n2xnyn‾H=\ell^2(\NN;\mathbb K),\qquad V=\left\{x\in H:\sum_{n=1}^\infty n^2|x_n|^2<\infty\right\},\qquad \langle x,y\rangle_V=\sum_{n=1}^\infty n^2x_n\overline{y_n}

とする。このとき、形式q(x,y)=⟨x,y⟩Vq(x,y)=\langle x,y\rangle_Vと外力f=0f=0に対し、任意のu0∈Hu_0\in Hから出る変分解は

un(t)=e−n2t(u0)nu_n(t)=e^{-n^2t}(u_0)_n

で与えられ、

∫0T∥u(t)∥V2 dt=12∑n=1∞(1−e−2n2T)∣(u0)n∣2≤12∥u0∥H2\int_0^T\|u(t)\|_V^2\,dt =\frac12\sum_{n=1}^\infty(1-e^{-2n^2T})|(u_0)_n|^2 \leq\frac12\|u_0\|_H^2

を満たす。特にu0=(1/n)n≥1u_0=(1/n)_{n\geq1}はH∖VH\setminus Vに属するが、その解はW(0,T)W(0,T)に属する。

証明.x↦(nxn)nx\mapsto(nx_n)_nはVVからℓ2\ell^2への線形等長全単射であるから、VVは可分 Hilbert 空間である。∥x∥H≤∥x∥V\|x\|_H\leq\|x\|_Vであり、有限支持列はVVとHHで稠密である。形式qqはM=α=1M=\alpha=1、λ=0\lambda=0について有界性と強圧性の条件を満たすので、定理 4.2により変分解が存在する。標準単位ベクトルen∈Ve_n\in Vで方程式を試験し、命題 2.3をスカラー時間関数とene_nの積に適用すると、各成分は

un′+n2un=0,un(0)=(u0)nu_n'+n^2u_n=0,\qquad u_n(0)=(u_0)_n

を弱微分の意味で満たす。§E10.48 定理 4.2と§E13.7 系 4.2の一意性からun(t)=e−n2t(u0)nu_n(t)=e^{-n^2t}(u_0)_nである。非負の級数を積分すると表示された評価を得る。∑n−2<∞\sum n^{-2}<\inftyかつ∑n2∣1/n∣2=∞\sum n^2|1/n|^2=\inftyであるから、最後の初期値はH∖VH\setminus Vに属する。▨

例 6.2.H,VH,Vを例 6.1の三つ組とし、T>0T>0、u0=(1/n)n≥1u_0=(1/n)_{n\geq1}、f=0f=0とする。このとき、零形式a(t;v,w)=0a(t;v,w)=0に対してu(0)=u0u(0)=u_0、u′=0u'=0を満たすu∈W(0,T)u\in W(0,T)は存在しない。したがって、正のα\alphaをもつ Gårding の不等式を除くと、HHの任意の初期値に対する変分解の存在は成り立たない。

証明.u∈W(0,T)u\in W(0,T)がu′=0u'=0とu(0)=u0u(0)=u_0を満たすと仮定する。§E10.48 定理 4.2により、uuのV×V^\times値連続代表は定数Iu0Iu_0である。定理 2.2のHH値連続代表も、IIによる像がこの代表に一致するので、IIの単射性からすべてのttでu(t)=u0u(t)=u_0となる。しかしu0∉Vu_0\notin Vであるから、この定数関数はL2(0,T;V)L^2(0,T;V)に属さず、u∈W(0,T)u\in W(0,T)に反する。

零形式に正のα\alphaと有限のλ\lambdaによる Gårding の不等式が成り立つとすれば、v=env=e_nを代入して0≥αn2−λ0\geq\alpha n^2-\lambdaがすべてのnnで成り立つことになる。n→∞n\to\inftyとするとこの不等式は成立しない。▨

例 6.3.T>0T>0、V=H=RV=H=\Rとし、b(t)=1+1(T/2,T)(t)b(t)=1+\mathbf1_{(T/2,T)}(t)、a(t;v,w)=b(t)vwa(t;v,w)=b(t)vwとする。このとき、初期値u(0)=1u(0)=1と外力f=0f=0に対する変分解は

u(t)=exp⁡(−∫0tb(s) ds)={e−t,0≤t≤T/2,e−2t+T/2,T/2≤t≤Tu(t)=\exp\left(-\int_0^t b(s)\,ds\right) =\begin{cases} e^{-t},&0\leq t\leq T/2,\\ e^{-2t+T/2},&T/2\leq t\leq T \end{cases}

である。uuはH1(0,T)H^1(0,T)に属するが、t=T/2t=T/2では微分可能でない。

証明.1≤b≤21\leq b\leq2であるから、aaはM=2M=2、α=1\alpha=1、λ=0\lambda=0について可測性、有界性、強圧性の条件を満たす。表示された関数はT/2T/2で連続であり、各側でu′=−buu'=-buを満たす。両側で積分した等式を接続するとu(t)=1−∫0tb(s)u(s) dsu(t)=1-\int_0^t b(s)u(s)\,dsとなる。bu∈L2(0,T)bu\in L^2(0,T)であるから、§E10.48 定理 4.2によりu∈H1(0,T)u\in H^1(0,T)かつu′=−buu'=-buである。定理 4.2により、これは唯一の変分解である。T/2T/2における左右の導関数は−e−T/2-e^{-T/2}と−2e−T/2-2e^{-T/2}であり、一致しない。▨

注意 6.4.

非線形の発展方程式では、L2(0,T;V)L^2(0,T;V)における弱収束だけから非線形項の極限を定めることができるとは限らない。空間のコンパクトな埋め込みと時間微分の評価を組み合わせて強収束を得る方法に、Aubin–Lions 型のコンパクト性がある。その定式化と証明はここでは扱わない。

前提記事

13 本の記事・単元を表示