§E13.8パラメータに関する微分可能性

最終更新

常微分方程式の初期値問題の解は、初期値や右辺に含まれるパラメータが変わるとともに変化し、その存在区間も変わることがある。たとえばy′=λ+y2y'=\lambda+y^2、y(0)=0y(0)=0の解は、λ>0\lambda>0ではλtan⁡(λ t)\sqrt\lambda\tan(\sqrt\lambda\,t)であり、存在区間は∣t∣<π/(2λ)\lvert t\rvert<\pi/(2\sqrt\lambda)に限られる。一方λ≤0\lambda\le0では存在区間はR\R全体である。したがって、解が初期値とパラメータの関数としてどのような正則性をもつかを述べるには、まず関数の定義域を明確にしなければならない。本記事では、解写像の定義域を定め、初期値とパラメータについての解の微分可能性を変分方程式によって述べる。

1 パラメータ付き Carathéodory 条件

定義 1.1.d,p∈N≥1d,p\in\NNとし、Rd\R^dとRp\R^pにノルム∥⋅∥\lVert\cdot\rVertを固定する。Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^pを開集合とし、f:Ω→Rdf:\Omega\to\R^dを写像とする。ffがパラメータ付き Carathéodory 条件 (parametric Carathéodory conditions) を局所的に満たすとは、長さが正の有限閉区間JJと半径が正の閉球B⊆RdB\subseteq\R^d、C⊆RpC\subseteq\R^pによる閉矩形J×B×CJ\times B\times CでΩ\Omegaに含まれるものそれぞれに対して、Lebesgue 零集合N⊆JN\subseteq Jと非負関数m,ℓ∈L1(J)m,\ell\in L^1(J)が存在し、次の条件が成り立つことをいう。

  1. すべての(u,λ)∈B×C(\boldsymbol u,\lambda)\in B\times Cに対して、t↦f(t,u,λ)t\mapsto f(t,\boldsymbol u,\lambda)はJJ上で Lebesgue 可測である。
  2. すべてのt∈J∖Nt\in J\setminus Nに対して、(u,λ)↦f(t,u,λ)(\boldsymbol u,\lambda)\mapsto f(t,\boldsymbol u,\lambda)はB×CB\times C上で連続である。
  3. すべてのt∈J∖Nt\in J\setminus N、u,v∈B\boldsymbol u,\boldsymbol v\in B、λ∈C\lambda\in Cに対して ∥f(t,u,λ)∥≤m(t),∥f(t,u,λ)−f(t,v,λ)∥≤ℓ(t)∥u−v∥\lVert f(t,\boldsymbol u,\lambda)\rVert\le m(t),\qquad \lVert f(t,\boldsymbol u,\lambda)-f(t,\boldsymbol v,\lambda)\rVert\le\ell(t)\lVert\boldsymbol u-\boldsymbol v\rVert が成り立つ。

パラメータ付き Carathéodory 条件は、§E13.4 定義 2.1の Carathéodory 条件と異なり、時間について可積分な Lipschitz 係数ℓ\ellによるu\boldsymbol uについての Lipschitz 条件を含む。

補題 1.2.n,k∈N≥1n,k\in\NNとし、Rn\R^nにノルム∥⋅∥\lVert\cdot\rVertを固定する。Ω⊆R×Rn\Omega\subseteq\R\times\R^nを開集合とする。このとき、次が成り立つ。

  1. 任意のコンパクト集合L⊆ΩL\subseteq\Omegaに対して、あるε>0\varepsilon>0が存在し、(t,v)∈L(t,\boldsymbol v)\in L、∣s−t∣≤ε\lvert s-t\rvert\le\varepsilon、∥w−v∥≤ε\lVert\boldsymbol w-\boldsymbol v\rVert\le\varepsilonならば(s,w)∈Ω(s,\boldsymbol w)\in\Omegaである。
  2. 写像g:Ω→Rkg:\Omega\to\R^kが次を満たすとする。長さが正の有限閉区間JJと半径が正の閉球B⊆RnB\subseteq\R^nについてJ×B⊆ΩJ\times B\subseteq\Omegaならば、各w∈B\boldsymbol w\in Bでt↦g(t,w)t\mapsto g(t,\boldsymbol w)はJJ上で Lebesgue 可測であり、ある Lebesgue 零集合N⊆JN\subseteq Jが存在し、各t∈J∖Nt\in J\setminus Nでw↦g(t,w)\boldsymbol w\mapsto g(t,\boldsymbol w)はBB上で連続である。このとき、a<ba<bと、すべてのt∈[a,b]t\in[a,b]で(t,γ(t))∈Ω(t,\boldsymbol\gamma(t))\in\Omegaを満たす連続写像γ:[a,b]→Rn\boldsymbol\gamma:[a,b]\to\R^nに対して、t↦g(t,γ(t))t\mapsto g(t,\boldsymbol\gamma(t))は[a,b][a,b]上で Lebesgue 可測である。

証明.主張 (1)を示す。R×Rn\R\times\R^nの点の距離dist∞\mathrm{dist}_\inftyと閉ρ\rho近傍LρL_\rhoを§E13.6 補題 3.1のとおりに定める。§E13.6 補題 3.1 (1)をLLに適用してρL>0\rho_L>0を取り、ε=ρL\varepsilon=\rho_Lとおく。(t,v)∈L(t,\boldsymbol v)\in L、∣s−t∣≤ε\lvert s-t\rvert\le\varepsilon、∥w−v∥≤ε\lVert\boldsymbol w-\boldsymbol v\rVert\le\varepsilonならばdist∞((s,w),L)≤dist∞((s,w),(t,v))≤ε\mathrm{dist}_\infty((s,\boldsymbol w),L)\le\mathrm{dist}_\infty((s,\boldsymbol w),(t,\boldsymbol v))\le\varepsilonであるから、(s,w)∈Lε⊆Ω(s,\boldsymbol w)\in L_\varepsilon\subseteq\Omegaである。

主張 (2)を示す。L={(t,γ(t)):a≤t≤b}L=\{(t,\boldsymbol\gamma(t)):a\le t\le b\}は連続写像によるコンパクト集合の像であるからコンパクトであり、Ω\Omegaに含まれる。主張 (1)のε\varepsilonを取る。γ\boldsymbol\gammaは一様連続であるから、∣s−t∣≤δ\lvert s-t\rvert\le\deltaならば∥γ(s)−γ(t)∥≤ε\lVert\boldsymbol\gamma(s)-\boldsymbol\gamma(t)\rVert\le\varepsilonとなるδ>0\delta>0が存在する。分点a=s0<s1<⋯<sj=ba=s_0<s_1<\cdots<s_j=bをsi−si−1≤min⁡{δ,ε}s_i-s_{i-1}\le\min\{\delta,\varepsilon\}となるように取り、Ji=[si−1,si]J_i=[s_{i-1},s_i]、Bi={w:∥w−γ(si−1)∥≤ε}B_i=\{\boldsymbol w:\lVert\boldsymbol w-\boldsymbol\gamma(s_{i-1})\rVert\le\varepsilon\}とおく。(si−1,γ(si−1))∈L(s_{i-1},\boldsymbol\gamma(s_{i-1}))\in LであるからJi×Bi⊆ΩJ_i\times B_i\subseteq\Omegaであり、γ(Ji)⊆Bi\boldsymbol\gamma(J_i)\subseteq B_iである。ggのJi×BiJ_i\times B_iへの制限は、d=nd=nとして§E13.4 補題 2.3の仮定を満たし、γ∣Ji\boldsymbol\gamma|_{J_i}はBiB_iに値をとる連続写像であるから可測である。同補題によりt↦g(t,γ(t))t\mapsto g(t,\boldsymbol\gamma(t))は各JiJ_i上で可測であり、有限個のJiJ_iの和である[a,b][a,b]上で可測である。▨

補題 1.3.d,p∈N≥1d,p\in\NNとし、Rd\R^dとRp\R^pにノルム∥⋅∥\lVert\cdot\rVertを固定する。開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dがパラメータ付き Carathéodory 条件を局所的に満たすとする。このとき、次が成り立つ。

  1. 各λ∈Rp\lambda\in\R^pについて、Ωλ={(t,u):(t,u,λ)∈Ω}\Omega^\lambda=\{(t,\boldsymbol u):(t,\boldsymbol u,\lambda)\in\Omega\}はR×Rd\R\times\R^dの開集合であり、(t,u)↦f(t,u,λ)(t,\boldsymbol u)\mapsto f(t,\boldsymbol u,\lambda)はΩλ\Omega^\lambda上で Carathéodory 条件と時間について可積分な Lipschitz 係数の条件を局所的に満たす。
  2. a<ba<bとし、γ:[a,b]→Rd×Rp\boldsymbol\gamma:[a,b]\to\R^d\times\R^pを、すべてのt∈[a,b]t\in[a,b]で(t,γ(t))∈Ω(t,\boldsymbol\gamma(t))\in\Omegaを満たす連続写像とする。このときt↦f(t,γ(t))t\mapsto f(t,\boldsymbol\gamma(t))は[a,b][a,b]上で Lebesgue 可測である。
  3. a<ba<b、λ0∈Rp\lambda_0\in\R^pとし、y:[a,b]→Rd\boldsymbol y:[a,b]\to\R^dを、すべてのt∈[a,b]t\in[a,b]で(t,y(t),λ0)∈Ω(t,\boldsymbol y(t),\lambda_0)\in\Omegaを満たす連続写像とする。このとき、ρ>0\rho>0、Lebesgue 零集合N⊆[a,b]N\subseteq[a,b]、非負関数m,ℓ∈L1([a,b])m,\ell\in L^1([a,b])が存在し、 Kρ={(t,u,λ):a≤t≤b, ∥u−y(t)∥≤ρ, ∥λ−λ0∥≤ρ}K_\rho=\{(t,\boldsymbol u,\lambda):a\le t\le b,\ \lVert\boldsymbol u-\boldsymbol y(t)\rVert\le\rho,\ \lVert\lambda-\lambda_0\rVert\le\rho\} はΩ\Omegaに含まれ、すべてのt∈[a,b]∖Nt\in[a,b]\setminus Nについて次が成り立つ。(u,λ)↦f(t,u,λ)(\boldsymbol u,\lambda)\mapsto f(t,\boldsymbol u,\lambda)は{(u,λ):(t,u,λ)∈Kρ}\{(\boldsymbol u,\lambda):(t,\boldsymbol u,\lambda)\in K_\rho\}上で連続であり、(t,u,λ),(t,v,λ)∈Kρ(t,\boldsymbol u,\lambda),(t,\boldsymbol v,\lambda)\in K_\rhoならば ∥f(t,u,λ)∥≤m(t),∥f(t,u,λ)−f(t,v,λ)∥≤ℓ(t)∥u−v∥\lVert f(t,\boldsymbol u,\lambda)\rVert\le m(t),\qquad \lVert f(t,\boldsymbol u,\lambda)-f(t,\boldsymbol v,\lambda)\rVert\le\ell(t)\lVert\boldsymbol u-\boldsymbol v\rVert である。

証明.Rd×Rp\R^d\times\R^pにノルム∥(u,λ)∥=max⁡{∥u∥,∥λ∥}\lVert(\boldsymbol u,\lambda)\rVert=\max\{\lVert\boldsymbol u\rVert,\lVert\lambda\rVert\}を入れる。このノルムの半径r>0r>0の閉球は、半径rrの閉球B⊆RdB\subseteq\R^dとC⊆RpC\subseteq\R^pの直積B×CB\times Cである。

主張 (1)を示す。Ωλ\Omega^\lambdaは連続写像(t,u)↦(t,u,λ)(t,\boldsymbol u)\mapsto(t,\boldsymbol u,\lambda)によるΩ\Omegaの逆像であるから開集合である。Ωλ\Omega^\lambdaに含まれる閉矩形J×BJ\times Bを取る。J×B×{λ}J\times B\times\{\lambda\}はコンパクトでΩ\Omegaに含まれるから、補題 1.2 (1)のε\varepsilonを取りC={λ′:∥λ′−λ∥≤ε}C=\{\lambda':\lVert\lambda'-\lambda\rVert\le\varepsilon\}とおくと、J×B×C⊆ΩJ\times B\times C\subseteq\Omegaである。この閉矩形に対する定義 1.1のN,m,ℓN,m,\ellを取り、λ\lambdaを固定すると、(t,u)↦f(t,u,λ)(t,\boldsymbol u)\mapsto f(t,\boldsymbol u,\lambda)はJ×BJ\times B上で§E13.4 定義 2.1の条件と Lipschitz 係数ℓ\ellの条件を満たす。

主張 (2)を示す。上のノルムの閉球はB×CB\times Cの形であるから、ffはn=d+pn=d+p、k=dk=dとして補題 1.2 (2)の仮定を満たす。同補題から結論を得る。

主張 (3)を示す。n=d+pn=d+pとし、Rn=Rd×Rp\R^n=\R^d\times\R^pに上のノルムを入れ、距離dist∞\mathrm{dist}_\inftyと閉ρ\rho近傍を§E13.6 補題 3.1のとおりに定める。L={(t,y(t),λ0):a≤t≤b}L=\{(t,\boldsymbol y(t),\lambda_0):a\le t\le b\}は連続写像によるコンパクト集合の像であるからコンパクトであり、Ω\Omegaに含まれる。§E13.6 補題 3.1 (1)をLLに適用してρ=ρL>0\rho=\rho_L>0を取ると、Lρ⊆ΩL_\rho\subseteq\Omegaである。(t,u,λ)∈Kρ(t,\boldsymbol u,\lambda)\in K_\rhoならばdist∞((t,u,λ),(t,y(t),λ0))=max⁡{∥u−y(t)∥,∥λ−λ0∥}≤ρ\mathrm{dist}_\infty((t,\boldsymbol u,\lambda),(t,\boldsymbol y(t),\lambda_0))=\max\{\lVert\boldsymbol u-\boldsymbol y(t)\rVert,\lVert\lambda-\lambda_0\rVert\}\le\rhoであるから、Kρ⊆Lρ⊆ΩK_\rho\subseteq L_\rho\subseteq\Omegaである。

上のノルムの半径r>0r>0の閉球はB×CB\times Cの形であるから、定義 1.1 条件 (b)と定義 1.1 条件 (c)により、ffはk=dk=dとして§E13.6 補題 3.1 (2)の仮定と§E13.6 補題 3.1 (3)の仮定を満たし、閉矩形ごとの零集合と非負関数はパラメータ付き Carathéodory 条件のN,m,ℓN,m,\ellである。§E13.6 補題 3.1 (3)をLLとρ\rhoに適用して、零集合NL⊆RN_L\subseteq\Rと非負関数mL,ℓL∈L1(R)m_L,\ell_L\in L^1(\R)を取る。N=NL∩[a,b]N=N_L\cap[a,b]、m=mL∣[a,b]m=m_L|_{[a,b]}、ℓ=ℓL∣[a,b]\ell=\ell_L|_{[a,b]}とおく。t∈[a,b]∖Nt\in[a,b]\setminus Nでは、KρK_\rhoの時刻ttにおける切り口はLρL_\rhoの時刻ttにおける切り口に含まれるから、同補題の結論により、(u,λ)↦f(t,u,λ)(\boldsymbol u,\lambda)\mapsto f(t,\boldsymbol u,\lambda)はKρK_\rhoの切り口上で連続であり、(t,u,λ),(t,v,λ)∈Kρ(t,\boldsymbol u,\lambda),(t,\boldsymbol v,\lambda)\in K_\rhoならば∥f(t,u,λ)∥≤m(t)\lVert f(t,\boldsymbol u,\lambda)\rVert\le m(t)、∥f(t,u,λ)−f(t,v,λ)∥≤ℓ(t)∥u−v∥\lVert f(t,\boldsymbol u,\lambda)-f(t,\boldsymbol v,\lambda)\rVert\le\ell(t)\lVert\boldsymbol u-\boldsymbol v\rVertである。▨

2 解写像の連続性

定義 2.1.d,p∈N≥1d,p\in\NNとし、開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dがパラメータ付き Carathéodory 条件を局所的に満たすとする。(t0,η,λ)∈Ω(t_0,\boldsymbol\eta,\lambda)\in\Omegaに対して、Ωλ={(t,u):(t,u,λ)∈Ω}\Omega^\lambda=\{(t,\boldsymbol u):(t,\boldsymbol u,\lambda)\in\Omega\}上の初期値問題

y′=f(t,y,λ),y(t0)=η\boldsymbol y'=f(t,\boldsymbol y,\lambda),\qquad \boldsymbol y(t_0)=\boldsymbol\eta

は一意な最大絶対連続解をもつ(補題 1.3 (1)、§E13.6 定理 2.1)。その最大存在区間をI(t0,η,λ)I(t_0,\boldsymbol\eta,\lambda)と書き、

W={(t,t0,η,λ)∈R×R×Rd×Rp:(t0,η,λ)∈Ω, t∈I(t0,η,λ)}W=\{(t,t_0,\boldsymbol\eta,\lambda)\in\R\times\R\times\R^d\times\R^p:(t_0,\boldsymbol\eta,\lambda)\in\Omega,\ t\in I(t_0,\boldsymbol\eta,\lambda)\}

とおく。(t,t0,η,λ)∈W(t,t_0,\boldsymbol\eta,\lambda)\in Wに、この最大絶対連続解の時刻ttにおける値φ(t,t0,η,λ)\varphi(t,t_0,\boldsymbol\eta,\lambda)を対応させる写像φ:W→Rd\varphi:W\to\R^dを、ffの解写像 (solution map) という。

定理 2.2.d,p∈N≥1d,p\in\NNとし、Rd\R^dとRp\R^pにノルム∥⋅∥\lVert\cdot\rVertを固定する。開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dがパラメータ付き Carathéodory 条件を局所的に満たすとし、φ:W→Rd\varphi:W\to\R^dをその解写像、I(t0,η,λ)I(t_0,\boldsymbol\eta,\lambda)を最大存在区間とする。このとき、次が成り立つ。

  1. (t0,η,λ)∈Ω(t_0,\boldsymbol\eta,\lambda)\in\Omegaと、a<t0<ba<t_0<bかつ[a,b]⊆I(t0,η,λ)[a,b]\subseteq I(t_0,\boldsymbol\eta,\lambda)を満たすa,ba,bを取る。このとき(t0,η,λ)(t_0,\boldsymbol\eta,\lambda)の開近傍V⊆ΩV\subseteq\Omegaが存在し、すべての(t0′,η′,λ′)∈V(t_0',\boldsymbol\eta',\lambda')\in Vについて[a,b]⊆I(t0′,η′,λ′)[a,b]\subseteq I(t_0',\boldsymbol\eta',\lambda')であり、(t0′,η′,λ′)→(t0,η,λ)(t_0',\boldsymbol\eta',\lambda')\to(t_0,\boldsymbol\eta,\lambda)のとき max⁡a≤t≤b∥φ(t,t0′,η′,λ′)−φ(t,t0,η,λ)∥⟶0\max_{a\le t\le b}\lVert\varphi(t,t_0',\boldsymbol\eta',\lambda')-\varphi(t,t_0,\boldsymbol\eta,\lambda)\rVert\longrightarrow0 である。
  2. WWはR×R×Rd×Rp\R\times\R\times\R^d\times\R^pの開集合であり、φ\varphiはWW上で連続である。

証明.主張 (1)を示す。y=φ(⋅,t0,η,λ)∣[a,b]\boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)|_{[a,b]}とおく。補題 1.3 (3)をy\boldsymbol yとλ0=λ\lambda_0=\lambdaに適用し、ρ>0\rho>0、零集合NN、非負関数m,ℓ∈L1([a,b])m,\ell\in L^1([a,b])とKρ⊆ΩK_\rho\subseteq\Omegaを取る。∥λ′−λ∥≤ρ\lVert\lambda'-\lambda\rVert\le\rhoを満たすλ′\lambda'に対して、t∈[a,b]∖Nt\in[a,b]\setminus Nでは

rλ′(t)=∥f(t,y(t),λ′)−f(t,y(t),λ)∥r_{\lambda'}(t)=\lVert f(t,\boldsymbol y(t),\lambda')-f(t,\boldsymbol y(t),\lambda)\rVert

とし、t∈Nt\in Nではrλ′(t)=0r_{\lambda'}(t)=0とおく。補題 1.3 (2)によりrλ′r_{\lambda'}は可測であり、rλ′≤2mr_{\lambda'}\le2mである。t∈[a,b]∖Nt\in[a,b]\setminus NではKρK_\rhoの切り口上の連続性から、λ′→λ\lambda'\to\lambdaのときrλ′(t)→0r_{\lambda'}(t)\to0である。λ\lambdaへ収束する任意の点列に§E10.11 定理 3.1を適用すると、λ′→λ\lambda'\to\lambdaのとき∫abrλ′→0\int_a^br_{\lambda'}\to0である。

E=exp⁡(∫abℓ)E=\exp(\int_a^b\ell)とおき、(t0′,η′,λ′)∈Ω(t_0',\boldsymbol\eta',\lambda')\in\Omegaに対して

Q(t0′,η′,λ′)=E(∥η′−y(t0′)∥+∫abrλ′(s) ds)Q(t_0',\boldsymbol\eta',\lambda')=E\left(\lVert\boldsymbol\eta'-\boldsymbol y(t_0')\rVert+\int_a^br_{\lambda'}(s)\,ds\right)

とおく(a≤t0′≤ba\le t_0'\le b、∥λ′−λ∥≤ρ\lVert\lambda'-\lambda\rVert\le\rhoの範囲で)。∥η′−y(t0′)∥≤∥η′−η∥+∥y(t0)−y(t0′)∥\lVert\boldsymbol\eta'-\boldsymbol y(t_0')\rVert\le\lVert\boldsymbol\eta'-\boldsymbol\eta\rVert+\lVert\boldsymbol y(t_0)-\boldsymbol y(t_0')\rVertとy\boldsymbol yの連続性と前段により、あるε∈(0,ρ]\varepsilon\in(0,\rho]が存在し、V={(t0′,η′,λ′):∣t0′−t0∣<ε, ∥η′−η∥<ε, ∥λ′−λ∥<ε}V=\{(t_0',\boldsymbol\eta',\lambda'):\lvert t_0'-t_0\rvert<\varepsilon,\ \lVert\boldsymbol\eta'-\boldsymbol\eta\rVert<\varepsilon,\ \lVert\lambda'-\lambda\rVert<\varepsilon\}の各点でa<t0′<ba<t_0'<bかつQ(t0′,η′,λ′)<ρQ(t_0',\boldsymbol\eta',\lambda')<\rhoとなり、(t0′,η′,λ′)→(t0,η,λ)(t_0',\boldsymbol\eta',\lambda')\to(t_0,\boldsymbol\eta,\lambda)のときQ(t0′,η′,λ′)→0Q(t_0',\boldsymbol\eta',\lambda')\to0である。(t0′,η′,λ′)∈V(t_0',\boldsymbol\eta',\lambda')\in Vならば∥η′−y(t0′)∥<ρ\lVert\boldsymbol\eta'-\boldsymbol y(t_0')\rVert<\rhoであるから、(t0′,η′,λ′)∈Kρ⊆Ω(t_0',\boldsymbol\eta',\lambda')\in K_\rho\subseteq\Omegaである。

(t0′,η′,λ′)∈V(t_0',\boldsymbol\eta',\lambda')\in Vを固定する。Ω′=Ωλ∩Ωλ′\Omega'=\Omega^\lambda\cap\Omega^{\lambda'}は開集合であり、補題 1.3 (1)によりf(⋅,⋅,λ′)f(\cdot,\cdot,\lambda')はΩ′\Omega'上で Carathéodory 条件と時間について可積分な Lipschitz 係数の条件を局所的に満たす。K′={(t,u):a≤t≤b, ∥u−y(t)∥≤ρ}K'=\{(t,\boldsymbol u):a\le t\le b,\ \lVert\boldsymbol u-\boldsymbol y(t)\rVert\le\rho\}は、Kρ⊆ΩK_\rho\subseteq\OmegaによりΩ′\Omega'に含まれるコンパクト集合である。t∈[a,b]∖Nt\in[a,b]\setminus Nと(t,u),(t,v)∈K′(t,\boldsymbol u),(t,\boldsymbol v)\in K'について∥f(t,u,λ′)−f(t,v,λ′)∥≤ℓ(t)∥u−v∥\lVert f(t,\boldsymbol u,\lambda')-f(t,\boldsymbol v,\lambda')\rVert\le\ell(t)\lVert\boldsymbol u-\boldsymbol v\rVertである。§E13.7 定理 3.1をΩ′\Omega'上の写像f(⋅,⋅,λ)f(\cdot,\cdot,\lambda)とg=f(⋅,⋅,λ′)g=f(\cdot,\cdot,\lambda')、K=K′K=K'、r=rλ′r=r_{\lambda'}、t∗=t0′t_*=t_0'、初期値η′\boldsymbol\eta'に適用すると、z′=f(t,z,λ′)\boldsymbol z'=f(t,\boldsymbol z,\lambda')、z(t0′)=η′\boldsymbol z(t_0')=\boldsymbol\eta'の[a,b][a,b]上の絶対連続解z\boldsymbol zが存在し、max⁡[a,b]∥z−y∥≤Q(t0′,η′,λ′)<ρ\max_{[a,b]}\lVert\boldsymbol z-\boldsymbol y\rVert\le Q(t_0',\boldsymbol\eta',\lambda')<\rhoである。

(α′,β′)=I(t0′,η′,λ′)(\alpha',\beta')=I(t_0',\boldsymbol\eta',\lambda')とし、ψ=φ(⋅,t0′,η′,λ′)\boldsymbol\psi=\varphi(\cdot,t_0',\boldsymbol\eta',\lambda')とする。z∣(a,b)\boldsymbol z|_{(a,b)}はt0′t_0'を含む開区間上の絶対連続解であるから、Ωλ′\Omega^{\lambda'}上で§E13.6 定理 2.1を適用すると(a,b)⊆(α′,β′)(a,b)\subseteq(\alpha',\beta')かつ(a,b)(a,b)上でψ=z\boldsymbol\psi=\boldsymbol zである。β′=b\beta'=bならば、[t0′,b)[t_0',b)上のψ\boldsymbol\psiのグラフはコンパクト集合K′⊆Ωλ′K'\subseteq\Omega^{\lambda'}に含まれる。これは§E13.6 系 4.3に反する。従ってβ′>b\beta'>bであり、同様にα′<a\alpha'<aである。[a,b]⊆I(t0′,η′,λ′)[a,b]\subseteq I(t_0',\boldsymbol\eta',\lambda')であり、連続性により[a,b][a,b]上でψ=z\boldsymbol\psi=\boldsymbol zであるから、

max⁡a≤t≤b∥φ(t,t0′,η′,λ′)−φ(t,t0,η,λ)∥≤Q(t0′,η′,λ′)\max_{a\le t\le b}\lVert\varphi(t,t_0',\boldsymbol\eta',\lambda')-\varphi(t,t_0,\boldsymbol\eta,\lambda)\rVert\le Q(t_0',\boldsymbol\eta',\lambda')

であり、右辺は(t0′,η′,λ′)→(t0,η,λ)(t_0',\boldsymbol\eta',\lambda')\to(t_0,\boldsymbol\eta,\lambda)のとき零へ収束する。

主張 (2)を示す。(t,t0,η,λ)∈W(t,t_0,\boldsymbol\eta,\lambda)\in Wを取る。I(t0,η,λ)I(t_0,\boldsymbol\eta,\lambda)はt,t0t,t_0を含む開区間であるから、a<min⁡{t,t0}a<\min\{t,t_0\}、max⁡{t,t0}<b\max\{t,t_0\}<b、[a,b]⊆I(t0,η,λ)[a,b]\subseteq I(t_0,\boldsymbol\eta,\lambda)を満たすa,ba,bが存在する。主張 (1)のVVを取ると、開集合(a,b)×V(a,b)\times Vは(t,t0,η,λ)(t,t_0,\boldsymbol\eta,\lambda)を含みWWに含まれる。(t′,t0′,η′,λ′)∈(a,b)×V(t',t_0',\boldsymbol\eta',\lambda')\in(a,b)\times Vについて

∥φ(t′,t0′,η′,λ′)−φ(t,t0,η,λ)∥≤max⁡a≤s≤b∥φ(s,t0′,η′,λ′)−φ(s,t0,η,λ)∥+∥φ(t′,t0,η,λ)−φ(t,t0,η,λ)∥\lVert\varphi(t',t_0',\boldsymbol\eta',\lambda')-\varphi(t,t_0,\boldsymbol\eta,\lambda)\rVert \le\max_{a\le s\le b}\lVert\varphi(s,t_0',\boldsymbol\eta',\lambda')-\varphi(s,t_0,\boldsymbol\eta,\lambda)\rVert+\lVert\varphi(t',t_0,\boldsymbol\eta,\lambda)-\varphi(t,t_0,\boldsymbol\eta,\lambda)\rVert

であり、右辺の第一項は主張 (1)により、第二項はφ(⋅,t0,η,λ)\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)の連続性により、(t′,t0′,η′,λ′)→(t,t0,η,λ)(t',t_0',\boldsymbol\eta',\lambda')\to(t,t_0,\boldsymbol\eta,\lambda)のとき零へ収束する。▨

3 変分方程式と一階の微分可能性

補題 3.1.K∈{R,C}\mathbb K\in\{\R,\C\}、k∈N≥1k\in\NNとし、Kk\mathbb K^kにノルム∥⋅∥\lVert\cdot\rVertを固定し、k×kk\times k行列にはその作用素ノルムを用いる。

  1. a<ba<b、J=[a,b]J=[a,b]、A∈L1(J;Kk×k)A\in L^1(J;\mathbb K^{k\times k})、b∈L1(J;Kk)\boldsymbol b\in L^1(J;\mathbb K^k)、t∗∈Jt_*\in J、ξ∈Kk\boldsymbol\xi\in\mathbb K^kとする。このとき、JJのほとんど至る所でx′=A(t)x+b(t)\boldsymbol x'=A(t)\boldsymbol x+\boldsymbol b(t)を満たし、x(t∗)=ξ\boldsymbol x(t_*)=\boldsymbol\xiを満たす絶対連続関数x:J→Kk\boldsymbol x:J\to\mathbb K^kが一意に存在する。
  2. 主張 (1)のJ,A,b,t∗,ξ,xJ,A,\boldsymbol b,t_*,\boldsymbol\xi,\boldsymbol xと、An∈L1(J;Kk×k)A_n\in L^1(J;\mathbb K^{k\times k})、bn∈L1(J;Kk)\boldsymbol b_n\in L^1(J;\mathbb K^k)、tn∈Jt_n\in J、ξn∈Kk\boldsymbol\xi_n\in\mathbb K^k(n∈N≥1n\in\NN)が ∫J∥An−A∥⟶0,∫J∥bn−b∥⟶0,tn⟶t∗,ξn⟶ξ\int_J\lVert A_n-A\rVert\longrightarrow0,\qquad\int_J\lVert\boldsymbol b_n-\boldsymbol b\rVert\longrightarrow0,\qquad t_n\longrightarrow t_*,\qquad\boldsymbol\xi_n\longrightarrow\boldsymbol\xi を満たすとし、xn\boldsymbol x_nをAn,bn,tn,ξnA_n,\boldsymbol b_n,t_n,\boldsymbol\xi_nに対する主張 (1)の解とする。このときmax⁡t∈J∥xn(t)−x(t)∥→0\max_{t\in J}\lVert\boldsymbol x_n(t)-\boldsymbol x(t)\rVert\to0である。
  3. I⊆RI\subseteq\Rを内部が空でない区間、A∈Lloc1(I;Kk×k)A\in L^1_{\mathrm{loc}}(I;\mathbb K^{k\times k})、b∈Lloc1(I;Kk)\boldsymbol b\in L^1_{\mathrm{loc}}(I;\mathbb K^k)、t∗∈It_*\in I、ξ∈Kk\boldsymbol\xi\in\mathbb K^kとする。このとき、IIのほとんど至る所でx′=A(t)x+b(t)\boldsymbol x'=A(t)\boldsymbol x+\boldsymbol b(t)を満たし、x(t∗)=ξ\boldsymbol x(t_*)=\boldsymbol\xiを満たすx∈ACloc(I;Kk)\boldsymbol x\in AC_{\mathrm{loc}}(I;\mathbb K^k)が一意に存在する。

q∈N≥1q\in\NNとし、Kk×q\mathbb K^{k\times q}に値をとる未知関数XXと、L1(J;Kk×q)L^1(J;\mathbb K^{k\times q})またはLloc1(I;Kk×q)L^1_{\mathrm{loc}}(I;\mathbb K^{k\times q})に属するBBに対する方程式X′=A(t)X+B(t)X'=A(t)X+B(t)についても、主張 (1)、主張 (2)、主張 (3)の結論が成り立つ。

証明.主張 (3)を示す。§E13.7 系 4.2をII、AA、b\boldsymbol b、t∗t_*、η=ξ\boldsymbol\eta=\boldsymbol\xiに適用すると、結論を得る。

主張 (1)を示す。L1(J)L^1(J)の関数はJJ上で局所可積分であるから、主張 (3)をI=JI=Jに適用すると、JJのほとんど至る所で方程式を満たしx(t∗)=ξ\boldsymbol x(t_*)=\boldsymbol\xiを満たすx∈ACloc(J;Kk)\boldsymbol x\in AC_{\mathrm{loc}}(J;\mathbb K^k)が一意に存在する。JJは有限閉区間であるから、ACloc(J;Kk)AC_{\mathrm{loc}}(J;\mathbb K^k)はJJ上の絶対連続関数の全体である。

主張 (2)を示す。K=C\mathbb K=\Cのときは、Ck\C^kを実部と虚部によりR2k\R^{2k}と同一視する。固定したノルムはR2k\R^{2k}のノルムであり、各ttについてu↦A(t)u\boldsymbol u\mapsto A(t)\boldsymbol uは実線形であり、絶対連続性は実座標成分ごとに判定される(§E13.2 補題 1.2 (3))。n∈N≥1n\in\NNを固定し、§E13.7 定理 2.1をS=J×KkS=J\times\mathbb K^k、二つの写像(t,u)↦A(t)u+b(t)(t,\boldsymbol u)\mapsto A(t)\boldsymbol u+\boldsymbol b(t)と(t,u)↦An(t)u+bn(t)(t,\boldsymbol u)\mapsto A_n(t)\boldsymbol u+\boldsymbol b_n(t)、二つの解x,xn\boldsymbol x,\boldsymbol x_n、時刻tnt_n、

ℓn=∥An∥,rn=∥An−A∥max⁡J∥x∥+∥bn−b∥\ell_n=\lVert A_n\rVert,\qquad r_n=\lVert A_n-A\rVert\max_J\lVert\boldsymbol x\rVert+\lVert\boldsymbol b_n-\boldsymbol b\rVert

に適用すると、

max⁡t∈J∥xn(t)−x(t)∥≤exp⁡(∫J∥An∥)(∥ξn−x(tn)∥+∫Jrn)\max_{t\in J}\lVert\boldsymbol x_n(t)-\boldsymbol x(t)\rVert\le\exp\left(\int_J\lVert A_n\rVert\right)\left(\lVert\boldsymbol\xi_n-\boldsymbol x(t_n)\rVert+\int_Jr_n\right)

を得る。∫J∥An∥≤∫J∥A∥+∫J∥An−A∥\int_J\lVert A_n\rVert\le\int_J\lVert A\rVert+\int_J\lVert A_n-A\rVertは有界であり、∥ξn−x(tn)∥≤∥ξn−ξ∥+∥x(t∗)−x(tn)∥→0\lVert\boldsymbol\xi_n-\boldsymbol x(t_n)\rVert\le\lVert\boldsymbol\xi_n-\boldsymbol\xi\rVert+\lVert\boldsymbol x(t_*)-\boldsymbol x(t_n)\rVert\to0、∫Jrn→0\int_Jr_n\to0であるから、右辺は零へ収束する。

行列値の方程式X′=A(t)X+B(t)X'=A(t)X+B(t)は、XXとBBの各列についての方程式xj′=A(t)xj+B(t)ej\boldsymbol x_j'=A(t)\boldsymbol x_j+B(t)\boldsymbol e_jの組であるから、各列に以上の結論を適用すればよい。▨

定義 3.2.d,p,r∈N≥1d,p,r\in\NNとし、Rd\R^dとRp\R^pにノルム∥⋅∥\lVert\cdot\rVertを固定する。Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^pを開集合とし、f:Ω→Rdf:\Omega\to\R^dを写像とする。w=(u,λ)∈Rd×Rp\boldsymbol w=(\boldsymbol u,\lambda)\in\R^d\times\R^pの座標をw1,…,wd+pw_1,\ldots,w_{d+p}と書き、wiw_iについての偏微分を∂i\partial_iと書く。t∈Rt\in\Rに対してΩt={w:(t,w)∈Ω}\Omega_t=\{\boldsymbol w:(t,\boldsymbol w)\in\Omega\}とおく。ffが状態とパラメータについてCrC^r級の Carathéodory 条件 (Carathéodory conditions of class C^r in the state and the parameter) を満たすとは、次の条件が成り立つことをいう。

  1. 各w∈Rd×Rp\boldsymbol w\in\R^d\times\R^pについて、t↦f(t,w)t\mapsto f(t,\boldsymbol w)は開集合{t:(t,w)∈Ω}\{t:(t,\boldsymbol w)\in\Omega\}上で Lebesgue 可測である。
  2. ある Lebesgue 零集合N⊆RN\subseteq\Rが存在し、すべてのt∈R∖Nt\in\R\setminus Nについてw↦f(t,w)\boldsymbol w\mapsto f(t,\boldsymbol w)はΩt\Omega_t上でCrC^r級である。
  3. 各コンパクト集合K⊆ΩK\subseteq\Omegaに対して非負関数mK∈L1(R)m_K\in L^1(\R)が存在し、(t,w)∈K(t,\boldsymbol w)\in Kかつt∉Nt\notin Nならば、0≤j≤r0\le j\le rとすべての(i1,…,ij)∈{1,…,d+p}j(i_1,\ldots,i_j)\in\{1,\ldots,d+p\}^jについて ∥∂ij⋯∂i1f(t,w)∥≤mK(t)\lVert\partial_{i_j}\cdots\partial_{i_1}f(t,\boldsymbol w)\rVert\le m_K(t) が成り立つ。j=0j=0の項はf(t,w)f(t,\boldsymbol w)とする。

t∈Nt\in Nでは、ffの偏導関数の値を00と定める。

補題 3.3.d,p,r∈N≥1d,p,r\in\NNとし、開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dが、零集合NNと支配関数mKm_Kをもって状態とパラメータについてCrC^r級の Carathéodory 条件を満たすとする。このとき、次が成り立つ。

  1. ffはパラメータ付き Carathéodory 条件を局所的に満たす。さらに、Ω\Omegaに含まれる各閉矩形R=J×B×CR=J\times B\times Cについて、定義 1.1 条件 (c)の関数としてm=mR∣Jm=m_R|_J、ℓ=c mR∣J\ell=c\,m_R|_Jを取ることができる。ここでccはddとRd\R^dのノルムだけで決まる定数である。
  2. 1≤j≤r1\le j\le r、(i1,…,ij)∈{1,…,d+p}j(i_1,\ldots,i_j)\in\{1,\ldots,d+p\}^j、w∈Rd×Rp\boldsymbol w\in\R^d\times\R^pについて、t↦∂ij⋯∂i1f(t,w)t\mapsto\partial_{i_j}\cdots\partial_{i_1}f(t,\boldsymbol w)は開集合{t:(t,w)∈Ω}\{t:(t,\boldsymbol w)\in\Omega\}上で Lebesgue 可測である。
  3. ffと各∂if\partial_ifは、n=d+pn=d+pとして補題 1.2 (2)のggの仮定を満たす。r≥2r\ge2ならば、各iiについて∂if\partial_ifは状態とパラメータについてCr−1C^{r-1}級の Carathéodory 条件を満たす。
  4. ffと、1≤j≤r1\le j\le rのすべての順序付き偏導関数∂ij⋯∂i1f\partial_{i_j}\cdots\partial_{i_1}fがΩ\Omega上で存在して連続ならば、ffはN=∅N=\varnothingとして状態とパラメータについてCrC^r級の Carathéodory 条件を満たす。

証明.主張 (1)を示す。Ω\Omegaに含まれる閉矩形R=J×B×CR=J\times B\times Cを取る。RRはコンパクトである。各(u,λ)∈B×C(\boldsymbol u,\lambda)\in B\times CについてJ⊆{t:(t,u,λ)∈Ω}J\subseteq\{t:(t,\boldsymbol u,\lambda)\in\Omega\}であるから、定義 3.2 条件 (a)によりt↦f(t,u,λ)t\mapsto f(t,\boldsymbol u,\lambda)はJJ上で可測である。t∈J∖Nt\in J\setminus Nではf(t,⋅)f(t,\cdot)はCrC^r級であり、特にB×CB\times C上で連続である。定義 3.2 条件 (c)をK=RK=Rに適用すると、t∈J∖Nt\in J\setminus Nと(u,λ)∈B×C(\boldsymbol u,\lambda)\in B\times Cについて∥f(t,u,λ)∥≤mR(t)\lVert f(t,\boldsymbol u,\lambda)\rVert\le m_R(t)である。Euclid ノルムを∣⋅∣\lvert\cdot\rvertと書き、すべてのx∈Rd\boldsymbol x\in\R^dで∣x∣≤c1∥x∥\lvert\boldsymbol x\rvert\le c_1\lVert\boldsymbol x\rVertかつ∥x∥≤c1∣x∣\lVert\boldsymbol x\rVert\le c_1\lvert\boldsymbol x\rvertを満たすc1≥1c_1\ge1を取る。t∈J∖Nt\in J\setminus N、λ∈C\lambda\in C、u,v∈B\boldsymbol u,\boldsymbol v\in Bとする。開集合U={u′:(t,u′,λ)∈Ω}U=\{\boldsymbol u':(t,\boldsymbol u',\lambda)\in\Omega\}は凸集合BBを含み、Φ(u′)=f(t,u′,λ)\Phi(\boldsymbol u')=f(t,\boldsymbol u',\lambda)はUU上でC1C^1級である。u′∈B\boldsymbol u'\in Bでは、d×dd\times d行列MMの Euclid ノルムに関する作用素ノルムが∑i∣Mei∣\sum_i\lvert M\boldsymbol e_i\rvert以下であることから

∥DΦ(u′)∥op≤∑i=1d∣∂if(t,u′,λ)∣≤c1d mR(t)\lVert D\Phi(\boldsymbol u')\rVert_{\mathrm{op}}\le\sum_{i=1}^d\lvert\partial_if(t,\boldsymbol u',\lambda)\rvert\le c_1d\,m_R(t)

である。§E4.3 定理 2.1により∣Φ(u)−Φ(v)∣≤c1d mR(t)∣u−v∣\lvert\Phi(\boldsymbol u)-\Phi(\boldsymbol v)\rvert\le c_1d\,m_R(t)\lvert\boldsymbol u-\boldsymbol v\rvertであり、c=c13dc=c_1^3dとおくと∥f(t,u,λ)−f(t,v,λ)∥≤c mR(t)∥u−v∥\lVert f(t,\boldsymbol u,\lambda)-f(t,\boldsymbol v,\lambda)\rVert\le c\,m_R(t)\lVert\boldsymbol u-\boldsymbol v\rVertを得る。

主張 (2)を示す。1≤j≤r1\le j\le rとし、g=∂ij−1⋯∂i1fg=\partial_{i_{j-1}}\cdots\partial_{i_1}f(j=1j=1ではg=fg=f)について、各w′\boldsymbol w'でt↦g(t,w′)t\mapsto g(t,\boldsymbol w')がOw′={t:(t,w′)∈Ω}O_{\boldsymbol w'}=\{t:(t,\boldsymbol w')\in\Omega\}上で可測であると仮定する。j=1j=1ではこの仮定は定義 3.2 条件 (a)である。w\boldsymbol wを固定し、e\boldsymbol eをwijw_{i_j}方向の単位ベクトルとする。n∈N≥1n\in\NNに対してOwO_{\boldsymbol w}上の関数hnh_nを、t∈(Ow∩Ow+e/n)∖Nt\in(O_{\boldsymbol w}\cap O_{\boldsymbol w+\boldsymbol e/n})\setminus Nでは

hn(t)=n(g(t,w+e/n)−g(t,w))h_n(t)=n\bigl(g(t,\boldsymbol w+\boldsymbol e/n)-g(t,\boldsymbol w)\bigr)

とし、その他のttではhn(t)=0h_n(t)=0と定める。Ow+e/nO_{\boldsymbol w+\boldsymbol e/n}は開集合であるからhnh_nは可測である。t∈Ow∖Nt\in O_{\boldsymbol w}\setminus Nならば、Ω\Omegaが開集合であることから十分大きいnnでt∈Ow+e/nt\in O_{\boldsymbol w+\boldsymbol e/n}であり、g(t,⋅)g(t,\cdot)はΩt\Omega_t上でCr−j+1C^{r-j+1}級であるからhn(t)→∂ijg(t,w)h_n(t)\to\partial_{i_j}g(t,\boldsymbol w)である。t∈Nt\in Nではhn(t)=0h_n(t)=0である。従ってhnh_nはOwO_{\boldsymbol w}の全点で∂ijg(⋅,w)\partial_{i_j}g(\cdot,\boldsymbol w)へ収束し、§E10.7 定理 5.1を各成分に適用すると、極限は可測である。j=1,2,…,rj=1,2,\ldots,rの順にこの議論を適用して結論を得る。

主張 (3)を示す。Ω\Omegaに含まれるJ×BJ\times B(B⊆Rd×RpB\subseteq\R^d\times\R^pは閉球)について、ffと∂if\partial_ifの時間についての可測性は定義 3.2 条件 (a)と主張 (2)をJJへ制限して得られ、t∈J∖Nt\in J\setminus NでのBB上の連続性は定義 3.2 条件 (b)による。r≥2r\ge2とする。主張 (2)により∂if\partial_ifは定義 3.2 条件 (a)を満たす。t∉Nt\notin Nでは∂if(t,⋅)\partial_if(t,\cdot)はΩt\Omega_t上でCr−1C^{r-1}級である。∂if\partial_ifのr−1r-1階以下の順序付き偏導関数はffのrr階以下の順序付き偏導関数であるから、同じmKm_Kで支配される。

主張 (4)を示す。各w\boldsymbol wについて、t↦f(t,w)t\mapsto f(t,\boldsymbol w)は開集合{t:(t,w)∈Ω}\{t:(t,\boldsymbol w)\in\Omega\}上で連続であるから可測である。各ttについてf(t,⋅)f(t,\cdot)のrr階以下の順序付き偏導関数はΩt\Omega_t上で存在して連続であるから、§E4.4 定理 1.7によりf(t,⋅)f(t,\cdot)はCrC^r級である。コンパクト集合K⊆ΩK\subseteq\Omegaに対して、KK上の∥∂ij⋯∂i1f∥\lVert\partial_{i_j}\cdots\partial_{i_1}f\rVert(0≤j≤r0\le j\le r)の最大値をMKM_Kとし、KKの時間成分への射影π(K)\pi(K)を用いてmK=MK1π(K)m_K=M_K\mathbf 1_{\pi(K)}とおく。π(K)\pi(K)はコンパクトであるからmK∈L1(R)m_K\in L^1(\R)である。▨

定義 3.4.d,p∈N≥1d,p\in\NNとし、開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dが状態とパラメータについてC1C^1級の Carathéodory 条件を満たすとする。φ\varphiをffの解写像とし、(t0,η,λ)∈Ω(t_0,\boldsymbol\eta,\lambda)\in\Omega、I=I(t0,η,λ)I=I(t_0,\boldsymbol\eta,\lambda)、y=φ(⋅,t0,η,λ)\boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)とする。DufD_{\boldsymbol u}f、DλfD_\lambda fを、u\boldsymbol uの座標とλ\lambdaの座標についてのffの偏導関数を列とする行列とし、t∈It\in Iに対して

A(t)=Duf(t,y(t),λ)∈Rd×d,B(t)=Dλf(t,y(t),λ)∈Rd×pA(t)=D_{\boldsymbol u}f(t,\boldsymbol y(t),\lambda)\in\R^{d\times d},\qquad B(t)=D_\lambda f(t,\boldsymbol y(t),\lambda)\in\R^{d\times p}

とおく。Rd×d\R^{d\times d}に値をとるYYとRd×p\R^{d\times p}に値をとるZZについての線形方程式

Y′=A(t)Y,Z′=A(t)Z+B(t)Y'=A(t)Y,\qquad Z'=A(t)Z+B(t)

を、y\boldsymbol yに沿う変分方程式 (variational equation) といい、A,BA,Bをその係数という。変分方程式の解とは、II上で局所絶対連続であり、IIのほとんど至る所でこれらの方程式を満たすY,ZY,Zをいう。

補題 3.5.d,p∈N≥1d,p\in\NNとし、開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dが状態とパラメータについてC1C^1級の Carathéodory 条件を満たすとし、φ\varphiをその解写像とする。(t0,η,λ)∈Ω(t_0,\boldsymbol\eta,\lambda)\in\Omega、I=I(t0,η,λ)I=I(t_0,\boldsymbol\eta,\lambda)とし、A,BA,Bをy=φ(⋅,t0,η,λ)\boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)に沿う変分方程式の係数とする。このとき、A∈Lloc1(I;Rd×d)A\in L^1_{\mathrm{loc}}(I;\R^{d\times d})、B∈Lloc1(I;Rd×p)B\in L^1_{\mathrm{loc}}(I;\R^{d\times p})であり、任意のt∗∈It_*\in I、Y∗∈Rd×dY_*\in\R^{d\times d}、Z∗∈Rd×pZ_*\in\R^{d\times p}に対して、Y(t∗)=Y∗Y(t_*)=Y_*、Z(t∗)=Z∗Z(t_*)=Z_*を満たすy\boldsymbol yに沿う変分方程式の解Y,ZY,Zが一意に存在する。

証明.NNを定義 3.2 条件 (b)の零集合とする。a<ba<b、[a,b]⊆I[a,b]\subseteq Iを取る。補題 3.3 (3)によりffの各偏導関数は補題 1.2 (2)の仮定を満たし、t↦(y(t),λ)t\mapsto(\boldsymbol y(t),\lambda)は連続で(t,y(t),λ)∈Ω(t,\boldsymbol y(t),\lambda)\in\Omegaであるから、A,BA,Bの各列は[a,b][a,b]上で可測である。K={(t,y(t),λ):a≤t≤b}K=\{(t,\boldsymbol y(t),\lambda):a\le t\le b\}はコンパクトでΩ\Omegaに含まれるから、定義 3.2 条件 (c)によりt∈[a,b]∖Nt\in[a,b]\setminus Nで各列のノルムはmK(t)m_K(t)以下である。従ってA,BA,Bは[a,b][a,b]上で可積分である。[a,b][a,b]は任意であるからA∈Lloc1(I;Rd×d)A\in L^1_{\mathrm{loc}}(I;\R^{d\times d})、B∈Lloc1(I;Rd×p)B\in L^1_{\mathrm{loc}}(I;\R^{d\times p})である。解の一意存在は補題 3.1 (3)を行列値の方程式に適用して得る。▨

補題 3.6.d,p∈N≥1d,p\in\NNとし、Rd\R^dとRp\R^pにノルム∥⋅∥\lVert\cdot\rVertを固定し、Rd×Rp\R^d\times\R^pにノルム∥(u,μ)∥=max⁡{∥u∥,∥μ∥}\lVert(\boldsymbol u,\mu)\rVert=\max\{\lVert\boldsymbol u\rVert,\lVert\mu\rVert\}を入れ、d×dd\times d行列とd×pd\times p行列には作用素ノルムを用いる。開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dが、零集合NNをもって状態とパラメータについてC1C^1級の Carathéodory 条件を満たすとする。a<ba<b、λ0∈Rp\lambda_0\in\R^pとし、y:[a,b]→Rd\boldsymbol y:[a,b]\to\R^dを、すべてのs∈[a,b]s\in[a,b]で(s,y(s),λ0)∈Ω(s,\boldsymbol y(s),\lambda_0)\in\Omegaを満たす連続写像とする。このとき、ρ>0\rho>0と非負関数m∈L1([a,b])m\in L^1([a,b])が存在し、

Kρ={(s,u,μ):a≤s≤b, ∥u−y(s)∥≤ρ, ∥μ−λ0∥≤ρ}K_\rho=\{(s,\boldsymbol u,\mu):a\le s\le b,\ \lVert\boldsymbol u-\boldsymbol y(s)\rVert\le\rho,\ \lVert\mu-\lambda_0\rVert\le\rho\}

について次が成り立つ。

  1. Kρ⊆ΩK_\rho\subseteq\Omegaであり、s∈[a,b]∖Ns\in[a,b]\setminus Nかつ(s,u,μ)∈Kρ(s,\boldsymbol u,\mu)\in K_\rhoならば、∥Duf(s,u,μ)∥≤m(s)\lVert D_{\boldsymbol u}f(s,\boldsymbol u,\mu)\rVert\le m(s)かつ∥Dλf(s,u,μ)∥≤m(s)\lVert D_\lambda f(s,\boldsymbol u,\mu)\rVert\le m(s)である。
  2. 写像γ:[a,b]×[0,1]→Rd×Rp\boldsymbol\gamma:[a,b]\times[0,1]\to\R^d\times\R^pは、各θ∈[0,1]\theta\in[0,1]についてs↦γ(s,θ)s\mapsto\boldsymbol\gamma(s,\theta)が連続であり、各s∈[a,b]s\in[a,b]についてθ↦γ(s,θ)\theta\mapsto\boldsymbol\gamma(s,\theta)が連続であり、すべての(s,θ)(s,\theta)で(s,γ(s,θ))∈Kρ(s,\boldsymbol\gamma(s,\theta))\in K_\rhoを満たすとする。s∈[a,b]∖Ns\in[a,b]\setminus Nでは Aγ(s)=∫01Duf(s,γ(s,θ)) dθ,Bγ(s)=∫01Dλf(s,γ(s,θ)) dθA_{\boldsymbol\gamma}(s)=\int_0^1D_{\boldsymbol u}f(s,\boldsymbol\gamma(s,\theta))\,d\theta,\qquad B_{\boldsymbol\gamma}(s)=\int_0^1D_\lambda f(s,\boldsymbol\gamma(s,\theta))\,d\theta とし、s∈[a,b]∩Ns\in[a,b]\cap NではAγ(s)=0A_{\boldsymbol\gamma}(s)=0、Bγ(s)=0B_{\boldsymbol\gamma}(s)=0とする。このとき、AγA_{\boldsymbol\gamma}とBγB_{\boldsymbol\gamma}は[a,b][a,b]上で可測であり、s∈[a,b]∖Ns\in[a,b]\setminus Nでは∥Aγ(s)∥≤m(s)\lVert A_{\boldsymbol\gamma}(s)\rVert\le m(s)かつ∥Bγ(s)∥≤m(s)\lVert B_{\boldsymbol\gamma}(s)\rVert\le m(s)である。特に、AγA_{\boldsymbol\gamma}とBγB_{\boldsymbol\gamma}は[a,b][a,b]上で可積分である。
  3. γ0(s,θ)=(y(s),λ0)\boldsymbol\gamma_0(s,\theta)=(\boldsymbol y(s),\lambda_0)とおく。γn\boldsymbol\gamma_n(n∈N≥1n\in\NN)は主張 (2)の仮定を満たし、すべてのs∈[a,b]s\in[a,b]についてn→∞n\to\inftyのときsup⁡0≤θ≤1∥γn(s,θ)−(y(s),λ0)∥→0\sup_{0\le\theta\le1}\lVert\boldsymbol\gamma_n(s,\theta)-(\boldsymbol y(s),\lambda_0)\rVert\to0であるとする。このとき、n→∞n\to\inftyのとき∫ab∥Aγn−Aγ0∥→0\int_a^b\lVert A_{\boldsymbol\gamma_n}-A_{\boldsymbol\gamma_0}\rVert\to0かつ∫ab∥Bγn−Bγ0∥→0\int_a^b\lVert B_{\boldsymbol\gamma_n}-B_{\boldsymbol\gamma_0}\rVert\to0である。

証明.主張 (1)を示す。補題 3.3 (1)によりffはパラメータ付き Carathéodory 条件を局所的に満たすから、補題 1.3 (3)をy\boldsymbol yとλ0\lambda_0に適用すると、Kρ⊆ΩK_\rho\subseteq\Omegaを満たすρ>0\rho>0が存在する。KρK_\rhoは有界な閉集合であるからコンパクトであり、定義 3.2 条件 (c)の支配関数mKρm_{K_\rho}が存在する。ノルムの同値性により、d×dd\times d行列とd×pd\times p行列の作用素ノルムがその列のノルムの最大値のcc倍以下となる定数c>0c>0が存在する。DufD_{\boldsymbol u}fとDλfD_\lambda fの列はffの偏導関数であるから、m=c mKρ∣[a,b]m=c\,m_{K_\rho}|_{[a,b]}とおくと結論を得る。

主張 (2)を示す。GGをDufD_{\boldsymbol u}fまたはDλfD_\lambda fとし、s∈[a,b]∖Ns\in[a,b]\setminus NではGγ(s)=∫01G(s,γ(s,θ)) dθG_{\boldsymbol\gamma}(s)=\int_0^1G(s,\boldsymbol\gamma(s,\theta))\,d\theta、s∈[a,b]∩Ns\in[a,b]\cap NではGγ(s)=0G_{\boldsymbol\gamma}(s)=0とおく。k∈N≥1k\in\NNとs∈[a,b]s\in[a,b]に対して

Rk(s)=1k∑i=1kG(s,γ(s,i/k))R_k(s)=\frac1k\sum_{i=1}^kG\bigl(s,\boldsymbol\gamma(s,i/k)\bigr)

とおく。s∈[a,b]∖Ns\in[a,b]\setminus Nでは、G(s,⋅)G(s,\cdot)はΩs\Omega_s上で連続であり、θ↦γ(s,θ)\theta\mapsto\boldsymbol\gamma(s,\theta)はΩs\Omega_sに値をとる連続写像であるから、被積分関数はθ\thetaについて連続であり、Gγ(s)G_{\boldsymbol\gamma}(s)は Riemann 和Rk(s)R_k(s)のk→∞k\to\inftyでの極限である。主張 (1)により∥Rk(s)∥≤m(s)\lVert R_k(s)\rVert\le m(s)であるから、∥Gγ(s)∥≤m(s)\lVert G_{\boldsymbol\gamma}(s)\rVert\le m(s)である。s∈[a,b]∩Ns\in[a,b]\cap Nでは、偏導関数の値の規約によりRk(s)=0=Gγ(s)R_k(s)=0=G_{\boldsymbol\gamma}(s)である。各iiについてs↦γ(s,i/k)s\mapsto\boldsymbol\gamma(s,i/k)は連続であるから、補題 3.3 (3)と補題 1.2 (2)によりs↦G(s,γ(s,i/k))s\mapsto G(s,\boldsymbol\gamma(s,i/k))の各列は[a,b][a,b]上で可測であり、RkR_kは可測である。RkR_kは[a,b][a,b]の各点でGγG_{\boldsymbol\gamma}へ収束するから、§E10.7 定理 5.1を各成分に適用すると、GγG_{\boldsymbol\gamma}は可測である。GγG_{\boldsymbol\gamma}は可測でありm∈L1([a,b])m\in L^1([a,b])で支配されるから、[a,b][a,b]上で可積分である。

主張 (3)を示す。GGを前段のとおりとする。γ0\boldsymbol\gamma_0は主張 (2)の仮定を満たし、s∈[a,b]∖Ns\in[a,b]\setminus NではGγ0(s)=G(s,y(s),λ0)G_{\boldsymbol\gamma_0}(s)=G(s,\boldsymbol y(s),\lambda_0)である。s∈[a,b]∖Ns\in[a,b]\setminus Nとε>0\varepsilon>0を取る。G(s,⋅)G(s,\cdot)は(y(s),λ0)(\boldsymbol y(s),\lambda_0)で連続であるから、w∈Ωs\boldsymbol w\in\Omega_sかつ∥w−(y(s),λ0)∥≤δ\lVert\boldsymbol w-(\boldsymbol y(s),\lambda_0)\rVert\le\deltaならば∥G(s,w)−G(s,y(s),λ0)∥≤ε\lVert G(s,\boldsymbol w)-G(s,\boldsymbol y(s),\lambda_0)\rVert\le\varepsilonとなるδ>0\delta>0が存在する。十分大きいnnではsup⁡θ∥γn(s,θ)−(y(s),λ0)∥≤δ\sup_\theta\lVert\boldsymbol\gamma_n(s,\theta)-(\boldsymbol y(s),\lambda_0)\rVert\le\deltaであり、Gγn(s)−Gγ0(s)G_{\boldsymbol\gamma_n}(s)-G_{\boldsymbol\gamma_0}(s)はθ↦G(s,γn(s,θ))−G(s,y(s),λ0)\theta\mapsto G(s,\boldsymbol\gamma_n(s,\theta))-G(s,\boldsymbol y(s),\lambda_0)の Riemann 和の極限で、各 Riemann 和のノルムはε\varepsilon以下であるから、∥Gγn(s)−Gγ0(s)∥≤ε\lVert G_{\boldsymbol\gamma_n}(s)-G_{\boldsymbol\gamma_0}(s)\rVert\le\varepsilonである。従って、[a,b]∖N[a,b]\setminus Nの各点で∥Gγn−Gγ0∥→0\lVert G_{\boldsymbol\gamma_n}-G_{\boldsymbol\gamma_0}\rVert\to0である。主張 (2)によりs↦∥Gγn(s)−Gγ0(s)∥s\mapsto\lVert G_{\boldsymbol\gamma_n}(s)-G_{\boldsymbol\gamma_0}(s)\rVertは可測であり、[a,b]∖N[a,b]\setminus Nで2m2m以下であるから、§E10.11 定理 3.1により∫ab∥Gγn−Gγ0∥→0\int_a^b\lVert G_{\boldsymbol\gamma_n}-G_{\boldsymbol\gamma_0}\rVert\to0である。▨

補題 3.7.d,p∈N≥1d,p\in\NNとし、開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dが状態とパラメータについてC1C^1級の Carathéodory 条件を満たすとし、φ:W→Rd\varphi:W\to\R^dをその解写像とする。(t,t0,η,λ)∈W(t,t_0,\boldsymbol\eta,\lambda)\in Wに対して、φ(⋅,t0,η,λ)\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)に沿う変分方程式の解でY(t0)=IdY(t_0)=I_d、Z(t0)=0Z(t_0)=0を満たすものの時刻ttにおける値をY(t,t0,η,λ)Y(t,t_0,\boldsymbol\eta,\lambda)、Z(t,t0,η,λ)Z(t,t_0,\boldsymbol\eta,\lambda)と書く。このときYYとZZはWW上で連続である。

証明.WWの点列(tn,t0,n,ηn,λn)(t_n,t_{0,n},\boldsymbol\eta_n,\lambda_n)が(t,t0,η,λ)∈W(t,t_0,\boldsymbol\eta,\lambda)\in Wへ収束するとする。a<min⁡{t,t0}a<\min\{t,t_0\}、max⁡{t,t0}<b\max\{t,t_0\}<b、[a,b]⊆I(t0,η,λ)[a,b]\subseteq I(t_0,\boldsymbol\eta,\lambda)を満たすa,ba,bを取る。y=φ(⋅,t0,η,λ)\boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)とおき、A,BA,Bをy\boldsymbol yに沿う変分方程式の係数とする。補題 3.6をy∣[a,b]\boldsymbol y|_{[a,b]}とλ0=λ\lambda_0=\lambdaに適用してρ>0\rho>0を取る。定理 2.2 (1)により、十分大きいnnではtn,t0,n∈(a,b)t_n,t_{0,n}\in(a,b)、[a,b]⊆I(t0,n,ηn,λn)[a,b]\subseteq I(t_{0,n},\boldsymbol\eta_n,\lambda_n)であり、yn=φ(⋅,t0,n,ηn,λn)\boldsymbol y_n=\varphi(\cdot,t_{0,n},\boldsymbol\eta_n,\lambda_n)は[a,b][a,b]上でy\boldsymbol yへ一様収束する。番号を付け直して、すべてのnnでこれらとmax⁡[a,b]∥yn−y∥≤ρ\max_{[a,b]}\lVert\boldsymbol y_n-\boldsymbol y\rVert\le\rho、∥λn−λ∥≤ρ\lVert\lambda_n-\lambda\rVert\le\rhoが成り立つとしてよい。

An,BnA_n,B_nをyn\boldsymbol y_nに沿う変分方程式の係数とする。γn(s,θ)=(yn(s),λn)\boldsymbol\gamma_n(s,\theta)=(\boldsymbol y_n(s),\lambda_n)は補題 3.6 (2)の仮定を満たし、偏導関数の値の規約により、[a,b][a,b]上でAγn=AnA_{\boldsymbol\gamma_n}=A_n、Bγn=BnB_{\boldsymbol\gamma_n}=B_n、Aγ0=AA_{\boldsymbol\gamma_0}=A、Bγ0=BB_{\boldsymbol\gamma_0}=Bである。n→∞n\to\inftyのときsup⁡θ∥γn(s,θ)−(y(s),λ)∥≤max⁡{max⁡[a,b]∥yn−y∥,∥λn−λ∥}→0\sup_\theta\lVert\boldsymbol\gamma_n(s,\theta)-(\boldsymbol y(s),\lambda)\rVert\le\max\{\max_{[a,b]}\lVert\boldsymbol y_n-\boldsymbol y\rVert,\lVert\lambda_n-\lambda\rVert\}\to0であるから、補題 3.6 (3)により∫ab∥An−A∥→0\int_a^b\lVert A_n-A\rVert\to0、∫ab∥Bn−B∥→0\int_a^b\lVert B_n-B\rVert\to0である。

Yn,ZnY_n,Z_nをyn\boldsymbol y_nに沿う変分方程式の解でYn(t0,n)=IdY_n(t_{0,n})=I_d、Zn(t0,n)=0Z_n(t_{0,n})=0を満たすものとすると、その[a,b][a,b]への制限は補題 3.1 (1)の解である。補題 3.1 (2)を行列値の方程式に適用すると、Yn→Y(⋅,t0,η,λ)Y_n\to Y(\cdot,t_0,\boldsymbol\eta,\lambda)、Zn→Z(⋅,t0,η,λ)Z_n\to Z(\cdot,t_0,\boldsymbol\eta,\lambda)は[a,b][a,b]上で一様収束する。従って

∥Yn(tn)−Y(t,t0,η,λ)∥≤max⁡[a,b]∥Yn−Y(⋅,t0,η,λ)∥+∥Y(tn,t0,η,λ)−Y(t,t0,η,λ)∥⟶0\lVert Y_n(t_n)-Y(t,t_0,\boldsymbol\eta,\lambda)\rVert\le\max_{[a,b]}\lVert Y_n-Y(\cdot,t_0,\boldsymbol\eta,\lambda)\rVert+\lVert Y(t_n,t_0,\boldsymbol\eta,\lambda)-Y(t,t_0,\boldsymbol\eta,\lambda)\rVert\longrightarrow0

であり、ZZについても同様である。▨

定理 3.8.d,p∈N≥1d,p\in\NNとし、開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dが状態とパラメータについてC1C^1級の Carathéodory 条件を満たすとし、φ:W→Rd\varphi:W\to\R^dをその解写像とする。このとき、各(t,t0)∈R2(t,t_0)\in\R^2について、Wt,t0={(η,λ):(t,t0,η,λ)∈W}W_{t,t_0}=\{(\boldsymbol\eta,\lambda):(t,t_0,\boldsymbol\eta,\lambda)\in W\}はRd×Rp\R^d\times\R^pの開集合であり、(η,λ)↦φ(t,t0,η,λ)(\boldsymbol\eta,\lambda)\mapsto\varphi(t,t_0,\boldsymbol\eta,\lambda)はWt,t0W_{t,t_0}上でC1C^1級であって、

Dηφ(t,t0,η,λ)=Y(t),Dλφ(t,t0,η,λ)=Z(t)D_{\boldsymbol\eta}\varphi(t,t_0,\boldsymbol\eta,\lambda)=Y(t),\qquad D_\lambda\varphi(t,t_0,\boldsymbol\eta,\lambda)=Z(t)

が成り立つ。ここでY,ZY,Zはφ(⋅,t0,η,λ)\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)に沿う変分方程式の解で、Y(t0)=IdY(t_0)=I_d、Z(t0)=0Z(t_0)=0を満たすものである。(t,t0,η,λ)↦Dηφ(t,t0,η,λ)(t,t_0,\boldsymbol\eta,\lambda)\mapsto D_{\boldsymbol\eta}\varphi(t,t_0,\boldsymbol\eta,\lambda)と(t,t0,η,λ)↦Dλφ(t,t0,η,λ)(t,t_0,\boldsymbol\eta,\lambda)\mapsto D_\lambda\varphi(t,t_0,\boldsymbol\eta,\lambda)はWW上で連続である。

証明.Wt,t0W_{t,t_0}は開集合WW(定理 2.2 (2))の切り口であるから開集合である。(t,t0,η,λ)∈W(t,t_0,\boldsymbol\eta,\lambda)\in Wとe=(e′,e′′)∈Rd×Rp\boldsymbol e=(\boldsymbol e',\boldsymbol e'')\in\R^d\times\R^pを取り、y=φ(⋅,t0,η,λ)\boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)とおき、A,BA,Bをy\boldsymbol yに沿う変分方程式の係数、NNを定義 3.2 条件 (b)の零集合とする。a<min⁡{t,t0}a<\min\{t,t_0\}、max⁡{t,t0}<b\max\{t,t_0\}<b、[a,b]⊆I(t0,η,λ)[a,b]\subseteq I(t_0,\boldsymbol\eta,\lambda)を満たすa,ba,bを取る。補題 3.6をy∣[a,b]\boldsymbol y|_{[a,b]}とλ0=λ\lambda_0=\lambdaに適用してρ>0\rho>0を取り、同補題のKρK_\rhoとRd×Rp\R^d\times\R^pのノルムを用いる。定理 2.2 (1)により、あるh0>0h_0>0が存在し、0<∣h∣≤h00<\lvert h\rvert\le h_0ならば[a,b]⊆I(t0,η+he′,λ+he′′)[a,b]\subseteq I(t_0,\boldsymbol\eta+h\boldsymbol e',\lambda+h\boldsymbol e'')であり、yh=φ(⋅,t0,η+he′,λ+he′′)\boldsymbol y_h=\varphi(\cdot,t_0,\boldsymbol\eta+h\boldsymbol e',\lambda+h\boldsymbol e'')はmax⁡[a,b]∥yh−y∥≤ρ\max_{[a,b]}\lVert\boldsymbol y_h-\boldsymbol y\rVert\le\rho、∣h∣∥e′′∥≤ρ\lvert h\rvert\lVert\boldsymbol e''\rVert\le\rhoを満たし、h→0h\to0のときmax⁡[a,b]∥yh−y∥→0\max_{[a,b]}\lVert\boldsymbol y_h-\boldsymbol y\rVert\to0である。

0<∣h∣≤h00<\lvert h\rvert\le h_0、s∈[a,b]s\in[a,b]、θ∈[0,1]\theta\in[0,1]に対して

γh(s,θ)=(y(s)+θ(yh(s)−y(s)), λ+θhe′′)\boldsymbol\gamma_h(s,\theta)=\bigl(\boldsymbol y(s)+\theta(\boldsymbol y_h(s)-\boldsymbol y(s)),\ \lambda+\theta h\boldsymbol e''\bigr)

とおくと、(s,γh(s,θ))∈Kρ(s,\boldsymbol\gamma_h(s,\theta))\in K_\rhoであり、γh\boldsymbol\gamma_hは補題 3.6 (2)の仮定を満たす。同項のAγhA_{\boldsymbol\gamma_h}、BγhB_{\boldsymbol\gamma_h}をAhA_h、BhB_hと書くと、同項によりAh,BhA_h,B_hは[a,b][a,b]上で可積分である。s∈[a,b]∖Ns\in[a,b]\setminus Nでは、§E4.3 定理 1.1によりθ↦f(s,γh(s,θ))\theta\mapsto f(s,\boldsymbol\gamma_h(s,\theta))は[0,1][0,1]上でC1C^1級であり、その導関数はDuf(s,γh(s,θ))(yh(s)−y(s))+hDλf(s,γh(s,θ))e′′D_{\boldsymbol u}f(s,\boldsymbol\gamma_h(s,\theta))(\boldsymbol y_h(s)-\boldsymbol y(s))+hD_\lambda f(s,\boldsymbol\gamma_h(s,\theta))\boldsymbol e''である。これを成分ごとに[0,1][0,1]で積分すると

f(s,yh(s),λ+he′′)−f(s,y(s),λ)=Ah(s)(yh(s)−y(s))+hBh(s)e′′f(s,\boldsymbol y_h(s),\lambda+h\boldsymbol e'')-f(s,\boldsymbol y(s),\lambda)=A_h(s)\bigl(\boldsymbol y_h(s)-\boldsymbol y(s)\bigr)+hB_h(s)\boldsymbol e''

を得る。

qh=(yh−y)/h\boldsymbol q_h=(\boldsymbol y_h-\boldsymbol y)/hは[a,b][a,b]上で絶対連続であり、qh(t0)=e′\boldsymbol q_h(t_0)=\boldsymbol e'である。yh\boldsymbol y_hとy\boldsymbol yがともに方程式を満たしs∉Ns\notin Nであるssでは、上の等式をhhで割ると

qh′(s)=Ah(s)qh(s)+Bh(s)e′′\boldsymbol q_h'(s)=A_h(s)\boldsymbol q_h(s)+B_h(s)\boldsymbol e''

である。従ってqh\boldsymbol q_hは、係数AhA_h、非斉次項Bhe′′B_h\boldsymbol e''、初期時刻t0t_0、初期値e′\boldsymbol e'に対する補題 3.1 (1)の解である。偏導関数の値の規約により、γ0(s,θ)=(y(s),λ)\boldsymbol\gamma_0(s,\theta)=(\boldsymbol y(s),\lambda)について[a,b][a,b]上でAγ0=AA_{\boldsymbol\gamma_0}=A、Bγ0=BB_{\boldsymbol\gamma_0}=Bである。s∈[a,b]s\in[a,b]では、h→0h\to0のときmax⁡θ∈[0,1]∥γh(s,θ)−(y(s),λ)∥≤max⁡{max⁡[a,b]∥yh−y∥,∣h∣∥e′′∥}→0\max_{\theta\in[0,1]}\lVert\boldsymbol\gamma_h(s,\theta)-(\boldsymbol y(s),\lambda)\rVert\le\max\{\max_{[a,b]}\lVert\boldsymbol y_h-\boldsymbol y\rVert,\lvert h\rvert\lVert\boldsymbol e''\rVert\}\to0である。零へ収束する任意の点列hnh_n(0<∣hn∣≤h00<\lvert h_n\rvert\le h_0)について、補題 3.6 (3)をγhn\boldsymbol\gamma_{h_n}に適用すると、∫ab∥Ahn−A∥→0\int_a^b\lVert A_{h_n}-A\rVert\to0であり、∫ab∥Bhne′′−Be′′∥≤∥e′′∥∫ab∥Bhn−B∥→0\int_a^b\lVert B_{h_n}\boldsymbol e''-B\boldsymbol e''\rVert\le\lVert\boldsymbol e''\rVert\int_a^b\lVert B_{h_n}-B\rVert\to0である。補題 3.1 (2)により、qhn\boldsymbol q_{h_n}は[a,b][a,b]上で、q′=Aq+Be′′\boldsymbol q'=A\boldsymbol q+B\boldsymbol e''、q(t0)=e′\boldsymbol q(t_0)=\boldsymbol e'の解q\boldsymbol qへ一様収束する。Ye′+Ze′′Y\boldsymbol e'+Z\boldsymbol e''はこの線形方程式の[a,b][a,b]上の解で時刻t0t_0にe′\boldsymbol e'をとるから、一意性によりq=Ye′+Ze′′\boldsymbol q=Y\boldsymbol e'+Z\boldsymbol e''である。点列hnh_nは任意であるから

lim⁡h→0φ(t,t0,η+he′,λ+he′′)−φ(t,t0,η,λ)h=Y(t)e′+Z(t)e′′\lim_{h\to0}\frac{\varphi(t,t_0,\boldsymbol\eta+h\boldsymbol e',\lambda+h\boldsymbol e'')-\varphi(t,t_0,\boldsymbol\eta,\lambda)}{h}=Y(t)\boldsymbol e'+Z(t)\boldsymbol e''

である。e\boldsymbol eをRd×Rp\R^d\times\R^pの標準基底の各ベクトルとすると、(η,λ)↦φ(t,t0,η,λ)(\boldsymbol\eta,\lambda)\mapsto\varphi(t,t_0,\boldsymbol\eta,\lambda)の偏導関数はY(t)Y(t)とZ(t)Z(t)の列である。補題 3.7により、これらの偏導関数はWW上で連続であり、特にWt,t0W_{t,t_0}上で連続である。§E4.2 系 2.2により、(η,λ)↦φ(t,t0,η,λ)(\boldsymbol\eta,\lambda)\mapsto\varphi(t,t_0,\boldsymbol\eta,\lambda)はWt,t0W_{t,t_0}上でC1C^1級であり、Dηφ=Y(t)D_{\boldsymbol\eta}\varphi=Y(t)、Dλφ=Z(t)D_\lambda\varphi=Z(t)である。▨

例 3.9.Ω=R×R×R\Omega=\R\times\R\times\R上のf(t,u,λ)=∣u∣f(t,u,\lambda)=\lvert u\rvertは連続であり、∣f(t,u,λ)−f(t,v,λ)∣≤∣u−v∣\lvert f(t,u,\lambda)-f(t,v,\lambda)\rvert\le\lvert u-v\rvertを満たすから、パラメータ付き Carathéodory 条件を局所的に満たす。η∈R\eta\in\Rに対して

y(t)={ηet,η≥0,ηe−t,η<0y(t)=\begin{cases}\eta e^{t},&\eta\ge0,\\ \eta e^{-t},&\eta<0\end{cases}

はR\R上でC1C^1級であり、η≥0\eta\ge0ではy≥0y\ge0、η<0\eta<0ではy<0y<0であるからy′=∣y∣y'=\lvert y\rvert、y(0)=ηy(0)=\etaを満たす。従ってI(0,η,λ)=RI(0,\eta,\lambda)=\Rであり、φ(t,0,η,λ)=y(t)\varphi(t,0,\eta,\lambda)=y(t)である。t≠0t\ne0を固定すると、η↦φ(t,0,η,λ)\eta\mapsto\varphi(t,0,\eta,\lambda)のη=0\eta=0における右微分係数はete^{t}、左微分係数はe−te^{-t}であり、一致しない。定理 3.8の結論は、状態についてのC1C^1級の仮定を Lipschitz 条件へ弱めると成り立たない。一方、定理 2.2 (2)により、このφ\varphiはWW上で連続である。

例 3.10.Ω=R×R×R\Omega=\R\times\R\times\R上で、(t,λ)≠(0,0)(t,\lambda)\ne(0,0)ならばf(t,u,λ)=λ/t2+λ2f(t,u,\lambda)=\lambda/\sqrt{t^2+\lambda^2}、f(0,u,0)=0f(0,u,0)=0とおく。ffはuuによらず、すべての点で∣f(t,u,λ)∣≤1\lvert f(t,u,\lambda)\rvert\le1である。各(u,λ)(u,\lambda)についてt↦f(t,u,λ)t\mapsto f(t,u,\lambda)はR∖{0}\R\setminus\{0\}上で連続であるから Lebesgue 可測であり、t≠0t\ne0では(u,λ)↦f(t,u,λ)(u,\lambda)\mapsto f(t,u,\lambda)はR2\R^2上でC∞C^\infty級である。従って、ffは各閉矩形J×B×CJ\times B\times Cについて零集合{0}∩J\{0\}\cap J、m=1m=1、ℓ=0\ell=0をもってパラメータ付き Carathéodory 条件を局所的に満たし、N={0}N=\{0\}として定義 3.2 条件 (a)と定義 3.2 条件 (b)をr=1r=1で満たす。t≠0t\ne0では∂uf(t,u,λ)=0\partial_uf(t,u,\lambda)=0であり

∂λf(t,u,λ)=t2(t2+λ2)3/2,∂λf(t,0,0)=1∣t∣\partial_\lambda f(t,u,\lambda)=\frac{t^2}{(t^2+\lambda^2)^{3/2}},\qquad\partial_\lambda f(t,0,0)=\frac1{\lvert t\rvert}

である。コンパクト集合K=[−1,1]×{0}×{0}K=[-1,1]\times\{0\}\times\{0\}について、どの零集合N′N'を取っても、t∈[−1,1]∖N′t\in[-1,1]\setminus N'でmK(t)≥1/∣t∣m_K(t)\ge1/\lvert t\rvertを満たすmKm_KはL1(R)L^1(\R)に属さない。従ってffは定義 3.2 条件 (c)を満たさず、状態とパラメータについてC1C^1級の Carathéodory 条件を満たさない。ffとその状態についての偏導関数は支配され、満たされないのはパラメータについての偏導関数の支配だけである。

(t0,η,λ)∈R3(t_0,\eta,\lambda)\in\R^3を取る。ffはuuによらないから、R\R上で定義されたy(t)=η+∫t0tf(s,0,λ) dsy(t)=\eta+\int_{t_0}^tf(s,0,\lambda)\,dsは、§E13.4 補題 2.4によりy′=f(t,y,λ)y'=f(t,y,\lambda)、y(t0)=ηy(t_0)=\etaのR\R上の絶対連続解である。従ってW=R4W=\R^4であり、φ(t,t0,η,λ)=y(t)\varphi(t,t_0,\eta,\lambda)=y(t)である。λ≠0\lambda\ne0ではddsarsinh⁡(s/∣λ∣)=1/s2+λ2\frac{d}{ds}\operatorname{arsinh}(s/\lvert\lambda\rvert)=1/\sqrt{s^2+\lambda^2}であるから

φ(t,0,η,λ)=η+λarsinh⁡t∣λ∣\varphi(t,0,\eta,\lambda)=\eta+\lambda\operatorname{arsinh}\frac{t}{\lvert\lambda\rvert}

であり、φ(t,0,η,0)=η\varphi(t,0,\eta,0)=\etaである。t>0t>0を固定すると、λ→0\lambda\to0のとき

φ(t,0,η,λ)−φ(t,0,η,0)λ=arsinh⁡t∣λ∣⟶∞\frac{\varphi(t,0,\eta,\lambda)-\varphi(t,0,\eta,0)}{\lambda}=\operatorname{arsinh}\frac{t}{\lvert\lambda\rvert}\longrightarrow\infty

であるから、λ↦φ(t,0,η,λ)\lambda\mapsto\varphi(t,0,\eta,\lambda)はλ=0\lambda=0で微分可能でない。定理 3.8の結論は、パラメータについての偏導関数の可積分な支配の仮定を外すと成り立たない。一方、∣λarsinh⁡(t/∣λ∣)∣→0\lvert\lambda\operatorname{arsinh}(t/\lvert\lambda\rvert)\rvert\to0であり、定理 2.2 (2)のとおりφ\varphiはWW上で連続である。

例 3.11.Ω=R×R×R\Omega=\R\times\R\times\R上のf(t,u,λ)=u2f(t,u,\lambda)=u^2は、補題 3.3 (4)によりすべてのr∈N≥1r\in\NNについて状態とパラメータについてCrC^r級の Carathéodory 条件を満たす。η∈R\eta\in\Rに対して、φ(t,0,η,λ)=η/(1−ηt)\varphi(t,0,\eta,\lambda)=\eta/(1-\eta t)であり、I(0,η,λ)={t:ηt<1}I(0,\eta,\lambda)=\{t:\eta t<1\}である。実際、関数y(t)=η/(1−ηt)y(t)=\eta/(1-\eta t)は{t:ηt<1}\{t:\eta t<1\}上でy′=y2y'=y^2、y(0)=ηy(0)=\etaを満たし、η≠0\eta\ne0ではt→1/ηt\to1/\etaのとき∣y(t)∣→∞\lvert y(t)\rvert\to\inftyであるから、これ以上延長することができない。t>0t>0を固定するとWt,0={(η,λ):η<1/t}W_{t,0}=\{(\eta,\lambda):\eta<1/t\}であり、Wt,0W_{t,0}はR2\R^2全体ではない。

Dηφ(t,0,η,λ)=(1−ηt)−2D_\eta\varphi(t,0,\eta,\lambda)=(1-\eta t)^{-2}であり、ηt↑1\eta t\uparrow1のとき非有界である。y=φ(⋅,0,η,λ)y=\varphi(\cdot,0,\eta,\lambda)に沿う変分方程式の係数はA(t)=2y(t)=2η/(1−ηt)A(t)=2y(t)=2\eta/(1-\eta t)、B(t)=0B(t)=0であり、Y(t)=(1−ηt)−2Y(t)=(1-\eta t)^{-2}は

Y′(t)=2η(1−ηt)3=A(t)Y(t),Y(0)=1Y'(t)=\frac{2\eta}{(1-\eta t)^3}=A(t)Y(t),\qquad Y(0)=1

を満たす。Z=0Z=0とDλφ=0D_\lambda\varphi=0も一致し、定理 3.8の等式が成り立つ。

4 高階の微分可能性と自律系の解写像

定理 4.1.d,p,r∈N≥1d,p,r\in\NNとし、開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dが状態とパラメータについてCrC^r級の Carathéodory 条件を満たすとし、φ:W→Rd\varphi:W\to\R^dをその解写像とする。このとき、各(t,t0)∈R2(t,t_0)\in\R^2について、(η,λ)↦φ(t,t0,η,λ)(\boldsymbol\eta,\lambda)\mapsto\varphi(t,t_0,\boldsymbol\eta,\lambda)は開集合Wt,t0={(η,λ):(t,t0,η,λ)∈W}W_{t,t_0}=\{(\boldsymbol\eta,\lambda):(t,t_0,\boldsymbol\eta,\lambda)\in W\}上でCrC^r級であり、(η,λ)(\boldsymbol\eta,\lambda)の座標についてのrr階以下の順序付き偏導関数はすべてWW上で連続である。

証明.r=1r=1の場合は定理 2.2 (2)と定理 3.8である。r≥2r\ge2とし、階数r−1r-1の場合の主張がすべてのd,pd,p、開集合と右辺について成り立つと仮定する。NNをffの零集合とする。

Ω^={(t,u,U,V,λ)∈R×Rd×Rd×d×Rd×p×Rp:(t,u,λ)∈Ω}\widehat\Omega=\{(t,\boldsymbol u,U,V,\lambda)\in\R\times\R^d\times\R^{d\times d}\times\R^{d\times p}\times\R^p:(t,\boldsymbol u,\lambda)\in\Omega\}

上で

F(t,u,U,V,λ)=(f(t,u,λ), Duf(t,u,λ)U, Duf(t,u,λ)V+Dλf(t,u,λ))F(t,\boldsymbol u,U,V,\lambda)=\bigl(f(t,\boldsymbol u,\lambda),\ D_{\boldsymbol u}f(t,\boldsymbol u,\lambda)U,\ D_{\boldsymbol u}f(t,\boldsymbol u,\lambda)V+D_\lambda f(t,\boldsymbol u,\lambda)\bigr)

と定める。状態を(u,U,V)(\boldsymbol u,U,V)、パラメータをλ\lambdaとみると、FFは状態とパラメータについてCr−1C^{r-1}級の Carathéodory 条件を満たす。この確認は演習とする(問題 6.2)。FFの解写像をφ^:W^→Rd×Rd×d×Rd×p\widehat\varphi:\widehat W\to\R^d\times\R^{d\times d}\times\R^{d\times p}とする。

(t0,η,λ)∈Ω(t_0,\boldsymbol\eta,\lambda)\in\Omega、I=I(t0,η,λ)I=I(t_0,\boldsymbol\eta,\lambda)、y=φ(⋅,t0,η,λ)\boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)とし、Y,ZY,Zをy\boldsymbol yに沿う変分方程式の解でY(t0)=IdY(t_0)=I_d、Z(t0)=0Z(t_0)=0を満たすものとする。(y,Y,Z)(\boldsymbol y,Y,Z)はII上で局所絶対連続であり、IIのほとんど至る所でFFの方程式を満たし、時刻t0t_0で(η,Id,0)(\boldsymbol\eta,I_d,0)をとる。§E13.6 定理 2.1により、(y,Y,Z)(\boldsymbol y,Y,Z)は初期値(η,Id,0)(\boldsymbol\eta,I_d,0)と時刻t0t_0、パラメータλ\lambdaに対するFFの最大絶対連続解y^\widehat{\boldsymbol y}の制限であり、y^\widehat{\boldsymbol y}の最大存在区間I^\widehat IはIIを含む。y^\widehat{\boldsymbol y}の第一成分はI^\widehat I上のy′=f(t,y,λ)\boldsymbol y'=f(t,\boldsymbol y,\lambda)の絶対連続解で時刻t0t_0にη\boldsymbol\etaをとるから、同定理によりy\boldsymbol yの制限であり、I^⊆I\widehat I\subseteq Iである。従ってI^=I\widehat I=Iであり、定理 3.8により、すべてのt∈It\in Iで

φ^(t,t0,(η,Id,0),λ)=(φ(t,t0,η,λ), Dηφ(t,t0,η,λ), Dλφ(t,t0,η,λ))\widehat\varphi\bigl(t,t_0,(\boldsymbol\eta,I_d,0),\lambda\bigr)=\bigl(\varphi(t,t_0,\boldsymbol\eta,\lambda),\ D_{\boldsymbol\eta}\varphi(t,t_0,\boldsymbol\eta,\lambda),\ D_\lambda\varphi(t,t_0,\boldsymbol\eta,\lambda)\bigr)

である。特に、(t,t0,η,λ)∈W(t,t_0,\boldsymbol\eta,\lambda)\in Wと(t,t0,(η,Id,0),λ)∈W^(t,t_0,(\boldsymbol\eta,I_d,0),\lambda)\in\widehat Wは同値である。

帰納法の仮定をFFに適用すると、φ^\widehat\varphiの初期値(η,U,V)(\boldsymbol\eta,U,V)とパラメータλ\lambdaの座標についてのr−1r-1階以下の順序付き偏導関数は存在してW^\widehat W上で連続である。U=IdU=I_d、V=0V=0に固定した写像の(η,λ)(\boldsymbol\eta,\lambda)の座標についての偏微分は、φ^\widehat\varphiの同じ座標についての偏微分の値を与える。従ってDηφD_{\boldsymbol\eta}\varphiとDλφD_\lambda\varphiの各成分の(η,λ)(\boldsymbol\eta,\lambda)の座標についてのr−1r-1階以下の順序付き偏導関数は存在し、連続写像(t,t0,η,λ)↦(t,t0,(η,Id,0),λ)(t,t_0,\boldsymbol\eta,\lambda)\mapsto(t,t_0,(\boldsymbol\eta,I_d,0),\lambda)との合成としてWW上で連続である。2≤j≤r2\le j\le rについて、φ\varphiのjj階の順序付き偏導関数は、一階の偏導関数であるDηφD_{\boldsymbol\eta}\varphi、DλφD_\lambda\varphiの成分のj−1j-1階の順序付き偏導関数であるから、存在してWW上で連続である。一階の偏導関数は定理 3.8により、φ\varphi自身は定理 2.2 (2)によりWW上で連続である。§E4.4 定理 1.7により、(η,λ)↦φ(t,t0,η,λ)(\boldsymbol\eta,\lambda)\mapsto\varphi(t,t_0,\boldsymbol\eta,\lambda)はWt,t0W_{t,t_0}上でCrC^r級である。▨

系 4.2.d,p,r∈N≥1d,p,r\in\NNとし、開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dと、(u,λ)(\boldsymbol u,\lambda)の座標についてのffのrr階以下のすべての順序付き偏導関数が、Ω\Omega上で存在して連続であるとする。補題 3.3 (4)と補題 3.3 (1)によりffはパラメータ付き Carathéodory 条件を局所的に満たすので、その解写像をφ:W→Rd\varphi:W\to\R^dとする。このとき、次が成り立つ。

  1. φ∈C1(W)\varphi\in C^1(W)であり、(η,λ)(\boldsymbol\eta,\lambda)の座標についてのrr階以下の順序付き偏導関数はWW上で連続である。すべての(t,t0,η,λ)∈W(t,t_0,\boldsymbol\eta,\lambda)\in Wで ∂tφ(t,t0,η,λ)=f(t,φ(t,t0,η,λ),λ),∂t0φ(t,t0,η,λ)=−Dηφ(t,t0,η,λ) f(t0,η,λ)\partial_t\varphi(t,t_0,\boldsymbol\eta,\lambda)=f\bigl(t,\varphi(t,t_0,\boldsymbol\eta,\lambda),\lambda\bigr),\qquad \partial_{t_0}\varphi(t,t_0,\boldsymbol\eta,\lambda)=-D_{\boldsymbol\eta}\varphi(t,t_0,\boldsymbol\eta,\lambda)\,f(t_0,\boldsymbol\eta,\lambda) が成り立つ。
  2. f∈Cr(Ω)f\in C^r(\Omega)ならばφ∈Cr(W)\varphi\in C^r(W)である。f∈C∞(Ω)f\in C^\infty(\Omega)ならばφ∈C∞(W)\varphi\in C^\infty(W)である。

証明.主張 (1)を示す。定理 4.1により、(η,λ)(\boldsymbol\eta,\lambda)の座標についてのrr階以下の順序付き偏導関数はWW上で連続である。(t0,η,λ)∈Ω(t_0,\boldsymbol\eta,\lambda)\in\Omegaに対してy=φ(⋅,t0,η,λ)\boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)とおくと、§E13.4 補題 2.4によりy(t)=η+∫t0tf(s,y(s),λ) ds\boldsymbol y(t)=\boldsymbol\eta+\int_{t_0}^tf(s,\boldsymbol y(s),\lambda)\,dsである。f(⋅,⋅,λ)f(\cdot,\cdot,\lambda)はΩλ={(t,u):(t,u,λ)∈Ω}\Omega^\lambda=\{(t,\boldsymbol u):(t,\boldsymbol u,\lambda)\in\Omega\}上で連続であるから、§E13.4 補題 2.5によりy\boldsymbol yは古典解であり、∂tφ(t,t0,η,λ)=f(t,φ(t,t0,η,λ),λ)\partial_t\varphi(t,t_0,\boldsymbol\eta,\lambda)=f(t,\varphi(t,t_0,\boldsymbol\eta,\lambda),\lambda)がすべての点で成り立つ。右辺は定理 2.2 (2)によりWW上で連続である。

(t,t0,η,λ)∈W(t,t_0,\boldsymbol\eta,\lambda)\in Wを取る。(t0,t0,η,λ)(t_0,t_0,\boldsymbol\eta,\lambda)と(t,t0,η,λ)(t,t_0,\boldsymbol\eta,\lambda)は開集合WWに属するから、あるh1>0h_1>0が存在し、∣h∣<h1\lvert h\rvert<h_1ならば(t0,t0+h,η,λ)∈W(t_0,t_0+h,\boldsymbol\eta,\lambda)\in W、(t,t0+h,η,λ)∈W(t,t_0+h,\boldsymbol\eta,\lambda)\in Wである。このときηh=φ(t0,t0+h,η,λ)\boldsymbol\eta_h=\varphi(t_0,t_0+h,\boldsymbol\eta,\lambda)とおく。φ(⋅,t0+h,η,λ)\varphi(\cdot,t_0+h,\boldsymbol\eta,\lambda)とφ(⋅,t0,ηh,λ)\varphi(\cdot,t_0,\boldsymbol\eta_h,\lambda)はともに、時刻t0t_0でηh\boldsymbol\eta_hをとる最大絶対連続解であるから、§E13.6 定理 2.1の一意性により一致する。従って

φ(t,t0+h,η,λ)=φ(t,t0,ηh,λ)\varphi(t,t_0+h,\boldsymbol\eta,\lambda)=\varphi(t,t_0,\boldsymbol\eta_h,\lambda)

である。古典解の積分表示により

ηh−ηh=−1h∫t0t0+hf(s,φ(s,t0+h,η,λ),λ) ds\frac{\boldsymbol\eta_h-\boldsymbol\eta}{h}=-\frac1h\int_{t_0}^{t_0+h}f\bigl(s,\varphi(s,t_0+h,\boldsymbol\eta,\lambda),\lambda\bigr)\,ds

であり、(s,h)↦f(s,φ(s,t0+h,η,λ),λ)(s,h)\mapsto f(s,\varphi(s,t_0+h,\boldsymbol\eta,\lambda),\lambda)は定理 2.2 (2)により(t0,0)(t_0,0)で連続であるから、h→0h\to0のとき右辺は−f(t0,η,λ)-f(t_0,\boldsymbol\eta,\lambda)へ収束する。η′↦φ(t,t0,η′,λ)\boldsymbol\eta'\mapsto\varphi(t,t_0,\boldsymbol\eta',\lambda)は定理 3.8によりη\boldsymbol\etaの近傍でC1C^1級であるから、§E4.3 定理 1.1により

∂t0φ(t,t0,η,λ)=ddh∣h=0φ(t,t0,ηh,λ)=−Dηφ(t,t0,η,λ) f(t0,η,λ)\partial_{t_0}\varphi(t,t_0,\boldsymbol\eta,\lambda)=\frac{d}{dh}\Big|_{h=0}\varphi(t,t_0,\boldsymbol\eta_h,\lambda)=-D_{\boldsymbol\eta}\varphi(t,t_0,\boldsymbol\eta,\lambda)\,f(t_0,\boldsymbol\eta,\lambda)

である。右辺は定理 3.8によりWW上で連続である。φ\varphiの(t,t0,η,λ)(t,t_0,\boldsymbol\eta,\lambda)の座標についての一階の偏導関数はすべて存在してWW上で連続であるから、§E4.2 系 2.2によりφ∈C1(W)\varphi\in C^1(W)である。

主張 (2)を示す。f∈Cr(Ω)f\in C^r(\Omega)とする。

Ω~={(σ,w,τ,s,λ)∈R×Rd×R×R×Rp:(τ,w,λ)∈Ω}\widetilde\Omega=\{(\sigma,\boldsymbol w,\tau,s,\lambda)\in\R\times\R^d\times\R\times\R\times\R^p:(\tau,\boldsymbol w,\lambda)\in\Omega\}

上で、状態を(w,τ)(\boldsymbol w,\tau)、パラメータを(s,λ)(s,\lambda)とする写像

H(σ,w,τ,s,λ)=(s f(τ,w,λ), s)H(\sigma,\boldsymbol w,\tau,s,\lambda)=\bigl(s\,f(\tau,\boldsymbol w,\lambda),\ s\bigr)

を考える。HHはσ\sigmaによらず、ssについて一次であるから、ssについて二回以上偏微分したHHの順序付き偏導関数は零である。従って、(w,τ,s,λ)(\boldsymbol w,\tau,s,\lambda)の座標についてのHHの順序付き偏導関数は、ffの順序付き偏導関数にssを掛けたもの、ffの順序付き偏導関数、定数、零のいずれかである。従ってrr階以下の順序付き偏導関数はΩ~\widetilde\Omega上で連続であり、補題 3.3 (4)によりHHは状態とパラメータについてCrC^r級の Carathéodory 条件を満たす。HHの解写像をχ:W~→Rd×R\chi:\widetilde W\to\R^d\times\Rとする。

(t,t0,η,λ)∈W(t,t_0,\boldsymbol\eta,\lambda)\in Wを取り、s=t−t0s=t-t_0とおく。T={σ∈R:t0+sσ∈I(t0,η,λ)}T=\{\sigma\in\R:t_0+s\sigma\in I(t_0,\boldsymbol\eta,\lambda)\}は00と11を含む開区間であり、TT上で

w^(σ)=(φ(t0+sσ,t0,η,λ), t0+sσ)\widehat{\boldsymbol w}(\sigma)=\bigl(\varphi(t_0+s\sigma,t_0,\boldsymbol\eta,\lambda),\ t_0+s\sigma\bigr)

とおく。主張 (1)の∂tφ\partial_t\varphiの式により、w^\widehat{\boldsymbol w}はC1C^1級でw^′=H(σ,w^,s,λ)\widehat{\boldsymbol w}'=H(\sigma,\widehat{\boldsymbol w},s,\lambda)、w^(0)=(η,t0)\widehat{\boldsymbol w}(0)=(\boldsymbol\eta,t_0)を満たす。§E13.6 定理 2.1によりw^\widehat{\boldsymbol w}はHHの最大絶対連続解の制限であるから、(1,0,(η,t0),(s,λ))∈W~(1,0,(\boldsymbol\eta,t_0),(s,\lambda))\in\widetilde Wであり

χ(1,0,(η,t0),(t−t0,λ))=(φ(t,t0,η,λ), t)\chi\bigl(1,0,(\boldsymbol\eta,t_0),(t-t_0,\lambda)\bigr)=\bigl(\varphi(t,t_0,\boldsymbol\eta,\lambda),\ t\bigr)

である。線形写像Λ(t,t0,η,λ)=((η,t0),(t−t0,λ))\Lambda(t,t_0,\boldsymbol\eta,\lambda)=((\boldsymbol\eta,t_0),(t-t_0,\lambda))はWWを開集合W~1,0={(η^,μ):(1,0,η^,μ)∈W~}\widetilde W_{1,0}=\{(\widehat{\boldsymbol\eta},\mu):(1,0,\widehat{\boldsymbol\eta},\mu)\in\widetilde W\}へ写し、φ\varphiはχ(1,0,⋅)∘Λ\chi(1,0,\cdot)\circ\Lambdaの第一成分である。定理 4.1によりχ(1,0,⋅)\chi(1,0,\cdot)はW~1,0\widetilde W_{1,0}上でCrC^r級である。Λ\Lambdaは線形であるから、§E4.3 定理 1.1を反復して適用すると、φ\varphiのrr階以下の順序付き偏導関数はχ(1,0,⋅)\chi(1,0,\cdot)の同じ階数の順序付き偏導関数とΛ\Lambdaの合成の定数係数の一次結合であり、WW上で連続である。§E4.4 定理 1.7によりφ∈Cr(W)\varphi\in C^r(W)である。f∈C∞(Ω)f\in C^\infty(\Omega)ならば、各r∈N≥1r\in\NNに適用してφ∈C∞(W)\varphi\in C^\infty(W)を得る。▨

注意 4.3. 状態とパラメータについてCrC^r級の Carathéodory 条件の下では、解写像は初期時刻について微分可能であるとは限らない。実際、§E13.4 例 5.1の係数aa(t<0t<0でa(t)=1a(t)=1、t≥0t\ge0でa(t)=2a(t)=2)を用いて、Ω=R×R×R\Omega=\R\times\R\times\R上でf(t,u,λ)=a(t)uf(t,u,\lambda)=a(t)uとおく。ffはN=∅N=\varnothingと、コンパクト集合KKに対するmK=2max⁡(t,u,λ)∈K(1+∣u∣) 1π(K)m_K=2\max_{(t,u,\lambda)\in K}(1+\lvert u\rvert)\,\mathbf 1_{\pi(K)}(π(K)\pi(K)はKKの時間成分への射影)をもって、すべてのr∈N≥1r\in\NNについて状態とパラメータについてCrC^r級の Carathéodory 条件を満たす。(t0,η,λ)∈R3(t_0,\eta,\lambda)\in\R^3に対して、y(t)=ηexp⁡∫t0ta(s) dsy(t)=\eta\exp\int_{t_0}^ta(s)\,dsは局所 Lipschitz 連続であり、t≠0t\ne0でy′(t)=a(t)y(t)y'(t)=a(t)y(t)を満たすから、初期値問題y′=a(t)yy'=a(t)y、y(t0)=ηy(t_0)=\etaのR\R上の絶対連続解である。従ってW=R4W=\R^4であり

φ(t,t0,η,λ)=ηexp⁡∫t0ta(s) ds\varphi(t,t_0,\eta,\lambda)=\eta\exp\int_{t_0}^ta(s)\,ds

である。η≠0\eta\ne0ならば、t0↦φ(t,t0,η,λ)t_0\mapsto\varphi(t,t_0,\eta,\lambda)のt0=0t_0=0における左微分係数は−φ(t,0,η,λ)-\varphi(t,0,\eta,\lambda)、右微分係数は−2φ(t,0,η,λ)-2\varphi(t,0,\eta,\lambda)であり、両者は異なる。系 4.2 (1)の式の右辺−Dηφ(t,t0,η,λ)f(t0,η,λ)=−a(t0)φ(t,t0,η,λ)-D_\eta\varphi(t,t_0,\eta,\lambda)f(t_0,\eta,\lambda)=-a(t_0)\varphi(t,t_0,\eta,\lambda)は、この例ではt0=0t_0=0の左右で異なる極限をもつ。

系 4.4.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、F:U→RnF:U\to\R^nを局所 Lipschitz 連続な写像とする。ξ∈U\boldsymbol\xi\in Uに対して、x′=F(x)\boldsymbol x'=F(\boldsymbol x)、x(0)=ξ\boldsymbol x(0)=\boldsymbol\xiの最大解をφ(⋅,ξ)\varphi(\cdot,\boldsymbol\xi)、その最大存在区間をI(ξ)I(\boldsymbol\xi)とし、D={(t,ξ):ξ∈U, t∈I(ξ)}D=\{(t,\boldsymbol\xi):\boldsymbol\xi\in U,\ t\in I(\boldsymbol\xi)\}とおく。このとき、次が成り立つ。

  1. DDはR×Rn\R\times\R^nの開集合であり、(t,ξ)↦φ(t,ξ)(t,\boldsymbol\xi)\mapsto\varphi(t,\boldsymbol\xi)はDD上で連続である。
  2. r∈N≥1r\in\NNとし、F∈Cr(U)F\in C^r(U)ならば、φ∈Cr(D)\varphi\in C^r(D)であり、∂tφ(t,ξ)=F(φ(t,ξ))\partial_t\varphi(t,\boldsymbol\xi)=F(\varphi(t,\boldsymbol\xi))である。さらに、各ξ∈U\boldsymbol\xi\in Uについてt↦Dξφ(t,ξ)t\mapsto D_{\boldsymbol\xi}\varphi(t,\boldsymbol\xi)はY′=DF(φ(t,ξ))YY'=DF(\varphi(t,\boldsymbol\xi))Y、Y(0)=InY(0)=I_nのI(ξ)I(\boldsymbol\xi)上の解である。F∈C∞(U)F\in C^\infty(U)ならばφ∈C∞(D)\varphi\in C^\infty(D)である。

証明.Ω=R×U×R\Omega=\R\times U\times\R上でf(t,u,λ)=F(u)f(t,\boldsymbol u,\lambda)=F(\boldsymbol u)とおく(d=nd=n、p=1p=1)。Ω\Omegaに含まれる閉矩形J×B×CJ\times B\times CについてJ×B⊆R×UJ\times B\subseteq\R\times Uであり、§E13.6 補題 5.1を(t,u)↦F(u)(t,\boldsymbol u)\mapsto F(\boldsymbol u)に適用すると、J×BJ\times B上の零集合、支配関数と Lipschitz 係数が得られる。これらはλ\lambdaによらないから、ffはパラメータ付き Carathéodory 条件を局所的に満たす。ffの解写像をφf:W→Rn\varphi_f:W\to\R^nとする。同補題によりΩ0=R×U\Omega^0=\R\times U上の最大絶対連続解は最大解であるから、φ(t,ξ)=φf(t,0,ξ,0)\varphi(t,\boldsymbol\xi)=\varphi_f(t,0,\boldsymbol\xi,0)であり、DDは連続写像ι(t,ξ)=(t,0,ξ,0)\iota(t,\boldsymbol\xi)=(t,0,\boldsymbol\xi,0)によるWWの逆像である。

主張 (1)を示す。定理 2.2 (2)によりWWは開集合でφf\varphi_fは連続であるから、D=ι−1(W)D=\iota^{-1}(W)は開集合であり、φ=φf∘ι\varphi=\varphi_f\circ\iotaは連続である。

主張 (2)を示す。F∈Cr(U)F\in C^r(U)ならばf∈Cr(Ω)f\in C^r(\Omega)であり、系 4.2 (2)によりφf∈Cr(W)\varphi_f\in C^r(W)である。ι\iotaは座標の一部を定数に固定する写像であるから、φ=φf∘ι\varphi=\varphi_f\circ\iotaの順序付き偏導関数はφf\varphi_fの対応する順序付き偏導関数とι\iotaの合成であり、§E4.4 定理 1.7によりφ∈Cr(D)\varphi\in C^r(D)である。系 4.2 (1)により∂tφ(t,ξ)=f(t,φ(t,ξ),0)=F(φ(t,ξ))\partial_t\varphi(t,\boldsymbol\xi)=f(t,\varphi(t,\boldsymbol\xi),0)=F(\varphi(t,\boldsymbol\xi))である。定理 3.8によりDξφ(t,ξ)=Dηφf(t,0,ξ,0)D_{\boldsymbol\xi}\varphi(t,\boldsymbol\xi)=D_{\boldsymbol\eta}\varphi_f(t,0,\boldsymbol\xi,0)は、φ(⋅,ξ)\varphi(\cdot,\boldsymbol\xi)に沿う変分方程式Y′=Duf(t,φ(t,ξ),0)Y=DF(φ(t,ξ))YY'=D_{\boldsymbol u}f(t,\varphi(t,\boldsymbol\xi),0)Y=DF(\varphi(t,\boldsymbol\xi))Yの解でY(0)=InY(0)=I_nを満たすものである。F∈C∞(U)F\in C^\infty(U)ならば、各rrに適用してφ∈C∞(D)\varphi\in C^\infty(D)を得る。▨

5 線形系の複素パラメータについての正則性

命題 5.1.d∈N≥1d\in\NNとし、Cd\C^dにノルム∥⋅∥\lVert\cdot\rVertを固定し、d×dd\times d複素行列にはその作用素ノルムを用いる。I⊆RI\subseteq\Rを内部が空でない区間とし、A0,A1∈Lloc1(I;Cd×d)A_0,A_1\in L^1_{\mathrm{loc}}(I;\C^{d\times d})、b0,b1∈Lloc1(I;Cd)\boldsymbol b_0,\boldsymbol b_1\in L^1_{\mathrm{loc}}(I;\C^d)、t0∈It_0\in I、η∈Cd\boldsymbol\eta\in\C^dとする。このとき、各λ∈C\lambda\in\Cについて、IIのほとんど至る所で

x′=(A0(t)+λA1(t))x+b0(t)+λb1(t)\boldsymbol x'=\bigl(A_0(t)+\lambda A_1(t)\bigr)\boldsymbol x+\boldsymbol b_0(t)+\lambda\boldsymbol b_1(t)

を満たしx(t0)=η\boldsymbol x(t_0)=\boldsymbol\etaを満たすx(⋅,λ)∈ACloc(I;Cd)\boldsymbol x(\cdot,\lambda)\in AC_{\mathrm{loc}}(I;\C^d)が一意に存在する。Picard 反復を、t∈It\in I、λ∈C\lambda\in\Cに対してx0(t,λ)=η\boldsymbol x_0(t,\lambda)=\boldsymbol\etaと

xk+1(t,λ)=η+∫t0t((A0(s)+λA1(s))xk(s,λ)+b0(s)+λb1(s)) ds\boldsymbol x_{k+1}(t,\lambda)=\boldsymbol\eta+\int_{t_0}^t\Bigl(\bigl(A_0(s)+\lambda A_1(s)\bigr)\boldsymbol x_k(s,\lambda)+\boldsymbol b_0(s)+\lambda\boldsymbol b_1(s)\Bigr)\,ds

によって定める。このとき、次が成り立つ。

  1. (t,λ)↦x(t,λ)(t,\lambda)\mapsto\boldsymbol x(t,\lambda)はI×CI\times\C上で連続である。
  2. 各t∈It\in I、λ∈C\lambda\in\Cについて、複素微分係数∂λx(t,λ)=lim⁡h→0(x(t,λ+h)−x(t,λ))/h\partial_\lambda\boldsymbol x(t,\lambda)=\lim_{h\to0}(\boldsymbol x(t,\lambda+h)-\boldsymbol x(t,\lambda))/h(h∈C∖{0}h\in\C\setminus\{0\})が存在する。各λ\lambdaについて、t↦∂λx(t,λ)t\mapsto\partial_\lambda\boldsymbol x(t,\lambda)は z′=(A0(t)+λA1(t))z+A1(t)x(t,λ)+b1(t),z(t0)=0\boldsymbol z'=\bigl(A_0(t)+\lambda A_1(t)\bigr)\boldsymbol z+A_1(t)\boldsymbol x(t,\lambda)+\boldsymbol b_1(t),\qquad\boldsymbol z(t_0)=\boldsymbol 0 のII上の局所絶対連続な解である。
  3. 各t∈It\in Iについて、λ↦x(t,λ)\lambda\mapsto\boldsymbol x(t,\lambda)は整関数である。
  4. IIに含まれる任意の有限閉区間JJとR>0R>0について、xk\boldsymbol x_kはJ×{λ∈C:∣λ∣≤R}J\times\{\lambda\in\C:\lvert\lambda\rvert\le R\}上でx\boldsymbol xへ一様収束する。
  5. A0,A1,b0,b1A_0,A_1,\boldsymbol b_0,\boldsymbol b_1が実数値でありη∈Rd\boldsymbol\eta\in\R^dならば、すべてのt∈It\in Iとλ∈R\lambda\in\Rについてx(t,λ)∈Rd\boldsymbol x(t,\lambda)\in\R^dである。

証明.補題 3.1 (3)をK=C\mathbb K=\C、係数A0+λA1A_0+\lambda A_1、非斉次項b0+λb1\boldsymbol b_0+\lambda\boldsymbol b_1に適用すると、x(⋅,λ)\boldsymbol x(\cdot,\lambda)は一意に存在する。t0t_0を含みIIに含まれる長さが正の有限閉区間JJ上で、x(⋅,λ)∣J\boldsymbol x(\cdot,\lambda)|_Jは補題 3.1 (1)の解である。

主張 (1)を示す。(t,λ)∈I×C(t,\lambda)\in I\times\Cを取り、t0t_0を含み、IIにおけるttの近傍を含むJJを取る。λn→λ\lambda_n\to\lambdaならば

∫J∥(A0+λnA1)−(A0+λA1)∥=∣λn−λ∣∫J∥A1∥→0\int_J\bigl\lVert(A_0+\lambda_nA_1)-(A_0+\lambda A_1)\bigr\rVert=\lvert\lambda_n-\lambda\rvert\int_J\lVert A_1\rVert\to0

であり、非斉次項についても同様であるから、補題 3.1 (2)によりx(⋅,λn)→x(⋅,λ)\boldsymbol x(\cdot,\lambda_n)\to\boldsymbol x(\cdot,\lambda)はJJ上で一様収束する。(tn,λn)→(t,λ)(t_n,\lambda_n)\to(t,\lambda)ならば、十分大きいnnでtn∈Jt_n\in Jであり、∥x(tn,λn)−x(t,λ)∥≤max⁡J∥x(⋅,λn)−x(⋅,λ)∥+∥x(tn,λ)−x(t,λ)∥→0\lVert\boldsymbol x(t_n,\lambda_n)-\boldsymbol x(t,\lambda)\rVert\le\max_J\lVert\boldsymbol x(\cdot,\lambda_n)-\boldsymbol x(\cdot,\lambda)\rVert+\lVert\boldsymbol x(t_n,\lambda)-\boldsymbol x(t,\lambda)\rVert\to0である。

主張 (2)を示す。λ∈C\lambda\in\Cを固定し、z\boldsymbol zを主張の方程式のII上の解とする(A1x(⋅,λ)+b1∈Lloc1A_1\boldsymbol x(\cdot,\lambda)+\boldsymbol b_1\in L^1_{\mathrm{loc}}であるから補題 3.1 (3)により存在する)。h∈C∖{0}h\in\C\setminus\{0\}に対してqh=(x(⋅,λ+h)−x(⋅,λ))/h\boldsymbol q_h=(\boldsymbol x(\cdot,\lambda+h)-\boldsymbol x(\cdot,\lambda))/hとおくと、qh(t0)=0\boldsymbol q_h(t_0)=\boldsymbol 0であり、ほとんど至る所で

qh′=(A0+λA1)qh+A1x(⋅,λ+h)+b1\boldsymbol q_h'=\bigl(A_0+\lambda A_1\bigr)\boldsymbol q_h+A_1\boldsymbol x(\cdot,\lambda+h)+\boldsymbol b_1

である。JJをt0t_0を含みIIに含まれる長さが正の有限閉区間とする。h→0h\to0のとき、主張 (1)の証明によりmax⁡J∥x(⋅,λ+h)−x(⋅,λ)∥→0\max_J\lVert\boldsymbol x(\cdot,\lambda+h)-\boldsymbol x(\cdot,\lambda)\rVert\to0であるから、∫J∥A1(x(⋅,λ+h)−x(⋅,λ))∥→0\int_J\lVert A_1(\boldsymbol x(\cdot,\lambda+h)-\boldsymbol x(\cdot,\lambda))\rVert\to0である。零へ収束する任意の点列hnh_nに補題 3.1 (2)を適用すると、qhn→z\boldsymbol q_{h_n}\to\boldsymbol zはJJ上で一様収束する。JJと点列は任意であるから、各t∈It\in Iで∂λx(t,λ)\partial_\lambda\boldsymbol x(t,\lambda)が存在してz(t)\boldsymbol z(t)に等しい。

主張 (3)を示す。主張 (2)により、各ttについてλ↦x(t,λ)\lambda\mapsto\boldsymbol x(t,\lambda)の各成分はC\Cの各点で複素微分可能である。

主張 (4)を示す。JJとt0t_0を含むIIの最小の閉区間をJ′J'とする。帰納法により、各xk(⋅,λ)\boldsymbol x_k(\cdot,\lambda)はII上で連続であり、Picard 反復の被積分関数は局所可積分である。J′J'上でα=∥A0∥+R∥A1∥\alpha=\lVert A_0\rVert+R\lVert A_1\rVert、β=∥b0∥+R∥b1∥\beta=\lVert\boldsymbol b_0\rVert+R\lVert\boldsymbol b_1\rVert、H(t)=∣∫t0tα∣H(t)=\lvert\int_{t_0}^t\alpha\rvert、C0=∥η∥∫J′α+∫J′βC_0=\lVert\boldsymbol\eta\rVert\int_{J'}\alpha+\int_{J'}\betaとおく。t∈J′t\in J'、∣λ∣≤R\lvert\lambda\rvert\le R、k≥0k\ge0について

∥xk+1(t,λ)−xk(t,λ)∥≤C0H(t)kk!\lVert\boldsymbol x_{k+1}(t,\lambda)-\boldsymbol x_k(t,\lambda)\rVert\le C_0\frac{H(t)^k}{k!}

が成り立つ。実際、k=0k=0では、§E13.2 補題 1.3 (2)をCd\C^dの固定したノルムに適用すると、左辺は∣∫t0t(α∥η∥+β)∣≤C0\lvert\int_{t_0}^t(\alpha\lVert\boldsymbol\eta\rVert+\beta)\rvert\le C_0以下である。kkで成り立つとすると、xk+2−xk+1=∫t0t(A0+λA1)(xk+1−xk) ds\boldsymbol x_{k+2}-\boldsymbol x_{k+1}=\int_{t_0}^t(A_0+\lambda A_1)(\boldsymbol x_{k+1}-\boldsymbol x_k)\,dsと同じ評価により、左辺はC0∣∫t0tα(s)H(s)k/k! ds∣C_0\lvert\int_{t_0}^t\alpha(s)H(s)^k/k!\,ds\rvert以下である。t≥t0t\ge t_0の側ではHHは絶対連続でほとんど至る所H′=αH'=\alphaであり、§E13.2 定理 3.1 (1)と§E13.2 定理 2.5により∫t0tαHk/k!=H(t)k+1/(k+1)!\int_{t_0}^t\alpha H^k/k!=H(t)^{k+1}/(k+1)!である。t≤t0t\le t_0の側も同様である。H∗=∫J′αH_*=\int_{J'}\alphaとおくと∑kC0H∗k/k!<∞\sum_kC_0H_*^k/k!<\inftyであるから、(xk)(\boldsymbol x_k)はJ′×{∣λ∣≤R}J'\times\{\lvert\lambda\rvert\le R\}上で一様に Cauchy 列であり、一様収束の極限x~\widetilde{\boldsymbol x}をもつ。∣λ∣≤R\lvert\lambda\rvert\le Rについて

∥∫t0t(A0+λA1)(xk(s,λ)−x~(s,λ)) ds∥≤H∗max⁡J′∥xk(⋅,λ)−x~(⋅,λ)∥→0\left\lVert\int_{t_0}^t\bigl(A_0+\lambda A_1\bigr)\bigl(\boldsymbol x_k(s,\lambda)-\widetilde{\boldsymbol x}(s,\lambda)\bigr)\,ds\right\rVert\le H_*\max_{J'}\lVert\boldsymbol x_k(\cdot,\lambda)-\widetilde{\boldsymbol x}(\cdot,\lambda)\rVert\to0

が§E13.2 補題 1.3 (2)により成り立つから、Picard 反復の定義式でk→∞k\to\inftyとすると、連続関数x~(⋅,λ)\widetilde{\boldsymbol x}(\cdot,\lambda)はすべてのt∈J′t\in J'で

x~(t,λ)=η+∫t0t((A0(s)+λA1(s))x~(s,λ)+b0(s)+λb1(s)) ds\widetilde{\boldsymbol x}(t,\lambda)=\boldsymbol\eta+\int_{t_0}^t\Bigl(\bigl(A_0(s)+\lambda A_1(s)\bigr)\widetilde{\boldsymbol x}(s,\lambda)+\boldsymbol b_0(s)+\lambda\boldsymbol b_1(s)\Bigr)\,ds

を満たす。§E13.2 系 2.6によりx~(⋅,λ)\widetilde{\boldsymbol x}(\cdot,\lambda)はJ′J'上で絶対連続であり、ほとんど至る所で方程式を満たし、x~(t0,λ)=η\widetilde{\boldsymbol x}(t_0,\lambda)=\boldsymbol\etaである。補題 3.1 (1)の一意性によりx~(⋅,λ)=x(⋅,λ)∣J′\widetilde{\boldsymbol x}(\cdot,\lambda)=\boldsymbol x(\cdot,\lambda)|_{J'}である。

主張 (5)を示す。A0,A1,b0,b1A_0,A_1,\boldsymbol b_0,\boldsymbol b_1が実数値であり、η∈Rd\boldsymbol\eta\in\R^d、λ∈R\lambda\in\Rであるとする。補題 3.1 (3)をK=R\mathbb K=\R、Rd\R^dの任意のノルム、実係数A0+λA1A_0+\lambda A_1、b0+λb1\boldsymbol b_0+\lambda\boldsymbol b_1に適用すると、IIのほとんど至る所で方程式を満たし時刻t0t_0にη\boldsymbol\etaをとるx~∈ACloc(I;Rd)\widetilde{\boldsymbol x}\in AC_{\mathrm{loc}}(I;\R^d)が存在する。x~\widetilde{\boldsymbol x}はCd\C^dに値をとる同じ初期値問題の局所絶対連続な解でもあるから、x(⋅,λ)\boldsymbol x(\cdot,\lambda)の一意性によりx(⋅,λ)=x~\boldsymbol x(\cdot,\lambda)=\widetilde{\boldsymbol x}であり、すべてのt∈It\in Iでx(t,λ)∈Rd\boldsymbol x(t,\lambda)\in\R^dである。▨

例 5.2.λ∈C\lambda\in\Cに対して、初期値問題y′′+λy=0y''+\lambda y=0、y(0)=1y(0)=1、y′(0)=0y'(0)=0を、x=(y,y′)\boldsymbol x=(y,y')についての一階の系

x′=(A0+λA1)x,x(0)=(1,0),A0=(0100),A1=(00−10)\boldsymbol x'=\bigl(A_0+\lambda A_1\bigr)\boldsymbol x,\qquad\boldsymbol x(0)=(1,0),\qquad A_0=\begin{pmatrix}0&1\\0&0\end{pmatrix},\quad A_1=\begin{pmatrix}0&0\\-1&0\end{pmatrix}

として扱う。d=2d=2、I=RI=\R、t0=0t_0=0、b0=b1=0\boldsymbol b_0=\boldsymbol b_1=\boldsymbol 0として命題 5.1を適用し、その解をx(t,λ)=(x1(t,λ),x2(t,λ))\boldsymbol x(t,\lambda)=(x_1(t,\lambda),x_2(t,\lambda))と書く。κ>0\kappa>0とする。このとき、R\R上で連続微分可能な写像

t↦(cos⁡κt, −κsin⁡κt),t↦(cosh⁡κt, κsinh⁡κt)t\mapsto(\cos\kappa t,\ -\kappa\sin\kappa t),\qquad t\mapsto(\cosh\kappa t,\ \kappa\sinh\kappa t)

は、それぞれλ=κ2\lambda=\kappa^2、λ=−κ2\lambda=-\kappa^2の系をすべてのttで満たし、時刻00で(1,0)(1,0)をとる。定数写像t↦(1,0)t\mapsto(1,0)はλ=0\lambda=0の系を満たす。§E13.2 系 2.7 (2)によりこれらは局所絶対連続であるから、同命題の解の一意性により、すべてのt∈Rt\in\Rで

x1(t,λ)={cos⁡(λ t),λ>0,1,λ=0,cosh⁡(−λ t),λ<0x_1(t,\lambda)=\begin{cases}\cos(\sqrt\lambda\,t),&\lambda>0,\\1,&\lambda=0,\\\cosh(\sqrt{-\lambda}\,t),&\lambda<0\end{cases}

である。実軸上のこの式は∣λ∣\sqrt{\lvert\lambda\rvert}を通じてλ\lambdaに依存するが、命題 5.1 (3)により、各t∈Rt\in\Rについてλ↦x1(t,λ)\lambda\mapsto x_1(t,\lambda)は整関数である。

命題 5.1 (2)により、t↦∂λx(t,0)t\mapsto\partial_\lambda\boldsymbol x(t,0)はz′=A0z+A1x(t,0)\boldsymbol z'=A_0\boldsymbol z+A_1\boldsymbol x(t,0)、z(0)=0\boldsymbol z(0)=\boldsymbol 0のR\R上の局所絶対連続な解である。A1x(t,0)=(0,−1)A_1\boldsymbol x(t,0)=(0,-1)であり、連続微分可能な写像z(t)=(−t2/2, −t)\boldsymbol z(t)=(-t^2/2,\,-t)はすべてのttでz′(t)=(−t,−1)=A0z(t)+(0,−1)\boldsymbol z'(t)=(-t,-1)=A_0\boldsymbol z(t)+(0,-1)とz(0)=0\boldsymbol z(0)=\boldsymbol 0を満たす。§E13.2 系 2.7 (2)と補題 3.1 (3)の一意性により、∂λx1(t,0)=−t2/2\partial_\lambda x_1(t,0)=-t^2/2である。

6 演習

問題 6.1.Ω=R×R×R\Omega=\R\times\R\times\R上のf(t,u,λ)=λ+u2f(t,u,\lambda)=\lambda+u^2の解写像をφ\varphiとし、初期値問題y′=λ+y2y'=\lambda+y^2、y(0)=0y(0)=0を考える。

  1. 定理 3.8を用いて、すべてのt∈Rt\in\Rについて∂λφ(t,0,0,0)=t\partial_\lambda\varphi(t,0,0,0)=tであることを示せ。
  2. λ>0\lambda>0ではφ(t,0,0,λ)=λtan⁡(λ t)\varphi(t,0,0,\lambda)=\sqrt\lambda\tan(\sqrt\lambda\,t)(∣t∣<π/(2λ)\lvert t\rvert<\pi/(2\sqrt\lambda))、λ<0\lambda<0ではφ(t,0,0,λ)=−−λtanh⁡(−λ t)\varphi(t,0,0,\lambda)=-\sqrt{-\lambda}\tanh(\sqrt{-\lambda}\,t)(t∈Rt\in\R)であることを示し、ttを固定してλ→0\lambda\to0のときφ(t,0,0,λ)/λ→t\varphi(t,0,0,\lambda)/\lambda\to tとなることを直接確かめよ。
解答.

(1)を示す。ffはC∞C^\infty級であるから、補題 3.3 (4)により状態とパラメータについてC1C^1級の Carathéodory 条件を満たす。λ=0\lambda=0ではy=0y=0がR\R上の解であるから、I(0,0,0)=RI(0,0,0)=\Rでありφ(t,0,0,0)=0\varphi(t,0,0,0)=0である。この解に沿う変分方程式の係数はA(t)=∂uf(t,0,0)=0A(t)=\partial_uf(t,0,0)=0、B(t)=∂λf(t,0,0)=1B(t)=\partial_\lambda f(t,0,0)=1である。Z(t)=tZ(t)=tはZ′=AZ+B=1Z'=A Z+B=1、Z(0)=0Z(0)=0を満たし、補題 3.5によりこの方程式の解は一意であるから、零解に沿う変分方程式の解でZ(0)=0Z(0)=0を満たすものはZ(t)=tZ(t)=tである。定理 3.8により∂λφ(t,0,0,0)=Z(t)=t\partial_\lambda\varphi(t,0,0,0)=Z(t)=tである。

(2)を示す。λ=κ2\lambda=\kappa^2(κ>0\kappa>0)とし、y(t)=κtan⁡(κt)y(t)=\kappa\tan(\kappa t)とおくと、∣t∣<π/(2κ)\lvert t\rvert<\pi/(2\kappa)でy′(t)=κ2(1+tan⁡2(κt))=λ+y(t)2y'(t)=\kappa^2(1+\tan^2(\kappa t))=\lambda+y(t)^2、y(0)=0y(0)=0である。∣t∣→π/(2κ)\lvert t\rvert\to\pi/(2\kappa)のとき∣y(t)∣→∞\lvert y(t)\rvert\to\inftyであるから、この解はそれ以上延長することができず、φ(t,0,0,λ)=y(t)\varphi(t,0,0,\lambda)=y(t)である。λ=−κ2\lambda=-\kappa^2(κ>0\kappa>0)とし、y(t)=−κtanh⁡(κt)y(t)=-\kappa\tanh(\kappa t)とおくと、y′(t)=−κ2(1−tanh⁡2(κt))=λ+y(t)2y'(t)=-\kappa^2(1-\tanh^2(\kappa t))=\lambda+y(t)^2、y(0)=0y(0)=0がR\R上で成り立つから、φ(t,0,0,λ)=y(t)\varphi(t,0,0,\lambda)=y(t)である。

t∈Rt\in\Rを固定する。t=0t=0ではφ(0,0,0,λ)=0\varphi(0,0,0,\lambda)=0である。t≠0t\ne0とする。0<λ<π2/(4t2)0<\lambda<\pi^2/(4t^2)ならば∣t∣<π/(2λ)\lvert t\rvert<\pi/(2\sqrt\lambda)であり、x=λ t→0x=\sqrt\lambda\,t\to0とすると

φ(t,0,0,λ)λ=tan⁡(λ t)λ=t⋅tan⁡xx⟶t\frac{\varphi(t,0,0,\lambda)}{\lambda}=\frac{\tan(\sqrt\lambda\,t)}{\sqrt\lambda}=t\cdot\frac{\tan x}{x}\longrightarrow t

である。λ<0\lambda<0では、x=−λ t→0x=\sqrt{-\lambda}\,t\to0とすると

φ(t,0,0,λ)λ=−−λtanh⁡(−λ t)λ=tanh⁡(−λ t)−λ=t⋅tanh⁡xx⟶t\frac{\varphi(t,0,0,\lambda)}{\lambda}=\frac{-\sqrt{-\lambda}\tanh(\sqrt{-\lambda}\,t)}{\lambda}=\frac{\tanh(\sqrt{-\lambda}\,t)}{\sqrt{-\lambda}}=t\cdot\frac{\tanh x}{x}\longrightarrow t

である。φ(t,0,0,0)=0\varphi(t,0,0,0)=0であるから、λ=0\lambda=0におけるλ↦φ(t,0,0,λ)\lambda\mapsto\varphi(t,0,0,\lambda)の微分係数はttであり、(1)と一致する。▨

問題 6.2.d,p∈N≥1d,p\in\NN、r≥2r\ge2とし、開集合Ω⊆R×Rd×Rp\Omega\subseteq\R\times\R^d\times\R^p上の写像f:Ω→Rdf:\Omega\to\R^dが、零集合NNと支配関数mKm_Kをもって状態とパラメータについてCrC^r級の Carathéodory 条件を満たすとする(t∈Nt\in Nでは偏導関数を00とする)。Rd×d\R^{d\times d}、Rd×p\R^{d\times p}をRd2\R^{d^2}、Rdp\R^{dp}と同一視してノルムを固定する。

Ω^={(t,u,U,V,λ):(t,u,λ)∈Ω}⊆R×Rd×Rd×d×Rd×p×Rp\widehat\Omega=\{(t,\boldsymbol u,U,V,\lambda):(t,\boldsymbol u,\lambda)\in\Omega\}\subseteq\R\times\R^d\times\R^{d\times d}\times\R^{d\times p}\times\R^p

上の写像

F(t,u,U,V,λ)=(f(t,u,λ), Duf(t,u,λ)U, Duf(t,u,λ)V+Dλf(t,u,λ))F(t,\boldsymbol u,U,V,\lambda)=\bigl(f(t,\boldsymbol u,\lambda),\ D_{\boldsymbol u}f(t,\boldsymbol u,\lambda)U,\ D_{\boldsymbol u}f(t,\boldsymbol u,\lambda)V+D_\lambda f(t,\boldsymbol u,\lambda)\bigr)

が、状態(u,U,V)(\boldsymbol u,U,V)、パラメータλ\lambdaについてCr−1C^{r-1}級の Carathéodory 条件を満たすことを示せ。

解答.

FFの成分は、ffの成分fif_iと、UUの成分UjkU_{jk}、VVの成分VjkV_{jk}を用いた

∑j=1d∂ujfi Ujk,∑j=1d∂ujfi Vjk+∂λkfi\sum_{j=1}^d\partial_{u_j}f_i\,U_{jk},\qquad\sum_{j=1}^d\partial_{u_j}f_i\,V_{jk}+\partial_{\lambda_k}f_i

である。

(u,U,V,λ)(\boldsymbol u,U,V,\lambda)を固定すると、{t:(t,u,U,V,λ)∈Ω^}={t:(t,u,λ)∈Ω}\{t:(t,\boldsymbol u,U,V,\lambda)\in\widehat\Omega\}=\{t:(t,\boldsymbol u,\lambda)\in\Omega\}であり、定義 3.2 条件 (a)と補題 3.3 (2)によりfif_i、∂ujfi\partial_{u_j}f_i、∂λkfi\partial_{\lambda_k}f_iはこの開集合上でttについて可測である。定数Ujk,VjkU_{jk},V_{jk}との積の有限和も可測である。

g(t,u,λ)g(t,\boldsymbol u,\lambda)をffのrr階以下の順序付き偏導関数の一つの成分とし、ggにUjkU_{jk}を掛けた関数gUjkgU_{jk}の(u,U,V,λ)(\boldsymbol u,U,V,\lambda)の座標についての順序付き偏導関数を考える。各段階の関数は、ggの(u,λ)(\boldsymbol u,\lambda)の座標についての順序付き偏導関数hhを用いてhUjkhU_{jk}、hh、00のいずれかの形であり、偏微分を一回行っても同じ形に留まる。実際、hUjkhU_{jk}を(u,λ)(\boldsymbol u,\lambda)の座標で偏微分すると(∂h)Ujk(\partial h)U_{jk}、UjkU_{jk}で偏微分するとhh、他のU,VU,Vの座標で偏微分すると00であり、hhをU,VU,Vの座標で偏微分すると00である。∂ujfi\partial_{u_j}f_iと∂λkfi\partial_{\lambda_k}f_iはt∉Nt\notin NでCr−1C^{r-1}級であるから、FFの成分のr−1r-1階以下の順序付き偏導関数は、ffのrr階以下の順序付き偏導関数の成分にUjkU_{jk}、VjkV_{jk}、11、00のいずれかを掛けたものの有限和であり、Ω^t\widehat\Omega_t上で存在して連続である。§E4.4 定理 1.7により、F(t,⋅)F(t,\cdot)はΩ^t\widehat\Omega_t上でCr−1C^{r-1}級である。

コンパクト集合K^⊆Ω^\widehat K\subseteq\widehat\Omegaを取り、その(t,u,λ)(t,\boldsymbol u,\lambda)成分への射影をKKとする。KKは連続写像によるコンパクト集合の像であるからコンパクトであり、Ω\Omegaに含まれる。M=max⁡K^max⁡j,k{1,∣Ujk∣,∣Vjk∣}M=\max_{\widehat K}\max_{j,k}\{1,\lvert U_{jk}\rvert,\lvert V_{jk}\rvert\}とおく。ノルムの同値性により、ベクトルのノルムと成分の絶対値の和を互いに定数倍で抑える定数が存在する。前段の形により、(t,u,U,V,λ)∈K^(t,\boldsymbol u,U,V,\lambda)\in\widehat K、t∉Nt\notin Nで、FFのr−1r-1階以下の各順序付き偏導関数のノルムは、ffのrr階以下の順序付き偏導関数のノルムの高々(d+1)(d+1)個の和にMMと定数を掛けたもので抑えられる。定義 3.2 条件 (c)によりそれらはmK(t)m_K(t)以下であるから、d,p,rd,p,r、ノルムとMMだけで決まる定数CCを用いてm^K^=C mK∈L1(R)\widehat m_{\widehat K}=C\,m_K\in L^1(\R)が支配関数となる。▨

参考文献

  1. Earl A. Coddington and Norman Levinson, Theory of Ordinary Differential Equations, McGraw Hill, 1955.Carathéodory 型の条件の下での解のパラメータについての連続性と微分可能性、およびパラメータについての解析性を参考にした。
  2. Philip Hartman, Ordinary Differential Equations, 2nd ed., Classics in Applied Mathematics, Society for Industrial and Applied Mathematics, Philadelphia, 2002.初期値とパラメータについての微分可能性と変分方程式を参考にした。

前提記事

11 本の記事・単元を表示