古塔的变形-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坐标iY<<古塔测量的数据y坐标X中心点的x坐标(1i108)iX<<古塔测量的数据z坐标Y中心点的y坐标(1i108)iZ中心点的z坐标N 古塔每层的测量点的个数α塔的倾斜角五、 模型的建立与求解5.1 问题1模型的建立与求解正八边形的重心等于中心，所以可以用均值法求每个面的中心点，公式如下：根据每个面内点的坐标(X i ,Y i ,Z i )，可求得平面的中心坐标：1Nii XX N ==∑ 1Nii YY N ==∑ 1Nii ZZ N ==∑ 由于每个面都有八个测量点，所以在这里N=8。

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

下面以求1996年古塔的中心点为例。

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

1352.83Z =Matlab curvefitting 软件对古塔X -Z 坐标系12层中心点和一个塔顶坐标进行拟合，结果如下图：Figure 1Linear 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 =同上用Matlab curvefitting软件对古塔Y-Z坐标系12层中心点和一个塔顶坐标进行拟合，结果如下图：Figure 2Linear model Poly1:f(x) = p1*x + p2Coefficients (with 95% confidence bounds):p1 = -114.5 (-133, -96.11)p2 = 5.987e+04 (5.025e+04, 6.95e+04)拟合公式：Z= -114.5*Y+59870将1352.83Z=代入拟合公式求得13522.2280Y=依据此方法，求出1986年的13层中心坐标。

5.2问题2模型的建立对于第二个问题，可以分成三个问题，倾角问题，弯曲问题，和扭曲问题。

5.2.1古塔的倾斜角图1图2图3图4图5图6图7图 8从上述八幅图中，不难看出，古塔已经向着X 轴的正方向，Y 轴的负方向发生了倾141cos Sα=进而求得倾斜角α的值：141arccos z zSα-=5.2.2古塔的弯曲Matlab curvefitting软件对古塔进行X-Z坐标系中心坐标曲线拟合，Y-Z坐标系中心坐标曲线拟合，结果如下图。

Figure 3General model Gauss1:f(x) = a1*exp(-((x-b1)/c1)^2)Coefficients (with 95% confidence bounds):a1 = 54.24 (50.25, 58.22)b1 = 567.3 (567.2, 567.3)c1 = 0.4135 (0.3546, 0.4723)拟合公式：2567.3()0.41554.24*XZ e--=然后代入曲率公式''3'22(1)ZKZ=+Figure 4General model Gauss1:f(x) = a1*exp(-((x-b1)/c1)^2) Coefficients (with 95% confidence bounds):a1 = 55.34 (53.18, 57.49) b1 = 522.3 (522.3, 522.3) c1 = 0.2665 (0.2421, 0.291)拟合公式：2522.3()0.266555.34*YZ e--=然后代入曲率公式''3'22(1)ZKZ=+5.2.3古塔的扭曲5.3分析古塔的变形趋势图9如图9所示，古塔的右下方的等高线比较密集，古塔左上方的等高线比较稀疏，说明古塔在X轴的正方向，Y轴的负方向已经有倾斜，在未来忽略不可抗力，古塔会一直沿着这个方向倾斜。

六、模型的推广与改进本模型简单易懂。

本模型解决了，古塔的各层中心点的确定，古塔倾斜角的求解。

改进建议，通过对倾斜角的取值，应该可以预测古塔的倒塌时间。

七、参考文献1、石宁刘竞刘青桂高等数学中国水利水电出版社2010年7月2、梁国业廖建平数学建模冶金工业出版社2004年9月八、附录1、求古塔倾斜角matlab程序function angle = angleacos(x,y,z)%UNTITLED5 Summary of this function goes here% Detailed explanation goes heres=sqrt((x(14)-x(1))^2+(y(14)-y(1))^2+(z(14)-z(1))^2);angle=acos((z(14)-z(1))/s);angle=180/pi*angle;end2、matlab画出古塔的3维图程序hold on;plot3(tower1986x, tower1986y, tower1986z);plot3(tower1986x, tower1986y, tower1986z,'g.');plot3(towercentre1986x, towercentre1986y, towercentre1986z,'r*');xlabel('X');ylabel('Y');zlabel('Z');grid on;hold off怎样写作数学建模竞赛论文一如何建立数学模型—建立数学模型的涉骤和方法建立数学模型没有固定的模式，通常它与实际问题的性质、建模的目的等有关。

当然，建模的过程也有共性，一般说来大致可以分以下几个步骤：1. 形成问题要建立现实问题的数学模型，首先要对所要解决的问题有一个十分明晰的提法。

只有明确问题的背景，尽量弄清对象的特征，掌握有关的数据，确切地了解建立数学模型要达到的目的，才能形成一个比较明晰的“问题”。

2. 假设和简化根据对象的特征和建模的目的，对问题进行必要的、合理的假设和简化。