next up previous
Next: 円盤領域における熱方程式 Up: A. Bessel 関数早分かり Previous: A..5 Fourier-Bessel 展開

A..6 応用1: 円盤領域のラプラシアンの固有値問題

$ \Omega=\{(x,y)\in\R^2; x^2+y^2<1\}$ におけるラプラシアンの固有値問題

(6)     $\displaystyle -\Laplacian u=\lambda u$   $\displaystyle \mbox{(in $\Omega$)}$$\displaystyle ,$
(7)     $\displaystyle u=0$   $\displaystyle \mbox{(on $\rd\Omega$)}$$\displaystyle ,$
(8)     $\displaystyle u\not\equiv 0$

を考える。 すなわちこの3つの方程式をみたす関数 $ u$ と 定数 $ \lambda$ を求める。$ \lambda$ を固有値、 $ u$$ \lambda$ に属する固有関数とよぶ。

極座標変換

$\displaystyle x=r\cos\theta,\quad
y=r\sin\theta,\quad
U(r,\theta)=u(x,y)
$

を導入すると
(9)     $\displaystyle U_{rr}+\frac{1}{r}U_r+\frac{1}{r^2}U_{\theta\theta}=-\lambda U$   $\displaystyle \mbox{($0<r<1$, $\theta\in[0,2\pi]$)}$$\displaystyle ,$
(10)     $\displaystyle U(1,\theta)=0$   $\displaystyle \mbox{($\theta\in[0,2\pi]$)}$$\displaystyle ,$
(11)     $\displaystyle U\not\equiv 0.$

変数分離解を求める。すなわち

$\displaystyle U(r,\theta)=R(r)\Theta(\theta)
$

の形をしているものを求める。 (9) に代入して $ R(r)\Theta(\theta)$ で割ると

$\displaystyle \frac{R''}{R}+\frac{1}{r}\frac{R'}{R}
+\frac{1}{r^2}\frac{\Theta''}{\Theta}
=-\lambda.
$

移項して

$\displaystyle \frac{r^2 R''+r R'+\lambda r^2 R}{R}
=-\frac{\Theta''}{\Theta}.
$

明らかにこの等式の値は定数である。それを $ \mu$ とおくと、

$\displaystyle \Theta''=-\mu\Theta.
$

$ \Theta$ が周期 $ 2\pi$ の関数であることを考えると、

$\displaystyle (\exists n\in\N\cup\{0\})\quad
(\exists A,B\in\C)$   s.t.$\displaystyle \quad
\mu=n^2,\quad\Theta(\theta)=A\cos n\theta+B\sin n\theta.
$

$ \mu=n^2$ のとき、$ R(r)$ については

$\displaystyle r^2 R''+r R'+(\lambda r^2-n^2)R=0.
$

すなわち

$\displaystyle R''+\frac{1}{r}R'+\left(\lambda-\frac{n^2}{r^2}\right)R=0.
$

ここで $ \lambda>0$ を仮定して

$\displaystyle z=\sqrt{\lambda}r,\quad
R(r)=f(z)
$

とおくと、

$\displaystyle \lambda f''\left(\sqrt{\lambda}r\right)
+\frac{\sqrt{\lambda}}{r}...
...\right)
+\left(\lambda-\frac{n^2}{r^2}\right)
f\left(\sqrt{\lambda}r\right)=0.
$

ゆえに

$\displaystyle f''(z)+\frac{1}{z}f'(z)+\left(1-\frac{n^2}{z^2}\right)f(z)=0.
$

これは Bessel の微分方程式である。 $ f(z)$$ z=0$ で有界であるから、

$\displaystyle \exists C\in\C$   s.t.$\displaystyle \quad
f(z)=C J_n(z).
$

すなわち

$\displaystyle R(r)=f(z)=f\left(\sqrt{\lambda}r\right)=C J_n(\sqrt{\lambda}r).
$

ところで $ R(1)=0$ であるから、

$\displaystyle C J_n(\sqrt{\lambda})=0.
$

$ C\ne 0$ であるから $ \sqrt{\lambda}$$ J_n$ の零点である。

$\displaystyle \exists m\in\N$   s.t.$\displaystyle \quad \sqrt{\lambda}=\mu_{nm}.
$

これから $ \lambda=\mu_{nm}^2$.


\begin{jremark}[$\lambda>0$ を仮定したことについて]
$-\Laplacian u=\lambda u$\...
...r
一方、計算だけで $\lambda>0$ を示すこともできる
(工事中)。 \qed
\end{jremark}

以上から、 $ \lambda:=\mu_{nm}$ は円盤領域におけるラプラシアンの固有値であり、

$\displaystyle u_{nm}(r,\theta):=
J_n(\mu_{nm}r)(A'\cos n\theta+B'\sin n\theta)
$

$ \lambda$ に属する固有関数であることが分かる。

実はこれ以外に固有値はなく、 固有関数もここに現われるものだけで十分であることが分かる。

(工事中)




next up previous
Next: 円盤領域における熱方程式 Up: A. Bessel 関数早分かり Previous: A..5 Fourier-Bessel 展開
Masashi Katsurada
平成18年11月21日