《ch绪论人口模型与混沌实用.pptx》由会员分享，可在线阅读，更多相关《ch绪论人口模型与混沌实用.pptx（30页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 巴西的一只蝴蝶扇动翅膀会引起 洛伦兹 数学的伟大使命在于从混沌中发现秩序。倍尔在那个混沌的体制中，结构上的微小差异几乎都会造成行为方式上的巨大变化，可控制的行为似乎已被排除。斯图尔特.考夫曼明年在得克萨斯的大风暴吗？第1页/共30页背背景景 年 1625 1830 1930 1960 1974 1987 1999人口(亿)5 10 20 30 40 50 60世界人口增长概况中国人口增长概况 年 1908 1933 1953 1964 1982 1990 1995 2000人口(亿)3.0 4.7 6.0 7.2 10.3 11.3 12.0 13.0研究人口变化规律控制人口过快增长如何预报
2、人口的增长第2页/共30页x(t)-时刻t t的人口基本假设 :人口(相对)增长率 r 是常数今年人口 x0,年增长率 rk年后人口随着时间增加，人口按指数规律无限增长建立数建立数学模型学模型1、指数增长模型(Malthus模型）第3页/共30页口数成正比,从而xn+1 xn r xn xn+1 a xn其中 a=r+1.设xn是某人类群体在第n个时间段(例如年)末时的总数，若在单位时间段内人口相对增长率为r（出生率与死亡率之差）,那么人口增长数与原人即是常用的计算公式常用的计算公式第4页/共30页记 g(x)=a x，则是函数迭代 xn=a xn1=a2xn2=an x0 于是Malthus
3、的结论：人口增长呈几何级数 据统计1951年-1961年人口增长率为2%，则a=1.02,1.0235 1.99989955 与近年统计结果有误差，由a 1,xn趋向无穷，模型在人口长期预测方面必定是失效的.xn1 g(xn)容易得到 所以,人口约35年增加一倍，则与17001961年世界人口统计结果一致。以此方法，到2626年世界人口将是1961年的19番，地球表面积为1.731014m2的话,每平方米就有9个人。第5页/共30页指数增长模型的应用及局限性 与1919世纪以前欧洲一些地区人口统计数据吻合 适用于1919世纪后迁往加拿大的欧洲移民后代 可用于短期人口增长预测 不符合1919世纪
4、后多数地区人口增长规律 不能预测较长期的人口增长过程1919世纪后人口数据人口增长率r r不是常数(逐渐下降)第6页/共30页2、阻滞增长模型(Logistic模型)人口增长到一定数量后，增长率下降的原因：资源、环境等因素对人口增长的阻滞作用且阻滞作用随人口数量增加而变大假设假设r固有增长率(x很小时)xm人口容量（资源、环境能容纳的最大数量）r是x的减函数第7页/共30页dx/dtx0 xmxm/2x(t)S形曲线,x增加先快后慢tx0 xmx0 xm/2第8页/共30页参数估参数估计计用指数增长模型或阻滞增长模型作人口预报，必须先估计模型参数 r r 或 r,xr,xm m 利用统计数据用
5、最小二乘法作拟合例：美国人口数据（单位-百万）1860 1870 1880 1960 1970 1980 1990 31.4 38.6 50.2 179.3 204.0 226.5 251.4阻滞增长模型(Logistic模型)r=0.2557,xm=392.1第9页/共30页连续形连续形式式t t,x,xN,x=NN,x=N是稳定平衡点(与r r大小无关)离散形离散形式式x(t):某种群 t 时刻的数量(人口)yk：某种群第k k代的数量(人口)若若yk=N,则则yk+1,yk+2,=N讨论平衡点的稳定性，即讨论平衡点的稳定性，即k,ykN？y*=N 是平衡点是平衡点第10页/共30页 生存
6、资源是重要的因素，修改模型为：yn+1 yn=r yn b yn2 b yn2为竞争(约束)项，r、b 称生命系数,则 yn+1=a yn-byn2 ，（a=r+1）这是一个非线性映射的迭代 f1(x)=ax-bx2数据观察（利用Mathematica）a=1.029;b=1.48654*10(-11);f1x_:=a*x-b*x2;gg=1979,x1979/108;Forn=1979;x1979=9.7542*108,n,xn+1/108;gg=Appendgg,n+1,xn+1/108;sawline=Linegg;sawgraph=ShowGraphicssawline,Axes-Tr
7、ue与统计数字接近第11页/共30页1980 9.89564 1981 10.0371982 10.17841983 10.31951984 10.46051985 10.60121986 10.74161987 10.88151988 11.02111989 11.16011990 11.29861991 11.43651992 11.57381993 11.71031994 11.84601995 11.98091996 12.11501997 12.24821998 12.38031999 12.51152000 12.64172001 12.77012002 12.89862003 1
8、3.0245第12页/共30页LogisticLogistic映射映射 (Robert.May的研究）f(x)=a x(1-x)，x 在0,1内变化xn+1=f(xn)从0,1内点x0出发，由Logistic映射的迭代形成xn=f n(x0),n=0,1,2,序列xn称为x0的轨道种群数的模型简化：相应的迭代为了一个序列，即什么是混沌?第13页/共30页数值迭代数值迭代 （a 逐渐增加，迭代会有何结果）1倍周期分叉现象 当0a 1时，由于0 xnaxn+1 当1a3时，任何(0,1)中初始值的轨道趋于 事实上，由方程 ax(1-x)=x 可以解出两个不动点x1*=1-1/a,x2*=0,一个稳
9、定(吸引),另一个不稳定，轨道xn趋向稳定点，xn 0 食物不够物种逐渐灭亡x*=1-1/a其中x*是方程f(x)=x的解，为映射f 的不动点（周期1点)例：a=1.5时 xn 1/3.种群数量趋于稳定第14页/共30页 这两个数满足 当3a1+61/2时，xn 绕着两个数 x3*，x4*振动,x2k-1 0.799455 x2k 0.513045 当1+61/2a3.5440903506时,从任意的点x0出 x4k 0.44391661 x4k+1 0.84768002 x4k+2 0.44596756 x4k+3 0.85242774也称为周期2点,对应轨道称周期2轨道.(原来周期例 a=
10、3.2点失稳)发的轨道将逐渐沿着四个数值振动例 a=3.45第15页/共30页这四个数满足称为周期4点，对应轨道称周期4轨道(原有周期点 若a再增大，周期4点又会失稳，而产生新的稳定分叉值如何求?又失稳)周期8点，这个周期不断加倍的过程将重复无限次，会依次出现周期16点，周期32点.，(请考虑什么是周期n)这种过程称为倍周期分叉.相应的分叉值c1=3,c2=1+61/2构成一个单调增加的数列ck.其极限值为c*=3.569945557391。第16页/共30页a=3.2;f1x_:=a*x-a*x2;gg=0,x0;Forn=0;x0=0.3,n TrueMATHEMATICAMATHEMAT
11、ICA程程序序第17页/共30页分叉图分叉图第18页/共30页混沌的特点混沌的特点 当c*a Identity;AppendTopointlist,x0,0;Fori=1,i,fvar;var=fvar;fvar=u*var*(1-var);p2=ListPlotpointlist,PlotJoined-True,DisplayFunction-Identity;Showp1,p2,DisplayFunction-$DisplayFunction;hh=Takepointlist,-20;MatrixFormhhIterGeo3.9,0.302第22页/共30页11xnxn1x0 x1x1x2第23页/共30页 1a3 从(0,1)中任何初值出发的轨道趋向不动点(周期1点)第24页/共30页 3a61/2+1 从任何初值出发的轨道趋向周期2点第25页/共30页61/2+1a 3.54409035从任何初值出发的轨道趋向周期4点第26页/共30页 a=3.58轨道进入浑沌状态第27页/共30页 a=4 轨道的浑沌性表现充分第28页/共30页出现在各个领域的一种现象:数学、物理、“侏罗纪公园”中的恐龙重现生物、金融、经济、管理等等：宇宙的起源 龙卷风的产生、厄尔尼诺现象 东南亚金融危机爆发 第29页/共30页感谢您的欣赏！第30页/共30页