数学建模笔记之马尔萨斯人口模型(一)

重光兰衣

数学建模笔记之马尔萨斯人口模型(一)
% C! b5 S+ s: ^* ~9 e[size=0.8]本文使用Latex公式排版,而简书不支持Latex排版,为获得更好的阅读体验,请移步个人博客原文地址:数学建模笔记之马尔萨斯人口模型(一)
[size=0.8]人口按照几何增长趋势发展( 按照指数函数增长的趋势 ) 而实物只有算术增长的趋势( 按照线性函数增长的趋势 ) 结论:控制人口增长
[size=0.8]现在我们使用 $P(t)$ 表示t时刻的人口数量,用$r$表示人口增长率 现在看做一个连续模型:变化在随时发生,也即人的生老病死随时在发生, 则有: $$P(t+\Delta t)-P(t)=rP(t)\Delta t$$ $$P(t+dt)-P(t)=rP(t)dt$$ $$\dfrac{dP(t)}{dt}=rP(t)$$ $$P(t_0)=P_0$$ 则我们可以根据以上公式推断出当t为某个具体值时, 我们得到的人口数目为 $P(t)=P_0e^{r(t-t_0)}$, 当$t$趋近于无穷大的时候我们发现人口确实成指数增长
[size=0.8]从这里我们可以看到数学建模研究的一个大致思路:
( _6 V  g2 d; T6 r4 s4 S$ r' p[size=0.8]数模思路
! L5 K0 j7 Z" R- c) C1 X7 M/ s[size=0.8]数模思路
2 ?' q5 V  f0 W# W0 f[size=0.8]在上述马尔萨斯模型中, 人口增长率是基于一定时间段的结果, 例如基于当时的工业农业现有人口现状是有一定意义的,即是说当时的情况下人口增长率是可以在在短期内保持一个特定的数值 但是在当今情况下, 人口增长率必定是一个关于时间的函数 $$r(t)=r(P(t))=r(1-\dfrac{P(t)}{K})$$这里的K表示我们所研究的生态系统最多可以容纳支撑的人口数量(即生物学上所说的最大容纳量) 则人口数量关于时间的积分有: $$\frac{dN(t)}{dt}=r(1-\frac{N(t)}{K})N(t)$$ $$N(t_0)=N_0$$ 则我们取积分后可以得到: $$N(t)=\frac{K}{1+Ce^{-r(t-t_0)}}$$ 其中C表示一个参数: $$C=\frac{K-P_0}{P_0}$$ 则我们考察$t$趋近于无穷大时的极限,结果与马尔萨斯模型完全不同 我们可以绘制出人口与时间关系的图形:
[size=0.8]logistic模型图
+ }) T0 R, H' K9 C2 e6 \[size=0.8]logistic模型图
[size=0.8]与马尔萨斯模型相比较,这里增加了一个参数$K$, 而且这个参数不容易计算或估计($K$表示环境最大容纳量) 则我们对模型 $\frac{dN}{dt}=rN(1-\frac{N}{K})$ 考虑离散化: $$\frac{\Delta N}{\Delta t}=rN(1-\frac{N}{K})$$ 则有: $$N_{t+1}-N_t=rN_t(1-\frac{N_t}{K})$$ 这里的时间离散步为1,每一代就是一个时间步 $$N_{t+1}=(1+r)N_t-\frac{r}{K}N_t^2$$ 然后我们取定参数$K$,考虑不同的参数$r$: 有:
' G2 P0 U+ v8 W' B9 u9 {7 Y& u[size=0.8]一周期解
7 Y4 _3 n: b+ p- v- ?[size=0.8]一周期解
7 s1 S: C" T9 E; R5 n[size=0.8]![二周期解] ](http://upload-images.jianshu.io/upload_images/3020761-2ad6a7cd778385cc.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/1240)
9 ?% w9 f9 V% O& a8 \[size=0.8]四周期的解
[size=0.8]四周期的解
: D% O- m1 X" Q( v& w" z" Y. r9 B[size=0.8]八周期的解
[size=0.8]八周期的解
0 k9 x- |8 s  Q' P3 p  F* _* E[size=0.8]倍周期示例
5 r& Z  ]; t; m[size=0.8]倍周期示例
[size=0.8]混沌
[size=0.8]混沌
1 n" T( a1 K$ i$ F' ~" d. K- p[size=0.8]后面的关于时间$t$的离散讨论已经超越这个模型本身了,仅作简单介绍
4 ]: g# M' k- Y, Y" D
& l0 i8 S9 y0 K
9 C" h  B+ q) s