- 1、本文档共29页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
Leslie模型数学建模
提出问题: 我们首先考虑Malthus 模型: x(t)为人口总数,r为自然增长率; 于是可以得出: x(t)=x0er(t-t0) 改进的模型 设地球能容纳的总人数为k,随着人口的增长,出生率必然会下降,于是r与x存在着一定的关系。基于上述假设,我们选择一种简单的函数。 r(x)=r0(1-x/k) r0为特定的常数 定义 一些定义: n为人类的年龄上限 F(x)=x岁的妇女所生的婴儿数/x岁的总人口数 S(x)=x岁人的存活率 P(x)=初始时x岁的总人口数 Nt(x)=距离初始t年时x岁的总人口数 数学表达: 建立模型: 一些定义: n为人类的年龄上限 a(x)=x岁的妇女所生的婴数/x岁的妇女总数 b(x)=x岁人的存活率 h(x)= x岁的妇女迁移数/x岁的妇女总数 Nt(x)=距离初始t年时x岁的总人口数 建立模型: * * 关于建立人口增长模型,我们考虑了两条主要思路: 一.以微分方程为主要手段: 二.以高等代数为主要手段: 解得: x(t)=k/[1+(k/x0-1)e-r(t-t0)] 分析以上两个模型: 每个个体的出生率与死亡率是相同的。但实际上不同年龄的年的生育率与死亡率有很大的不同。 基于这种考虑,下面将建立一个人口按年龄分布的模型 r表示年龄,函数F(r,t)为t时刻年龄小于r的人口总数,称其为人口分布函数 令p(r,t)= F/ r p(r,t)为年龄密度函数 则t时刻年龄处在[r,r+dr)的人口总数为p(r,t)dr 设μ(r,t)为t时刻年龄为r的人的死亡率,t时刻年龄在[r,r+dr)单位时间死亡的人数为μ(r,t)p(r,t)dr 分析: 下面考虑从t到t+dt这一过程的人口变化: 年龄处在[r,r+dr)到t+dt时刻活着的人的年龄变为[r+dt,r+dr+dt)而这一时刻死亡的人数为μ(r,t)p(r,t)drdt 则p(r,t)dr-p(r+dt,t+dt)dr= μ(r,t)p(r,t)drdt p/ r + p/ t=- μ(r,t)p(r,t) p(r,0)=p0(r) p(0,t)=f(t) p0(r-t)e f(t-r)e 在社会比较安定的情况下,死亡率大致与时间无关. μ(r,t)=μ(r) p(r,t)= 0≤t≤r tr 分析: 1.当tr时,p(r,t)完全由年龄为r-t的人口的初始密度及这些人的死亡率决定。 2.tr 时,p(r,t)完全由未来的生育状况f(t-r)及死亡率决定。 两个重要模型: Keyfitz Leslie K= …… P(0) P(1) P(n) I(t)= …… Nt(0) Nt(1) Nt(n) 第一年新生儿的总数: F(0)?P(0)+ F(1)?P(1)+ ??? +F(n)?P(n) 第一年x岁人口总数: N1(x)=S(x-1)?P(x-1) 第一年末人口总数:F(0)?P(0)+ F(1)?P(1)+ ??? +F(n)?P(n)+ S(0)?P(0)+ S(1)?P(1)+ ??? +S(n-1)?P(n-1) 构造n+1阶方阵 M= F(0) F(1) F(2) ? ? ? F(n) S(0) S(1) S(2) ? ? ? S(n-1) 那么I (1)=MK I (t)=MtK 考虑到在一段稳定的时间段内:总的女性人口数比上总的男性人口数为一个近似为1的定值.为了更为确切地分析女性个体数量的分布对总人口数的影响,我们单独把女性人口数作为研究对象. 另外在这个模型中我们还加上了人口迁移对起其总数的影响. K= I(0) I(t)= …… Nt(0) Nt(1) Nt(n) H(t)= h(0) h(1) ….. h(n) 数学表达: 第一年新生女婴的总数: a(0)?Nt(0)+ a(1)?Nt(1)+ ??? +a(n)? Nt(n) 第一年x岁女性人口总数: N1(x)=b(x-1)?Nt(x-1)- h(x-1)?b(x-1)?Nt(x-1)=(1- h(x-1) )?b(x-1)?Nt(x-1) 第一年末
文档评论(0)