当前位置：文档之家› 2012年全国大学生数学建模竞赛江西科技学院朱文强洪振晖庄朝舜论文

2012年全国大学生数学建模竞赛江西科技学院朱文强洪振晖庄朝舜论文

题目： C

参赛队员：洪振晖朱文强庄朝舜指导教师：教练组

单位：江西科技学院

承诺书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。

我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： C

我们的参赛报名号为（如果赛区设置报名号的话）：

所属学校（请填写完整的全名）：江西科技学院

参赛队员(打印并签名) ：1. 洪振晖

2. 朱文强

3. 庄朝舜

指导教师或指导教师组负责人(打印并签名)：

日期： 2012 年 09 月 10日赛区评阅编号（由赛区组委会评阅前进行编号）：

编号专用页

赛区评阅编号（由赛区组委会评阅前进行编号）：

全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：

脑卒中发病环境因素分析及干预模型

摘要

本文主要对来源于中国某城市各家医院提供的脑卒病人基本信息以及相应期间当地的逐日气象资料进行分析并统计，通过考虑当地环境因素气压、温度和相对湿度对该疾病发病率，建立合理的数学模型，能很好地掌握疾病发病率的规律。

针对问题一，主要是运用excel，matlab等软件对题目所提供的四年期间脑卒中发病病例相关信息进行编程分析与统计并汇成三维散点图，从而分析图象得出结论：1. 年龄是影响脑卒中发病率最重要且独立的危险因素之一，且脑卒中发病率随会年龄的增加而增高。2. 男性脑卒中发病率明显高于女性脑卒中发病率。3.脑卒中发病率也受职业影响。

针对问题二，根据问题一中分析结果，脑卒发病率与环境因素之间存在密切的关系。利用四年的数据进行数据拟合模型分别得到脑卒发病率与气压、温度和相对湿度之间一元关系。利用气压、温度和相对湿度已知条件下，通过建立灰色模型来预测发病率，与真实值进行比较，可知预测的结果的精确度是比较准确的。再通过对灰色模型的改进建立了多元线性回归模型，通过四年的平均发病率、气压、温度和相对湿度建立了一个多元的关系，通过对模型的检验，建立的发病率与气压、温度和相对湿度多元线性关系比较准确。

针对问题三中对高危人群患病特征的分析，得出脑卒中主要是因为心血管疾病及各种不良因素长期作用于脑血管导致而成的。通过对脑卒中发病的关键指标的查找与确定，结合问题一和问题二所得结论，对脑卒中高危人群以不同性别、不同职业状况等条件以及天气状况和脑卒中的关系，得出男性发病率高于女性、各职业人群发病率区别的重要原因、严寒酷暑温差大对发病率的重要影响，从而进行预警和干预，最终向其提供合适的建议方案。

关键词：脑卒中、matlab、数据拟合、灰色模型、多元线性回归、预警和干预

1.问题重述

1.1背景

脑卒中（俗称脑中风）是目前威胁人类生命的严重疾病之一，它的发生是一个漫长的过程，一旦得病就很难逆转。这种疾病的诱发已经被证实与环境因素，包括气温和湿度之间存在密切的关系。对脑卒中的发病环境因素进行分析，其目的是为了进行疾病的风险评估，对脑卒中高危人群能够及时采取干预措施，也让尚未得病的健康人，或者亚健康人了解自己得脑卒中风险程度，进行自我保护。同时，通过数据模型的建立，掌握疾病发病率的规律，对于卫生行政部门和医疗机构合理调配医务力量、改善就诊治疗环境、配置床位和医疗药物等都具有实际的指导意义。

1.2相关情况

数据（见Appendix-C1）来源于中国某城市各家医院2007年1月至2010年12月的脑卒中发病病例信息以及相应期间当地的逐日气象资料（Appendix-C2）。请你们根据题目提供的数据，回答以下问题：

1.3问题的提出

问题一：根据病人基本信息，对发病人群进行统计描述。

问题二：建立数学模型研究脑卒中发病率与气温、气压、相对湿度间的关系。

问题三：查阅和收集文献中有关脑卒中高位人群的重要特征和关键指标结合1、2中所得结论，对高危人群提出预警和干预的建议方案。

2.符号说明

Y：代表发病率；

T：代表温度；

P：代表气压；

H：代表相对湿度；

3.模型假设

（1）假设在条件允许的情况下，题目所提供的所有数据表中存在着的一些不同的时间格式或者错误以及部分缺失和不符合现实的数据均暂不纳入考虑的范围

内。

（2）假设某城市各家医院统计的数据无较大误差，不会影响预测结果。

（3）假设2007-2010年这四年间中国某城市的气候条件正常。

（4）假设每年每月气压、气温和相对湿度有着相似的变化。

4.问题一的解答

4.1问题一分析

题目的附件（Appendix-C1）中给出了与性别、年龄和职业有关联的人群发病的病例信息表，根据题目要求我们应用了MATLAB绘制了2007-2010年的人群脑卒中发病病例信息图，以此来进行统计与描述

首先，以X轴表示性别、Y轴表示年龄、Z轴表示职业，作出了关于这三方面的三维散点图进行观测。由于题目所提供的所有数据表中存在的一些格式上或者数据上的误差以及部分数据的缺失，故在数据受损方面我们暂不将其纳入考虑范围之内。因此，我们作出了以下的脑卒中发病病例信息图。

4.2 问题一模型的建立与求解：对问题一分析后，将附件(Appendix-C1)中数据文件data1.xls -data4.xls的表格，按年份重新分类，依次为2007-2010年的，分别重命名为data5-data8,再将每年表格的性别，年龄，职业用MATLAB绘成三维图，如下所示：

图1 2007年脑卒中发病病例信息图图2 2008年脑卒中发病病例信息图

图3 2009年脑卒中发病病例信息图

图4 2010年脑卒中发病病例信息图

上述图中表明，假设X轴上面的1代表的是男性，2代表的是女性。Y轴代表的是

年龄从0-100岁的各阶段人群。Z轴上的数字代表的是各职业的代号。由图可知：年龄是影响脑卒中发病率最重要且独立的危险因素之一。脑卒中发病人群大多数趋向于中老年人身上，青少年发病率则相对较少，由此可以得出结论:脑卒中发病率，随年龄的增加而越发而增高。

同时从以上四张图中也可以看出，Z5（渔民）、Z6（医务人员）在这四年中的脑卒中发病率明显低于其它各职业的。这是由于Z5人群长期在海面上生存，大胆勇敢地在海上作业，对于恶劣天气的适应能力强于其它各职业人群，另外，海上渔民喜好饮酒，最主要的原因是在炎夏烈日当空时饮酒可解暑，严冬霜雪时饮酒能抗寒，因此对于脑卒中这种与天气密切相关的疾病（问题二会给出具体的分析方案）对渔民的危害明显小于其它职业人群。Z6人群则是因为长期在医院、卫生所等地工作，对于脑卒中的预防以及防止相较于其它职业人群更胜一筹，因此发病率也相对较低。Z1（农民）、Z2（工人）、Z3（退休人员）在这四年中的发病率相对于其它职业人群高出许多，这是由于Z1、Z2都是处于社会的最底层，随着社会的不断发展，来自生活各方面的压力越来越巨大，使得他们不得不面对繁琐沉重的高强度工作，日出而作日落而息，严寒酷暑劳碌奔波，这就给脑卒中疾病的发生提供了巨大的可能性。Z3人群均为退休人员，即年龄均在60岁以上，故发病率高。从表中散点图的分布看来，Z6（医务人员）、Z7（职工）、Z4(教师)这三部分职业人群在脑卒中发病率男性明显大于女性，而其他职业人群男女发病比例大致相同。由上面的四张图中我们可以大致对比出，男性脑卒中发病率明显高于女性。

5.问题二的解答

5.1问题二分析

在问题一的基础上，从Appendix-C2中所提供的2007年到2010年这四年间的当地的逐日气象资料，可以预测到脑卒中与气温、气压和相对湿度具有相关性。绘制发病人数、平均气压、温度和相对湿度的对比图，可以很明显反映出它们之间存在着很直接关系。

为了建立脑卒中发病率与气温、气压、相对湿度间的关系，利用2007年到2010年数据，通过数据拟合建模可以得出发病率与气压、发病率与温度、发病率与相对湿度的函数关系。这样当已知气压、温度和相对湿度任何一个数据可以预测近似的发病率。

但由于发病率与气压、温度和相对湿度存在密切的关系。在综合考虑气压、温度和相对湿度情况下，利用灰色模型对发病率进行预测，并与实际值进行对比与残差检验。观测发病率是否可以预测。

如果可以预测那么可以通过多元线性回归确定发病率与气压、温度和相对湿度之间的一个近似多元线性函数，使得在已知气压、温度和相对湿度情况下，能更加近似地得到发病率，因而更好地掌握疾病发病率的规律。

5.2 问题二的模型建立与求解

5.2.1 建立数据对比图

根据题目所给的Appenddix-C2以及问题一中的统计出来的每一天的发病人数做出了2007年和2008年每一天的发病人数与平均气温、气压和相对湿度的对比图。

图5 2007年发病人数、平均气压、温度和相对湿度的对比图

图6 2008年发病人数、平均气压、温度和相对湿度的对比图从以上两图观察可知，平均气压与发病人数仅有较为细微的关联，即在2007年和2008年中，平均气压略微成曲线的状态延伸，与发病人数在这两年中的曲线状态延伸大

致一致，由此可预测，平均气压与发病人数存有较为细微的关联。平均温度与发病人数大致呈负相关关系，可以看出在寒冷的冬季发病率最高，春季气候转暖，脑卒中发病率也开始下降，夏季发病率达到最低，秋季后又开始上升。平均相对湿度与发病人数成正相关。综上所述，气温和相对湿度是影响脑卒中发病的重要因素。

由于气温对脑卒的发病率有着密切的关系，从每年以月进行对发病率与气温进行对比。

(a) (b)

图7 2007-2010年四年的温差影响发病率（以月计图）根据以上四张温差图我们可以明显的看出，入秋及早春时节也易发生脑卒中。因为此时节最易受寒流侵袭、寒流袭击、气温骤降、日内最高气温与最低气温之差太大。对于刚刚度过的炎热夏季或刚刚感受到春天融融暖意的人们，突来的气候变化会使之不适，对老年人和具备脑卒中高危因素的人来说会是一个噩梦的时季。

同时根据以上四张图易得出结论：温差与脑卒中的发病率大致呈正相关，即温差变化越大，脑卒中发病率越高；温差变化越小，脑卒中发病率也随之降低。

5.2.2 问题二模型的建立与分析

根据发病人数、平均气压、温度和相对湿度的数据对比图，可知它们之间存在密切的关系。为了建立脑卒中发病率与气压、气温和相对湿度间的关系模型，由于每年每月气压、温度和相对湿度变化相似，可以利用2007年到2010年四年的气压、气温和相对湿度进行数据的拟合[1]，拟合出近似的函数关系。

将2007-2010年的数据汇总，提取出四年48个月的月平均气温的数据点，由MATLAB 编程可得，脑卒中发病率与气温之间存在着密切的关系。从而通过数据拟合得出发病率与气温之间的近似函数：

()

432

=-+-+÷

0.0020.01540.3726 3.667519.6530100

Y T T T T

图8 2007-2010年每月的平均温度与发病率的关系图由图中可得，四年中每个月的平均温度在10℃~25℃的时候每月的平均发病率较为稳定，大致为7%左右。当每月平均温度偏低（＜10℃）或偏高（＞25℃）的时候，每月平均发病率明显变高。

同样将2007-2010年的数据汇总，提取出四年48个月的平均气压的数据点，由MATLAB编程可得，脑卒中发病率与气压之间存在着一定的关系。从而通过数据拟合得出气压与发病率之间近似函数：

()

Y P

=-÷

0.082576.1797100

图9 2007-2010年每月的平均气压与发病率的关系图由图中可得，四年中每个月的平均气压与发病率呈正相关关系，即气压越低，发病率越低；反之，气压越高，脑卒中发病率就越高。整体呈平稳向上缓慢递增的趋势。

同样将2007-2010年的数据汇总，提取出四年48个月的平均相对湿度的是据点，由MA TLAB 编程可得，脑卒中发病率与相对湿度之间存在着密切的关系，从而通过数据拟合得出发病率与相对湿度之间近似函数：

()320.00210.436730.3765708.9656100Y H H H =-+-+÷

图10 2007-2010年每月的平均相对湿度与发病率的关系图

由图可得，四年中每个月的平均相对湿度与发病率呈负相关关系，即相对湿度越高，脑卒中发病率越低；反之，相对湿度越低，脑卒中发病率越高。其中，当四年中每月平均相对湿度小于60%的时候，脑卒中的发病率高于8%；当平均相对湿度大于60%小于75%的时候，脑卒中的发病率则大致稳定在8%；当相对湿度小于75%的时候，脑卒中发病率低于8%。

5.2.3 问题二建立灰色模型[2]

但由于发病率与气压、温度和相对湿度存在密切的关系。在综合考虑气压、温度和相对湿度情况下，利用灰色模型对发病率进行了预测。

把2007年至2010年每月的发病率看成一个数列)0(x

，然后建立)1,1(GM 模型：令 ))(),2(),1(()0()0()0()0(n x x x

x ， = 作1-A G O ， ∑==k m m x k x

1)0()1()()( 有

))()1(,),2()1(),1(())(),2(),1(()0()1()0()1()0()1()1()1()0(n x n x x x x n x x x x +-+== ， )1(x 可以建立下述白化形式的方程：

b ax dt

dx =+)1()

这是一阶一个变量的微分方程模型。

记参数列为∧a ，则??

???=∧b a a ，

按最小二乘法解，求∧a ，

N T T y B B B a 1)

(-∧

= ?????????

????

???+--+-+-=1)()1((211)3()2((21

1)2()1((21

)1()1()1()1()1()1(n x n x x x x x B T

N n x x x y ))(,),3(),2(()0()0()0( = 白化形式微分方程的解为：

(1)(0)(1)((1))ak b

x k x e a a ∧

-+=-+

还原数列：

)1()()()

1()1()0(--=∧∧∧k x k x k x

则可以通过灰色模型预测每年的发病率，利用预测与真实值进行对比与残差检验。建立GM （1，1）模型，得2007年发病率预测值：

52.1897752.19865)1(0433.0-=+-k e k x

表1 2007年每个月的模型检验结果

2007年月份实际值预测值

残差检验 1月 888 888.00

0 2月 719 879

-0.22253 3月 1004 918

0.085657 4月 1054 959

0.090133 5月 1043 1001

0.040268 6月 1012 1046

-0.0336 7月 1000 1092

-0.092 8月 1175 1140

0.029787 9月 1210 1191

0.015702 10月 1356 1243

0.083333 11月 1198 1298

-0.08347 12月 1349 1356

-0.00519

建立GM （1，1）模型，得2008年发病率预测值：

??????=??????=∧1.20800565.0b a a 93.3681593.35022)1(0565.0+-=+k e

k x 表2 2008年每个月的模型检验结果

??????-=??????=∧7267.8210433.0b a a

2008年月份实际值预测值

残差检验 1月 1793 1793

0 2月 1925 1924

0.000519 3月 1891 1818

0.038604 4月 1718 1718

0 5月 1707 1624

0.048623 6月 1387 1535

-0.10671 7月 1327 1451

-0.09344 8月 1336 1371

-0.0262 9月 1242 1296

-0.04348 10月 1401 1225

0.125625 11月 1282 1157

0.097504 12月 1005 1094

-0.08856

建立GM （1，1）模型，得2009年发病率预测值：

??????=??????=∧62.8960196.0b a a 03.4574603.44883)1(0196.0+-=+k e k x

表3 2009年每个月的模型检验结果

2009年月份实际值预测值

残差检验 1月 863 863

0 2月 828 871

-0.05193 3月 793 854

-0.07692 4月 851 838

0.015276 5月 860 822

0.044186 6月 774 806

-0.04134 7月 874 790

0.09611 8月 905 775

0.143646 9月 774 760

0.018088 10月 748 745

0.004011 11月 620 731

-0.17903 12月 679 716

-0.05449

建立GM （1，1）模型，得2010年发病率预测值：

??????=??????=∧6.145830164.0b a a

02.8893902.87515)1(0164.0+-=+k e

k x 表4 2010年每个月的模型检验结果 2010年月份

实际值预测值残差检验 1月

1424 1424 0 2月

1224 1424 -0.1634 3月

1420 1400 0.014085 4月

1355 1377 -0.01624 5月

1509 1355 0.102054

6月

1277 1333 -0.04385 7月

1435 1311 0.086411 8月

1367 1290 0.056328 9月

1298 1269 0.022342 10月

1410 1248 0.114894 11月

1304 1228 0.058282 12月

840 1208 -0.4381

根据四年验差检验:

1、计算原始序列标准差：

1S =2、计算绝对误差序列的标准差

:2S =

3、计算方差比：2

1S C S = 得出：0.0021C =

计算小误差概率：()()(){}

0010.7645P P i S =?-?< 得出：1p =。结论：通过残差分析、后验差检验，发现，模型的预测值与实际值相差不大，说明发病率在气压、温度和相对湿度已知条件下，发病率可以被预测，可以看出模型预测的结果的精确度是挺高的。

5.3 问题二的改进模型多元线性回归[3]

一元线性回归是一个主要影响因素作为自变量来解释因变量的变化，在现实问题研究中，因变量的变化往往受几个重要因素的影响，此时就需要用两个或两个以上的影响因素作为自变量来解释因变量的变化，这就是多元回归亦称多重回归。由于发病率与气压、温度和相对湿度存在密切的关系，为了更好地掌握疾病发病率的规律，对单一变量的数据拟合进行改进，考虑多变量的下对发病率的影响。

设y 为因变量，12,,k x x x 为自变量，并且自变量与因变量之间为线性关系时，则

多元线性回归模型为：

01122k k y b b x b x b x e

=+++++ 其中，0b 为常数项，12,,k b b b 为回归系数，1b 为23,k x x x 固定时，1x 每增加一个单位对y 的效应，即1x 对y 的偏回归系数；同理2b 为12,,,k x x x 固定时，2x 每增加一个

单位对y 的效应，即，2x 对y 的偏回归系数等。

多元性回归模型的参数估计，同一元线性回归方程一样，也是在要求误差平方和(2e ∑)为最小的前提下，用最小二乘法求解参数。以二线性回归模型为例，求解回归参数的标准方程组为

01122210111212

220211222y nb b x b x x y b x b x b x x x y b x b x x b x ?=++?=++??=++?∑∑∑∑∑∑∑∑∑∑∑

解此方程可求得02b b b 、

1、的数值。亦可用下列矩阵法求得 21()()e

b x x x y -''=?=∑ 即

0b ?