演習で使う数学的ツール - 数理生物学演習

この演習で必要となる数学的なツールを整理する．

変数分離で微分方程式を解く¶

変数分離法は，常微分方程式の最も基本的な解法の一つ．

$\frac{dx}{dt} = f(x) \cdot g(t)$

の形に書ける微分方程式に対して， $x$ を含む項と $t$ を含む項をそれぞれ左辺・右辺に分離し，両辺を積分する．

\frac{dx}{dt} = rx

変数分離すると，

\frac{1}{x} dx = r \, dt

両辺を積分して，

\int \frac{1}{x} dN = \int r \, dt

\log x = rt + C

x(t) = x_0 e^{rt}

ここで $x_0 = e^C$ は初期条件 $x(0) = x_0$ から決まる．

多変数の連立微分方程式系

\frac{dx_1}{dt} = f_1(x_1, x_2, \ldots, x_n), \quad \frac{dx_2}{dt} = f_2(x_1, x_2, \ldots, x_n), \quad \ldots

に対して，各関数の各変数に関する偏微分を並べた行列

J = \begin{pmatrix} \frac{\partial f_1}{\partial x_1} & \frac{\partial f_1}{\partial x_2} & \cdots \\ \frac{\partial f_2}{\partial x_1} & \frac{\partial f_2}{\partial x_2} & \cdots \\ \vdots & \vdots & \ddots \end{pmatrix}

をヤコビ行列という．

平衡点 $(x_1^*, x_2^*, \ldots, x_n^*)$ でこのヤコビ行列を評価し，その固有値から平衡点の局所安定性を判定する．

正方行列 $A$ に対して，

A \mathbf{v} = \lambda \mathbf{v}

を満たすスカラー $\lambda$ を固有値，ベクトル $\mathbf{v}$ を固有ベクトルと呼ぶ．

固有値は特性方程式

\det(A - \lambda I) = 0

を解くことで求められる．

A = \begin{pmatrix} 1 & 2 \\ 3 & 4 \end{pmatrix}

特性方程式は $(1-\lambda)(4-\lambda) - 6 = 0$ ，すなわち $\lambda^2 - 5\lambda - 2 = 0$ である．

平衡点におけるヤコビ行列の固有値 $\lambda_1, \lambda_2, \ldots$ から，平衡点の安定性がわかる．