2013高教社杯全国大学生数学建模竞赛c题论文-古塔变形

《2013高教社杯全国大学生数学建模竞赛c题论文-古塔变形》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、古塔变形古塔变形一、问题分析：通过对数据的初步观察，知道个问题其实就是一个通过古塔的观测数据找到古塔变形的规律进而推测其未来发展的问题。我认为题目的关键是对这些数据进行拟合，然后根据函数的变化规律来进行预测。先对古塔大致形态进行分析，我先用Matlab将原始数据用三维图的形式描绘出来，这样有助于我对问题进一步分析与求解：由于前两次数据有缺失，故选择2009年的数据：古塔侧视图：25古塔变形古塔俯视图：通过对这个大致图形的观察，我掌握了古塔的一些基本信息，这是一个十三层的八角古塔，每一层都是由正八边形（由于测量误差或变形可能不

2、那么正，但古人肯定是这么设计的），于是可以大胆假设古塔的各层都是正八边形。二、模型假设：1）假设古塔的各层质量分布均匀，即重心即为古塔几何中心2）假设古塔的各层都是正八边形，即便有些数据误差，不过可以忽略二、符号申明：25古塔变形三、模型求解：由于1986年和1996年的数据有缺失（如下图13层第5点）：这对于利用几何法求中心是有阻碍的，于是先将它们补齐。根据假设2：古塔各层都是正八边形，故这里用几何法进行推算:由于是正八边形，如图，E为所求的点：由几何关系应该有：用坐标表示出来：即为：25古塔变形纵坐标可直接利用其余各点平

3、均值求得。利用Matlab求解这个方程，分别将1986和1996的数据带入方程，解得缺失的中心点的坐标分别为：年份XYZ1986568.1644519.625352.834291996568.1704519.619152.83于是现在每一次数据都是齐全的了，下面来求各层的中心坐标点：首先给出各层中心点求法的通用方法：首先，先将各层都拟合到一个平面，由于竖直面上各层差异不大，于是采用对各层各点z坐标取均值的方法将古塔的每一层拟合到同一个平面上,然后再求中心坐标，根据假设1：塔的各层分布均匀，故几何中心即为重心，在网上找到一篇论

4、文就是研究关于求多边形几何中心的方法，这里不妨直接拿来用一下：25古塔变形在这里，只需将n设置为8，利用Matlab求解公式即可得到各层的中心坐标。得到的结果，如下：1986年：1996年：2009年：25古塔变形2011年：散点连线图：25古塔变形将四年的放在一起看：25古塔变形现在大致可以知道1996到2009年这个时间段古塔变形较大，下面进行进一步量化分析：（一）、倾斜程度：对中心点进行线性拟合，求出拟合直线与竖坐标轴的夹角，夹角的余弦值可作为倾斜程度的标识：25古塔变形即：对各个年份的中心点进行拟合，采用线性插值法得

5、：25古塔变形放大后：解得的数据为：年份余弦值角度值正切值余弦值19860.1.93.852230.19960.1.92.804840.20090.1.86.353880.20110.1.86.170690.有表格可以看出，塔身与地面的夹角从1.变化到1.，逐渐靠近地面，下面利用Matlab对其拟合分析：25古塔变形原始散点图：年份-角度用三次样条插值后：25古塔变形由模拟的函数图可知，古塔的倾斜度每年都有所增加，但是在2000年到2010年这个时间段内变化特别明显，这说明在这个时间段内自然环境可能有比较剧烈的变化，联想到2

6、008年我国的汶川地震，可以知道那段时间我国的地壳运动是比较活跃的，由此猜测也许与地壳活动剧烈有关，到了2010年以后，变化速度急剧变慢，回到了90年代的速度，由此我大胆的做出一个与本题无关的推测，趋势变缓意味着地壳运动变缓，也就是说在未来的一段日子里我国应该不会发生强烈的地震灾害，所以大家可以安稳的生活。同时，我对古塔未来一段时间的倾斜度变化做出一个预测：变化存在，但很缓慢，与地面的夹角逐渐在变小，由图形走向来看，可能趋向停止（如果没有不可抗力因素，如地震狂风的影响）。（二）、弯曲程度：由于弯曲程度是指中心轴线的弯曲度，将

7、各个中心点连接起来，依次求出相邻点的之间的夹角，然后取平均值，作为弯曲度的度量：用MatLab解得为：年份余弦角度（弧度）1986-9.2E-051.1996-0.000191.2009-0.000291.2011-0.00041.25古塔变形用三次样条拟合后得到的图像：弯曲程度在1986到2007年之前都基本保持线性，但是2008到2010突然变为指数级增长，估计为自然灾害影响。所以可以预测，2011年以后，若无自然灾害或其他不可抗力的因素的影响，弯曲程度将与1986年左右持平，每年大约弯曲,其实这个只是很小的，至少要等上

8、千年才能看到较明显的变化。分析这个问题我还采用了另外一种做法，就是将各层坐标投影到xoz平面上，然后通过观测曲线曲率的变化来分析弯曲程度的变化：具体做法是，分别对每年都做一次三次样条插值，然后在曲线上抽样一些点，取这些点的曲率的平均值作为该年度塔的弯曲度的特征值，然后通过对改特征值的分析来