地震预测问题数学建模论文解析.doc

2015年北京交通大学数学建模竞赛 承 诺 书 我们仔细阅读了数学建模竞赛的竞赛规则。 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。 我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 日期： 年 月 日 赛区评阅编号（由赛区组委会评阅前进行编号）： 目录 一、摘要 2 二、问题的重述 3 三、模型的假设与符号说明 6 （1）、模型的假设 6 （2）、符号说明 6 四、模型的建立 2 （1）. 由原始数据计算一次累加数据 2 （2）. 建立矩阵 2 （3）. 求u和a 2 （4）. 时间响应方程 2 五、模型的求解 3 （1）. 矩阵求解 3 （2）.求u和a 3 （3）. 时间响应方程 4 六、模型检验或误差分析 6 （1）、残差检验： 6 （2）、后验差检验： 6 （3）、结论 7 （4）、误差分析 7 七、模型评价 7 八、参考文献 7 地 震 预 报 问 题 一、摘要 本文通过建立模型，解决了地震预报问题，简述如何利用正确的建模，从几组数据中预测下一次地震的时间与地点。 20世纪以来发生在中国8级以上大地震共四次，数据较少，很难找到规律，用神经网络等建模形式来“训练数据”较为困难，且到目前为止，地震的发生尚无规律可言，即使是相邻的两次地震，彼此之间也很难说有什么影响。对于地震来说，筛选出的四个数据，它们的时间地点没有必然的联系，所以可以将时间和地点的预测分开进行建模。 在时间上，由于数据较少，内部数据规律不能很好地表示出来，可以选择灰色预测，它是一种对含有不确定因素的系统进行预测的方法，可以弥补地震事件中的不确定性。通过关联分析几次地震时间，并对原始数据进行生成处理来寻找系统变动的规律，生成有较强规律性的数据序列，然后建立相应的微分方程模型，从而预测事物未来发展趋势的状况，预测未来某一时刻的特征量，或达到某一特征量的时间。最终运用matlab进行编程得出结果。 在地点上，根据板块漂移假说和地震带理论，发生在地震带上的可能性较高。对于地点，主要通过经纬度来确定，假设经度和纬度之间没有关系，通过分别建模预测下次地震的经度和纬度。由于数据较少，可以运用曲线拟合法，通过图像的类比得到近似的函数。最小化问题的精度，依赖于所选择的函数模型，用q个相关变量或p个附加的相关变量去拟和。通常人们将一个可能的、对不相关变量的构成都无困难的函数类型充作函数模型(本题中运用正弦函数)。求出参数后，用参数使所选择的函数模型同观测值Y相匹配。最终得到可预测下次地震地点的函数。在求解方面，利用了mathematica软件求解。 经过建模求解，并进行误差分析之后，预测下次8级及以上的地震时间是2087年10月15日，地点是N36.55?,E106.76?（甘肃省内） 论文的末尾给出了模型优缺点的分析和评价，并提出了改进方向，如果数据再多一些，或者是要进行7级以上地震预测的话，会更准确一些。 关键词： 地震预报 时间大地震的破坏性是众所周知的，为了减少大地震带来的灾难，人们提出了各种预报地震的方法，以求减少大地震产生的破坏。本赛题请大家用数学建模的方式预报下一次大地震发生的时间和地点。为了减少问题的复杂性，请根据20世纪以来我国发生8级及以上大地震的时间和地点预报下一次我国发生8级及以上大地震的时间和地点。对于时间预测模型所涉及： X(0) ：实际时间 X(1) ：累加之后的时间 u、a : 时间响应方程中的常量 X’(1)（K+1）：时间预测值 E（k）：模型计算值与实际值之间的残差 E(k) ：残差与实际值之间的相对误差