數(shù)學(xué)建模:人口增長(zhǎng)模型

模型目標(biāo):
通過給定的一組人口增長(zhǎng)數(shù)據(jù),預(yù)測(cè)后續(xù)的人口增長(zhǎng)情況.

一、指數(shù)增長(zhǎng)模型

假設(shè)增長(zhǎng)率不變:

若已知人口年增長(zhǎng)率為r,今年人口為 $x_0$,預(yù)測(cè)k年后的人口可以用簡(jiǎn)單的公式得到:
$$x_k=x_0(1+r{)}^k$$
*以美國(guó)人口為例,數(shù)據(jù)點(diǎn)取下表:

用matlab輸入好數(shù)據(jù):

p = [3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76 92 105.7 122.8 131.7 150.7 179.3 203.2 226.5 248.7 281.4];

該公式得到的人口增長(zhǎng)模型如下圖(matlab作圖):

選取的增長(zhǎng)率0.2與實(shí)際的人口數(shù)據(jù)點(diǎn)貼合度差,于是進(jìn)一步建立模型,找到偏差更小的增長(zhǎng)率r

精確化增長(zhǎng)率r

把人口看作關(guān)于時(shí)間的可微函數(shù) $x(t)$,記 $x(0)=x_0$
$r{x}_0$ 即為單位時(shí)間 $x(t)$ 的增長(zhǎng)量
所以得到微分方程:
$$\frac{dx}{dt}=rx,x(0)=x_0$$
解得:
$$x(t)=x_0{e}^{rt}$$
上述公式即位指數(shù)增長(zhǎng)模型

數(shù)據(jù)擬合:

接下來(lái)對(duì)指數(shù)增長(zhǎng)模型的參數(shù)進(jìn)行估計(jì),即數(shù)據(jù)擬合

·法一

用人口數(shù)據(jù)和線性最小二乘法,對(duì)上式取對(duì)數(shù):
$$y=rt+a,y=lnx,a=ln{x}_0$$
用matlab進(jìn)行編程擬合:

%取1790年為t=0,數(shù)據(jù)點(diǎn)取到2000年 t=0:10:210; lnp = log(p); cftools

解得:
$$r=0.202,x_0=e^{1.8}=6.049$$
帶入得$x(t)=6.049e^{0.0202t}$

·法2
對(duì)人口數(shù)據(jù)做數(shù)值微分,計(jì)算平均值r‘,x0直接選用原始數(shù)據(jù).
函數(shù)在各點(diǎn)的近似導(dǎo)數(shù)值為(數(shù)值微分中點(diǎn)公式):
$$\begin{gathered} x^{\prime }(t_k)=\frac{x_{k+1}?x_{k?1}}{2\Delta t}, \\ (k=1,2,3,...,n?1) \end{gathered}$$
$$x^{\prime }(0)=\frac{4x_1?3x_0?x_2}{2\Delta t},x^{\prime }(n)=\frac{?4x_{n?1}+3x_n+x_{n?2}}{2\Delta t}$$
那么增長(zhǎng)率為: $$r(t_k)=\frac{x^{\prime }(t_k)}{x(t_k)}$$

*公式相關(guān)推導(dǎo)可參考數(shù)值計(jì)算方法 第六章 數(shù)值積分和數(shù)值微分

增長(zhǎng)率 $r_k=\frac{x^{\prime }(t_k)}{x(t_k)}$再取平均值得到r = 0.0205

改進(jìn)的指數(shù)增長(zhǎng)模型

上述模型對(duì)于增長(zhǎng)率r不變的假設(shè)導(dǎo)致預(yù)測(cè)曲線與實(shí)際偏差較大
所以改進(jìn)模型中假設(shè) r 為 t 的函數(shù) $r(t)$ , 根據(jù)上述法2的 $x^{\prime }(t)$ 畫r—t圖像:

r=[]; for i=1:22if i == 1r(i)=(4*p(i+1)-3*p(i)-p(i+2))/(20*p(i));elseif i == 22r(i)=(-4*p(i-1)+3*p(i)+p(i-2))/(20*p(i));elser(i)=(p(i+1)-p(i-1))/(20*p(i));end end plot(t,r,'.','MarkerSize',20); ylim ([0.005,0.04]); xlabel('t'); ylabel('增長(zhǎng)率r');

根據(jù)散點(diǎn)圖假設(shè)$$r(t)=r_0?r_{1}t$$的線性函數(shù),用最小二乘法線性擬合得到:
$$r_0=0.03252,r_1=0.0001143$$
根據(jù)微分方程:
$$dx/dt=r(t)x$$
$$?x(t)=x_0{e}^{(r_{0}t?r_1{t}^2/2)}$$
把擬合后的參數(shù)帶入:

xt2 = 3.9.*exp(0.03252.*t-0.0001143.*t.^2./2); plot(t,xt2,'LineWidth',1);

最終根據(jù)改進(jìn)模型預(yù)測(cè)的2010年人口為290million,與實(shí)際數(shù)據(jù)281.4吻合度較高.顯然改進(jìn)后的模型優(yōu)于前兩者.

二、logistic模型

改進(jìn)的指數(shù)增長(zhǎng)模型中增長(zhǎng)率線性下降,但沒有體現(xiàn)其下降的相關(guān)影響因素,只是以時(shí)間為變量.logistic模型考慮了自然資源,環(huán)境等對(duì)人口增長(zhǎng)的阻滯作用.

未完…

總結(jié)

以上是生活随笔為你收集整理的数学建模:人口增长模型的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。