人口増加とロジスティック曲線

　今日は軽い話題を一つ述べましょう。

　全世界,または比較的出入りの少ない閉じた地域の現在の人口をN 人とし,Δｔ 年間に人口 N 人に比例して(ｋΔｔ)Ｎ人だけ人口が増加するとします。

　今の時刻（年）を ｔ として,年間の人口増加率ｋが ｋ ＝ｂ－ｄ （ｂは出生率,ｄは死亡率）で与えられる単純なモデルを想定するわけです。

　ｋは一定であると仮定すると,(ｔ＋Δｔ)年の人口:N₁＝ N＋ΔNはΔN /Ｎ＝ｋΔｔにより,Ｎ₁＝Ｎ (１＋ｋΔｔ）となりますから,結局,時刻（ｔ＋nΔｔ)年での人口N_n人はN_n＝N（１＋ｋΔｔ)ⁿになると予測されます。

　ｋ＞0 であれば,正に人口はネズミ算的に増えてゆきますね。

　Δｔ が無限小:ｄｔであるとすると,ΔN /Ｎ＝ｋΔｔはｄＮ/ｄｔ＝ｋＮ となります。

　この微分方程式を解けば,ｔ＝0 での人口をＮ(ｔ＝0)＝Ｎ₀人として時刻:ｔにおける人口はＮ＝Ｎ(ｔ)＝N₀ｅｘｐ(ｋｔ )で与えられるということになります。

　これを見ると,ｋ＞0 なら ｔ → ∞ ではＮ → ∞ ですが, 逆にｋ＜0 なら ｔ → ∞ で Ｎ → 0  なのでやがて絶滅してしまいます。

　しかし,実際には Δｔ の間にはいろいろな災害や環境の変化などあって,人口増加率 ｋ は一定ではなくかなりの変化を受けると考えられます。

　一般に人間をも含む生物個体の増加は個体総数Ｎが増えれば増えるほど妨げられる傾向がありますから,それは増加率がｋ＝(一定)から ｋ（１－αＮ）(α＞0 )となるような効果で表わすことができます。

　このモデルをロジスティックモデル(logistic model)といいます。

　増加率:ｋ（１－αＮ）(α＞0 )においてαＮ＞１なら人口(個体総数)Ｎは増加し,逆ならＮは減少しますね。

　このモデルは,Ｎに対する微分方程式の形ではｄＮ/ｄｔ＝ｋＮ（１－αＮ）という非線型微分方程式になります。

　具体的にこれを解くと,ｄＮ/Ｎ（１－αＮ）＝αｄＮ/(αＮ)＋αｄＮ/(1－αＮ)}＝ｋｄｔより,ln(αN)－ln(1－αN)＝ｋｔ＋CですからαN/(１－αN)＝Aexp(ｋｔ),つまりαN＝Aexp(ｋｔ)/{1＋Aexp(ｋｔ)},またはN＝Aα^-1/{A＋exp(－ｋｔ)}です。

　そしてαＮ₀＝A/(1＋A)によりA＝αＮ₀/(1－αＮ₀)ですから,N＝Ｎ₀/{αＮ₀＋(1－αＮ₀)exp(－ｋｔ)},結局,Ｎ＝Ｎ(ｔ)＝Ｎ₀/[αＮ₀{1－ｅｘｐ(- ｋｔ)}＋exp(－ｋｔ)]が得られます。

　あるいは,N(ｔ)＝(1/α)/[1＋{1/(αＮ₀)－1}exp(－ｋｔ)]です。

　これの描く(Ｎ-ｔ)曲線をロジスティック曲線(logistic-curve)と呼びます。

　これを見ると, ｔ → ∞ の極限ではＮ → １/αとなって,人口(個体総数)Ｎは,あ

る一定の極限値に到達します。それ以上は増加も減少もしません。

　ロジスティックモデルは実際に生態学(ecology)において個体の増加減少の履歴と一致する例が多々あり,人口にもこれが適用できると考えられます。

　これは,正に「増え過ぎた生物は抑制される。」という自然の摂理(神の摂理)を体現するモデルになっています。

　人類は天敵がいないことや医学の進歩,そして軍縮などによる戦争の減少？etc.によってこの摂理を破壊し,結果的に生態系を破壊しつつあります。

　やがては,この神の摂理の破壊の報いを受けるかも知れません。

　ところで,ロジスティック微分方程式のｄｔ＝Δｔ の刻みを調節して中心差分の差分方程式として離散化すると ,ｋ の値によってはｔ が大きいところで不安定な人口増減の振動をするカオス現象を起こすことが知られています。

　この不安定性は数値解析の目的で「離散化＝差分化」を行なったために生じたものですが,現実の現象のモデルとしては時間刻みが無限小の微分方程式よりも時間刻み有限の差分方程式の方が適切かも知れません。

　カオス(chaos)の例としては,ロジスティック模型:ｘ_n+1＝ａｘ_n(1－ｘ_n)は典型的なものであり,上式で時間刻みに相当するａの値によっては「リー・ヨーク(Li-Yorke)の定理」でのカオスになるべき条件を満たします。

※(注):ロジスティック差分方程式:N＋ΔＮ＝Ｎ＋ｋNΔｔ（１－αＮ）＝Ｎ(1＋ｋΔｔ)（１－αＮ）において,ｘ_n＝αN, ｘ_n+1＝α(N＋ΔＮ)と置いてパラメータをａ＝1＋ｋΔｔとして簡単化すれば,ｘ_n+1＝ａｘ_n(1－ｘ_n)に帰着します。※

参考文献；山口昌也 著「カオス入門」（朝倉書店）,山口昌也 編著「数値解析と非線型現象」（日本評論社）

http://fphys.nifty.com/(ニフティ「物理フォーラム」サブマネージャー）                                       TOSHI

http://blog.with2.net/link.php?269343(ブログ・ランキングの投票）

↑ここをクリックして投票をお願いします。

ここをクリックして投票をお願いします。