2013数学建模-古塔变形

承诺书

我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。

我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）：

我们的参赛报名号为（如果赛区设置报名号的话）：

所属学校（请填写完整的全名）：

参赛队员(打印并签名) ：1.

2.

3.

指导教师或指导教师组负责人(打印并签名)：

（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。）

日期： 2013 年 09 月 13 日赛区评阅编号（由赛区组委会评阅前进行编号）：

编号专用页

赛区评阅编号（由赛区组委会评阅前进行编号）：

全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：

古塔的变形

摘要

中国古塔---是中国五千年文明史的载体之一。古塔为祖国城市山林增光添彩，被佛教界人士尊为佛塔。矗立在大江南北的古塔，被誉为中国古代杰出的高层建筑。但是，由于长时间承受自重、气温、风力等各种作用，偶然还要受到地震等自然灾害的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。

现在已知某测绘公司先后四次对某古塔进行测量的数据，要求通过这些数据来分析古塔变形的情况以及变形的趋势。

针对问题一，根据古塔每层已知的八个测量点，本文利用最小二乘法，通过Mathematica数学软件中的拟合曲线fit[]函数将这八个点的x、y坐标拟合成一个圆，则圆心就是各层的平面中心。最后对各测量点的z坐标求平均值，即可得到各层的空间中心。

针对问题二,倾斜程度可以通过古塔塔尖坐标相对于底层坐标在竖直方向偏移的角度来判断；弯曲度通过古塔各相邻两层之间的中心轴线的夹角大小反映；扭曲度通过古塔顶层各测量点相对底层相对应的测量点在水平方向偏移的角度大小反映。

针对问题三，通过问题二得出的数据分析古塔各方面变形程度的强弱。进而判断该古塔变形的趋势。最后制定出相应的保护措施，防止古塔继续形变、遭到破坏。

关键词: 最小二乘法 Mathematica Matlab 数据拟合

一、问题重述

因为长期承受自重、气温、风力等各种作用，偶然还要受到地震、飓风的影响。所以古塔会产生各种形变。出于对古塔保护的目的，文物部门委托测绘公司对某塔适时进行测量。然后通过对测量数据的分析，总结出古塔变形的情况以及变形的趋势。最后制定出相应的保护措施。防止古塔继续形变，遭到破坏。

二、问题分析

2.1 问题一分析

利用Matlab软件将附件1中的坐标画在坐标轴上，然后从俯视图中我们可以清楚得看出，该塔为八边形（如图1）。又因为在测量数据中，1986年和1996年该塔第13层的第5个数据空白，而2009年和2011年该数据有测量值。则分析可知，该塔在1986年到1996年之间，第13层的第5个角遭到破坏，1996年到2009年之间对该塔进行过维修，并将被破坏的角修好。这也很好的解释了附件1中1986年和1996年第13层第5个数据空白，而2009年和2011年该数据有测量值。同时可判断，对该塔是以每层的八个角为测量点进行测量。于是我们利用Mathmatic软件拟合曲线中的fit[]函数将每层的八个点拟合为圆。然后再联系圆的圆心和八个角的纵坐标的平均值，得出每层的中心坐标。

图1 1986古塔坐标俯视图

2.2 问题二分析

○1古塔长时间受到各种因素的影响，使其主轴线发生偏移。从上往下投影，塔顶与塔底的中心不在一个点上。于是可以利用坐标求出塔顶与塔底中心在水平面上的距离L，然后结合塔顶的高度H构建直角三角形。最后在三角形中解得古塔中轴线与垂线的夹角θ，即该塔的倾斜度。如图2所示：

图2 古塔中轴线倾斜度示意图

○2根据第一问得出的中心坐标，将相邻三个坐标相连得到一个三角形。如∆kgf，然后利用matlab求出三角形中间角度的补角，即θ。则可知道上一层塔相对于下一层塔偏移的角度。即这两层塔弯曲的程度。如图3所示：

图3 古塔各层弯曲示意图

○3利用matlab将古塔的坐标画在坐标轴上，从俯视图中可以看见，同一方位各层的坐标相连是一条曲线。说明古塔已有扭曲形变。那么将同一方位古塔最底层坐标、最高层坐标、底层中轴线坐标相连得到一个三角形，如∆cba。然后利用matlab求出角度θ，即可知道顶层相对于底层扭曲了多少。如图4所示：

图4 古塔扭曲示意图

三、模型假设

1.假设古塔每层以八个角为测量点；

2.假设每次测量选择的坐标原点为定点；

3.假设古塔不受到人为破坏。

四、符号说明

H：塔高；

D：塔顶投影到塔底距离塔底中心的距离；

：古塔倾斜角度；

x i：第i次测量古塔顶层x坐标；

y i：第i次测量古塔顶层y坐标；

z i：第i次测量古塔顶层z坐标；

x i，：第i次测量古塔底层x坐标；

y i，：第i次测量古塔底层y坐标；

z i，：第i次测量古塔底层z坐标。

五、模型的建立和求解

○1由于古塔是八边形，每层的八个坐标都已知。所以可以直接利用Mathematica的fit[]函数将古塔每层八个点的x、y坐标拟合成圆。则圆心就是各层水平面的中心位置。然后再用平均法求出各层竖直方向的中心位置。由此即可确定古塔各层的中心位置。以计算1986年第1层中心位置为例。程序如下:

t=List[{565.454,528.012},{562.058,525.544},{561.395,521.447},{56 3.782,518.108},{567.941,517.407},{571.255,519.857},{571.938,523.