1、精選優(yōu)質(zhì)文檔-傾情為你奉上數(shù)學(xué)建模人口增長(zhǎng)模型摘要:人口的增長(zhǎng)是當(dāng)前世界上引起普遍關(guān)注的問題作為世界上人口最多的國(guó)家，我國(guó)的人口問題是十分突出的由于人口基數(shù)大盡管我國(guó)已經(jīng)實(shí)行了20多年的計(jì)劃生育政策人口的增長(zhǎng)依然很快，巨大人口壓力會(huì)給我國(guó)的社會(huì) 政治經(jīng)濟(jì)醫(yī)療就業(yè)等帶來了一系列的問題。因此研究和解決人口問題在我國(guó)顯得尤為重要。我們經(jīng)常在報(bào)刊上看見關(guān)于人口增長(zhǎng)預(yù)報(bào)，說到本世紀(jì)，或下世紀(jì)中葉，全世界的人口將達(dá)到多少億。你可能注意到不同報(bào)刊對(duì)同一時(shí)間人口的預(yù)報(bào)在數(shù)字商場(chǎng)有較大的區(qū)別，這顯然是由于用了不同的人口整張模型計(jì)算出來的結(jié)果。人類社會(huì)進(jìn)入20世紀(jì)以來，在科學(xué)和技術(shù)和生產(chǎn)力飛速發(fā)展的同時(shí)世界人口也

2、以空前的規(guī)模增長(zhǎng)。人口每增加十億的時(shí)間，有一百年縮短為十幾年。我們賴以生存的地球已經(jīng)攜帶著他的60億子民踏入下一個(gè)世紀(jì)。長(zhǎng)期以來，人類的繁殖一直在自然地進(jìn)行著，只是由于人口數(shù)量的迅速膨脹和環(huán)境質(zhì)量的急劇惡化，人們才猛然醒悟，開始研究人類和自然的關(guān)系、人口數(shù)量的變化規(guī)律以及如何驚醒人口控制等問題。 本論文中有兩個(gè)模型： （1） ：中國(guó)人口的指數(shù)增長(zhǎng)模型，并用該模型進(jìn)行預(yù)測(cè)，與實(shí)際人口數(shù)據(jù)進(jìn)行比較。（2）：中國(guó)人口的Logistic圖形，標(biāo)出中國(guó)人口的實(shí)際統(tǒng)計(jì)數(shù)據(jù)進(jìn)行比較。 而且利用MATLAB圖形 ，標(biāo)出中國(guó)人口的實(shí)際統(tǒng)計(jì)數(shù)據(jù)，并畫出兩種模型的預(yù)測(cè)曲線。 關(guān)鍵字：人口預(yù)測(cè)；Malthus模型；L

3、ogistic模型；MATLAB軟件專心-專注-專業(yè)1、 問題背景及重述1.1問題的背景中國(guó)是一個(gè)人口大國(guó)，人口問題始終是制約我國(guó)發(fā)展的關(guān)鍵因素之一。我國(guó)自1973年全面推行計(jì)劃生育以來，生育率迅速下降，取得了舉世矚目的成就，但全面建設(shè)小康社會(huì)仍面臨著人口的形勢(shì)和嚴(yán)峻挑戰(zhàn)。隨著我國(guó)經(jīng)濟(jì)的發(fā)展、國(guó)家人口政策的實(shí)施，未來我國(guó)人口高峰期到底有多少人口，專家學(xué)者們的預(yù)測(cè)結(jié)果不一。因此，根據(jù)已有數(shù)據(jù)，運(yùn)用數(shù)學(xué)建模的方法，對(duì)中國(guó)人口做出分析和預(yù)測(cè)是一個(gè)重要問題。 1.2 問題的重述下表列出了中國(guó)19821998年的人口統(tǒng)計(jì)數(shù)據(jù)，取1982年為起始年(t=0)，1982年的人口萬(wàn)人，人口自然增長(zhǎng)率為14，以

4、36億作為我國(guó)人口的容納量，試建立一個(gè)較好的人口數(shù)學(xué)模型并給出相應(yīng)的算法和程序，并與實(shí)際人口進(jìn)行比較。時(shí)間(年)198219831984198519861987人口(萬(wàn)人)時(shí)間(年)198819891990199119921993人口(萬(wàn)人)時(shí)間(年)19941995199619971998人口(萬(wàn)人)2、 問題分析對(duì)于人口增長(zhǎng)的問題，其影響因素有很多，比如：人口基數(shù)，出生率，死亡率，人口男女比例，人口年齡結(jié)構(gòu)的組成，人口的遷入率和遷出率，人口的生育率和生育模式，國(guó)家的醫(yī)療發(fā)展情況，國(guó)家的政治策略等眾多的因素。如果把這些因素都要考慮進(jìn)去，則該問題根本無從下手。因此，應(yīng)該根據(jù)中國(guó)人口自身發(fā)展的特

5、點(diǎn)，選取相應(yīng)的能夠體現(xiàn)我國(guó)人口發(fā)展特點(diǎn)的模型。人口發(fā)展模型有連續(xù)形式和離散形式，因?yàn)轭}目所給的 圖1 1790-1980年間美國(guó)的人口數(shù)據(jù)圖數(shù)據(jù)是每個(gè)年份的具體數(shù)據(jù)，可以將這些數(shù)據(jù)視為連續(xù)的。根據(jù)表格中的數(shù)據(jù)，我們使用MATLAB編程（附錄1）畫出散點(diǎn)圖。圖1中國(guó)19821998年的人口數(shù)據(jù)圖從圖中我們可以看到人口數(shù)在 19821998 年是呈增長(zhǎng)趨勢(shì)的，且增長(zhǎng)趨勢(shì)類似于指數(shù)型增長(zhǎng)，因此，我們可以先建立一個(gè)指數(shù)增長(zhǎng)模型（Malthus模型）。但是，由于地球上的資源是有限的，它只能提供一定數(shù)量的生命生存所需的條件，因此人口不可能無限制增加。隨著人口數(shù)量的增加，自然資源，環(huán)境條件等對(duì)人口再增長(zhǎng)的限

6、制作用將越來越顯著。于是我們假設(shè)在人口較少時(shí)，可以把人口增長(zhǎng)率看成常數(shù)，但隨著人口的增加，我們應(yīng)該把人口增長(zhǎng)率視為一個(gè)隨著人口增加而減小的量，從而我們可以將模型一（Malthus模型）優(yōu)化為一個(gè)阻滯增長(zhǎng)模型（Logistic模型）。3、 模型假設(shè)（1） Malthus模型假設(shè)我國(guó)人口的增長(zhǎng)符合人口指數(shù)增長(zhǎng)的規(guī)律，即滿足Malthus模型的兩個(gè)前提：第一， 食物是人類生存所必需的；第二， 兩性間的情欲是必然的，而且?guī)缀鯐?huì)保持現(xiàn)狀。從這兩個(gè)“人類本性的固定法則” 出發(fā)，可以得出一個(gè)最基本的經(jīng)濟(jì)比例：食物或生活資料的增長(zhǎng)與人口的增殖之間的關(guān)系。馬爾薩斯說， 人口的增殖比生活資料增長(zhǎng)的要快，人口是按幾

7、何級(jí)數(shù)增長(zhǎng)的，而生活資料則只按算術(shù)級(jí)數(shù)增長(zhǎng)。但是，馬爾薩斯并不認(rèn)為這兩個(gè)級(jí)數(shù)就是人口規(guī)律的反映，他提出，保持兩個(gè)級(jí)數(shù)平衡的唯一出路就是抑制人口的增長(zhǎng)。他把所謂支配人類命運(yùn)的永恒的人口自然法則， 歸納成以下三個(gè)定理。三個(gè)定理：第一點(diǎn)是人口的制約原理， 說明人口與生活資料之間必然存在某種正常的比例， 即“人口的增長(zhǎng)， 必然要受到生活資料的限制”；第二點(diǎn)是人口的增殖原理， 即“生活資料增加， 人口也常隨著增加”；第三點(diǎn)是人口原理的核心， 稱之為人口的均衡原理， 即“占優(yōu)勢(shì)的人口繁殖力為貧困和罪惡所抑制，因而使現(xiàn)實(shí)的人口得以與生活資料保持平衡”。這個(gè)原理與前兩個(gè)原理是緊密相連的，它說明人口與生活資料之

8、間最終將實(shí)現(xiàn)均衡， 但是這種均衡不是自然實(shí)現(xiàn)的，而是種種“抑制”的產(chǎn)物。所以，Malthus模型假設(shè)條件如下：1 設(shè)P(t)表示t時(shí)刻的人口數(shù)，且P(t)連續(xù)可微。2 人口的增長(zhǎng)率r是常數(shù)(增長(zhǎng)率=出生率-死亡率)。3 人口數(shù)量的變化是封閉，即人口數(shù)量的增加與減少只取決于人口中個(gè)體的生育與死亡，且每一個(gè)都具有同樣的的生育能力與死亡率。 （二）Logistic模型由于地球上的資源有限，當(dāng)人口數(shù)量發(fā)展到一定階段后，會(huì)產(chǎn)生一系列問題，如食物短缺、居住和交通擁擠等。另外，隨著人口密度的增加，疾病將會(huì)增多，死亡率會(huì)上升，因此，人口的增長(zhǎng)率不會(huì)是Malthus所假設(shè)的是一個(gè)常數(shù)不改變，而是會(huì)隨著人口數(shù)量增

9、加而減少。假設(shè)增長(zhǎng)率r表示P(t)的函數(shù)r(p)，且r(p)為p的減函數(shù)。1. 設(shè)r(p)為p的線性函數(shù)，r(p)=r-kp。2. 自然資源與環(huán)境條件所能容納的最大人口數(shù)為Pm，即當(dāng)P=Pm時(shí)，增長(zhǎng)率r(p)=0。4、 變量說明符號(hào)表示意義P人口數(shù)量t年份r人口自然增長(zhǎng)率人口最大容納量起始年人口5、 模型建立與模型求解5.1 Malthus模型由假設(shè)一，時(shí)刻到時(shí)刻人口增量為于是可得 由分離變量法解得模型的解為對(duì)該模型兩邊同時(shí)取對(duì)數(shù)得到一次線性擬合函數(shù)取表中1982到1998年的數(shù)據(jù)在MATLAB中M文件（附錄2）進(jìn)行線性最小二乘擬合可得出:f =0. t - 14.5121所以可知r=0.,p

10、(t)=*exp(0.*(t-1982)用MATLAB進(jìn)行指數(shù)擬合得到下圖圖2可以看出擬合曲線基本吻合，但是隨著時(shí)間t的增加其誤差逐漸加大，所以需要對(duì)其修正。5.2 Logistic模型由假設(shè)二可知，記p(t)是第t年的人口數(shù)量，人口增長(zhǎng)率r(p)是p的線性函數(shù)，r(p)=r-kp。最大人口容納量為Pm。即當(dāng)P=Pm時(shí)，增長(zhǎng)率r(p)=0。所以，(5.2.1)同樣利用分離變量法求得其解(5.2.2)根據(jù)（5.2.1）式作出的曲線圖（圖1）以及由（5.2.2）式作出p-t曲線圖（圖2）O圖1 曲線圖圖2 p-t曲線圖從上述曲線圖以及表達(dá)式中，我們可以總結(jié)出如下規(guī)律：，它表明不管人口初始狀態(tài)是什么

11、樣，人口總數(shù)最終都將趨于最大人口容納量。當(dāng)p(t)>pm時(shí)，<0；當(dāng)p(t)<pm時(shí)，>0。它表明當(dāng)人口數(shù)量超過最大人口容納量時(shí)，人口數(shù)量將減少，當(dāng)人口數(shù)量小于最大人口容納量時(shí)，人口數(shù)量將增加。人口變化率在時(shí)取到最大值，即人口總數(shù)達(dá)到極限值一半之前是加速生長(zhǎng)的，經(jīng)過此點(diǎn)后，增長(zhǎng)率會(huì)逐漸減小至0。采用非線性最小二乘估計(jì)法對(duì)參數(shù)r和pm進(jìn)行估計(jì)，通過使用matlab編寫程序（附錄4）可得：r =0.01137，pm =3.7465e+04用MATLAB擬合圖像如下圖36、 模型檢驗(yàn)及結(jié)果分析經(jīng)過前面模型建立的工作，已建立出Malthus模型和Logisic模型?，F(xiàn)在根據(jù)所建

12、立的模型預(yù)測(cè)相關(guān)年份的人口數(shù)量，并與實(shí)際人口數(shù)量相比較以檢驗(yàn)?zāi)Ｐ偷膬?yōu)劣性。Malthus模型與Logistic模型對(duì)我國(guó)人數(shù)據(jù)的擬合結(jié)果年份實(shí)際人口/萬(wàn)計(jì)算人口P1計(jì)算人口P21982198319841985104201986198719881989199019911992199319941995199619971998對(duì)表中數(shù)據(jù)進(jìn)行分析可知：對(duì)于短期預(yù)測(cè)，這兩個(gè)模型基本一致，但使用模型一更簡(jiǎn)單；對(duì)于中長(zhǎng)期預(yù)測(cè)，模型二要強(qiáng)于模型一。7、 模型評(píng)價(jià)與推廣1、 優(yōu)點(diǎn)：首先我們采用圖表結(jié)合法，比較直觀地表達(dá)出題中所給的信息，并據(jù)此得出了人口增長(zhǎng)的基本規(guī)律。根據(jù)所給出的數(shù)據(jù)，對(duì)其進(jìn)行分析得出了人口增長(zhǎng)

13、率與人口總數(shù)的線性關(guān)系，從而建立了人口阻滯增長(zhǎng)模型，對(duì)未來人口數(shù)的預(yù)測(cè)作出了較為準(zhǔn)確的判斷。模型一是依據(jù)英國(guó)神父T·Malthus的發(fā)現(xiàn)建立了指數(shù)型增長(zhǎng)模型，經(jīng)過我們實(shí)際數(shù)據(jù)的檢驗(yàn)，發(fā)現(xiàn)其人口早期的增長(zhǎng)情況與Malthus模型的預(yù)測(cè)基本相符，然而隨著時(shí)間的增加，該模型的預(yù)測(cè)結(jié)果明顯出現(xiàn)了不合理性。其原因就是我們將人口增長(zhǎng)率視為常數(shù)，因此需要對(duì)r進(jìn)行修正。所以，我們將r表示為p的減函數(shù)，從而推導(dǎo)建立了模型二2、 缺點(diǎn)：本文對(duì)模型一中的參數(shù)只做了線性估計(jì)，所以其計(jì)算結(jié)果與實(shí)際誤差較大模型二中僅考慮了r與p的關(guān)系是線性的，沒有考慮非線性關(guān)系8、 參考文獻(xiàn)1.司守奎，孫兆亮，孫璽菁，周剛，仲

14、維杰，康淑瑰.數(shù)學(xué)建模算法與應(yīng)用（第二版）.國(guó)防工業(yè)出版社，2016年2.姜啟源，謝金星，葉俊.數(shù)學(xué)模型（第四版）.高等教育出版社·北京.2011年3.儲(chǔ)昌木，沈長(zhǎng)春.數(shù)學(xué)建模及其應(yīng)用.西南交通大學(xué)出版社·成都.2015年4 胡守信，李柏年.基于MATLAB的數(shù)學(xué)實(shí)驗(yàn)M.北京科學(xué)出版社.2004年6月5 揚(yáng)啟帆，康旭升，等.數(shù)學(xué)模型M.北京：高等教育出版社.2006年5月6 于學(xué)軍.中國(guó)人口科學(xué)2000年第2期，時(shí)間：2000-4-6，中國(guó)人口信息網(wǎng).附錄：1) syms x yx0=1982:1:1998y0= xlabel('x')ylabel('

15、;y')plot(x0,y0,'*')2) t=1982:1:1998;y=log() log() log() log() log() log() log() log() log() log() log() log() log() log() log() log() log()p1=polyfit(t,y,1);f=poly2str(p1,'t')3) syms x y px0=1982:1:1998y0= xlabel('x')ylabel('y')plot(x0,y0,'*')hold ont=1982:1:1998p=.*exp(0.*(t-1982)plot(x,p,'r','LineWidth',0.5);legend('原始數(shù)據(jù)散點(diǎn)圖','指數(shù)擬合曲線');grid o