数学建模作业6指数增长模型和Logistic模型.docxVIP

佛山科学技术学院 上机报告 课程名称 数学建模 上机项目 指数增长模型和Logistic模型 专业班级 姓名 学号 一、问题提出 人口问题是当前世界上人们最关心的问题之一。认识人口数量的变化规律，作出较准确的 预报，是有效控制人口增长的前提。 要求：分别建立并求解两个最基本的人口模型，即：指数增长模型和Logistic模型，并利 川表1给出的近两百年的人口统计数据，価出图形拟合数据，对模型做出检验，最后用它预报 2000年的人口。 表1人口统计数据 年（公元） 1790 1800 181() 1820 183() 1840 1850 人口（百万） 3.9 5.3 7.2 9.6 12.9 17.1 23.2 年（公元） 1860 1870 1880 1890 1900 1910 1920 人口（百万） 31.4 36 50.2 62.9 76.0 92.0 106.5 年（公兀） 1930 1940 1950 1960 1970 1980 1990 人口（百万） 123.2 131.7 150.7 179.3 204.0 226.5 251.4 模型一：指数增长（Malthus）模型: （1）模型假设 常用的计算公式：今年人口工0,年增长率厂R年后人口为无二兀。（1 +厂） 假设：人口壇长率r是常数（或单位时间内人口的壇长量与当时的人口成正比）。 符号说明： 兀。~匸0时的人口数，垃）~时刻1的人口数 模型建立 (显示模型函数的构造过程) 由于量大，兀⑴口J看作连续、口J微函数，r到f +△/时间内人口的增量为 = 于是兀⑴满足微分方程 dx dt (1) x(0) = x0 模型求解 (显示模型的求解方法、步骤及运算程序、结果) 解微分方程(1)得 x(r) = xQerl (2) 当/Too时，兀⑴Too(厂0),即随着时间增加，人口按指数规律无限增长。 (4)模型的参数估计 要用模型的结果(2)式来预报人口，必须对常数r进行估计，可以用表1的数 据通过拟合得到。取兀。=3.9,通过(2)式以及表中1790-1990的数据进行最小二 乘法拟合得r=0. 2169.程序如下: 模型求解：取初始值x(0)=3.9 Matlab 程序： 建立M文件volum.m,如下： function y=volum(beta,t) y=3?9*exp(beta(1)*t)； ? Editor - C:\Users\Administrator\Documents\MATLAB\volum.m File Edit Text Go Cell Tools Debug Desktop Window Help 第■唔町e ◎勢打鼻??已|闵▼档老fin :?回厂目 目b回 ■ l.o + Li x 癖疵级 function y=volum(bett) 一 y=3. 9*exp (beta(l)*t): 3 再建立rl?m程序，如下： t = 0:1:2 0; x= [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 251.4]; beta0= [0.01]； [beta,r,J]=nlinfit(t1,x1,1volum1,betaO)； beta Editor ? C:\User5\Administrator\Documents\MATLAB\rl. File Edit 「ext ( jo Cell Tool s Debug Desktop Window Help Base ▼ | fx :?目「目 ?回b目 1.0 + T 1.1 x囁囁Q ? t=0:l:20: - x=[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 - bet a0= [0.01]： - [beta, r, J]=nlinfit (t,x?/ volm ;betaO); - beta; A 150.7 179.3 204.0 226. 5 251.4]; 结果为 beta =0.2169 J MATLAB 7.10.0 (R2010a) File Edit Debug Parallel Desktop Window Help ° 1 晶匕 为 5 G 西 Ml ? Current Folder: C:\Users\Administrator\Documents\MA Shortcuts 4 How to Add 上]Whats New Current Folder 计口冇 X Command W