原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
2014年数模实验报告《100米最好成绩的灰色预测》
实验五：100米最好成绩的灰色预测
2014-11-25
问题表述
利用互联网查询2014年以前10年的100米最好成绩（分男子和女子两种情况），然后利用灰色理论预测2015年和2016年的最好成绩，并作精度检验。
二．实验过程与结果（含程序代码）
本实验我们只讨论2014年以前男子100米的最好成绩。
第一步：经查询后2014年以前10年的男子100米最好成绩如下表所示：
年 份 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 男子最好成绩 9.78 9.77 9.76 9.74 9.69 9.58 9.72 9.79 9.63 9.77
原始序列为：
即 =（9.78,9.77,9.76,9.74,9.69,9.58,9.72,9.79,9.63,9.77）
第二步：用MATLAB作1-AGO生成序列
第三步：再作的一阶均值生成，得
即：其中： k= 2, 3…..，n 即构成了灰色模块，
第四步：计算矩阵B,Y 。
则可得出
第五步：由最小二乘法求得：，
利用MATLAB计算出，程序如下：
B=[-14.665,1;-24.43,1;-34.18,1;-43.895,1;-53.53,1;-63.18,1;-72.935,1;-82.645,1;-92.345,1]
Y=[9.77,9.76,9.74,9.69,9.58,9.72,9.79,9.63,9.77]
inv(B*B)*B*Y
ans =0.0004
9.7406
得出
第六步：即可建立灰色模型，GM（1，1）模型的一般式为：
解此微分方程得： ；
第七步：取 ，再利用MATLAB计算出，程序如下：
function y=eg1_1f(k)
for k=1:10
y(k)=(-24341.72)*exp(-0.0004*(k-1))+24351.5
End
eg1_1f
ans =
Columns 1 through 9
9.7800 19.5147 29.2456 38.9725 48.6956 58.4148 68.1301 77.8415 87.5490
Column 10
97.2526
当k=0时，=9.78 所以
第八步：用X(1)的模拟值累减还原求出X(0)的模拟值,取。
（k=2,...,10）得：
= 为其GM(1,1)时间响应式的累减还原值，则当时，称为模型模拟值；当tn时，称为模型预测值。用MATLAB计算预测值（）程序如下：
function y=eg1_1f(k)
for k=1:12
y(k)=(1-exp(0.0004))*(-24341.72)*exp(-0.0004*(k-1));
end
ans = 9.7386 9.7347 9.7308 9.7270 9.7231 9.7192 9.7153 9.7114 9.7075 9.7036 9.6998 9.6959
第九步：用残差检验。MATLAB计算程序如下：
x0=[9.78,9.77,9.76,9.74,9.69,9.58,9.72,9.79,9.63,9.77]
x0 = 9.7800 9.7700 9.7600 9.7400 9.6900 9.5800 9.7200 9.7900 9.6300 9.7700
x1=[9.78,9.7347,9.7309,9.7269,9.7231,9.7192,9.7153,9.7114,9.7075,9.7036]
x1 = 9.7800 9.7347 9.7309 9.7269 9.7231 9.7192 9.7153 9.7114 9.7075 9.7036
q0=x0-x1
q0 = 0 0.0353 0.0291 0.0131 -0.0331 -0.1392 0.0047 0.0786 -0.0775 0.0664
q1=sum(q0)/10
q1 = -0.0023
x2=sum(x0)/10
x2 = 9.7230
s1=sum((x0-x2).^2)/10
s1 = 0.0044
s2=sum((q0-q1).^2)/10
s2 = 0.0039
C=sqrt(s2/s1)
C = 0.2450
参照下列P,C表
好
文档评论（0）