QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 9516|回复: 0
打印 上一主题 下一主题

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

[复制链接]
字体大小: 正常 放大

600

主题

29

听众

6868

积分

  • TA的每日心情
    奋斗
    2023-5-24 09:14
  • 签到天数: 119 天

    [LV.6]常住居民II

    群组2018高中组美赛 课堂

    群组2018国赛冲刺

    群组2018 夏令营面授课堂

    群组2016美赛交流群组

    跳转到指定楼层
    1#
    发表于 2018-11-2 08:37 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    数学建模笔记之马尔萨斯人口模型(一)
    ( t2 K0 r7 Y- x3 Y[size=0.8]本文使用Latex公式排版,而简书不支持Latex排版,为获得更好的阅读体验,请移步个人博客原文地址:数学建模笔记之马尔萨斯人口模型(一)
    1 V" O$ w- j& q' \+ A5 w) D1 g[size=0.8]人口按照几何增长趋势发展( 按照指数函数增长的趋势 ) 而实物只有算术增长的趋势( 按照线性函数增长的趋势 ) 结论:控制人口增长
    7 n, y) W) e" l* X* p, j[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$趋近于无穷大的时候我们发现人口确实成指数增长) w0 m, S/ K& s
    [size=0.8]从这里我们可以看到数学建模研究的一个大致思路:
    7 _. }+ k& O4 j0 r, |" }[size=0.8]数模思路
    : C# `( M, P4 W: v3 x[size=0.8]数模思路* L  K) U- U2 |: H2 j
    [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$趋近于无穷大时的极限,结果与马尔萨斯模型完全不同 我们可以绘制出人口与时间关系的图形:
    " C" S, F9 E6 |6 s[size=0.8]logistic模型图. p' ~! F2 K. b6 X& h1 J- d) N
    [size=0.8]logistic模型图% [* Q0 X! g$ P! M
    [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$: 有:+ J, M3 E$ H* o' r/ e; _! _
    [size=0.8]一周期解1 p0 S9 Q- m7 H$ b, N/ [3 A
    [size=0.8]一周期解; \7 }7 Q; p! a" w5 m# u8 T$ Q0 q
    [size=0.8]![二周期解] ](http://upload-images.jianshu.io/upload_images/3020761-2ad6a7cd778385cc.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/1240)
    5 G4 L. X$ Q" q  t: B, a# z[size=0.8]四周期的解9 _, ]# l3 Q7 {
    [size=0.8]四周期的解
    8 G1 }9 R, d8 @9 B[size=0.8]八周期的解- }1 D0 _/ O6 C! ^; F- }9 ?( ?8 A
    [size=0.8]八周期的解" a' c) h# v8 P( R- @2 z
    [size=0.8]倍周期示例5 m( w; [. C" O; d1 E8 u8 I* @9 r
    [size=0.8]倍周期示例: R, c& i  N. x& K2 T$ u0 X' c
    [size=0.8]混沌
    5 I; G- ^0 A8 G& i" ]7 N3 o: F[size=0.8]混沌
    * g9 ^3 B0 O) t# @' ~[size=0.8]后面的关于时间$t$的离散讨论已经超越这个模型本身了,仅作简单介绍5 M: q% }" p( T6 ], @0 Z; ^

    $ v7 J% _, A9 a) L( x2 N" D
    ( d; ~. G/ N1 Y. C; D8 ~, I) {
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-12 12:43 , Processed in 0.272122 second(s), 49 queries .

    回顶部