奥运会奖牌榜预测问题

1、2016年中国矿业大学数学建模模拟竞赛参赛基本信息专用页承诺书我们仔细阅读了大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。参赛学生基本信息题目里约奥运会奖牌榜预测评阅编号（由组委会统一填写）姓名（打印）班级电

2、话本人签名日期：年月日里约奥运会奖牌榜预测摘要本文主要探讨了奥运会奖牌榜的预测问题，通过建立模型对2016年里约奥运会奖牌榜进行预测。本文主要通过建立比较时间序列模型、多远线性模型喝灰色预测模型，使用MATLAB和SPSS软件对历届奥运会奖牌榜进行分析并最终得到最优预测模型。在模型一的建立中，我们假设奥运会奖牌榜之语往届奥运会奖牌榜有关，使用时间序列法来进行预测，通过对数据拟合建立回归模型。由于模型一得预测结果与实际奖牌榜出入较大，在模型二的建立中，我们综合考虑了国家GDP、人口数量、经济发展水平、人口结构和东道主效应等因素建立多元回归模型。关键词：一、问题的重述奥运会是竞技体育顶级盛会，其所

3、获奖牌数及国家排名在一定程度上是衡量一个国家体育水平的尺度，也是一个国家的经济、政治和综合实力的具体体现之一。2016年8月6日早上7:00在巴西的里约热内卢隆重开幕（以下简称里约奥运会），对于正在举行的里约奥运会，大家普遍关心的问题就是奖牌榜的排名，奥运会奖牌榜成了大家关心的热点问题。现请查阅资料，并根据以往各国奖牌榜排名情况，以及各国经济发展、人口体质、政府政策等各种能影响到奖牌榜的因素，建立数学模型，预测2016里约奥运会的奖牌榜前十名。并据此对各国体育水平进行分类。二、问题的分析三、模型假设俄罗斯队因为禁赛原因结果取预测值的2/3世界银行和奥委会提供的资料真实可靠今年奥运会共颁发302

4、块金牌，960块奖牌四、符号说明i国家编号t奥运会届数M每个国家的金牌总数y每个国家奖牌总数POP国家的人口总数GDP国民生产总值Host奥运会主办国Sys社会制度Gen青少年在总人口总所占比例五、建立模型及求解5.1、模型一得建立与求解5.1.1、模型一的建立：最简单的增长模型就是马尔萨斯的指数增长模型，设某年的奖牌数为M(),n年后的奖牌数为M()，届增长率为常数k，则有：M(n)二M()(1+k)n(1)虽然该模型可以较好的预当参赛国家较少时的奥运奖牌榜，但在实际中增长率k并不是常数且由于各种原因增长率差别较大，所以此模型无法有效的预测本届奥运会奖牌榜。时间序列模型是描述时间序列统一性的

5、一种常用方法，在预测未来事件时，通过该事件的历史数据揭示事假未来的发展规律，从而对该事件的未来做出预测。时间序列模型强调时间因素在预测中的主体作用，暂不考虑外界因素的影响。所以在模型的建立中我们假设奥运会奖牌榜只与历届奥运会成绩相关，暂不考虑其他因素的影响。因为每一届奥运会比赛项目都在发生变化，奖牌总数也随之变化，为了使任意两届奥运会奖牌数目具有可比性，本文采取奖牌数目的相对值来进行计算，从而消除奖牌数目变化带来影响，即使用Ms=M一来计算。iYMii=1,2,3在模型的建立当中，我们使用线性回归模型来预测2016年奥运会奖牌榜。奥运奖牌总数和时间之间的一元线性回归模型为：y=a+at+u(2

6、)1其中t为时间，a、a为待定参数，u为随机变量，用MATLAB拟合后得到总奖牌数目拟1合曲线。5.1.2、模型一的求解以德国为例进行2016年奥运会奖牌榜排名：表一：1992-2012年德国奥运会奖牌数年份金牌数金牌占比奖牌总数奖牌总数占比19884819.92%14217.91%19923312.69%821.6%199627.94%658.82%2000134.36%566.7%2004144.64%485.16%2008165.3%414.28%2012113.64%444.57%由于1988年德国分为联邦德国和民主德国进行参赛，数据不具有参考价值，拟合时不予考虑。德国历届奥运会奖牌总

7、数拟合结果：图一：德国历届奖牌总数线性回归拟合曲线回归系数回归系数值估计a00.394a1-0.012R2=0.877F=28.649p=0.006届数总奖牌占比011111124252627282930312186Jo-0-0ooO04a20a表二：奖牌总数线性回归结果：德国奖牌总数线性回归结果为：y=0.0121+0.394使用SPSS软件预测德国队在今年奥运会中拿到的总奖牌占比为0.0649。使用同样方法对其他国家进行线性拟合结果如表三所示：、表三：各国奖牌总数占比拟合结果：国家拟合方程美国y=0.003t+0.232中国y=0.108t0.407俄罗斯y=0.0161+0.550英国y

8、=0.0050.123澳大利亚y=0.002t0.032法国y=0.001t+3x10人（-6）日本y=0.003t0.052德国y=0.0121+0.394意大利y=0.001t+0.005韩国y=0.012t+0.394荷兰y=0.00071+0.00108匈牙利y=0.003t+0.019通过对各个国家历届奖牌数目占比进行线性拟合，我们得到了每个国家金牌占比M和总奖牌占比y与奥运会举办届数t的线性关系，榜排名。表四：第31届奥运会奖牌榜预测并由此得到第26界奥运会的奖牌国家奖牌总数奖牌总数排名金牌数金牌数排名美国991421中国972312英国503124澳大利亚464115俄罗斯425

9、183法国376117日本377711德国378106意大利349108韩国25101095.2、模型二的建立与求解5.2.1、模型二的建立：由于在模型一中仅考虑时间因素，线性拟合结果不理想。考虑到是否为东道主对比赛结果特别是金牌数量，具有显著影响，因此在模型二中，采用了灰色GM(1,1)模型与东道主因素相互结合的方法进行预测，假设奖牌榜上各成员成绩只与是否为东道主和在往届的表现有关。客观世界中既有大量已知信息又有大量未知信息，相互交织。如果将已知信息称为白色信息，未知或非确知信息称为黑色信息，既含有已知信息又含有未知或非确知信息的系统称为灰色系统。这些系统仅具有以时间为序列的数据，因此在研究

10、时可将时间序列转化为微分方程，建立抽象系统发展变化的动态模型，这就是灰色DM或GM模型。在模型二中，为排除东道主所产生的显著影响，首先对原始数据进行了处理，将东道主当年奖牌与金牌占有量修正为上下两届的平均数。同时考虑到各国的体育水平更多地由金牌(或奖牌)占有量体现，而不是由纯粹的数量决定，因此为了使任意两届奥运会奖牌数目具有可比性，本文采取奖牌数目的相对值来进行计算，从而消除奖牌数目变化带6)7)i=1i=1i=1来影响，即使用Ms=理一来计算。特别的，为简化运算假设本届将有金牌量为302i工Mii=l,2,3个，奖牌量为960个，因此M二302DMs，y二960QMs。iiiiGM(n,h)

11、模型是指n序列h阶线性动态模型，在模型二中我们采取单序列一阶线性动态模型，即GM(1,1)模型。假设按照时间排列的处理后的金牌榜数据为M(o)=M(o)(k)1k二1,2,3利用累加生成法得到新数列M，其中：3)M(k)=丈M(o)(i)i=1生成的数列M(i)(k)为一阶一个参数的灰色微分方程，记为GM(1,1)，即dMdt+aM=u4)其中a，u为待定常数，t为时间变量。该微分方程的解为:MM(k+1)=M(1)_-aueak+a5)该式的参数列为M=(a,u)t用最小二乘法求出a=(BtB)-1BtHN由灰色GM(1，1)模型，有B=MM(1)(2)+M(1)(3)M(1)(n)H=1N

12、M(0)(2),M(0)(3)M(0)(n)T根据线性代数的相关知识可以求解、化简得到：m=(n-1)艺(b)2-(艺b)2iii=1i=1a=(n一1)艺b匚h艺b艺h一miiii8)(b)2艺hiii=1由微分方程的解的表达式可求得M(k+1)因此第n项的估测值为:M(o)(k+1)=M(i)(k+1)-M(i)(k)(9)通过以上计算得到的拟合值、预测值，与给定数列存在一定误差，称为残差。构造残差数列E(0)(k)=M(0)(k)-MM(0)(k)(10)残差数列可以揭示数据拟合值与实际值的差额，也可利用此序列再次应用GM(1,1)模型进行同样的计算，得到预测值的修正值。将金牌榜奖牌榜数

13、据代入，可得到预测的金牌榜奖牌榜数据及排名。奖牌预测同理。5.2.2、模型二的求解以韩国的金牌数预测为例，对模型二进行求解。首先，对韩国历届获得金牌数量排除东道主因素的影响，并利用韩国每一届的金牌占比得到经过修正的金牌数，如下表所示：表五韩国历届奥运会金牌获得情况年份金牌数金牌占比修正金牌数1988120.04979253113.938461992120.04615384613.93846199670.0277777788.388889200080.0268456388.107383200490.02980132592008130.043046358132012130.04304635813利

14、用式(3)得到待处理数据：数列M(0)(k)、M(1)(k)，如下表所示表六灰色模型待处理数据年份1988199219962000200420082012序号k1234567M(0)(k)13.9384613.938468.3888898.10738391313M(1)(k)13.9384627.8769236.2658144.3731953.3731966.3731979.37319为求出a、u由MM(i)(n)H=M（o）（2）,M（o）（3）M（o）（n）TN两个式子可计算出矩阵B、H中各元素，即h（i）与b（i）的值，具体数值如下表所示表七灰色模型计算过程量123456Sumh(i)1

15、3.948.398.119.0013.0013.0065.43b(i)-20.91-32.07-40.32-48.87-59.87-72.87-274.92b2(i)437.131028.571625.662388.593584.805310.5014375.26b(i)Di(i)-291.42-269.04-326.89-439.86-778.35-947.35-3052.91由上表数据，以及式（6）、（7）、（8）可以求出a、u的值，由式（5）可知所构造的一阶微分方程的解为M（k+1）=M（1）（1）e-ak+_aa将k依次代入,得到韩国各届奥运金牌的预测值，并根据式（10）求出残差。所求

16、数据如下表所示：表八韩国各届奥运金牌预测值及残差年份19881992199620002004200820122016序号k12345678衍(0)(k)13.9410.0810.3910.7211.0511.4011.7612.12M(1)(k)13.9424.0234.4145.1356.1867.5879.3491.46残差E（k）0.00-3.862.012.612.05-1.60-1.24拟合值与原值曲线如下图二韩国各届奥运金牌数与预测金牌数曲线国家金牌金牌名次奖牌奖牌名次中国4611011美国402982英国283593澳大利亚124445韩国1252710日本126406俄罗斯11

17、7444德国108377法国99358意大利9102995.3、模型三的建立与求解：5.2.1、模型三的建立：在实际问题中，不同国家性质与经济发展水平都极大的影响奥运会奖牌榜排名，在模型二的建立中我们引入一国国民经济总收入（GDP）、人口总数（POP）、青少年在总人口中所占比例（Gen）主办国（Host）和政治体系（Sys）等因素进行综合分析预测。基于历届奥运会奖牌榜，以各个国家奖牌榜占比作为因变量，各影响因素作为自变量，其中国家体制和东道主效应是虚拟变量，建立线性模型如下：M二B+PGDP+PPOP+PGen+pHost+pSys+（3）01S2S345其中GDP、POP、Gen均取奥运会举

18、办前一年的数据，奥运会东道主国家取1,非东道主国家取0；社会主义国家取1，非社会主义国家取0。5.2.2、模型三的求解：5.2.2.1、变量的统计描述和标准化：由于世界发展的不断变化，影响奥运会奖牌榜的各个因素在不同的时期差异颇大，为使各个因素在不同时期具有可比性，首先对影响奥运会各个参赛国奖牌榜排名的的因素进行标准化处理，计算了每个国家各个指标的份额，即X(j)=,X表示影响变量，工X(j)j表示每一项指标。5.2.2.1、模型三的求解：以美国为例求解模型三：表五：1988-2015年每届奥运会美国各项指标年份届数金牌占比奖牌占比GDP占比人口占比青年占比GDP增长率i-、东道主19882414.94%12.72%37.82%4.65%21.6