身近な自然と科学

微分方程式解法入門 部分分数分解

例として、非線形微分方程式 $\dfrac{dy}{dx}=ay-by^{2}$ を解きます。
ただし、$a$と$b$は定数です。
この式は、ロジスティクス微分方程式と呼ばれ、1838年ベルハルストが人口論を解いたときに使ったものです。
人口は時間の経過と共に増えるがやがて急ブレーキが掛かって飽和するというものです。 人口が増えれば食料などの資源の供給が間に合わなくなって養えなくなりますから人口は増えなくなります。 人の世界で考えると嫌な気持ちになりますが、細菌や雑草でも限られた空間内では無限に増え続けることは出来ないので、当たり前の事を言っている微分方程式です。
この $\dfrac{dy}{dx}=ay-by^{2}$ の微分方程式の解法は易しい部類に属しますが、数式に弱い私には何をやっているのか判らなくなるほど複雑でした。
微分方程式と解く前に、解くのに使うテクニックの部分分数分解について(部分分数分解は、元の式を簡単にして微積分をしやすくする役割があります。)

部分分数分解とは

たとえば、分数 $\dfrac{1}{x^{2}-3x+2}$ を $\dfrac{-1}{x-1}+\dfrac{1}{x-2}$ の形に変えることです。
この例の部分分数分解は簡単なのですが、複雑な式になると才能に恵まれた方はよいのですが凡才では、数式の型毎に効率よく分解する方法を習熟するしかないです。
簡単な例ですが
$$\displaystyle \frac{1}{x^{2}-3x+2}=\displaystyle \frac{1}{(x-1)(x-2)} \tag{・・・式1}$$
式1の右辺を下 式2の形にしたい
$$\displaystyle \frac{a}{x-1}+\frac{b}{x-2} \tag{・・・式2}$$
式2を通分し
$$\displaystyle \frac{ax+bx-2a-b}{(x-1)(x-2)}=\frac{(a+b)x-(2a+b)}{(x-1)(x-2)} \tag{・・・式3}$$
式3の分子は$1$なので
$$(a+b)x-(2a+b)=1$$ $$a+b=0$$ $$2a+b=-1$$ $$a=-1$$ $$b=1$$
よって、
$$\displaystyle \frac{1}{x^{2}-3x+2}=\frac{-1}{x-1}+\frac{1}{x-2}$$

ロジスティクス微分方程式の解法

最も難解なのは、部分分数分解です。教えて貰えば「なるほど簡単」ですが
式5が解くロジスティクス微分方程式です。
$$\displaystyle \frac{dy}{dx}=ay-by^{2} \tag{・・・式5}$$
変数分離の方法で解くために右辺を整理してみます。
$$\displaystyle \frac{dy}{dx}=ay(1-\frac{b}{a}y) \tag{・・・式6}$$
少しでも見やすいように$M=\displaystyle \frac{b}{a}$として
$$\displaystyle \frac{dy}{dx}=ay(1-My) \tag{・・・式7}$$ $$\displaystyle \frac{1}{y(1-My)}\frac{dy}{dx}=a\tag{・・・式8}$$
ここで、式8の左辺を部分分数分解します。
$$(\displaystyle \frac{1}{y}+\frac{M}{1-My})\frac{dy}{dx}=a\tag{・・・式9}$$ $$\displaystyle \int(\frac{1}{y}+\frac{M}{1-My})\frac{dy}{dx}dx=\int adx\tag{・・・式10}$$ $$\displaystyle \int\frac{1}{y}dy+M\int\frac{1}{1-My}dy=\int adx\tag{・・・式11}$$
式11の左辺2項目は置換積分法を使います
$$ln(y)+C_{1}-ln(1-My)+C_{2}=ax+C_{3}\tag{・・・式12}$$
対数同士の引き算は割り算になるので次の様に変形できます
積分定数は$C_{4}$にまとめました
$$ln(\displaystyle \frac{y}{1-My})=ax+C_{4}\tag{・・・式13}$$
対数を指数に直します
$$e^{ax+c_{4}}=\displaystyle \frac{y}{1-My}\tag{・・・式14}$$ $$e^{C_{4}}e^{ax}=\displaystyle \frac{y}{1-My} \tag{・・・式15}$$
$e^{C_{4}}$は定数なので、積分定数$C$と置き換えます
$$Ce^{ax}=\displaystyle \frac{y}{1-My}\tag{・・・式16}$$
式16を$y$について解き
$$y=\displaystyle \frac{Ce^{ax}}{1+CMe^{ax}}\tag{・・・式17}$$
$M$を元に戻してロジスティクス微分方程式の一般解が求められました
$$y=\displaystyle \frac{Ce^{ax}}{1+C\frac{b}{a}e^{ax}}\tag{・・・式18}$$