§E13.25非線形境界値問題

最終更新

区間[a,b][a,b]上の二階の方程式y′′=f(x,y,y′)y''=f(x,y,y')に対して、両端の値y(a)=Ay(a)=A、y(b)=By(b)=Bを課す二点境界値問題は、弦の平衡や定常的な拡散など、空間的に拡がった系の釣り合いを記述する基本的な問題である。初期値問題では、右辺が局所 Lipschitz 連続であれば、初期時刻を含む小さな区間の上で解がただ一つ存在する。一方、境界値問題では、右辺が局所 Lipschitz 連続であっても、区間[a,b][a,b]上の解が存在するとは限らず、存在しても一つとは限らない。たとえば、右辺が−ω2y-\omega^2y(ω=π/(b−a)\omega=\pi/(b-a))である線形の場合でも、y(a)=y(b)=0y(a)=y(b)=0を満たす解はcsin⁡ω(x−a)c\sin\omega(x-a)(ccは任意の実数)であり、無限に多く存在する。このため、境界値問題では、右辺に課す条件と解の存在・一意性との関係を、初期値問題とは別に調べる必要がある。本記事では、非線形の二点境界値問題の解の存在と一意性について、基本的な結果や代表的な例を解説する。

1 解の定義と Dirichlet 問題の Green 関数

定義 1.1.a<ba<bとし、f:[a,b]×R2→Rf:[a,b]\times\R^2\to\Rを写像とする。R2\R^2のノルムを∥(u,v)∥=∣u∣+∣v∣\|(u,v)\|=|u|+|v|と定め、ρ>0\rho>0に対してBρ={(u,v)∈R2:∥(u,v)∥≤ρ}B_\rho=\{(u,v)\in\R^2:\|(u,v)\|\leq\rho\}とおく。ffがCarathéodory 型 (of Carathéodory type) であるとは、次の三条件を満たす Lebesgue 零集合N⊆[a,b]N\subseteq[a,b]が存在することをいう。

  1. すべての(u,v)∈R2(u,v)\in\R^2に対して、x↦f(x,u,v)x\mapsto f(x,u,v)は Lebesgue 可測である。
  2. すべてのx∈[a,b]∖Nx\in[a,b]\setminus Nに対して、(u,v)↦f(x,u,v)(u,v)\mapsto f(x,u,v)はR2\R^2上で連続である。
  3. すべてのρ>0\rho>0に対して、非負関数mρ∈L1(a,b)m_\rho\in L^1(a,b)が存在し、すべてのx∈[a,b]∖Nx\in[a,b]\setminus Nと(u,v)∈Bρ(u,v)\in B_\rhoに対して∣f(x,u,v)∣≤mρ(x)|f(x,u,v)|\leq m_\rho(x)が成り立つ。

f(x,u,v)f(x,u,v)がvvに依らない場合はf(x,u)f(x,u)と書く。f:[a,b]×R→Rf:[a,b]\times\R\to\Rに対して(x,u,v)↦f(x,u)(x,u,v)\mapsto f(x,u)が Carathéodory 型であるとき、ffは Carathéodory 型であるという。

A,B∈RA,B\in\Rとする。関数y∈C1([a,b];R)y\in C^1([a,b];\R)がy′∈AC([a,b])y'\in AC([a,b])を満たし、ほとんどすべてのx∈[a,b]x\in[a,b]で

y′′(x)=f(x,y(x),y′(x)),y''(x)=f(x,y(x),y'(x)),

かつy(a)=Ay(a)=A、y(b)=By(b)=Bを満たすとき、yyを二点境界値問題

y′′=f(x,y,y′),y(a)=A,y(b)=By''=f(x,y,y'),\qquad y(a)=A,\quad y(b)=B

の境界値問題の解 (solution of the boundary value problem) という。

補題 1.2.a<ba<b、g∈L1(a,b)g\in L^1(a,b)とし、∥g∥1=∫ab∣g(s)∣ ds\|g\|_1=\int_a^b|g(s)|\,ds、連続関数yyに対して∥y∥∞=max⁡[a,b]∣y∣\|y\|_\infty=\max_{[a,b]}|y|とおく。[a,b]2[a,b]^2上の関数GGを

G(x,s)={−(x−a)(b−s)b−a,a≤x≤s≤b,−(s−a)(b−x)b−a,a≤s≤x≤bG(x,s)= \begin{cases} -\dfrac{(x-a)(b-s)}{b-a},&a\leq x\leq s\leq b,\\[3mm] -\dfrac{(s-a)(b-x)}{b-a},&a\leq s\leq x\leq b \end{cases}

と定める。このとき、GGは連続であり、∣G∣≤b−a4|G|\leq\frac{b-a}4を満たす。また、

y′′=gほとんど至る所,y(a)=y(b)=0y''=g\quad\text{ほとんど至る所},\qquad y(a)=y(b)=0

を満たすy∈C1([a,b];R)y\in C^1([a,b];\R)、y′∈AC([a,b])y'\in AC([a,b])はただ一つ存在し、y(x)=∫abG(x,s)g(s) dsy(x)=\int_a^bG(x,s)g(s)\,dsと表され、

∥y∥∞≤b−a4∥g∥1,∥y′∥∞≤∥g∥1\|y\|_\infty\leq\frac{b-a}4\|g\|_1,\qquad\|y'\|_\infty\leq\|g\|_1

を満たす。

証明.P=Q=0P=Q=0、(αa,βa)=(αb,βb)=(1,0)(\alpha_a,\beta_a)=(\alpha_b,\beta_b)=(1,0)の分離型二点境界値問題を考える。§E13.24 補題 3.2のu,vu,vはu(x)=a−xu(x)=a-x、v(x)=b−xv(x)=b-xであり、W=uv′−u′v=b−a≠0W=uv'-u'v=b-a\neq0である。同補題の三条件の同値により、この問題の境界核は自明である。§E13.24 定義 3.3の Green 関数は上のGGに一致し、§E13.24 定理 3.4により、GGは連続であり、x≠sx\neq sに対して

∂xG(x,s)={−b−sb−a,x<s,s−ab−a,s<x\partial_xG(x,s)= \begin{cases} -\dfrac{b-s}{b-a},&x<s,\\[3mm] \dfrac{s-a}{b-a},&s<x \end{cases}

である。同定理により、y′′=gy''=g、y(a)=y(b)=0y(a)=y(b)=0を満たすy∈C1([a,b];C)y\in C^1([a,b];\C)、y′∈AC([a,b];C)y'\in AC([a,b];\C)はただ一つ存在し、すべてのx∈[a,b]x\in[a,b]で

y(x)=∫abG(x,s)g(s) ds,y′(x)=∫ab∂xG(x,s)g(s) dsy(x)=\int_a^bG(x,s)g(s)\,ds,\qquad y'(x)=\int_a^b\partial_xG(x,s)g(s)\,ds

が成り立つ。ggは実数値なのでyyは実数値である。

x,s∈[a,b]x,s\in[a,b]に対して、x≤sx\leq sならば(x−a)(b−s)≤(s−a)(b−s)(x-a)(b-s)\leq(s-a)(b-s)、s≤xs\leq xならば(s−a)(b−x)≤(s−a)(b−s)(s-a)(b-x)\leq(s-a)(b-s)であるから、(s−a)(b−s)≤(b−a2)2(s-a)(b-s)\leq\bigl(\frac{b-a}2\bigr)^2と合わせて

∣G(x,s)∣≤(s−a)(b−s)b−a≤b−a4|G(x,s)|\leq\frac{(s-a)(b-s)}{b-a}\leq\frac{b-a}4

である。従って∥y∥∞≤b−a4∥g∥1\|y\|_\infty\leq\frac{b-a}4\|g\|_1である。x≠sx\neq sでは0≤b−s≤b−a0\leq b-s\leq b-a、0≤s−a≤b−a0\leq s-a\leq b-aにより∣∂xG(x,s)∣≤1|\partial_xG(x,s)|\leq1であり、∥y′∥∞≤∥g∥1\|y'\|_\infty\leq\|g\|_1である。▨

2 射撃法による存在

定理 2.1.a<ba<bとし、f:[a,b]×R2→Rf:[a,b]\times\R^2\to\Rは Carathéodory 型であるとする。あるm∈L1(a,b)m\in L^1(a,b)と Lebesgue 零集合N′⊆[a,b]N'\subseteq[a,b]が存在し、すべてのx∈[a,b]∖N′x\in[a,b]\setminus N'とすべての(u,v)∈R2(u,v)\in\R^2に対して∣f(x,u,v)∣≤m(x)|f(x,u,v)|\leq m(x)が成り立つとする。また、各ρ>0\rho>0に対して非負関数ℓρ∈L1(a,b)\ell_\rho\in L^1(a,b)が存在し、すべてのx∈[a,b]∖N′x\in[a,b]\setminus N'と∥(u,v)∥≤ρ\|(u,v)\|\leq\rho、∥(uˉ,vˉ)∥≤ρ\|(\bar u,\bar v)\|\leq\rhoを満たす(u,v),(uˉ,vˉ)(u,v),(\bar u,\bar v)に対して

∣f(x,u,v)−f(x,uˉ,vˉ)∣≤ℓρ(x)(∣u−uˉ∣+∣v−vˉ∣)|f(x,u,v)-f(x,\bar u,\bar v)|\leq\ell_\rho(x)\bigl(|u-\bar u|+|v-\bar v|\bigr)

が成り立つとする。A∈RA\in\Rとする。

このとき、任意のs∈Rs\in\Rに対して、y(a)=Ay(a)=A、y′(a)=sy'(a)=sを満たすy′′=f(x,y,y′)y''=f(x,y,y')の解y∈C1([a,b];R)y\in C^1([a,b];\R)、y′∈AC([a,b])y'\in AC([a,b])はただ一つ存在する。これをy(⋅ ;s)y(\cdot\,;s)と書くと、写像s↦y(b;s)s\mapsto y(b;s)は連続であり、

∣y(b;s)−A−s(b−a)∣≤(b−a)∥m∥1\bigl|y(b;s)-A-s(b-a)\bigr|\leq(b-a)\|m\|_1

を満たす。従って、任意のB∈RB\in\Rに対してy′′=f(x,y,y′)y''=f(x,y,y')、y(a)=Ay(a)=A、y(b)=By(b)=Bの解が存在する。

証明.F:[a,b]×R2→R2F:[a,b]\times\R^2\to\R^2をF(x,u,v)=(v,f(x,u,v))F(x,u,v)=(v,f(x,u,v))で定め、R2\R^2には定義 1.1のノルム∥(u,v)∥=∣u∣+∣v∣\|(u,v)\|=|u|+|v|を用いる。NNをffの Carathéodory 型の零集合とし、N0=N∪N′N_0=N\cup N'とおく。各(u,v)∈R2(u,v)\in\R^2に対してx↦F(x,u,v)x\mapsto F(x,u,v)は可測であり、x∈[a,b]∖N0x\in[a,b]\setminus N_0に対して(u,v)↦F(x,u,v)(u,v)\mapsto F(x,u,v)は連続である。ρ>0\rho>0、x∈[a,b]∖N0x\in[a,b]\setminus N_0、(u,v),(uˉ,vˉ)∈Bρ(u,v),(\bar u,\bar v)\in B_\rhoに対して∥F(x,u,v)∥≤ρ+∣m(x)∣\|F(x,u,v)\|\leq\rho+|m(x)|であり、

∥F(x,u,v)−F(x,uˉ,vˉ)∥≤∣v−vˉ∣+ℓρ(x)(∣u−uˉ∣+∣v−vˉ∣)≤(1+ℓρ(x))∥(u,v)−(uˉ,vˉ)∥\|F(x,u,v)-F(x,\bar u,\bar v)\|\leq|v-\bar v|+\ell_\rho(x)\bigl(|u-\bar u|+|v-\bar v|\bigr)\leq\bigl(1+\ell_\rho(x)\bigr)\|(u,v)-(\bar u,\bar v)\|

である。R2\R^2の任意の閉球はあるBρB_\rhoに含まれるので、[a,b][a,b]に含まれる長さが正の各閉区間JJとR2\R^2の各閉球BBについて、FFのJ×BJ\times Bへの制限は§E13.4 定義 2.1の Carathéodory 条件を満たし、時間について可積分な Lipschitz 係数をもつ。また、x∈[a,b]∖N0x\in[a,b]\setminus N_0とすべての(u,v)∈R2(u,v)\in\R^2に対して∥F(x,u,v)∥≤∣m(x)∣+∥(u,v)∥\|F(x,u,v)\|\leq|m(x)|+\|(u,v)\|である。

s∈Rs\in\Rとし、y∈C1([a,b];R)y\in C^1([a,b];\R)、y′∈AC([a,b])y'\in AC([a,b])がy(a)=Ay(a)=A、y′(a)=sy'(a)=sを満たすy′′=f(x,y,y′)y''=f(x,y,y')の解であることと、Y=(y,y′)Y=(y,y')が[a,b][a,b]上の絶対連続関数としてY′=F(x,Y)Y'=F(x,Y)をほとんど至る所で満たしY(a)=(A,s)Y(a)=(A,s)を満たすことは同値である。実際、Y=(u,v)Y=(u,v)が後者を満たすならば、u′=vu'=vがほとんど至る所で成り立ち、§E13.2 系 2.6によりu(x)=A+∫axv(t) dtu(x)=A+\int_a^xv(t)\,dtである。vvは連続なのでu∈C1u\in C^1、u′=vu'=vが全点で成り立ち、v=u′∈ACv=u'\in AC、u′′=v′=f(x,u,u′)u''=v'=f(x,u,u')である。逆は、y∈C1([a,b];R)y\in C^1([a,b];\R)が Lipschitz 連続であり絶対連続であることから従う。

§E13.7 定理 4.1をd=2d=2、I=[a,b]I=[a,b]、t∗=at_*=a、初期値(A,s)(A,s)、p=∣m∣p=|m|、q=1q=1に適用すると、Y(a)=(A,s)Y(a)=(A,s)を満たし[a,b][a,b]のほとんど至る所でY′=F(x,Y)Y'=F(x,Y)を満たすY∈AC([a,b];R2)Y\in AC([a,b];\R^2)がただ一つ存在する。前段の同値により、y(a)=Ay(a)=A、y′(a)=sy'(a)=sを満たすy′′=f(x,y,y′)y''=f(x,y,y')の解y(⋅ ;s)y(\cdot\,;s)がただ一つ存在する。

§E13.2 系 2.6により、a≤x≤ba\leq x\leq bに対して

y′(x;s)=s+∫axf(t,y(t;s),y′(t;s)) dt,y(x;s)=A+∫axy′(t;s) dty'(x;s)=s+\int_a^xf\bigl(t,y(t;s),y'(t;s)\bigr)\,dt,\qquad y(x;s)=A+\int_a^xy'(t;s)\,dt

であるから、∣y′(x;s)−s∣≤∥m∥1|y'(x;s)-s|\leq\|m\|_1、∣y(x;s)∣≤∣A∣+(∣s∣+∥m∥1)(b−a)|y(x;s)|\leq|A|+(|s|+\|m\|_1)(b-a)である。特に

∣y(b;s)−A−s(b−a)∣=∣∫ab(y′(x;s)−s) dx∣≤(b−a)∥m∥1\bigl|y(b;s)-A-s(b-a)\bigr|=\left|\int_a^b\bigl(y'(x;s)-s\bigr)\,dx\right|\leq(b-a)\|m\|_1

である。

s0∈Rs_0\in\Rとし、ρ0=∣A∣+(∣s0∣+1+∥m∥1)(b−a+1)\rho_0=|A|+(|s_0|+1+\|m\|_1)(b-a+1)とおく。∣s−s0∣≤1|s-s_0|\leq1ならば、前段の評価によりYs=(y(⋅ ;s),y′(⋅ ;s))Y_s=(y(\cdot\,;s),y'(\cdot\,;s))はすべてのx∈[a,b]x\in[a,b]で∥Ys(x)∥≤ρ0\|Y_s(x)\|\leq\rho_0を満たす。∣s−s0∣≤1|s-s_0|\leq1とし、§E13.7 定理 2.1をd=2d=2、S=[a,b]×Bρ0S=[a,b]\times B_{\rho_0}、二つの右辺F,FF,F、二解Ys0,YsY_{s_0},Y_s、t0=at_0=a、ℓ=1+ℓρ0\ell=1+\ell_{\rho_0}、r=0r=0に適用すると、Λ=exp⁡(∫ab(1+ℓρ0(x)) dx)\Lambda=\exp\bigl(\int_a^b(1+\ell_{\rho_0}(x))\,dx\bigr)に対して

max⁡x∈[a,b]∥Ys(x)−Ys0(x)∥≤Λ∣s−s0∣\max_{x\in[a,b]}\|Y_s(x)-Y_{s_0}(x)\|\leq\Lambda|s-s_0|

である。従って∣y(b;s)−y(b;s0)∣≤Λ∣s−s0∣|y(b;s)-y(b;s_0)|\leq\Lambda|s-s_0|であり、s↦y(b;s)s\mapsto y(b;s)はs0s_0で連続である。

B∈RB\in\Rとし、C=(b−a)∥m∥1C=(b-a)\|m\|_1、s1=B−A+Cb−as_1=\frac{B-A+C}{b-a}、s2=B−A−Cb−as_2=\frac{B-A-C}{b-a}とおく。上の評価からy(b;s1)≥A+s1(b−a)−C=By(b;s_1)\geq A+s_1(b-a)-C=B、y(b;s2)≤A+s2(b−a)+C=By(b;s_2)\leq A+s_2(b-a)+C=Bである。C=0C=0のときはs1=s2s_1=s_2であり、y(b;s1)=By(b;s_1)=Bである。C>0C>0のとき、s↦y(b;s)s\mapsto y(b;s)は[s2,s1][s_2,s_1]上で連続なので、§D1.12 系 1.2によりy(b;s)=By(b;s)=Bとなるs∈[s2,s1]s\in[s_2,s_1]が存在する。y(⋅ ;s)y(\cdot\,;s)は求める解である。▨

注意 2.2. 初期勾配ssをパラメータとする初期値問題の解y(⋅ ;s)y(\cdot\,;s)の右端の値y(b;s)y(b;s)のssに関する連続性と増大度を評価し、中間値の定理によってy(b;s)=By(b;s)=Bを満たすssを選ぶ論法を、射撃法という。

3 縮小写像による一意存在

定理 3.1.a<ba<bとし、f:[a,b]×R2→Rf:[a,b]\times\R^2\to\Rは Carathéodory 型であるとする。非負関数ℓ0,ℓ1∈L1(a,b)\ell_0,\ell_1\in L^1(a,b)と Lebesgue 零集合N′⊆[a,b]N'\subseteq[a,b]が存在し、すべてのx∈[a,b]∖N′x\in[a,b]\setminus N'とすべての(u,v),(uˉ,vˉ)∈R2(u,v),(\bar u,\bar v)\in\R^2に対して

∣f(x,u,v)−f(x,uˉ,vˉ)∣≤ℓ0(x)∣u−uˉ∣+ℓ1(x)∣v−vˉ∣|f(x,u,v)-f(x,\bar u,\bar v)|\leq\ell_0(x)|u-\bar u|+\ell_1(x)|v-\bar v|

が成り立つとする。さらに

κ=b−a4∥ℓ0∥1+∥ℓ1∥1<1\kappa=\frac{b-a}4\|\ell_0\|_1+\|\ell_1\|_1<1

とする。このとき、任意のA,B∈RA,B\in\Rに対して境界値問題y′′=f(x,y,y′)y''=f(x,y,y')、y(a)=Ay(a)=A、y(b)=By(b)=Bの解はただ一つ存在する。

証明.X=C1([a,b];R)X=C^1([a,b];\R)に、ノルム

∥y∥∗=max⁡{4b−a∥y∥∞, ∥y′∥∞}\|y\|_*=\max\Bigl\{\frac4{b-a}\|y\|_\infty,\ \|y'\|_\infty\Bigr\}

を入れる。(yn)(y_n)が∥⋅∥∗\|\cdot\|_*について Cauchy 列ならば、yny_nとyn′y_n'は[a,b][a,b]上で一様収束し、極限をyy、vvとすると、yn(x)=yn(a)+∫axyn′(t) dty_n(x)=y_n(a)+\int_a^xy_n'(t)\,dtの極限からy(x)=y(a)+∫axv(t) dty(x)=y(a)+\int_a^xv(t)\,dtを得る。従ってy∈C1y\in C^1、y′=vy'=vであり、yn→yy_n\to yである。以上より、XXは完備である。

h(x)=A+B−Ab−a(x−a)h(x)=A+\frac{B-A}{b-a}(x-a)とおく。y∈Xy\in Xに対してgy(s)=f(s,y(s),y′(s))g_y(s)=f(s,y(s),y'(s))とおく。ρ=∥y∥∞+∥y′∥∞+1\rho=\|y\|_\infty+\|y'\|_\infty+1とし、§E13.4 補題 2.3をd=2d=2、J=[a,b]J=[a,b]、閉球BρB_\rhoと、可測写像s↦(y(s),y′(s))s\mapsto(y(s),y'(s))に適用すると、gyg_yは可測であり、∣gy∣≤mρ|g_y|\leq m_\rhoなのでgy∈L1(a,b)g_y\in L^1(a,b)である。

(Ty)(x)=h(x)+∫abG(x,s)gy(s) ds(Ty)(x)=h(x)+\int_a^bG(x,s)g_y(s)\,ds

と定める。補題 1.2により、Ty∈XTy\in X、(Ty)′∈AC(Ty)'\in AC、(Ty)′′=gy(Ty)''=g_yがほとんど至る所で成り立ち、(Ty)(a)=A(Ty)(a)=A、(Ty)(b)=B(Ty)(b)=Bである。

yyが境界値問題の解であることと、Ty=yTy=yであることは同値である。実際、yyが解ならばgy∈L1g_y\in L^1であり、y−hy-hは(y−h)′′=gy(y-h)''=g_y、(y−h)(a)=(y−h)(b)=0(y-h)(a)=(y-h)(b)=0を満たすので、補題 1.2の一意性によりy−h=∫Ggyy-h=\int G g_y、すなわちTy=yTy=yである。逆は、Ty=yTy=yならばy′′=gyy''=g_yがほとんど至る所で成り立ち、y(a)=Ay(a)=A、y(b)=By(b)=Bであることから従う。

y,z∈Xy,z\in Xとし、w=y−zw=y-zとおく。s∉N′s\notin N'について∣gy(s)−gz(s)∣≤ℓ0(s)∣w(s)∣+ℓ1(s)∣w′(s)∣|g_y(s)-g_z(s)|\leq\ell_0(s)|w(s)|+\ell_1(s)|w'(s)|である。∥w∥∞≤b−a4∥w∥∗\|w\|_\infty\leq\frac{b-a}4\|w\|_*、∥w′∥∞≤∥w∥∗\|w'\|_\infty\leq\|w\|_*により

∥gy−gz∥1≤∥ℓ0∥1b−a4∥w∥∗+∥ℓ1∥1∥w∥∗=κ∥w∥∗\|g_y-g_z\|_1\leq\|\ell_0\|_1\frac{b-a}4\|w\|_*+\|\ell_1\|_1\|w\|_*=\kappa\|w\|_*

である。補題 1.2をg=gy−gzg=g_y-g_zに適用すると、Ty−TzTy-Tzはgy−gzg_y-g_zに対応する Dirichlet 問題の解であるから、∥Ty−Tz∥∞≤b−a4∥gy−gz∥1\|Ty-Tz\|_\infty\leq\frac{b-a}4\|g_y-g_z\|_1、∥(Ty−Tz)′∥∞≤∥gy−gz∥1\|(Ty-Tz)'\|_\infty\leq\|g_y-g_z\|_1であり、∥Ty−Tz∥∗≤κ∥y−z∥∗\|Ty-Tz\|_*\leq\kappa\|y-z\|_*を得る。κ<1\kappa<1なので、§E2.7 定理 2.2によりTTはただ一つの不動点をもつ。▨

4 片側条件による一意性

補題 4.1.a<ba<b、L=b−aL=b-aとし、z∈C1([a,b];R)z\in C^1([a,b];\R)、z′∈AC([a,b])z'\in AC([a,b])、z(a)=z(b)=0z(a)=z(b)=0とする。このとき

∫abz(x)2 dx≤(Lπ)2∫abz′(x)2 dx\int_a^bz(x)^2\,dx\leq\Bigl(\frac L\pi\Bigr)^2\int_a^bz'(x)^2\,dx

が成り立つ。

証明.k=π/Lk=\pi/L、φ(x)=sin⁡k(x−a)\varphi(x)=\sin k(x-a)とおく。k2k^2は、§E13.24 例 2.2の第一固有値λ1\lambda_1(長さLLの区間、n=1n=1)であり、φ\varphiはその固有関数である。x∈(a,b)x\in(a,b)ならば0<k(x−a)<π0<k(x-a)<\piであるからφ(x)>0\varphi(x)>0であり、g=φ′/φg=\varphi'/\varphiは(a,b)(a,b)上でC1C^1級である。φ′′=−k2φ\varphi''=-k^2\varphiにより

g′=φ′′φ−(φ′φ)2=−k2−g2g'=\frac{\varphi''}\varphi-\Bigl(\frac{\varphi'}\varphi\Bigr)^2=-k^2-g^2

であるから、(a,b)(a,b)上で

(z′−gz)2+(gz2)′=z′2−2gzz′+g2z2+g′z2+2gzz′=z′2−k2z2(z'-gz)^2+(gz^2)'=z'^2-2gzz'+g^2z^2+g'z^2+2gzz'=z'^2-k^2z^2

が成り立つ。

M=max⁡[a,b]∣z′∣M=\max_{[a,b]}|z'|とおく。z(a)=z(b)=0z(a)=z(b)=0と平均値の定理により、x∈[a,b]x\in[a,b]に対して∣z(x)∣≤M(x−a)|z(x)|\leq M(x-a)、∣z(x)∣≤M(b−x)|z(x)|\leq M(b-x)である。φ(a)=0\varphi(a)=0、φ′(a)=k>0\varphi'(a)=k>0であるから、あるδ∈(0,L/2)\delta\in(0,L/2)が存在し、a<x≤a+δa<x\leq a+\deltaならばφ(x)≥k2(x−a)\varphi(x)\geq\frac k2(x-a)である。φ(x)=sin⁡k(b−x)\varphi(x)=\sin k(b-x)でもあるから、同じδ\deltaについてb−δ≤x<bb-\delta\leq x<bならばφ(x)≥k2(b−x)\varphi(x)\geq\frac k2(b-x)である。∣φ′∣≤k|\varphi'|\leq kと合わせると、a<x≤a+δa<x\leq a+\deltaならば

∣g(x)z(x)2∣=∣φ′(x)∣z(x)2φ(x)≤kM2(x−a)2k2(x−a)=2M2(x−a)|g(x)z(x)^2|=\frac{|\varphi'(x)|z(x)^2}{\varphi(x)}\leq\frac{kM^2(x-a)^2}{\frac k2(x-a)}=2M^2(x-a)

であり、b−δ≤x<bb-\delta\leq x<bならば∣g(x)z(x)2∣≤2M2(b−x)|g(x)z(x)^2|\leq2M^2(b-x)である。

0<ε≤δ0<\varepsilon\leq\deltaとする。gz2gz^2は[a+ε,b−ε][a+\varepsilon,b-\varepsilon]上でC1C^1級であり、従って Lipschitz 連続で絶対連続であるから、§E13.2 系 2.6により

∫a+εb−ε(z′2−k2z2) dx=∫a+εb−ε(z′−gz)2 dx+(gz2)(b−ε)−(gz2)(a+ε)≥−4M2ε\int_{a+\varepsilon}^{b-\varepsilon}\bigl(z'^2-k^2z^2\bigr)\,dx=\int_{a+\varepsilon}^{b-\varepsilon}(z'-gz)^2\,dx+(gz^2)(b-\varepsilon)-(gz^2)(a+\varepsilon)\geq-4M^2\varepsilon

である。z′2−k2z2z'^2-k^2z^2は[a,b][a,b]上で連続であるから、ε→0\varepsilon\to0とすると∫ab(z′2−k2z2) dx≥0\int_a^b\bigl(z'^2-k^2z^2\bigr)\,dx\geq0を得る。k−2=(L/π)2k^{-2}=(L/\pi)^2であるから、これは主張の不等式である。▨

定理 4.2.a<ba<bとし、f:[a,b]×R→Rf:[a,b]\times\R\to\Rは Carathéodory 型であるとする。あるK<(πb−a)2K<\bigl(\frac\pi{b-a}\bigr)^2と Lebesgue 零集合N′⊆[a,b]N'\subseteq[a,b]が存在し、すべてのx∈[a,b]∖N′x\in[a,b]\setminus N'とすべてのu,v∈Ru,v\in\Rに対して

(f(x,u)−f(x,v))(u−v)≥−K(u−v)2\bigl(f(x,u)-f(x,v)\bigr)(u-v)\geq-K(u-v)^2

が成り立つとする。このとき、任意のA,B∈RA,B\in\Rに対して境界値問題y′′=f(x,y)y''=f(x,y)、y(a)=Ay(a)=A、y(b)=By(b)=Bの解は高々一つである。

証明.y,yˉy,\bar yを二つの解とし、z=y−yˉz=y-\bar yとおく。z∈C1z\in C^1、z′∈ACz'\in AC、z(a)=z(b)=0z(a)=z(b)=0であり、ほとんど至る所でz′′=f(x,y)−f(x,yˉ)z''=f(x,y)-f(x,\bar y)である。§E13.2 定理 3.4をz′z'とzzに適用すると、z(a)=z(b)=0z(a)=z(b)=0により

∫abz′′z dx=−∫abz′2 dx\int_a^bz''z\,dx=-\int_a^bz'^2\,dx

である。仮定により、ほとんどすべてのxxでz′′(x)z(x)≥−Kz(x)2z''(x)z(x)\geq-Kz(x)^2であるから、K+=max⁡{K,0}K_+=\max\{K,0\}とおくと

∫abz′2 dx≤K∫abz2 dx≤K+∫abz2 dx\int_a^bz'^2\,dx\leq K\int_a^bz^2\,dx\leq K_+\int_a^bz^2\,dx

である。補題 4.1により右辺はK+(b−aπ)2∫abz′2 dxK_+\bigl(\frac{b-a}\pi\bigr)^2\int_a^bz'^2\,dx以下であり、K+(b−aπ)2<1K_+\bigl(\frac{b-a}\pi\bigr)^2<1なので∫abz′2 dx=0\int_a^bz'^2\,dx=0である。z′z'は連続なのでz′=0z'=0であり、z(a)=0z(a)=0によりz=0z=0である。▨

例 4.3.a<ba<b、ω=πb−a\omega=\frac\pi{b-a}とし、f(x,u)=−ω2uf(x,u)=-\omega^2uとする。このとき、任意のc∈Rc\in\Rに対してyc(x)=csin⁡ω(x−a)y_c(x)=c\sin\omega(x-a)はy′′=−ω2yy''=-\omega^2y、y(a)=y(b)=0y(a)=y(b)=0の解である。従って境界値問題は無限に多くの解をもつ。また、y(a)=0y(a)=0、y(b)=By(b)=B(B≠0B\neq0)を課す問題は、一般解がc1cos⁡ω(x−a)+c2sin⁡ω(x−a)c_1\cos\omega(x-a)+c_2\sin\omega(x-a)でありy(a)=0y(a)=0からc1=0c_1=0、y(b)=c2sin⁡π=0y(b)=c_2\sin\pi=0となるので、解をもたない。

前半の解は、§E13.24 例 2.2をL=b−aL=b-a、n=1n=1として平行移動x↦x−ax\mapsto x-aの下で適用した結果と一致する。すなわち、ω2=(πb−a)2\omega^2=\bigl(\frac\pi{b-a}\bigr)^2は区間[a,b][a,b]上の Dirichlet 問題y′′+λy=0y''+\lambda y=0、y(a)=y(b)=0y(a)=y(b)=0の第一固有値であり、固有関数はsin⁡ω(x−a)\sin\omega(x-a)の零でない定数倍である。

このffは(f(x,u)−f(x,v))(u−v)=−ω2(u−v)2(f(x,u)-f(x,v))(u-v)=-\omega^2(u-v)^2を満たすので、定理 4.2の不等式を満たすKKはK≥ω2K\geq\omega^2に限られ、同定理の仮定は成り立たない。また、定理 3.1の Lipschitz 評価を満たすℓ0\ell_0は、ほとんど至る所でℓ0≥ω2\ell_0\geq\omega^2を満たす。したがってκ≥b−a4ω2(b−a)=π24>1\kappa\geq\frac{b-a}4\omega^2(b-a)=\frac{\pi^2}4>1であり、同定理の仮定も成り立たない。定理 4.2の閾値ω2\omega^2は、KKの許容範囲としてこれ以上広げることができない。

5 最大値原理と下解・上解

補題 5.1.α<β\alpha<\betaとし、w∈C1([α,β];R)w\in C^1([\alpha,\beta];\R)、w′∈AC([α,β])w'\in AC([\alpha,\beta])、w(α)=w(β)=0w(\alpha)=w(\beta)=0とする。このとき、ほとんど至る所でw′′≤0w''\leq0ならば、[α,β][\alpha,\beta]上でw≥0w\geq0である。

証明.§E13.2 系 2.6により、α≤s<t≤β\alpha\leq s<t\leq\betaに対してw′(t)−w′(s)=∫stw′′(r) dr≤0w'(t)-w'(s)=\int_s^tw''(r)\,dr\leq0である。従ってw′w'は単調非増加である。w(x)<0w(x)<0となるx∈(α,β)x\in(\alpha,\beta)が存在すると仮定する。平均値の定理により、ξ1∈(α,x)\xi_1\in(\alpha,x)とξ2∈(x,β)\xi_2\in(x,\beta)が存在して

w′(ξ1)=w(x)−w(α)x−α<0,w′(ξ2)=w(β)−w(x)β−x>0w'(\xi_1)=\frac{w(x)-w(\alpha)}{x-\alpha}<0,\qquad w'(\xi_2)=\frac{w(\beta)-w(x)}{\beta-x}>0

となる。ξ1<ξ2\xi_1<\xi_2なのでw′(ξ1)≥w′(ξ2)w'(\xi_1)\geq w'(\xi_2)でなければならないが、これはw′(ξ1)<0<w′(ξ2)w'(\xi_1)<0<w'(\xi_2)と両立しない。▨

補題 5.2.a<ba<bとし、w∈C1([a,b];R)w\in C^1([a,b];\R)、w′∈AC([a,b])w'\in AC([a,b])、w(a)≥0w(a)\geq0、w(b)≥0w(b)\geq0とする。ほとんどすべてのx∈[a,b]x\in[a,b]に対して、w(x)<0w(x)<0ならばw′′(x)≤0w''(x)\leq0であるとする。このとき、[a,b][a,b]上でw≥0w\geq0である。

証明.U={x∈[a,b]:w(x)<0}U=\{x\in[a,b]:w(x)<0\}が空でないと仮定する。w(a),w(b)≥0w(a),w(b)\geq0なのでU⊆(a,b)U\subseteq(a,b)であり、UUは開集合である。UUの連結成分(α,β)(\alpha,\beta)を取る。α∉U\alpha\notin Uなのでw(α)≥0w(\alpha)\geq0であり、wwは連続で(α,β)(\alpha,\beta)上でw<0w<0なのでw(α)≤0w(\alpha)\leq0である。従ってw(α)=0w(\alpha)=0であり、同様にw(β)=0w(\beta)=0である。ほとんどすべてのx∈(α,β)x\in(\alpha,\beta)でw′′(x)≤0w''(x)\leq0なので、補題 5.1により[α,β][\alpha,\beta]上でw≥0w\geq0である。これは(α,β)(\alpha,\beta)上でw<0w<0であることと両立しない。▨

補題 5.3.a<ba<bとし、c∈L1(a,b)c\in L^1(a,b)はc≥0c\geq0をほとんど至る所で満たすとする。このとき、任意のg∈L1(a,b)g\in L^1(a,b)とA,B∈RA,B\in\Rに対して、境界値問題

−y′′+cy=gほとんど至る所,y(a)=A,y(b)=B-y''+cy=g\quad\text{ほとんど至る所},\qquad y(a)=A,\quad y(b)=B

の解y∈C1([a,b];R)y\in C^1([a,b];\R)、y′∈AC([a,b])y'\in AC([a,b])はただ一つ存在する。この解をSA,BgS_{A,B}gと書く。

a,b,ca,b,cだけで定まる連続関数Γ:[a,b]2→R\Gamma:[a,b]^2\to\Rが存在し、任意のg,A,Bg,A,Bに対して

SA,Bg(x)=φA,B(x)+∫abΓ(x,s)g(s) dsS_{A,B}g(x)=\varphi_{A,B}(x)+\int_a^b\Gamma(x,s)g(s)\,ds

が成り立つ。ここでφA,B=SA,B0\varphi_{A,B}=S_{A,B}0である。

さらに、g≥0g\geq0がほとんど至る所で成り立ちA,B≥0A,B\geq0ならば、SA,Bg≥0S_{A,B}g\geq0である。従ってg1≤g2g_1\leq g_2がほとんど至る所で成り立てばS0,0g1≤S0,0g2S_{0,0}g_1\leq S_{0,0}g_2である。

証明.w∈C1([a,b];C)w\in C^1([a,b];\C)、w′∈ACw'\in ACがw′′=cww''=cwをほとんど至る所で満たし、w(a)=w(b)=0w(a)=w(b)=0であるとする。§E13.2 定理 3.4の複素形をw′w'とwwに適用すると、∫abw′′w‾ dx=−∫ab∣w′∣2 dx\int_a^bw''\overline w\,dx=-\int_a^b|w'|^2\,dxであり、左辺は∫abc∣w∣2 dx\int_a^bc|w|^2\,dxに等しい。従って∫ab∣w′∣2 dx+∫abc∣w∣2 dx=0\int_a^b|w'|^2\,dx+\int_a^bc|w|^2\,dx=0であり、w′=0w'=0がほとんど至る所で成り立つ。§E13.2 系 2.6によりwwは定数であり、w(a)=0w(a)=0によりw=0w=0である。従って、P=0P=0、Q=−cQ=-c、(αa,βa)=(αb,βb)=(1,0)(\alpha_a,\beta_a)=(\alpha_b,\beta_b)=(1,0)の Dirichlet 問題の境界核は自明である。

§E13.24 定理 3.4を適用して得られる連続な核をGcG_cとし、Γ=−Gc\Gamma=-G_cとおく。h(x)=A+B−Ab−a(x−a)h(x)=A+\frac{B-A}{b-a}(x-a)とし、g~=g−ch∈L1(a,b)\tilde g=g-ch\in L^1(a,b)とおく。同定理によりw′′−cw=−g~w''-cw=-\tilde g、w(a)=w(b)=0w(a)=w(b)=0の解w=∫Γg~w=\int\Gamma\tilde gがただ一つ存在する。y=h+wy=h+wとおくと、h′′=0h''=0なので−y′′+cy=−w′′+cw+ch=g~+ch=g-y''+cy=-w''+cw+ch=\tilde g+ch=g、y(a)=Ay(a)=A、y(b)=By(b)=Bである。二つの解の差は上の斉次問題の解であるから零であり、解は一意である。y=h−∫Γch+∫Γgy=h-\int\Gamma ch+\int\Gamma gであり、h−∫Γchh-\int\Gamma chはggに依らず、g=0g=0とおくとφA,B\varphi_{A,B}に一致する。

g≥0g\geq0、A,B≥0A,B\geq0とし、y=SA,Bgy=S_{A,B}gとする。y(a)=A≥0y(a)=A\geq0、y(b)=B≥0y(b)=B\geq0であり、ほとんどすべてのxxに対して、y(x)<0y(x)<0ならばy′′(x)=c(x)y(x)−g(x)≤0y''(x)=c(x)y(x)-g(x)\leq0である。従って補題 5.2によりy≥0y\geq0である。

g1≤g2g_1\leq g_2ならば、線形性によりS0,0g2−S0,0g1=S0,0(g2−g1)≥0S_{0,0}g_2-S_{0,0}g_1=S_{0,0}(g_2-g_1)\geq0である。▨

定義 5.4.a<ba<bとし、f:[a,b]×R→Rf:[a,b]\times\R\to\Rは Carathéodory 型、A,B∈RA,B\in\Rとする。関数α∈C1([a,b];R)\alpha\in C^1([a,b];\R)がα′∈AC([a,b])\alpha'\in AC([a,b])を満たし、ほとんど至る所でα′′≥f(x,α)\alpha''\geq f(x,\alpha)、かつα(a)≤A\alpha(a)\leq A、α(b)≤B\alpha(b)\leq Bを満たすとき、α\alphaを境界値問題y′′=f(x,y)y''=f(x,y)、y(a)=Ay(a)=A、y(b)=By(b)=Bの下解 (lower solution of the boundary value problem) という。関数β∈C1([a,b];R)\beta\in C^1([a,b];\R)がβ′∈AC([a,b])\beta'\in AC([a,b])を満たし、ほとんど至る所でβ′′≤f(x,β)\beta''\leq f(x,\beta)、かつβ(a)≥A\beta(a)\geq A、β(b)≥B\beta(b)\geq Bを満たすとき、β\betaを同じ問題の上解 (upper solution of the boundary value problem) という。

定理 5.5.a<ba<bとし、f:[a,b]×R→Rf:[a,b]\times\R\to\Rは Carathéodory 型であるとする。A,B∈RA,B\in\Rとし、α\alphaを下解、β\betaを上解とし、α≤β\alpha\leq\betaが[a,b][a,b]上で成り立つとする。非負関数c∈L1(a,b)c\in L^1(a,b)と Lebesgue 零集合N′⊆[a,b]N'\subseteq[a,b]が存在し、すべてのx∈[a,b]∖N′x\in[a,b]\setminus N'とα(x)≤v≤u≤β(x)\alpha(x)\leq v\leq u\leq\beta(x)を満たすu,vu,vに対して

f(x,u)−f(x,v)≤c(x)(u−v)f(x,u)-f(x,v)\leq c(x)(u-v)

が成り立つとする。F(x,u)=c(x)u−f(x,u)F(x,u)=c(x)u-f(x,u)とし、補題 5.3のSA,BS_{A,B}を用いて

α0=α,αk+1=SA,BF(⋅,αk),β0=β,βk+1=SA,BF(⋅,βk)\alpha_0=\alpha,\quad\alpha_{k+1}=S_{A,B}F(\cdot,\alpha_k),\qquad\beta_0=\beta,\quad\beta_{k+1}=S_{A,B}F(\cdot,\beta_k)

と定める。このとき、すべてのk≥0k\geq0に対して

α≤αk≤αk+1≤βk+1≤βk≤β\alpha\leq\alpha_k\leq\alpha_{k+1}\leq\beta_{k+1}\leq\beta_k\leq\beta

であり、(αk)(\alpha_k)と(βk)(\beta_k)は各点で収束する。極限y∗=lim⁡αky_*=\lim\alpha_k、y∗=lim⁡βky^*=\lim\beta_kは境界値問題y′′=f(x,y)y''=f(x,y)、y(a)=Ay(a)=A、y(b)=By(b)=Bの解であり、y∗≤y∗y_*\leq y^*である。α≤y≤β\alpha\leq y\leq\betaを満たす任意の解yyはy∗≤y≤y∗y_*\leq y\leq y^*を満たす。

証明. 連続関数qqに対して、§E13.4 補題 2.3によりF(⋅,q(⋅))F(\cdot,q(\cdot))は可測であり、∣F(x,q(x))∣≤c(x)ρ+mρ(x)|F(x,q(x))|\leq c(x)\rho+m_\rho(x)(ρ=∥q∥∞+1\rho=\|q\|_\infty+1)をほとんど至る所で満たす。従ってF(⋅,q)∈L1F(\cdot,q)\in L^1であり、αk,βk\alpha_k,\beta_kは定義される。x∉N′x\notin N'とα(x)≤v≤u≤β(x)\alpha(x)\leq v\leq u\leq\beta(x)に対してF(x,u)−F(x,v)=c(x)(u−v)−(f(x,u)−f(x,v))≥0F(x,u)-F(x,v)=c(x)(u-v)-(f(x,u)-f(x,v))\geq0であるから、F(x,⋅)F(x,\cdot)は[α(x),β(x)][\alpha(x),\beta(x)]上で単調非減少である。

k=0k=0のとき、α≤α0≤β0≤β\alpha\leq\alpha_0\leq\beta_0\leq\betaは仮定である。α≤αk≤βk≤β\alpha\leq\alpha_k\leq\beta_k\leq\betaが成り立つとする。このとき、以下の比較ではF(x,⋅)F(x,\cdot)を[α(x),β(x)][\alpha(x),\beta(x)]の中の値でだけ比べる。d=αk+1−αd=\alpha_{k+1}-\alphaとおくと、α′′≥f(x,α)\alpha''\geq f(x,\alpha)により−α′′+cα≤F(x,α)-\alpha''+c\alpha\leq F(x,\alpha)なので

−d′′+cd≥F(x,αk)−F(x,α)≥0,d(a)=A−α(a)≥0,d(b)=B−α(b)≥0-d''+cd\geq F(x,\alpha_k)-F(x,\alpha)\geq0,\qquad d(a)=A-\alpha(a)\geq0,\quad d(b)=B-\alpha(b)\geq0

であり、補題 5.3によりd≥0d\geq0である。同様にd=β−βk+1d=\beta-\beta_{k+1}は−d′′+cd≥F(x,β)−F(x,βk)≥0-d''+cd\geq F(x,\beta)-F(x,\beta_k)\geq0、d(a),d(b)≥0d(a),d(b)\geq0を満たし、βk+1≤β\beta_{k+1}\leq\betaである。d=βk+1−αk+1d=\beta_{k+1}-\alpha_{k+1}は−d′′+cd=F(x,βk)−F(x,αk)≥0-d''+cd=F(x,\beta_k)-F(x,\alpha_k)\geq0、d(a)=d(b)=0d(a)=d(b)=0を満たし、αk+1≤βk+1\alpha_{k+1}\leq\beta_{k+1}である。以上より、α≤αk+1≤βk+1≤β\alpha\leq\alpha_{k+1}\leq\beta_{k+1}\leq\betaが成り立つ。

k=0k=0のときα1≥α=α0\alpha_1\geq\alpha=\alpha_0とβ1≤β=β0\beta_1\leq\beta=\beta_0は上で示した。k≥1k\geq1でαk−1≤αk\alpha_{k-1}\leq\alpha_kとすると、d=αk+1−αkd=\alpha_{k+1}-\alpha_kは−d′′+cd=F(x,αk)−F(x,αk−1)≥0-d''+cd=F(x,\alpha_k)-F(x,\alpha_{k-1})\geq0、d(a)=d(b)=0d(a)=d(b)=0を満たし、αk+1≥αk\alpha_{k+1}\geq\alpha_kである。βk+1≤βk\beta_{k+1}\leq\beta_kも同様である。

従って各点でα≤αk↑y∗≤β\alpha\leq\alpha_k\uparrow y_*\leq\beta、α≤βk↓y∗≤β\alpha\leq\beta_k\downarrow y^*\leq\betaであり、αk≤βk\alpha_k\leq\beta_kからy∗≤y∗y_*\leq y^*である。y∗y_*は連続関数の各点極限なので可測である。補題 5.3の表示により

αk+1(x)=φA,B(x)+∫abΓ(x,s)F(s,αk(s)) ds\alpha_{k+1}(x)=\varphi_{A,B}(x)+\int_a^b\Gamma(x,s)F(s,\alpha_k(s))\,ds

である。s∉N∪N′s\notin N\cup N'(NNはffの Carathéodory 型の零集合)について、F(s,⋅)F(s,\cdot)は連続なのでF(s,αk(s))→F(s,y∗(s))F(s,\alpha_k(s))\to F(s,y_*(s))であり、∣F(s,αk(s))∣≤c(s)ρ+mρ(s)|F(s,\alpha_k(s))|\leq c(s)\rho+m_\rho(s)(ρ=max⁡{∥α∥∞,∥β∥∞}+1\rho=\max\{\|\alpha\|_\infty,\|\beta\|_\infty\}+1)である。Γ\Gammaは有界で、§E10.11 定理 3.1により

y∗(x)=φA,B(x)+∫abΓ(x,s)F(s,y∗(s)) dsy_*(x)=\varphi_{A,B}(x)+\int_a^b\Gamma(x,s)F(s,y_*(s))\,ds

である。§E13.4 補題 2.3によりF(⋅,y∗)∈L1F(\cdot,y_*)\in L^1であり、補題 5.3により右辺はSA,BF(⋅,y∗)S_{A,B}F(\cdot,y_*)である。従ってy∗∈C1y_*\in C^1、y∗′∈ACy_*'\in AC、−y∗′′+cy∗=cy∗−f(x,y∗)-y_*''+cy_*=cy_*-f(x,y_*)、すなわちy∗′′=f(x,y∗)y_*''=f(x,y_*)がほとんど至る所で成り立ち、y∗(a)=Ay_*(a)=A、y∗(b)=By_*(b)=Bである。y∗y^*も同様に解である。

α≤y≤β\alpha\leq y\leq\betaを満たす解yyを取る。−y′′+cy=F(x,y)-y''+cy=F(x,y)とy(a)=Ay(a)=A、y(b)=By(b)=Bにより、補題 5.3の一意性からy=SA,BF(⋅,y)y=S_{A,B}F(\cdot,y)である。k=0k=0のとき、α0≤y≤β0\alpha_0\leq y\leq\beta_0は仮定である。αk≤y\alpha_k\leq yならば、F(x,αk)≤F(x,y)F(x,\alpha_k)\leq F(x,y)(α≤αk≤y≤β\alpha\leq\alpha_k\leq y\leq\beta)なので、線形性と補題 5.3によりy−αk+1=S0,0(F(⋅,y)−F(⋅,αk))≥0y-\alpha_{k+1}=S_{0,0}\bigl(F(\cdot,y)-F(\cdot,\alpha_k)\bigr)\geq0である。y≤βk+1y\leq\beta_{k+1}も同様である。極限を取るとy∗≤y≤y∗y_*\leq y\leq y^*である。▨

6 事前評価による存在

定理 6.1.a<ba<bとし、f:[a,b]×R→Rf:[a,b]\times\R\to\Rは Carathéodory 型であるとする。A,B∈RA,B\in\R、R≥max⁡{∣A∣,∣B∣}R\geq\max\{|A|,|B|\}とし、Lebesgue 零集合N′⊆[a,b]N'\subseteq[a,b]が存在して、すべてのx∈[a,b]∖N′x\in[a,b]\setminus N'に対して

f(x,u)≥0 (u≥R),f(x,u)≤0 (u≤−R)f(x,u)\geq0\ (u\geq R),\qquad f(x,u)\leq0\ (u\leq-R)

が成り立つとする。このとき、境界値問題y′′=f(x,y)y''=f(x,y)、y(a)=Ay(a)=A、y(b)=By(b)=Bの任意の解yyは∣y∣≤R|y|\leq Rを満たす。

さらに、各ρ>0\rho>0に対して非負関数ℓρ∈L1(a,b)\ell_\rho\in L^1(a,b)が存在し、すべてのx∈[a,b]∖N′x\in[a,b]\setminus N'と∣u∣,∣uˉ∣≤ρ|u|,|\bar u|\leq\rhoに対して∣f(x,u)−f(x,uˉ)∣≤ℓρ(x)∣u−uˉ∣|f(x,u)-f(x,\bar u)|\leq\ell_\rho(x)|u-\bar u|が成り立つならば、この境界値問題の解が存在する。

証明.yyを解とする。w=R−yw=R-yとおくと、w∈C1w\in C^1、w′∈ACw'\in AC、w(a)=R−A≥0w(a)=R-A\geq0、w(b)=R−B≥0w(b)=R-B\geq0である。EEを、y′′(x)=f(x,y(x))y''(x)=f(x,y(x))が成り立たないx∈[a,b]x\in[a,b]全体の集合とすると、EEは Lebesgue 零集合である。x∉N′∪Ex\notin N'\cup Eかつw(x)<0w(x)<0、すなわちy(x)>Ry(x)>Rならばw′′(x)=−f(x,y(x))≤0w''(x)=-f(x,y(x))\leq0である。従って、ほとんどすべてのx∈[a,b]x\in[a,b]に対して、w(x)<0w(x)<0ならばw′′(x)≤0w''(x)\leq0である。補題 5.2によりw≥0w\geq0、すなわちy≤Ry\leq Rである。w=y+Rw=y+Rとおくと、w(a)=A+R≥0w(a)=A+R\geq0、w(b)=B+R≥0w(b)=B+R\geq0であり、x∉N′∪Ex\notin N'\cup Eかつw(x)<0w(x)<0、すなわちy(x)<−Ry(x)<-Rならばw′′(x)=f(x,y(x))≤0w''(x)=f(x,y(x))\leq0である。同じ補題によりw≥0w\geq0、すなわちy≥−Ry\geq-Rである。

ℓρ\ell_\rhoの仮定の下で、fR(x,u,v)=f(x,max⁡{−R,min⁡{R,u}})f_R(x,u,v)=f\bigl(x,\max\{-R,\min\{R,u\}\}\bigr)とおく。fRf_Rは Carathéodory 型であり、∣fR(x,u,v)∣≤mR+1(x)|f_R(x,u,v)|\leq m_{R+1}(x)(mR+1m_{R+1}はffに対するρ=R+1\rho=R+1の支配関数)、∣fR(x,u,v)−fR(x,uˉ,vˉ)∣≤ℓR+1(x)∣u−uˉ∣≤ℓR+1(x)(∣u−uˉ∣+∣v−vˉ∣)|f_R(x,u,v)-f_R(x,\bar u,\bar v)|\leq\ell_{R+1}(x)|u-\bar u|\leq\ell_{R+1}(x)\bigl(|u-\bar u|+|v-\bar v|\bigr)をすべての(u,v),(uˉ,vˉ)(u,v),(\bar u,\bar v)について満たす。従ってfRf_Rは定理 2.1の仮定を満たし、y′′=fR(x,y)y''=f_R(x,y)、y(a)=Ay(a)=A、y(b)=By(b)=Bの解yyが存在する。fR(x,⋅)f_R(x,\cdot)はf(x,⋅)f(x,\cdot)と同じ符号条件を満たすので、前半の議論をfRf_Rに適用すると∣y∣≤R|y|\leq Rである。∣y∣≤R|y|\leq Rの下でfR(x,y(x))=f(x,y(x))f_R(x,y(x))=f(x,y(x))なので、yyは元の問題の解である。▨

例 6.2 (Bratu 問題).λ>π2/e\lambda>\pi^2/eとする。このとき、f(x,u)=−λeuf(x,u)=-\lambda e^uは Carathéodory 型であり、境界値問題

y′′+λey=0,y(0)=y(1)=0y''+\lambda e^y=0,\qquad y(0)=y(1)=0

は解をもたない。このffはすべての(x,u)(x,u)でf(x,u)<0f(x,u)<0を満たすので、どのR≥0R\geq0についても定理 6.1の符号条件を満たさない。

解yyが存在すると仮定する。ほとんど至る所でy′′=−λey<0y''=-\lambda e^y<0なので、補題 5.1により[0,1][0,1]上でy≥0y\geq0である。φ(x)=sin⁡πx\varphi(x)=\sin\pi xとおく。φ\varphiは§E13.24 例 2.2のL=1L=1、n=1n=1の固有関数であり、固有値はπ2\pi^2である。§E13.2 定理 3.4をy′y'とφ\varphi、次にyyとφ′\varphi'に適用し、φ(0)=φ(1)=0\varphi(0)=\varphi(1)=0、y(0)=y(1)=0y(0)=y(1)=0を用いると

∫01y′′φ dx=−∫01y′φ′ dx=∫01yφ′′ dx=−π2∫01yφ dx\int_0^1y''\varphi\,dx=-\int_0^1y'\varphi'\,dx=\int_0^1y\varphi''\,dx=-\pi^2\int_0^1y\varphi\,dx

である。左辺は−λ∫01eyφ dx-\lambda\int_0^1e^y\varphi\,dxに等しいので、λ∫01eyφ dx=π2∫01yφ dx\lambda\int_0^1e^y\varphi\,dx=\pi^2\int_0^1y\varphi\,dxである。I=∫01yφ dxI=\int_0^1y\varphi\,dxとおく。y≥0y\geq0、φ≥0\varphi\geq0であるからI≥0I\geq0である。また、すべてのt∈Rt\in\Rに対してet≥ete^t\geq etであるから、π2I≥λeI\pi^2I\geq\lambda eIである。λe>π2\lambda e>\pi^2なのでI≤0I\leq0、従ってI=0I=0である。yφy\varphiは連続で非負なので、(0,1)(0,1)上でyφ=0y\varphi=0、すなわちy=0y=0であり、yyは[0,1][0,1]上で恒等的に零である。このときy′′=0y''=0であるが、y′′=−λe0=−λ≠0y''=-\lambda e^0=-\lambda\neq0がほとんど至る所で成り立つので、両立しない。

7 演習

問題 7.1. 境界値問題y′′=eyy''=e^y、y(0)=y(1)=0y(0)=y(1)=0の解がただ一つ存在し、その解はx(x−1)2≤y(x)≤0\frac{x(x-1)}2\leq y(x)\leq0を満たすことを示せ。

解答.

f(x,u)=euf(x,u)=e^uは Carathéodory 型である。α(x)=x(x−1)2\alpha(x)=\frac{x(x-1)}2、β=0\beta=0とおく。α∈C2\alpha\in C^2、α′′=1\alpha''=1であり、α≤0\alpha\leq0によりeα(x)≤1=α′′(x)e^{\alpha(x)}\leq1=\alpha''(x)である。α(0)=α(1)=0\alpha(0)=\alpha(1)=0なので、α\alphaは下解である。β′′=0≤e0=f(x,β)\beta''=0\leq e^0=f(x,\beta)、β(0)=β(1)=0\beta(0)=\beta(1)=0なので、β\betaは上解であり、α≤β\alpha\leq\betaである。c=1c=1とおく。α(x)≤v≤u≤0\alpha(x)\leq v\leq u\leq0に対して、平均値の定理によりeu−ev=eξ(u−v)e^u-e^v=e^\xi(u-v)(v<ξ<uv<\xi<uまたはu=vu=v)であり、eξ≤1e^\xi\leq1なのでeu−ev≤u−ve^u-e^v\leq u-vである。定理 5.5により、α≤y≤β\alpha\leq y\leq\betaを満たす解y∗y_*が存在する。

任意のu,vu,vに対して(eu−ev)(u−v)≥0(e^u-e^v)(u-v)\geq0である。従ってK=0<π2K=0<\pi^2として定理 4.2の仮定が成り立ち、解は高々一つである。存在した解y∗y_*が唯一の解であり、x(x−1)2≤y∗≤0\frac{x(x-1)}2\leq y_*\leq0を満たす。▨

問題 7.2.λ∈R\lambda\in\R、A,B∈RA,B\in\Rとし、境界値問題y′′=λsin⁡yy''=\lambda\sin y、y(0)=Ay(0)=A、y(1)=By(1)=Bを考える。この問題の解は、任意のλ\lambdaに対して存在し、∣λ∣<π2|\lambda|<\pi^2ならばただ一つであることを示せ。

解答.

f(x,u,v)=λsin⁡uf(x,u,v)=\lambda\sin uは連続で、∣f∣≤∣λ∣|f|\leq|\lambda|、∣f(x,u,v)−f(x,uˉ,vˉ)∣≤∣λ∣∣u−uˉ∣≤∣λ∣(∣u−uˉ∣+∣v−vˉ∣)|f(x,u,v)-f(x,\bar u,\bar v)|\leq|\lambda||u-\bar u|\leq|\lambda|(|u-\bar u|+|v-\bar v|)を満たす。従ってffは Carathéodory 型であり、m=∣λ∣m=|\lambda|、ℓρ=∣λ∣\ell_\rho=|\lambda|として定理 2.1の仮定を満たす。同定理により、任意のA,BA,Bに対して解が存在する。

∣sin⁡u−sin⁡v∣≤∣u−v∣|\sin u-\sin v|\leq|u-v|なので、(f(x,u)−f(x,v))(u−v)≥−∣λ∣(u−v)2(f(x,u)-f(x,v))(u-v)\geq-|\lambda|(u-v)^2である。∣λ∣<π2=(π1−0)2|\lambda|<\pi^2=\bigl(\frac\pi{1-0}\bigr)^2ならば、K=∣λ∣K=|\lambda|として定理 4.2の仮定が成り立ち、解は高々一つである。存在と合わせて、解はただ一つである。▨

参考文献

  1. Philip Hartman, Ordinary Differential Equations, 2nd ed., Classics in Applied Mathematics, Society for Industrial and Applied Mathematics, Philadelphia, 2002.射撃法、下解と上解による単調反復、事前評価による存在の定式化を参考にした。

前提記事

10 本の記事・単元を表示