2013数学建模全真模拟2.docVIP

2013数学建模全真模拟2 承 诺 书 我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛 赛区评阅编号（由赛区组委会评阅前进行编号）： 2013数学建模全真模拟2 编 号 专 用 页 赛区评阅编号（由赛区组委会评阅前进行编号）： 赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注 全国统一编号（由赛区组委会送交全国前编号）： 全国评阅编号（由全国组委会评阅前进行编号）： 地应力拓展映射关系 摘要 我们通过应用统计分析的方法即最小二乘法, 在残差平方和最小的前提下, 求得回归方程中各自变量系数的最优解, 用回归系数乘以每一种待定因素下的应力值, 即得全场的地应力反演计算值。找出最优解，以达到最优拟合的目的，求得工程区域初始应力场。为了证明我们找出的解是最优解，我们用将实测数据与拟合数据的三个力运用模糊的数学方法赋予权重，计算出最优解在进行拟合。 关键词： 统计分析方法 方差 残差平方和 matlab 一 问题的重述 二 模型的建立 模型的分析： 1 地应力的分布规律受很多因素影响，在本次分析中我们只考虑这四个因素：自重（V方向）、水平（H向）单位位移和水平（h向）单位位移和剪切位移，其他因素不予考虑。 2 实测数据与四种工况的数据都是埋深在500米以上数据的，所以水平地应力多数大于垂直地应力（布朗和霍克的结论）。 3 地应力是存在于地层中的未受工程扰动的天然应力，也称岩体初始应力、绝对应力或原岩应力。它是引起开挖工程变形和破坏的根本作用力，是进行围岩稳定性分析，实现岩石工程开挖设计和决策科学化的必要前提条件。所以要通过有效的计算过分析方法，进行第应力场的反演分析，以获得更为准确的适应范围较大的地应力场。 4 为了合并数据，计算出方差值，如下表： 表1 表2 5. 建立回归模型，并对结果进行检验 模型为： 至于检验，R、R2和F值可以直接从回归结果中读出。实际上，，检验通过。有了R值，F值和t值均可计算出来。F值的计算公式和结果为： 自重（V方向）、水平（H向）单位位移和水平（h向） 图1 图2 依照上面的方法得到下面的两个图： 图3 图4 二次拟合数据表： 序号 埋深/m 工况值 实测值 1 83.91 2.09897 3.142 2 102.11 2.4156 3.084 3 152.9 3.33952 4.727 4 162.23 3.66366 4.904 5 164.13 3.76854 4.939 6 146.9 2.8177 5.238 7 159.7 3.9334 5.073 8 188.4 4.53185 5.478 9 219.71 6.1359 5.533 10 255.51 7.38073 6.688 11 265.71 7.57351 7.933 12 273.91 8.0032 7.338 13 283.72 8.50917 7.618 14 284.72 8.66075 7.788 15 237.53 4.0124 10.606 16 293.93 5.14051 12.451 17 334.66 5.55701 11.756 18 367.01 5.95395 14.201 19 376.16 6.26202 13.723 表3 根据上表做出来的图： 图5 图6 观察散点图，判断点列分布是否具有线性趋势。只有当数据具有线性分布特征时，才能采用线性回归分析方法。从图中可以看出，本例数据具有线性分布趋势，可以进行线性回归。 回归结果如下： 表4 图7 总结：根据二次拟合曲线可以看出我们找出的解是最优解，使实测值与拟合值达到最优拟合。 附表一 题目中给出的数据表： 表1 实测数据 钻孔号 序号 埋深/m 实测值 /MPa /MPa /MPa ZK-1 1 83.91 4.43 3.34 2.29 2 102.11 3.61 3.03 2.79 3 152.9 6.12 4.03 4.17 4 162.23 6.21 4.13 4.43 5 164.13 6.23 4.15 4.48 ZK-2 6 146.9 7.53 4.97 3.97 7 159.7 6.66 4.6 4.31 8 188.4 6.85 4.39 5.09 9 219.71 5.76 4.2 5.93 10 255.51 7.42 5.06 6.9 11 265.71 9.72 7.16 7.17 12 273.91 8.3 5.74 7.4 13 283.72