人口指数增长模型及Logistic模型.docVIP

范文范例 值得参考 word完美 整理版 根据美国人口从1790年到1990年间的人口数据（如下表），确定人口指数增长模型和Logistic模型中的待定参数，估计出美国2010年的人口，同时画出拟合效果的图形。 表1 美国人口统计数据 年 份 1790 1800 1810 1820 1830 1840 1850 人口(×106) 3.9 5.3 7.2 9.6 12.9 17.1 23.2 年 份 1860 1870 1880 1890 1900 1910 1920 人口(×106) 31.4 38.6 50.2 62.9 76.0 92.0 106.5 年 份 1930 1940 1950 1960 1970 1980 人口(×106) 123.2 131.7 150.7 179.3 204.0 226.5 提示： 指数增长模型： Logistic模型： 解：模型一:指数增长模型。Malthus 模型的基本假设下，人口的增长率为常数，记为r，记时刻t的人口为 ，（即为模型的状态变量）且初始时刻的人口为，因为由假设可知 经拟合得到： 程序： t=1790:10:1980; x(t)=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 ]; y=log(x(t));a=polyfit(t,y,1) r=a(1),x0=exp(a(2)) x1=x0.*exp(r.*t); plot(t,x(t),r,t,x1,b) 结果：a = 0.0214 -36.6198 r= 0.0214 x0= 1.2480e-016 所以得到人口关于时间的函数为：，其中x0 = 1.2480e-016， 输入：t=2010; x0 = 1.2480e-016; x(t)=x0*exp(0.0214*t) 得到x(t)= 598.3529。即在此模型下到2010年人口大约为598.3529 。 模型二：阻滞增长模型（或 Logistic 模型） 由于资源、环境等因素对人口增长的阻滞作用，人口增长到一定数量后，增长率会下降，假设人口的增长率为 x 的减函数，如设，其中 r 为固有增长率 (x 很小时 ) ，为人口容量（资源、环境能容纳的最大数量）， 于是得到如下微分方程： 建立函数文件curvefit_fun2.m function f=curvefit_fun2 (a,t) f=a(1)./(1+(a(1)/3.9-1)*exp(-a(2)*(t-1790))); 在命令文件main.m中调用函数文件curvefit_fun2.m % 定义向量（数组） x=1790:10:1990; y=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76 ... 92 106.5 123.2 131.7 150.7 179.3 204 226.5 251.4]; plot(x,y,*,x,y); % 画点，并且画一直线把各点连起来 hold on; a0=[0.001,1]; % 初值 % 最重要的函数，第1个参数是函数名（一个同名的m文件定义），第2个参数是初值，第3、4个参数是已知数据点 a=lsqcurvefit(curvefit_fun2,a0,x,y); disp([a= num2str(a)]); % 显示结果 % 画图检验结果 xi=1790:5:2020; yi=curvefit_fun2(a,xi); plot(xi,yi,r); % 预测2010年的数据 x1=2010; y1=curvefit_fun2(a,x1) hold off 运行结果： a=311.9531 0y1 =267.1947 其中a(1)、a(2)分别表示中的和，y1则是对美国美国2010年的人口的估计。 第二题： 问题重述： 一垂钓俱乐部鼓励垂钓者将钓上的鱼放生，打算按照放生的鱼的重量给与鼓励，俱乐部只准备了一把软尺用于测量，请你设计按照测量的长度估计鱼的重量的方法。假定鱼池中只有一种鲈鱼，并且得到8条鱼的如下数据（胸围指鱼身的最大周长）： 身长(cm) 36.8 31.8 43.8 36.8 32.1 45.1 35.9 32.1 重量(g) 765 482 1162 737 482 1389 652 454 胸围(cm) 24.8 21.3 27.9 24.8 21.6 31