2012數(shù)學(xué)建模深圳杯A答案

答卷編號(hào)（參賽學(xué)校填寫(xiě)）：答卷編號(hào)（競(jìng)賽組委會(huì)填寫(xiě)）：論文題目：深圳人口與醫(yī)療需求預(yù)測(cè)（A）組別：本科生參賽學(xué)校：報(bào)名序號(hào)：參賽隊(duì)員信息(必填)： 答卷編號(hào)（競(jìng)賽組委會(huì)填寫(xiě)）：評(píng)閱情況（省賽評(píng)閱專家填寫(xiě)）：省賽評(píng)閱1：省賽評(píng)閱2：省賽評(píng)閱3：省賽評(píng)閱4：省賽評(píng)閱5：深圳市人口與醫(yī)療需求預(yù)測(cè)模型摘要：人口與醫(yī)療問(wèn)題是關(guān)系到國(guó)計(jì)民生的大問(wèn)題，能夠合理而準(zhǔn)確地預(yù)測(cè)就顯得非常重要。但不同城市有不同的人口特點(diǎn)，本文在吸取前人經(jīng)驗(yàn)的基礎(chǔ)上，以深圳的人口為依托提出了一些新的簡(jiǎn)單而實(shí)用方法，希望能為政府決策提供幫助。針對(duì)深圳市人口結(jié)構(gòu)中非戶籍人口比重大，流動(dòng)人口多這一特點(diǎn)，我們采用了灰色GM(1,1)模型，通過(guò)matlab對(duì)深圳市自2001至2010年的數(shù)據(jù)進(jìn)行擬合，發(fā)現(xiàn)其人口變化近似呈線性增長(zhǎng)，線性相關(guān)系數(shù)高達(dá)0.99，我們就此認(rèn)定其為線性相關(guān)并給出線性方程。同理，針對(duì)其非戶籍人口，我們進(jìn)行matlab擬合發(fā)現(xiàn)，其為非線性相關(guān)，并得出相關(guān)函數(shù)。通過(guò)模擬出的常住人口與非戶籍人口的函數(shù)，我們可以很容易的得出深圳市的人口數(shù)量變化情況，同時(shí)我們以非戶籍人口與常住人口的函數(shù)之比作為深圳市人口結(jié)構(gòu)的變化，通過(guò)作圖發(fā)現(xiàn)，深圳市非戶籍人口正逐年下降，這正與官方以及媒體報(bào)道深圳市產(chǎn)業(yè)轉(zhuǎn)型相對(duì)應(yīng)。由于深圳市人口結(jié)構(gòu)中外來(lái)人口比例接近76%，而且外來(lái)人口中以青壯年居多，可以認(rèn)為在較短時(shí)間內(nèi)（十年內(nèi)）外來(lái)人口年齡結(jié)構(gòu)近似不變，同時(shí)當(dāng)?shù)貞艏丝谝驗(yàn)槭軞v史條件影響，人口年齡結(jié)構(gòu)在短期內(nèi)也不會(huì)發(fā)生較大變化，所以16.——胃癌患病人數(shù)；17.——肺癌在市級(jí)醫(yī)院的病床需求；18.——肺癌在區(qū)級(jí)醫(yī)院的病床需求；——胃癌在市級(jí)醫(yī)院的病床需求；——胃癌在區(qū)級(jí)醫(yī)院的病床需求；模型建立及求解問(wèn)題=1\*ROMANI：深圳市最近十年常住人口、非常住人口變化特征與未來(lái)十年發(fā)展情況。首先，我們采用灰色GM(1,1)模型.采用原因：灰色模型適用于小樣本、貧信息、內(nèi)在規(guī)律未充分外露的系統(tǒng)，按適當(dāng)辦法處理原始數(shù)據(jù)后得到規(guī)律性較強(qiáng)的生成函數(shù).本題給出的常住人口、非常住人口數(shù)據(jù)受到難以區(qū)分的多重因素影響，且數(shù)據(jù)量較小，適用于灰色模型.由于常住人口數(shù)量受歷史影響較大，不易發(fā)生較大變化，且在數(shù)據(jù)處理中發(fā)現(xiàn)了較強(qiáng)的線性關(guān)系，我們之后采用了一元線性擬合來(lái)簡(jiǎn)化模型；而非常住人口受各方面因素影響較大，仍保持灰色模型不變。最后我們得出了與常住/非常住人口相關(guān)的人口結(jié)構(gòu)變化規(guī)律?；疑Ｐ筒糠?目前使用最廣泛的灰色預(yù)測(cè)模型就是關(guān)于數(shù)列預(yù)測(cè)的一個(gè)變量、一階微分的GM(1,1)模型.它是基于隨機(jī)的原始時(shí)間序列，經(jīng)按時(shí)間累加后所形成的新的時(shí)間序列呈現(xiàn)的規(guī)律可用一階線性微分方程的解來(lái)逼近.經(jīng)證明，經(jīng)一階線性微分方程的解逼近所揭示的原始時(shí)間序列呈指數(shù)變化規(guī)律.因此，當(dāng)原始時(shí)間序列隱含著指數(shù)變化規(guī)律時(shí)，灰色模型GM(1,1)的預(yù)測(cè)是非常成功的.GM(1,1)的定義設(shè)為個(gè)元素的數(shù)列，的AGO生成數(shù)列為，其中.則定義的灰導(dǎo)數(shù)為，令為數(shù)列的緊鄰均值數(shù)列，即，則.于是定義GM(1,1)灰微分方程模型為，即(a)其中稱為灰導(dǎo)數(shù)，稱為發(fā)展系數(shù)，稱為白化背景值，稱為灰作用量.將時(shí)刻代入(a)式中有令，，，稱為數(shù)據(jù)向量，為數(shù)據(jù)矩陣，為參數(shù)向量，則GM(1,1)模型可以表示為矩陣方程.由最小二乘法可以求得GM(1,1)的白化型對(duì)于GM(1,1)的灰微分方程（7），如果將的時(shí)刻視為連續(xù)的變量，則數(shù)列就可以視為時(shí)間的函數(shù)，記為，并讓灰導(dǎo)數(shù)對(duì)應(yīng)于導(dǎo)數(shù)，背景值對(duì)應(yīng)于.于是得到GM(1,1)的灰微分方程對(duì)應(yīng)的白微分方程為稱之為GM(1,1)的白化型.模型建立此預(yù)測(cè)模型是擬合參數(shù)模型，通過(guò)原始數(shù)據(jù)累加生成，得到規(guī)律性較強(qiáng)的序列，用函數(shù)曲線擬合得到預(yù)測(cè)值.建立過(guò)程如下：設(shè)原始數(shù)據(jù)序列有n個(gè)觀察值，，通過(guò)累加生成新序列，利用新生成的序列擬合函數(shù)曲線.利用擬合出的函數(shù)求出新生序列的預(yù)測(cè)值序列.利用累減還原，得到灰色預(yù)測(cè)值序列（共n+m個(gè)，m個(gè)未來(lái)預(yù)測(cè)值）.將序列分為和，其中反映的確定性增長(zhǎng)趨勢(shì)，反映的平穩(wěn)周期變化趨勢(shì).對(duì)序列的確定增長(zhǎng)趨勢(shì)進(jìn)行預(yù)測(cè).模型求解整理得深圳市2001年~2010年常住人口數(shù)，見(jiàn)下表.表一：深圳市2001~2010年年末常住人口數(shù)根據(jù)上述數(shù)據(jù)建立含有10個(gè)觀察值的原始數(shù)據(jù)序列：使用Matlab軟件對(duì)進(jìn)行一次累加，得到新數(shù)列，見(jiàn)表二.表二：GM(1,1)算法擬合值及誤差序號(hào)年份模型值殘差相對(duì)誤差級(jí)比偏差2001724.57002002739.770.00920.123%-0.01212003771.490.00870.112%-0.00042004804.56-0.00470.059%-0.01352005839.06-0.01340.1650%0.00892006875.03-0.00450.0518%0.00902007912.54-0.00020.0020%0.00432008951.660.00270.0287%0.00292009992.460.00260.0257%-0.000220101035.010.00210.0204%-0.0005擬合函數(shù)：由殘差、相對(duì)誤差、級(jí)比偏差可知此模型精度較高，可用于預(yù)測(cè).預(yù)測(cè)值見(jiàn)下表.表三：2011年~2020年深圳市常住人口預(yù)測(cè)人數(shù)同理，整理得2001年~2010年非戶籍人口數(shù)如下表.表四：2001年~2010年非戶籍人口數(shù)表五：GM(1,1)算法擬合值及誤差表六：2011~2020年非戶籍人口預(yù)測(cè)值擬合函數(shù)：一元線性擬合部分由表一作出年份-常住人口數(shù)（單位：萬(wàn)人）曲線如下圖：圖一由圖可見(jiàn)數(shù)據(jù)的線性關(guān)系很強(qiáng)，且在一段時(shí)間內(nèi)仍保持線性增長(zhǎng)趨勢(shì)。由最小二乘法得擬合函數(shù)為，相關(guān)系數(shù)為。非常住人口擬合函數(shù)為。設(shè)為人口結(jié)構(gòu)變化率，則。使用EXCEL作圖如下：2001200220032004200520062007200820092010年份%2001200220032004200520062007200820092010年份%圖二由圖可見(jiàn)，深圳市非戶籍人口所占比例逐年下降，并在2009年左右達(dá)到穩(wěn)定值，約為76%。問(wèn)題二：未來(lái)全市和各區(qū)醫(yī)療床位需求。模型建立的前提：未來(lái)十年深圳市各年齡段比例不變。原因：由上題結(jié)果可知，深圳市非戶籍人口所占比例逐年下降，并趨近于76%的穩(wěn)定值，且深圳市非戶籍人口以年輕人為主，在較短時(shí)間內(nèi)年齡結(jié)構(gòu)基本不變；而常住人口受歷史影響較大，年齡結(jié)構(gòu)在短期內(nèi)基本不變。因此，可以認(rèn)為未來(lái)十年深圳市各年齡段比例保持不變。查找資料得病床需求量公式如下：[3]設(shè)病床總需求為Q，為第年（以2011年為第1年）深圳人口總數(shù)，為年齡段所占比例（0~4歲為1，5~14歲為2，依此類(lèi)推，），為年齡段住院率，平均住院天數(shù)為，平均年床開(kāi)放日數(shù)為。則(b)各年齡段住院率見(jiàn)下圖：‰‰圖三由于自然規(guī)律，各年齡段發(fā)病住院率基本不變，因此可用圖中數(shù)據(jù)計(jì)算。深圳市各年齡段比例見(jiàn)下圖：圖四查找資料得，。由(b)式結(jié)合圖三、圖四可得預(yù)測(cè)結(jié)果如下表。模型的進(jìn)一步處理：醫(yī)院在實(shí)際配置病床時(shí)，除考慮需求量外，還需考慮潛在需求量和流動(dòng)人口對(duì)病床的需求量。由文獻(xiàn)可知城市病床潛在需求量為實(shí)際利用量的20%；根據(jù)第一次國(guó)家衛(wèi)生服務(wù)調(diào)查和流動(dòng)人口調(diào)查結(jié)構(gòu)估算，流動(dòng)人口醫(yī)院病床需要數(shù)相當(dāng)于本地居民需要的10%?？紤]上述因素后，醫(yī)院病床配置的推薦公式為：[4]設(shè)醫(yī)院病床推薦配置數(shù)為，病床潛在需求量為，流動(dòng)人口病床需要數(shù)為。因此(c)根據(jù)(b)可求出未來(lái)十年醫(yī)院病床推薦配置數(shù)如下表。表七：未來(lái)十年醫(yī)院病床推薦配置數(shù)根據(jù)各區(qū)人數(shù)與全市人口的比例可求出各區(qū)病床推薦配置數(shù)。圖五表八：各區(qū)推薦床位數(shù)（以2020年為例）問(wèn)題=2\*ROMANII：不同疾病在不同醫(yī)療機(jī)構(gòu)的床位需求——以肺癌和胃癌為例考慮到模型的實(shí)用情況，我們選取了在人群中發(fā)病率較高的肺癌和胃癌作為預(yù)測(cè)對(duì)象。查閱文獻(xiàn)可得不同年齡段肺癌和胃癌發(fā)病率如下表：表九：不同年齡段肺癌發(fā)病率分組（歲）調(diào)查人數(shù)患病人數(shù)患病率（%）44歲以下82100.0045~54109110.0955~64174390.5265~741857241.2975以上18942.12表十“不同年齡段胃癌發(fā)病率分組 （歲）調(diào)查人數(shù)患病人數(shù)患病率（%）44歲以下8210045~5410910055~64174330.1765~741857140.7575以上18963.18針對(duì)以上數(shù)據(jù)，由第一問(wèn)得出深圳未來(lái)人口變化情況為，設(shè)為肺癌患病人數(shù)，為胃癌患病人數(shù)。因?yàn)槟撤N疾病的患病人數(shù)等于該所有年齡段的人數(shù)與它們所對(duì)應(yīng)的患病率乘積之和，則由圖4的人口年齡結(jié)構(gòu)得每種疾病的患病人數(shù)如下：（單位：萬(wàn)）(單位：萬(wàn)）經(jīng)查閱資料，不同疾病在不同類(lèi)型的醫(yī)院的住院天數(shù)如表十一所示：表十一.肺癌和胃癌在不同級(jí)別醫(yī)院住院天數(shù)[2]病種市級(jí)區(qū)級(jí)住院天數(shù)術(shù)后住院天數(shù)住院天數(shù)術(shù)后住院天數(shù)肺癌36204441胃癌30204735則由上表可知，肺癌在市級(jí)醫(yī)院平均住院天數(shù)為56天，在區(qū)級(jí)醫(yī)院平均住院天數(shù)為85天；胃癌在市級(jí)醫(yī)院平均住院天數(shù)為50天，在區(qū)級(jí)醫(yī)院平均住院天數(shù)為82天。設(shè)肺癌在市級(jí)醫(yī)院的病床需求為，在區(qū)級(jí)醫(yī)院的病床需求為，則由第一問(wèn)得知，理想的病床需求為：[3]在本模型中假設(shè)得肺癌與胃癌的人都住院，則(注：查閱資料知大城市年床開(kāi)放日數(shù)為317天，區(qū)級(jí)則為237天，C為深圳市得肺癌的人數(shù)），。則，（注：為年份，結(jié)果取整）設(shè)胃癌在市級(jí)醫(yī)院的病床需求為，在區(qū)級(jí)醫(yī)院的病床需求為，則同理可得，，，則，(注：為年份，結(jié)果取整）模型評(píng)價(jià)與改進(jìn)方向模型建造中沒(méi)有考慮到年齡結(jié)構(gòu)的變化，而深圳市年齡結(jié)構(gòu)是趨向老齡化方向的，老年人體質(zhì)弱，易生病住院，因此實(shí)際得出的床位數(shù)小于實(shí)際需求的床位數(shù)。模型沒(méi)有體現(xiàn)出深圳不同區(qū)的人口變化函數(shù)，而只是假設(shè)其人口總量對(duì)深圳總?cè)丝诹勘３植蛔?，?shí)際上每個(gè)區(qū)的發(fā)展不同，人口變化也不同。本模型只適應(yīng)于預(yù)測(cè)短期，不適應(yīng)于長(zhǎng)期。本模型還是比較合理的，相比其它同類(lèi)模型而言，本模型立意新穎，結(jié)構(gòu)簡(jiǎn)單，易于理解，同時(shí)具有很強(qiáng)的操作性，值得政府部門(mén)參考。本模型巧妙的忽略了問(wèn)題的次要矛盾，即年齡結(jié)構(gòu)，各區(qū)所占比例，雖然喪失了一些精確度，但是相對(duì)于其它繁瑣的模型而言，還是利大于弊的從提高準(zhǔn)確性與普遍性的角度考慮，還可以引入關(guān)于年齡結(jié)構(gòu)的擬合函數(shù)以及各區(qū)獨(dú)立的人口發(fā)展函數(shù)。參考文獻(xiàn)唐麗芳、賈冬青、孟慶鵬，2008.《用MATLAB實(shí)現(xiàn)灰色預(yù)測(cè)GM(1,1)模型》.滄州師范?？茖W(xué)校學(xué)報(bào)，24:35上海醫(yī)科大學(xué)陳興寶、鄭恩群、陳潔上海市不同級(jí)別醫(yī)院平均住院天數(shù)分析，饒克勤、陳育德關(guān)于制定衛(wèi)生資源配置標(biāo)準(zhǔn)的幾點(diǎn)建議，1999年03期39頁(yè)饒克勤、陳育德關(guān)于制定衛(wèi)生資源配置標(biāo)準(zhǔn)的幾點(diǎn)建議，1999年03期40頁(yè)八、附件GM(1,1)灰色模型常住人口：clearx0=[724.57746.62778.27800.8827.75871.1912.37954.28995.011037.2]n=length(x0);lamda=x0(1:n-1)./x0(2:n)range=minmax(lamda);x1=cumsum(x0);fori=2:nz(i)=0.5*(x1(i)+x1(i-1));endB=[-z(2:n)',ones(n-1,1)];Y=x0(2:n)';u=B\Yx=dsolve('Dx+a*x=b','x(0)=x0');x=subs(x,{'a','b','x0'},{u(1),u(2),x1(1)});yuce1=subs(x,'t',[0:n-1]);digits(6),y=vpa(x)yuce=[x0(1),diff(yuce1)]epsilon=(x0-yuce)./x0%求殘差delta=abs(epsilon./x0)%求相對(duì)誤差rho=1-(1-0.5*u(1))/(1+0.5*u(1))*lamda%求級(jí)比誤差yuce1=subs(x,'t',[0:19]);%共20年預(yù)測(cè)值非戶籍人口：x0=[592.53607.17627.34635.67645.82674.27699.99726.21753.56786.17];n=length(x0);lamda=x0(1:n-1)./x0(2:n)range=minmax(lamda);x1=cumsum(x0);fori=2:nz(i)=0.5*(x1(i)+x1(i-1));endB=[-z(2:n)',ones(n-1,1)