深圳杯数学建模A题答案
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
摘要
作为中国经济发展的重点城市,人口与医疗问题已经成为我们的焦点话题,是一个复杂的系统工程。
本文针对地区人口年龄分布情况,外来务工人员的数量,从实际出发,在基于一些合理简化假设的基础上,建立数学模型,并充分利用matlab 等软件简化计算,对相关问题进行了有针对性的求解。
在预测未来十年常住人口时,我们运用了matlab 一元线性回归对近十年的数据进行了多次拟合,并对这些拟合进行了比较得出常住人口模型公式为:2() 1.00050.00838.1671Q x e x x =+-+,通过拟合预测出了未来十年市常住人口的数量,同时在网上2000年到2010年的人口结构的数据,通过Leslie 矩阵预测出了未来十年人口结构的分布。
通过分析近人口数量和人口结构的变化,预测未来十年市人口数量和结构的发展趋势,以此为基础预测未来全市和各区医疗床位需求呈线性递增趋势。
同时选取了高血压,脑出血,癌症这三种疾病进行预测,运用matlab 最小二乘法散点拟合,得出这三种疾病的发展趋势,由此预测出未来十年这三种疾病的就医的床位需求。
关键词:matlab 、一元线性回归、Leslie 、最小二乘法、床位需求
一、问题重述
从的人口的结构来看,显著的特点是流动人口远远超过户籍人口,且年轻人口占主绝对优势。
流动人口主要从事第二、三产业的企业一线工人等。
年轻人身
体好,发病少,导致目前人均医疗设施低于全国类似城市平均水平,但仍能满足现有人口的就医需求。
然而,政策的调整与世界的推移会使市老年人增加。
产业结构的变化也会影流动人口的数量。
直接会导致市未来的医疗需求的变化。
现有人口社会发展模型在面对情况时,难以满足人口和医疗预测的要求。
为了解决此问题,请根据人口发展变化态势以及全社会医疗卫生资源投入情况(医疗设施、医护人员结构等方面)收集数据、建立针对具体情况的数学模型,预测未来的人口增长和医疗需求,解决下面几个问题:
1.分析近十年常住人口、非常住人口变化特征,预测未来十年市人口数量和结构的发展趋势,以此为基础预测未来全市和各区医疗床位需求;
2.根据市人口的年龄结构和患病情况及所收集的数据,对几种病进行预测,在不同类型的医疗机构就医的床位需求。
二、问题分析
2.1背景分析
作为我国的经济重镇,经济迅猛发展,带动人口发生了极大变化,大量的人才需求使外来人口大量增加。
劳动力的需求使年轻人占据的的主要地位。
年轻人身体健壮,发病较少,弥补了医疗稍差的缺陷。
然而,由于政府的各项政策(如计划生育等)使得人口结构发生了变化,市统计局12日公布了全市第六次全国人口普查主要数据,显示特区在2000年至2010年的10年中人口增长率近“50%”,人口密度大幅提高。
政府部门需要更详细的人口数量与人口结构的发展趋势,以此为基础来满足市各区几种病的床位需求。
近些年来,对人口结构的分析预测仅限于粗线条分析,
只能预测年龄与性别的大致分布围。
随着人们对健康要求的提高,床位的需求逐渐受到重视,这就是人口与医疗需求的预测。
2.2 问题的分析
题目中所给的两个问题都属于预测的数学问题。
其中问题一需要通过对人口数量极其人口结构进行预测,以此为基础预测未来全市和各区医疗床位需求。
为了解决此问题,我们首先要对近十年的常住人口与非常住人口进行分析,其次再对人口数量和结构进行分析,通过对这些已知数据的分析和统计,在预测未来十年常住人口时,我们运用了matlab 对近十年的数据进行了多次拟合,并对这些拟合进行了比较得出常住人口模型公式为:2() 1.00050.00838.1671Q x e x x =+-+,通过这个模型对未来十年常住人口进行预测。
接而得出未来十年,即得到市2011年到2020年每年的人口增长率,得出市未来十年的人口数量发展趋势。
通过按照年龄来划分儿童、青壮年、老年三个年龄层,求出三个年龄层的比例模型,通过得出关系函数在计算得出未来十年的结构发展趋势。
通过如下关系:年龄结构和患病率相关,患病率和住院率相关,住院人口数和床位有关,建立数学模型,预测得出未来十年的床位需求数。
对于问题二,要求预测不同类型的医疗机构就医的床位需求,根据问题一中得到的全市人口年龄结构和患病情况,对高血压,癌症,脑出血三种病症在不同类型的医疗机构就医床位需求.按照规模大小划分市的医院类别,再通过各等级医院的床位需求与某种病的患病人数和同一等级医院的数量,可治疗这种病的医院总个数的关系得出不同医疗机构就医床位需求。
三、模型假设
1、假设题目所给的数据真实可靠。
2、假设在政府政策的稳定前提下,生育和死亡率都比较稳定。
3、不考虑战争,瘟疫,大规模流行病对人口的影响。
4、假设市人口为年末常住人口。
5、假设同一年龄段的人死亡率相同,同一年龄段的育龄女性生育率相同。
6、假设当地人们的生育观念不发生太大变化。
7、假设人们生病时都能支付起医疗费。
8、假设各区域的患病者不相互交换,即各区域是相互独立的。
四、定义符号与说明
见文中标注
五、模型的建立与求解
一、问题一的分析
(一)市常住人口的预测
(1)利用现有数据(表一)分析从1979年到2010年的年末常住人口数变化规律。
运用Excel软件画出1979年到2010年的人口数量折线图(图一):
表1:1979—2010年年末常住人口数
图1:1979—2010年年末常住人口数
年末常住人口数
0200
400
600
800
1000
1200
1979
1982
1985
1988
1991
1994
1997
2000
2003
2006
2009
年份
人口数(单位:万人)年末常住人口数
(2)通过现有的数据及其折线图,可以很明显地观察出常住人口数从1980到1992的人口处于缓慢增长,呈线性增长。
但随着高速的发展,优质的社会公共资源对流动人口形成了强大的吸引力,因此外来人口的迁入增多导致从1994年到2010年年末常住人口数的增长率相对以前增大,但也基本保持一次函数的增长。
(3)模型的建立
我们通过运用matlab 软件对这一组数据进行多次拟合,其根本思想就是:观测散点走势来确定拟合函数,利用散点但又不拘泥于散点。
他的整体思路与我们的数据分析非常相似。
并对这些拟合进行了比较得出常住人口模型公式为:2() 1.00050.00838.1671Q x e x x =+-+,拟合结果如下图(图二):
图2:常住人口的拟合结果图
(二)流动人口的预测
从的人口的结构来看,显著的特点是流动人口远远超过户籍人口,因此对流动人口的预测对整个及各区医疗床位需求的预测中起到至关重要的作用。
(1)流动人口定义
流动人口是相对于某地的常住人口而言的, 指离开常住户籍所在地, 跨越一定的行政辖区围, 在某一地区滞留的人口。
其包括:
1、进入城镇务工、经商、和从事劳动服务的暂住人口;
2、为探亲访友、旅游、求学、治病等而外出的人员;
3、无职业、无收入、无暂住证的三无人员即盲流人口。
为此我们可得:123Q Q Q Q =++非
其中:
Q 非——非常住人口总和;
1Q ——进入城镇务工、经商、和从事劳动服务的暂住人口总和;
2Q ——为探亲访友、旅游、求学、治病等而外出的人员;
3Q ——无职业、无收入、无暂住证的三无人员即盲流人口。
(2)求解进入城镇务工、经商、和从事劳动服务的暂住人口:
显然对于1Q ,它是市经济发展主要的带动者,因此与市GDP 有很大的关系,GDP 越多,则市外来人口就越多。
为此我们假设1Q 与外来人口所产生的GDP 成正比例关系,由此我们可得:
1()tGDP Q p X X b =-+1
其中:
p ——比例因素;
tGDP X ——市t 当年GDP 总量;
X ——常住人口GDP 值;
b 1——进入城镇务工、经商、和从事劳动服务的暂住人口总和的初始值;
对于一个非平稳序列来说,其数字特征,如均值,方差和协方差等是随着时间的变化而变化的。
也就是说,非平稳序列在各个时间点上的随机规律是不同的,难以通过序列已知的信息去掌握序列整体上的随机性。
而GDP 时间序列都是非平稳的,为此我们采用ARIMA 模型求解:ARIMA 模型使用包括自回归项(AR 项) , 单整项和MA 移动平均项三种形式对扰动项进行建模分析, 使模型同时综合考虑了预测变量的过去值, 当前值和误差值, 从而有效地提高了模型的预测精度。
(1)ARIMA 模型的形式:
考虑序列t y ,若其能通过d 次差分后变为平稳序列, 即~()t y I d , 则 (1)d d t t t u y B y =∆=-
t u 为平稳序列, 即~(0)t u I , 于是可建立ARIMA (,)p q 模型:
1111t t p t p t t q t q u c u u φφεθεθε----=+++++++
经d 阶差分后的ARIMA (,)p q 模型称为ARIMA (,,)p d q 模型。
其中p 为自回归模型的阶数,q 为移动平均的阶数,t ε为一个白噪声过程。
(2)建立ARIMA 模型的一般方法:
1)检验原序列的平稳性 检验的标准方法是单位根检验, 若序列不满足平稳性条件, 则可通过数学方法, 如差分变换或者对数差分变换使其满足平稳性条件;
2)通过计算能够描述序列特征的一些统计量, 如自相关(ACP )系数和偏自相关
p q模型的阶数p和q,并根据一定的准则, 如ATC (PACP)系数来确定ARIMA(,)
准则或SC准则等综合考虑来确定模型的参数;
3)估计模型的未知参数[2], 并通过参数的 统计量检验其显著性, 以及模型的合理性;
4)进行诊断分析, 检验模型的拟合值和实际值的残差序列是否为一个白噪声序列。
(3)数据的来源与描述:
从《统计年鉴》各卷统计出1979 至2006 年国生产总值, 见表5:
X, 具有长期上升趋势, 非水平平稳。
并按此数据作图1从中可以粗略地看出
t
表2:1979——2006年国生产总值统计表(亿元)
图3 图4
(4)序列的平稳性处理:
X,进行平稳性检验(ADF检验) ,结果如表2 :
对
t
表3:序列A D F 检验结果
由表7可知其不平稳。
为了消除原始数据序列的不平稳性, 使数据更为平稳, 本文采用对国生产
∆,总值序列取对数形式, 记为ln t X,序列ln t X一阶差分后的序列记为ln t X
∆,按二阶差分后数据作序列图2 , 可见时间趋势二阶差分后的序列记为2ln t X
基本消除, 可认为是平稳序列但序列图只能粗略地判断序列具有平稳性, 理论上应用单位根检验方法检验。
∆, 进行平稳性检验(A D F 检验) , 结果如表3 :
对2ln t X
表4:序列ADF检验结果
由表7可知其平稳,说明GDP 序列为2 阶单整序列, 即2
ln ~(2)t X I ∆ 模型的识别与建立
由以上对序列ln ~(2)t X I ∆, 的A D F 检验, 我们可确定(,,)ARIMA p d q , 模型中的d 应取为2为了确定模型中的p 和q , 作出序列2ln t X ∆直至滞后16 阶的自相关(ACP )图和偏自相关(PACP) 图, 分别见图3 和图4.
由图7和图8可看出, 少In Xt 序列的自相关图与偏自相关图都是拖尾的, 因此可建立:
图5 图6
ARIMA 模型。
经反复计算比较, 最终取1p =,
2q =, 建立如下(1,2,2)ARIMA 模型: (括号中的数据为对应估计值的T 检验统计量)
2ln (0.031188,(1)0.19417,(2) 2.087428)t X c AR MA ∆==-==-
( 6.899257)-( 4.005350)-( 4.247829)-
图7 .0.050796S E = 3.009846AIC =- 2.863581SC =-
即:
2212ln 0.0311880.19417ln 2.087428t t t t X X εε--∆=--∆+-
20.842R = 2.23DW =
由模型(1 ) , 对其进行回归拟合, 模型中的残差序列(Residual) 以及过
ln t X ∆的实际值(Actual)和拟合值(Fitted )的序列图见图9:
从图9可以看出, 模型的拟合值和实际值的变动具有较好的一致性。
其次, 模型的残差值较小,消除了线性或者指数趋势, 表现得较为平稳, 说明模型通过了适应性检验, 所以该模型还是比较
理想的。
为了进一步检验该模型的效果, 记ˆt u
为该模型的残差序列, 对其进行DF 检验, 得:1ˆˆ1.118299t t u
u -∆=-,DF 的值为-5.3921而在1%显著水平下,DF 的临界值为-2.6649,因此,残差序列ˆt u
, 即误差项序列能在1 %显著水平下被看作白噪声过程,这说明2
ln t X ∆的拟合值是实际值的无偏估计, 模型具有较好的拟合
效果。
作出残差序列ˆt u
前16 阶的自相关(ACP)和偏自相关(PACP)图, 分别见图10和图11。
从两图我们也可看出, 自相关函数和偏自相关函数均落在置信区间, 残差序列应为白噪声过程, 这与上面D F 检验的结果一致。
图8: 自相关(ACP)图图9:偏自相关(PACP)图 (5)模型的预测: 由(1,2,2)ARIMA 模型得:
2212ln 0.0311880.19417ln 2.087428t t t t X X εε--∆=--∆+-
又因为:2
12ln ln 2ln ln t t t t X X X X --∆=-+ 可得ln t X 的预测公式为:
21212ln 2ln ln 0.0311880.19417ln 2.087428t t t t t t X X X X εε----=---∆+-
因此得序列t X 的预测公式为:
21212
2ln ln 0.0311880.19417ln 2.087428t t t t t X X X t X e
εε-------∆+-=
用(1,2,2)ARIMA 模型对国生产总值作预测, 结果见表4
表5:实际值与ARIMA 模型预测值比较衰(亿元)
为此,我们可以求出p 和b 1的值: 由1()tGDP Q p X X b =-+1可得:
()tGDP Q p X X b =-+出
通过1979年初始可知31.41b =,X 出几乎可以忽略不计,则: 通过上面数据求出p 的平均值为:
=14.871%p
由此可得:
14.871%31.41tGDP Q X =+
则:Q Q Q =-1
常
可得下表:
表6:Q 与时间关系表
3)求解为探亲访友、旅游、求学、治病等而外出的人员:
对于2Q ,探亲访友与市现有人口总数成正比,旅游人数可以通过市旅游人口数情况可直接求的;求学人数同样可以通过市教育机构统计数求解,但考虑到未来市不断在发展,所以求学人数也不断在上升,然而整个国家已经入老龄化社会,而且据国家统计局统计年签表明,我国学生数量在不断下降,这两因素一综合,我们假设外来求学人数为恒定不变的;对于外来治病人数,显然与市公有医院服务水平有很大关系,我们假设成正比关系,因此我们可得:
22++Q Q Q Q Q b ++2医旅学=p
其中:
2p ——探亲访友人数的概率;
Q 旅——旅行人数; Q 学——外来学习人数; Q 医——外来求医人数;
2b ——其他人数。
求解2p ,对于探亲访友人数应该和在该地区中人口成正比,在1979年,刚开放,以此那时没有几乎没有其它外来人员,为此我们可得:
20.15
0.00423631.41
p =
= 可得:
表7:访友人数表
求解Q
旅,根据现有的资料,我们查的市南山区2008年统计年签旅行情况可得:
如下表:
表8:09年6月旅游者接待情况统计
为此,根据上表我们求解出该区每天平均每天接待人数和同比增长率如下表所示:
表9:每天平均每天接待人数和同比增长率
由上表可知近几年来旅游增长幅度不大,而且旅游是非常住人口的一小部分,为了减少计算难度,我们忽略的这种增长。
对于共有7个区,为了简化计算,我们假设旅客到每一区去旅行都是随机的,去每区每年平均每天接待人数为3.425万可得:
7 3.42523.927Q =⨯=旅万
根据资料可得:目前有35万左右的义务教育阶段非户籍学生。
这数字占了义务阶段学生的一半。
同理可以求得:
35Q =学万
对于外来求医人数,市公有医院服务水平有很大关系,我们假设与公有医院的等级成正比,与公有医院的总数成正比关系,因此我们可得:
2=Q p S 医
其中:
2p ——公有医院等级因数;
S ——公有医院总数;
但是根据题意:此目前人均医疗设施虽然低于全国类似城市平均水平,但仍
能满足现有人口的就医需求。
可知,对于医疗水平,相对于其它如、等一些大城市相比,医疗水平很弱,因此为了简化模型,我们将Q 医直接放到其它人口2b 中考虑。
3)求解三无人口数目Q 三无:
三无人口定义:无职业、无收入、无暂住证的三无人员即盲流人口。
由此我们可以得到该三无人口出现的概率非常小,几乎可以忽略不计,为此我们也将他归为其它人口。
综上所述:
23=p()++GDP Q W W b Q Q Q Q b Q -+++++12
医旅学非p =()58.927GDP p W W p Q b -+++2
综上所述求解Q 非可得:
表10:非常住人口表
(三)市未来十年人口结构的预测
根据2010 年人口总数是1037.2万,按照每五岁为一个年龄组,把0~99 岁划分成20 个年龄组,即0~4 岁为第1 个年龄组,5~9 岁为第2 个年龄组,10~14 岁为第3 个年龄组, ⋯,95~99 岁组第20 个年龄组,100 岁以上为第21 个年龄组,并设各年龄组人口构成的初始人口列向量为X(0) = [ x1 (0) ,x2 (0) ,x3
(0) ,⋯,x21 (0) ] T ; 第5t 年各年龄组人口构成的人口列向量为X(t) = [ x1 (t) ,x2 (t) ,x3 (t) , ⋯,x21 (t) ] T ,称X(t) 为人口状态向量。
如果设所有年龄组女性人口占同一组总人口比例的系数向量为C = [c1 ,c2 ,c3 , ⋯,c21 ] T ,那么在5t 年时,女性人口的列向量应为C ·X( t ) = [ c1x1 ( t ) ,c2x2 ( t ) ,c3x3 ( t ) , ⋯,c21x21 (t) ] T 。
各年龄组妇女在五年的平均生育率向量为B = [ b1 ,b1 ,b2 , ⋯,b21 ] T ;由于在2000 年以后,随着独生子女群体结婚高峰的到来,按照我国现行计划生育政策,这一群体允许生育第二胎,因此育龄妇女的生育率将会上升,其上升幅度现在很难准确估计,但总和生育率R 应满足不等式:1 < R < 2 , (即平均一对夫妇终生只能生育R 个孩子) 。
如果2000 年以后按2000年总和生育率(1 125 ‰) 的a (0. 9 < a < 1. 3)倍进行估算,那么可取B = a[ b1 ,b1 ,b2 , ⋯, b21 ] T 。
若把t 阶段存活的全部新生儿划分到第t + 1 阶段的第一年龄组,并设各年龄组人口在五年期的自然存活率向量为S = [ s1 ,s2 ,s3 , ⋯, s21 ] T 。
由于第t 阶段k - 1 年龄组的人存活到第t + 1 阶段就是k 年龄组的人, (k = 2 ,3 ,4 , ⋯,20) ,且第21 年龄组(即100 岁以上) 的老年人五年后存活下来的仍然属于第21 年龄组。
由此可得人口系统状态X( t ) 关于离散时间变量t ( t = 1 ,2 ,3 , ⋯,n ,⋯) 的状态转移方程组
()()
t a t x b c x k
k
k k
∑==+21
1
1
1
()()t t x s x k k k
1
1
1--=+
x 21( t + 1) = s 20x 20 ( t) + s 21x 21( t) (1)
引进系数矩阵:
⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫
⎝
⎛=s s s b c b c b c s s s b c b c b c a a a a a a a 212019
2121202019193
213322110000000000000000000
000000000000000000
则方程组(1) 可用矩阵形式表示成X(t + 1) = AX(t) t = 0 ,1 ,2 ,3 , ⋯(2) 矩阵A 为Leslie 矩阵[2 ] ,以A 为系数矩阵的人口状态向量X(t) 的转移方程(2) ,就是人口增长的动力学模型。
若以2010 年的人口向量为初始向量X(0) ,把X(0) 代入方程(2) 可依次求得2015年、2020年等以后第5t 年的人口向量X(t) 的预测值。
由于方程(2) 以五年为一个时间单位,故应根据表 2 中的数据计算出五年各年龄组的死亡率与生育率。
假设第k 组人口年平均死亡率为λk ,则由于单位时
间
dt 的死亡人数与人口总数()t x k 成正比,即有
)
()
(t dt
t d
x x
k k k
λ-=,解此微分
方程可得五年的人口存活率为
()
21,.....3,2,15==-k e s k k λ。
但当第k 组育龄妇女的
年平均生育率为f k 时,五年的平均生育率就是bk = 5fk (1 ,2 ,3 , ⋯,21) 。
经计算可得以五年为一个单位时间时这两组数据组成的向量分别为
S = [ 0. 970 009 , 0. 997 553 , 0. 997 802 , 0. 996 357 ,0. 994 068 , 0. 993 372 , 0. 992 578 , 0. 991 189 , 0. 987 973 ,0. 982 161 , 0. 971 999 , 0. 955 042 , 0. 924 641 , 0. 872 406 ,0. 774 103 , 0. 672 032 , 0. 505
554 , 0. 396 135 , 0. 275 891 ,0. 313 627 , 0. 301 194 ]T ;
B = [ 0,0,0,0.0034,0.06014,0.09029,0.036,0.0093,0.0018,0.00048,0,0,0,0,0, 0,0,0,0,0,0 ] T ;
C=[0.4520,0.4169,0.4162,0.5829,0.5435,0,4727,0.4542,0.4298,0.4270,0.4 207,0.4268,0.5065,0.5275,0.4967,0.5107,0.5277,0.5672,0,7708,0.7442,0. 9091,1]
经计算得到未来十年的人口人口结构:
表11
由表格,对比2010年的年龄结构,可以看出,市将面临老龄化严重的问题,并且中青年人口也在总人口中的比重降低。
(三)预测未来全市和各区医疗床位需求
假设光明新区和坪山新区是在2010年时新增加的两个区.并设定Y 总为全市医疗床位总需求量,n Q ()1,2,
,8n =为各区医疗床位需求量,不妨令81n
n Y Q ==∑总()1,2,,8n =.
设D 为各区年龄结构比例,W 为全区总人数,R 为全市总床位数(见[附件]),P 总
为全市总人数,Z 为住院率,Q 为区床位需求,则 区床位需求=各区年龄结构⨯全区总人数⨯(全市总床位数/全市总人数),可得Q D W Z =⨯⨯即(/)G D W R P =⨯⨯.
由于所给数据有限,我们只得到了2000年和2010年的各区人数和各区中各个年龄层的人口数量分布,运用matlab 最小二乘法拟合散点,得出2000--2010年各区床位需求大致走向是呈正向发展趋势,如图所示:
图10:2000--2010年各区床位需求曲线图
进而得出了各个区的床位需求量与年份的函数关系式:
罗湖区:189.8t 1140Q =+,
福田区:2156.6t 1337Q =+,
南山区:3134t 1060Q =+,
宝安区:4483.5t 4011Q =+,
龙岗区:5191t 2581Q =+,
盐田区:646.3t 224Q =+.(取t =11,12, ,20)
则全市总床位需求621)(Q Q Q t R +++= ,以此预测出未来十年的各区和全市医疗床位需求。
二、问题二的分析
表13:2005年患病人数表
表14:2010年患病人数表
根据拟合得出:
高血压人口的发展趋势为
1246.6t1083
y=+,
癌症人口的发展趋势为:
21067.36657.2
y t
=+,
脑出血人口的发展趋势为:
3100.5804.5
y t
=+.
按照病情的严重性和人群就医的一般规则,癌症和脑出血应在综合医院和专科医院救治,高血压可在综合医院,专科医院,街道(镇)医院进行治疗,癌症
六、模型的评价
优点:本文在预测问题一中未来十年年末常住人口时,通过matlab的一元线性回归对已知数据进行多次拟合,由此得出了未来十年年末常住人口的函数模型,通过多次的拟合,得到了较高的拟合度。
在预测人口结构时,使用了Leslie 矩阵,该模型最大的优点就是即使已知的数据很少,也可以通过建立模型对未来的数据进行预测。
通过对常住人口数的预测和人口结构的预测,加上题目提供的各区的人口数据,我们运用matlab最小二乘法拟合散点得出2000到2010年各区床位需求大致走向是呈正向发展趋势,而且在应用最小二乘法拟合数据实可以根据需要进行多次数的拟合,5次的精确率高达95%,所以我们的数据误差较小准确性较高。
我们应用了几乎所有数据更能体现出数据的可靠性和真实性,更加准确地预测出未来全市和各区医疗床位需求。
问题二则是根据市人口的年龄结构和患病情况,及所收集的数据,计算预测出高血压、癌症和脑出血在不同类型医疗机构就医的床位需求,具体如下:3种病症的患者人数在未来一段时间大致接近一次函数发展趋势,并且得出各医疗机构床位需求函数:
= 同一等级医院的数量
各级医院床位需求某种病的患病人数
可治疗这种病的医院总个数
.
由于原题中未给出相应的各等级医院数和各医院床位数,因此,我们无法得出具体结果。
缺点:首先我们计算得出的预测数据是基于较多的假设的基础之上得出的,从客观上讲,较为理想化,有很多客观或主观因素不能全面的考虑进去,只能得到在一定程度上更接近正确的合理性数据;其次由于数据有限,未能充分考虑影响人口增长的因素,因此所建模型不全面,所得的结果与实际有一定的出入。
并且由
于有限的时间,我们只能将最关键的情况考虑进去。
改进:在第一个问题的解决中我们没有考虑户籍人口和非户籍人口的关系,直接使用了常住人口来拟合数据。
虽然有很多地方都出现了问题,但是计算方法和计算的结果绝对具有真实性。
该模型在构造时有个很重要的假设,即考虑床位需求和患病人数成正比,但实际情况要复杂的多,要考虑患病高峰期、环境变化、医疗结构和医疗投入等诸多因素的影响,因此我们可以通过修正住院率的值,将住院率(总床位数/总人数)作为变化函数带入到模型中使模型更加准确等等。
七、全文总结
人口预测是一个极其复杂的过程,当中考虑的因素极为复杂,有人口迁移,经济变化,自然灾害等等。
是一个高速发展的城市,经济体制,产业体制,政府政策等等都会影响人口预测,对此造成极大的影响。
本文人口预测是建立在的产业结构,政府政策等等不发生太大变化,允许存在小围的变动,运用一元线性及Leslie矩阵,模拟估计。
并以此为依据对医疗需求进行了相关的预测。
由于作者能力有限,知识水平有限,故模型可能存在问题,希望各位给予指出,进行指导,本人定会对模型进行完善修改。
参考文献:
[1] 卫国.MATLAB程序设计与应用(第二版)[M].:高等教
育.2002.
[2] 启帆.数学建模.[M].:高等教育,2004.
[3] 谭永基,蔡志杰,数学模型,:复旦大学,2005.2.
[4]强,人口系统模型及人口状况分析,中国优秀硕士论文,2004.9.
[5]王晓军,蔡正高,死亡率预测模型的新进展[J].统计研究,2008,(25).。