古塔的变形-2013年数学建模大赛C题

承诺书

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

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

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

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

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

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

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

所属学校（请填写完整的全名）：石家庄职业技术学院

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

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

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

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

编号专用页

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

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

C题：古塔的变形

摘要

古塔由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。

对于第一个问题，求中心点坐标，采用的是均值法，由于前两次测量中第13层第5个点没有数据，要是采用均值法求中心坐标，会产生较大的误差，所以在求第13层中心坐标，采用的是拟合法。

对于第二个问题，分析古塔倾斜、弯曲、扭曲等变形情况。这个问题可以分三个小问题考虑。

1、分析古塔的倾斜情况，先用Matlab软件绘制出，古塔的俯视图，观察古塔

的倾斜情况，大致的倾斜方向，再用三角函数求出古塔的倾斜角度，再把

四次算的倾斜角，做一下比较，观察古塔的倾斜状况。

2、分析古塔的弯曲情况，首先观察X-Z坐标系中心点坐标，用Matlab软件把

X-Z坐标系中的中心坐标拟合成一条曲线，求出这条曲线的曲率，然后按照

上述方法求出Y-Z坐标系中心点坐标的曲线方程，求出这条线的曲率，分

别观察古塔在X轴方向的弯曲情况，和Y轴方向的弯曲情况。

3、分析古塔的扭曲情况，由于时间关系，没有分析古塔的扭曲。

对于第三个问题，分析古塔的变形趋势，可以根据第二问中的倾斜角，弯曲情况，进行简单的分析。

关键词：Matlab拟合，Matlab绘图，均值法，Matlab curve fitting软件，Matlab编程

一、问题重述

古塔由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。

某古塔已有上千年历史，是我国重点保护文物。管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。

请你们根据附件1提供的4次观测数据，讨论以下问题：

1. 给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。

2. 分析该塔倾斜、弯曲、扭曲等变形情况。

3. 分析该塔的变形趋势。

二、问题分析

本文我们是根据相关人员对古塔的观测的数据来分析该古塔的倾斜与弯曲的程度，并且分析出未来古塔的趋势走向。首先，我们利用均值法求出各层的中心点并且拟合出图形，然后我们对该塔构建了三角形，利用三角函数求出该塔的倾斜角度，并且利用曲率算出弯曲的程度。最后，我们利用所求出的数据以及图表进行分析得到该塔未来的发展趋势。

三、模型假设

1、假设古塔每层都是正八边形。

2、假设题目中提供的数据真实可靠。

3、假设地面平整。

4、假设每层的测量点在一个平面内。

四、符号说明

S 塔身长度(1i108)

X<<古塔测量的数据x坐标

i

Y<<古塔测量的数据y坐标X中心点的x坐标(1i108)

i

X<<古塔测量的数据z坐标Y中心点的y坐标(1i108)

i

Z中心点的z坐标N 古塔每层的测量点的个数

α塔的倾斜角

五、 模型的建立与求解

5.1 问题1模型的建立与求解

正八边形的重心等于中心，所以可以用均值法求每个面的中心点，公式如下：

根据每个面内点的坐标(X i ,Y i ,Z i )，可求得平面的中心坐标：

1

N

i

i X

X N ==

∑ 1

N

i

i Y

Y N ==

∑ 1

N

i

i Z

Z N ==∑ 由于每个面都有八个测量点，所以在这里N=8。

在求第13层中心点时，由于缺失数据，用均值法得出的中心坐标有很大的偏差，所以在求13层中心点改用拟合法。下面以求1996年古塔的中心点为例。

因为古塔的每层测量点都在一个平面内，所以13层的Z 轴坐标为7个测量点Z 轴坐标的平均值。

1352.83

Z =

Matlab curvefitting 软件对古塔X -Z 坐标系12层中心点和一个塔顶坐标进行拟合，结果如下图：

Figure 1

Linear model Poly1:

Coefficients (with 95% confidence bounds): p1 = 93.06 (90.11, 96.01)

p2 = -5.273e+04 (-5.44e+04, -5.106e+04) 拟合的公式为：Z = 93.06*X -52730

将1352.83Z =代入拟合公式求得13567.1951X =