コンテンツにスキップ

変数分離形

対数積分では dyy=logy+C\int \frac{dy}{y}=\log|y|+C です。ただし指数化後の正負は任意定数 AA に吸収できます。

個体数 yy、成長率 r>0r>0、環境収容力 K>0K>0 に対して y=ry(1yK)y'=ry\left(1-\frac yK\right) を考えます。定数解は y=0,Ky=0,K。それ以外では部分分数分解により Ky(Ky)=1y+1Ky\frac{K}{y(K-y)}=\frac1y+\frac1{K-y} だから

logyKy=rx+C,y=K1+Aerx.\log\left|\frac{y}{K-y}\right|=rx+C,\qquad y=\frac{K}{1+Ae^{-rx}}.

初期値 y(0)=y0y(0)=y_0y00y_0\ne0)なら A=Ky0y0A=\frac{K-y_0}{y_0} です。分母が0にならない、初期点を含む区間が存在区間です。特に 0<y0<K0<y_0<K なら全実数で定義され、xx\to\inftyKK へ近づきます。

積分後に F(y)=G(x)+CF(y)=G(x)+C のような陰関数形しか得られないこともあります。それでも微分して元の方程式へ戻せれば立派な解の表現です。

Q基本:定数解★☆☆
y=x(y2)y'=x(y-2) の定数解を求め、非定数解の一般形も求めてください。
HINTヒントを見る
y2y-2 で割る前に確認します。
ANS解答を見る
定数解は y=2y=2。非定数解は dyy2=xdx\frac{dy}{y-2}=x\,dx より logy2=x22+C\log|y-2|=\frac{x^2}{2}+C、したがって y=2+Aex22y=2+Ae^{\frac{x^2}{2}} です。A=0A=0 とすれば定数解も含められます。
Q標準:ロジスティック方程式★★☆
y=y(1y),y(0)=12y'=y(1-y), y(0)=\frac{1}{2} を解き、長時間挙動を答えてください。
HINTヒントを見る
K=r=1K=r=1 の公式を使えます。
ANS解答を見る
A=11212=1A=\frac{1-\frac{1}{2}}{\frac{1}{2}}=1 より y=11+exy=\frac{1}{1+e^{-x}}。代入すると y=ex(1+ex)2=y(1y)y'=\frac{e^{-x}}{(1+e^{-x})^2}=y(1-y)。全実数上で定義され、xx\to\infty で1へ近づきます。
Q発展:有限時間発散★★★
y=2xy2,y(0)=1y'=2xy^2, y(0)=1 を解き、最大存在区間を求めてください。
HINTヒントを見る
y2dy=2xdxy^{-2}dy=2x\,dx
ANS解答を見る
1y=x2+C-\frac{1}{y}=x^2+C、初期条件から C=1C=-1。よって y=11x2y=\frac{1}{1-x^2}。微分すると 2x(1x2)2=2xy2\frac{2x}{(1-x^2)^2}=2xy^2。特異点は x=±1x=\pm1 なので、初期点0を含む最大存在区間は (1,1)(-1,1) です。