微分方程式解法入門 部分分数分解
例として、非線形微分方程式 $\dfrac{dy}{dx}=ay-by^{2}$ を解きます。
ただし、$a$と$b$は定数です。
この式は、ロジスティクス微分方程式と呼ばれ、１８３８年ベルハルストが人口論を解いたときに使ったものです。
人口は時間の経過と共に増えるがやがて急ブレーキが掛かって飽和するというものです。 人口が増えれば食料などの資源の供給が間に合わなくなって養えなくなりますから人口は増えなくなります。 人の世界で考えると嫌な気持ちになりますが、細菌や雑草でも限られた空間内では無限に増え続けることは出来ないので、当たり前の事を言っている微分方程式です。
この $\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{・・・式１}$$
式１の右辺を下 式２の形にしたい $$\displaystyle \frac{a}{x-1}+\frac{b}{x-2} \tag{・・・式２}$$
式２を通分し $$\displaystyle \frac{ax+bx-2a-b}{(x-1)(x-2)}=\frac{(a+b)x-(2a+b)}{(x-1)(x-2)} \tag{・・・式３}$$
式３の分子は$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}$$
ロジスティクス微分方程式の解法
最も難解なのは、部分分数分解です。教えて貰えば「なるほど簡単」ですが
式５が解くロジスティクス微分方程式です。
積分定数は$C_{4}$にまとめました
式５が解くロジスティクス微分方程式です。
$$\displaystyle \frac{dy}{dx}=ay-by^{2} \tag{・・・式５}$$
変数分離の方法で解くために右辺を整理してみます。 $$\displaystyle \frac{dy}{dx}=ay(1-\frac{b}{a}y) \tag{・・・式６}$$
少しでも見やすいように$M=\displaystyle \frac{b}{a}$として $$\displaystyle \frac{dy}{dx}=ay(1-My) \tag{・・・式７}$$ $$\displaystyle \frac{1}{y(1-My)}\frac{dy}{dx}=a\tag{・・・式８}$$
ここで、式８の左辺を部分分数分解します。 $$(\displaystyle \frac{1}{y}+\frac{M}{1-My})\frac{dy}{dx}=a\tag{・・・式９}$$ $$\displaystyle \int(\frac{1}{y}+\frac{M}{1-My})\frac{dy}{dx}dx=\int adx\tag{・・・式１０}$$ $$\displaystyle \int\frac{1}{y}dy+M\int\frac{1}{1-My}dy=\int adx\tag{・・・式１１}$$
式１１の左辺２項目は置換積分法を使います $$ln(y)+C_{1}-ln(1-My)+C_{2}=ax+C_{3}\tag{・・・式１２}$$
対数同士の引き算は割り算になるので次の様に変形できます積分定数は$C_{4}$にまとめました
$$ln(\displaystyle \frac{y}{1-My})=ax+C_{4}\tag{・・・式１３}$$
対数を指数に直します $$e^{ax+c_{4}}=\displaystyle \frac{y}{1-My}\tag{・・・式１４}$$ $$e^{C_{4}}e^{ax}=\displaystyle \frac{y}{1-My} \tag{・・・式１５}$$
$e^{C_{4}}$は定数なので、積分定数$C$と置き換えます $$Ce^{ax}=\displaystyle \frac{y}{1-My}\tag{・・・式１６}$$
式１６を$y$について解き $$y=\displaystyle \frac{Ce^{ax}}{1+CMe^{ax}}\tag{・・・式１７}$$
$M$を元に戻してロジスティクス微分方程式の一般解が求められました $$y=\displaystyle \frac{Ce^{ax}}{1+C\frac{b}{a}e^{ax}}\tag{・・・式１８}$$