1 パラメータ付き Carathéodory 条件
定義 1.1. d , p ∈ N ≥ 1 d,p\in\NN d , p ∈ N ≥ 1 とし、R d \R^d R d とR p \R^p R p にノルム∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ を固定する。Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p を開集合とし、f : Ω → R d f:\Omega\to\R^d f : Ω → R d を写像とする。f f f がパラメータ付き Carathéodory 条件 (parametric Carathéodory conditions ) を局所的に満たすとは、長さが正の有限閉区間J J J と半径が正の閉球B ⊆ R d B\subseteq\R^d B ⊆ R d 、C ⊆ R p C\subseteq\R^p C ⊆ R p による閉矩形J × B × C J\times B\times C J × B × C でΩ \Omega Ω に含まれるものそれぞれに対して、Lebesgue 零集合N ⊆ J N\subseteq J N ⊆ J と非負関数m , ℓ ∈ L 1 ( J ) m,\ell\in L^1(J) m , ℓ ∈ L 1 ( J ) が存在し、次の条件が成り立つことをいう。
すべての( u , λ ) ∈ B × C (\boldsymbol u,\lambda)\in B\times C ( u , λ ) ∈ B × C に対して、t ↦ f ( t , u , λ ) t\mapsto f(t,\boldsymbol u,\lambda) t ↦ f ( t , u , λ ) はJ J J 上で Lebesgue 可測である。
すべてのt ∈ J ∖ N t\in J\setminus N t ∈ J ∖ N に対して、( u , λ ) ↦ f ( t , u , λ ) (\boldsymbol u,\lambda)\mapsto f(t,\boldsymbol u,\lambda) ( u , λ ) ↦ f ( t , u , λ ) はB × C B\times C B × C 上で連続である。
すべてのt ∈ J ∖ N t\in J\setminus N t ∈ J ∖ N 、u , v ∈ B \boldsymbol u,\boldsymbol v\in B u , v ∈ B 、λ ∈ C \lambda\in C λ ∈ 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 ∥ f ( t , u , λ )∥ ≤ m ( t ) , ∥ f ( t , u , λ ) − f ( t , v , λ )∥ ≤ ℓ ( t ) ∥ u − v ∥
が成り立つ。
パラメータ付き Carathéodory 条件は、§E13.4 定義 2.1 の Carathéodory 条件と異なり、時間について可積分な Lipschitz 係数ℓ \ell ℓ によるu \boldsymbol u u についての Lipschitz 条件を含む。
補題 1.2. n , k ∈ N ≥ 1 n,k\in\NN n , k ∈ N ≥ 1 とし、R n \R^n R n にノルム∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ を固定する。Ω ⊆ R × R n \Omega\subseteq\R\times\R^n Ω ⊆ R × R n を開集合とする。このとき、次が成り立つ。
任意のコンパクト集合L ⊆ Ω L\subseteq\Omega L ⊆ Ω に対して、あるε > 0 \varepsilon>0 ε > 0 が存在し、( t , v ) ∈ L (t,\boldsymbol v)\in L ( t , v ) ∈ L 、∣ s − t ∣ ≤ ε \lvert s-t\rvert\le\varepsilon ∣ s − t ∣ ≤ ε 、∥ w − v ∥ ≤ ε \lVert\boldsymbol w-\boldsymbol v\rVert\le\varepsilon ∥ w − v ∥ ≤ ε ならば( s , w ) ∈ Ω (s,\boldsymbol w)\in\Omega ( s , w ) ∈ Ω である。
写像g : Ω → R k g:\Omega\to\R^k g : Ω → R k が次を満たすとする。長さが正の有限閉区間J J J と半径が正の閉球B ⊆ R n B\subseteq\R^n B ⊆ R n についてJ × B ⊆ Ω J\times B\subseteq\Omega J × B ⊆ Ω ならば、各w ∈ B \boldsymbol w\in B w ∈ B でt ↦ g ( t , w ) t\mapsto g(t,\boldsymbol w) t ↦ g ( t , w ) はJ J J 上で Lebesgue 可測であり、ある Lebesgue 零集合N ⊆ J N\subseteq J N ⊆ J が存在し、各t ∈ J ∖ N t\in J\setminus N t ∈ J ∖ N でw ↦ g ( t , w ) \boldsymbol w\mapsto g(t,\boldsymbol w) w ↦ g ( t , w ) はB B B 上で連続である。このとき、a < b a<b a < b と、すべてのt ∈ [ a , b ] t\in[a,b] t ∈ [ a , b ] で( t , γ ( t ) ) ∈ Ω (t,\boldsymbol\gamma(t))\in\Omega ( t , γ ( t )) ∈ Ω を満たす連続写像γ : [ a , b ] → R n \boldsymbol\gamma:[a,b]\to\R^n γ : [ a , b ] → R n に対して、t ↦ g ( t , γ ( t ) ) t\mapsto g(t,\boldsymbol\gamma(t)) t ↦ g ( t , γ ( t )) は[ a , b ] [a,b] [ a , b ] 上で Lebesgue 可測である。
証明. 主張 (1) を示す。R × R n \R\times\R^n R × R n の点の距離d i s t ∞ \mathrm{dist}_\infty dist ∞ と閉ρ \rho ρ 近傍L ρ L_\rho L ρ を§E13.6 補題 3.1 のとおりに定める。§E13.6 補題 3.1 (1) をL L L に適用してρ L > 0 \rho_L>0 ρ L > 0 を取り、ε = ρ L \varepsilon=\rho_L ε = ρ L とおく。( t , v ) ∈ L (t,\boldsymbol v)\in L ( t , v ) ∈ L 、∣ s − t ∣ ≤ ε \lvert s-t\rvert\le\varepsilon ∣ s − t ∣ ≤ ε 、∥ w − v ∥ ≤ ε \lVert\boldsymbol w-\boldsymbol v\rVert\le\varepsilon ∥ w − v ∥ ≤ ε ならばd i s t ∞ ( ( s , w ) , L ) ≤ d i s t ∞ ( ( s , w ) , ( t , v ) ) ≤ ε \mathrm{dist}_\infty((s,\boldsymbol w),L)\le\mathrm{dist}_\infty((s,\boldsymbol w),(t,\boldsymbol v))\le\varepsilon dist ∞ (( s , w ) , L ) ≤ dist ∞ (( s , w ) , ( t , v )) ≤ ε であるから、( s , w ) ∈ L ε ⊆ Ω (s,\boldsymbol w)\in L_\varepsilon\subseteq\Omega ( s , w ) ∈ L ε ⊆ Ω である。
主張 (2) を示す。L = { ( t , γ ( t ) ) : a ≤ t ≤ b } L=\{(t,\boldsymbol\gamma(t)):a\le t\le b\} L = {( t , γ ( t )) : a ≤ t ≤ b } は連続写像によるコンパクト集合の像であるからコンパクトであり、Ω \Omega Ω に含まれる。主張 (1) のε \varepsilon ε を取る。γ \boldsymbol\gamma γ は一様連続であるから、∣ s − t ∣ ≤ δ \lvert s-t\rvert\le\delta ∣ s − t ∣ ≤ δ ならば∥ γ ( s ) − γ ( t ) ∥ ≤ ε \lVert\boldsymbol\gamma(s)-\boldsymbol\gamma(t)\rVert\le\varepsilon ∥ γ ( s ) − γ ( t )∥ ≤ ε となるδ > 0 \delta>0 δ > 0 が存在する。分点a = s 0 < s 1 < ⋯ < s j = b a=s_0<s_1<\cdots<s_j=b a = s 0 < s 1 < ⋯ < s j = b をs i − s i − 1 ≤ min { δ , ε } s_i-s_{i-1}\le\min\{\delta,\varepsilon\} s i − s i − 1 ≤ min { δ , ε } となるように取り、J i = [ s i − 1 , s i ] J_i=[s_{i-1},s_i] J i = [ s i − 1 , s i ] 、B i = { w : ∥ w − γ ( s i − 1 ) ∥ ≤ ε } B_i=\{\boldsymbol w:\lVert\boldsymbol w-\boldsymbol\gamma(s_{i-1})\rVert\le\varepsilon\} B i = { w : ∥ w − γ ( s i − 1 )∥ ≤ ε } とおく。( s i − 1 , γ ( s i − 1 ) ) ∈ L (s_{i-1},\boldsymbol\gamma(s_{i-1}))\in L ( s i − 1 , γ ( s i − 1 )) ∈ L であるからJ i × B i ⊆ Ω J_i\times B_i\subseteq\Omega J i × B i ⊆ Ω であり、γ ( J i ) ⊆ B i \boldsymbol\gamma(J_i)\subseteq B_i γ ( J i ) ⊆ B i である。g g g のJ i × B i J_i\times B_i J i × B i への制限は、d = n d=n d = n として§E13.4 補題 2.3 の仮定を満たし、γ ∣ J i \boldsymbol\gamma|_{J_i} γ ∣ J i はB i B_i B i に値をとる連続写像であるから可測である。同補題によりt ↦ g ( t , γ ( t ) ) t\mapsto g(t,\boldsymbol\gamma(t)) t ↦ g ( t , γ ( t )) は各J i J_i J i 上で可測であり、有限個のJ i J_i J i の和である[ a , b ] [a,b] [ a , b ] 上で可測である。▨
補題 1.3. d , p ∈ N ≥ 1 d,p\in\NN d , p ∈ N ≥ 1 とし、R d \R^d R d とR p \R^p R p にノルム∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ を固定する。開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d がパラメータ付き Carathéodory 条件を局所的に満たすとする。このとき、次が成り立つ。
各λ ∈ R p \lambda\in\R^p λ ∈ R p について、Ω λ = { ( t , u ) : ( t , u , λ ) ∈ Ω } \Omega^\lambda=\{(t,\boldsymbol u):(t,\boldsymbol u,\lambda)\in\Omega\} Ω λ = {( t , u ) : ( t , u , λ ) ∈ Ω } はR × R d \R\times\R^d R × R d の開集合であり、( t , u ) ↦ f ( t , u , λ ) (t,\boldsymbol u)\mapsto f(t,\boldsymbol u,\lambda) ( t , u ) ↦ f ( t , u , λ ) はΩ λ \Omega^\lambda Ω λ 上で Carathéodory 条件と時間について可積分な Lipschitz 係数の条件を局所的に満たす。
a < b a<b a < b とし、γ : [ a , b ] → R d × R p \boldsymbol\gamma:[a,b]\to\R^d\times\R^p γ : [ a , b ] → R d × R p を、すべてのt ∈ [ a , b ] t\in[a,b] t ∈ [ a , b ] で( t , γ ( t ) ) ∈ Ω (t,\boldsymbol\gamma(t))\in\Omega ( t , γ ( t )) ∈ Ω を満たす連続写像とする。このときt ↦ f ( t , γ ( t ) ) t\mapsto f(t,\boldsymbol\gamma(t)) t ↦ f ( t , γ ( t )) は[ a , b ] [a,b] [ a , b ] 上で Lebesgue 可測である。
a < b a<b a < b 、λ 0 ∈ R p \lambda_0\in\R^p λ 0 ∈ R p とし、y : [ a , b ] → R d \boldsymbol y:[a,b]\to\R^d y : [ a , b ] → R d を、すべてのt ∈ [ a , b ] t\in[a,b] t ∈ [ a , b ] で( t , y ( t ) , λ 0 ) ∈ Ω (t,\boldsymbol y(t),\lambda_0)\in\Omega ( t , y ( t ) , λ 0 ) ∈ Ω を満たす連続写像とする。このとき、ρ > 0 \rho>0 ρ > 0 、Lebesgue 零集合N ⊆ [ a , b ] N\subseteq[a,b] N ⊆ [ a , b ] 、非負関数m , ℓ ∈ L 1 ( [ a , b ] ) m,\ell\in L^1([a,b]) m , ℓ ∈ 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\} K ρ = {( t , u , λ ) : a ≤ t ≤ b , ∥ u − y ( t )∥ ≤ ρ , ∥ λ − λ 0 ∥ ≤ ρ }
はΩ \Omega Ω に含まれ、すべてのt ∈ [ a , b ] ∖ N t\in[a,b]\setminus N t ∈ [ a , b ] ∖ N について次が成り立つ。( u , λ ) ↦ f ( t , u , λ ) (\boldsymbol u,\lambda)\mapsto f(t,\boldsymbol u,\lambda) ( u , λ ) ↦ f ( t , u , λ ) は{ ( u , λ ) : ( t , u , λ ) ∈ K ρ } \{(\boldsymbol u,\lambda):(t,\boldsymbol u,\lambda)\in K_\rho\} {( u , λ ) : ( t , u , λ ) ∈ K ρ } 上で連続であり、( t , u , λ ) , ( t , v , λ ) ∈ K ρ (t,\boldsymbol u,\lambda),(t,\boldsymbol v,\lambda)\in K_\rho ( t , u , λ ) , ( t , v , λ ) ∈ K ρ ならば
∥ 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 ∥ f ( t , u , λ )∥ ≤ m ( t ) , ∥ f ( t , u , λ ) − f ( t , v , λ )∥ ≤ ℓ ( t ) ∥ u − v ∥
である。
証明. R d × R p \R^d\times\R^p R d × R p にノルム∥ ( u , λ ) ∥ = max { ∥ u ∥ , ∥ λ ∥ } \lVert(\boldsymbol u,\lambda)\rVert=\max\{\lVert\boldsymbol u\rVert,\lVert\lambda\rVert\} ∥( u , λ )∥ = max {∥ u ∥ , ∥ λ ∥} を入れる。このノルムの半径r > 0 r>0 r > 0 の閉球は、半径r r r の閉球B ⊆ R d B\subseteq\R^d B ⊆ R d とC ⊆ R p C\subseteq\R^p C ⊆ R p の直積B × C B\times C B × C である。
主張 (1) を示す。Ω λ \Omega^\lambda Ω λ は連続写像( t , u ) ↦ ( t , u , λ ) (t,\boldsymbol u)\mapsto(t,\boldsymbol u,\lambda) ( t , u ) ↦ ( t , u , λ ) によるΩ \Omega Ω の逆像であるから開集合である。Ω λ \Omega^\lambda Ω λ に含まれる閉矩形J × B J\times B J × B を取る。J × B × { λ } J\times B\times\{\lambda\} J × B × { λ } はコンパクトでΩ \Omega Ω に含まれるから、補題 1.2 (1) のε \varepsilon ε を取りC = { λ ′ : ∥ λ ′ − λ ∥ ≤ ε } C=\{\lambda':\lVert\lambda'-\lambda\rVert\le\varepsilon\} C = { λ ′ : ∥ λ ′ − λ ∥ ≤ ε } とおくと、J × B × C ⊆ Ω J\times B\times C\subseteq\Omega J × B × C ⊆ Ω である。この閉矩形に対する定義 1.1 のN , m , ℓ N,m,\ell N , m , ℓ を取り、λ \lambda λ を固定すると、( t , u ) ↦ f ( t , u , λ ) (t,\boldsymbol u)\mapsto f(t,\boldsymbol u,\lambda) ( t , u ) ↦ f ( t , u , λ ) はJ × B J\times B J × B 上で§E13.4 定義 2.1 の条件と Lipschitz 係数ℓ \ell ℓ の条件を満たす。
主張 (2) を示す。上のノルムの閉球はB × C B\times C B × C の形であるから、f f f はn = d + p n=d+p n = d + p 、k = d k=d k = d として補題 1.2 (2) の仮定を満たす。同補題から結論を得る。
主張 (3) を示す。n = d + p n=d+p n = d + p とし、R n = R d × R p \R^n=\R^d\times\R^p R n = R d × R p に上のノルムを入れ、距離d i s t ∞ \mathrm{dist}_\infty dist ∞ と閉ρ \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\} L = {( t , y ( t ) , λ 0 ) : a ≤ t ≤ b } は連続写像によるコンパクト集合の像であるからコンパクトであり、Ω \Omega Ω に含まれる。§E13.6 補題 3.1 (1) をL L L に適用してρ = ρ L > 0 \rho=\rho_L>0 ρ = ρ L > 0 を取ると、L ρ ⊆ Ω L_\rho\subseteq\Omega L ρ ⊆ Ω である。( t , u , λ ) ∈ K ρ (t,\boldsymbol u,\lambda)\in K_\rho ( t , u , λ ) ∈ K ρ ならばd i s t ∞ ( ( 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 dist ∞ (( t , u , λ ) , ( t , y ( t ) , λ 0 )) = max {∥ u − y ( t )∥ , ∥ λ − λ 0 ∥} ≤ ρ であるから、K ρ ⊆ L ρ ⊆ Ω K_\rho\subseteq L_\rho\subseteq\Omega K ρ ⊆ L ρ ⊆ Ω である。
上のノルムの半径r > 0 r>0 r > 0 の閉球はB × C B\times C B × C の形であるから、定義 1.1 条件 (b) と定義 1.1 条件 (c) により、f f f はk = d k=d k = d として§E13.6 補題 3.1 (2) の仮定と§E13.6 補題 3.1 (3) の仮定を満たし、閉矩形ごとの零集合と非負関数はパラメータ付き Carathéodory 条件のN , m , ℓ N,m,\ell N , m , ℓ である。§E13.6 補題 3.1 (3) をL L L とρ \rho ρ に適用して、零集合N L ⊆ R N_L\subseteq\R N L ⊆ R と非負関数m L , ℓ L ∈ L 1 ( R ) m_L,\ell_L\in L^1(\R) m L , ℓ L ∈ L 1 ( R ) を取る。N = N L ∩ [ a , b ] N=N_L\cap[a,b] N = N L ∩ [ a , b ] 、m = m L ∣ [ a , b ] m=m_L|_{[a,b]} m = m L ∣ [ a , b ] 、ℓ = ℓ L ∣ [ a , b ] \ell=\ell_L|_{[a,b]} ℓ = ℓ L ∣ [ a , b ] とおく。t ∈ [ a , b ] ∖ N t\in[a,b]\setminus N t ∈ [ a , b ] ∖ N では、K ρ K_\rho K ρ の時刻t t t における切り口はL ρ L_\rho L ρ の時刻t t t における切り口に含まれるから、同補題の結論により、( u , λ ) ↦ f ( t , u , λ ) (\boldsymbol u,\lambda)\mapsto f(t,\boldsymbol u,\lambda) ( u , λ ) ↦ f ( t , u , λ ) はK ρ K_\rho K ρ の切り口上で連続であり、( t , u , λ ) , ( t , v , λ ) ∈ K ρ (t,\boldsymbol u,\lambda),(t,\boldsymbol v,\lambda)\in K_\rho ( t , u , λ ) , ( t , v , λ ) ∈ K ρ ならば∥ f ( t , u , λ ) ∥ ≤ m ( t ) \lVert f(t,\boldsymbol u,\lambda)\rVert\le m(t) ∥ f ( t , u , λ )∥ ≤ 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 ∥ f ( t , u , λ ) − f ( t , v , λ )∥ ≤ ℓ ( t ) ∥ u − v ∥ である。▨
2 解写像の連続性
定義 2.1. d , p ∈ N ≥ 1 d,p\in\NN d , p ∈ N ≥ 1 とし、開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d がパラメータ付き Carathéodory 条件を局所的に満たすとする。( t 0 , η , λ ) ∈ Ω (t_0,\boldsymbol\eta,\lambda)\in\Omega ( t 0 , η , λ ) ∈ Ω に対して、Ω λ = { ( t , u ) : ( t , u , λ ) ∈ Ω } \Omega^\lambda=\{(t,\boldsymbol u):(t,\boldsymbol u,\lambda)\in\Omega\} Ω λ = {( t , u ) : ( t , u , λ ) ∈ Ω } 上の初期値問題
y ′ = f ( t , y , λ ) , y ( t 0 ) = η \boldsymbol y'=f(t,\boldsymbol y,\lambda),\qquad \boldsymbol y(t_0)=\boldsymbol\eta y ′ = f ( t , y , λ ) , y ( t 0 ) = η は一意な最大絶対連続解をもつ(補題 1.3 (1) 、§E13.6 定理 2.1 )。その最大存在区間をI ( t 0 , η , λ ) I(t_0,\boldsymbol\eta,\lambda) I ( t 0 , η , λ ) と書き、
W = { ( t , t 0 , η , λ ) ∈ R × R × R d × R p : ( t 0 , η , λ ) ∈ Ω , t ∈ I ( t 0 , η , λ ) } 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)\} W = {( t , t 0 , η , λ ) ∈ R × R × R d × R p : ( t 0 , η , λ ) ∈ Ω , t ∈ I ( t 0 , η , λ )} とおく。( t , t 0 , η , λ ) ∈ W (t,t_0,\boldsymbol\eta,\lambda)\in W ( t , t 0 , η , λ ) ∈ W に、この最大絶対連続解の時刻t t t における値φ ( t , t 0 , η , λ ) \varphi(t,t_0,\boldsymbol\eta,\lambda) φ ( t , t 0 , η , λ ) を対応させる写像φ : W → R d \varphi:W\to\R^d φ : W → R d を、f f f の解写像 (solution map ) という。
定理 2.2. d , p ∈ N ≥ 1 d,p\in\NN d , p ∈ N ≥ 1 とし、R d \R^d R d とR p \R^p R p にノルム∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ を固定する。開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d がパラメータ付き Carathéodory 条件を局所的に満たすとし、φ : W → R d \varphi:W\to\R^d φ : W → R d をその解写像、I ( t 0 , η , λ ) I(t_0,\boldsymbol\eta,\lambda) I ( t 0 , η , λ ) を最大存在区間とする。このとき、次が成り立つ。
( t 0 , η , λ ) ∈ Ω (t_0,\boldsymbol\eta,\lambda)\in\Omega ( t 0 , η , λ ) ∈ Ω と、a < t 0 < b a<t_0<b a < t 0 < b かつ[ a , b ] ⊆ I ( t 0 , η , λ ) [a,b]\subseteq I(t_0,\boldsymbol\eta,\lambda) [ a , b ] ⊆ I ( t 0 , η , λ ) を満たすa , b a,b a , b を取る。このとき( t 0 , η , λ ) (t_0,\boldsymbol\eta,\lambda) ( t 0 , η , λ ) の開近傍V ⊆ Ω V\subseteq\Omega V ⊆ Ω が存在し、すべての( t 0 ′ , η ′ , λ ′ ) ∈ V (t_0',\boldsymbol\eta',\lambda')\in V ( t 0 ′ , η ′ , λ ′ ) ∈ V について[ a , b ] ⊆ I ( t 0 ′ , η ′ , λ ′ ) [a,b]\subseteq I(t_0',\boldsymbol\eta',\lambda') [ a , b ] ⊆ I ( t 0 ′ , η ′ , λ ′ ) であり、( t 0 ′ , η ′ , λ ′ ) → ( t 0 , η , λ ) (t_0',\boldsymbol\eta',\lambda')\to(t_0,\boldsymbol\eta,\lambda) ( t 0 ′ , η ′ , λ ′ ) → ( t 0 , η , λ ) のとき
max a ≤ t ≤ b ∥ φ ( t , t 0 ′ , η ′ , λ ′ ) − φ ( t , t 0 , η , λ ) ∥ ⟶ 0 \max_{a\le t\le b}\lVert\varphi(t,t_0',\boldsymbol\eta',\lambda')-\varphi(t,t_0,\boldsymbol\eta,\lambda)\rVert\longrightarrow0 a ≤ t ≤ b max ∥ φ ( t , t 0 ′ , η ′ , λ ′ ) − φ ( t , t 0 , η , λ )∥ ⟶ 0
である。
W W W はR × R × R d × R p \R\times\R\times\R^d\times\R^p R × R × R d × R p の開集合であり、φ \varphi φ はW W W 上で連続である。
証明. 主張 (1) を示す。y = φ ( ⋅ , t 0 , η , λ ) ∣ [ a , b ] \boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda)|_{[a,b]} y = φ ( ⋅ , t 0 , η , λ ) ∣ [ a , b ] とおく。補題 1.3 (3) をy \boldsymbol y y とλ 0 = λ \lambda_0=\lambda λ 0 = λ に適用し、ρ > 0 \rho>0 ρ > 0 、零集合N N N 、非負関数m , ℓ ∈ L 1 ( [ a , b ] ) m,\ell\in L^1([a,b]) m , ℓ ∈ L 1 ([ a , b ]) とK ρ ⊆ Ω K_\rho\subseteq\Omega K ρ ⊆ Ω を取る。∥ λ ′ − λ ∥ ≤ ρ \lVert\lambda'-\lambda\rVert\le\rho ∥ λ ′ − λ ∥ ≤ ρ を満たすλ ′ \lambda' λ ′ に対して、t ∈ [ a , b ] ∖ N t\in[a,b]\setminus N t ∈ [ a , b ] ∖ 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 r λ ′ ( t ) = ∥ f ( t , y ( t ) , λ ′ ) − f ( t , y ( t ) , λ )∥ とし、t ∈ N t\in N t ∈ N ではr λ ′ ( t ) = 0 r_{\lambda'}(t)=0 r λ ′ ( t ) = 0 とおく。補題 1.3 (2) によりr λ ′ r_{\lambda'} r λ ′ は可測であり、r λ ′ ≤ 2 m r_{\lambda'}\le2m r λ ′ ≤ 2 m である。t ∈ [ a , b ] ∖ N t\in[a,b]\setminus N t ∈ [ a , b ] ∖ N ではK ρ K_\rho K ρ の切り口上の連続性から、λ ′ → λ \lambda'\to\lambda λ ′ → λ のときr λ ′ ( t ) → 0 r_{\lambda'}(t)\to0 r λ ′ ( t ) → 0 である。λ \lambda λ へ収束する任意の点列に§E10.11 定理 3.1 を適用すると、λ ′ → λ \lambda'\to\lambda λ ′ → λ のとき∫ a b r λ ′ → 0 \int_a^br_{\lambda'}\to0 ∫ a b r λ ′ → 0 である。
E = exp ( ∫ a b ℓ ) E=\exp(\int_a^b\ell) E = exp ( ∫ a b ℓ ) とおき、( t 0 ′ , η ′ , λ ′ ) ∈ Ω (t_0',\boldsymbol\eta',\lambda')\in\Omega ( t 0 ′ , η ′ , λ ′ ) ∈ Ω に対して
Q ( t 0 ′ , η ′ , λ ′ ) = E ( ∥ η ′ − y ( t 0 ′ ) ∥ + ∫ a b r λ ′ ( s ) d s ) Q(t_0',\boldsymbol\eta',\lambda')=E\left(\lVert\boldsymbol\eta'-\boldsymbol y(t_0')\rVert+\int_a^br_{\lambda'}(s)\,ds\right) Q ( t 0 ′ , η ′ , λ ′ ) = E ( ∥ η ′ − y ( t 0 ′ )∥ + ∫ a b r λ ′ ( s ) d s ) とおく(a ≤ t 0 ′ ≤ b a\le t_0'\le b a ≤ t 0 ′ ≤ b 、∥ λ ′ − λ ∥ ≤ ρ \lVert\lambda'-\lambda\rVert\le\rho ∥ λ ′ − λ ∥ ≤ ρ の範囲で)。∥ η ′ − y ( t 0 ′ ) ∥ ≤ ∥ η ′ − η ∥ + ∥ y ( t 0 ) − y ( t 0 ′ ) ∥ \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 ( t 0 ′ )∥ ≤ ∥ η ′ − η ∥ + ∥ y ( t 0 ) − y ( t 0 ′ )∥ とy \boldsymbol y y の連続性と前段により、あるε ∈ ( 0 , ρ ] \varepsilon\in(0,\rho] ε ∈ ( 0 , ρ ] が存在し、V = { ( t 0 ′ , η ′ , λ ′ ) : ∣ t 0 ′ − t 0 ∣ < ε , ∥ η ′ − η ∥ < ε , ∥ λ ′ − λ ∥ < ε } 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\} V = {( t 0 ′ , η ′ , λ ′ ) : ∣ t 0 ′ − t 0 ∣ < ε , ∥ η ′ − η ∥ < ε , ∥ λ ′ − λ ∥ < ε } の各点でa < t 0 ′ < b a<t_0'<b a < t 0 ′ < b かつQ ( t 0 ′ , η ′ , λ ′ ) < ρ Q(t_0',\boldsymbol\eta',\lambda')<\rho Q ( t 0 ′ , η ′ , λ ′ ) < ρ となり、( t 0 ′ , η ′ , λ ′ ) → ( t 0 , η , λ ) (t_0',\boldsymbol\eta',\lambda')\to(t_0,\boldsymbol\eta,\lambda) ( t 0 ′ , η ′ , λ ′ ) → ( t 0 , η , λ ) のときQ ( t 0 ′ , η ′ , λ ′ ) → 0 Q(t_0',\boldsymbol\eta',\lambda')\to0 Q ( t 0 ′ , η ′ , λ ′ ) → 0 である。( t 0 ′ , η ′ , λ ′ ) ∈ V (t_0',\boldsymbol\eta',\lambda')\in V ( t 0 ′ , η ′ , λ ′ ) ∈ V ならば∥ η ′ − y ( t 0 ′ ) ∥ < ρ \lVert\boldsymbol\eta'-\boldsymbol y(t_0')\rVert<\rho ∥ η ′ − y ( t 0 ′ )∥ < ρ であるから、( t 0 ′ , η ′ , λ ′ ) ∈ K ρ ⊆ Ω (t_0',\boldsymbol\eta',\lambda')\in K_\rho\subseteq\Omega ( t 0 ′ , η ′ , λ ′ ) ∈ K ρ ⊆ Ω である。
( t 0 ′ , η ′ , λ ′ ) ∈ V (t_0',\boldsymbol\eta',\lambda')\in V ( t 0 ′ , η ′ , λ ′ ) ∈ V を固定する。Ω ′ = Ω λ ∩ Ω λ ′ \Omega'=\Omega^\lambda\cap\Omega^{\lambda'} Ω ′ = Ω λ ∩ Ω λ ′ は開集合であり、補題 1.3 (1) によりf ( ⋅ , ⋅ , λ ′ ) f(\cdot,\cdot,\lambda') f ( ⋅ , ⋅ , λ ′ ) はΩ ′ \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 ′ = {( t , u ) : a ≤ t ≤ b , ∥ u − y ( t )∥ ≤ ρ } は、K ρ ⊆ Ω K_\rho\subseteq\Omega K ρ ⊆ Ω によりΩ ′ \Omega' Ω ′ に含まれるコンパクト集合である。t ∈ [ a , b ] ∖ N t\in[a,b]\setminus N t ∈ [ a , b ] ∖ N と( t , u ) , ( t , v ) ∈ K ′ (t,\boldsymbol u),(t,\boldsymbol v)\in K' ( t , u ) , ( t , v ) ∈ 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 ∥ f ( t , u , λ ′ ) − f ( t , v , λ ′ )∥ ≤ ℓ ( t ) ∥ u − v ∥ である。§E13.7 定理 3.1 をΩ ′ \Omega' Ω ′ 上の写像f ( ⋅ , ⋅ , λ ) f(\cdot,\cdot,\lambda) f ( ⋅ , ⋅ , λ ) とg = f ( ⋅ , ⋅ , λ ′ ) g=f(\cdot,\cdot,\lambda') g = f ( ⋅ , ⋅ , λ ′ ) 、K = K ′ K=K' K = K ′ 、r = r λ ′ r=r_{\lambda'} r = r λ ′ 、t ∗ = t 0 ′ t_*=t_0' t ∗ = t 0 ′ 、初期値η ′ \boldsymbol\eta' η ′ に適用すると、z ′ = f ( t , z , λ ′ ) \boldsymbol z'=f(t,\boldsymbol z,\lambda') z ′ = f ( t , z , λ ′ ) 、z ( t 0 ′ ) = η ′ \boldsymbol z(t_0')=\boldsymbol\eta' z ( t 0 ′ ) = η ′ の[ a , b ] [a,b] [ a , b ] 上の絶対連続解z \boldsymbol z z が存在し、max [ a , b ] ∥ z − y ∥ ≤ Q ( t 0 ′ , η ′ , λ ′ ) < ρ \max_{[a,b]}\lVert\boldsymbol z-\boldsymbol y\rVert\le Q(t_0',\boldsymbol\eta',\lambda')<\rho max [ a , b ] ∥ z − y ∥ ≤ Q ( t 0 ′ , η ′ , λ ′ ) < ρ である。
( α ′ , β ′ ) = I ( t 0 ′ , η ′ , λ ′ ) (\alpha',\beta')=I(t_0',\boldsymbol\eta',\lambda') ( α ′ , β ′ ) = I ( t 0 ′ , η ′ , λ ′ ) とし、ψ = φ ( ⋅ , t 0 ′ , η ′ , λ ′ ) \boldsymbol\psi=\varphi(\cdot,t_0',\boldsymbol\eta',\lambda') ψ = φ ( ⋅ , t 0 ′ , η ′ , λ ′ ) とする。z ∣ ( a , b ) \boldsymbol z|_{(a,b)} z ∣ ( a , b ) はt 0 ′ t_0' t 0 ′ を含む開区間上の絶対連続解であるから、Ω λ ′ \Omega^{\lambda'} Ω λ ′ 上で§E13.6 定理 2.1 を適用すると( a , b ) ⊆ ( α ′ , β ′ ) (a,b)\subseteq(\alpha',\beta') ( a , b ) ⊆ ( α ′ , β ′ ) かつ( a , b ) (a,b) ( a , b ) 上でψ = z \boldsymbol\psi=\boldsymbol z ψ = z である。β ′ = b \beta'=b β ′ = b ならば、[ t 0 ′ , b ) [t_0',b) [ t 0 ′ , b ) 上のψ \boldsymbol\psi ψ のグラフはコンパクト集合K ′ ⊆ Ω λ ′ K'\subseteq\Omega^{\lambda'} K ′ ⊆ Ω λ ′ に含まれる。これは§E13.6 系 4.3 に反する。従ってβ ′ > b \beta'>b β ′ > b であり、同様にα ′ < a \alpha'<a α ′ < a である。[ a , b ] ⊆ I ( t 0 ′ , η ′ , λ ′ ) [a,b]\subseteq I(t_0',\boldsymbol\eta',\lambda') [ a , b ] ⊆ I ( t 0 ′ , η ′ , λ ′ ) であり、連続性により[ a , b ] [a,b] [ a , b ] 上でψ = z \boldsymbol\psi=\boldsymbol z ψ = z であるから、
max a ≤ t ≤ b ∥ φ ( t , t 0 ′ , η ′ , λ ′ ) − φ ( t , t 0 , η , λ ) ∥ ≤ Q ( t 0 ′ , η ′ , λ ′ ) \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') a ≤ t ≤ b max ∥ φ ( t , t 0 ′ , η ′ , λ ′ ) − φ ( t , t 0 , η , λ )∥ ≤ Q ( t 0 ′ , η ′ , λ ′ ) であり、右辺は( t 0 ′ , η ′ , λ ′ ) → ( t 0 , η , λ ) (t_0',\boldsymbol\eta',\lambda')\to(t_0,\boldsymbol\eta,\lambda) ( t 0 ′ , η ′ , λ ′ ) → ( t 0 , η , λ ) のとき零へ収束する。
主張 (2) を示す。( t , t 0 , η , λ ) ∈ W (t,t_0,\boldsymbol\eta,\lambda)\in W ( t , t 0 , η , λ ) ∈ W を取る。I ( t 0 , η , λ ) I(t_0,\boldsymbol\eta,\lambda) I ( t 0 , η , λ ) はt , t 0 t,t_0 t , t 0 を含む開区間であるから、a < min { t , t 0 } a<\min\{t,t_0\} a < min { t , t 0 } 、max { t , t 0 } < b \max\{t,t_0\}<b max { t , t 0 } < b 、[ a , b ] ⊆ I ( t 0 , η , λ ) [a,b]\subseteq I(t_0,\boldsymbol\eta,\lambda) [ a , b ] ⊆ I ( t 0 , η , λ ) を満たすa , b a,b a , b が存在する。主張 (1) のV V V を取ると、開集合( a , b ) × V (a,b)\times V ( a , b ) × V は( t , t 0 , η , λ ) (t,t_0,\boldsymbol\eta,\lambda) ( t , t 0 , η , λ ) を含みW W W に含まれる。( t ′ , t 0 ′ , η ′ , λ ′ ) ∈ ( a , b ) × V (t',t_0',\boldsymbol\eta',\lambda')\in(a,b)\times V ( t ′ , t 0 ′ , η ′ , λ ′ ) ∈ ( a , b ) × V について
∥ φ ( t ′ , t 0 ′ , η ′ , λ ′ ) − φ ( t , t 0 , η , λ ) ∥ ≤ max a ≤ s ≤ b ∥ φ ( s , t 0 ′ , η ′ , λ ′ ) − φ ( s , t 0 , η , λ ) ∥ + ∥ φ ( t ′ , t 0 , η , λ ) − φ ( t , t 0 , η , λ ) ∥ \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 ∥ φ ( t ′ , t 0 ′ , η ′ , λ ′ ) − φ ( t , t 0 , η , λ )∥ ≤ a ≤ s ≤ b max ∥ φ ( s , t 0 ′ , η ′ , λ ′ ) − φ ( s , t 0 , η , λ )∥ + ∥ φ ( t ′ , t 0 , η , λ ) − φ ( t , t 0 , η , λ )∥ であり、右辺の第一項は主張 (1) により、第二項はφ ( ⋅ , t 0 , η , λ ) \varphi(\cdot,t_0,\boldsymbol\eta,\lambda) φ ( ⋅ , t 0 , η , λ ) の連続性により、( t ′ , t 0 ′ , η ′ , λ ′ ) → ( t , t 0 , η , λ ) (t',t_0',\boldsymbol\eta',\lambda')\to(t,t_0,\boldsymbol\eta,\lambda) ( t ′ , t 0 ′ , η ′ , λ ′ ) → ( t , t 0 , η , λ ) のとき零へ収束する。▨
3 変分方程式と一階の微分可能性
補題 3.1. K ∈ { R , C } \mathbb K\in\{\R,\C\} K ∈ { R , C } 、k ∈ N ≥ 1 k\in\NN k ∈ N ≥ 1 とし、K k \mathbb K^k K k にノルム∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ を固定し、k × k k\times k k × k 行列にはその作用素ノルムを用いる。
a < b a<b a < b 、J = [ a , b ] J=[a,b] J = [ a , b ] 、A ∈ L 1 ( J ; K k × k ) A\in L^1(J;\mathbb K^{k\times k}) A ∈ L 1 ( J ; K k × k ) 、b ∈ L 1 ( J ; K k ) \boldsymbol b\in L^1(J;\mathbb K^k) b ∈ L 1 ( J ; K k ) 、t ∗ ∈ J t_*\in J t ∗ ∈ J 、ξ ∈ K k \boldsymbol\xi\in\mathbb K^k ξ ∈ K k とする。このとき、J J J のほとんど至る所でx ′ = A ( t ) x + b ( t ) \boldsymbol x'=A(t)\boldsymbol x+\boldsymbol b(t) x ′ = A ( t ) x + b ( t ) を満たし、x ( t ∗ ) = ξ \boldsymbol x(t_*)=\boldsymbol\xi x ( t ∗ ) = ξ を満たす絶対連続関数x : J → K k \boldsymbol x:J\to\mathbb K^k x : J → K k が一意に存在する。
主張 (1) のJ , A , b , t ∗ , ξ , x J,A,\boldsymbol b,t_*,\boldsymbol\xi,\boldsymbol x J , A , b , t ∗ , ξ , x と、A n ∈ L 1 ( J ; K k × k ) A_n\in L^1(J;\mathbb K^{k\times k}) A n ∈ L 1 ( J ; K k × k ) 、b n ∈ L 1 ( J ; K k ) \boldsymbol b_n\in L^1(J;\mathbb K^k) b n ∈ L 1 ( J ; K k ) 、t n ∈ J t_n\in J t n ∈ J 、ξ n ∈ K k \boldsymbol\xi_n\in\mathbb K^k ξ n ∈ K k (n ∈ N ≥ 1 n\in\NN n ∈ N ≥ 1 )が
∫ J ∥ A n − A ∥ ⟶ 0 , ∫ J ∥ b n − b ∥ ⟶ 0 , t n ⟶ 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 ∫ J ∥ A n − A ∥ ⟶ 0 , ∫ J ∥ b n − b ∥ ⟶ 0 , t n ⟶ t ∗ , ξ n ⟶ ξ
を満たすとし、x n \boldsymbol x_n x n をA n , b n , t n , ξ n A_n,\boldsymbol b_n,t_n,\boldsymbol\xi_n A n , b n , t n , ξ n に対する主張 (1) の解とする。このときmax t ∈ J ∥ x n ( t ) − x ( t ) ∥ → 0 \max_{t\in J}\lVert\boldsymbol x_n(t)-\boldsymbol x(t)\rVert\to0 max t ∈ J ∥ x n ( t ) − x ( t )∥ → 0 である。
I ⊆ R I\subseteq\R I ⊆ R を内部が空でない区間、A ∈ L l o c 1 ( I ; K k × k ) A\in L^1_{\mathrm{loc}}(I;\mathbb K^{k\times k}) A ∈ L loc 1 ( I ; K k × k ) 、b ∈ L l o c 1 ( I ; K k ) \boldsymbol b\in L^1_{\mathrm{loc}}(I;\mathbb K^k) b ∈ L loc 1 ( I ; K k ) 、t ∗ ∈ I t_*\in I t ∗ ∈ I 、ξ ∈ K k \boldsymbol\xi\in\mathbb K^k ξ ∈ K k とする。このとき、I I I のほとんど至る所でx ′ = A ( t ) x + b ( t ) \boldsymbol x'=A(t)\boldsymbol x+\boldsymbol b(t) x ′ = A ( t ) x + b ( t ) を満たし、x ( t ∗ ) = ξ \boldsymbol x(t_*)=\boldsymbol\xi x ( t ∗ ) = ξ を満たすx ∈ A C l o c ( I ; K k ) \boldsymbol x\in AC_{\mathrm{loc}}(I;\mathbb K^k) x ∈ A C loc ( I ; K k ) が一意に存在する。
q ∈ N ≥ 1 q\in\NN q ∈ N ≥ 1 とし、K k × q \mathbb K^{k\times q} K k × q に値をとる未知関数X X X と、L 1 ( J ; K k × q ) L^1(J;\mathbb K^{k\times q}) L 1 ( J ; K k × q ) またはL l o c 1 ( I ; K k × q ) L^1_{\mathrm{loc}}(I;\mathbb K^{k\times q}) L loc 1 ( I ; K k × q ) に属するB B B に対する方程式X ′ = A ( t ) X + B ( t ) X'=A(t)X+B(t) X ′ = A ( t ) X + B ( t ) についても、主張 (1) 、主張 (2) 、主張 (3) の結論が成り立つ。
証明. 主張 (3) を示す。§E13.7 系 4.2 をI I I 、A A A 、b \boldsymbol b b 、t ∗ t_* t ∗ 、η = ξ \boldsymbol\eta=\boldsymbol\xi η = ξ に適用すると、結論を得る。
主張 (1) を示す。L 1 ( J ) L^1(J) L 1 ( J ) の関数はJ J J 上で局所可積分であるから、主張 (3) をI = J I=J I = J に適用すると、J J J のほとんど至る所で方程式を満たしx ( t ∗ ) = ξ \boldsymbol x(t_*)=\boldsymbol\xi x ( t ∗ ) = ξ を満たすx ∈ A C l o c ( J ; K k ) \boldsymbol x\in AC_{\mathrm{loc}}(J;\mathbb K^k) x ∈ A C loc ( J ; K k ) が一意に存在する。J J J は有限閉区間であるから、A C l o c ( J ; K k ) AC_{\mathrm{loc}}(J;\mathbb K^k) A C loc ( J ; K k ) はJ J J 上の絶対連続関数の全体である。
主張 (2) を示す。K = C \mathbb K=\C K = C のときは、C k \C^k C k を実部と虚部によりR 2 k \R^{2k} R 2 k と同一視する。固定したノルムはR 2 k \R^{2k} R 2 k のノルムであり、各t t t についてu ↦ A ( t ) u \boldsymbol u\mapsto A(t)\boldsymbol u u ↦ A ( t ) u は実線形であり、絶対連続性は実座標成分ごとに判定される(§E13.2 補題 1.2 (3) )。n ∈ N ≥ 1 n\in\NN n ∈ N ≥ 1 を固定し、§E13.7 定理 2.1 をS = J × K k S=J\times\mathbb K^k S = J × K k 、二つの写像( t , u ) ↦ A ( t ) u + b ( t ) (t,\boldsymbol u)\mapsto A(t)\boldsymbol u+\boldsymbol b(t) ( t , u ) ↦ A ( t ) u + b ( t ) と( t , u ) ↦ A n ( t ) u + b n ( t ) (t,\boldsymbol u)\mapsto A_n(t)\boldsymbol u+\boldsymbol b_n(t) ( t , u ) ↦ A n ( t ) u + b n ( t ) 、二つの解x , x n \boldsymbol x,\boldsymbol x_n x , x n 、時刻t n t_n t n 、
ℓ n = ∥ A n ∥ , r n = ∥ A n − A ∥ max J ∥ x ∥ + ∥ b n − 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 ℓ n = ∥ A n ∥ , r n = ∥ A n − A ∥ J max ∥ x ∥ + ∥ b n − b ∥ に適用すると、
max t ∈ J ∥ x n ( t ) − x ( t ) ∥ ≤ exp ( ∫ J ∥ A n ∥ ) ( ∥ ξ n − x ( t n ) ∥ + ∫ J r n ) \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) t ∈ J max ∥ x n ( t ) − x ( t )∥ ≤ exp ( ∫ J ∥ A n ∥ ) ( ∥ ξ n − x ( t n )∥ + ∫ J r n ) を得る。∫ J ∥ A n ∥ ≤ ∫ J ∥ A ∥ + ∫ J ∥ A n − A ∥ \int_J\lVert A_n\rVert\le\int_J\lVert A\rVert+\int_J\lVert A_n-A\rVert ∫ J ∥ A n ∥ ≤ ∫ J ∥ A ∥ + ∫ J ∥ A n − A ∥ は有界であり、∥ ξ n − x ( t n ) ∥ ≤ ∥ ξ n − ξ ∥ + ∥ x ( t ∗ ) − x ( t n ) ∥ → 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 ∥ ξ n − x ( t n )∥ ≤ ∥ ξ n − ξ ∥ + ∥ x ( t ∗ ) − x ( t n )∥ → 0 、∫ J r n → 0 \int_Jr_n\to0 ∫ J r n → 0 であるから、右辺は零へ収束する。
行列値の方程式X ′ = A ( t ) X + B ( t ) X'=A(t)X+B(t) X ′ = A ( t ) X + B ( t ) は、X X X とB B B の各列についての方程式x j ′ = A ( t ) x j + B ( t ) e j \boldsymbol x_j'=A(t)\boldsymbol x_j+B(t)\boldsymbol e_j x j ′ = A ( t ) x j + B ( t ) e j の組であるから、各列に以上の結論を適用すればよい。▨
定義 3.2. d , p , r ∈ N ≥ 1 d,p,r\in\NN d , p , r ∈ N ≥ 1 とし、R d \R^d R d とR p \R^p R p にノルム∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ を固定する。Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p を開集合とし、f : Ω → R d f:\Omega\to\R^d f : Ω → R d を写像とする。w = ( u , λ ) ∈ R d × R p \boldsymbol w=(\boldsymbol u,\lambda)\in\R^d\times\R^p w = ( u , λ ) ∈ R d × R p の座標をw 1 , … , w d + p w_1,\ldots,w_{d+p} w 1 , … , w d + p と書き、w i w_i w i についての偏微分を∂ i \partial_i ∂ i と書く。t ∈ R t\in\R t ∈ R に対してΩ t = { w : ( t , w ) ∈ Ω } \Omega_t=\{\boldsymbol w:(t,\boldsymbol w)\in\Omega\} Ω t = { w : ( t , w ) ∈ Ω } とおく。f f f が状態とパラメータについてC r C^r C r 級の Carathéodory 条件 (Carathéodory conditions of class C^r in the state and the parameter ) を満たすとは、次の条件が成り立つことをいう。
各w ∈ R d × R p \boldsymbol w\in\R^d\times\R^p w ∈ R d × R p について、t ↦ f ( t , w ) t\mapsto f(t,\boldsymbol w) t ↦ f ( t , w ) は開集合{ t : ( t , w ) ∈ Ω } \{t:(t,\boldsymbol w)\in\Omega\} { t : ( t , w ) ∈ Ω } 上で Lebesgue 可測である。
ある Lebesgue 零集合N ⊆ R N\subseteq\R N ⊆ R が存在し、すべてのt ∈ R ∖ N t\in\R\setminus N t ∈ R ∖ N についてw ↦ f ( t , w ) \boldsymbol w\mapsto f(t,\boldsymbol w) w ↦ f ( t , w ) はΩ t \Omega_t Ω t 上でC r C^r C r 級である。
各コンパクト集合K ⊆ Ω K\subseteq\Omega K ⊆ Ω に対して非負関数m K ∈ L 1 ( R ) m_K\in L^1(\R) m K ∈ L 1 ( R ) が存在し、( t , w ) ∈ K (t,\boldsymbol w)\in K ( t , w ) ∈ K かつt ∉ N t\notin N t ∈ / N ならば、0 ≤ j ≤ r 0\le j\le r 0 ≤ j ≤ r とすべての( i 1 , … , i j ) ∈ { 1 , … , d + p } j (i_1,\ldots,i_j)\in\{1,\ldots,d+p\}^j ( i 1 , … , i j ) ∈ { 1 , … , d + p } j について
∥ ∂ i j ⋯ ∂ i 1 f ( t , w ) ∥ ≤ m K ( t ) \lVert\partial_{i_j}\cdots\partial_{i_1}f(t,\boldsymbol w)\rVert\le m_K(t) ∥ ∂ i j ⋯ ∂ i 1 f ( t , w )∥ ≤ m K ( t )
が成り立つ。j = 0 j=0 j = 0 の項はf ( t , w ) f(t,\boldsymbol w) f ( t , w ) とする。
t ∈ N t\in N t ∈ N では、f f f の偏導関数の値を0 0 0 と定める。
補題 3.3. d , p , r ∈ N ≥ 1 d,p,r\in\NN d , p , r ∈ N ≥ 1 とし、開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d が、零集合N N N と支配関数m K m_K m K をもって状態とパラメータについてC r C^r C r 級の Carathéodory 条件を満たすとする。このとき、次が成り立つ。
f f f はパラメータ付き Carathéodory 条件を局所的に満たす。さらに、Ω \Omega Ω に含まれる各閉矩形R = J × B × C R=J\times B\times C R = J × B × C について、定義 1.1 条件 (c) の関数としてm = m R ∣ J m=m_R|_J m = m R ∣ J 、ℓ = c m R ∣ J \ell=c\,m_R|_J ℓ = c m R ∣ J を取ることができる。ここでc c c はd d d とR d \R^d R d のノルムだけで決まる定数である。
1 ≤ j ≤ r 1\le j\le r 1 ≤ j ≤ r 、( i 1 , … , i j ) ∈ { 1 , … , d + p } j (i_1,\ldots,i_j)\in\{1,\ldots,d+p\}^j ( i 1 , … , i j ) ∈ { 1 , … , d + p } j 、w ∈ R d × R p \boldsymbol w\in\R^d\times\R^p w ∈ R d × R p について、t ↦ ∂ i j ⋯ ∂ i 1 f ( t , w ) t\mapsto\partial_{i_j}\cdots\partial_{i_1}f(t,\boldsymbol w) t ↦ ∂ i j ⋯ ∂ i 1 f ( t , w ) は開集合{ t : ( t , w ) ∈ Ω } \{t:(t,\boldsymbol w)\in\Omega\} { t : ( t , w ) ∈ Ω } 上で Lebesgue 可測である。
f f f と各∂ i f \partial_if ∂ i f は、n = d + p n=d+p n = d + p として補題 1.2 (2) のg g g の仮定を満たす。r ≥ 2 r\ge2 r ≥ 2 ならば、各i i i について∂ i f \partial_if ∂ i f は状態とパラメータについてC r − 1 C^{r-1} C r − 1 級の Carathéodory 条件を満たす。
f f f と、1 ≤ j ≤ r 1\le j\le r 1 ≤ j ≤ r のすべての順序付き偏導関数∂ i j ⋯ ∂ i 1 f \partial_{i_j}\cdots\partial_{i_1}f ∂ i j ⋯ ∂ i 1 f がΩ \Omega Ω 上で存在して連続ならば、f f f はN = ∅ N=\varnothing N = ∅ として状態とパラメータについてC r C^r C r 級の Carathéodory 条件を満たす。
証明. 主張 (1) を示す。Ω \Omega Ω に含まれる閉矩形R = J × B × C R=J\times B\times C R = J × B × C を取る。R R R はコンパクトである。各( u , λ ) ∈ B × C (\boldsymbol u,\lambda)\in B\times C ( u , λ ) ∈ B × C についてJ ⊆ { t : ( t , u , λ ) ∈ Ω } J\subseteq\{t:(t,\boldsymbol u,\lambda)\in\Omega\} J ⊆ { t : ( t , u , λ ) ∈ Ω } であるから、定義 3.2 条件 (a) によりt ↦ f ( t , u , λ ) t\mapsto f(t,\boldsymbol u,\lambda) t ↦ f ( t , u , λ ) はJ J J 上で可測である。t ∈ J ∖ N t\in J\setminus N t ∈ J ∖ N ではf ( t , ⋅ ) f(t,\cdot) f ( t , ⋅ ) はC r C^r C r 級であり、特にB × C B\times C B × C 上で連続である。定義 3.2 条件 (c) をK = R K=R K = R に適用すると、t ∈ J ∖ N t\in J\setminus N t ∈ J ∖ N と( u , λ ) ∈ B × C (\boldsymbol u,\lambda)\in B\times C ( u , λ ) ∈ B × C について∥ f ( t , u , λ ) ∥ ≤ m R ( t ) \lVert f(t,\boldsymbol u,\lambda)\rVert\le m_R(t) ∥ f ( t , u , λ )∥ ≤ m R ( t ) である。Euclid ノルムを∣ ⋅ ∣ \lvert\cdot\rvert ∣ ⋅ ∣ と書き、すべてのx ∈ R d \boldsymbol x\in\R^d x ∈ R d で∣ x ∣ ≤ c 1 ∥ x ∥ \lvert\boldsymbol x\rvert\le c_1\lVert\boldsymbol x\rVert ∣ x ∣ ≤ c 1 ∥ x ∥ かつ∥ x ∥ ≤ c 1 ∣ x ∣ \lVert\boldsymbol x\rVert\le c_1\lvert\boldsymbol x\rvert ∥ x ∥ ≤ c 1 ∣ x ∣ を満たすc 1 ≥ 1 c_1\ge1 c 1 ≥ 1 を取る。t ∈ J ∖ N t\in J\setminus N t ∈ J ∖ N 、λ ∈ C \lambda\in C λ ∈ C 、u , v ∈ B \boldsymbol u,\boldsymbol v\in B u , v ∈ B とする。開集合U = { u ′ : ( t , u ′ , λ ) ∈ Ω } U=\{\boldsymbol u':(t,\boldsymbol u',\lambda)\in\Omega\} U = { u ′ : ( t , u ′ , λ ) ∈ Ω } は凸集合B B B を含み、Φ ( u ′ ) = f ( t , u ′ , λ ) \Phi(\boldsymbol u')=f(t,\boldsymbol u',\lambda) Φ ( u ′ ) = f ( t , u ′ , λ ) はU U U 上でC 1 C^1 C 1 級である。u ′ ∈ B \boldsymbol u'\in B u ′ ∈ B では、d × d d\times d d × d 行列M M M の Euclid ノルムに関する作用素ノルムが∑ i ∣ M e i ∣ \sum_i\lvert M\boldsymbol e_i\rvert ∑ i ∣ M e i ∣ 以下であることから
∥ D Φ ( u ′ ) ∥ o p ≤ ∑ i = 1 d ∣ ∂ i f ( t , u ′ , λ ) ∣ ≤ c 1 d m R ( 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) ∥ D Φ ( u ′ ) ∥ op ≤ i = 1 ∑ d ∣ ∂ i f ( t , u ′ , λ )∣ ≤ c 1 d m R ( t ) である。§E4.3 定理 2.1 により∣ Φ ( u ) − Φ ( v ) ∣ ≤ c 1 d m R ( t ) ∣ u − v ∣ \lvert\Phi(\boldsymbol u)-\Phi(\boldsymbol v)\rvert\le c_1d\,m_R(t)\lvert\boldsymbol u-\boldsymbol v\rvert ∣ Φ ( u ) − Φ ( v )∣ ≤ c 1 d m R ( t ) ∣ u − v ∣ であり、c = c 1 3 d c=c_1^3d c = c 1 3 d とおくと∥ f ( t , u , λ ) − f ( t , v , λ ) ∥ ≤ c m R ( 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 ∥ f ( t , u , λ ) − f ( t , v , λ )∥ ≤ c m R ( t ) ∥ u − v ∥ を得る。
主張 (2) を示す。1 ≤ j ≤ r 1\le j\le r 1 ≤ j ≤ r とし、g = ∂ i j − 1 ⋯ ∂ i 1 f g=\partial_{i_{j-1}}\cdots\partial_{i_1}f g = ∂ i j − 1 ⋯ ∂ i 1 f (j = 1 j=1 j = 1 ではg = f g=f g = f )について、各w ′ \boldsymbol w' w ′ でt ↦ g ( t , w ′ ) t\mapsto g(t,\boldsymbol w') t ↦ g ( t , w ′ ) がO w ′ = { t : ( t , w ′ ) ∈ Ω } O_{\boldsymbol w'}=\{t:(t,\boldsymbol w')\in\Omega\} O w ′ = { t : ( t , w ′ ) ∈ Ω } 上で可測であると仮定する。j = 1 j=1 j = 1 ではこの仮定は定義 3.2 条件 (a) である。w \boldsymbol w w を固定し、e \boldsymbol e e をw i j w_{i_j} w i j 方向の単位ベクトルとする。n ∈ N ≥ 1 n\in\NN n ∈ N ≥ 1 に対してO w O_{\boldsymbol w} O w 上の関数h n h_n h n を、t ∈ ( O w ∩ O w + e / n ) ∖ N t\in(O_{\boldsymbol w}\cap O_{\boldsymbol w+\boldsymbol e/n})\setminus N t ∈ ( O w ∩ O w + e / n ) ∖ N では
h n ( 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) h n ( t ) = n ( g ( t , w + e / n ) − g ( t , w ) ) とし、その他のt t t ではh n ( t ) = 0 h_n(t)=0 h n ( t ) = 0 と定める。O w + e / n O_{\boldsymbol w+\boldsymbol e/n} O w + e / n は開集合であるからh n h_n h n は可測である。t ∈ O w ∖ N t\in O_{\boldsymbol w}\setminus N t ∈ O w ∖ N ならば、Ω \Omega Ω が開集合であることから十分大きいn n n でt ∈ O w + e / n t\in O_{\boldsymbol w+\boldsymbol e/n} t ∈ O w + e / n であり、g ( t , ⋅ ) g(t,\cdot) g ( t , ⋅ ) はΩ t \Omega_t Ω t 上でC r − j + 1 C^{r-j+1} C r − j + 1 級であるからh n ( t ) → ∂ i j g ( t , w ) h_n(t)\to\partial_{i_j}g(t,\boldsymbol w) h n ( t ) → ∂ i j g ( t , w ) である。t ∈ N t\in N t ∈ N ではh n ( t ) = 0 h_n(t)=0 h n ( t ) = 0 である。従ってh n h_n h n はO w O_{\boldsymbol w} O w の全点で∂ i j g ( ⋅ , w ) \partial_{i_j}g(\cdot,\boldsymbol w) ∂ i j g ( ⋅ , w ) へ収束し、§E10.7 定理 5.1 を各成分に適用すると、極限は可測である。j = 1 , 2 , … , r j=1,2,\ldots,r j = 1 , 2 , … , r の順にこの議論を適用して結論を得る。
主張 (3) を示す。Ω \Omega Ω に含まれるJ × B J\times B J × B (B ⊆ R d × R p B\subseteq\R^d\times\R^p B ⊆ R d × R p は閉球)について、f f f と∂ i f \partial_if ∂ i f の時間についての可測性は定義 3.2 条件 (a) と主張 (2) をJ J J へ制限して得られ、t ∈ J ∖ N t\in J\setminus N t ∈ J ∖ N でのB B B 上の連続性は定義 3.2 条件 (b) による。r ≥ 2 r\ge2 r ≥ 2 とする。主張 (2) により∂ i f \partial_if ∂ i f は定義 3.2 条件 (a) を満たす。t ∉ N t\notin N t ∈ / N では∂ i f ( t , ⋅ ) \partial_if(t,\cdot) ∂ i f ( t , ⋅ ) はΩ t \Omega_t Ω t 上でC r − 1 C^{r-1} C r − 1 級である。∂ i f \partial_if ∂ i f のr − 1 r-1 r − 1 階以下の順序付き偏導関数はf f f のr r r 階以下の順序付き偏導関数であるから、同じm K m_K m K で支配される。
主張 (4) を示す。各w \boldsymbol w w について、t ↦ f ( t , w ) t\mapsto f(t,\boldsymbol w) t ↦ f ( t , w ) は開集合{ t : ( t , w ) ∈ Ω } \{t:(t,\boldsymbol w)\in\Omega\} { t : ( t , w ) ∈ Ω } 上で連続であるから可測である。各t t t についてf ( t , ⋅ ) f(t,\cdot) f ( t , ⋅ ) のr r r 階以下の順序付き偏導関数はΩ t \Omega_t Ω t 上で存在して連続であるから、§E4.4 定理 1.7 によりf ( t , ⋅ ) f(t,\cdot) f ( t , ⋅ ) はC r C^r C r 級である。コンパクト集合K ⊆ Ω K\subseteq\Omega K ⊆ Ω に対して、K K K 上の∥ ∂ i j ⋯ ∂ i 1 f ∥ \lVert\partial_{i_j}\cdots\partial_{i_1}f\rVert ∥ ∂ i j ⋯ ∂ i 1 f ∥ (0 ≤ j ≤ r 0\le j\le r 0 ≤ j ≤ r )の最大値をM K M_K M K とし、K K K の時間成分への射影π ( K ) \pi(K) π ( K ) を用いてm K = M K 1 π ( K ) m_K=M_K\mathbf 1_{\pi(K)} m K = M K 1 π ( K ) とおく。π ( K ) \pi(K) π ( K ) はコンパクトであるからm K ∈ L 1 ( R ) m_K\in L^1(\R) m K ∈ L 1 ( R ) である。▨
定義 3.4. d , p ∈ N ≥ 1 d,p\in\NN d , p ∈ N ≥ 1 とし、開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d が状態とパラメータについてC 1 C^1 C 1 級の Carathéodory 条件を満たすとする。φ \varphi φ をf f f の解写像とし、( t 0 , η , λ ) ∈ Ω (t_0,\boldsymbol\eta,\lambda)\in\Omega ( t 0 , η , λ ) ∈ Ω 、I = I ( t 0 , η , λ ) I=I(t_0,\boldsymbol\eta,\lambda) I = I ( t 0 , η , λ ) 、y = φ ( ⋅ , t 0 , η , λ ) \boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda) y = φ ( ⋅ , t 0 , η , λ ) とする。D u f D_{\boldsymbol u}f D u f 、D λ f D_\lambda f D λ f を、u \boldsymbol u u の座標とλ \lambda λ の座標についてのf f f の偏導関数を列とする行列とし、t ∈ I t\in I t ∈ I に対して
A ( t ) = D u f ( t , y ( t ) , λ ) ∈ R d × d , B ( t ) = D λ f ( t , y ( t ) , λ ) ∈ R d × p A(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} A ( t ) = D u f ( t , y ( t ) , λ ) ∈ R d × d , B ( t ) = D λ f ( t , y ( t ) , λ ) ∈ R d × p とおく。R d × d \R^{d\times d} R d × d に値をとるY Y Y とR d × p \R^{d\times p} R d × p に値をとるZ Z Z についての線形方程式
Y ′ = A ( t ) Y , Z ′ = A ( t ) Z + B ( t ) Y'=A(t)Y,\qquad Z'=A(t)Z+B(t) Y ′ = A ( t ) Y , Z ′ = A ( t ) Z + B ( t ) を、y \boldsymbol y y に沿う変分方程式 (variational equation ) といい、A , B A,B A , B をその係数という。変分方程式の解とは、I I I 上で局所絶対連続であり、I I I のほとんど至る所でこれらの方程式を満たすY , Z Y,Z Y , Z をいう。
補題 3.5. d , p ∈ N ≥ 1 d,p\in\NN d , p ∈ N ≥ 1 とし、開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d が状態とパラメータについてC 1 C^1 C 1 級の Carathéodory 条件を満たすとし、φ \varphi φ をその解写像とする。( t 0 , η , λ ) ∈ Ω (t_0,\boldsymbol\eta,\lambda)\in\Omega ( t 0 , η , λ ) ∈ Ω 、I = I ( t 0 , η , λ ) I=I(t_0,\boldsymbol\eta,\lambda) I = I ( t 0 , η , λ ) とし、A , B A,B A , B をy = φ ( ⋅ , t 0 , η , λ ) \boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda) y = φ ( ⋅ , t 0 , η , λ ) に沿う変分方程式の係数とする。このとき、A ∈ L l o c 1 ( I ; R d × d ) A\in L^1_{\mathrm{loc}}(I;\R^{d\times d}) A ∈ L loc 1 ( I ; R d × d ) 、B ∈ L l o c 1 ( I ; R d × p ) B\in L^1_{\mathrm{loc}}(I;\R^{d\times p}) B ∈ L loc 1 ( I ; R d × p ) であり、任意のt ∗ ∈ I t_*\in I t ∗ ∈ I 、Y ∗ ∈ R d × d Y_*\in\R^{d\times d} Y ∗ ∈ R d × d 、Z ∗ ∈ R d × p Z_*\in\R^{d\times p} Z ∗ ∈ R d × p に対して、Y ( t ∗ ) = Y ∗ Y(t_*)=Y_* Y ( t ∗ ) = Y ∗ 、Z ( t ∗ ) = Z ∗ Z(t_*)=Z_* Z ( t ∗ ) = Z ∗ を満たすy \boldsymbol y y に沿う変分方程式の解Y , Z Y,Z Y , Z が一意に存在する。
証明. N N N を定義 3.2 条件 (b) の零集合とする。a < b a<b a < b 、[ a , b ] ⊆ I [a,b]\subseteq I [ a , b ] ⊆ I を取る。補題 3.3 (3) によりf f f の各偏導関数は補題 1.2 (2) の仮定を満たし、t ↦ ( y ( t ) , λ ) t\mapsto(\boldsymbol y(t),\lambda) t ↦ ( y ( t ) , λ ) は連続で( t , y ( t ) , λ ) ∈ Ω (t,\boldsymbol y(t),\lambda)\in\Omega ( t , y ( t ) , λ ) ∈ Ω であるから、A , B A,B A , B の各列は[ a , b ] [a,b] [ a , b ] 上で可測である。K = { ( t , y ( t ) , λ ) : a ≤ t ≤ b } K=\{(t,\boldsymbol y(t),\lambda):a\le t\le b\} K = {( t , y ( t ) , λ ) : a ≤ t ≤ b } はコンパクトでΩ \Omega Ω に含まれるから、定義 3.2 条件 (c) によりt ∈ [ a , b ] ∖ N t\in[a,b]\setminus N t ∈ [ a , b ] ∖ N で各列のノルムはm K ( t ) m_K(t) m K ( t ) 以下である。従ってA , B A,B A , B は[ a , b ] [a,b] [ a , b ] 上で可積分である。[ a , b ] [a,b] [ a , b ] は任意であるからA ∈ L l o c 1 ( I ; R d × d ) A\in L^1_{\mathrm{loc}}(I;\R^{d\times d}) A ∈ L loc 1 ( I ; R d × d ) 、B ∈ L l o c 1 ( I ; R d × p ) B\in L^1_{\mathrm{loc}}(I;\R^{d\times p}) B ∈ L loc 1 ( I ; R d × p ) である。解の一意存在は補題 3.1 (3) を行列値の方程式に適用して得る。▨
補題 3.6. d , p ∈ N ≥ 1 d,p\in\NN d , p ∈ N ≥ 1 とし、R d \R^d R d とR p \R^p R p にノルム∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ を固定し、R d × R p \R^d\times\R^p R d × R p にノルム∥ ( u , μ ) ∥ = max { ∥ u ∥ , ∥ μ ∥ } \lVert(\boldsymbol u,\mu)\rVert=\max\{\lVert\boldsymbol u\rVert,\lVert\mu\rVert\} ∥( u , μ )∥ = max {∥ u ∥ , ∥ μ ∥} を入れ、d × d d\times d d × d 行列とd × p d\times p d × p 行列には作用素ノルムを用いる。開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d が、零集合N N N をもって状態とパラメータについてC 1 C^1 C 1 級の Carathéodory 条件を満たすとする。a < b a<b a < b 、λ 0 ∈ R p \lambda_0\in\R^p λ 0 ∈ R p とし、y : [ a , b ] → R d \boldsymbol y:[a,b]\to\R^d y : [ a , b ] → R d を、すべてのs ∈ [ a , b ] s\in[a,b] s ∈ [ a , b ] で( s , y ( s ) , λ 0 ) ∈ Ω (s,\boldsymbol y(s),\lambda_0)\in\Omega ( s , y ( s ) , λ 0 ) ∈ Ω を満たす連続写像とする。このとき、ρ > 0 \rho>0 ρ > 0 と非負関数m ∈ L 1 ( [ a , b ] ) m\in L^1([a,b]) m ∈ 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\} K ρ = {( s , u , μ ) : a ≤ s ≤ b , ∥ u − y ( s )∥ ≤ ρ , ∥ μ − λ 0 ∥ ≤ ρ } について次が成り立つ。
K ρ ⊆ Ω K_\rho\subseteq\Omega K ρ ⊆ Ω であり、s ∈ [ a , b ] ∖ N s\in[a,b]\setminus N s ∈ [ a , b ] ∖ N かつ( s , u , μ ) ∈ K ρ (s,\boldsymbol u,\mu)\in K_\rho ( s , u , μ ) ∈ K ρ ならば、∥ D u f ( s , u , μ ) ∥ ≤ m ( s ) \lVert D_{\boldsymbol u}f(s,\boldsymbol u,\mu)\rVert\le m(s) ∥ D u f ( s , u , μ )∥ ≤ m ( s ) かつ∥ D λ f ( s , u , μ ) ∥ ≤ m ( s ) \lVert D_\lambda f(s,\boldsymbol u,\mu)\rVert\le m(s) ∥ D λ f ( s , u , μ )∥ ≤ m ( s ) である。
写像γ : [ a , b ] × [ 0 , 1 ] → R d × R p \boldsymbol\gamma:[a,b]\times[0,1]\to\R^d\times\R^p γ : [ a , b ] × [ 0 , 1 ] → R d × R p は、各θ ∈ [ 0 , 1 ] \theta\in[0,1] θ ∈ [ 0 , 1 ] についてs ↦ γ ( s , θ ) s\mapsto\boldsymbol\gamma(s,\theta) s ↦ γ ( s , θ ) が連続であり、各s ∈ [ a , b ] s\in[a,b] s ∈ [ a , b ] についてθ ↦ γ ( s , θ ) \theta\mapsto\boldsymbol\gamma(s,\theta) θ ↦ γ ( s , θ ) が連続であり、すべての( s , θ ) (s,\theta) ( s , θ ) で( s , γ ( s , θ ) ) ∈ K ρ (s,\boldsymbol\gamma(s,\theta))\in K_\rho ( s , γ ( s , θ )) ∈ K ρ を満たすとする。s ∈ [ a , b ] ∖ N s\in[a,b]\setminus N s ∈ [ a , b ] ∖ N では
A γ ( s ) = ∫ 0 1 D u f ( s , γ ( s , θ ) ) d θ , B γ ( s ) = ∫ 0 1 D λ 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 A γ ( s ) = ∫ 0 1 D u f ( s , γ ( s , θ )) d θ , B γ ( s ) = ∫ 0 1 D λ f ( s , γ ( s , θ )) d θ
とし、s ∈ [ a , b ] ∩ N s\in[a,b]\cap N s ∈ [ a , b ] ∩ N ではA γ ( s ) = 0 A_{\boldsymbol\gamma}(s)=0 A γ ( s ) = 0 、B γ ( s ) = 0 B_{\boldsymbol\gamma}(s)=0 B γ ( s ) = 0 とする。このとき、A γ A_{\boldsymbol\gamma} A γ とB γ B_{\boldsymbol\gamma} B γ は[ a , b ] [a,b] [ a , b ] 上で可測であり、s ∈ [ a , b ] ∖ N s\in[a,b]\setminus N s ∈ [ a , b ] ∖ N では∥ A γ ( s ) ∥ ≤ m ( s ) \lVert A_{\boldsymbol\gamma}(s)\rVert\le m(s) ∥ A γ ( s )∥ ≤ m ( s ) かつ∥ B γ ( s ) ∥ ≤ m ( s ) \lVert B_{\boldsymbol\gamma}(s)\rVert\le m(s) ∥ B γ ( s )∥ ≤ m ( s ) である。特に、A γ A_{\boldsymbol\gamma} A γ とB γ B_{\boldsymbol\gamma} B γ は[ a , b ] [a,b] [ a , b ] 上で可積分である。
γ 0 ( s , θ ) = ( y ( s ) , λ 0 ) \boldsymbol\gamma_0(s,\theta)=(\boldsymbol y(s),\lambda_0) γ 0 ( s , θ ) = ( y ( s ) , λ 0 ) とおく。γ n \boldsymbol\gamma_n γ n (n ∈ N ≥ 1 n\in\NN n ∈ N ≥ 1 )は主張 (2) の仮定を満たし、すべてのs ∈ [ a , b ] s\in[a,b] s ∈ [ a , b ] についてn → ∞ n\to\infty n → ∞ のとき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 sup 0 ≤ θ ≤ 1 ∥ γ n ( s , θ ) − ( y ( s ) , λ 0 )∥ → 0 であるとする。このとき、n → ∞ n\to\infty n → ∞ のとき∫ a b ∥ A γ n − A γ 0 ∥ → 0 \int_a^b\lVert A_{\boldsymbol\gamma_n}-A_{\boldsymbol\gamma_0}\rVert\to0 ∫ a b ∥ A γ n − A γ 0 ∥ → 0 かつ∫ a b ∥ B γ n − B γ 0 ∥ → 0 \int_a^b\lVert B_{\boldsymbol\gamma_n}-B_{\boldsymbol\gamma_0}\rVert\to0 ∫ a b ∥ B γ n − B γ 0 ∥ → 0 である。
証明. 主張 (1) を示す。補題 3.3 (1) によりf f f はパラメータ付き Carathéodory 条件を局所的に満たすから、補題 1.3 (3) をy \boldsymbol y y とλ 0 \lambda_0 λ 0 に適用すると、K ρ ⊆ Ω K_\rho\subseteq\Omega K ρ ⊆ Ω を満たすρ > 0 \rho>0 ρ > 0 が存在する。K ρ K_\rho K ρ は有界な閉集合であるからコンパクトであり、定義 3.2 条件 (c) の支配関数m K ρ m_{K_\rho} m K ρ が存在する。ノルムの同値性により、d × d d\times d d × d 行列とd × p d\times p d × p 行列の作用素ノルムがその列のノルムの最大値のc c c 倍以下となる定数c > 0 c>0 c > 0 が存在する。D u f D_{\boldsymbol u}f D u f とD λ f D_\lambda f D λ f の列はf f f の偏導関数であるから、m = c m K ρ ∣ [ a , b ] m=c\,m_{K_\rho}|_{[a,b]} m = c m K ρ ∣ [ a , b ] とおくと結論を得る。
主張 (2) を示す。G G G をD u f D_{\boldsymbol u}f D u f またはD λ f D_\lambda f D λ f とし、s ∈ [ a , b ] ∖ N s\in[a,b]\setminus N s ∈ [ a , b ] ∖ N ではG γ ( s ) = ∫ 0 1 G ( s , γ ( s , θ ) ) d θ G_{\boldsymbol\gamma}(s)=\int_0^1G(s,\boldsymbol\gamma(s,\theta))\,d\theta G γ ( s ) = ∫ 0 1 G ( s , γ ( s , θ )) d θ 、s ∈ [ a , b ] ∩ N s\in[a,b]\cap N s ∈ [ a , b ] ∩ N ではG γ ( s ) = 0 G_{\boldsymbol\gamma}(s)=0 G γ ( s ) = 0 とおく。k ∈ N ≥ 1 k\in\NN k ∈ N ≥ 1 とs ∈ [ a , b ] s\in[a,b] s ∈ [ a , b ] に対して
R k ( s ) = 1 k ∑ i = 1 k G ( s , γ ( s , i / k ) ) R_k(s)=\frac1k\sum_{i=1}^kG\bigl(s,\boldsymbol\gamma(s,i/k)\bigr) R k ( s ) = k 1 i = 1 ∑ k G ( s , γ ( s , i / k ) ) とおく。s ∈ [ a , b ] ∖ N s\in[a,b]\setminus N s ∈ [ a , b ] ∖ N では、G ( s , ⋅ ) G(s,\cdot) G ( s , ⋅ ) はΩ s \Omega_s Ω s 上で連続であり、θ ↦ γ ( s , θ ) \theta\mapsto\boldsymbol\gamma(s,\theta) θ ↦ γ ( s , θ ) はΩ s \Omega_s Ω s に値をとる連続写像であるから、被積分関数はθ \theta θ について連続であり、G γ ( s ) G_{\boldsymbol\gamma}(s) G γ ( s ) は Riemann 和R k ( s ) R_k(s) R k ( s ) のk → ∞ k\to\infty k → ∞ での極限である。主張 (1) により∥ R k ( s ) ∥ ≤ m ( s ) \lVert R_k(s)\rVert\le m(s) ∥ R k ( s )∥ ≤ m ( s ) であるから、∥ G γ ( s ) ∥ ≤ m ( s ) \lVert G_{\boldsymbol\gamma}(s)\rVert\le m(s) ∥ G γ ( s )∥ ≤ m ( s ) である。s ∈ [ a , b ] ∩ N s\in[a,b]\cap N s ∈ [ a , b ] ∩ N では、偏導関数の値の規約によりR k ( s ) = 0 = G γ ( s ) R_k(s)=0=G_{\boldsymbol\gamma}(s) R k ( s ) = 0 = G γ ( s ) である。各i i i についてs ↦ γ ( s , i / k ) s\mapsto\boldsymbol\gamma(s,i/k) s ↦ γ ( s , i / k ) は連続であるから、補題 3.3 (3) と補題 1.2 (2) によりs ↦ G ( s , γ ( s , i / k ) ) s\mapsto G(s,\boldsymbol\gamma(s,i/k)) s ↦ G ( s , γ ( s , i / k )) の各列は[ a , b ] [a,b] [ a , b ] 上で可測であり、R k R_k R k は可測である。R k R_k R k は[ a , b ] [a,b] [ a , b ] の各点でG γ G_{\boldsymbol\gamma} G γ へ収束するから、§E10.7 定理 5.1 を各成分に適用すると、G γ G_{\boldsymbol\gamma} G γ は可測である。G γ G_{\boldsymbol\gamma} G γ は可測でありm ∈ L 1 ( [ a , b ] ) m\in L^1([a,b]) m ∈ L 1 ([ a , b ]) で支配されるから、[ a , b ] [a,b] [ a , b ] 上で可積分である。
主張 (3) を示す。G G G を前段のとおりとする。γ 0 \boldsymbol\gamma_0 γ 0 は主張 (2) の仮定を満たし、s ∈ [ a , b ] ∖ N s\in[a,b]\setminus N s ∈ [ a , b ] ∖ N ではG γ 0 ( s ) = G ( s , y ( s ) , λ 0 ) G_{\boldsymbol\gamma_0}(s)=G(s,\boldsymbol y(s),\lambda_0) G γ 0 ( s ) = G ( s , y ( s ) , λ 0 ) である。s ∈ [ a , b ] ∖ N s\in[a,b]\setminus N s ∈ [ a , b ] ∖ N とε > 0 \varepsilon>0 ε > 0 を取る。G ( s , ⋅ ) G(s,\cdot) G ( s , ⋅ ) は( y ( s ) , λ 0 ) (\boldsymbol y(s),\lambda_0) ( y ( s ) , λ 0 ) で連続であるから、w ∈ Ω s \boldsymbol w\in\Omega_s w ∈ Ω s かつ∥ w − ( y ( s ) , λ 0 ) ∥ ≤ δ \lVert\boldsymbol w-(\boldsymbol y(s),\lambda_0)\rVert\le\delta ∥ w − ( y ( s ) , λ 0 )∥ ≤ δ ならば∥ G ( s , w ) − G ( s , y ( s ) , λ 0 ) ∥ ≤ ε \lVert G(s,\boldsymbol w)-G(s,\boldsymbol y(s),\lambda_0)\rVert\le\varepsilon ∥ G ( s , w ) − G ( s , y ( s ) , λ 0 )∥ ≤ ε となるδ > 0 \delta>0 δ > 0 が存在する。十分大きいn n n ではsup θ ∥ γ n ( s , θ ) − ( y ( s ) , λ 0 ) ∥ ≤ δ \sup_\theta\lVert\boldsymbol\gamma_n(s,\theta)-(\boldsymbol y(s),\lambda_0)\rVert\le\delta sup θ ∥ γ n ( s , θ ) − ( y ( s ) , λ 0 )∥ ≤ δ であり、G γ n ( s ) − G γ 0 ( s ) G_{\boldsymbol\gamma_n}(s)-G_{\boldsymbol\gamma_0}(s) G γ n ( s ) − G γ 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) θ ↦ G ( s , γ n ( s , θ )) − G ( s , y ( s ) , λ 0 ) の Riemann 和の極限で、各 Riemann 和のノルムはε \varepsilon ε 以下であるから、∥ G γ n ( s ) − G γ 0 ( s ) ∥ ≤ ε \lVert G_{\boldsymbol\gamma_n}(s)-G_{\boldsymbol\gamma_0}(s)\rVert\le\varepsilon ∥ G γ n ( s ) − G γ 0 ( s )∥ ≤ ε である。従って、[ a , b ] ∖ N [a,b]\setminus N [ a , b ] ∖ N の各点で∥ G γ n − G γ 0 ∥ → 0 \lVert G_{\boldsymbol\gamma_n}-G_{\boldsymbol\gamma_0}\rVert\to0 ∥ G γ n − G γ 0 ∥ → 0 である。主張 (2) によりs ↦ ∥ G γ n ( s ) − G γ 0 ( s ) ∥ s\mapsto\lVert G_{\boldsymbol\gamma_n}(s)-G_{\boldsymbol\gamma_0}(s)\rVert s ↦ ∥ G γ n ( s ) − G γ 0 ( s )∥ は可測であり、[ a , b ] ∖ N [a,b]\setminus N [ a , b ] ∖ N で2 m 2m 2 m 以下であるから、§E10.11 定理 3.1 により∫ a b ∥ G γ n − G γ 0 ∥ → 0 \int_a^b\lVert G_{\boldsymbol\gamma_n}-G_{\boldsymbol\gamma_0}\rVert\to0 ∫ a b ∥ G γ n − G γ 0 ∥ → 0 である。▨
補題 3.7. d , p ∈ N ≥ 1 d,p\in\NN d , p ∈ N ≥ 1 とし、開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d が状態とパラメータについてC 1 C^1 C 1 級の Carathéodory 条件を満たすとし、φ : W → R d \varphi:W\to\R^d φ : W → R d をその解写像とする。( t , t 0 , η , λ ) ∈ W (t,t_0,\boldsymbol\eta,\lambda)\in W ( t , t 0 , η , λ ) ∈ W に対して、φ ( ⋅ , t 0 , η , λ ) \varphi(\cdot,t_0,\boldsymbol\eta,\lambda) φ ( ⋅ , t 0 , η , λ ) に沿う変分方程式の解でY ( t 0 ) = I d Y(t_0)=I_d Y ( t 0 ) = I d 、Z ( t 0 ) = 0 Z(t_0)=0 Z ( t 0 ) = 0 を満たすものの時刻t t t における値をY ( t , t 0 , η , λ ) Y(t,t_0,\boldsymbol\eta,\lambda) Y ( t , t 0 , η , λ ) 、Z ( t , t 0 , η , λ ) Z(t,t_0,\boldsymbol\eta,\lambda) Z ( t , t 0 , η , λ ) と書く。このときY Y Y とZ Z Z はW W W 上で連続である。
証明. W W W の点列( t n , t 0 , n , η n , λ n ) (t_n,t_{0,n},\boldsymbol\eta_n,\lambda_n) ( t n , t 0 , n , η n , λ n ) が( t , t 0 , η , λ ) ∈ W (t,t_0,\boldsymbol\eta,\lambda)\in W ( t , t 0 , η , λ ) ∈ W へ収束するとする。a < min { t , t 0 } a<\min\{t,t_0\} a < min { t , t 0 } 、max { t , t 0 } < b \max\{t,t_0\}<b max { t , t 0 } < b 、[ a , b ] ⊆ I ( t 0 , η , λ ) [a,b]\subseteq I(t_0,\boldsymbol\eta,\lambda) [ a , b ] ⊆ I ( t 0 , η , λ ) を満たすa , b a,b a , b を取る。y = φ ( ⋅ , t 0 , η , λ ) \boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda) y = φ ( ⋅ , t 0 , η , λ ) とおき、A , B A,B A , B をy \boldsymbol y y に沿う変分方程式の係数とする。補題 3.6 をy ∣ [ a , b ] \boldsymbol y|_{[a,b]} y ∣ [ a , b ] とλ 0 = λ \lambda_0=\lambda λ 0 = λ に適用してρ > 0 \rho>0 ρ > 0 を取る。定理 2.2 (1) により、十分大きいn n n ではt n , t 0 , n ∈ ( a , b ) t_n,t_{0,n}\in(a,b) t n , t 0 , n ∈ ( a , b ) 、[ a , b ] ⊆ I ( t 0 , n , η n , λ n ) [a,b]\subseteq I(t_{0,n},\boldsymbol\eta_n,\lambda_n) [ a , b ] ⊆ I ( t 0 , n , η n , λ n ) であり、y n = φ ( ⋅ , t 0 , n , η n , λ n ) \boldsymbol y_n=\varphi(\cdot,t_{0,n},\boldsymbol\eta_n,\lambda_n) y n = φ ( ⋅ , t 0 , n , η n , λ n ) は[ a , b ] [a,b] [ a , b ] 上でy \boldsymbol y y へ一様収束する。番号を付け直して、すべてのn n n でこれらとmax [ a , b ] ∥ y n − y ∥ ≤ ρ \max_{[a,b]}\lVert\boldsymbol y_n-\boldsymbol y\rVert\le\rho max [ a , b ] ∥ y n − y ∥ ≤ ρ 、∥ λ n − λ ∥ ≤ ρ \lVert\lambda_n-\lambda\rVert\le\rho ∥ λ n − λ ∥ ≤ ρ が成り立つとしてよい。
A n , B n A_n,B_n A n , B n をy n \boldsymbol y_n y n に沿う変分方程式の係数とする。γ n ( s , θ ) = ( y n ( s ) , λ n ) \boldsymbol\gamma_n(s,\theta)=(\boldsymbol y_n(s),\lambda_n) γ n ( s , θ ) = ( y n ( s ) , λ n ) は補題 3.6 (2) の仮定を満たし、偏導関数の値の規約により、[ a , b ] [a,b] [ a , b ] 上でA γ n = A n A_{\boldsymbol\gamma_n}=A_n A γ n = A n 、B γ n = B n B_{\boldsymbol\gamma_n}=B_n B γ n = B n 、A γ 0 = A A_{\boldsymbol\gamma_0}=A A γ 0 = A 、B γ 0 = B B_{\boldsymbol\gamma_0}=B B γ 0 = B である。n → ∞ n\to\infty n → ∞ のときsup θ ∥ γ n ( s , θ ) − ( y ( s ) , λ ) ∥ ≤ max { max [ a , b ] ∥ y n − 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 sup θ ∥ γ n ( s , θ ) − ( y ( s ) , λ )∥ ≤ max { max [ a , b ] ∥ y n − y ∥ , ∥ λ n − λ ∥} → 0 であるから、補題 3.6 (3) により∫ a b ∥ A n − A ∥ → 0 \int_a^b\lVert A_n-A\rVert\to0 ∫ a b ∥ A n − A ∥ → 0 、∫ a b ∥ B n − B ∥ → 0 \int_a^b\lVert B_n-B\rVert\to0 ∫ a b ∥ B n − B ∥ → 0 である。
Y n , Z n Y_n,Z_n Y n , Z n をy n \boldsymbol y_n y n に沿う変分方程式の解でY n ( t 0 , n ) = I d Y_n(t_{0,n})=I_d Y n ( t 0 , n ) = I d 、Z n ( t 0 , n ) = 0 Z_n(t_{0,n})=0 Z n ( t 0 , n ) = 0 を満たすものとすると、その[ a , b ] [a,b] [ a , b ] への制限は補題 3.1 (1) の解である。補題 3.1 (2) を行列値の方程式に適用すると、Y n → Y ( ⋅ , t 0 , η , λ ) Y_n\to Y(\cdot,t_0,\boldsymbol\eta,\lambda) Y n → Y ( ⋅ , t 0 , η , λ ) 、Z n → Z ( ⋅ , t 0 , η , λ ) Z_n\to Z(\cdot,t_0,\boldsymbol\eta,\lambda) Z n → Z ( ⋅ , t 0 , η , λ ) は[ a , b ] [a,b] [ a , b ] 上で一様収束する。従って
∥ Y n ( t n ) − Y ( t , t 0 , η , λ ) ∥ ≤ max [ a , b ] ∥ Y n − Y ( ⋅ , t 0 , η , λ ) ∥ + ∥ Y ( t n , t 0 , η , λ ) − Y ( t , t 0 , η , λ ) ∥ ⟶ 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 ∥ Y n ( t n ) − Y ( t , t 0 , η , λ )∥ ≤ [ a , b ] max ∥ Y n − Y ( ⋅ , t 0 , η , λ )∥ + ∥ Y ( t n , t 0 , η , λ ) − Y ( t , t 0 , η , λ )∥ ⟶ 0 であり、Z Z Z についても同様である。▨
定理 3.8. d , p ∈ N ≥ 1 d,p\in\NN d , p ∈ N ≥ 1 とし、開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d が状態とパラメータについてC 1 C^1 C 1 級の Carathéodory 条件を満たすとし、φ : W → R d \varphi:W\to\R^d φ : W → R d をその解写像とする。このとき、各( t , t 0 ) ∈ R 2 (t,t_0)\in\R^2 ( t , t 0 ) ∈ R 2 について、W t , t 0 = { ( η , λ ) : ( t , t 0 , η , λ ) ∈ W } W_{t,t_0}=\{(\boldsymbol\eta,\lambda):(t,t_0,\boldsymbol\eta,\lambda)\in W\} W t , t 0 = {( η , λ ) : ( t , t 0 , η , λ ) ∈ W } はR d × R p \R^d\times\R^p R d × R p の開集合であり、( η , λ ) ↦ φ ( t , t 0 , η , λ ) (\boldsymbol\eta,\lambda)\mapsto\varphi(t,t_0,\boldsymbol\eta,\lambda) ( η , λ ) ↦ φ ( t , t 0 , η , λ ) はW t , t 0 W_{t,t_0} W t , t 0 上でC 1 C^1 C 1 級であって、
D η φ ( t , t 0 , η , λ ) = Y ( t ) , D λ φ ( t , t 0 , η , λ ) = 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) D η φ ( t , t 0 , η , λ ) = Y ( t ) , D λ φ ( t , t 0 , η , λ ) = Z ( t ) が成り立つ。ここでY , Z Y,Z Y , Z はφ ( ⋅ , t 0 , η , λ ) \varphi(\cdot,t_0,\boldsymbol\eta,\lambda) φ ( ⋅ , t 0 , η , λ ) に沿う変分方程式の解で、Y ( t 0 ) = I d Y(t_0)=I_d Y ( t 0 ) = I d 、Z ( t 0 ) = 0 Z(t_0)=0 Z ( t 0 ) = 0 を満たすものである。( t , t 0 , η , λ ) ↦ D η φ ( t , t 0 , η , λ ) (t,t_0,\boldsymbol\eta,\lambda)\mapsto D_{\boldsymbol\eta}\varphi(t,t_0,\boldsymbol\eta,\lambda) ( t , t 0 , η , λ ) ↦ D η φ ( t , t 0 , η , λ ) と( t , t 0 , η , λ ) ↦ D λ φ ( t , t 0 , η , λ ) (t,t_0,\boldsymbol\eta,\lambda)\mapsto D_\lambda\varphi(t,t_0,\boldsymbol\eta,\lambda) ( t , t 0 , η , λ ) ↦ D λ φ ( t , t 0 , η , λ ) はW W W 上で連続である。
証明. W t , t 0 W_{t,t_0} W t , t 0 は開集合W W W (定理 2.2 (2) )の切り口であるから開集合である。( t , t 0 , η , λ ) ∈ W (t,t_0,\boldsymbol\eta,\lambda)\in W ( t , t 0 , η , λ ) ∈ W とe = ( e ′ , e ′ ′ ) ∈ R d × R p \boldsymbol e=(\boldsymbol e',\boldsymbol e'')\in\R^d\times\R^p e = ( e ′ , e ′′ ) ∈ R d × R p を取り、y = φ ( ⋅ , t 0 , η , λ ) \boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda) y = φ ( ⋅ , t 0 , η , λ ) とおき、A , B A,B A , B をy \boldsymbol y y に沿う変分方程式の係数、N N N を定義 3.2 条件 (b) の零集合とする。a < min { t , t 0 } a<\min\{t,t_0\} a < min { t , t 0 } 、max { t , t 0 } < b \max\{t,t_0\}<b max { t , t 0 } < b 、[ a , b ] ⊆ I ( t 0 , η , λ ) [a,b]\subseteq I(t_0,\boldsymbol\eta,\lambda) [ a , b ] ⊆ I ( t 0 , η , λ ) を満たすa , b a,b a , b を取る。補題 3.6 をy ∣ [ a , b ] \boldsymbol y|_{[a,b]} y ∣ [ a , b ] とλ 0 = λ \lambda_0=\lambda λ 0 = λ に適用してρ > 0 \rho>0 ρ > 0 を取り、同補題のK ρ K_\rho K ρ とR d × R p \R^d\times\R^p R d × R p のノルムを用いる。定理 2.2 (1) により、あるh 0 > 0 h_0>0 h 0 > 0 が存在し、0 < ∣ h ∣ ≤ h 0 0<\lvert h\rvert\le h_0 0 < ∣ h ∣ ≤ h 0 ならば[ a , b ] ⊆ I ( t 0 , η + h e ′ , λ + h e ′ ′ ) [a,b]\subseteq I(t_0,\boldsymbol\eta+h\boldsymbol e',\lambda+h\boldsymbol e'') [ a , b ] ⊆ I ( t 0 , η + h e ′ , λ + h e ′′ ) であり、y h = φ ( ⋅ , t 0 , η + h e ′ , λ + h e ′ ′ ) \boldsymbol y_h=\varphi(\cdot,t_0,\boldsymbol\eta+h\boldsymbol e',\lambda+h\boldsymbol e'') y h = φ ( ⋅ , t 0 , η + h e ′ , λ + h e ′′ ) はmax [ a , b ] ∥ y h − y ∥ ≤ ρ \max_{[a,b]}\lVert\boldsymbol y_h-\boldsymbol y\rVert\le\rho max [ a , b ] ∥ y h − y ∥ ≤ ρ 、∣ h ∣ ∥ e ′ ′ ∥ ≤ ρ \lvert h\rvert\lVert\boldsymbol e''\rVert\le\rho ∣ h ∣ ∥ e ′′ ∥ ≤ ρ を満たし、h → 0 h\to0 h → 0 のときmax [ a , b ] ∥ y h − y ∥ → 0 \max_{[a,b]}\lVert\boldsymbol y_h-\boldsymbol y\rVert\to0 max [ a , b ] ∥ y h − y ∥ → 0 である。
0 < ∣ h ∣ ≤ h 0 0<\lvert h\rvert\le h_0 0 < ∣ h ∣ ≤ h 0 、s ∈ [ a , b ] s\in[a,b] s ∈ [ a , b ] 、θ ∈ [ 0 , 1 ] \theta\in[0,1] θ ∈ [ 0 , 1 ] に対して
γ h ( s , θ ) = ( y ( s ) + θ ( y h ( s ) − y ( s ) ) , λ + θ h e ′ ′ ) \boldsymbol\gamma_h(s,\theta)=\bigl(\boldsymbol y(s)+\theta(\boldsymbol y_h(s)-\boldsymbol y(s)),\ \lambda+\theta h\boldsymbol e''\bigr) γ h ( s , θ ) = ( y ( s ) + θ ( y h ( s ) − y ( s )) , λ + θ h e ′′ ) とおくと、( s , γ h ( s , θ ) ) ∈ K ρ (s,\boldsymbol\gamma_h(s,\theta))\in K_\rho ( s , γ h ( s , θ )) ∈ K ρ であり、γ h \boldsymbol\gamma_h γ h は補題 3.6 (2) の仮定を満たす。同項のA γ h A_{\boldsymbol\gamma_h} A γ h 、B γ h B_{\boldsymbol\gamma_h} B γ h をA h A_h A h 、B h B_h B h と書くと、同項によりA h , B h A_h,B_h A h , B h は[ a , b ] [a,b] [ a , b ] 上で可積分である。s ∈ [ a , b ] ∖ N s\in[a,b]\setminus N s ∈ [ a , b ] ∖ N では、§E4.3 定理 1.1 によりθ ↦ f ( s , γ h ( s , θ ) ) \theta\mapsto f(s,\boldsymbol\gamma_h(s,\theta)) θ ↦ f ( s , γ h ( s , θ )) は[ 0 , 1 ] [0,1] [ 0 , 1 ] 上でC 1 C^1 C 1 級であり、その導関数はD u f ( s , γ h ( s , θ ) ) ( y h ( s ) − y ( s ) ) + h D λ 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'' D u f ( s , γ h ( s , θ )) ( y h ( s ) − y ( s )) + h D λ f ( s , γ h ( s , θ )) e ′′ である。これを成分ごとに[ 0 , 1 ] [0,1] [ 0 , 1 ] で積分すると
f ( s , y h ( s ) , λ + h e ′ ′ ) − f ( s , y ( s ) , λ ) = A h ( s ) ( y h ( s ) − y ( s ) ) + h B h ( 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'' f ( s , y h ( s ) , λ + h e ′′ ) − f ( s , y ( s ) , λ ) = A h ( s ) ( y h ( s ) − y ( s ) ) + h B h ( s ) e ′′ を得る。
q h = ( y h − y ) / h \boldsymbol q_h=(\boldsymbol y_h-\boldsymbol y)/h q h = ( y h − y ) / h は[ a , b ] [a,b] [ a , b ] 上で絶対連続であり、q h ( t 0 ) = e ′ \boldsymbol q_h(t_0)=\boldsymbol e' q h ( t 0 ) = e ′ である。y h \boldsymbol y_h y h とy \boldsymbol y y がともに方程式を満たしs ∉ N s\notin N s ∈ / N であるs s s では、上の等式をh h h で割ると
q h ′ ( s ) = A h ( s ) q h ( s ) + B h ( s ) e ′ ′ \boldsymbol q_h'(s)=A_h(s)\boldsymbol q_h(s)+B_h(s)\boldsymbol e'' q h ′ ( s ) = A h ( s ) q h ( s ) + B h ( s ) e ′′ である。従ってq h \boldsymbol q_h q h は、係数A h A_h A h 、非斉次項B h e ′ ′ B_h\boldsymbol e'' B h e ′′ 、初期時刻t 0 t_0 t 0 、初期値e ′ \boldsymbol e' e ′ に対する補題 3.1 (1) の解である。偏導関数の値の規約により、γ 0 ( s , θ ) = ( y ( s ) , λ ) \boldsymbol\gamma_0(s,\theta)=(\boldsymbol y(s),\lambda) γ 0 ( s , θ ) = ( y ( s ) , λ ) について[ a , b ] [a,b] [ a , b ] 上でA γ 0 = A A_{\boldsymbol\gamma_0}=A A γ 0 = A 、B γ 0 = B B_{\boldsymbol\gamma_0}=B B γ 0 = B である。s ∈ [ a , b ] s\in[a,b] s ∈ [ a , b ] では、h → 0 h\to0 h → 0 のときmax θ ∈ [ 0 , 1 ] ∥ γ h ( s , θ ) − ( y ( s ) , λ ) ∥ ≤ max { max [ a , b ] ∥ y h − 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 max θ ∈ [ 0 , 1 ] ∥ γ h ( s , θ ) − ( y ( s ) , λ )∥ ≤ max { max [ a , b ] ∥ y h − y ∥ , ∣ h ∣ ∥ e ′′ ∥} → 0 である。零へ収束する任意の点列h n h_n h n (0 < ∣ h n ∣ ≤ h 0 0<\lvert h_n\rvert\le h_0 0 < ∣ h n ∣ ≤ h 0 )について、補題 3.6 (3) をγ h n \boldsymbol\gamma_{h_n} γ h n に適用すると、∫ a b ∥ A h n − A ∥ → 0 \int_a^b\lVert A_{h_n}-A\rVert\to0 ∫ a b ∥ A h n − A ∥ → 0 であり、∫ a b ∥ B h n e ′ ′ − B e ′ ′ ∥ ≤ ∥ e ′ ′ ∥ ∫ a b ∥ B h n − 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 ∫ a b ∥ B h n e ′′ − B e ′′ ∥ ≤ ∥ e ′′ ∥ ∫ a b ∥ B h n − B ∥ → 0 である。補題 3.1 (2) により、q h n \boldsymbol q_{h_n} q h n は[ a , b ] [a,b] [ a , b ] 上で、q ′ = A q + B e ′ ′ \boldsymbol q'=A\boldsymbol q+B\boldsymbol e'' q ′ = A q + B e ′′ 、q ( t 0 ) = e ′ \boldsymbol q(t_0)=\boldsymbol e' q ( t 0 ) = e ′ の解q \boldsymbol q q へ一様収束する。Y e ′ + Z e ′ ′ Y\boldsymbol e'+Z\boldsymbol e'' Y e ′ + Z e ′′ はこの線形方程式の[ a , b ] [a,b] [ a , b ] 上の解で時刻t 0 t_0 t 0 にe ′ \boldsymbol e' e ′ をとるから、一意性によりq = Y e ′ + Z e ′ ′ \boldsymbol q=Y\boldsymbol e'+Z\boldsymbol e'' q = Y e ′ + Z e ′′ である。点列h n h_n h n は任意であるから
lim h → 0 φ ( t , t 0 , η + h e ′ , λ + h e ′ ′ ) − φ ( t , t 0 , η , λ ) 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'' h → 0 lim h φ ( t , t 0 , η + h e ′ , λ + h e ′′ ) − φ ( t , t 0 , η , λ ) = Y ( t ) e ′ + Z ( t ) e ′′ である。e \boldsymbol e e をR d × R p \R^d\times\R^p R d × R p の標準基底の各ベクトルとすると、( η , λ ) ↦ φ ( t , t 0 , η , λ ) (\boldsymbol\eta,\lambda)\mapsto\varphi(t,t_0,\boldsymbol\eta,\lambda) ( η , λ ) ↦ φ ( t , t 0 , η , λ ) の偏導関数はY ( t ) Y(t) Y ( t ) とZ ( t ) Z(t) Z ( t ) の列である。補題 3.7 により、これらの偏導関数はW W W 上で連続であり、特にW t , t 0 W_{t,t_0} W t , t 0 上で連続である。§E4.2 系 2.2 により、( η , λ ) ↦ φ ( t , t 0 , η , λ ) (\boldsymbol\eta,\lambda)\mapsto\varphi(t,t_0,\boldsymbol\eta,\lambda) ( η , λ ) ↦ φ ( t , t 0 , η , λ ) はW t , t 0 W_{t,t_0} W t , t 0 上でC 1 C^1 C 1 級であり、D η φ = Y ( t ) D_{\boldsymbol\eta}\varphi=Y(t) D η φ = Y ( t ) 、D λ φ = Z ( t ) D_\lambda\varphi=Z(t) D λ φ = Z ( t ) である。▨
例 3.9. Ω = R × R × R \Omega=\R\times\R\times\R Ω = R × R × R 上のf ( t , u , λ ) = ∣ u ∣ f(t,u,\lambda)=\lvert u\rvert f ( t , u , λ ) = ∣ u ∣ は連続であり、∣ f ( t , u , λ ) − f ( t , v , λ ) ∣ ≤ ∣ u − v ∣ \lvert f(t,u,\lambda)-f(t,v,\lambda)\rvert\le\lvert u-v\rvert ∣ f ( t , u , λ ) − f ( t , v , λ )∣ ≤ ∣ u − v ∣ を満たすから、パラメータ付き Carathéodory 条件を局所的に満たす。η ∈ R \eta\in\R η ∈ R に対して
y ( t ) = { η e t , η ≥ 0 , η e − t , η < 0 y(t)=\begin{cases}\eta e^{t},&\eta\ge0,\\ \eta e^{-t},&\eta<0\end{cases} y ( t ) = { η e t , η e − t , η ≥ 0 , η < 0 はR \R R 上でC 1 C^1 C 1 級であり、η ≥ 0 \eta\ge0 η ≥ 0 ではy ≥ 0 y\ge0 y ≥ 0 、η < 0 \eta<0 η < 0 ではy < 0 y<0 y < 0 であるからy ′ = ∣ y ∣ y'=\lvert y\rvert y ′ = ∣ y ∣ 、y ( 0 ) = η y(0)=\eta y ( 0 ) = η を満たす。従ってI ( 0 , η , λ ) = R I(0,\eta,\lambda)=\R I ( 0 , η , λ ) = R であり、φ ( t , 0 , η , λ ) = y ( t ) \varphi(t,0,\eta,\lambda)=y(t) φ ( t , 0 , η , λ ) = y ( t ) である。t ≠ 0 t\ne0 t = 0 を固定すると、η ↦ φ ( t , 0 , η , λ ) \eta\mapsto\varphi(t,0,\eta,\lambda) η ↦ φ ( t , 0 , η , λ ) のη = 0 \eta=0 η = 0 における右微分係数はe t e^{t} e t 、左微分係数はe − t e^{-t} e − t であり、一致しない。定理 3.8 の結論は、状態についてのC 1 C^1 C 1 級の仮定を Lipschitz 条件へ弱めると成り立たない。一方、定理 2.2 (2) により、このφ \varphi φ はW W W 上で連続である。
例 3.10. Ω = R × R × R \Omega=\R\times\R\times\R Ω = R × R × R 上で、( t , λ ) ≠ ( 0 , 0 ) (t,\lambda)\ne(0,0) ( t , λ ) = ( 0 , 0 ) ならばf ( t , u , λ ) = λ / t 2 + λ 2 f(t,u,\lambda)=\lambda/\sqrt{t^2+\lambda^2} f ( t , u , λ ) = λ / t 2 + λ 2 、f ( 0 , u , 0 ) = 0 f(0,u,0)=0 f ( 0 , u , 0 ) = 0 とおく。f f f はu u u によらず、すべての点で∣ f ( t , u , λ ) ∣ ≤ 1 \lvert f(t,u,\lambda)\rvert\le1 ∣ f ( t , u , λ )∣ ≤ 1 である。各( u , λ ) (u,\lambda) ( u , λ ) についてt ↦ f ( t , u , λ ) t\mapsto f(t,u,\lambda) t ↦ f ( t , u , λ ) はR ∖ { 0 } \R\setminus\{0\} R ∖ { 0 } 上で連続であるから Lebesgue 可測であり、t ≠ 0 t\ne0 t = 0 では( u , λ ) ↦ f ( t , u , λ ) (u,\lambda)\mapsto f(t,u,\lambda) ( u , λ ) ↦ f ( t , u , λ ) はR 2 \R^2 R 2 上でC ∞ C^\infty C ∞ 級である。従って、f f f は各閉矩形J × B × C J\times B\times C J × B × C について零集合{ 0 } ∩ J \{0\}\cap J { 0 } ∩ J 、m = 1 m=1 m = 1 、ℓ = 0 \ell=0 ℓ = 0 をもってパラメータ付き Carathéodory 条件を局所的に満たし、N = { 0 } N=\{0\} N = { 0 } として定義 3.2 条件 (a) と定義 3.2 条件 (b) をr = 1 r=1 r = 1 で満たす。t ≠ 0 t\ne0 t = 0 では∂ u f ( t , u , λ ) = 0 \partial_uf(t,u,\lambda)=0 ∂ u f ( t , u , λ ) = 0 であり
∂ λ f ( t , u , λ ) = t 2 ( t 2 + λ 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} ∂ λ f ( t , u , λ ) = ( t 2 + λ 2 ) 3/2 t 2 , ∂ λ f ( t , 0 , 0 ) = ∣ t ∣ 1 である。コンパクト集合K = [ − 1 , 1 ] × { 0 } × { 0 } K=[-1,1]\times\{0\}\times\{0\} K = [ − 1 , 1 ] × { 0 } × { 0 } について、どの零集合N ′ N' N ′ を取っても、t ∈ [ − 1 , 1 ] ∖ N ′ t\in[-1,1]\setminus N' t ∈ [ − 1 , 1 ] ∖ N ′ でm K ( t ) ≥ 1 / ∣ t ∣ m_K(t)\ge1/\lvert t\rvert m K ( t ) ≥ 1/ ∣ t ∣ を満たすm K m_K m K はL 1 ( R ) L^1(\R) L 1 ( R ) に属さない。従ってf f f は定義 3.2 条件 (c) を満たさず、状態とパラメータについてC 1 C^1 C 1 級の Carathéodory 条件を満たさない。f f f とその状態についての偏導関数は支配され、満たされないのはパラメータについての偏導関数の支配だけである。
( t 0 , η , λ ) ∈ R 3 (t_0,\eta,\lambda)\in\R^3 ( t 0 , η , λ ) ∈ R 3 を取る。f f f はu u u によらないから、R \R R 上で定義されたy ( t ) = η + ∫ t 0 t f ( s , 0 , λ ) d s y(t)=\eta+\int_{t_0}^tf(s,0,\lambda)\,ds y ( t ) = η + ∫ t 0 t f ( s , 0 , λ ) d s は、§E13.4 補題 2.4 によりy ′ = f ( t , y , λ ) y'=f(t,y,\lambda) y ′ = f ( t , y , λ ) 、y ( t 0 ) = η y(t_0)=\eta y ( t 0 ) = η のR \R R 上の絶対連続解である。従ってW = R 4 W=\R^4 W = R 4 であり、φ ( t , t 0 , η , λ ) = y ( t ) \varphi(t,t_0,\eta,\lambda)=y(t) φ ( t , t 0 , η , λ ) = y ( t ) である。λ ≠ 0 \lambda\ne0 λ = 0 ではd d s arsinh ( s / ∣ λ ∣ ) = 1 / s 2 + λ 2 \frac{d}{ds}\operatorname{arsinh}(s/\lvert\lambda\rvert)=1/\sqrt{s^2+\lambda^2} d s d arsinh ( s / ∣ λ ∣) = 1/ s 2 + λ 2 であるから
φ ( t , 0 , η , λ ) = η + λ arsinh t ∣ λ ∣ \varphi(t,0,\eta,\lambda)=\eta+\lambda\operatorname{arsinh}\frac{t}{\lvert\lambda\rvert} φ ( t , 0 , η , λ ) = η + λ arsinh ∣ λ ∣ t であり、φ ( t , 0 , η , 0 ) = η \varphi(t,0,\eta,0)=\eta φ ( t , 0 , η , 0 ) = η である。t > 0 t>0 t > 0 を固定すると、λ → 0 \lambda\to0 λ → 0 のとき
φ ( 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 , η , λ ) − φ ( t , 0 , η , 0 ) = arsinh ∣ λ ∣ t ⟶ ∞ であるから、λ ↦ φ ( t , 0 , η , λ ) \lambda\mapsto\varphi(t,0,\eta,\lambda) λ ↦ φ ( t , 0 , η , λ ) はλ = 0 \lambda=0 λ = 0 で微分可能でない。定理 3.8 の結論は、パラメータについての偏導関数の可積分な支配の仮定を外すと成り立たない。一方、∣ λ arsinh ( t / ∣ λ ∣ ) ∣ → 0 \lvert\lambda\operatorname{arsinh}(t/\lvert\lambda\rvert)\rvert\to0 ∣ λ arsinh ( t / ∣ λ ∣)∣ → 0 であり、定理 2.2 (2) のとおりφ \varphi φ はW W W 上で連続である。
例 3.11. Ω = R × R × R \Omega=\R\times\R\times\R Ω = R × R × R 上のf ( t , u , λ ) = u 2 f(t,u,\lambda)=u^2 f ( t , u , λ ) = u 2 は、補題 3.3 (4) によりすべてのr ∈ N ≥ 1 r\in\NN r ∈ N ≥ 1 について状態とパラメータについてC r C^r C r 級の Carathéodory 条件を満たす。η ∈ R \eta\in\R η ∈ R に対して、φ ( t , 0 , η , λ ) = η / ( 1 − η t ) \varphi(t,0,\eta,\lambda)=\eta/(1-\eta t) φ ( t , 0 , η , λ ) = η / ( 1 − η t ) であり、I ( 0 , η , λ ) = { t : η t < 1 } I(0,\eta,\lambda)=\{t:\eta t<1\} I ( 0 , η , λ ) = { t : η t < 1 } である。実際、関数y ( t ) = η / ( 1 − η t ) y(t)=\eta/(1-\eta t) y ( t ) = η / ( 1 − η t ) は{ t : η t < 1 } \{t:\eta t<1\} { t : η t < 1 } 上でy ′ = y 2 y'=y^2 y ′ = y 2 、y ( 0 ) = η y(0)=\eta y ( 0 ) = η を満たし、η ≠ 0 \eta\ne0 η = 0 ではt → 1 / η t\to1/\eta t → 1/ η のとき∣ y ( t ) ∣ → ∞ \lvert y(t)\rvert\to\infty ∣ y ( t )∣ → ∞ であるから、これ以上延長することができない。t > 0 t>0 t > 0 を固定するとW t , 0 = { ( η , λ ) : η < 1 / t } W_{t,0}=\{(\eta,\lambda):\eta<1/t\} W t , 0 = {( η , λ ) : η < 1/ t } であり、W t , 0 W_{t,0} W t , 0 はR 2 \R^2 R 2 全体ではない。
D η φ ( t , 0 , η , λ ) = ( 1 − η t ) − 2 D_\eta\varphi(t,0,\eta,\lambda)=(1-\eta t)^{-2} D η φ ( t , 0 , η , λ ) = ( 1 − η t ) − 2 であり、η t ↑ 1 \eta t\uparrow1 η t ↑ 1 のとき非有界である。y = φ ( ⋅ , 0 , η , λ ) y=\varphi(\cdot,0,\eta,\lambda) y = φ ( ⋅ , 0 , η , λ ) に沿う変分方程式の係数はA ( t ) = 2 y ( t ) = 2 η / ( 1 − η t ) A(t)=2y(t)=2\eta/(1-\eta t) A ( t ) = 2 y ( t ) = 2 η / ( 1 − η t ) 、B ( t ) = 0 B(t)=0 B ( t ) = 0 であり、Y ( t ) = ( 1 − η t ) − 2 Y(t)=(1-\eta t)^{-2} Y ( t ) = ( 1 − η t ) − 2 は
Y ′ ( t ) = 2 η ( 1 − η t ) 3 = A ( t ) Y ( t ) , Y ( 0 ) = 1 Y'(t)=\frac{2\eta}{(1-\eta t)^3}=A(t)Y(t),\qquad Y(0)=1 Y ′ ( t ) = ( 1 − η t ) 3 2 η = A ( t ) Y ( t ) , Y ( 0 ) = 1 を満たす。Z = 0 Z=0 Z = 0 とD λ φ = 0 D_\lambda\varphi=0 D λ φ = 0 も一致し、定理 3.8 の等式が成り立つ。
4 高階の微分可能性と自律系の解写像
定理 4.1. d , p , r ∈ N ≥ 1 d,p,r\in\NN d , p , r ∈ N ≥ 1 とし、開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d が状態とパラメータについてC r C^r C r 級の Carathéodory 条件を満たすとし、φ : W → R d \varphi:W\to\R^d φ : W → R d をその解写像とする。このとき、各( t , t 0 ) ∈ R 2 (t,t_0)\in\R^2 ( t , t 0 ) ∈ R 2 について、( η , λ ) ↦ φ ( t , t 0 , η , λ ) (\boldsymbol\eta,\lambda)\mapsto\varphi(t,t_0,\boldsymbol\eta,\lambda) ( η , λ ) ↦ φ ( t , t 0 , η , λ ) は開集合W t , t 0 = { ( η , λ ) : ( t , t 0 , η , λ ) ∈ W } W_{t,t_0}=\{(\boldsymbol\eta,\lambda):(t,t_0,\boldsymbol\eta,\lambda)\in W\} W t , t 0 = {( η , λ ) : ( t , t 0 , η , λ ) ∈ W } 上でC r C^r C r 級であり、( η , λ ) (\boldsymbol\eta,\lambda) ( η , λ ) の座標についてのr r r 階以下の順序付き偏導関数はすべてW W W 上で連続である。
証明. r = 1 r=1 r = 1 の場合は定理 2.2 (2) と定理 3.8 である。r ≥ 2 r\ge2 r ≥ 2 とし、階数r − 1 r-1 r − 1 の場合の主張がすべてのd , p d,p d , p 、開集合と右辺について成り立つと仮定する。N N N をf f f の零集合とする。
Ω ^ = { ( t , u , U , V , λ ) ∈ R × R d × R d × d × R d × p × R p : ( 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\} Ω = {( t , u , U , V , λ ) ∈ R × R d × R d × d × R d × p × R p : ( t , u , λ ) ∈ Ω } 上で
F ( t , u , U , V , λ ) = ( f ( t , u , λ ) , D u f ( t , u , λ ) U , D u f ( 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) F ( t , u , U , V , λ ) = ( f ( t , u , λ ) , D u f ( t , u , λ ) U , D u f ( t , u , λ ) V + D λ f ( t , u , λ ) ) と定める。状態を( u , U , V ) (\boldsymbol u,U,V) ( u , U , V ) 、パラメータをλ \lambda λ とみると、F F F は状態とパラメータについてC r − 1 C^{r-1} C r − 1 級の Carathéodory 条件を満たす。この確認は演習とする(問題 6.2 )。F F F の解写像をφ ^ : W ^ → R d × R d × d × R d × p \widehat\varphi:\widehat W\to\R^d\times\R^{d\times d}\times\R^{d\times p} φ : W → R d × R d × d × R d × p とする。
( t 0 , η , λ ) ∈ Ω (t_0,\boldsymbol\eta,\lambda)\in\Omega ( t 0 , η , λ ) ∈ Ω 、I = I ( t 0 , η , λ ) I=I(t_0,\boldsymbol\eta,\lambda) I = I ( t 0 , η , λ ) 、y = φ ( ⋅ , t 0 , η , λ ) \boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda) y = φ ( ⋅ , t 0 , η , λ ) とし、Y , Z Y,Z Y , Z をy \boldsymbol y y に沿う変分方程式の解でY ( t 0 ) = I d Y(t_0)=I_d Y ( t 0 ) = I d 、Z ( t 0 ) = 0 Z(t_0)=0 Z ( t 0 ) = 0 を満たすものとする。( y , Y , Z ) (\boldsymbol y,Y,Z) ( y , Y , Z ) はI I I 上で局所絶対連続であり、I I I のほとんど至る所でF F F の方程式を満たし、時刻t 0 t_0 t 0 で( η , I d , 0 ) (\boldsymbol\eta,I_d,0) ( η , I d , 0 ) をとる。§E13.6 定理 2.1 により、( y , Y , Z ) (\boldsymbol y,Y,Z) ( y , Y , Z ) は初期値( η , I d , 0 ) (\boldsymbol\eta,I_d,0) ( η , I d , 0 ) と時刻t 0 t_0 t 0 、パラメータλ \lambda λ に対するF F F の最大絶対連続解y ^ \widehat{\boldsymbol y} y の制限であり、y ^ \widehat{\boldsymbol y} y の最大存在区間I ^ \widehat I I はI I I を含む。y ^ \widehat{\boldsymbol y} y の第一成分はI ^ \widehat I I 上のy ′ = f ( t , y , λ ) \boldsymbol y'=f(t,\boldsymbol y,\lambda) y ′ = f ( t , y , λ ) の絶対連続解で時刻t 0 t_0 t 0 にη \boldsymbol\eta η をとるから、同定理によりy \boldsymbol y y の制限であり、I ^ ⊆ I \widehat I\subseteq I I ⊆ I である。従ってI ^ = I \widehat I=I I = I であり、定理 3.8 により、すべてのt ∈ I t\in I t ∈ I で
φ ^ ( t , t 0 , ( η , I d , 0 ) , λ ) = ( φ ( t , t 0 , η , λ ) , D η φ ( t , t 0 , η , λ ) , D λ φ ( t , t 0 , η , λ ) ) \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 , t 0 , ( η , I d , 0 ) , λ ) = ( φ ( t , t 0 , η , λ ) , D η φ ( t , t 0 , η , λ ) , D λ φ ( t , t 0 , η , λ ) ) である。特に、( t , t 0 , η , λ ) ∈ W (t,t_0,\boldsymbol\eta,\lambda)\in W ( t , t 0 , η , λ ) ∈ W と( t , t 0 , ( η , I d , 0 ) , λ ) ∈ W ^ (t,t_0,(\boldsymbol\eta,I_d,0),\lambda)\in\widehat W ( t , t 0 , ( η , I d , 0 ) , λ ) ∈ W は同値である。
帰納法の仮定をF F F に適用すると、φ ^ \widehat\varphi φ の初期値( η , U , V ) (\boldsymbol\eta,U,V) ( η , U , V ) とパラメータλ \lambda λ の座標についてのr − 1 r-1 r − 1 階以下の順序付き偏導関数は存在してW ^ \widehat W W 上で連続である。U = I d U=I_d U = I d 、V = 0 V=0 V = 0 に固定した写像の( η , λ ) (\boldsymbol\eta,\lambda) ( η , λ ) の座標についての偏微分は、φ ^ \widehat\varphi φ の同じ座標についての偏微分の値を与える。従ってD η φ D_{\boldsymbol\eta}\varphi D η φ とD λ φ D_\lambda\varphi D λ φ の各成分の( η , λ ) (\boldsymbol\eta,\lambda) ( η , λ ) の座標についてのr − 1 r-1 r − 1 階以下の順序付き偏導関数は存在し、連続写像( t , t 0 , η , λ ) ↦ ( t , t 0 , ( η , I d , 0 ) , λ ) (t,t_0,\boldsymbol\eta,\lambda)\mapsto(t,t_0,(\boldsymbol\eta,I_d,0),\lambda) ( t , t 0 , η , λ ) ↦ ( t , t 0 , ( η , I d , 0 ) , λ ) との合成としてW W W 上で連続である。2 ≤ j ≤ r 2\le j\le r 2 ≤ j ≤ r について、φ \varphi φ のj j j 階の順序付き偏導関数は、一階の偏導関数であるD η φ D_{\boldsymbol\eta}\varphi D η φ 、D λ φ D_\lambda\varphi D λ φ の成分のj − 1 j-1 j − 1 階の順序付き偏導関数であるから、存在してW W W 上で連続である。一階の偏導関数は定理 3.8 により、φ \varphi φ 自身は定理 2.2 (2) によりW W W 上で連続である。§E4.4 定理 1.7 により、( η , λ ) ↦ φ ( t , t 0 , η , λ ) (\boldsymbol\eta,\lambda)\mapsto\varphi(t,t_0,\boldsymbol\eta,\lambda) ( η , λ ) ↦ φ ( t , t 0 , η , λ ) はW t , t 0 W_{t,t_0} W t , t 0 上でC r C^r C r 級である。▨
系 4.2. d , p , r ∈ N ≥ 1 d,p,r\in\NN d , p , r ∈ N ≥ 1 とし、開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d と、( u , λ ) (\boldsymbol u,\lambda) ( u , λ ) の座標についてのf f f のr r r 階以下のすべての順序付き偏導関数が、Ω \Omega Ω 上で存在して連続であるとする。補題 3.3 (4) と補題 3.3 (1) によりf f f はパラメータ付き Carathéodory 条件を局所的に満たすので、その解写像をφ : W → R d \varphi:W\to\R^d φ : W → R d とする。このとき、次が成り立つ。
φ ∈ C 1 ( W ) \varphi\in C^1(W) φ ∈ C 1 ( W ) であり、( η , λ ) (\boldsymbol\eta,\lambda) ( η , λ ) の座標についてのr r r 階以下の順序付き偏導関数はW W W 上で連続である。すべての( t , t 0 , η , λ ) ∈ W (t,t_0,\boldsymbol\eta,\lambda)\in W ( t , t 0 , η , λ ) ∈ W で
∂ t φ ( t , t 0 , η , λ ) = f ( t , φ ( t , t 0 , η , λ ) , λ ) , ∂ t 0 φ ( t , t 0 , η , λ ) = − D η φ ( t , t 0 , η , λ ) f ( t 0 , η , λ ) \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) ∂ t φ ( t , t 0 , η , λ ) = f ( t , φ ( t , t 0 , η , λ ) , λ ) , ∂ t 0 φ ( t , t 0 , η , λ ) = − D η φ ( t , t 0 , η , λ ) f ( t 0 , η , λ )
が成り立つ。
f ∈ C r ( Ω ) f\in C^r(\Omega) f ∈ C r ( Ω ) ならばφ ∈ C r ( W ) \varphi\in C^r(W) φ ∈ C r ( W ) である。f ∈ C ∞ ( Ω ) f\in C^\infty(\Omega) f ∈ C ∞ ( Ω ) ならばφ ∈ C ∞ ( W ) \varphi\in C^\infty(W) φ ∈ C ∞ ( W ) である。
証明. 主張 (1) を示す。定理 4.1 により、( η , λ ) (\boldsymbol\eta,\lambda) ( η , λ ) の座標についてのr r r 階以下の順序付き偏導関数はW W W 上で連続である。( t 0 , η , λ ) ∈ Ω (t_0,\boldsymbol\eta,\lambda)\in\Omega ( t 0 , η , λ ) ∈ Ω に対してy = φ ( ⋅ , t 0 , η , λ ) \boldsymbol y=\varphi(\cdot,t_0,\boldsymbol\eta,\lambda) y = φ ( ⋅ , t 0 , η , λ ) とおくと、§E13.4 補題 2.4 によりy ( t ) = η + ∫ t 0 t f ( s , y ( s ) , λ ) d s \boldsymbol y(t)=\boldsymbol\eta+\int_{t_0}^tf(s,\boldsymbol y(s),\lambda)\,ds y ( t ) = η + ∫ t 0 t f ( s , y ( s ) , λ ) d s である。f ( ⋅ , ⋅ , λ ) f(\cdot,\cdot,\lambda) f ( ⋅ , ⋅ , λ ) はΩ λ = { ( t , u ) : ( t , u , λ ) ∈ Ω } \Omega^\lambda=\{(t,\boldsymbol u):(t,\boldsymbol u,\lambda)\in\Omega\} Ω λ = {( t , u ) : ( t , u , λ ) ∈ Ω } 上で連続であるから、§E13.4 補題 2.5 によりy \boldsymbol y y は古典解であり、∂ t φ ( t , t 0 , η , λ ) = f ( t , φ ( t , t 0 , η , λ ) , λ ) \partial_t\varphi(t,t_0,\boldsymbol\eta,\lambda)=f(t,\varphi(t,t_0,\boldsymbol\eta,\lambda),\lambda) ∂ t φ ( t , t 0 , η , λ ) = f ( t , φ ( t , t 0 , η , λ ) , λ ) がすべての点で成り立つ。右辺は定理 2.2 (2) によりW W W 上で連続である。
( t , t 0 , η , λ ) ∈ W (t,t_0,\boldsymbol\eta,\lambda)\in W ( t , t 0 , η , λ ) ∈ W を取る。( t 0 , t 0 , η , λ ) (t_0,t_0,\boldsymbol\eta,\lambda) ( t 0 , t 0 , η , λ ) と( t , t 0 , η , λ ) (t,t_0,\boldsymbol\eta,\lambda) ( t , t 0 , η , λ ) は開集合W W W に属するから、あるh 1 > 0 h_1>0 h 1 > 0 が存在し、∣ h ∣ < h 1 \lvert h\rvert<h_1 ∣ h ∣ < h 1 ならば( t 0 , t 0 + h , η , λ ) ∈ W (t_0,t_0+h,\boldsymbol\eta,\lambda)\in W ( t 0 , t 0 + h , η , λ ) ∈ W 、( t , t 0 + h , η , λ ) ∈ W (t,t_0+h,\boldsymbol\eta,\lambda)\in W ( t , t 0 + h , η , λ ) ∈ W である。このときη h = φ ( t 0 , t 0 + h , η , λ ) \boldsymbol\eta_h=\varphi(t_0,t_0+h,\boldsymbol\eta,\lambda) η h = φ ( t 0 , t 0 + h , η , λ ) とおく。φ ( ⋅ , t 0 + h , η , λ ) \varphi(\cdot,t_0+h,\boldsymbol\eta,\lambda) φ ( ⋅ , t 0 + h , η , λ ) とφ ( ⋅ , t 0 , η h , λ ) \varphi(\cdot,t_0,\boldsymbol\eta_h,\lambda) φ ( ⋅ , t 0 , η h , λ ) はともに、時刻t 0 t_0 t 0 でη h \boldsymbol\eta_h η h をとる最大絶対連続解であるから、§E13.6 定理 2.1 の一意性により一致する。従って
φ ( t , t 0 + h , η , λ ) = φ ( t , t 0 , η h , λ ) \varphi(t,t_0+h,\boldsymbol\eta,\lambda)=\varphi(t,t_0,\boldsymbol\eta_h,\lambda) φ ( t , t 0 + h , η , λ ) = φ ( t , t 0 , η h , λ ) である。古典解の積分表示により
η h − η h = − 1 h ∫ t 0 t 0 + h f ( s , φ ( s , t 0 + h , η , λ ) , λ ) d s \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 h η h − η = − h 1 ∫ t 0 t 0 + h f ( s , φ ( s , t 0 + h , η , λ ) , λ ) d s であり、( s , h ) ↦ f ( s , φ ( s , t 0 + h , η , λ ) , λ ) (s,h)\mapsto f(s,\varphi(s,t_0+h,\boldsymbol\eta,\lambda),\lambda) ( s , h ) ↦ f ( s , φ ( s , t 0 + h , η , λ ) , λ ) は定理 2.2 (2) により( t 0 , 0 ) (t_0,0) ( t 0 , 0 ) で連続であるから、h → 0 h\to0 h → 0 のとき右辺は− f ( t 0 , η , λ ) -f(t_0,\boldsymbol\eta,\lambda) − f ( t 0 , η , λ ) へ収束する。η ′ ↦ φ ( t , t 0 , η ′ , λ ) \boldsymbol\eta'\mapsto\varphi(t,t_0,\boldsymbol\eta',\lambda) η ′ ↦ φ ( t , t 0 , η ′ , λ ) は定理 3.8 によりη \boldsymbol\eta η の近傍でC 1 C^1 C 1 級であるから、§E4.3 定理 1.1 により
∂ t 0 φ ( t , t 0 , η , λ ) = d d h ∣ h = 0 φ ( t , t 0 , η h , λ ) = − D η φ ( t , t 0 , η , λ ) f ( t 0 , η , λ ) \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) ∂ t 0 φ ( t , t 0 , η , λ ) = d h d h = 0 φ ( t , t 0 , η h , λ ) = − D η φ ( t , t 0 , η , λ ) f ( t 0 , η , λ ) である。右辺は定理 3.8 によりW W W 上で連続である。φ \varphi φ の( t , t 0 , η , λ ) (t,t_0,\boldsymbol\eta,\lambda) ( t , t 0 , η , λ ) の座標についての一階の偏導関数はすべて存在してW W W 上で連続であるから、§E4.2 系 2.2 によりφ ∈ C 1 ( W ) \varphi\in C^1(W) φ ∈ C 1 ( W ) である。
主張 (2) を示す。f ∈ C r ( Ω ) f\in C^r(\Omega) f ∈ C r ( Ω ) とする。
Ω ~ = { ( σ , w , τ , s , λ ) ∈ R × R d × R × R × R p : ( τ , 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 , τ , s , λ ) ∈ R × R d × R × R × R p : ( τ , w , λ ) ∈ Ω } 上で、状態を( w , τ ) (\boldsymbol w,\tau) ( w , τ ) 、パラメータを( s , λ ) (s,\lambda) ( s , λ ) とする写像
H ( σ , w , τ , s , λ ) = ( s f ( τ , w , λ ) , s ) H(\sigma,\boldsymbol w,\tau,s,\lambda)=\bigl(s\,f(\tau,\boldsymbol w,\lambda),\ s\bigr) H ( σ , w , τ , s , λ ) = ( s f ( τ , w , λ ) , s ) を考える。H H H はσ \sigma σ によらず、s s s について一次であるから、s s s について二回以上偏微分したH H H の順序付き偏導関数は零である。従って、( w , τ , s , λ ) (\boldsymbol w,\tau,s,\lambda) ( w , τ , s , λ ) の座標についてのH H H の順序付き偏導関数は、f f f の順序付き偏導関数にs s s を掛けたもの、f f f の順序付き偏導関数、定数、零のいずれかである。従ってr r r 階以下の順序付き偏導関数はΩ ~ \widetilde\Omega Ω 上で連続であり、補題 3.3 (4) によりH H H は状態とパラメータについてC r C^r C r 級の Carathéodory 条件を満たす。H H H の解写像をχ : W ~ → R d × R \chi:\widetilde W\to\R^d\times\R χ : W → R d × R とする。
( t , t 0 , η , λ ) ∈ W (t,t_0,\boldsymbol\eta,\lambda)\in W ( t , t 0 , η , λ ) ∈ W を取り、s = t − t 0 s=t-t_0 s = t − t 0 とおく。T = { σ ∈ R : t 0 + s σ ∈ I ( t 0 , η , λ ) } T=\{\sigma\in\R:t_0+s\sigma\in I(t_0,\boldsymbol\eta,\lambda)\} T = { σ ∈ R : t 0 + s σ ∈ I ( t 0 , η , λ )} は0 0 0 と1 1 1 を含む開区間であり、T T T 上で
w ^ ( σ ) = ( φ ( t 0 + s σ , t 0 , η , λ ) , t 0 + s σ ) \widehat{\boldsymbol w}(\sigma)=\bigl(\varphi(t_0+s\sigma,t_0,\boldsymbol\eta,\lambda),\ t_0+s\sigma\bigr) w ( σ ) = ( φ ( t 0 + s σ , t 0 , η , λ ) , t 0 + s σ ) とおく。主張 (1) の∂ t φ \partial_t\varphi ∂ t φ の式により、w ^ \widehat{\boldsymbol w} w はC 1 C^1 C 1 級でw ^ ′ = H ( σ , w ^ , s , λ ) \widehat{\boldsymbol w}'=H(\sigma,\widehat{\boldsymbol w},s,\lambda) w ′ = H ( σ , w , s , λ ) 、w ^ ( 0 ) = ( η , t 0 ) \widehat{\boldsymbol w}(0)=(\boldsymbol\eta,t_0) w ( 0 ) = ( η , t 0 ) を満たす。§E13.6 定理 2.1 によりw ^ \widehat{\boldsymbol w} w はH H H の最大絶対連続解の制限であるから、( 1 , 0 , ( η , t 0 ) , ( s , λ ) ) ∈ W ~ (1,0,(\boldsymbol\eta,t_0),(s,\lambda))\in\widetilde W ( 1 , 0 , ( η , t 0 ) , ( s , λ )) ∈ W であり
χ ( 1 , 0 , ( η , t 0 ) , ( t − t 0 , λ ) ) = ( φ ( t , t 0 , η , λ ) , t ) \chi\bigl(1,0,(\boldsymbol\eta,t_0),(t-t_0,\lambda)\bigr)=\bigl(\varphi(t,t_0,\boldsymbol\eta,\lambda),\ t\bigr) χ ( 1 , 0 , ( η , t 0 ) , ( t − t 0 , λ ) ) = ( φ ( t , t 0 , η , λ ) , t ) である。線形写像Λ ( t , t 0 , η , λ ) = ( ( η , t 0 ) , ( t − t 0 , λ ) ) \Lambda(t,t_0,\boldsymbol\eta,\lambda)=((\boldsymbol\eta,t_0),(t-t_0,\lambda)) Λ ( t , t 0 , η , λ ) = (( η , t 0 ) , ( t − t 0 , λ )) はW W W を開集合W ~ 1 , 0 = { ( η ^ , μ ) : ( 1 , 0 , η ^ , μ ) ∈ W ~ } \widetilde W_{1,0}=\{(\widehat{\boldsymbol\eta},\mu):(1,0,\widehat{\boldsymbol\eta},\mu)\in\widetilde W\} W 1 , 0 = {( η , μ ) : ( 1 , 0 , η , μ ) ∈ W } へ写し、φ \varphi φ はχ ( 1 , 0 , ⋅ ) ∘ Λ \chi(1,0,\cdot)\circ\Lambda χ ( 1 , 0 , ⋅ ) ∘ Λ の第一成分である。定理 4.1 によりχ ( 1 , 0 , ⋅ ) \chi(1,0,\cdot) χ ( 1 , 0 , ⋅ ) はW ~ 1 , 0 \widetilde W_{1,0} W 1 , 0 上でC r C^r C r 級である。Λ \Lambda Λ は線形であるから、§E4.3 定理 1.1 を反復して適用すると、φ \varphi φ のr r r 階以下の順序付き偏導関数はχ ( 1 , 0 , ⋅ ) \chi(1,0,\cdot) χ ( 1 , 0 , ⋅ ) の同じ階数の順序付き偏導関数とΛ \Lambda Λ の合成の定数係数の一次結合であり、W W W 上で連続である。§E4.4 定理 1.7 によりφ ∈ C r ( W ) \varphi\in C^r(W) φ ∈ C r ( W ) である。f ∈ C ∞ ( Ω ) f\in C^\infty(\Omega) f ∈ C ∞ ( Ω ) ならば、各r ∈ N ≥ 1 r\in\NN r ∈ N ≥ 1 に適用してφ ∈ C ∞ ( W ) \varphi\in C^\infty(W) φ ∈ C ∞ ( W ) を得る。▨
系 4.4. n ∈ N ≥ 1 n\in\NN n ∈ N ≥ 1 とし、U ⊆ R n U\subseteq\R^n U ⊆ R n を開集合、F : U → R n F:U\to\R^n F : U → R n を局所 Lipschitz 連続な写像とする。ξ ∈ U \boldsymbol\xi\in U ξ ∈ U に対して、x ′ = F ( x ) \boldsymbol x'=F(\boldsymbol x) x ′ = F ( x ) 、x ( 0 ) = ξ \boldsymbol x(0)=\boldsymbol\xi x ( 0 ) = ξ の最大解をφ ( ⋅ , ξ ) \varphi(\cdot,\boldsymbol\xi) φ ( ⋅ , ξ ) 、その最大存在区間をI ( ξ ) I(\boldsymbol\xi) I ( ξ ) とし、D = { ( t , ξ ) : ξ ∈ U , t ∈ I ( ξ ) } D=\{(t,\boldsymbol\xi):\boldsymbol\xi\in U,\ t\in I(\boldsymbol\xi)\} D = {( t , ξ ) : ξ ∈ U , t ∈ I ( ξ )} とおく。このとき、次が成り立つ。
D D D はR × R n \R\times\R^n R × R n の開集合であり、( t , ξ ) ↦ φ ( t , ξ ) (t,\boldsymbol\xi)\mapsto\varphi(t,\boldsymbol\xi) ( t , ξ ) ↦ φ ( t , ξ ) はD D D 上で連続である。
r ∈ N ≥ 1 r\in\NN r ∈ N ≥ 1 とし、F ∈ C r ( U ) F\in C^r(U) F ∈ C r ( U ) ならば、φ ∈ C r ( D ) \varphi\in C^r(D) φ ∈ C r ( D ) であり、∂ t φ ( t , ξ ) = F ( φ ( t , ξ ) ) \partial_t\varphi(t,\boldsymbol\xi)=F(\varphi(t,\boldsymbol\xi)) ∂ t φ ( t , ξ ) = F ( φ ( t , ξ )) である。さらに、各ξ ∈ U \boldsymbol\xi\in U ξ ∈ U についてt ↦ D ξ φ ( t , ξ ) t\mapsto D_{\boldsymbol\xi}\varphi(t,\boldsymbol\xi) t ↦ D ξ φ ( t , ξ ) はY ′ = D F ( φ ( t , ξ ) ) Y Y'=DF(\varphi(t,\boldsymbol\xi))Y Y ′ = D F ( φ ( t , ξ )) Y 、Y ( 0 ) = I n Y(0)=I_n Y ( 0 ) = I n のI ( ξ ) I(\boldsymbol\xi) I ( ξ ) 上の解である。F ∈ C ∞ ( U ) F\in C^\infty(U) F ∈ C ∞ ( U ) ならばφ ∈ C ∞ ( D ) \varphi\in C^\infty(D) φ ∈ C ∞ ( D ) である。
証明. Ω = R × U × R \Omega=\R\times U\times\R Ω = R × U × R 上でf ( t , u , λ ) = F ( u ) f(t,\boldsymbol u,\lambda)=F(\boldsymbol u) f ( t , u , λ ) = F ( u ) とおく(d = n d=n d = n 、p = 1 p=1 p = 1 )。Ω \Omega Ω に含まれる閉矩形J × B × C J\times B\times C J × B × C についてJ × B ⊆ R × U J\times B\subseteq\R\times U J × B ⊆ R × U であり、§E13.6 補題 5.1 を( t , u ) ↦ F ( u ) (t,\boldsymbol u)\mapsto F(\boldsymbol u) ( t , u ) ↦ F ( u ) に適用すると、J × B J\times B J × B 上の零集合、支配関数と Lipschitz 係数が得られる。これらはλ \lambda λ によらないから、f f f はパラメータ付き Carathéodory 条件を局所的に満たす。f f f の解写像をφ f : W → R n \varphi_f:W\to\R^n φ f : W → R n とする。同補題によりΩ 0 = R × U \Omega^0=\R\times U Ω 0 = R × U 上の最大絶対連続解は最大解であるから、φ ( t , ξ ) = φ f ( t , 0 , ξ , 0 ) \varphi(t,\boldsymbol\xi)=\varphi_f(t,0,\boldsymbol\xi,0) φ ( t , ξ ) = φ f ( t , 0 , ξ , 0 ) であり、D D D は連続写像ι ( t , ξ ) = ( t , 0 , ξ , 0 ) \iota(t,\boldsymbol\xi)=(t,0,\boldsymbol\xi,0) ι ( t , ξ ) = ( t , 0 , ξ , 0 ) によるW W W の逆像である。
主張 (1) を示す。定理 2.2 (2) によりW W W は開集合でφ f \varphi_f φ f は連続であるから、D = ι − 1 ( W ) D=\iota^{-1}(W) D = ι − 1 ( W ) は開集合であり、φ = φ f ∘ ι \varphi=\varphi_f\circ\iota φ = φ f ∘ ι は連続である。
主張 (2) を示す。F ∈ C r ( U ) F\in C^r(U) F ∈ C r ( U ) ならばf ∈ C r ( Ω ) f\in C^r(\Omega) f ∈ C r ( Ω ) であり、系 4.2 (2) によりφ f ∈ C r ( W ) \varphi_f\in C^r(W) φ f ∈ C r ( W ) である。ι \iota ι は座標の一部を定数に固定する写像であるから、φ = φ f ∘ ι \varphi=\varphi_f\circ\iota φ = φ f ∘ ι の順序付き偏導関数はφ f \varphi_f φ f の対応する順序付き偏導関数とι \iota ι の合成であり、§E4.4 定理 1.7 によりφ ∈ C r ( D ) \varphi\in C^r(D) φ ∈ 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)) ∂ t φ ( t , ξ ) = f ( t , φ ( t , ξ ) , 0 ) = F ( φ ( t , ξ )) である。定理 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) D ξ φ ( t , ξ ) = D η φ f ( t , 0 , ξ , 0 ) は、φ ( ⋅ , ξ ) \varphi(\cdot,\boldsymbol\xi) φ ( ⋅ , ξ ) に沿う変分方程式Y ′ = D u f ( t , φ ( t , ξ ) , 0 ) Y = D F ( φ ( t , ξ ) ) Y Y'=D_{\boldsymbol u}f(t,\varphi(t,\boldsymbol\xi),0)Y=DF(\varphi(t,\boldsymbol\xi))Y Y ′ = D u f ( t , φ ( t , ξ ) , 0 ) Y = D F ( φ ( t , ξ )) Y の解でY ( 0 ) = I n Y(0)=I_n Y ( 0 ) = I n を満たすものである。F ∈ C ∞ ( U ) F\in C^\infty(U) F ∈ C ∞ ( U ) ならば、各r r r に適用してφ ∈ C ∞ ( D ) \varphi\in C^\infty(D) φ ∈ C ∞ ( D ) を得る。▨
5 線形系の複素パラメータについての正則性
命題 5.1. d ∈ N ≥ 1 d\in\NN d ∈ N ≥ 1 とし、C d \C^d C d にノルム∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ を固定し、d × d d\times d d × d 複素行列にはその作用素ノルムを用いる。I ⊆ R I\subseteq\R I ⊆ R を内部が空でない区間とし、A 0 , A 1 ∈ L l o c 1 ( I ; C d × d ) A_0,A_1\in L^1_{\mathrm{loc}}(I;\C^{d\times d}) A 0 , A 1 ∈ L loc 1 ( I ; C d × d ) 、b 0 , b 1 ∈ L l o c 1 ( I ; C d ) \boldsymbol b_0,\boldsymbol b_1\in L^1_{\mathrm{loc}}(I;\C^d) b 0 , b 1 ∈ L loc 1 ( I ; C d ) 、t 0 ∈ I t_0\in I t 0 ∈ I 、η ∈ C d \boldsymbol\eta\in\C^d η ∈ C d とする。このとき、各λ ∈ C \lambda\in\C λ ∈ C について、I I I のほとんど至る所で
x ′ = ( A 0 ( t ) + λ A 1 ( t ) ) x + b 0 ( t ) + λ b 1 ( t ) \boldsymbol x'=\bigl(A_0(t)+\lambda A_1(t)\bigr)\boldsymbol x+\boldsymbol b_0(t)+\lambda\boldsymbol b_1(t) x ′ = ( A 0 ( t ) + λ A 1 ( t ) ) x + b 0 ( t ) + λ b 1 ( t ) を満たしx ( t 0 ) = η \boldsymbol x(t_0)=\boldsymbol\eta x ( t 0 ) = η を満たすx ( ⋅ , λ ) ∈ A C l o c ( I ; C d ) \boldsymbol x(\cdot,\lambda)\in AC_{\mathrm{loc}}(I;\C^d) x ( ⋅ , λ ) ∈ A C loc ( I ; C d ) が一意に存在する。Picard 反復を、t ∈ I t\in I t ∈ I 、λ ∈ C \lambda\in\C λ ∈ C に対してx 0 ( t , λ ) = η \boldsymbol x_0(t,\lambda)=\boldsymbol\eta x 0 ( t , λ ) = η と
x k + 1 ( t , λ ) = η + ∫ t 0 t ( ( A 0 ( s ) + λ A 1 ( s ) ) x k ( s , λ ) + b 0 ( s ) + λ b 1 ( s ) ) d s \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 x k + 1 ( t , λ ) = η + ∫ t 0 t ( ( A 0 ( s ) + λ A 1 ( s ) ) x k ( s , λ ) + b 0 ( s ) + λ b 1 ( s ) ) d s によって定める。このとき、次が成り立つ。
( t , λ ) ↦ x ( t , λ ) (t,\lambda)\mapsto\boldsymbol x(t,\lambda) ( t , λ ) ↦ x ( t , λ ) はI × C I\times\C I × C 上で連続である。
各t ∈ I t\in I t ∈ I 、λ ∈ C \lambda\in\C λ ∈ 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 ∂ λ x ( t , λ ) = lim h → 0 ( x ( t , λ + h ) − x ( t , λ )) / h (h ∈ C ∖ { 0 } h\in\C\setminus\{0\} h ∈ C ∖ { 0 } )が存在する。各λ \lambda λ について、t ↦ ∂ λ x ( t , λ ) t\mapsto\partial_\lambda\boldsymbol x(t,\lambda) t ↦ ∂ λ x ( t , λ ) は
z ′ = ( A 0 ( t ) + λ A 1 ( t ) ) z + A 1 ( t ) x ( t , λ ) + b 1 ( t ) , z ( t 0 ) = 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 z ′ = ( A 0 ( t ) + λ A 1 ( t ) ) z + A 1 ( t ) x ( t , λ ) + b 1 ( t ) , z ( t 0 ) = 0
のI I I 上の局所絶対連続な解である。
各t ∈ I t\in I t ∈ I について、λ ↦ x ( t , λ ) \lambda\mapsto\boldsymbol x(t,\lambda) λ ↦ x ( t , λ ) は整関数である。
I I I に含まれる任意の有限閉区間J J J とR > 0 R>0 R > 0 について、x k \boldsymbol x_k x k はJ × { λ ∈ C : ∣ λ ∣ ≤ R } J\times\{\lambda\in\C:\lvert\lambda\rvert\le R\} J × { λ ∈ C : ∣ λ ∣ ≤ R } 上でx \boldsymbol x x へ一様収束する。
A 0 , A 1 , b 0 , b 1 A_0,A_1,\boldsymbol b_0,\boldsymbol b_1 A 0 , A 1 , b 0 , b 1 が実数値でありη ∈ R d \boldsymbol\eta\in\R^d η ∈ R d ならば、すべてのt ∈ I t\in I t ∈ I とλ ∈ R \lambda\in\R λ ∈ R についてx ( t , λ ) ∈ R d \boldsymbol x(t,\lambda)\in\R^d x ( t , λ ) ∈ R d である。
証明. 補題 3.1 (3) をK = C \mathbb K=\C K = C 、係数A 0 + λ A 1 A_0+\lambda A_1 A 0 + λ A 1 、非斉次項b 0 + λ b 1 \boldsymbol b_0+\lambda\boldsymbol b_1 b 0 + λ b 1 に適用すると、x ( ⋅ , λ ) \boldsymbol x(\cdot,\lambda) x ( ⋅ , λ ) は一意に存在する。t 0 t_0 t 0 を含みI I I に含まれる長さが正の有限閉区間J J J 上で、x ( ⋅ , λ ) ∣ J \boldsymbol x(\cdot,\lambda)|_J x ( ⋅ , λ ) ∣ J は補題 3.1 (1) の解である。
主張 (1) を示す。( t , λ ) ∈ I × C (t,\lambda)\in I\times\C ( t , λ ) ∈ I × C を取り、t 0 t_0 t 0 を含み、I I I におけるt t t の近傍を含むJ J J を取る。λ n → λ \lambda_n\to\lambda λ n → λ ならば
∫ J ∥ ( A 0 + λ n A 1 ) − ( A 0 + λ A 1 ) ∥ = ∣ λ n − λ ∣ ∫ J ∥ A 1 ∥ → 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 ∫ J ( A 0 + λ n A 1 ) − ( A 0 + λ A 1 ) = ∣ λ n − λ ∣ ∫ J ∥ A 1 ∥ → 0 であり、非斉次項についても同様であるから、補題 3.1 (2) によりx ( ⋅ , λ n ) → x ( ⋅ , λ ) \boldsymbol x(\cdot,\lambda_n)\to\boldsymbol x(\cdot,\lambda) x ( ⋅ , λ n ) → x ( ⋅ , λ ) はJ J J 上で一様収束する。( t n , λ n ) → ( t , λ ) (t_n,\lambda_n)\to(t,\lambda) ( t n , λ n ) → ( t , λ ) ならば、十分大きいn n n でt n ∈ J t_n\in J t n ∈ J であり、∥ x ( t n , λ n ) − x ( t , λ ) ∥ ≤ max J ∥ x ( ⋅ , λ n ) − x ( ⋅ , λ ) ∥ + ∥ x ( t n , λ ) − 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 ∥ x ( t n , λ n ) − x ( t , λ )∥ ≤ max J ∥ x ( ⋅ , λ n ) − x ( ⋅ , λ )∥ + ∥ x ( t n , λ ) − x ( t , λ )∥ → 0 である。
主張 (2) を示す。λ ∈ C \lambda\in\C λ ∈ C を固定し、z \boldsymbol z z を主張の方程式のI I I 上の解とする(A 1 x ( ⋅ , λ ) + b 1 ∈ L l o c 1 A_1\boldsymbol x(\cdot,\lambda)+\boldsymbol b_1\in L^1_{\mathrm{loc}} A 1 x ( ⋅ , λ ) + b 1 ∈ L loc 1 であるから補題 3.1 (3) により存在する)。h ∈ C ∖ { 0 } h\in\C\setminus\{0\} h ∈ C ∖ { 0 } に対してq h = ( x ( ⋅ , λ + h ) − x ( ⋅ , λ ) ) / h \boldsymbol q_h=(\boldsymbol x(\cdot,\lambda+h)-\boldsymbol x(\cdot,\lambda))/h q h = ( x ( ⋅ , λ + h ) − x ( ⋅ , λ )) / h とおくと、q h ( t 0 ) = 0 \boldsymbol q_h(t_0)=\boldsymbol 0 q h ( t 0 ) = 0 であり、ほとんど至る所で
q h ′ = ( A 0 + λ A 1 ) q h + A 1 x ( ⋅ , λ + h ) + b 1 \boldsymbol q_h'=\bigl(A_0+\lambda A_1\bigr)\boldsymbol q_h+A_1\boldsymbol x(\cdot,\lambda+h)+\boldsymbol b_1 q h ′ = ( A 0 + λ A 1 ) q h + A 1 x ( ⋅ , λ + h ) + b 1 である。J J J をt 0 t_0 t 0 を含みI I I に含まれる長さが正の有限閉区間とする。h → 0 h\to0 h → 0 のとき、主張 (1) の証明によりmax J ∥ x ( ⋅ , λ + h ) − x ( ⋅ , λ ) ∥ → 0 \max_J\lVert\boldsymbol x(\cdot,\lambda+h)-\boldsymbol x(\cdot,\lambda)\rVert\to0 max J ∥ x ( ⋅ , λ + h ) − x ( ⋅ , λ )∥ → 0 であるから、∫ J ∥ A 1 ( x ( ⋅ , λ + h ) − x ( ⋅ , λ ) ) ∥ → 0 \int_J\lVert A_1(\boldsymbol x(\cdot,\lambda+h)-\boldsymbol x(\cdot,\lambda))\rVert\to0 ∫ J ∥ A 1 ( x ( ⋅ , λ + h ) − x ( ⋅ , λ ))∥ → 0 である。零へ収束する任意の点列h n h_n h n に補題 3.1 (2) を適用すると、q h n → z \boldsymbol q_{h_n}\to\boldsymbol z q h n → z はJ J J 上で一様収束する。J J J と点列は任意であるから、各t ∈ I t\in I t ∈ I で∂ λ x ( t , λ ) \partial_\lambda\boldsymbol x(t,\lambda) ∂ λ x ( t , λ ) が存在してz ( t ) \boldsymbol z(t) z ( t ) に等しい。
主張 (3) を示す。主張 (2) により、各t t t についてλ ↦ x ( t , λ ) \lambda\mapsto\boldsymbol x(t,\lambda) λ ↦ x ( t , λ ) の各成分はC \C C の各点で複素微分可能である。
主張 (4) を示す。J J J とt 0 t_0 t 0 を含むI I I の最小の閉区間をJ ′ J' J ′ とする。帰納法により、各x k ( ⋅ , λ ) \boldsymbol x_k(\cdot,\lambda) x k ( ⋅ , λ ) はI I I 上で連続であり、Picard 反復の被積分関数は局所可積分である。J ′ J' J ′ 上でα = ∥ A 0 ∥ + R ∥ A 1 ∥ \alpha=\lVert A_0\rVert+R\lVert A_1\rVert α = ∥ A 0 ∥ + R ∥ A 1 ∥ 、β = ∥ b 0 ∥ + R ∥ b 1 ∥ \beta=\lVert\boldsymbol b_0\rVert+R\lVert\boldsymbol b_1\rVert β = ∥ b 0 ∥ + R ∥ b 1 ∥ 、H ( t ) = ∣ ∫ t 0 t α ∣ H(t)=\lvert\int_{t_0}^t\alpha\rvert H ( t ) = ∣ ∫ t 0 t α ∣ 、C 0 = ∥ η ∥ ∫ J ′ α + ∫ J ′ β C_0=\lVert\boldsymbol\eta\rVert\int_{J'}\alpha+\int_{J'}\beta C 0 = ∥ η ∥ ∫ J ′ α + ∫ J ′ β とおく。t ∈ J ′ t\in J' t ∈ J ′ 、∣ λ ∣ ≤ R \lvert\lambda\rvert\le R ∣ λ ∣ ≤ R 、k ≥ 0 k\ge0 k ≥ 0 について
∥ x k + 1 ( t , λ ) − x k ( t , λ ) ∥ ≤ C 0 H ( t ) k k ! \lVert\boldsymbol x_{k+1}(t,\lambda)-\boldsymbol x_k(t,\lambda)\rVert\le C_0\frac{H(t)^k}{k!} ∥ x k + 1 ( t , λ ) − x k ( t , λ )∥ ≤ C 0 k ! H ( t ) k が成り立つ。実際、k = 0 k=0 k = 0 では、§E13.2 補題 1.3 (2) をC d \C^d C d の固定したノルムに適用すると、左辺は∣ ∫ t 0 t ( α ∥ η ∥ + β ) ∣ ≤ C 0 \lvert\int_{t_0}^t(\alpha\lVert\boldsymbol\eta\rVert+\beta)\rvert\le C_0 ∣ ∫ t 0 t ( α ∥ η ∥ + β )∣ ≤ C 0 以下である。k k k で成り立つとすると、x k + 2 − x k + 1 = ∫ t 0 t ( A 0 + λ A 1 ) ( x k + 1 − x k ) d s \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 x k + 2 − x k + 1 = ∫ t 0 t ( A 0 + λ A 1 ) ( x k + 1 − x k ) d s と同じ評価により、左辺はC 0 ∣ ∫ t 0 t α ( s ) H ( s ) k / k ! d s ∣ C_0\lvert\int_{t_0}^t\alpha(s)H(s)^k/k!\,ds\rvert C 0 ∣ ∫ t 0 t α ( s ) H ( s ) k / k ! d s ∣ 以下である。t ≥ t 0 t\ge t_0 t ≥ t 0 の側ではH H H は絶対連続でほとんど至る所H ′ = α H'=\alpha H ′ = α であり、§E13.2 定理 3.1 (1) と§E13.2 定理 2.5 により∫ t 0 t α H k / k ! = H ( t ) k + 1 / ( k + 1 ) ! \int_{t_0}^t\alpha H^k/k!=H(t)^{k+1}/(k+1)! ∫ t 0 t α H k / k ! = H ( t ) k + 1 / ( k + 1 )! である。t ≤ t 0 t\le t_0 t ≤ t 0 の側も同様である。H ∗ = ∫ J ′ α H_*=\int_{J'}\alpha H ∗ = ∫ J ′ α とおくと∑ k C 0 H ∗ k / k ! < ∞ \sum_kC_0H_*^k/k!<\infty ∑ k C 0 H ∗ k / k ! < ∞ であるから、( x k ) (\boldsymbol x_k) ( x k ) はJ ′ × { ∣ λ ∣ ≤ R } J'\times\{\lvert\lambda\rvert\le R\} J ′ × {∣ λ ∣ ≤ R } 上で一様に Cauchy 列であり、一様収束の極限x ~ \widetilde{\boldsymbol x} x をもつ。∣ λ ∣ ≤ R \lvert\lambda\rvert\le R ∣ λ ∣ ≤ R について
∥ ∫ t 0 t ( A 0 + λ A 1 ) ( x k ( s , λ ) − x ~ ( s , λ ) ) d s ∥ ≤ H ∗ max J ′ ∥ x k ( ⋅ , λ ) − 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 ∫ t 0 t ( A 0 + λ A 1 ) ( x k ( s , λ ) − x ( s , λ ) ) d s ≤ H ∗ J ′ max ∥ x k ( ⋅ , λ ) − x ( ⋅ , λ )∥ → 0 が§E13.2 補題 1.3 (2) により成り立つから、Picard 反復の定義式でk → ∞ k\to\infty k → ∞ とすると、連続関数x ~ ( ⋅ , λ ) \widetilde{\boldsymbol x}(\cdot,\lambda) x ( ⋅ , λ ) はすべてのt ∈ J ′ t\in J' t ∈ J ′ で
x ~ ( t , λ ) = η + ∫ t 0 t ( ( A 0 ( s ) + λ A 1 ( s ) ) x ~ ( s , λ ) + b 0 ( s ) + λ b 1 ( s ) ) d s \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 x ( t , λ ) = η + ∫ t 0 t ( ( A 0 ( s ) + λ A 1 ( s ) ) x ( s , λ ) + b 0 ( s ) + λ b 1 ( s ) ) d s を満たす。§E13.2 系 2.6 によりx ~ ( ⋅ , λ ) \widetilde{\boldsymbol x}(\cdot,\lambda) x ( ⋅ , λ ) はJ ′ J' J ′ 上で絶対連続であり、ほとんど至る所で方程式を満たし、x ~ ( t 0 , λ ) = η \widetilde{\boldsymbol x}(t_0,\lambda)=\boldsymbol\eta x ( t 0 , λ ) = η である。補題 3.1 (1) の一意性によりx ~ ( ⋅ , λ ) = x ( ⋅ , λ ) ∣ J ′ \widetilde{\boldsymbol x}(\cdot,\lambda)=\boldsymbol x(\cdot,\lambda)|_{J'} x ( ⋅ , λ ) = x ( ⋅ , λ ) ∣ J ′ である。
主張 (5) を示す。A 0 , A 1 , b 0 , b 1 A_0,A_1,\boldsymbol b_0,\boldsymbol b_1 A 0 , A 1 , b 0 , b 1 が実数値であり、η ∈ R d \boldsymbol\eta\in\R^d η ∈ R d 、λ ∈ R \lambda\in\R λ ∈ R であるとする。補題 3.1 (3) をK = R \mathbb K=\R K = R 、R d \R^d R d の任意のノルム、実係数A 0 + λ A 1 A_0+\lambda A_1 A 0 + λ A 1 、b 0 + λ b 1 \boldsymbol b_0+\lambda\boldsymbol b_1 b 0 + λ b 1 に適用すると、I I I のほとんど至る所で方程式を満たし時刻t 0 t_0 t 0 にη \boldsymbol\eta η をとるx ~ ∈ A C l o c ( I ; R d ) \widetilde{\boldsymbol x}\in AC_{\mathrm{loc}}(I;\R^d) x ∈ A C loc ( I ; R d ) が存在する。x ~ \widetilde{\boldsymbol x} x はC d \C^d C d に値をとる同じ初期値問題の局所絶対連続な解でもあるから、x ( ⋅ , λ ) \boldsymbol x(\cdot,\lambda) x ( ⋅ , λ ) の一意性によりx ( ⋅ , λ ) = x ~ \boldsymbol x(\cdot,\lambda)=\widetilde{\boldsymbol x} x ( ⋅ , λ ) = x であり、すべてのt ∈ I t\in I t ∈ I でx ( t , λ ) ∈ R d \boldsymbol x(t,\lambda)\in\R^d x ( t , λ ) ∈ R d である。▨
例 5.2. λ ∈ C \lambda\in\C λ ∈ C に対して、初期値問題y ′ ′ + λ y = 0 y''+\lambda y=0 y ′′ + λ y = 0 、y ( 0 ) = 1 y(0)=1 y ( 0 ) = 1 、y ′ ( 0 ) = 0 y'(0)=0 y ′ ( 0 ) = 0 を、x = ( y , y ′ ) \boldsymbol x=(y,y') x = ( y , y ′ ) についての一階の系
x ′ = ( A 0 + λ A 1 ) x , x ( 0 ) = ( 1 , 0 ) , A 0 = ( 0 1 0 0 ) , A 1 = ( 0 0 − 1 0 ) \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} x ′ = ( A 0 + λ A 1 ) x , x ( 0 ) = ( 1 , 0 ) , A 0 = ( 0 0 1 0 ) , A 1 = ( 0 − 1 0 0 ) として扱う。d = 2 d=2 d = 2 、I = R I=\R I = R 、t 0 = 0 t_0=0 t 0 = 0 、b 0 = b 1 = 0 \boldsymbol b_0=\boldsymbol b_1=\boldsymbol 0 b 0 = b 1 = 0 として命題 5.1 を適用し、その解をx ( t , λ ) = ( x 1 ( t , λ ) , x 2 ( t , λ ) ) \boldsymbol x(t,\lambda)=(x_1(t,\lambda),x_2(t,\lambda)) x ( t , λ ) = ( x 1 ( t , λ ) , x 2 ( t , λ )) と書く。κ > 0 \kappa>0 κ > 0 とする。このとき、R \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) t ↦ ( cos κ t , − κ sin κ t ) , t ↦ ( cosh κ t , κ sinh κ t ) は、それぞれλ = κ 2 \lambda=\kappa^2 λ = κ 2 、λ = − κ 2 \lambda=-\kappa^2 λ = − κ 2 の系をすべてのt t t で満たし、時刻0 0 0 で( 1 , 0 ) (1,0) ( 1 , 0 ) をとる。定数写像t ↦ ( 1 , 0 ) t\mapsto(1,0) t ↦ ( 1 , 0 ) はλ = 0 \lambda=0 λ = 0 の系を満たす。§E13.2 系 2.7 (2) によりこれらは局所絶対連続であるから、同命題の解の一意性により、すべてのt ∈ R t\in\R t ∈ R で
x 1 ( t , λ ) = { cos ( λ t ) , λ > 0 , 1 , λ = 0 , cosh ( − λ t ) , λ < 0 x_1(t,\lambda)=\begin{cases}\cos(\sqrt\lambda\,t),&\lambda>0,\\1,&\lambda=0,\\\cosh(\sqrt{-\lambda}\,t),&\lambda<0\end{cases} x 1 ( t , λ ) = ⎩ ⎨ ⎧ cos ( λ t ) , 1 , cosh ( − λ t ) , λ > 0 , λ = 0 , λ < 0 である。実軸上のこの式は∣ λ ∣ \sqrt{\lvert\lambda\rvert} ∣ λ ∣ を通じてλ \lambda λ に依存するが、命題 5.1 (3) により、各t ∈ R t\in\R t ∈ R についてλ ↦ x 1 ( t , λ ) \lambda\mapsto x_1(t,\lambda) λ ↦ x 1 ( t , λ ) は整関数である。
命題 5.1 (2) により、t ↦ ∂ λ x ( t , 0 ) t\mapsto\partial_\lambda\boldsymbol x(t,0) t ↦ ∂ λ x ( t , 0 ) はz ′ = A 0 z + A 1 x ( t , 0 ) \boldsymbol z'=A_0\boldsymbol z+A_1\boldsymbol x(t,0) z ′ = A 0 z + A 1 x ( t , 0 ) 、z ( 0 ) = 0 \boldsymbol z(0)=\boldsymbol 0 z ( 0 ) = 0 のR \R R 上の局所絶対連続な解である。A 1 x ( t , 0 ) = ( 0 , − 1 ) A_1\boldsymbol x(t,0)=(0,-1) A 1 x ( t , 0 ) = ( 0 , − 1 ) であり、連続微分可能な写像z ( t ) = ( − t 2 / 2 , − t ) \boldsymbol z(t)=(-t^2/2,\,-t) z ( t ) = ( − t 2 /2 , − t ) はすべてのt t t でz ′ ( t ) = ( − t , − 1 ) = A 0 z ( t ) + ( 0 , − 1 ) \boldsymbol z'(t)=(-t,-1)=A_0\boldsymbol z(t)+(0,-1) z ′ ( t ) = ( − t , − 1 ) = A 0 z ( t ) + ( 0 , − 1 ) とz ( 0 ) = 0 \boldsymbol z(0)=\boldsymbol 0 z ( 0 ) = 0 を満たす。§E13.2 系 2.7 (2) と補題 3.1 (3) の一意性により、∂ λ x 1 ( t , 0 ) = − t 2 / 2 \partial_\lambda x_1(t,0)=-t^2/2 ∂ λ x 1 ( t , 0 ) = − t 2 /2 である。
6 演習
問題 6.1. Ω = R × R × R \Omega=\R\times\R\times\R Ω = R × R × R 上のf ( t , u , λ ) = λ + u 2 f(t,u,\lambda)=\lambda+u^2 f ( t , u , λ ) = λ + u 2 の解写像をφ \varphi φ とし、初期値問題y ′ = λ + y 2 y'=\lambda+y^2 y ′ = λ + y 2 、y ( 0 ) = 0 y(0)=0 y ( 0 ) = 0 を考える。
定理 3.8 を用いて、すべてのt ∈ R t\in\R t ∈ R について∂ λ φ ( t , 0 , 0 , 0 ) = t \partial_\lambda\varphi(t,0,0,0)=t ∂ λ φ ( t , 0 , 0 , 0 ) = t であることを示せ。
λ > 0 \lambda>0 λ > 0 ではφ ( t , 0 , 0 , λ ) = λ tan ( λ t ) \varphi(t,0,0,\lambda)=\sqrt\lambda\tan(\sqrt\lambda\,t) φ ( t , 0 , 0 , λ ) = λ tan ( λ t ) (∣ t ∣ < π / ( 2 λ ) \lvert t\rvert<\pi/(2\sqrt\lambda) ∣ t ∣ < π / ( 2 λ ) )、λ < 0 \lambda<0 λ < 0 ではφ ( t , 0 , 0 , λ ) = − − λ tanh ( − λ t ) \varphi(t,0,0,\lambda)=-\sqrt{-\lambda}\tanh(\sqrt{-\lambda}\,t) φ ( t , 0 , 0 , λ ) = − − λ tanh ( − λ t ) (t ∈ R t\in\R t ∈ R )であることを示し、t t t を固定してλ → 0 \lambda\to0 λ → 0 のときφ ( t , 0 , 0 , λ ) / λ → t \varphi(t,0,0,\lambda)/\lambda\to t φ ( t , 0 , 0 , λ ) / λ → t となることを直接確かめよ。
解答. (1) を示す。f f f はC ∞ C^\infty C ∞ 級であるから、補題 3.3 (4) により状態とパラメータについてC 1 C^1 C 1 級の Carathéodory 条件を満たす。λ = 0 \lambda=0 λ = 0 ではy = 0 y=0 y = 0 がR \R R 上の解であるから、I ( 0 , 0 , 0 ) = R I(0,0,0)=\R I ( 0 , 0 , 0 ) = R でありφ ( t , 0 , 0 , 0 ) = 0 \varphi(t,0,0,0)=0 φ ( t , 0 , 0 , 0 ) = 0 である。この解に沿う変分方程式の係数はA ( t ) = ∂ u f ( t , 0 , 0 ) = 0 A(t)=\partial_uf(t,0,0)=0 A ( t ) = ∂ u f ( t , 0 , 0 ) = 0 、B ( t ) = ∂ λ f ( t , 0 , 0 ) = 1 B(t)=\partial_\lambda f(t,0,0)=1 B ( t ) = ∂ λ f ( t , 0 , 0 ) = 1 である。Z ( t ) = t Z(t)=t Z ( t ) = t はZ ′ = A Z + B = 1 Z'=A Z+B=1 Z ′ = A Z + B = 1 、Z ( 0 ) = 0 Z(0)=0 Z ( 0 ) = 0 を満たし、補題 3.5 によりこの方程式の解は一意であるから、零解に沿う変分方程式の解でZ ( 0 ) = 0 Z(0)=0 Z ( 0 ) = 0 を満たすものはZ ( t ) = t Z(t)=t Z ( t ) = t である。定理 3.8 により∂ λ φ ( t , 0 , 0 , 0 ) = Z ( t ) = t \partial_\lambda\varphi(t,0,0,0)=Z(t)=t ∂ λ φ ( t , 0 , 0 , 0 ) = Z ( t ) = t である。
(2) を示す。λ = κ 2 \lambda=\kappa^2 λ = κ 2 (κ > 0 \kappa>0 κ > 0 )とし、y ( t ) = κ tan ( κ t ) y(t)=\kappa\tan(\kappa t) y ( t ) = κ tan ( κ t ) とおくと、∣ t ∣ < π / ( 2 κ ) \lvert t\rvert<\pi/(2\kappa) ∣ t ∣ < π / ( 2 κ ) でy ′ ( t ) = κ 2 ( 1 + tan 2 ( κ t ) ) = λ + y ( t ) 2 y'(t)=\kappa^2(1+\tan^2(\kappa t))=\lambda+y(t)^2 y ′ ( t ) = κ 2 ( 1 + tan 2 ( κ t )) = λ + y ( t ) 2 、y ( 0 ) = 0 y(0)=0 y ( 0 ) = 0 である。∣ t ∣ → π / ( 2 κ ) \lvert t\rvert\to\pi/(2\kappa) ∣ t ∣ → π / ( 2 κ ) のとき∣ y ( t ) ∣ → ∞ \lvert y(t)\rvert\to\infty ∣ y ( t )∣ → ∞ であるから、この解はそれ以上延長することができず、φ ( t , 0 , 0 , λ ) = y ( t ) \varphi(t,0,0,\lambda)=y(t) φ ( t , 0 , 0 , λ ) = y ( t ) である。λ = − κ 2 \lambda=-\kappa^2 λ = − κ 2 (κ > 0 \kappa>0 κ > 0 )とし、y ( t ) = − κ tanh ( κ t ) y(t)=-\kappa\tanh(\kappa t) y ( t ) = − κ tanh ( κ t ) とおくと、y ′ ( t ) = − κ 2 ( 1 − tanh 2 ( κ t ) ) = λ + y ( t ) 2 y'(t)=-\kappa^2(1-\tanh^2(\kappa t))=\lambda+y(t)^2 y ′ ( t ) = − κ 2 ( 1 − tanh 2 ( κ t )) = λ + y ( t ) 2 、y ( 0 ) = 0 y(0)=0 y ( 0 ) = 0 がR \R R 上で成り立つから、φ ( t , 0 , 0 , λ ) = y ( t ) \varphi(t,0,0,\lambda)=y(t) φ ( t , 0 , 0 , λ ) = y ( t ) である。
t ∈ R t\in\R t ∈ R を固定する。t = 0 t=0 t = 0 ではφ ( 0 , 0 , 0 , λ ) = 0 \varphi(0,0,0,\lambda)=0 φ ( 0 , 0 , 0 , λ ) = 0 である。t ≠ 0 t\ne0 t = 0 とする。0 < λ < π 2 / ( 4 t 2 ) 0<\lambda<\pi^2/(4t^2) 0 < λ < π 2 / ( 4 t 2 ) ならば∣ t ∣ < π / ( 2 λ ) \lvert t\rvert<\pi/(2\sqrt\lambda) ∣ t ∣ < π / ( 2 λ ) であり、x = λ t → 0 x=\sqrt\lambda\,t\to0 x = λ t → 0 とすると
φ ( t , 0 , 0 , λ ) λ = tan ( λ t ) λ = t ⋅ tan x x ⟶ t \frac{\varphi(t,0,0,\lambda)}{\lambda}=\frac{\tan(\sqrt\lambda\,t)}{\sqrt\lambda}=t\cdot\frac{\tan x}{x}\longrightarrow t λ φ ( t , 0 , 0 , λ ) = λ tan ( λ t ) = t ⋅ x tan x ⟶ t である。λ < 0 \lambda<0 λ < 0 では、x = − λ t → 0 x=\sqrt{-\lambda}\,t\to0 x = − λ t → 0 とすると
φ ( t , 0 , 0 , λ ) λ = − − λ tanh ( − λ t ) λ = tanh ( − λ t ) − λ = t ⋅ tanh x x ⟶ 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 , λ ) = λ − − λ tanh ( − λ t ) = − λ tanh ( − λ t ) = t ⋅ x tanh x ⟶ t である。φ ( t , 0 , 0 , 0 ) = 0 \varphi(t,0,0,0)=0 φ ( t , 0 , 0 , 0 ) = 0 であるから、λ = 0 \lambda=0 λ = 0 におけるλ ↦ φ ( t , 0 , 0 , λ ) \lambda\mapsto\varphi(t,0,0,\lambda) λ ↦ φ ( t , 0 , 0 , λ ) の微分係数はt t t であり、(1) と一致する。▨
問題 6.2. d , p ∈ N ≥ 1 d,p\in\NN d , p ∈ N ≥ 1 、r ≥ 2 r\ge2 r ≥ 2 とし、開集合Ω ⊆ R × R d × R p \Omega\subseteq\R\times\R^d\times\R^p Ω ⊆ R × R d × R p 上の写像f : Ω → R d f:\Omega\to\R^d f : Ω → R d が、零集合N N N と支配関数m K m_K m K をもって状態とパラメータについてC r C^r C r 級の Carathéodory 条件を満たすとする(t ∈ N t\in N t ∈ N では偏導関数を0 0 0 とする)。R d × d \R^{d\times d} R d × d 、R d × p \R^{d\times p} R d × p をR d 2 \R^{d^2} R d 2 、R d p \R^{dp} R d p と同一視してノルムを固定する。
Ω ^ = { ( t , u , U , V , λ ) : ( t , u , λ ) ∈ Ω } ⊆ R × R d × R d × d × R d × p × R p \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 Ω = {( t , u , U , V , λ ) : ( t , u , λ ) ∈ Ω } ⊆ R × R d × R d × d × R d × p × R p 上の写像
F ( t , u , U , V , λ ) = ( f ( t , u , λ ) , D u f ( t , u , λ ) U , D u f ( 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) F ( t , u , U , V , λ ) = ( f ( t , u , λ ) , D u f ( t , u , λ ) U , D u f ( t , u , λ ) V + D λ f ( t , u , λ ) ) が、状態( u , U , V ) (\boldsymbol u,U,V) ( u , U , V ) 、パラメータλ \lambda λ についてC r − 1 C^{r-1} C r − 1 級の Carathéodory 条件を満たすことを示せ。
解答. F F F の成分は、f f f の成分f i f_i f i と、U U U の成分U j k U_{jk} U j k 、V V V の成分V j k V_{jk} V j k を用いた
∑ j = 1 d ∂ u j f i U j k , ∑ j = 1 d ∂ u j f i V j k + ∂ λ k f i \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 j = 1 ∑ d ∂ u j f i U j k , j = 1 ∑ d ∂ u j f i V j k + ∂ λ k f i である。
( u , U , V , λ ) (\boldsymbol u,U,V,\lambda) ( u , U , V , λ ) を固定すると、{ 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\} { t : ( t , u , U , V , λ ) ∈ Ω } = { t : ( t , u , λ ) ∈ Ω } であり、定義 3.2 条件 (a) と補題 3.3 (2) によりf i f_i f i 、∂ u j f i \partial_{u_j}f_i ∂ u j f i 、∂ λ k f i \partial_{\lambda_k}f_i ∂ λ k f i はこの開集合上でt t t について可測である。定数U j k , V j k U_{jk},V_{jk} U j k , V j k との積の有限和も可測である。
g ( t , u , λ ) g(t,\boldsymbol u,\lambda) g ( t , u , λ ) をf f f のr r r 階以下の順序付き偏導関数の一つの成分とし、g g g にU j k U_{jk} U j k を掛けた関数g U j k gU_{jk} g U j k の( u , U , V , λ ) (\boldsymbol u,U,V,\lambda) ( u , U , V , λ ) の座標についての順序付き偏導関数を考える。各段階の関数は、g g g の( u , λ ) (\boldsymbol u,\lambda) ( u , λ ) の座標についての順序付き偏導関数h h h を用いてh U j k hU_{jk} h U j k 、h h h 、0 0 0 のいずれかの形であり、偏微分を一回行っても同じ形に留まる。実際、h U j k hU_{jk} h U j k を( u , λ ) (\boldsymbol u,\lambda) ( u , λ ) の座標で偏微分すると( ∂ h ) U j k (\partial h)U_{jk} ( ∂ h ) U j k 、U j k U_{jk} U j k で偏微分するとh h h 、他のU , V U,V U , V の座標で偏微分すると0 0 0 であり、h h h をU , V U,V U , V の座標で偏微分すると0 0 0 である。∂ u j f i \partial_{u_j}f_i ∂ u j f i と∂ λ k f i \partial_{\lambda_k}f_i ∂ λ k f i はt ∉ N t\notin N t ∈ / N でC r − 1 C^{r-1} C r − 1 級であるから、F F F の成分のr − 1 r-1 r − 1 階以下の順序付き偏導関数は、f f f のr r r 階以下の順序付き偏導関数の成分にU j k U_{jk} U j k 、V j k V_{jk} V j k 、1 1 1 、0 0 0 のいずれかを掛けたものの有限和であり、Ω ^ t \widehat\Omega_t Ω t 上で存在して連続である。§E4.4 定理 1.7 により、F ( t , ⋅ ) F(t,\cdot) F ( t , ⋅ ) はΩ ^ t \widehat\Omega_t Ω t 上でC r − 1 C^{r-1} C r − 1 級である。
コンパクト集合K ^ ⊆ Ω ^ \widehat K\subseteq\widehat\Omega K ⊆ Ω を取り、その( t , u , λ ) (t,\boldsymbol u,\lambda) ( t , u , λ ) 成分への射影をK K K とする。K K K は連続写像によるコンパクト集合の像であるからコンパクトであり、Ω \Omega Ω に含まれる。M = max K ^ max j , k { 1 , ∣ U j k ∣ , ∣ V j k ∣ } M=\max_{\widehat K}\max_{j,k}\{1,\lvert U_{jk}\rvert,\lvert V_{jk}\rvert\} M = max K max j , k { 1 , ∣ U j k ∣ , ∣ V j k ∣} とおく。ノルムの同値性により、ベクトルのノルムと成分の絶対値の和を互いに定数倍で抑える定数が存在する。前段の形により、( t , u , U , V , λ ) ∈ K ^ (t,\boldsymbol u,U,V,\lambda)\in\widehat K ( t , u , U , V , λ ) ∈ K 、t ∉ N t\notin N t ∈ / N で、F F F のr − 1 r-1 r − 1 階以下の各順序付き偏導関数のノルムは、f f f のr r r 階以下の順序付き偏導関数のノルムの高々( d + 1 ) (d+1) ( d + 1 ) 個の和にM M M と定数を掛けたもので抑えられる。定義 3.2 条件 (c) によりそれらはm K ( t ) m_K(t) m K ( t ) 以下であるから、d , p , r d,p,r d , p , r 、ノルムとM M M だけで決まる定数C C C を用いてm ^ K ^ = C m K ∈ L 1 ( R ) \widehat m_{\widehat K}=C\,m_K\in L^1(\R) m K = C m K ∈ L 1 ( R ) が支配関数となる。▨