关于偏微分方程人口模型中死亡率的确定
用偏微分方程求解人口分布问题

用偏微分方程求解人口分布问题陈喜林【摘要】随着物理科学的发展,偏微分方程的应用范围越来越广泛.而人口问题始终是制约我国发展的关键因素,预测我国人口的发展趋势,是我国人口发展战略的核心内容之一.应用一维迁移方程导出人口问题的数学模型,并得出偏微分方程的解.【期刊名称】《佳木斯大学学报(自然科学版)》【年(卷),期】2016(034)006【总页数】3页(P990-992)【关键词】一维迁移方程;一阶线性偏微分方程;人口【作者】陈喜林【作者单位】罗定职业技术学院教育系,广东罗东527200【正文语种】中文【中图分类】O175设人口的数量不仅与时间t有关,还和年龄x有关,作为一个连续模型(设人口总数很大),以u(x,t)表示人口在任意时刻x年龄段t岁的分布密度,其中t≥0,x∈[0,A],A为人的最大寿命。
这里所说的人口,并不一定限于人,可以是任何一种生物群体,只要有类似的性质,也可以应用一维迁移方程导出人口问题的数学模型。
考虑到死亡对人口分布的影响,假设某年龄人口的死亡数与该年龄人口数成正比,其比例系数为d=d(x),称为死亡率。
由于假定社会是稳定的,它与时间无关,以人口数作为迁移量,则其源密度为-d·u(x,t),于是便可导出人口分布密度u(x,t)所满足的偏微分方程。
为什么0岁人口的出生未考虑在源密度内呢?因为在任何时刻t,所有出生的婴儿都是0岁人口,其总数应是u(t,0),这便决定了x=0处的边值条件,设b(x)为出生率,它仅与年龄有关,而与时间t无关,记a>0为最小生育年龄,则应有它便是一个边值条件,然而它并非我们学过的三类边值条件(Dirichletn条件,Neumann条件,Robin条件)中的任何一个,这是一个非局部边值条件。
那么如何导出所满足的偏微分方程呢?为此将问题简化成:导出人口分布密度u(x,t),t≥0,x∈[0,A],A为人的最大寿命,d(x)为死亡率,b(x)为出生率,a<A为最小生育年龄,设u(x,0)=u0(x),问u0(x)应满足怎样与边值条件相容的条件?在解决物理学、力学及工程技术等实际问题的过程中,有许多被考察的量不停地在空间进行传递和迁移,如流体的质量、能量、物体的热量等等,一般地,统称这些量为被迁移量,被迁移量与其所在的空间的位置与时间是有关的,引入被迁移量密度φ(x,t),x∈R3为时刻t,点x处的被迁移量,取包含x点的区域上在时刻t的被迁移量的总和除以该区域之体积为在时刻t区域上的平均被迁移量。
偏微分方程在数学建模中的应用分析——评《偏微分方程》

书评•广告偏微分方程在数学建模中的应用分析—评《偏微分方程》蒋楠(上海工程技术大学高等职业技术学院,上海200439)偏微分方程是高校数学应用专业的专业课程,是大学理工科学生的必选科目之一,其学习的主要内容是数学方程在物理、化学、经济、生物等诸多领域中建模的实际应用。
由于该课程内容复杂,概念抽象,导致学生在学习过程中不易消化理解,长期以来这成为了高校工科大学生学习的一个难题。
偏微分方程作为数学理论与实际应用有机衔接的纽带,对于当前我国科技创新背景下的社会经济发展非常重要,对于探索我国工程教育领跑世界新模式、总结中国发展经验、促进创新驱动战略、助力科技教育强国具有重要意义。
如何创新偏微分方程课程的内容和教学形式,激发学生的学习兴趣,是当前基于高校教学改革的偏微分方程教学改革的重要研究课题之一。
由孔德兴编著、高等教育出版社于2010年出版的《偏微分方程》一书在这方面做了积极有益的探索。
第一,该书作者根据自己长期从事偏微分方程教学的实践经验,结合自己的科研工作,以全新的理论角度和实践创新经验向广大读者详细地介绍了偏微分方程的理论基础知识,阐述了当代数学有关偏微分方程的一些新概念,并结合具体事例论述了偏微分方程在工科领域中的应用。
偏微分方程是现代数学的重要分支,人们在偏微分方程的研究过程中取得了丰硕的成果,并将其应用到各个领域,使其成为解决这些领域中实际难题的重要工具。
第二,该书在向广大读者诠释传统微积分方程基础理论知识的过程中,还基于数学理论的系统性向读者详细地介绍了当代与偏微分方程有关的其他重要的数学概念,这不仅有利于读者学习和理解,更有利于相关人员进行应用实践。
另外,亥书作者为了把传统偏微分方程的基础理论知识讲得更通俗易懂,还结合偏微分方程建模应用系统地总结了当代与偏微分方程相关的一些新的数学基础知识,展示了有关方面的最新成果,这有利于读者融会贯通,进一步提升对偏微分方程的理解与应用,也增强了数学理论的趣味性,使得对偏微分方程的学习不再枯燥而抽象,能够有效提升学习者的学习兴趣。
-偏微分方程模型

交通管理部门
尽可能多的人安全地通过
集中参数法:
假设车流量是均匀分布 目标使车流密度保持在安全的范围之内,让司机尽 可能开得快些即可,必要时司机自己会刹车。
现实生活中可能吗?
车流密度和车速不可能是常数
分布参数法:
x轴表示公路,x轴正向表示车流方向。
如果采用连续模型,设u(t,x)为时刻t时车辆按x方向分布 的密度,再设q(t,x)为车辆通过x点的流通率。
初值条件:
h(0, x) u f
f uj
u0 (x)
将Greenshields的基本方程代入(3.41),利用复合函数求导法则
并注意到uf、uj均为常数,可得:
u t
(t, x) (u f
2u u
f j
u
)
u x
(t,
x)
0
令
h uf
2u f uj
u ,方程可简化为: h (t, x) h h (t, x) 0
t
x
令 p(t, x) 为 t 时刻年龄为 x 的人口密度,则 t 时人口总数为:
A
P(t) 0 P(t, x)dx
其中A为人的最大寿命。
设t时刻年龄为x的人的死亡率为d(t,x),则有:
p(t dt, x)dx p(t, x dt)dx d (t, x dt) p(t, x)dxdt
P(t )
P(t )
此即Malthus模型
B(t)、D(t)分别为t时刻的生育率和死亡率。则有:
dP (B(t) D(t))P(t) dt 若B(t)、D(t)与t无关,则可得:
•
dP
(B
数学建模讲座机理分析方法及例子1

不稳定,轨道{xn}趋向稳定点
■ 当3<a<1+61/2时, xn 绕着两个数 x3*,x4*振动,
例 a =3.2
x2k-1 →0.799455
x2k →o.513045
这两个数满足
x f 2 ( x), x f ( x)
也称为周期2点,对应轨道称周期2轨道.(原来周期
n = 0,1,2,…
● 数值迭代( a 逐渐增加,迭代会有何结果)
1.倍周期分叉现象
■ 当0<a <1时,由于0<xn<axn+1
xn →0
物种逐渐灭亡
■ 当1<a<3时,任何(0,1)中初始值的轨道趋于
x*=1-1/a 其中x*是方程f(x)=x的解,为映射f 的不动点
(周期1点)例:a =1.5时 xn → 1/3.
~总和生育率
f
(t )
(t) r2 r1
h(r , t )k
(r,
t)
p(r , t )dr
人口发展方程和生育率
f
(t)
(t) r2 r1
h(r , t )k
(r,t)
p(r,
t)dr
(t) ~总和生育率——控制生育的多少
h(r, t ) ~生育模式——控制生育的早晚和疏密
p(r,t)
p0
约35年增加一倍,与1700-1961年世界人 统口计结果一致
与近年统计结果有误差,由a >1,xn趋向无穷, 模型在人口长期预测方面必定是失效的.
● Logistic模型
.
生存资源是重要的因素,修改模型为:
xn+1 - xn= r xn- b xn2 - b xn2为竞争(约束)项,r、b 称生命系数,则
4.3偏微分方程模型

§4.3 偏微分方程模型如果一个微分方程中出现多元函数的偏导数,或者说如果未知函数和几个变量有关,而且方程中出现未知函数对几个变量的导数,那么这种微分方程就是偏微分方程。
本节以人口增长模型和扩散模型为例说明偏微分方程的建模过程以及相应的数值解法。
4.3.1 人口增长模型统计数据表明,世界人口在1800年达到10亿,1930年达到20亿,1960年达到30亿,1974年达到40亿,1987年达到50亿,1999年达到60亿,2011年10月31日突破70亿。
可以看出,人口每增加10亿的时间由100多年缩短为10余年。
人口的剧增导致资源消费量增加,引起资源蓄积量减少甚至枯竭,出现诸如过度开垦土地、沙漠化日益严重、不合理地砍伐森林、绿色空间缩小、能源紧张等问题。
人口剧增还会带来空气污染,引起全球气候变化异常等环境问题,造成全球性生态平衡失调。
而且,这么多数量的人口空间分布极其不均衡。
全球45个发达国家的生育率都低于人口平均增长率。
在世界出生率最低的25个国家中,有22个在欧洲。
人口数量的减少成为这些国家最大的危机,对经济发展和国家安全带来严峻挑战。
同时,世界上人口增长率最高的都是一些最不发达的国家,如阿富汗、布隆迪、刚果、利比里亚等,而发展速度较快的发展中国家,如中国、印度、埃及等,也身负人口增加给经济和环境带来的巨大压力。
中国是世界上人口最多的国家,根据2010年第六次人口普查登记的全国总人口为13.3972亿(不包括港澳台地区),其中,男性人口6.8685亿,女性5287亿;60岁以上人口为1.7765亿,占总人口的13.26%;城市人口为6.6558亿,农村人口为6.7415亿。
老龄化问题、男女比例失调、城镇化建设加速等问题成为我国人口问题的一些新特点,直接影响着我国人口的发展趋势[1]。
准确地对人口进行预测,有效地控制人口增长并制定合理的人口政策,是全面落实科学发展观、实现适当生育水平、提高人口素质、改善人口结构、引导人口合理分布、保障人口安全、促进人口与经济社会资源环境的协调和可持续发展的重要手段。
第二章微分方程模型-22人口问题模型

*** 模型修改
分析表明,以上这些现象的主要原因是 随着人口的增长,自然资源,环境条件等因 素对人口增长的限制作用越来越显著。人口 较少时,人口的自然增长率基本上是常数, 而当人口增加到一定数量以后,这个增长率 就要随着人口的增加而减少。因此,我们将 对指数模型关于净相对增长率是常数的基本 假设进行修改。
我国是世界第一人口大国,地球上每九 个人中就有二个中国人,在20世纪的一段 时间内我国人口的增长速度过快,如下表:
年
1908 1933 1953 1964 1982 1990 2000
人口(亿)3.0 4.7 6.0 7.2 10.3 11.3 12.95
2020年,13.4亿,2013年,13.5亿。有效地 控制人口的增长,不仅是使我国全面进入小 康社会、到21世纪中叶建成富强民主文明的 社会主义国家的需要,而且对于全人类社会 的美好理想来说,也是我们义不容辞的责任。
返回
*** 模型求解 分离变量得
dN
N
1
N Nm
=rdt 即
N
NmdN
Nm
N
=rdt
或
dN dN
N + Nm N =rdt 两边积分得
ln N -ln Nm-N =C1ert
代入边界条件,并整理得
N t=
Nm
1+
模型分析 容易看出
t 时,N t Nm.
陆地,每人只有9.3平方英尺的活动范围。
而到2670年,人口达到 361015亿,只有一个人
站到另一个的肩上了。 因此,Malthus人口模型是不完善的。从根本上 说是不完整的,必须修正。
问题在于:在上述模型中假设r是常数,从 而人口方程是线性常微分方程。这个模型 在群体总数不太大时才合理。而没有考虑 总数增大时,生物群体的各成员之间由于 有限的生存空间,有限的自然资源及食物 等原因,就要进行生存竞争。
数学建模习题中国人口增长预测

中国人口增长预测本题是一个人口发展预测的问题。
人口发展与一般种群增长一样,是由自然增长率决定的。
然而,人类个体是一种社会的个体,所以人口发展有自己的特点。
想到人口的迁移,性别比例,城镇化等。
同时,人口发展受政策的影响,例如计划生育;也要受到人们意识的影响,像生育意识等。
但是从社会层面上看,生育意识在整个社会上体现为妇女的生育模式,进而可以特别地去考虑。
思考方法:首先,数据的处理。
在经过EXCEL分析和验证后,适当修正题中的个别有误数据后,利用有效数据进行建模求解,在此过程中,我们提取出死亡率、生育率等感念,且把人的一生按年龄分为青年期、衰老期等阶段。
这是求解人口增长模型的必要过程和方法。
其次,模型建立。
和一般的预测模型一样,本模型也是个预测模型,所以考虑到用题目所给的五年的信息,来推测今后几十年的人口的总数和结构情况。
对此,我们选用差分方程模型和数据参数拟合等方法。
同时,将死亡率与出生率分开分别计算和拟合,通过五年的实际数据拟合出相应函数的参数,再利用此函数进行评估和预测。
最后,利用已有信息以及上述所求出的对应函数和方程,对中短期与长期进行估计和预测,进而得出人口结构、人口比例、人口数量等一系列的相关数据。
以下是解答过程:1.数据说明:x:表示最大的年龄;mi=1,2,3,4,5,6 其中1表示市男性,2表示市女性,3表示镇男性,4表示镇女性,5表示乡男性,6表示乡女性;A :表示婴儿性别比例矩阵;* :表示点乘;P(x,t):表示t时刻年龄为x的人口数量;ibir(x,t):表示t时刻年龄为x的出生率;i)(,i dea x t:表示t时刻年龄为x的死亡率;)(i t k:表示t时刻婴儿的死亡率;tra(x,t):表示t时刻年龄为x的人口迁出率;i2.假设条件1. 假设国内社会环境稳定,无异常大量死亡或出生情况发生,人口比例,人口总数不会出现突变状况; 2. 假设只存在乡向城镇迁出,不存在其他迁移方式,且不同年龄段迁移率相同; 3. 假设不考虑国家之间的迁入与迁出,把中国内部看为一个封闭的模型; 4. 对于90岁以上的人都按照90岁处理; 5. 假设只存在乡向城镇迁出,不存在其他迁移方式,且不同年龄段迁出率相同,按照0.6%均匀增长。
Logistic人口模型的检验与中国未来人口的预测与控制

数学建模大作业——人口模型班级:周一、周三姓名:石星宇学号:02123010学院:电子工程学院任课教师:李伟Logistic人口模型的检验与中国未来人口的预测与控制摘要本文利用已有的Logistic人口模型,对中国近年来的人口数据进行了分析。
用MATLAB数学软件对人口数据做曲线拟合,得到三次多项式如下:33269=-+⨯-⨯+⨯.y x x x0.7011 4.1738108.280910 5.475510±,可见该模型根据实际人口数和计算人口数的比较得知,最大的误差为0.5能够很好地预测人口发展趋势,由此验证了Logistic人口模型的正确性。
并用得到的三次多项式对今年年末的人口进行预测,计算结果为132850万人。
本文最后还对当前国内的人口政策做出了一定的讨论,论述了控制人口发张的一些方法及策略。
例如,提倡一对夫妇只生一个孩子、晚婚晚育等。
还对近年来提出的“单独二胎”、“双独二胎”政策做了相应的讨论。
得出结论:这些政策的实施将会导致婴儿出生率短期内显著回升,有望在5年内新增750万新生儿,在人口结构方面,政策调整将使2030年增加2200万劳动年龄人口。
关键词:Logistic人口模型曲线拟合人口预测与控制1、问题重述人口问题是当今世界上最令人关注的问题之一。
一些发展中国家的人口出生率过高,越来越严重的威胁着人类的正常生活,而有些发达国家的自然生长率趋于零,甚至变为负数,造成劳动力短缺,也是不容忽视的问题。
由于我国20世纪50~60年代人口政策方面的失误,不仅造成人口总数增长过快,而且年龄结构也不合理,使得对人口增长的严格控制会导致人口老化问题严重。
因此,在首先保证人口有限增长的前提下适当控制人口老化,把年龄结构调整到合适的水平,是一项长期而又艰巨的任务。
我们目前面对的问题有:(1)检验Logistic人口模型的正确性;(2)预测中国未来人口的发展状况;(3)评价中国现有的人口政策。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
- 〔exp ( - ∫0xλN (s) f2 (t - x + s) ds) 〕dx + ∫LtλN (x) φ珔 (x - t) ·
π (x) 〔exp ( - ∫xx - λt N (s) f1 (t - x + s) ds)
- exp ( - ∫xx - λt N (s) f2 (t - x + s) ds) 〕dx
ห้องสมุดไป่ตู้
≤D (t ; 0)
从而就有 :
B
(t) - P′(t) D (t ; 0)
-
1 ≤T 〔f〕 ≤B
(t) - P′(t) D (t ; M)
-
1
再由假设条件 A4 可知 : 0 ≤T 〔f〕 ≤M
因此 , T : CM →CM 记 : Λ (x) = ∫0xλN (s) ds , π (x)
= exp ( - ∫0xλN (s) ds) , φ珔 (x)
3 问题的解决
定理 : 如 果 假 设 条 件 A1 —A4 成 立 , 则 反 问 题 ( 211) — ( 214) 在 CM 中 存 在 唯 一 解 f , 其 中 CM = {f ∈C (0 , t) | 0 ≤f (t) ≤M}
证明 : 利用特征线法可知正问题 (211) — (213) 的解具有如下形式 :
1998 - 06 - 08 收稿 第一作者 郝新生 (1968 - ) , 男 , 山西武乡人 , 讲师 , 主要从事偏微分方程反问题及应用方面的研究 。
2000 郝新生 : 关于偏微分方程人口模型中死亡率的确定
ut (x ,t) + ux (x ,t) = - λN (x)〔1 + f ( P(t) ) 〕u (x ,t) , 0 ≤x ≤L ,0 ≤t ≤T u (x ,0) = φ(x) , 0 ≤x ≤L u (0 ,t) = B (t) , 0 ≤t ≤T
- P′(t) D (t ; f2)
|
·| D
(t ;
f1)
- D (t ; f2) |
≤| B
(t) D2
- P′(t) (t ; M)
||D
(t ; f1)
- D (t ; f2) | ≤珔C ‖f1 - f2 ‖〔0 ,T〕
其中 , 珔C =珔C
(t , λN ,
B , φ,
M,
p)
,
且 lim珔C = 0 t →0
∫0tλN (x) π (x) Λ (x) dx = - ∫0tπ′(x) Λ (x) dx = 1 - π (t) - Λ (t) e - Λ(t)
=1,π
∫LtλN (x) π (x) 〔Λ (x) - Λ (x - t) 〕dx =π (t) Λ (t) + ∫Lt 〔λN (x) - λN (x - t) 〕π (x) dx 令 : C1 (t , λN) = 1 - π (t) - Λ (t) e - Λ(t)
001382 J . Shanxi Agric. Univ. Vol. 20 No. 1 山西农业大学学报 A 1000 - 162X (2000) 01 - 0027 - 03
关于偏微分方程人口模型中死亡率的确定
郝新生
(山西农业大学 基础部 , 山西 太谷 030801)
摘 要 : 本文利用不动点方法讨论了一类描述人口变化规律的一阶双曲型方程确定死亡率 d (x , t) 的一个反问题 , 在已知函数的一些假设条件下 , 给出该问题解的存在唯一性定理 。 关键词 : 人口模型 ; 一阶双曲型方程 ; 反问题 ; 死亡率 中图分类号 :) 175. 2
其中 , C (t , B , φ, λN) = C1 (t , λN) ‖B ‖〔O ,T〕+ C2 (t , λN) ‖φ珔‖〔O ,t〕
显然 ,
C
(t , B , φ, λN)
在 〔0 , L〕上有界 ,
且 limC = 0 t →0
| T 〔f1〕 -
T 〔f1〕|
≤|
D
B (t
(t) ; f1)
φ =π
(x) (x)
,
Πf1 , f2 ∈CM
D (t ; f1) - D (t ; f2) = ∫0xλNB (t - x) π (x) 〔exp ( - ∫0xλN (s) f1 (t - x + s) ds) 〕
90
山 西 农 业 大 学 学 报 20 (1)
∫x
B (t - x) exp{ - λN (s)〔1 + f (t - x + s) 〕ds} ,x < t
0
∫ u (x ,t ;f) =
x
φ(x - t) exp{ - λN (x - t + s)〔1 + f (s) 〕ds} ,x ≥t
0
(3. 1)
Π (x , t) ∈〔0 , L〕 ×〔0 , T〕, 显然有 :
Determing Death - rate for Population Model of Partial Differential Equation Hao Xinsheng
( Department of Basic , Shanxi Agricultural University , Taigu , Shanxi , 030801) Abstract : Using Fixed - Point Method , this paper deals with an inverse problem of first - order hyperbolic equation of describing developing law of Population , under some assumed conditions of known functions , the existence and uniqueness of the solution of this inverse probem are demonstrated. Key words : Population Model ; First - order hyperboic equation ; Inverse problem ; Death - rate
参考文献 [ 1 ] W. Rundell , Determing the birth function for an age structured population [J ] . Mathematical Population studies 1989 (4) : 377~395 [2 ] 谷越豪 , 李大潜 , 沈玮熙. 应用偏微分方程 [M] , 北京 : 高等教育出版社 , 1993 [3 ] 宋 健编. 人口预测和人口控制 [M] , 北京 : 人民教育出版社 , 1981 [4 ] 姜礼尚 , 陈亚渐 1 数理方程讲义 [M] , 北京 : 高等教育出版社 , 1986
- ∫L0λN (x) (1 + f) u (x , t ; f) dx = B (t) - (1 + f) D (t ; f)
因此 ,
f =B
(t) D
- p′(t) (t ; f)
-
1
作算子 T T 〔f〕 = B
(t) D
- P′(t) (t ; f)
-1
由 0 ≤f ≤M , 可得 D (t , M) ≤D (t ; f)
其中 , L 为人口的生命极限 , P (t) 为 t 时刻的人口总数 , 从而问题 (1. 1) — (1. 4) 就构成了
一阶双曲型方程确定未知系数的反问题 。
2 问题的转化和数据假设
我们认为人口的死亡率不仅与年龄有关 , 而且还与当时的人口总数有关 , 因此我们可以把死亡率 d (x , t) 表示为λN (x) 〔1 + f (p (t) ) 〕, 即 d (x , t) =λN (x) 〔1 + f (p (t) ) 〕, 其中 , λN (x) 表示人 口的自然死亡率 , 为一已知函数 , f (p (t) ) 表示人口环境因素造成的死亡率 , 为一未知函数 。这样问 题就转化为如下确定函数偶 (u , f) 的反问题 :
∫L
P(t) = u (x ,t) dx , 0 ≤t ≤T 0
89
(2. 1) (2. 2) (2. 3) (2. 4)
下面对已知数据作适当假设 : A11λN (x) 为连续正值函数 , 且 lxi→mLλN (x) = ∞ A21φ (x) 为 〔0 , L〕上的非负连续函数 。 A31 B (t) 为 〔0 , T〕上连续正值函数 。 A41 P (t) 为 〔0 , T〕上的连续可微的正值增函数 , 并且对 M > 0 满足 : B (t) - (M + 1) D (t ; 0) ≤p′(t) ≤B (t) - D (t ; M) 其中 , D (t ; f) = ∫L0λN (x) u (x , t ; f) dx
(3. 3)
对 (313) 式应用不等式 : | e - x - e - y| ≤| x - y| 可得 :
| D ( t ; f1 ) - D ( t ; f2 ) | ≤ ‖ B ‖〔O ,T〕 ∫0tλN ( x ) π ( x ) Λ ( x ) dx ·
‖f1 - f2 ‖〔O ,T〕+ ‖φ珔〔O ,L〕‖∫LtλN (x)π (x) 〔Λ (x) - Λ (x - t) 〕dx·‖f1 - f2 ‖〔O ,T〕利用π (0) (L) = 0 , Λ (0) = 0 , π′ (x) = 2π (x) λN (x) 可得 :
1 问题的提出
一个稳定社会的人口发展规律可用下面的一阶双曲线方程描述 :
ut (x , t) + ux (x , t) = - d (x , t) u (x , t)