§E13.26Sturm–Liouville 問題

最終更新

両端を固定した長さL>0L>0の一様な弦の振動を表す方程式を変数分離すると、空間方向の関数は

−y′′=λy,y(0)=y(L)=0-y''=\lambda y,\qquad y(0)=y(L)=0

を満たす。この問題の解として現れる正弦関数は、振動数が異なるものどうしが直交する。係数や端点条件を変えたとき、どのような条件のもとで固有値は実数となり、異なる固有値に属する解は直交するのだろうか。

個々の方程式を解いて固有関数を求めるだけでは、固有値の実性と固有関数の直交性が多くの方程式に共通して成り立つ理由は明らかにならない。二階微分方程式に重みと境界条件を加え、固有値をもつ境界値問題として定式化した問題が Sturm–Liouville 問題である。Sturm–Liouville 問題の重要性は、解の具体式に依存せず、微分式と境界条件から得られる対称性を用いて固有値と固有関数の性質を導く点にある。

本記事では、二階線形方程式の Sturm–Liouville 形への変形、端点の型、Liouville の標準形を扱い、両端が正則端点である分離型固有値問題について、流束による Lagrange の恒等式から固有値の実性と固有関数の直交性を、端点での初期値問題の一意性から固有空間の一次元性を導く。

1 流束、端点の型、分離型固有値問題

定義 1.1.−∞≤a<b≤∞-\infty\le a<b\le\inftyとし、p,q,w:(a,b)→Rp,q,w:(a,b)\to\Rを可測関数で、ほとんどすべてのx∈(a,b)x\in(a,b)でp(x)>0p(x)>0かつw(x)>0w(x)>0を満たし、1/p,q,w1/p,q,wが(a,b)(a,b)上で局所可積分であるものとする。y∈ACloc((a,b);C)y\in AC_{\mathrm{loc}}((a,b);\C)について、py′py'とほとんど至る所で一致する局所絶対連続関数(a,b)→C(a,b)\to\Cが存在するとき、その関数をyyの流束 (flux) と呼び、y[1]y^{[1]}と書く。流束は存在すれば一つに限る。実際、二つの連続関数がほとんど至る所で一致すれば、値が異なる点の集合は測度零の開集合であるから空である。流束をもつyyに対して

ℓ[y]=−(y[1])′+qy\ell[y]=-\bigl(y^{[1]}\bigr)'+qy

とおく。右辺はほとんど至る所で定まる。λ∈C\lambda\in\Cと、wfwfが(a,b)(a,b)上で局所可積分である可測関数f:(a,b)→Cf:(a,b)\to\Cについて、方程式ℓ[y]=λwy+wf\ell[y]=\lambda wy+wfの解とは、流束をもつy∈ACloc((a,b);C)y\in AC_{\mathrm{loc}}((a,b);\C)であって、ℓ[y]=λwy+wf\ell[y]=\lambda wy+wfがほとんど至る所で成り立つものをいう。

定義 1.2.−∞≤a<b≤∞-\infty\le a<b\le\inftyとし、p,q,wp,q,wを(a,b)(a,b)上で定義 1.1の仮定を満たす関数とする。端点aaが正則端点 (regular endpoint) であるとは、a>−∞a>-\inftyであり、(a,c)(a,c)上で1/p,q,w1/p,q,wが可積分であるようなc∈(a,b)c\in(a,b)が存在することをいう。この条件はc∈(a,b)c\in(a,b)の取り方によらない。実際、c<c′c<c'ならば[c,c′][c,c']は(a,b)(a,b)に含まれる有限閉区間であり、1/p,q,w1/p,q,wはその上で可積分である。端点bbが正則端点であるとは、b<∞b<\inftyであり、(c,b)(c,b)上で1/p,q,w1/p,q,wが可積分であるようなc∈(a,b)c\in(a,b)が存在することをいう。正則端点でない端点を特異端点 (singular endpoint) という。

命題 1.3.−∞≤a<b≤∞-\infty\le a<b\le\inftyとし、p,q,wp,q,wを(a,b)(a,b)上で定義 1.1の仮定を満たす関数とする。λ∈C\lambda\in\Cとし、f:(a,b)→Cf:(a,b)\to\Cをwfwfが(a,b)(a,b)上で局所可積分である可測関数とする。

Aλ=(01/pq−λw0),b=(0−wf)A_\lambda=\begin{pmatrix}0&1/p\\ q-\lambda w&0\end{pmatrix},\qquad \boldsymbol b=\begin{pmatrix}0\\ -wf\end{pmatrix}

とおく。このとき、次が成り立つ。

  1. 関数y:(a,b)→Cy:(a,b)\to\Cがℓ[y]=λwy+wf\ell[y]=\lambda wy+wfの解であるための必要十分条件は、Y=(Y0,Y1)T∈ACloc((a,b);C2)\boldsymbol Y=(Y_0,Y_1)^{\mathsf T}\in AC_{\mathrm{loc}}((a,b);\C^2)であって、Y′=AλY+b\boldsymbol Y'=A_\lambda\boldsymbol Y+\boldsymbol bをほとんど至る所で満たし、Y0=yY_0=yとなるものが存在することである。このときY1=y[1]Y_1=y^{[1]}である。
  2. 任意のc∈(a,b)c\in(a,b)と(η0,η1)∈C2(\eta_0,\eta_1)\in\C^2に対して、y(c)=η0y(c)=\eta_0かつy[1](c)=η1y^{[1]}(c)=\eta_1を満たすℓ[y]=λwy+wf\ell[y]=\lambda wy+wfの解yyがただ一つ存在する。λ,η0,η1\lambda,\eta_0,\eta_1が実数であり、ffが実数値であるならば、この解は実数値である。
  3. u,vu,vがℓ[y]=λwy\ell[y]=\lambda wyの解ならば、W(u,v)=uv[1]−u[1]vW(u,v)=uv^{[1]}-u^{[1]}vは(a,b)(a,b)上で定数である。
  4. e∈{a,b}e\in\{a,b\}を正則端点とし、J=(a,b)∪{e}J=(a,b)\cup\{e\}とおく。eeとccを端点とする開区間の上でwfwfが可積分であるようなc∈(a,b)c\in(a,b)が存在するとする。このときℓ[y]=λwy+wf\ell[y]=\lambda wy+wfの任意の解yyについて、yyとy[1]y^{[1]}はJJ上の局所絶対連続関数に延びる。延長のeeでの値をy(e)y(e)、y[1](e)y^{[1]}(e)と書く。任意の(η0,η1)∈C2(\eta_0,\eta_1)\in\C^2に対して、y(e)=η0y(e)=\eta_0かつy[1](e)=η1y^{[1]}(e)=\eta_1を満たすℓ[y]=λwy+wf\ell[y]=\lambda wy+wfの解yyがただ一つ存在する。
  5. ppが(a,b)(a,b)上の正の連続関数ならば、流束をもつ任意のy∈ACloc((a,b);C)y\in AC_{\mathrm{loc}}((a,b);\C)はC1C^1級であり、すべてのx∈(a,b)x\in(a,b)でy[1](x)=p(x)y′(x)y^{[1]}(x)=p(x)y'(x)が成り立つ。とくにℓ[y]=λwy\ell[y]=\lambda wyの解u,vu,vについて、W(u,v)W(u,v)はppと Wronskianuv′−u′vuv'-u'vの積である。

証明.主張 (1)を示す。yyが解であるとし、Y=(y,y[1])T\boldsymbol Y=(y,y^{[1]})^{\mathsf T}とおく。Y\boldsymbol Yは局所絶対連続である。y[1]=py′y^{[1]}=py'がほとんど至る所で成り立ち、ほとんどすべての点でp>0p>0であるから、y′=y[1]/py'=y^{[1]}/pがほとんど至る所で成り立つ。ℓ[y]\ell[y]の定義から(y[1])′=(q−λw)y−wf(y^{[1]})'=(q-\lambda w)y-wfがほとんど至る所で成り立つ。したがってY′=AλY+b\boldsymbol Y'=A_\lambda\boldsymbol Y+\boldsymbol bがほとんど至る所で成り立つ。逆にY\boldsymbol Yが条件を満たすとすると、Y0′=Y1/pY_0'=Y_1/pがほとんど至る所で成り立つからpY0′=Y1pY_0'=Y_1がほとんど至る所で成り立ち、Y1Y_1は局所絶対連続であるからY1Y_1はY0Y_0の流束である。Y1′=(q−λw)Y0−wfY_1'=(q-\lambda w)Y_0-wfがほとんど至る所で成り立つから、ℓ[Y0]=−Y1′+qY0=λwY0+wf\ell[Y_0]=-Y_1'+qY_0=\lambda wY_0+wfがほとんど至る所で成り立つ。

主張 (2)を示す。1/p,q,w,wf1/p,q,w,wfは(a,b)(a,b)上で局所可積分であるから、Aλ∈Lloc1((a,b);C2×2)A_\lambda\in L^1_{\mathrm{loc}}((a,b);\C^{2\times2})、b∈Lloc1((a,b);C2)\boldsymbol b\in L^1_{\mathrm{loc}}((a,b);\C^2)である。§E13.9 定理 1.1をI=(a,b)I=(a,b)、K=C\mathbb K=\Cに適用すると、Y(c)=(η0,η1)T\boldsymbol Y(c)=(\eta_0,\eta_1)^{\mathsf T}を満たすY′=AλY+b\boldsymbol Y'=A_\lambda\boldsymbol Y+\boldsymbol bの局所絶対連続解がただ一つ存在する。主張 (1)により、ℓ[y]=λwy+wf\ell[y]=\lambda wy+wfの解yyと、この系の解Y\boldsymbol YはY=(y,y[1])T\boldsymbol Y=(y,y^{[1]})^{\mathsf T}によって一対一に対応するから、求める解はただ一つ存在する。λ,η0,η1\lambda,\eta_0,\eta_1が実数でffが実数値ならば、AλA_\lambdaとb\boldsymbol bは実数値であり、§E13.9 定理 1.1をK=R\mathbb K=\Rに適用して得られる実数値の解はK=C\mathbb K=\Cの解でもあるから、一意性により求める解に一致する。

主張 (3)を示す。行列値関数

Φ=(uvu[1]v[1])\Phi=\begin{pmatrix}u&v\\ u^{[1]}&v^{[1]}\end{pmatrix}

は、主張 (1)をf=0f=0に適用すると、Φ′=AλΦ\Phi'=A_\lambda\Phiをほとんど至る所で満たす局所絶対連続関数である。tr⁡Aλ=0\operatorname{tr}A_\lambda=0であるから、§E13.9 定理 3.1によりdet⁡Φ=uv[1]−u[1]v\det\Phi=uv^{[1]}-u^{[1]}vは定数である。

主張 (4)を示す。e=ae=aの場合、1/p,q,w,wf1/p,q,w,wfは(a,c)(a,c)上で可積分であり、JJに含まれる任意の有限閉区間KKについてK⊆[a,c′]K\subseteq[a,c']を満たすc′∈(c,b)c'\in(c,b)が存在するから、Aλ∈Lloc1(J;C2×2)A_\lambda\in L^1_{\mathrm{loc}}(J;\C^{2\times2})、b∈Lloc1(J;C2)\boldsymbol b\in L^1_{\mathrm{loc}}(J;\C^2)である。e=be=bの場合も、JJに含まれる任意の有限閉区間KKについてK⊆[c′,b]K\subseteq[c',b]を満たすc′∈(a,c)c'\in(a,c)が存在するから、同じ結論が成り立つ。yyを解とし、Y=(y,y[1])T\boldsymbol Y=(y,y^{[1]})^{\mathsf T}とおく。§E13.9 定理 1.1をI=JI=Jに適用すると、Z(c)=Y(c)\boldsymbol Z(c)=\boldsymbol Y(c)を満たすZ′=AλZ+b\boldsymbol Z'=A_\lambda\boldsymbol Z+\boldsymbol bの解Z∈ACloc(J;C2)\boldsymbol Z\in AC_{\mathrm{loc}}(J;\C^2)が存在する。Z\boldsymbol Zの(a,b)(a,b)への制限とY\boldsymbol Yは、§E13.9 定理 1.1をI=(a,b)I=(a,b)に適用したときの同じ初期値の解であるから一致する。したがってyyとy[1]y^{[1]}はJJ上の局所絶対連続関数Z0,Z1Z_0,Z_1に延びる。Z(e)=(η0,η1)T\boldsymbol Z(e)=(\eta_0,\eta_1)^{\mathsf T}を満たすJJ上の解Z\boldsymbol Zの(a,b)(a,b)への制限は、主張 (1)により解y=Z0y=Z_0を与え、y(e)=η0y(e)=\eta_0、y[1](e)=η1y^{[1]}(e)=\eta_1を満たす。y,y~y,\tilde yをy(e)=y~(e)y(e)=\tilde y(e)、y[1](e)=y~[1](e)y^{[1]}(e)=\tilde y^{[1]}(e)を満たす二つの解とすると、(y,y[1])(y,y^{[1]})と(y~,y~[1])(\tilde y,\tilde y^{[1]})のJJへの延長は、ともにJJ上でZ′=AλZ+b\boldsymbol Z'=A_\lambda\boldsymbol Z+\boldsymbol bをほとんど至る所で満たしeeで同じ値をとる局所絶対連続関数であるから、§E13.9 定理 1.1の一意性により一致する。

主張 (5)を示す。h=y[1]/ph=y^{[1]}/pは(a,b)(a,b)上の連続関数であり、y′=hy'=hがほとんど至る所で成り立つ。§E13.2 系 2.7 (1)をI=(a,b)I=(a,b)に適用すると、yyはC1C^1級であり、すべての点でy′=hy'=h、すなわちy[1]=py′y^{[1]}=py'が成り立つ。u,vu,vに適用するとW(u,v)=u pv′−pu′ v=p(uv′−u′v)W(u,v)=u\,pv'-pu'\,v=p(uv'-u'v)である。▨

定理 1.4.−∞≤a<b≤∞-\infty\le a<b\le\inftyとし、p,q,wp,q,wを(a,b)(a,b)上で定義 1.1の仮定を満たす関数とする。c,dc,dをa≤c<d≤ba\le c<d\le bを満たす実数とし、c=ac=aならばaaは正則端点、d=bd=bならばbbは正則端点であるとする。u,v∈ACloc((a,b);C)u,v\in AC_{\mathrm{loc}}((a,b);\C)は流束をもち、u,u[1],v,v[1]u,u^{[1]},v,v^{[1]}の(c,d)(c,d)への制限は[c,d][c,d]上の絶対連続関数に延びるとする。延長のc,dc,dでの値をu(c),u[1](c)u(c),u^{[1]}(c)などと書き、

H(x)=u(x)v[1](x)‾−u[1](x)v(x)‾(x∈[c,d])H(x)=u(x)\overline{v^{[1]}(x)}-u^{[1]}(x)\overline{v(x)}\qquad(x\in[c,d])

とおく。このときℓ[u]v‾\ell[u]\overline vとuℓ[v]‾u\overline{\ell[v]}は(c,d)(c,d)上で可積分であり、次が成り立つ。

  1. 次の等式が成り立つ。 ∫cd(ℓ[u]v‾−uℓ[v]‾) dx=H(d)−H(c).\int_c^d\bigl(\ell[u]\overline v-u\overline{\ell[v]}\bigr)\,dx=H(d)-H(c).
  2. e∈{c,d}e\in\{c,d\}とし、(α,β)∈R2∖{(0,0)}(\alpha,\beta)\in\R^2\setminus\{(0,0)\}についてαu(e)+βu[1](e)=0\alpha u(e)+\beta u^{[1]}(e)=0かつαv(e)+βv[1](e)=0\alpha v(e)+\beta v^{[1]}(e)=0が成り立つならば、H(e)=0H(e)=0である。
  3. ε∈{1,−1}\varepsilon\in\{1,-1\}について(u(d),u[1](d))=ε(u(c),u[1](c))(u(d),u^{[1]}(d))=\varepsilon(u(c),u^{[1]}(c))かつ(v(d),v[1](d))=ε(v(c),v[1](c))(v(d),v^{[1]}(d))=\varepsilon(v(c),v^{[1]}(c))が成り立つならば、H(d)−H(c)=0H(d)-H(c)=0である。

証明.c<k<dc<k<dを満たすkkを取る。c=ac=aならばaaは正則端点であるから定義 1.2によりqqは(c,k)(c,k)上で可積分であり、c>ac>aならば[c,k][c,k]は(a,b)(a,b)に含まれる有限閉区間であるからqqは(c,k)(c,k)上で可積分である。d=bd=bの場合とd<bd<bの場合に分けると、同じ理由でqqは(k,d)(k,d)上で可積分である。したがってqqは(c,d)(c,d)上で可積分である。u,vu,vの延長は[c,d][c,d]上で有界であるからquv‾qu\overline vとuqv‾u\overline{qv}は(c,d)(c,d)上で可積分である。(u[1])′(u^{[1]})'と(v[1])′(v^{[1]})'は§E13.2 定理 2.5により(c,d)(c,d)上で可積分であるから、ℓ[u]v‾\ell[u]\overline vとuℓ[v]‾u\overline{\ell[v]}は(c,d)(c,d)上で可積分である。

主張 (1)を示す。複素共役z↦z‾z\mapsto\overline zは実変数のC1C^1写像であるから、§E13.2 定理 3.1 (1)によりv‾,v[1]‾\overline v,\overline{v^{[1]}}は[c,d][c,d]上で絶対連続であり、その導関数はv′‾,(v[1])′‾\overline{v'},\overline{(v^{[1]})'}である。B((z0,z1),(ζ0,ζ1))=z0ζ1−z1ζ0B((z_0,z_1),(\zeta_0,\zeta_1))=z_0\zeta_1-z_1\zeta_0は双線形写像C2×C2→C\C^2\times\C^2\to\Cであり、H=B((u,u[1]),(v‾,v[1]‾))H=B\bigl((u,u^{[1]}),(\overline v,\overline{v^{[1]}})\bigr)であるから、§E13.2 定理 3.1 (2)によりHHは[c,d][c,d]上で絶対連続であり、ほとんど至る所で

H′=u′v[1]‾+u(v[1])′‾−(u[1])′v‾−u[1]v′‾H'=u'\overline{v^{[1]}}+u\overline{(v^{[1]})'}-(u^{[1]})'\overline v-u^{[1]}\overline{v'}

が成り立つ。u[1]=pu′u^{[1]}=pu'、v[1]=pv′v^{[1]}=pv'がほとんど至る所で成り立ちppは実数値であるから、u′v[1]‾=pu′v′‾=u[1]v′‾u'\overline{v^{[1]}}=pu'\overline{v'}=u^{[1]}\overline{v'}がほとんど至る所で成り立つ。qqは実数値であるからquv‾=uqv‾qu\overline v=u\overline{qv}である。したがってほとんど至る所で

H′=u(v[1])′‾−(u[1])′v‾=(ℓ[u]−qu)v‾−u(ℓ[v]−qv)‾=ℓ[u]v‾−uℓ[v]‾H'=u\overline{(v^{[1]})'}-(u^{[1]})'\overline v =\bigl(\ell[u]-qu\bigr)\overline v-u\overline{\bigl(\ell[v]-qv\bigr)} =\ell[u]\overline v-u\overline{\ell[v]}

が成り立つ。§E13.2 定理 2.5 (3)によりH(d)−H(c)=∫cdH′ dxH(d)-H(c)=\int_c^dH'\,dxである。

主張 (2)を示す。α,β\alpha,\betaは実数であるからαv(e)‾+βv[1](e)‾=αv(e)+βv[1](e)‾=0\alpha\overline{v(e)}+\beta\overline{v^{[1]}(e)}=\overline{\alpha v(e)+\beta v^{[1]}(e)}=0である。二つのベクトル(u(e),u[1](e))(u(e),u^{[1]}(e))と(v(e)‾,v[1](e)‾)(\overline{v(e)},\overline{v^{[1]}(e)})は、零でない線形形式(z0,z1)↦αz0+βz1(z_0,z_1)\mapsto\alpha z_0+\beta z_1の核に属する。この核はC2\C^2の一次元部分空間であるから、二つのベクトルは一次従属であり、行列式H(e)=u(e)v[1](e)‾−u[1](e)v(e)‾H(e)=u(e)\overline{v^{[1]}(e)}-u^{[1]}(e)\overline{v(e)}は零である。

主張 (3)を示す。ε\varepsilonは実数でε2=1\varepsilon^2=1であるから、H(d)=ε2(u(c)v[1](c)‾−u[1](c)v(c)‾)=H(c)H(d)=\varepsilon^2\bigl(u(c)\overline{v^{[1]}(c)}-u^{[1]}(c)\overline{v(c)}\bigr)=H(c)である。▨

補題 1.5.a<ba<bを実数とし、p,q,wp,q,wを(a,b)(a,b)上で定義 1.1の仮定を満たす関数で、a,ba,bがともに正則端点であるものとする。λ∈C\lambda\in\Cとする。このとき、次が成り立つ。

  1. 1/p,q,w1/p,q,wは(a,b)(a,b)上で可積分である。
  2. ℓ[y]=λwy\ell[y]=\lambda wyの任意の解yyについて、yyとy[1]y^{[1]}は[a,b][a,b]上の絶対連続関数に延びる。延長のa,ba,bでの値をy(a),y[1](a),y(b),y[1](b)y(a),y^{[1]}(a),y(b),y^{[1]}(b)と書く。
  3. 任意のc∈[a,b]c\in[a,b]と(η0,η1)∈C2(\eta_0,\eta_1)\in\C^2に対して、y(c)=η0y(c)=\eta_0かつy[1](c)=η1y^{[1]}(c)=\eta_1を満たすℓ[y]=λwy\ell[y]=\lambda wyの解yyがただ一つ存在する。λ,η0,η1\lambda,\eta_0,\eta_1が実数であるならば、この解は実数値である。
  4. ℓ[y]=λwy\ell[y]=\lambda wyの零でない解yyについて、すべてのx∈[a,b]x\in[a,b]で(y(x),y[1](x))≠(0,0)(y(x),y^{[1]}(x))\ne(0,0)である。
  5. ℓ[y]=λwy\ell[y]=\lambda wyの解u,vu,vについて、W(u,v)=uv[1]−u[1]vW(u,v)=uv^{[1]}-u^{[1]}vは[a,b][a,b]上で定数であり、u,vu,vが一次独立であることとW(u,v)≠0W(u,v)\ne0は同値である。

証明.主張 (1)を示す。c∈(a,b)c\in(a,b)を取る。定義 1.2により、aaが正則端点であるから1/p,q,w1/p,q,wは(a,c)(a,c)上で可積分であり、bbが正則端点であるから(c,b)(c,b)上で可積分である。したがって1/p,q,w1/p,q,wは(a,b)(a,b)上で可積分である。

主張 (2)を示す。yyを解とし、c∈(a,b)c\in(a,b)を取る。e∈{a,b}e\in\{a,b\}について、命題 1.3 (4)をf=0f=0に適用すると、yyとy[1]y^{[1]}は(a,b)∪{e}(a,b)\cup\{e\}上の局所絶対連続関数に延びる。e=ae=aの延長は[a,c][a,c]上で絶対連続であり、e=be=bの延長は[c,b][c,b]上で絶対連続であり、両者はccで一致する。[a,c][a,c]上でe=ae=aの延長に、[c,b][c,b]上でe=be=bの延長に等しい関数は、(a,b)(a,b)上でyyまたはy[1]y^{[1]}に一致し、§E13.2 命題 2.8により[a,b][a,b]上で絶対連続である。

主張 (3)を示す。c∈(a,b)c\in(a,b)のとき、命題 1.3 (2)をf=0f=0に適用すると、条件を満たす解がただ一つ存在する。c∈{a,b}c\in\{a,b\}のとき、命題 1.3 (4)をe=ce=c、f=0f=0に適用すると、条件を満たす解がただ一つ存在する。λ,η0,η1\lambda,\eta_0,\eta_1が実数であるとし、yyをこの解とする。複素共役z↦z‾z\mapsto\overline zは実変数のC1C^1写像であるから、§E13.2 定理 3.1 (1)によりy‾\overline yとy[1]‾\overline{y^{[1]}}は(a,b)(a,b)上で局所絶対連続であり、その導関数はほとんど至る所でy′‾\overline{y'}と(y[1])′‾\overline{(y^{[1]})'}である。p,q,w,λp,q,w,\lambdaは実数値であるから、ほとんど至る所でy[1]‾=py′‾\overline{y^{[1]}}=p\overline{y'}かつ−(y[1]‾)′+qy‾=ℓ[y]‾=λwy‾-\bigl(\overline{y^{[1]}}\bigr)'+q\overline y=\overline{\ell[y]}=\lambda w\overline yが成り立つ。したがってy‾\overline yは流束y[1]‾\overline{y^{[1]}}をもつℓ[y]=λwy\ell[y]=\lambda wyの解であり、y‾\overline yとy[1]‾\overline{y^{[1]}}の[a,b][a,b]への延長はyyとy[1]y^{[1]}の延長の複素共役である。y‾(c)=η0‾=η0\overline y(c)=\overline{\eta_0}=\eta_0、y[1]‾(c)=η1‾=η1\overline{y^{[1]}}(c)=\overline{\eta_1}=\eta_1であるから、一意性によりy‾=y\overline y=yであり、yyは実数値である。

主張 (4)を示す。零でない解yyとx∈[a,b]x\in[a,b]について(y(x),y[1](x))=(0,0)(y(x),y^{[1]}(x))=(0,0)であるとすると、零関数もxxで同じ値をとる解であるから、主張 (3)の一意性によりy=0y=0となり、yyが零でないことに反する。

主張 (5)を示す。流束は存在すれば一つに限るから、解u,vu,vとk1,k2∈Ck_1,k_2\in\Cについてk1u+k2vk_1u+k_2vは流束k1u[1]+k2v[1]k_1u^{[1]}+k_2v^{[1]}をもつ解である。命題 1.3 (3)によりW(u,v)W(u,v)は(a,b)(a,b)上で定数であり、主張 (2)の延長によって[a,b][a,b]上の連続関数に延びるから、[a,b][a,b]上で定数である。(k1,k2)≠(0,0)(k_1,k_2)\ne(0,0)についてk1u+k2v=0k_1u+k_2v=0ならば、その流束k1u[1]+k2v[1]k_1u^{[1]}+k_2v^{[1]}も零であるから、行列(uvu[1]v[1])\begin{pmatrix}u&v\\u^{[1]}&v^{[1]}\end{pmatrix}の列は各点で一次従属であり、W(u,v)=0W(u,v)=0である。逆にW(u,v)=0W(u,v)=0とし、c∈[a,b]c\in[a,b]を取ると、(k1,k2)≠(0,0)(k_1,k_2)\ne(0,0)が存在してk1u(c)+k2v(c)=0k_1u(c)+k_2v(c)=0かつk1u[1](c)+k2v[1](c)=0k_1u^{[1]}(c)+k_2v^{[1]}(c)=0が成り立つ。k1u+k2vk_1u+k_2vはccでの値が(0,0)(0,0)の解であるから、主張 (3)の一意性により零であり、u,vu,vは一次従属である。▨

定義 1.6.a<ba<bを実数とし、p,q,wp,q,wを(a,b)(a,b)上で定義 1.1の仮定を満たす関数で、a,ba,bがともに正則端点であるものとする。(αa,βa),(αb,βb)∈R2∖{(0,0)}(\alpha_a,\beta_a),(\alpha_b,\beta_b)\in\R^2\setminus\{(0,0)\}とする。λ∈C\lambda\in\Cについて、ℓ[y]=λwy\ell[y]=\lambda wyの零でない解yyであって、補題 1.5 (2)の延長の値について

αay(a)+βay[1](a)=0,αby(b)+βby[1](b)=0\alpha_ay(a)+\beta_ay^{[1]}(a)=0,\qquad \alpha_by(b)+\beta_by^{[1]}(b)=0

を満たすものが存在するとき、λ\lambdaをこの境界条件の下の固有値 (Sturm–Liouville eigenvalue)、yyをλ\lambdaに属する固有関数 (Sturm–Liouville eigenfunction) という。固有値λ\lambdaに属する固有関数の全体に零関数を加えた集合を、λ\lambdaの固有空間 (eigenspace) という。

命題 1.7.a<ba<bを実数とし、p,q,wp,q,wを(a,b)(a,b)上で定義 1.1の仮定を満たす関数で、a,ba,bがともに正則端点であるものとする。(αa,βa),(αb,βb)∈R2∖{(0,0)}(\alpha_a,\beta_a),(\alpha_b,\beta_b)\in\R^2\setminus\{(0,0)\}とし、固有値、固有関数、固有空間はこの境界条件の下の定義 1.6の意味のものとする。このとき、次が成り立つ。

  1. すべての固有値は実数である。
  2. λ≠μ\lambda\ne\muを固有値とし、uuをλ\lambdaに属する固有関数、vvをμ\muに属する固有関数とする。このとき、wuv‾wu\overline vは(a,b)(a,b)上で可積分であり、∫abuv‾w dx=0\displaystyle\int_a^bu\overline vw\,dx=0である。
  3. 各固有値の固有空間は、実数値の固有関数一つで張られる一次元の複素ベクトル空間である。

証明.λ,μ\lambda,\muを固有値、uuをλ\lambdaに属する固有関数、vvをμ\muに属する固有関数とする。補題 1.5 (2)によりu,u[1],v,v[1]u,u^{[1]},v,v^{[1]}は[a,b][a,b]上の絶対連続関数に延びるから、定理 1.4をc=ac=a、d=bd=bに適用することができる。補題 1.5 (1)によりwwは(a,b)(a,b)上で可積分であり、u,vu,vの延長は[a,b][a,b]上で有界であるから、wuv‾wu\overline vは(a,b)(a,b)上で可積分である。wwは実数値であるから、ほとんど至る所でℓ[u]v‾−uℓ[v]‾=(λ−μ‾)wuv‾\ell[u]\overline v-u\overline{\ell[v]}=(\lambda-\overline\mu)wu\overline vが成り立ち、定理 1.4 (1)により(λ−μ‾)∫abwuv‾ dx=H(b)−H(a)(\lambda-\overline\mu)\int_a^bwu\overline v\,dx=H(b)-H(a)である。e∈{a,b}e\in\{a,b\}について、uuとvvはともにαez(e)+βez[1](e)=0\alpha_ez(e)+\beta_ez^{[1]}(e)=0を満たすから、定理 1.4 (2)によりH(e)=0H(e)=0である。したがって

(λ−μ‾)∫abwuv‾ dx=0(\lambda-\overline\mu)\int_a^bwu\overline v\,dx=0

が成り立つ。

主張 (1)を示す。(λ−μ‾)∫abwuv‾ dx=0(\lambda-\overline\mu)\int_a^bwu\overline v\,dx=0にμ=λ\mu=\lambda、v=uv=uを代入すると、(λ−λ‾)∫abw∣u∣2 dx=0(\lambda-\overline\lambda)\int_a^bw\lvert u\rvert^2\,dx=0である。uuは連続で零でないから、∣u∣>0\lvert u\rvert>0がその上で成り立つ空でない開区間が存在し、w>0w>0がほとんど至る所で成り立つから、w∣u∣2w\lvert u\rvert^2は正の測度をもつ集合の上で正である。§E10.10 命題 3.4により∫abw∣u∣2 dx>0\int_a^bw\lvert u\rvert^2\,dx>0であり、λ=λ‾\lambda=\overline\lambdaである。

主張 (2)を示す。主張 (1)によりμ‾=μ\overline\mu=\muであるから、(λ−μ‾)∫abwuv‾ dx=0(\lambda-\overline\mu)\int_a^bwu\overline v\,dx=0は(λ−μ)∫abwuv‾ dx=0(\lambda-\mu)\int_a^bwu\overline v\,dx=0となる。したがって、λ≠μ\lambda\ne\muより積分は零である。

主張 (3)を示す。λ\lambdaを固有値とすると、主張 (1)によりλ\lambdaは実数である。(αa,βa)(\alpha_a,\beta_a)は零でない実数の組であるから、線形形式(z0,z1)↦αaz0+βaz1(z_0,z_1)\mapsto\alpha_az_0+\beta_az_1のC2\C^2における核は、零でない実ベクトル(−βa,αa)(-\beta_a,\alpha_a)で張られる一次元部分空間である。補題 1.5 (3)により、yλ(a)=−βay_\lambda(a)=-\beta_a、yλ[1](a)=αay_\lambda^{[1]}(a)=\alpha_aを満たすℓ[y]=λwy\ell[y]=\lambda wyの解yλy_\lambdaがただ一つ存在し、yλy_\lambdaは実数値で零でない。yyをλ\lambdaに属する固有関数とすると、(y(a),y[1](a))(y(a),y^{[1]}(a))は核に属するから、(y(a),y[1](a))=k(−βa,αa)(y(a),y^{[1]}(a))=k(-\beta_a,\alpha_a)を満たすk∈Ck\in\Cが存在する。y−kyλy-ky_\lambdaは流束y[1]−kyλ[1]y^{[1]}-ky_\lambda^{[1]}をもちaaでの値が(0,0)(0,0)の解であるから、補題 1.5 (3)の一意性によりy=kyλy=ky_\lambdaであり、y≠0y\ne0からk≠0k\ne0である。したがってyλ=k−1yy_\lambda=k^{-1}yはbbでの境界条件を満たし、λ\lambdaに属する実数値の固有関数である。固有空間の元はyλy_\lambdaの複素数倍であり、yλy_\lambdaの零でない複素数倍はすべてλ\lambdaに属する固有関数であるから、固有空間はCyλ\C y_\lambdaに等しい。▨

補題 1.8.IIを開区間とし、P,Q,R,W:I→RP,Q,R,W:I\to\Rを可測関数で、すべてのx∈Ix\in IでP(x)≠0P(x)\ne0を満たし、Q/P,R/P,W/PQ/P,R/P,W/PがII上で局所可積分であるものとする。c∈Ic\in Iを固定し、

p(x)=exp⁡(∫cxQ(s)P(s) ds),q=−pRP,w=pWPp(x)=\exp\left(\int_c^x\frac{Q(s)}{P(s)}\,ds\right),\qquad q=-\frac{pR}{P},\qquad w=\frac{pW}{P}

とおく。このとき、次が成り立つ。

  1. ppと1/p1/pはII上で局所絶対連続であり、すべての点でp>0p>0である。q,wq,wはII上で局所可積分である。ほとんどすべての点でW/P>0W/P>0ならば、p,q,wp,q,wはII上で定義 1.1の仮定を満たす。
  2. λ∈C\lambda\in\Cとし、y:I→Cy:I\to\Cを関数とする。y∈C1(I;C)y\in C^1(I;\C)、y′∈ACloc(I;C)y'\in AC_{\mathrm{loc}}(I;\C)であり y′′+QPy′+R+λWPy=0ほとんど至る所y''+\frac QPy'+\frac{R+\lambda W}Py=0\quad\text{ほとんど至る所} が成り立つことと、y∈ACloc(I;C)y\in AC_{\mathrm{loc}}(I;\C)であり、py′py'とほとんど至る所で一致する局所絶対連続関数y[1]y^{[1]}が存在して−(y[1])′+qy=λwy-(y^{[1]})'+qy=\lambda wyがほとんど至る所で成り立つことは同値である。
  3. ccをc′∈Ic'\in Iに替えると、p,q,wp,q,wはいずれも正の定数exp⁡(∫c′cQ(s)/P(s) ds)\exp\bigl(\int_{c'}^cQ(s)/P(s)\,ds\bigr)倍になる。

証明.主張 (1)を示す。g(x)=∫cxQ(s)/P(s) dsg(x)=\int_c^xQ(s)/P(s)\,dsとおくと、§E13.2 系 2.6によりggは局所絶対連続であり、g′=Q/Pg'=Q/Pがほとんど至る所で成り立つ。exp⁡\expはC1C^1写像であるから、§E13.2 定理 3.1 (1)によりp=exp⁡∘gp=\exp\circ gと1/p=exp⁡∘(−g)1/p=\exp\circ(-g)は局所絶対連続であり、p′=pQ/Pp'=pQ/Pがほとんど至る所で成り立つ。ppは正の連続関数であるからIIに含まれる各有限閉区間の上で有界であり、q=−pR/Pq=-pR/Pとw=pW/Pw=pW/Pは局所可積分である。1/p1/pは連続であるから局所可積分である。ほとんどすべての点でW/P>0W/P>0ならば、ほとんどすべての点でw>0w>0であり、定義 1.1の仮定が満たされる。

主張 (2)を示す。y∈C1(I;C)y\in C^1(I;\C)、y′∈ACloc(I;C)y'\in AC_{\mathrm{loc}}(I;\C)が方程式を満たすとする。yyはIIに含まれる各有限閉区間の上で導関数が有界であるから Lipschitz 連続であり、§E13.2 補題 1.2 (4)により局所絶対連続である。§E13.2 定理 3.1 (2)によりpy′py'は局所絶対連続であり、ほとんど至る所で

(py′)′=p′y′+py′′=p(y′′+QPy′)=−p R+λWPy=(q−λw)y(py')'=p'y'+py''=p\Bigl(y''+\frac QPy'\Bigr)=-p\,\frac{R+\lambda W}Py=(q-\lambda w)y

が成り立つ。したがってy[1]=py′y^{[1]}=py'について−(y[1])′+qy=λwy-(y^{[1]})'+qy=\lambda wyがほとんど至る所で成り立つ。

逆に、y∈ACloc(I;C)y\in AC_{\mathrm{loc}}(I;\C)とy[1]y^{[1]}が後者の条件を満たすとする。p,0,1p,0,1はII上で定義 1.1の仮定を満たし、y[1]y^{[1]}はその意味でのyyの流束であるから、命題 1.3 (5)によりy∈C1(I;C)y\in C^1(I;\C)であり、すべての点でy′=y[1]/py'=y^{[1]}/pである。§E13.2 定理 3.1 (2)によりy′=y[1]⋅(1/p)y'=y^{[1]}\cdot(1/p)は局所絶対連続であり、(1/p)′=−p′/p2=−(Q/P)/p(1/p)'=-p'/p^2=-(Q/P)/pと(y[1])′=(q−λw)y(y^{[1]})'=(q-\lambda w)yがほとんど至る所で成り立つから、ほとんど至る所で

y′′=(y[1])′p−QP y[1]p=(q−λw)yp−QPy′=−R+λWPy−QPy′y''=\frac{(y^{[1]})'}p-\frac QP\,\frac{y^{[1]}}p =\frac{(q-\lambda w)y}p-\frac QPy' =-\frac{R+\lambda W}Py-\frac QPy'

が成り立つ。

主張 (3)を示す。c′c'から定めた関数をpc′p_{c'}と書くと、∫c′xQ/P=∫c′cQ/P+∫cxQ/P\int_{c'}^xQ/P=\int_{c'}^cQ/P+\int_c^xQ/Pであるからpc′=exp⁡(∫c′cQ/P)pp_{c'}=\exp\bigl(\int_{c'}^cQ/P\bigr)pである。q,wq,wはppの定数倍について線形であるから、同じ定数倍になる。▨

例 1.9. 端点の型は、最高階係数PPの零点や∫Q/P\int Q/Pの発散では決まらない。

  1. 開区間(−1,1)(-1,1)上の方程式(1−x2)y′′−xy′+λy=0(1-x^2)y''-xy'+\lambda y=0ではP(x)=1−x2P(x)=1-x^2は±1\pm1で零になる。Q(x)=−xQ(x)=-x、R=0R=0、W=1W=1とし、補題 1.8をc=0c=0で適用すると、∫0xQ/P=12log⁡(1−x2)\int_0^xQ/P=\frac12\log(1-x^2)はx→±1x\to\pm1で−∞-\inftyに発散し、p(x)=(1−x2)1/2p(x)=(1-x^2)^{1/2}は±1\pm1で零に近づく。一方、q=0q=0であり、1/p=w=(1−x2)−1/21/p=w=(1-x^2)^{-1/2}は0≤x<10\le x<1で(1−x)−1/2(1-x)^{-1/2}以下である。0<ε<10<\varepsilon<1について∫01−ε(1−x)−1/2 dx=2−2ε1/2≤2\int_0^{1-\varepsilon}(1-x)^{-1/2}\,dx=2-2\varepsilon^{1/2}\le2であるから、§E10.11 定理 1.1により(1−x)−1/2(1-x)^{-1/2}は(0,1)(0,1)上で可積分である。x↦−xx\mapsto-xについての対称性により1/p,q,w1/p,q,wは(−1,0)(-1,0)上でも可積分であり、端点±1\pm1はともに正則端点である。
  2. 開区間(0,1)(0,1)上の方程式−y′′+x−2y=λy-y''+x^{-2}y=\lambda yではp=w=1p=w=1は[0,1][0,1]上の正の連続関数であるが、q(x)=x−2q(x)=x^{-2}は(0,1)(0,1)上で可積分でないから、端点00は特異端点である。λ=0\lambda=0のとき、r±=(1±5)/2r_\pm=(1\pm\sqrt5)/2はr2−r−1=0r^2-r-1=0の根であり、y±(x)=xr±y_\pm(x)=x^{r_\pm}はy±′′=r±(r±−1)xr±−2=x−2y±y_\pm''=r_\pm(r_\pm-1)x^{r_\pm-2}=x^{-2}y_\pmを満たす解である。r−<0r_-<0であるからy−y_-はx→0x\to0で有界でなく、00へ連続に延びない。したがって、この方程式の端点00では、λ=0\lambda=0について命題 1.3 (4)の結論が成り立たない。

注意 1.10. 有限閉区間[a,b][a,b]上の正則 Sturm–Liouville 式(定義 3.1)の係数p,q,wp,q,wを(a,b)(a,b)へ制限したものは定義 1.1の仮定を満たす。1/p,q,w1/p,q,wは[a,b][a,b]上の連続関数であるから、aaとbbはともに正則端点である。逆は成り立たない。実際、例 1.9 (1)の方程式では±1\pm1がともに正則端点であるが、p(x)=(1−x2)1/2p(x)=(1-x^2)^{1/2}は±1\pm1で零になり、[−1,1][-1,1]上の正則 Sturm–Liouville 式の係数ではない。正則 Sturm–Liouville 式であることは閉区間上の係数の滑らかさと正値性による条件であり、両端が正則端点であることは1/p,q,w1/p,q,wの可積分性だけによる条件である。

2 Liouville の標準形と古典的な方程式

定理 2.1.−∞≤a<b≤∞-\infty\le a<b\le\inftyとし、p,w:(a,b)→Rp,w:(a,b)\to\Rを正のC2C^2級関数、q:(a,b)→Rq:(a,b)\to\Rを局所可積分な可測関数とする。t:(a,b)→Rt:(a,b)\to\Rをw/p\sqrt{w/p}の原始関数の一つとし、ta=lim⁡x↓at(x)t_a=\lim_{x\downarrow a}t(x)、tb=lim⁡x↑bt(x)t_b=\lim_{x\uparrow b}t(x)を[−∞,∞][-\infty,\infty]の元とする。このとき、p,q,wp,q,wは(a,b)(a,b)上で定義 1.1の仮定を満たし、ttは(a,b)(a,b)から(ta,tb)(t_a,t_b)への狭義増加な全単射であり、その逆写像x:(ta,tb)→(a,b)x:(t_a,t_b)\to(a,b)はC2C^2級である。

m=(pw)1/4∘x,V=qw∘x+m′′mm=(pw)^{1/4}\circ x,\qquad V=\frac qw\circ x+\frac{m''}m

とおき、可測関数y:(a,b)→Cy:(a,b)\to\Cに対して(Uy)(τ)=m(τ) y(x(τ))(Uy)(\tau)=m(\tau)\,y(x(\tau))(τ∈(ta,tb)\tau\in(t_a,t_b))と定める。このとき、次が成り立つ。

  1. VVは(ta,tb)(t_a,t_b)上で局所可積分である。qqが連続ならばVVは連続である。
  2. UUはw dxw\,dxに関する二乗可積分関数の空間L2((a,b),w dx)L^2((a,b),w\,dx)からL2(ta,tb)L^2(t_a,t_b)へのユニタリ写像を定める。
  3. λ∈C\lambda\in\Cとし、y:(a,b)→Cy:(a,b)\to\Cを関数とする。yyがp,q,wp,q,wについての方程式ℓ[y]=λwy\ell[y]=\lambda wyの解であることと、u=Uyu=Uyがu∈C1((ta,tb);C)u\in C^1((t_a,t_b);\C)、u′∈ACloc((ta,tb);C)u'\in AC_{\mathrm{loc}}((t_a,t_b);\C)を満たし、−u′′+Vu=λu-u''+Vu=\lambda uがほとんど至る所で成り立つことは同値である。このときy[1]∘x=mu′−m′uy^{[1]}\circ x=mu'-m'uである。
  4. e∈{a,b}e\in\{a,b\}は有限であり、p,wp,wは(a,b)∪{e}(a,b)\cup\{e\}上の正のC2C^2級関数に延び、eeとccを端点とする開区間の上でqqが可積分であるようなc∈(a,b)c\in(a,b)が存在するとする。このときeeはp,q,wp,q,wについて正則端点であり、tet_eは有限であって、tet_eは1,V,11,V,1について正則端点である。mmとm′m'はtet_eへ連続に延び、m(te)>0m(t_e)>0である。主張 (3)のy,uy,uについてy(e)=u(te)/m(te)y(e)=u(t_e)/m(t_e)、y[1](e)=m(te)u′(te)−m′(te)u(te)y^{[1]}(e)=m(t_e)u'(t_e)-m'(t_e)u(t_e)である。したがって(α,β)∈R2∖{(0,0)}(\alpha,\beta)\in\R^2\setminus\{(0,0)\}について、αy(e)+βy[1](e)=0\alpha y(e)+\beta y^{[1]}(e)=0であることと (αm(te)−βm′(te))u(te)+βm(te)u′(te)=0\Bigl(\frac\alpha{m(t_e)}-\beta m'(t_e)\Bigr)u(t_e)+\beta m(t_e)u'(t_e)=0 であることは同値であり、後者の係数の組は(0,0)(0,0)でない実数の組で、β=0\beta=0のとき、かつそのときに限りu′(te)u'(t_e)の係数が零である。

証明.t′=w/pt'=\sqrt{w/p}は正のC2C^2級関数であるから、ttは狭義増加なC3C^3級関数である。τ∈(ta,tb)\tau\in(t_a,t_b)についてt(x1)<τ<t(x2)t(x_1)<\tau<t(x_2)を満たすx1<x2x_1<x_2を取り、§D1.12 系 1.2を[x1,x2][x_1,x_2]上のttに適用するとt(x)=τt(x)=\tauを満たすxxが存在するから、ttの像は(ta,tb)(t_a,t_b)である。§E4.7 定理 2.1をn=1n=1で各点に適用すると、逆写像xxはC1C^1級でありx′=1/(t′∘x)=p/w∘xx'=1/(t'\circ x)=\sqrt{p/w}\circ xである。右辺はC1C^1級であるからxxはC2C^2級である。p,wp,wは正の連続関数であるから1/p,w1/p,wは局所可積分であり、p,q,wp,q,wは定義 1.1の仮定を満たす。mmは正のC2C^2級関数(pw)1/4(pw)^{1/4}とxxの合成であるから正のC2C^2級関数であり、m2=pw∘xm^2=\sqrt{pw}\circ xである。

主張 2.1.1.J1,J2J_1,J_2を開区間、φ:J1→J2\varphi:J_1\to J_2をC1C^1級の全単射で逆写像もC1C^1級であり、すべての点でφ′>0\varphi'>0を満たすものとする。N⊆J2N\subseteq J_2が測度零ならばφ−1(N)\varphi^{-1}(N)は測度零である。F∈ACloc(J2;C)F\in AC_{\mathrm{loc}}(J_2;\C)ならばF∘φ∈ACloc(J1;C)F\circ\varphi\in AC_{\mathrm{loc}}(J_1;\C)であり、(F∘φ)′=(F′∘φ)φ′(F\circ\varphi)'=(F'\circ\varphi)\varphi'がほとんど至る所で成り立つ。

証明.§E10.44 定理 2.1 (2)をn=1n=1とC1C^1級微分同相写像φ−1:J2→J1\varphi^{-1}:J_2\to J_1に適用すると、φ−1(N)\varphi^{-1}(N)の測度はNN上の積分であるから零である。s0∈J1s_0\in J_1を固定し、s∈J1s\in J_1、s>s0s>s_0とする。φ\varphiの開区間(s0,s)(s_0,s)への制限は(φ(s0),φ(s))(\varphi(s_0),\varphi(s))へのC1C^1級微分同相写像であり、F′F'は(φ(s0),φ(s))(\varphi(s_0),\varphi(s))上で可積分であるから、§E10.44 定理 2.1 (4)と§E13.2 系 2.6により

F(φ(s))−F(φ(s0))=∫φ(s0)φ(s)F′(ξ) dξ=∫s0sF′(φ(σ))φ′(σ) dσF(\varphi(s))-F(\varphi(s_0))=\int_{\varphi(s_0)}^{\varphi(s)}F'(\xi)\,d\xi=\int_{s_0}^sF'(\varphi(\sigma))\varphi'(\sigma)\,d\sigma

が成り立ち、(F′∘φ)φ′(F'\circ\varphi)\varphi'は(s0,s)(s_0,s)上で可積分である。s<s0s<s_0についても向き付き積分で同じ等式が成り立つ。したがって§E13.2 系 2.6によりF∘φF\circ\varphiは局所絶対連続であり、その導関数はほとんど至る所で(F′∘φ)φ′(F'\circ\varphi)\varphi'に等しい。▨

主張 (1)を示す。[σ1,σ2]⊆(ta,tb)[\sigma_1,\sigma_2]\subseteq(t_a,t_b)とする。§E10.44 定理 2.1 (3)をxxの(σ1,σ2)(\sigma_1,\sigma_2)への制限とg=∣q∣/pwg=\lvert q\rvert/\sqrt{pw}に適用すると、(g∘x)x′=∣q/w∣∘x(g\circ x)x'=\lvert q/w\rvert\circ xであるから

∫σ1σ2∣qw∣∘x dτ=∫x(σ1)x(σ2)∣q∣pw dξ\int_{\sigma_1}^{\sigma_2}\Bigl\lvert\frac qw\Bigr\rvert\circ x\,d\tau=\int_{x(\sigma_1)}^{x(\sigma_2)}\frac{\lvert q\rvert}{\sqrt{pw}}\,d\xi

が成り立つ。1/pw1/\sqrt{pw}は[x(σ1),x(σ2)][x(\sigma_1),x(\sigma_2)]上で有界であり、qqはその上で可積分であるから、右辺は有限である。m′′/mm''/mは連続であるからVVは局所可積分である。qqが連続ならば(q/w)∘x(q/w)\circ xは連続であり、VVは連続である。

主張 (2)を示す。yyを可測関数とする。§E10.44 定理 2.1 (3)をC1C^1級微分同相写像xxとyyの実部・虚部の正の部分と負の部分に適用すると、y∘xy\circ xは可測であり、UyUyは可測である。同じ定理をg=w∣y∣2g=w\lvert y\rvert^2に適用すると、x′=p/w∘xx'=\sqrt{p/w}\circ xであるから

∫tatb∣Uy∣2 dτ=∫tatb(pw ∣y∣2)∘x dτ=∫tatb(w∣y∣2)∘x⋅x′ dτ=∫abw∣y∣2 dx\int_{t_a}^{t_b}\lvert Uy\rvert^2\,d\tau=\int_{t_a}^{t_b}\bigl(\sqrt{pw}\,\lvert y\rvert^2\bigr)\circ x\,d\tau=\int_{t_a}^{t_b}\bigl(w\lvert y\rvert^2\bigr)\circ x\cdot x'\,d\tau=\int_a^bw\lvert y\rvert^2\,dx

が成り立つ。主張 2.1.1により、ほとんど至る所で一致する関数はUUによってほとんど至る所で一致する関数に写る。したがってUUはL2((a,b),w dx)L^2((a,b),w\,dx)からL2(ta,tb)L^2(t_a,t_b)への線形な等長写像を定める。u∈L2(ta,tb)u\in L^2(t_a,t_b)に対し、y=(u/m)∘ty=(u/m)\circ tは、§E10.44 定理 2.1 (3)をttとu/mu/mの実部・虚部の正の部分と負の部分に適用すると可測であり、Uy=uUy=uを満たす。したがってUUは全射であり、ユニタリ写像である。

主張 (3)を示す。yyが解であるとする。命題 1.3 (5)によりyyはC1C^1級でありy[1]=py′y^{[1]}=py'である。z=y∘xz=y\circ xはC1C^1級でz′=(y′∘x) x′z'=(y'\circ x)\,x'であり、u=mzu=mzはC1C^1級でu′=m′z+mz′u'=m'z+mz'である。したがって

mu′−m′u=m2z′=(pw y′p/w)∘x=(py′)∘x=y[1]∘xmu'-m'u=m^2z'=\bigl(\sqrt{pw}\,y'\sqrt{p/w}\bigr)\circ x=(py')\circ x=y^{[1]}\circ x

である。G=y[1]∘xG=y^{[1]}\circ xとおくと、主張 2.1.1をφ=x\varphi=xに適用して、GGは局所絶対連続であり、ほとんど至る所で

G′=((y[1])′∘x)x′=((q−λw)yp/w)∘x=(qw∘x−λ)(pw∘x)um=(qw∘x−λ)muG'=\bigl((y^{[1]})'\circ x\bigr)x'=\bigl((q-\lambda w)y\sqrt{p/w}\bigr)\circ x=\Bigl(\frac qw\circ x-\lambda\Bigr)\bigl(\sqrt{pw}\circ x\bigr)\frac um=\Bigl(\frac qw\circ x-\lambda\Bigr)mu

が成り立つ。m,m′m,m'はC1C^1級であるから、§E13.2 定理 3.1 (2)によりu′=(G+m′u)/mu'=(G+m'u)/mは局所絶対連続である。G=mu′−m′uG=mu'-m'uをほとんど至る所で微分するとG′=mu′′−m′′uG'=mu''-m''uであるから、mu′′−m′′u=((q/w)∘x−λ)mumu''-m''u=((q/w)\circ x-\lambda)muがほとんど至る所で成り立ち、m>0m>0で割ると−u′′+Vu=λu-u''+Vu=\lambda uがほとんど至る所で成り立つ。

逆に、u=Uyu=Uyがu∈C1u\in C^1、u′∈AClocu'\in AC_{\mathrm{loc}}、−u′′+Vu=λu-u''+Vu=\lambda u(ほとんど至る所)を満たすとする。H=mu′−m′uH=mu'-m'uは§E13.2 定理 3.1 (2)により局所絶対連続であり、ほとんど至る所で

H′=mu′′−m′′u=m(V−λ)u−m′′u=(qw∘x−λ)muH'=mu''-m''u=m(V-\lambda)u-m''u=\Bigl(\frac qw\circ x-\lambda\Bigr)mu

が成り立つ。y=(u/m)∘ty=(u/m)\circ tはC1C^1級であり、y′=((u′m−um′)/m2)∘t⋅w/py'=\bigl((u'm-um')/m^2\bigr)\circ t\cdot\sqrt{w/p}であるから、m2∘t=pwm^2\circ t=\sqrt{pw}により

py′=pw (Hm2)∘t=H∘tpy'=\sqrt{pw}\,\Bigl(\frac{H}{m^2}\Bigr)\circ t=H\circ t

がすべての点で成り立つ。主張 2.1.1をφ=t\varphi=tに適用すると、H∘tH\circ tは局所絶対連続であるからyyの流束であり、H′H'の等式が成り立たない点の集合のttによる逆像は測度零であるから、ほとんど至る所で

(y[1])′=(H′∘t) t′=(qw−λ)(m∘t)(u∘t)wp=(qw−λ)pw y wp=(q−λw)y(y^{[1]})'=(H'\circ t)\,t'=\Bigl(\frac qw-\lambda\Bigr)(m\circ t)(u\circ t)\sqrt{\frac wp}=\Bigl(\frac qw-\lambda\Bigr)\sqrt{pw}\,y\,\sqrt{\frac wp}=(q-\lambda w)y

が成り立つ。yyはC1C^1級であるから局所絶対連続であり、ℓ[y]=λwy\ell[y]=\lambda wyの解である。y[1]∘x=H=mu′−m′uy^{[1]}\circ x=H=mu'-m'uである。

主張 (4)を示す。e=ae=aの場合を述べる。e=be=bの場合は、以下の[a,c][a,c]、(a,c)(a,c)、(ta,t(c))(t_a,t(c))、x↓ax\downarrow a、τ↓ta\tau\downarrow t_aを[c,b][c,b]、(c,b)(c,b)、(t(c),tb)(t(c),t_b)、x↑bx\uparrow b、τ↑tb\tau\uparrow t_bに置き換えて同じ議論が成り立つ。p,wp,wは[a,c][a,c]上の正の連続関数に延びるから、1/p1/pとwwは(a,c)(a,c)上で可積分であり、仮定によりqqも可積分であるから、aaは正則端点である。w/p\sqrt{w/p}は[a,c][a,c]上で有界であるから、t(c)−t(x)=∫xcw/p dξt(c)-t(x)=\int_x^c\sqrt{w/p}\,d\xiはx∈(a,c)x\in(a,c)について有界であり、tat_aは有限である。k0=(pw)1/4k_0=(pw)^{1/4}、k1=k0′p/wk_1=k_0'\sqrt{p/w}、k2=k1′p/wk_2=k_1'\sqrt{p/w}は[a,c][a,c]上の連続関数に延び、x′=p/w∘xx'=\sqrt{p/w}\circ xであるからm=k0∘xm=k_0\circ x、m′=k1∘xm'=k_1\circ x、m′′=k2∘xm''=k_2\circ xである。τ↓ta\tau\downarrow t_aとx(τ)↓ax(\tau)\downarrow aは同値であるから、mmとm′m'はtat_aへ連続に延び、m(ta)=k0(a)>0m(t_a)=k_0(a)>0であり、m′′/mm''/mは(ta,t(c))(t_a,t(c))上で有界である。§E10.44 定理 2.1 (3)をxxの(ta,t(c))(t_a,t(c))への制限とg=∣q∣/pwg=\lvert q\rvert/\sqrt{pw}に適用すると、主張 (1)の証明と同じ等式が(ta,t(c))(t_a,t(c))と(a,c)(a,c)の上で成り立ち、∣q∣/pw\lvert q\rvert/\sqrt{pw}は(a,c)(a,c)上で可積分であるから、(q/w)∘x(q/w)\circ xは(ta,t(c))(t_a,t(c))上で可積分である。したがってtat_aは1,V,11,V,1について正則端点である。命題 1.3 (4)と命題 1.3 (5)により、uuとu′u'はtat_aへ連続に延び、yyとy[1]y^{[1]}はaaへ連続に延びる。x↓ax\downarrow aとt(x)↓tat(x)\downarrow t_aは同値であるから、y=(u/m)∘ty=(u/m)\circ tとy[1]=(mu′−m′u)∘ty^{[1]}=(mu'-m'u)\circ tの極限をとってy(a)=u(ta)/m(ta)y(a)=u(t_a)/m(t_a)、y[1](a)=m(ta)u′(ta)−m′(ta)u(ta)y^{[1]}(a)=m(t_a)u'(t_a)-m'(t_a)u(t_a)を得る。これらを代入すると

αy(a)+βy[1](a)=(αm(ta)−βm′(ta))u(ta)+βm(ta)u′(ta)\alpha y(a)+\beta y^{[1]}(a)=\Bigl(\frac\alpha{m(t_a)}-\beta m'(t_a)\Bigr)u(t_a)+\beta m(t_a)u'(t_a)

である。m(ta)>0m(t_a)>0であるから、β≠0\beta\ne0ならばu′(ta)u'(t_a)の係数は零でなく、β=0\beta=0ならばα≠0\alpha\ne0でありu(ta)u(t_a)の係数α/m(ta)\alpha/m(t_a)は零でない。▨

命題 2.2.λ∈C\lambda\in\Cとする。このとき、次の各方程式について、y∈C1y\in C^1、y′∈AClocy'\in AC_{\mathrm{loc}}であって方程式がほとんど至る所で成り立つ関数yyの全体は、示した開区間上で示したp,q,wp,q,wについてのℓ[y]=λwy\ell[y]=\lambda wyの解の全体に一致する。p,q,wp,q,wは補題 1.8 (3)の正の定数倍を除いて定まり、定義 1.1の仮定を満たす。各端点の型は次のとおりである。

  1. α,β>−1\alpha,\beta>-1とする。Jacobi の方程式(1−x2)y′′+(β−α−(α+β+2)x)y′+λy=0(1-x^2)y''+\bigl(\beta-\alpha-(\alpha+\beta+2)x\bigr)y'+\lambda y=0は(−1,1)(-1,1)上でp(x)=(1−x)α+1(1+x)β+1p(x)=(1-x)^{\alpha+1}(1+x)^{\beta+1}、q=0q=0、w(x)=(1−x)α(1+x)βw(x)=(1-x)^\alpha(1+x)^\betaをもつ。端点11はα<0\alpha<0のとき正則端点、α≥0\alpha\ge0のとき特異端点である。端点−1-1はβ<0\beta<0のとき正則端点、β≥0\beta\ge0のとき特異端点である。
  2. Legendre の方程式(1−x2)y′′−2xy′+λy=0(1-x^2)y''-2xy'+\lambda y=0は(−1,1)(-1,1)上でp(x)=1−x2p(x)=1-x^2、q=0q=0、w=1w=1をもつ。端点±1\pm1はともに特異端点である。
  3. Chebyshev の方程式(1−x2)y′′−xy′+λy=0(1-x^2)y''-xy'+\lambda y=0は(−1,1)(-1,1)上でp(x)=(1−x2)1/2p(x)=(1-x^2)^{1/2}、q=0q=0、w(x)=(1−x2)−1/2w(x)=(1-x^2)^{-1/2}をもつ。端点±1\pm1はともに正則端点である。
  4. α>−1\alpha>-1とする。Laguerre の方程式xy′′+(α+1−x)y′+λy=0xy''+(\alpha+1-x)y'+\lambda y=0は(0,∞)(0,\infty)上でp(x)=xα+1e−xp(x)=x^{\alpha+1}e^{-x}、q=0q=0、w(x)=xαe−xw(x)=x^\alpha e^{-x}をもつ。端点00はα<0\alpha<0のとき正則端点、α≥0\alpha\ge0のとき特異端点であり、端点∞\inftyは特異端点である。
  5. Hermite の方程式y′′−2xy′+λy=0y''-2xy'+\lambda y=0はR\R上でp(x)=w(x)=e−x2p(x)=w(x)=e^{-x^2}、q=0q=0をもち、端点±∞\pm\inftyはともに特異端点である。u(x)=e−x2/2y(x)u(x)=e^{-x^2/2}y(x)とおくと、yyが解であることと、u∈C1(R;C)u\in C^1(\R;\C)、u′∈ACloc(R;C)u'\in AC_{\mathrm{loc}}(\R;\C)が−u′′+(x2−1)u=λu-u''+(x^2-1)u=\lambda u、すなわち−u′′+x2u=(λ+1)u-u''+x^2u=(\lambda+1)uをほとんど至る所で満たすことは同値であり、y↦uy\mapsto uはL2(R,e−x2dx)L^2(\R,e^{-x^2}dx)からL2(R)L^2(\R)へのユニタリ写像を定める。
  6. ν≥0\nu\ge0とする。方程式x2y′′+xy′+(λx2−ν2)y=0x^2y''+xy'+(\lambda x^2-\nu^2)y=0は(0,∞)(0,\infty)上でp(x)=xp(x)=x、q(x)=ν2/xq(x)=\nu^2/x、w(x)=xw(x)=xをもつ。端点00はすべてのν≥0\nu\ge0について特異端点であり、端点∞\inftyは特異端点である。同じp,q,wp,q,wを(0,1)(0,1)へ制限すると、端点11は正則端点である。u(x)=x1/2y(x)u(x)=x^{1/2}y(x)とおくと、yyが解であることと、u∈C1((0,∞);C)u\in C^1((0,\infty);\C)、u′∈ACloc((0,∞);C)u'\in AC_{\mathrm{loc}}((0,\infty);\C)が−u′′+(ν2−14)x−2u=λu-u''+(\nu^2-\tfrac14)x^{-2}u=\lambda uをほとんど至る所で満たすことは同値であり、y↦uy\mapsto uはL2((0,∞),x dx)L^2((0,\infty),x\,dx)からL2(0,∞)L^2(0,\infty)へのユニタリ写像を定める。
  7. T>0T>0とし、V:R→RV:\R\to\Rを周期TTの局所可積分関数とする。Hill の方程式−y′′+Vy=λy-y''+Vy=\lambda yはR\R上でp=w=1p=w=1、q=Vq=Vをもち、端点±∞\pm\inftyはともに特異端点である。同じp,q,wp,q,wを(0,T)(0,T)へ制限すると、端点0,T0,Tはともに正則端点である。h∈Rh\in\RについてV(x)=2hcos⁡2xV(x)=2h\cos2x、T=πT=\piとした場合が Mathieu の方程式である。

証明.

主張 2.2.1.γ∈R\gamma\in\R、δ>0\delta>0とする。s↦sγs\mapsto s^\gammaが(0,δ)(0,\delta)上で可積分であるための必要十分条件はγ>−1\gamma>-1である。

証明.0<ε<δ0<\varepsilon<\deltaについて、γ≠−1\gamma\ne-1ならば∫εδsγ ds=(δγ+1−εγ+1)/(γ+1)\int_\varepsilon^\delta s^\gamma\,ds=(\delta^{\gamma+1}-\varepsilon^{\gamma+1})/(\gamma+1)、γ=−1\gamma=-1ならば∫εδs−1 ds=log⁡δ−log⁡ε\int_\varepsilon^\delta s^{-1}\,ds=\log\delta-\log\varepsilonである。ε↓0\varepsilon\downarrow0のとき、この値はγ>−1\gamma>-1ならばδγ+1/(γ+1)\delta^{\gamma+1}/(\gamma+1)に収束し、γ≤−1\gamma\le-1ならば+∞+\inftyに発散する。§E10.11 定理 1.1により∫0δsγ ds\int_0^\delta s^\gamma\,dsはこの極限に等しい。▨

各方程式をPy′′+Qy′+(R+λW)y=0Py''+Qy'+(R+\lambda W)y=0の形に書くと、PPは区間上で零にならず、Q/P,R/P,W/PQ/P,R/P,W/Pは局所可積分であり、W/P>0W/P>0である。したがって補題 1.8 (1)によりp,q,wp,q,wは定義 1.1の仮定を満たし、解の全体の一致は補題 1.8 (2)による。κ>0\kappa>0についてκp,κq,κw\kappa p,\kappa q,\kappa wはp,q,wp,q,wと同じ解の全体と同じ端点の型をもつから、p,q,wp,q,wを正の定数倍で取り替えてよい。±∞\pm\inftyである端点は定義 1.2により特異端点である。

主張 (1)を示す。P=1−x2P=1-x^2、Q=β−α−(α+β+2)xQ=\beta-\alpha-(\alpha+\beta+2)x、R=0R=0、W=1W=1である。

QP=−(α+1)(1+x)+(β+1)(1−x)(1−x)(1+x)=−α+11−x+β+11+x\frac QP=\frac{-(\alpha+1)(1+x)+(\beta+1)(1-x)}{(1-x)(1+x)}=-\frac{\alpha+1}{1-x}+\frac{\beta+1}{1+x}

であるから、c=0c=0について∫0xQ/P=(α+1)log⁡(1−x)+(β+1)log⁡(1+x)\int_0^xQ/P=(\alpha+1)\log(1-x)+(\beta+1)\log(1+x)であり、p,q,wp,q,wは示したとおりである。0≤x<10\le x<1では(1+x)β+1(1+x)^{\beta+1}と(1+x)β(1+x)^\betaは正の定数で上下から押さえられるから、1/p1/pが(0,1)(0,1)上で可積分であることは(1−x)−α−1(1-x)^{-\alpha-1}が(0,1)(0,1)上で可積分であることと同値であり、s=1−xs=1-xと置き換えて主張 2.2.1をγ=−α−1\gamma=-\alpha-1、δ=1\delta=1で適用すると、−α−1>−1-\alpha-1>-1、すなわちα<0\alpha<0と同値である。同じ理由で、wwが(0,1)(0,1)上で可積分であることはα>−1\alpha>-1と同値であり、仮定により成り立つ。q=0q=0であるから、端点11はα<0\alpha<0のとき正則端点、α≥0\alpha\ge0のとき特異端点である。端点−1-1についてはx↦−xx\mapsto-xでα\alphaとβ\betaの役割が入れ替わり、同じ結論が成り立つ。

主張 (2)と主張 (3)を示す。Legendre の方程式は主張 (1)のα=β=0\alpha=\beta=0の場合であり、Chebyshev の方程式はα=β=−12\alpha=\beta=-\tfrac12の場合である。

主張 (4)を示す。P=xP=x、Q=α+1−xQ=\alpha+1-x、R=0R=0、W=1W=1であり、c=1c=1について∫1xQ/P=(α+1)log⁡x−x+1\int_1^xQ/P=(\alpha+1)\log x-x+1であるから、p=e xα+1e−xp=e\,x^{\alpha+1}e^{-x}、q=0q=0、w=e xαe−xw=e\,x^\alpha e^{-x}であり、正の定数eeを除いて示したとおりである。0<x<10<x<1ではe−1≤e−x≤1e^{-1}\le e^{-x}\le1であるから、主張 2.2.1により、1/p1/pが(0,1)(0,1)上で可積分であることはα<0\alpha<0と同値であり、wwはα>−1\alpha>-1により(0,1)(0,1)上で可積分である。

主張 (5)を示す。P=1P=1、Q=−2xQ=-2x、R=0R=0、W=1W=1であり、c=0c=0について∫0xQ/P=−x2\int_0^xQ/P=-x^2であるから、p=w=e−x2p=w=e^{-x^2}、q=0q=0である。定理 2.1を(a,b)=R(a,b)=\Rに適用する。w/p=1\sqrt{w/p}=1であるからt(x)=xt(x)=xとし、m(x)=(e−2x2)1/4=e−x2/2m(x)=(e^{-2x^2})^{1/4}=e^{-x^2/2}である。m′=−xmm'=-xm、m′′=(x2−1)mm''=(x^2-1)mであるからV=x2−1V=x^2-1であり、主張は定理 2.1 (2)と定理 2.1 (3)から従う。

主張 (6)を示す。P=x2P=x^2、Q=xQ=x、R=−ν2R=-\nu^2、W=x2W=x^2であり、c=1c=1について∫1xQ/P=log⁡x\int_1^xQ/P=\log xであるから、p=xp=x、q=−pR/P=ν2/xq=-pR/P=\nu^2/x、w=pW/P=xw=pW/P=xである。1/p=1/x1/p=1/xは主張 2.2.1により(0,1)(0,1)上で可積分でないから、端点00は特異端点である。1/p,q,w1/p,q,wは[12,1][\tfrac12,1]上で連続であるから、(0,1)(0,1)の端点11は正則端点である。定理 2.1を(a,b)=(0,∞)(a,b)=(0,\infty)に適用する。w/p=1\sqrt{w/p}=1であるからt(x)=xt(x)=xとし、m(x)=(x2)1/4=x1/2m(x)=(x^2)^{1/4}=x^{1/2}である。m′′(x)=−14x−3/2m''(x)=-\tfrac14x^{-3/2}であるからV(x)=ν2x−2−14x−2V(x)=\nu^2x^{-2}-\tfrac14x^{-2}であり、主張は定理 2.1 (2)と定理 2.1 (3)から従う。

主張 (7)を示す。P=−1P=-1、Q=0Q=0、R=VR=V、W=−1W=-1とすると方程式はPy′′+Qy′+(R+λW)y=0Py''+Qy'+(R+\lambda W)y=0の形であり、c=0c=0についてp=1p=1、q=−pR/P=Vq=-pR/P=V、w=pW/P=1w=pW/P=1である。VVは[0,T][0,T]上で可積分であるから、(0,T)(0,T)の端点0,T0,Tは正則端点である。▨

注意 2.3. 特異端点eeは、ℓ[y]=λwy\ell[y]=\lambda wyのすべての解がeeの近くでw dxw\,dxに関して二乗可積分であるようなλ∈C\lambda\in\Cが存在するとき極限円型、そうでないとき極限点型と呼ばれる。この区別はλ\lambdaの取り方によらない。命題 2.2のうち、Legendre、Laguerre、Hermite、Bessel、Hill の方程式の特異端点の型は次のとおりである。Legendre の方程式の±1\pm1は極限円型である。Laguerre の方程式の00は0≤α<10\le\alpha<1のとき極限円型、α≥1\alpha\ge1のとき極限点型であり、∞\inftyは極限点型である。Hermite の方程式の±∞\pm\inftyは極限点型であり、Liouville の標準形−u′′+(x2−1)u=λu-u''+(x^2-1)u=\lambda uの±∞\pm\inftyも極限点型である。Bessel の方程式の00は0≤ν<10\le\nu<1のとき極限円型、ν≥1\nu\ge1のとき極限点型であり、∞\inftyは極限点型である。Hill の方程式の±∞\pm\inftyは極限点型であり、V=0V=0の場合として、R\R上および(0,∞)(0,\infty)上の−u′′=λu-u''=\lambda uの無限遠の端点も極限点型である。これらの事実は本記事のどの証明にも用いない。

注意 2.4. 次の型の問題は本記事の定理の仮定を満たさない。

  1. 境界条件がλ\lambdaを含む問題。たとえば[0,1][0,1]上の−y′′=λy-y''=\lambda yにy(0)=0y(0)=0、y′(1)=λy(1)y'(1)=\lambda y(1)を課すと、異なるλ,μ\lambda,\muに属する解は端点11で異なる係数の組(λ,−1)(\lambda,-1)、(μ,−1)(\mu,-1)の条件を満たし、定理 1.4 (2)の仮定が成り立たない。
  2. 正の測度の集合の上でw<0w<0となる重みをもつ問題。∫w∣y∣2 dx\int w\lvert y\rvert^2\,dxは正定値でなく、重み付き内積を定めない。
  3. 実数値でないqqまたは実数でない境界係数をもつ問題。定理 1.4 (1)の証明はquv‾=uqv‾qu\overline v=u\overline{qv}を用い、定理 1.4 (2)の証明は境界係数が実数であることを用いる。問題 7.3は実数でない境界係数をもつ問題の例である。
  4. 四階以上の微分式の固有値問題。

3 正則 Sturm–Liouville 問題

定義 3.1.a<ba<bとする。実数値関数p,q,wp,q,wが

p∈C1([a,b]),q,w∈C([a,b]),p(x)>0,w(x)>0(a≤x≤b)p\in C^1([a,b]),\qquad q,w\in C([a,b]),\qquad p(x)>0,\quad w(x)>0 \quad(a\le x\le b)

を満たすとする。微分式

L[y]=−(p(x)y′(x))′+q(x)y(x)L[y]=-\bigl(p(x)y'(x)\bigr)'+q(x)y(x)(3.1.1)\text{(3.1.1)}

を正則 Sturm–Liouville 式 (regular Sturm–Liouville expression) という。

端点でppまたはwwが零になる係数は正則 Sturm–Liouville 式の係数ではないから、正則 Sturm–Liouville 式の係数を仮定する本節以降の固有値問題の定理をそのまま適用することはできない。このような係数でも、p,q,wp,q,wの(a,b)(a,b)への制限が定義 1.1の仮定を満たし、両端が正則端点であれば、命題 1.7を適用することができる。命題 2.2 (2)の Legendre の方程式の係数p(x)=1−x2p(x)=1-x^2、q=0q=0、w=1w=1は端点±1\pm1でppが零になる係数であり、端点±1\pm1はともに特異端点であるから、定理 1.4をc=−1c=-1、d=1d=1で適用することもできない。

定義 3.2.a<ba<bとし、p,q,wp,q,wを定義 3.1の係数とする。実数αa,βa,αb,βb\alpha_a,\beta_a,\alpha_b,\beta_bが

(αa,βa)≠(0,0),(αb,βb)≠(0,0)(\alpha_a,\beta_a)\ne(0,0),\qquad (\alpha_b,\beta_b)\ne(0,0)

を満たすとする。

αay(a)+βap(a)y′(a)=0,αby(b)+βbp(b)y′(b)=0\begin{aligned} \alpha_a y(a)+\beta_a p(a)y'(a)&=0,\\ \alpha_b y(b)+\beta_b p(b)y'(b)&=0 \end{aligned}(3.2.1)\text{(3.2.1)}

を実係数の斉次分離型境界条件 (real homogeneous separated boundary conditions) という。

式 (3.1.1)、式 (3.2.1)の下で

L[y]=λw(x)yL[y]=\lambda w(x)y(3.2.2)\text{(3.2.2)}

を正則 Sturm–Liouville 固有値問題 (regular Sturm–Liouville eigenvalue problem) という。p,q,wp,q,wの(a,b)(a,b)への制限は定義 1.1の仮定を満たし、1/p,q,w1/p,q,wは[a,b][a,b]上で連続であるから、定義 1.2によりa,ba,bはともに正則端点である。この問題の固有値、固有関数、固有空間は、p,q,wp,q,wの(a,b)(a,b)への制限と境界係数(αa,βa),(αb,βb)(\alpha_a,\beta_a),(\alpha_b,\beta_b)についての定義 1.6の固有値、固有関数、固有空間とし、固有関数は[a,b][a,b]上の連続な延長と同一視する。

βa=βb=0\beta_a=\beta_b=0の場合を Dirichlet 条件、αa=αb=0\alpha_a=\alpha_b=0の場合を Neumann 条件という。端点eeでαeβe≠0\alpha_e\beta_e\ne0である条件を、その端点での Robin 条件という。

補題 3.3.a<ba<bとし、p,q,wp,q,wを定義 3.1の係数、(αa,βa),(αb,βb)(\alpha_a,\beta_a),(\alpha_b,\beta_b)を定義 3.2の境界係数とし、λ∈C\lambda\in\Cとする。このとき、次が成り立つ。

  1. 連続関数y:[a,b]→Cy:[a,b]\to\Cについて、yyの(a,b)(a,b)への制限がp,q,wp,q,wの(a,b)(a,b)への制限についてのℓ[y]=λwy\ell[y]=\lambda wyの解であることと、y∈C2([a,b];C)y\in C^2([a,b];\C)であり[a,b][a,b]の各点でL[y]=λwyL[y]=\lambda wyが成り立つことは同値である。このとき、補題 1.5 (2)の延長について、[a,b][a,b]の各点でy[1]=py′y^{[1]}=py'である。
  2. λ\lambdaが定義 3.2の固有値であることと、式 (3.2.1)を満たすL[y]=λwyL[y]=\lambda wyの零でない解y∈C2([a,b];C)y\in C^2([a,b];\C)が存在することは同値であり、λ\lambdaに属する固有関数はこのようなyyに限る。

証明.主張 (1)を示す。y∈C2([a,b];C)y\in C^2([a,b];\C)が[a,b][a,b]上でL[y]=λwyL[y]=\lambda wyを満たすとする。py′∈C1([a,b];C)py'\in C^1([a,b];\C)であるから、§E13.2 系 2.7 (2)によりyyとpy′py'は[a,b][a,b]上で局所絶対連続である。したがってpy′py'の(a,b)(a,b)への制限はyyの流束であり、−(py′)′+qy=L[y]=λwy-(py')'+qy=L[y]=\lambda wyが(a,b)(a,b)の各点で成り立つから、yyの制限は解である。流束の[a,b][a,b]への連続な延長はpy′py'である。

逆に、yyの制限が解であるとする。定義 3.2によりa,ba,bはともに正則端点であるから、補題 1.5 (2)により、流束は[a,b][a,b]上の絶対連続関数Y1Y_1に延び、yyの制限の連続な延長はyy自身であるから、yyは[a,b][a,b]上で絶対連続である。ほとんど至る所でY1′=(q−λw)yY_1'=(q-\lambda w)yとy′=Y1/py'=Y_1/pが成り立ち、右辺はともに[a,b][a,b]上で連続である。§E13.2 系 2.7 (1)をI=[a,b]I=[a,b]に適用すると、Y1Y_1とyyは[a,b][a,b]上で連続微分可能であり、[a,b][a,b]の各点でY1′=(q−λw)yY_1'=(q-\lambda w)yかつy′=Y1/py'=Y_1/pが成り立つ。p∈C1([a,b])p\in C^1([a,b])は正であるから1/p∈C1([a,b])1/p\in C^1([a,b])であり、y′=Y1⋅(1/p)y'=Y_1\cdot(1/p)は連続微分可能である。したがってy∈C2([a,b];C)y\in C^2([a,b];\C)であり、py′=Y1py'=Y_1とL[y]=−Y1′+qy=λwyL[y]=-Y_1'+qy=\lambda wyが[a,b][a,b]の各点で成り立つ。

主張 (2)を示す。補題 1.5 (2)により、p,q,wp,q,wの制限についてのℓ[y]=λwy\ell[y]=\lambda wyの解は[a,b][a,b]上の連続関数にただ一つの方法で延び、主張 (1)により、その延長はL[y]=λwyL[y]=\lambda wyを満たすC2([a,b];C)C^2([a,b];\C)の元であって、e∈{a,b}e\in\{a,b\}でy[1](e)=p(e)y′(e)y^{[1]}(e)=p(e)y'(e)を満たす。逆に、L[y]=λwyL[y]=\lambda wyを満たすy∈C2([a,b];C)y\in C^2([a,b];\C)の制限は解である。この対応で零でない関数は零でない関数に対応し、定義 1.6の境界条件αey(e)+βey[1](e)=0\alpha_ey(e)+\beta_ey^{[1]}(e)=0は式 (3.2.1)の条件αey(e)+βep(e)y′(e)=0\alpha_ey(e)+\beta_ep(e)y'(e)=0に一致する。▨

4 Lagrange の恒等式

定理 4.1.a<ba<bとし、p:[a,b]→Rp:[a,b]\to\Rを連続関数でp(x)>0p(x)>0をすべてのx∈[a,b]x\in[a,b]で満たすもの、q:[a,b]→Rq:[a,b]\to\Rを可積分関数とする。u,v∈AC([a,b];C)u,v\in AC([a,b];\C)とし、あるU,V∈AC([a,b];C)U,V\in AC([a,b];\C)が存在して、ほとんど至る所で

u′=Up,v′=Vpu'=\frac{U}{p},\qquad v'=\frac{V}{p}

が成り立つと仮定する。このときu,v∈C1([a,b])u,v\in C^1([a,b])であり、すべてのx∈[a,b]x\in[a,b]に対してU(x)=p(x)u′(x)U(x)=p(x)u'(x)、V(x)=p(x)v′(x)V(x)=p(x)v'(x)が成り立つ。ほとんど至る所で

L[u]=−U′+qu,L[v]=−V′+qvL[u]=-U'+qu,\qquad L[v]=-V'+qv

と定めると、L[u],L[v]∈L1([a,b];C)L[u],L[v]\in L^1([a,b];\C)であり、次が成り立つ。

  1. 積分恒等式 ∫ab(L[u] v‾−u L[v]‾)dx=[uV‾−Uv‾]ab=[p(uv′‾−u′v‾)]ab\int_a^b\left(L[u]\,\overline v-u\,\overline{L[v]}\right)dx =\left[u\overline V-U\overline v\right]_a^b =\left[p\bigl(u\overline{v'}-u'\overline v\bigr)\right]_a^b が成り立つ。
  2. 実数αa,βa,αb,βb\alpha_a,\beta_a,\alpha_b,\beta_bが(αa,βa)≠(0,0)(\alpha_a,\beta_a)\ne(0,0)、(αb,βb)≠(0,0)(\alpha_b,\beta_b)\ne(0,0)を満たすとする。uuとvvがともに αay(a)+βap(a)y′(a)=0,αby(b)+βbp(b)y′(b)=0\begin{aligned} \alpha_a y(a)+\beta_a p(a)y'(a)&=0,\\ \alpha_b y(b)+\beta_b p(b)y'(b)&=0 \end{aligned} を満たすならば、 [p(uv′‾−u′v‾)]ab=0\left[p\bigl(u\overline{v'}-u'\overline v\bigr)\right]_a^b=0 であり、主張 (1)の積分は零である。

証明.ppは正の連続関数であり、UUとVVは連続であるから、U/pU/pとV/pV/pは[a,b][a,b]上で連続である。§E13.2 系 2.7 (1)をI=[a,b]I=[a,b]に適用すると、u,v∈C1([a,b])u,v\in C^1([a,b])であり、すべてのx∈[a,b]x\in[a,b]でu′=U/pu'=U/p、v′=V/pv'=V/p、すなわちU=pu′U=pu'、V=pv′V=pv'が成り立つ。また、U′,V′∈L1([a,b];C)U',V'\in L^1([a,b];\C)であり、qqは可積分、u,vu,vは有界であるから、L[u]L[u]とL[v]L[v]は可積分である。

1/p1/pと11は[a,b][a,b]上で連続であり、qqは[a,b][a,b]上で可積分であるから、p,q,1p,q,1の(a,b)(a,b)への制限は定義 1.1の仮定を満たし、1/p,q,11/p,q,1は(a,b)(a,b)上で可積分である。定義 1.2によりaaとbbはともに正則端点である。U,VU,Vの(a,b)(a,b)への制限は局所絶対連続であり、ほとんど至る所でpu′=Upu'=U、pv′=Vpv'=Vが成り立つから、それぞれu,vu,vの流束である。したがってほとんど至る所でℓ[u]=−U′+qu=L[u]\ell[u]=-U'+qu=L[u]、ℓ[v]=L[v]\ell[v]=L[v]が成り立つ。u,U,v,Vu,U,v,Vは[a,b][a,b]上で絶対連続であるから、定理 1.4をc=ac=a、d=bd=bに適用することができ、そこでのHHはuV‾−Uv‾u\overline V-U\overline vである。定理 1.4 (1)により主張 (1)の第一の等号が成り立つ。ppは実数値であるからV‾=pv′‾\overline V=p\overline{v'}であり、U=pu′U=pu'と合わせて第二の等号が成り立つ。

e∈{a,b}e\in\{a,b\}とし、その端点の境界係数を(αe,βe)(\alpha_e,\beta_e)と書く。U(e)=p(e)u′(e)U(e)=p(e)u'(e)、V(e)=p(e)v′(e)V(e)=p(e)v'(e)であるから、仮定はαeu(e)+βeU(e)=0\alpha_eu(e)+\beta_eU(e)=0かつαev(e)+βeV(e)=0\alpha_ev(e)+\beta_eV(e)=0である。定理 1.4 (2)によりH(a)=H(b)=0H(a)=H(b)=0であるから、[p(uv′‾−u′v‾)]ab=H(b)−H(a)=0\left[p\bigl(u\overline{v'}-u'\overline v\bigr)\right]_a^b=H(b)-H(a)=0であり、主張 (1)により積分は零である。▨

定理 4.2 (Lagrange の恒等式).a<ba<bとし、実数値関数p∈C1([a,b])p\in C^1([a,b])とq∈C([a,b])q\in C([a,b])がp(x)>0p(x)>0をすべてのx∈[a,b]x\in[a,b]で満たすとし、LLを式 (3.1.1)で定める。u,v∈C2([a,b];C)u,v\in C^2([a,b];\C)とする。このとき、[a,b][a,b]の各点で

uL[v]−vL[u]=ddx[p(x)(u′(x)v(x)−u(x)v′(x))]uL[v]-vL[u] =\frac{d}{dx} \left[ p(x)\bigl(u'(x)v(x)-u(x)v'(x)\bigr) \right]

が成り立ち、

∫ab(uL[v]−vL[u]) dx=[p(u′v−uv′)]ab\int_a^b\bigl(uL[v]-vL[u]\bigr)\,dx = \left[ p(u'v-uv') \right]_a^b(4.2.1)\text{(4.2.1)}

が成り立つ。

証明.pu′,pv′∈C1([a,b];C)pu',pv'\in C^1([a,b];\C)であるから、積の微分法則により[a,b][a,b]の各点で

ddx[p(u′v−uv′)]=(pu′)′v+pu′v′−pu′v′−u(pv′)′=(qu−L[u])v−u(qv−L[v])=uL[v]−vL[u]\frac{d}{dx}\bigl[p(u'v-uv')\bigr] =(pu')'v+pu'v'-pu'v'-u(pv')' =\bigl(qu-L[u]\bigr)v-u\bigl(qv-L[v]\bigr) =uL[v]-vL[u]

が成り立つ。p(u′v−uv′)p(u'v-uv')は[a,b][a,b]上で連続微分可能であるから、§E13.2 系 2.7 (2)をt0=at_0=a、t=bt=bで適用すると式 (4.2.1)を得る。▨

5 境界形式の消滅と対称性

定義 5.1.a<ba<bとし、w:[a,b]→Rw:[a,b]\to\Rを連続関数でw(x)>0w(x)>0をすべてのx∈[a,b]x\in[a,b]で満たすものとする。連続関数u,v:[a,b]→Cu,v:[a,b]\to\Cに対して、重み付き内積 (weighted inner product) を

⟨u,v⟩w=∫abu(x)v(x)‾w(x) dx\langle u,v\rangle_w =\int_a^b u(x)\overline{v(x)}w(x)\,dx(5.1.1)\text{(5.1.1)}

と定める。w>0w>0であるから、これはC([a,b];C)C([a,b];\C)上の正定値な内積である。

定理 5.2.a<ba<bとし、p,q,wp,q,wと境界係数を定義 3.2のとおりとする。u,v∈C2([a,b];C)u,v\in C^2([a,b];\C)がともに式 (3.2.1)を満たすとする。このとき

⟨w−1L[u],v⟩w=⟨u,w−1L[v]⟩w\left\langle w^{-1}L[u],v\right\rangle_w = \left\langle u,w^{-1}L[v]\right\rangle_w

が成り立つ。

証明.U=pu′U=pu'とV=pv′V=pv'は[a,b][a,b]上で連続微分可能であるから、§E13.2 系 2.7 (2)によりu,v,U,Vu,v,U,Vは[a,b][a,b]上で絶対連続であり、[a,b][a,b]の各点でu′=U/pu'=U/p、v′=V/pv'=V/pが成り立つ。ppは正の連続関数であり、qqは連続であるから[a,b][a,b]上で可積分である。したがって定理 4.1をu,vu,vに適用することができ、そこでのL[u],L[v]L[u],L[v]は式 (3.1.1)のL[u],L[v]L[u],L[v]に一致する。定理 4.1 (1)と定理 4.1 (2)により、式 (5.1.1)の定義から

⟨w−1L[u],v⟩w−⟨u,w−1L[v]⟩w=∫ab(L[u]v‾−uL[v]‾)dx=0\left\langle w^{-1}L[u],v\right\rangle_w -\left\langle u,w^{-1}L[v]\right\rangle_w = \int_a^b \left( L[u]\overline v-u\overline{L[v]} \right)dx =0

である。▨

6 固有値の実性と固有関数の直交性

系 6.1.a<ba<bとし、p,q,wp,q,wと境界係数を定義 3.2のとおりとする。このとき、式 (3.2.2)で定めた正則 Sturm–Liouville 問題の固有値は実数である。

証明.定義 3.2により、この問題の固有値は、両端が正則端点であるp,q,wp,q,wの(a,b)(a,b)への制限についての定義 1.6の固有値であるから、命題 1.7 (1)により実数である。▨

系 6.2.a<ba<bとし、p,q,wp,q,wと境界係数を定義 3.2のとおりとする。ym,yny_m,y_nを、それぞれ異なる固有値λm≠λn\lambda_m\ne\lambda_nに属する固有関数とする。このとき

∫abym(x)yn(x)‾w(x) dx=0\int_a^b y_m(x)\overline{y_n(x)}w(x)\,dx =0

が成り立つ。

証明.定義 3.2により、ym,yny_m,y_nの(a,b)(a,b)への制限は、両端が正則端点であるp,q,wp,q,wの(a,b)(a,b)への制限についての定義 1.6の意味で異なる固有値に属する固有関数であるから、主張は命題 1.7 (2)による。▨

例 6.3.b>0b>0とし、[0,b][0,b]上でp=w=1p=w=1、q=0q=0とすると、L[y]=−y′′L[y]=-y''であり、Dirichlet 条件y(0)=y(b)=0y(0)=y(b)=0の下のL[y]=λyL[y]=\lambda yは正則 Sturm–Liouville 固有値問題である。区間[0,b][0,b]について§E13.24 例 2.2で得た

λn=(nπb)2,yn(x)=sin⁡(nπxb)(n∈N≥1)\lambda_n=\left(\frac{n\pi}{b}\right)^2,\qquad y_n(x)=\sin\left(\frac{n\pi x}{b}\right)\qquad(n\in\NN)

は、補題 3.3 (2)によりこの問題の固有値と固有関数である。m≠nm\ne nならばλm≠λn\lambda_m\ne\lambda_nであるから、

∫0bsin⁡(mπxb)sin⁡(nπxb) dx=0(m≠n)\int_0^b \sin\left(\frac{m\pi x}{b}\right) \sin\left(\frac{n\pi x}{b}\right)\,dx=0 \qquad(m\ne n)

が系 6.2から従う。

注意 6.4. 正則 Sturm–Liouville 固有値問題の固有値の全体は、+∞+\inftyへ発散する狭義増加の実数列λ1<λ2<⋯\lambda_1<\lambda_2<\cdotsをなす。各λn\lambda_nに属する実数値の固有関数ene_nを⟨en,en⟩w=1\langle e_n,e_n\rangle_w=1となるように選ぶと、(en)n∈N≥1(e_n)_{n\in\NN}はL2([a,b],w(x) dx;C)L^2([a,b],w(x)\,dx;\C)の完全正規直交系である。これらの事実は §E13.27 固有関数展開 で証明し、本記事のどの証明にも用いない。

7 演習

問題 7.1 (Robin 条件の固有値). 区間[0,1][0,1]上でp=w=1p=w=1、q=0q=0とし、分離境界条件y(0)=0y(0)=0、y′(1)+y(1)=0y'(1)+y(1)=0の下で固有値問題−y′′=λy-y''=\lambda yを考える。固有値はすべて正であることを示せ。また、λ\lambdaが固有値であることと、k>0k>0とtan⁡k=−k\tan k=-kを満たすkkによってλ=k2\lambda=k^{2}と書くことができることが同値であることを示せ。さらに、各m∈N≥1m\in\NNについて区間((2m−1)π/2,mπ)\bigl((2m-1)\pi/2,m\pi\bigr)にtan⁡k=−k\tan k=-kの解がちょうど一つあり、(0,∞)(0,\infty)の解はこれらに限ることを示せ。

解答.

境界条件は実係数の分離型であるから、系 6.1により固有値は実数である。λ\lambdaを固有値、yyを対応する固有関数とすると、補題 3.3 (2)によりyyは零でないC2C^{2}級関数である。部分積分とy(0)=0y(0)=0、y′(1)=−y(1)y'(1)=-y(1)により

λ∫01∣y∣2 dx=∫01(−y′′)y‾ dx=[−y′y‾]01+∫01∣y′∣2 dx=∣y(1)∣2+∫01∣y′∣2 dx\lambda\int_0^1|y|^{2}\,dx=\int_0^1(-y'')\overline y\,dx=\bigl[-y'\overline y\bigr]_0^1+\int_0^1|y'|^{2}\,dx=|y(1)|^{2}+\int_0^1|y'|^{2}\,dx

である。右辺は00以上である。右辺が00ならばy′≡0y'\equiv0であり、yyは定数であり、y(0)=0y(0)=0からy≡0y\equiv0である。これはyyが零でないことに反する。したがって右辺は正であり、∫01∣y∣2>0\int_0^1|y|^{2}>0であるからλ>0\lambda>0である。

λ=k2\lambda=k^{2}(k>0k>0)とすると、−y′′=k2y-y''=k^{2}yの解はy=Acos⁡kx+Bsin⁡kxy=A\cos kx+B\sin kxである。y(0)=0y(0)=0からA=0A=0であり、y′(1)+y(1)=B(kcos⁡k+sin⁡k)y'(1)+y(1)=B(k\cos k+\sin k)である。非零解が存在することはkcos⁡k+sin⁡k=0k\cos k+\sin k=0と同値である。cos⁡k=0\cos k=0ならばsin⁡k=±1\sin k=\pm1となってこの等式は成り立たないので、等式はcos⁡k≠0\cos k\ne0かつtan⁡k=−k\tan k=-kと同値である。このときsin⁡kx\sin kxが固有関数である。前段と合わせて同値性を得る。

k>0k>0でcos⁡k≠0\cos k\ne0を満たす点は、区間(0,π/2)(0,\pi/2)と、m∈N≥1m\in\NNについての区間((2m−1)π/2,mπ)\bigl((2m-1)\pi/2,m\pi\bigr)と[mπ,(2m+1)π/2)\bigl[m\pi,(2m+1)\pi/2\bigr)に分かれる。(0,π/2)(0,\pi/2)と[mπ,(2m+1)π/2)\bigl[m\pi,(2m+1)\pi/2\bigr)ではtan⁡k≥0>−k\tan k\ge0>-kであるから解はない。Im=((2m−1)π/2,mπ)I_{m}=\bigl((2m-1)\pi/2,m\pi\bigr)の上でg(k)=tan⁡k+kg(k)=\tan k+kとおくと、g′(k)=1/cos⁡2k+1>0g'(k)=1/\cos^{2}k+1>0であるからggは狭義単調増加であり、k↓(2m−1)π/2k\downarrow(2m-1)\pi/2でg(k)→−∞g(k)\to-\infty、g(mπ)=mπ>0g(m\pi)=m\pi>0である。中間値の定理によりImI_{m}にggの零点があり、単調性によりただ一つである。▨

問題 7.2 (非一定の重みをもつ固有関数). 区間[0,π][0,\pi]上で

p(x)=w(x)=e2x,q(x)=0p(x)=w(x)=e^{2x},\qquad q(x)=0

とし、Dirichlet 条件y(0)=y(π)=0y(0)=y(\pi)=0の下の正則 Sturm–Liouville 固有値問題を考える。定理 2.1を用いて、固有値の全体が{n2+1:n∈N≥1}\{n^2+1:n\in\NN\}であり、固有値n2+1n^2+1に属する固有関数がyn(x)=e−xsin⁡(nx)y_n(x)=e^{-x}\sin(nx)の零でない複素数倍に限ることを示せ。さらに、m≠nm\ne nのとき⟨ym,yn⟩w=0\langle y_m,y_n\rangle_w=0であることを、定理 2.1 (2)のユニタリ写像を用いて計算によって確かめよ。

解答.

p,wp,wはR\R上の正のC∞C^\infty級関数e2xe^{2x}の制限である。定理 2.1を(a,b)=(0,π)(a,b)=(0,\pi)に適用する。w/p=1\sqrt{w/p}=1であるからt(x)=xt(x)=xと取ると、ta=0t_a=0、tb=πt_b=\pi、x(τ)=τx(\tau)=\tauであり、m(τ)=(e4τ)1/4=eτm(\tau)=(e^{4\tau})^{1/4}=e^\tau、m′′/m=1m''/m=1であるからV=q/w+m′′/m=1V=q/w+m''/m=1である。定理 2.1 (3)により、yyがℓ[y]=λwy\ell[y]=\lambda wyの解であることと、u=eτyu=e^\tau yがu∈C1u\in C^1、u′∈AClocu'\in AC_{\mathrm{loc}}を満たし−u′′+u=λu-u''+u=\lambda uがほとんど至る所で成り立つことは同値である。後者は、uuが係数1,0,11,0,1についての方程式ℓ[u]=(λ−1)u\ell[u]=(\lambda-1)uの、流束u′u'をもつ解であることと同じである。端点e∈{0,π}e\in\{0,\pi\}でp,wp,wは正のC2C^2級関数に延び、q=0q=0であるから、定理 2.1 (4)を(α,β)=(1,0)(\alpha,\beta)=(1,0)に適用すると、y(e)=0y(e)=0はu(e)=0u(e)=0と同値である。定義 3.2と定義 1.6により、λ\lambdaが固有値であることは、u(0)=u(π)=0u(0)=u(\pi)=0を満たすℓ[u]=(λ−1)u\ell[u]=(\lambda-1)u(係数1,0,11,0,1)の零でない解uuが存在することと同値であり、固有関数はy=e−xuy=e^{-x}uの形のものに限る。

系 6.1により固有値は実数であるから、μ=λ−1∈R\mu=\lambda-1\in\Rとする。係数1,0,11,0,1についてのℓ[u]=μu\ell[u]=\mu uの解は、補題 1.5 (3)により(u(0),u′(0))(u(0),u'(0))によってただ一つに定まる。μ>0\mu>0のときk=μk=\sqrt\muとおくと、cos⁡kx\cos kxとsin⁡kx\sin kxは解であるから、(u(0),u′(0))=(A,kB)(u(0),u'(0))=(A,kB)の解はAcos⁡kx+Bsin⁡kxA\cos kx+B\sin kxである。u(0)=0u(0)=0からA=0A=0であり、u(π)=Bsin⁡kπu(\pi)=B\sin k\piであるから、B≠0B\ne0である解がu(π)=0u(\pi)=0を満たすことはk∈N≥1k\in\NNと同値である。μ=0\mu=0のとき解はA+BxA+Bxの形であり、u(0)=u(π)=0u(0)=u(\pi)=0からA=B=0A=B=0である。μ<0\mu<0のときκ=−μ\kappa=\sqrt{-\mu}とおくと解はAcosh⁡κx+Bsinh⁡κxA\cosh\kappa x+B\sinh\kappa xの形であり、u(0)=0u(0)=0からA=0A=0、sinh⁡κπ>0\sinh\kappa\pi>0とu(π)=0u(\pi)=0からB=0B=0である。したがって固有値の全体は{n2+1:n∈N≥1}\{n^2+1:n\in\NN\}であり、n2+1n^2+1に属する固有関数はB≠0B\ne0についてBe−xsin⁡(nx)=BynBe^{-x}\sin(nx)=By_nに限る。

UUを定理 2.1 (2)のユニタリ写像とすると(Uyn)(τ)=eτe−τsin⁡(nτ)=sin⁡(nτ)(Uy_n)(\tau)=e^\tau e^{-\tau}\sin(n\tau)=\sin(n\tau)であり、UUは内積を保つから、m≠nm\ne nのとき

⟨ym,yn⟩w=∫0πsin⁡(mτ)sin⁡(nτ) dτ=12[sin⁡((m−n)τ)m−n−sin⁡((m+n)τ)m+n]0π=0\langle y_m,y_n\rangle_w =\int_0^\pi\sin(m\tau)\sin(n\tau)\,d\tau =\frac12\left[ \frac{\sin((m-n)\tau)}{m-n} -\frac{\sin((m+n)\tau)}{m+n} \right]_{0}^{\pi} =0

である。この値は系 6.2の結論と一致する。▨

問題 7.3 (複素境界条件による対称性の破れ). 区間[0,1][0,1]上でp=w=1p=w=1、q=0q=0とし、

L[y]=−y′′,y′(0)−(1+i)y(0)=0,y′(1)−(1+i)y(1)=0L[y]=-y'',\qquad y'(0)-(1+i)y(0)=0,\qquad y'(1)-(1+i)y(1)=0

とする。境界条件を満たす零でないy∈C2([0,1];C)y\in C^2([0,1];\C)と実数でないλ\lambdaでL[y]=λyL[y]=\lambda yを満たすものを一つ挙げ、そのyyについてu=v=yu=v=yとした境界形式[yy′‾−y′y‾]01\bigl[y\overline{y'}-y'\overline y\bigr]_0^1が零でないことを示せ。境界係数を実数とする仮定を外すと、定理 4.1 (2)と系 6.1の結論が成り立たないことがわかる。

解答.

y(x)=e(1+i)xy(x)=e^{(1+i)x}と置くと、両端で

y′(x)−(1+i)y(x)=0y'(x)-(1+i)y(x)=0

である。また、

L[y]=−(1+i)2y=−2iyL[y]=-(1+i)^2y=-2iy

であるから、yyは境界条件を満たす零でないC2C^2級関数であり、実数でないλ=−2i\lambda=-2iについてL[y]=λyL[y]=\lambda yを満たす。

さらに、u=v=yu=v=yとして境界形式を計算すると、

yy′‾−y′y‾=∣y∣2((1−i)−(1+i))=−2i∣y∣2y\overline{y'}-y'\overline y =|y|^2\bigl((1-i)-(1+i)\bigr) =-2i|y|^2

である。したがって、

[yy′‾−y′y‾]01=−2i(e2−1)≠0.\left[y\overline{y'}-y'\overline y\right]_0^1 =-2i(e^2-1)\ne0.

y‾\overline yはy‾′−(1+i)y‾=−2iy‾\overline y'-(1+i)\overline y=-2i\overline yを満たすから、端点での値の組(y(e)‾,y′(e)‾)(\overline{y(e)},\overline{y'(e)})(e∈{0,1}e\in\{0,1\})は線形形式(z0,z1)↦z1−(1+i)z0(z_0,z_1)\mapsto z_1-(1+i)z_0の核に属さず、定理 1.4 (2)の証明の一次従属性の議論は適用することができない。▨

参考文献

  1. Anton Zettl, Sturm-Liouville Theory, Mathematical Surveys and Monographs 121, American Mathematical Society, 2005.可積分係数の Sturm–Liouville 式と流束、正則端点と特異端点の判定、古典的な方程式の Sturm–Liouville 形、正則問題の定義、Lagrange の恒等式、分離型自己共役境界条件を参考にした。
  2. E. C. Titchmarsh, Eigenfunction Expansions Associated with Second-Order Differential Equations, 2nd ed., vol. I, Oxford University Press, Oxford, 1962.正則 Sturm–Liouville 問題における固有値の実性と固有関数の直交性を参考にした。

前提記事

11 本の記事・単元を表示