洪峰流量预测数学建模.doc

2013年河南科技大学数学建模选拔赛 承 诺 书 我们仔细阅读了数学建模选拔赛的规则. 我们完全明白，在做题期间不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人研究、讨论与选拔题有关的问题。 我们知道，抄袭别人的成果是违反选拔规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺，严格遵守选拔规则，以保证选拔的公正、公平性。如有违反选拔规则的行为，我们将受到严肃处理。 我们选择的题号是（从A/B/C中选择一项填写）： 队员签名 ： 日期：2013年08月23日 2013年河南科技大学数学建模竞赛选拔 编 号 专 用 页 评阅编号（评阅前进行编号）： 评阅记录（评阅时使用）： 评 阅 人 评 分 备 注 洪峰流量预测 摘要 三峡大坝建成之时号称“固若金汤，可抵万年之洪”；五年之后，能防千年之大涝；时至今日，抗洪之任则不可博弈于三峡。就近年来三峡大坝的非全能抗洪能力的问题，本文从理性的角度，结合某水文站最大洪峰的观测值，运用相关数学知识，评估大坝的抗洪能力，探讨“百年一遇”洪峰流量值并对未来洪峰流量做出预测。 对于问题一：“百年一遇”并非一百年出现一次，它是一个关于频率的概念。水文站“百年一遇”的问题亦属频率发生事件，首先我们用matlab画出历年洪峰流量折线图，然后对洪峰流量进行升序排列后进行数据分析猜测其可能满足正态分布。之后用spss做出Q-Q图对洪峰流量进行正态分析性检验，判断其数据是否服从正态分布，然后画出洪峰流量频率散点图，“百年一遇”即为图中频率为1%时的相应数据值。 对于问题二：对于百年一遇的流量的计算，即为对问题一的定量分析，先用spss画出已给数据频率散点图，然后用matlab画出皮尔逊Ⅲ型曲线进行拟合从而确定频率P为1%时的最大洪峰流量值，在运算皮尔逊Ⅲ型曲线时可以通过配线法计算相关变量,“百年一遇”最大洪峰流量值即可迎刃。 对于问题三：对于预测未来三年最高洪峰流量，属于预测问题，首先建立了三层BP神经网络预测模型，对预测的原理进行了有力的阐述，然后对所给的数据重新编排、整理，调整输入、输出数据，根据所建立的BP神经网络预测模型，使用Matlab编写计算程序，将整理的数据代入程序，即可预测出未来三年的最大洪峰流量值。 关键词：PⅢ型曲线、配线法、水文频率计算、BP神经网络预测模型 问题重述 洪峰流量预测 外界对三峡工程“万年一遇”“千年一遇”“百年一遇”等防洪标准说法不一提出质疑，并将相关报道整理如下： 现假设附件中是某水文站的每年的最大洪峰流量观测值。 有人说，百年一遇就是一百年内只出现过一次，能否这样理解？给出你对“百年一遇”的定义。 能否计算出附件中水文站百年一遇的洪峰流量值；如能，给出具体的计算过程及结果。 预测该水文站将来3年内每一年的最大洪峰流量。 年份 最大洪流量 年份 最大洪流量 1985 638.1 1999 521.35 1986 658.99 2000 535.82 1987 577.11 2001 348.53 1988 502.42 2002 463.45 1989 384.55 2003 502.04 1990 455.04 2004 597.32 1991 405.57 2005 472.09 1992 555.06 2006 433.74 1993 510.22 2007 610.67 1994 578.33 2008 512.23 1995 470.9 2009 465.35 1996 834 2010 555.56 1997 703.92 2011 485.52 1998 597.79 2012 672.23 符号说明 :设计洪峰流量； : 测量年数； : 实测系列在n中从大到小排列序号； : 样本频率； : 模比系数； : 均方差； : 变差系数； : 偏态系数； 、：连接权系数； 、：阀值； ：隐含层输出； ：实际输出； ：校正误差； 模型假设 所有属性变量均为有序列属性变量 不考虑属性变量和数值变量在模型中的交叉作用 不考虑模型对变量的多次回归相关性 、数据经转换后保留了原始数据的某些统计特征，即问题的可预测性 、本文就以该水文站为单一研