陰公式で最も扱いやすいのは
としたときで、
(8) は
![$\displaystyle -\alpha U_{j-1}^{n+1} +(1+2\alpha)U_j^{n+1} -\alpha U_{j+1}^{n+1} =U_j^n$](Diffusion_equation-img40.png) |
(11) |
と書き換えることができ、
(9)と(10)は
![$\displaystyle \begin{pmatrix}1+2\alpha & -\alpha & & & & 0 -\alpha & 1+2\al...
... U_2^{n} \vdots \vdots U_{J-2}^{n} U_{J-1}^{n} \end{pmatrix}$](Diffusion_equation-img41.png) |
(12) |
と書けるので、
時刻
での
が
から
まで与えられたとすれば、
(12)の両辺に大きな行列
の逆行列を掛けてやれば、
時刻
での
が求められることになる。
ここで、境界条件(4) は
、
として表している。
大きな行列の逆行列を求めるのは一般には必ずしも容易ではないし、
行列
の様に帯対角成分以外は零であるような行列をそのままひっくり返すのは、
経済的にも余り良い方法であるとは言えない。
ここでは行列
が三重対角行列となっているので、
漸化式(11)を用いて解くことを考える。
(11)は形式的に
![$\displaystyle P_j U_{j-1}^{n+1} +Q_j U_j^{n+1} +R_j U_{j+1}^{n+1} =U_j^n$](Diffusion_equation-img50.png) |
(13) |
と書ける。ここで
![$\displaystyle P_j = -\alpha , \quad Q_j =1+2\alpha , \quad R_j = -\alpha$](Diffusion_equation-img51.png) |
(14) |
であり、
![$\displaystyle P_1 =0 ,\quad R_{J-1} =0$](Diffusion_equation-img52.png) |
(15) |
とした。
(13)に於いて、
と
とした二つの式を使って
を消去すれば、
![$\displaystyle \left(Q_2 -\frac{P_2 R_1}{Q_1}\right) U_2^{n+1} +R_2 U_3^{n+1} =U_2^n -\frac{P_2U_1^n}{Q_1}$](Diffusion_equation-img55.png) |
(16) |
となり、または
![$\displaystyle \tilde{Q}_2 U_2^{n+1} + R_2 U_3^{n+1} =\tilde{U}_2^n$](Diffusion_equation-img56.png) |
(17) |
を得る。ここで
![$\displaystyle \tilde{Q}_2 = Q_2 -\frac{P_2 R_1}{Q_1} ,\quad \tilde{U}_2^n = U_2^n -\frac{P_2U_1^n}{Q_1}$](Diffusion_equation-img57.png) |
(18) |
である。これと、(13)の
とした式とを使って今度は
を消去すれば
![$\displaystyle \tilde{Q}_3 U_3^{n+1} + R_3 U_4^{n+1} =\tilde{U}_3^n$](Diffusion_equation-img60.png) |
(19) |
を得る。ここで
![$\displaystyle \tilde{Q}_3 = Q_3 -\frac{P_3 R_2}{\tilde{Q}_2} ,\quad \tilde{U}_3^n = U_3^n -\frac{P_3\tilde{U}_1^2}{\tilde{Q}_2}$](Diffusion_equation-img61.png) |
(20) |
である。同様の操作を繰り返せば
![$\displaystyle \tilde{Q}_j U_j^{n+1} + R_j U_{j+1}^{n+1} =\tilde{U}_j^n$](Diffusion_equation-img62.png) |
(21) |
![$\displaystyle \tilde{Q}_j = Q_j -\frac{P_j R_{j-1}}{\tilde{Q}_{j-1}} ,\quad \tilde{U}_j^n = U_j^n -\frac{P_j\tilde{U}_{j-1}^n}{\tilde{Q}_{j-1}}$](Diffusion_equation-img63.png) |
(22) |
を得る。これを
まで繰り返せば
![$\displaystyle \tilde{Q}_{J-2} U_{J-2}^{n+1} + R_{J-2} U_{J-1}^{n+1} =\tilde{U}_{J-2}^n$](Diffusion_equation-img65.png) |
(23) |
となるが、これと(13)で
とした式
![$\displaystyle P_{J-1} U_{J-2}^{n+1} + Q_{J-1}U_{J-1}^{n+1} =U_{J-1}^n$](Diffusion_equation-img66.png) |
(24) |
とを連立させて解けば、
![% latex2html id marker 1289
$\displaystyle \tilde{Q}_{J-1} U_{J-1}^{n+1} = \tild...
...^n \quad \therefore U_{J-1}^{n+1} = \frac{\tilde{U}_{J-1}^n}{\tilde{Q}_{J-1}}$](Diffusion_equation-img67.png) |
(25) |
を得る。
このように求められた
と漸化式(21)とを使えば、
![$\displaystyle U_j^{n+1} = \frac{\tilde{U}_j^n -R_j U_{j+1}^{n+1}}{\tilde{Q}_j}$](Diffusion_equation-img69.png) |
(26) |
より、
から
までの
を求めることができる。
fat-cat
平成16年11月30日