次へ: 7 数値計算 上へ: 糸の方程式の境界条件について 前へ: 5 ベッセルの方程式への帰着と境界条件 (PDF ファイル: lstring1.pdf)


6 一般の μ の場合

それでは、一般の $\mu $ の場合はどうなるのかを次に考えてみることにする。

5 節の $\mu =1$ の場合の議論を見れば、 $x=0$ での境界条件の不要性は、 結局 $w$ の $x=0$ での有界性から (20) の $C_4=0$ が 導かれることに起因することがわかる。 よって、一般の $\mu $ の場合にも 3 節、5 節と同様の変数分離法を用いて、 $x$ に関する方程式をベッセルの方程式に帰着させることで考察してみることにする。

方程式 (5) に対して $u=v(t)w(x)$ とすると、

\begin{displaymath} \frac{v''(t)}{v(t)}=\frac{xw''(x)+\mu w'(x)}{w(x)}=-\lambda^2 \end{displaymath}

となるので、
\begin{eqnarray*}&& v(t)=C_1\cos\lambda t+C_2\sin\lambda t\hspace{1zw}(t>0)\\ && xw''(x)+\mu w'(x)+\lambda^2 w(x)=0\hspace{1zw}(0<x<L)\end{eqnarray*}


となる。 5 節と同様に $w$ を (18) のように置いて変形すると、この場合は
\begin{displaymath} \hat{x}^2\hat{w}''+\frac{2\alpha+\beta+\mu-1}{\beta}\hat{x}... ...w}' +\frac{\alpha(\alpha+\mu-1)+\lambda^2 x}{\beta^2}\hat{w}=0\end{displaymath} (26)

となるので、よって
\begin{displaymath} \beta=\frac{1}{2},\hspace{1zw}\gamma=2\lambda,\hspace{1zw}\alpha=\frac{1-\mu}{2} \end{displaymath}

とすることで
\begin{displaymath} \frac{\lambda^2}{\beta^2}x=\hat{x}^2, \hspace{1zw}\frac{2\al... ...1, \hspace{1zw}\frac{\alpha(\alpha+\mu-1)}{\beta^2}=-(1-\mu)^2 \end{displaymath}

となって、よって (26) は
\begin{displaymath} \hat{x}^2\hat{w}''+\hat{x}\hat{w}'+\{\hat{x}^2-(1-\mu)^2\}\hat{w}=0 \end{displaymath}

の形、すなわち $\nu=\vert 1-\mu\vert$ のベッセルの微分方程式となる。

まず $\mu $ が整数でない場合は、$\nu=\vert 1-\mu\vert$ も整数ではなく、 よって $\hat{w}$ は

\begin{displaymath} \hat{w}=C_3J_\nu(\hat{x})+C_4J_{-\nu}(\hat{x}) \end{displaymath}

と書けるので、これにより $w(x)$ は
\begin{displaymath} w(x)=x^{(1-\mu)/2} \{C_3J_\nu(2\lambda\sqrt{x})+C_4J_{-\nu}(2\lambda\sqrt{x})\}\end{displaymath} (27)

となるが、(13) より $x^{(1-\mu)/2}J_\nu(2\lambda\sqrt{x})$ は
\begin{displaymath} x^{(1-\mu)/2}J_\nu(2\lambda\sqrt{x}) = x^{(1-\mu)/2}(\sqrt{x})^\nu O(1) = x^{(1-\mu+\vert 1-\mu\vert)/2}O(1) \end{displaymath}

となるので、これはどんな $\mu $ に対しても $x=0$ の近くで有界となる。 一方、(14) より $x^{(1-\mu)/2}J_{-\nu}(2\lambda\sqrt{x})$ は
\begin{displaymath} x^{(1-\mu)/2}J_{-\nu}(2\lambda\sqrt{x}) = x^{(1-\mu-\vert 1-\mu\vert)/2}O(1) \end{displaymath}

となるので、これは $\mu<1$ ならば有界であるが、$\mu>1$ ならば有界ではない。 つまり、$\mu<1$ ならば (27) の $C_4$ は残るが、 $\mu>1$ ならば $C_4=0$ ではなくてはならないことになる。

次に $\mu $ が整数の場合は、$\nu=\vert 1-\mu\vert$ も整数なので、

\begin{displaymath} \hat{w}=C_3J_\nu(\hat{x})+C_4Y_\nu(\hat{x}) \end{displaymath}

より、
\begin{displaymath} w(x)=x^{(1-\mu)/2} \{C_3J_\nu(2\lambda\sqrt{x})+C_4Y_\nu(2\lambda\sqrt{x})\}\end{displaymath} (28)

となる。 $x^{(1-\mu)/2}J_\nu(2\lambda\sqrt{x})$ は $\mu $ が整数でない場合と同様に $x=0$ で常に有界であるが、 $x^{(1-\mu)/2}Y_\nu(2\lambda\sqrt{x})$ は (15) より、 $\mu =1$ のときは
\begin{displaymath} x^0Y_0(2\lambda\sqrt{x})=O(\log x) \end{displaymath}

より有界ではなく、$\mu\neq 1$ のときは、
\begin{displaymath} x^{(1-\mu)/2}Y_\nu(2\lambda\sqrt{x}) = x^{(1-\mu)/2}(\sqrt{x})^{-\nu} O(1) = x^{(1-\mu-\vert 1-\mu\vert)/2}O(1) \end{displaymath}

となるので、やはり $\mu<1$ ならば有界、$\mu>1$ ならば有界ではない。 よって、$\mu $ が整数の場合も $\mu<1$ ならば $w(x)$ は有界なので (28) の $C_4$ は残るが、 $\mu\geq 1$ ならば有界ではないので $C_4=0$ でなくてはならない。

結局以上をまとめると、以下のようになる。

次へ: 7 数値計算 上へ: 糸の方程式の境界条件について 前へ: 5 ベッセルの方程式への帰着と境界条件
竹野茂治＠新潟工科大学
2009年6月22日