数学建模大赛C题

1、2013高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛

2、规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： C 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）：石家庄职业技术学院参赛队员 (打印并签名) ：1.魏鹏飞 2.邢磊 3.刘力恒 指导教师或指导教师组负责人 (打印并签名)：陈佩宁 （论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。） 日

3、期：2013年9月16日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：C题：古塔的变形摘要古塔由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。对于第一个问题，求中心点坐标，采用的是均值法，由于前两次测量中第13层第5个

4、点没有数据，要是采用均值法求中心坐标，会产生较大的误差，所以在求第13层中心坐标，采用的是拟合法。对于第二个问题，分析古塔倾斜、弯曲、扭曲等变形情况。这个问题可以分三个小问题考虑。1、 分析古塔的倾斜情况，先用Matlab软件绘制出，古塔的俯视图，观察古塔的倾斜情况，大致的倾斜方向，再用三角函数求出古塔的倾斜角度，再把四次算的倾斜角，做一下比较，观察古塔的倾斜状况。2、 分析古塔的弯曲情况，首先观察X-Z坐标系中心点坐标，用Matlab软件把X-Z坐标系中的中心坐标拟合成一条曲线，求出这条曲线的曲率，然后按照上述方法求出Y-Z坐标系中心点坐标的曲线方程，求出这条线的曲率，分别观察古塔在X轴方向

5、的弯曲情况，和Y轴方向的弯曲情况。3、 分析古塔的扭曲情况，由于时间关系，没有分析古塔的扭曲。对于第三个问题，分析古塔的变形趋势，可以根据第二问中的倾斜角，弯曲情况，进行简单的分析。关键词：Matlab拟合，Matlab绘图，均值法，Matlab curve fitting软件，Matlab编程一、 问题重述古塔由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。某古塔已有上千年历史，是我国重点保护文物。管理部门委托测绘公司先后于1986年7月、1996

6、年8月、2009年3月和2011年3月对该塔进行了4次观测。请你们根据附件1提供的4次观测数据，讨论以下问题：1. 给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。2. 分析该塔倾斜、弯曲、扭曲等变形情况。3. 分析该塔的变形趋势。二、 问题分析本文我们是根据相关人员对古塔的观测的数据来分析该古塔的倾斜与弯曲的程度，并且分析出未来古塔的趋势走向。首先，我们利用均值法求出各层的中心点并且拟合出图形，然后我们对该塔构建了三角形，利用三角函数求出该塔的倾斜角度，并且利用曲率算出弯曲的程度。最后，我们利用所求出的数据以及图表进行分析得到该塔未来的发展趋势。三、 模型假设1、

7、假设古塔每层都是正八边形。2、 假设题目中提供的数据真实可靠。3、 假设地面平整。4、 假设每层的测量点在一个平面内。四、 符号说明S 塔身长度 古塔测量的数据x坐标 中心点的x坐标 古塔测量的数据y坐标 中心点的y坐标 古塔测量的数据z坐标 中心点的z坐标 N 古塔每层的测量点的个数 塔的倾斜角五、 模型的建立与求解5.1 问题1模型的建立与求解正八边形的重心等于中心，所以可以用均值法求每个面的中心点，公式如下：根据每个面内点的坐标()，可求得平面的中心坐标：由于每个面都有八个测量点，所以在这里N=8。在求第13层中心点时，由于缺失数据，用均值法得出的中心坐标有很大的偏差，所以在求13层中心

8、点改用拟合法。下面以求1996年古塔的中心点为例。因为古塔的每层测量点都在一个平面内，所以13层的Z轴坐标为7个测量点Z轴坐标的平均值。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)拟合的公式为：= 93.06* -52730将代入拟合公式求得同上用Matlab cu

9、rvefitting软件对古塔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)拟合公式：= -114.5* +59870将代入拟合公式求得依据此方法，求出1986年的13层中心坐标。古塔各层的中心坐标，如下表：塔的各层中心点坐标年份层次X坐标Y坐标Z坐标年份层次X坐标Y坐标Z坐标1986

10、1566.6648522.71051.78737519961566.665522.71021.7832566.7196522.66847.320252566.7205522.66747.3146253566.7735522.627312.755253566.7751522.625612.750754566.8161522.594417.078254566.8183522.592217.075135566.8621522.559121.72055566.8649522.556321.7166566.9084522.524426.235136566.9118522.52126.22957566.9

11、468522.508129.836887566.9506522.504229.832258566.9843522.492433.350888566.9884522.488133.345389567.0218522.476436.854889567.0265522.471436.8482510567.0569522.462440.1721310567.062522.457240.1676311567.1045522.42344.4408811567.1102522.417344.4353812567.1518522.383648.7118812567.1578522.377548.7073813

12、567.204522.28652.8342913567.1951522.22852.83塔尖567.2473522.243855.12325塔尖567.2544522.236755.11975表格 1年份层次X坐标Y坐标Z坐标年份层次X坐标Y坐标Z坐标20091566.7268522.70151.764520111566.727522.70141.763252566.764522.66937.3092566.7642522.6697.29053566.8001522.638412.732253566.8004522.638712.726884566.8293522.613217.0697545

13、66.8297522.612717.0525566.8604522.586621.709385566.861522.58621.703886566.9471522.534226.2116566.9478522.533526.20457566.9792522.512329.824637566.98522.511529.8178567.0305522.479733.339888567.0313522.478833.336639567.0816522.446636.843759567.0825522.445736.8222510567.137522.393740.1611310567.1381522

14、.392640.1441311567.1799522.354744.4326311567.181522.353544.4248812567.2225522.31648.6997512567.2238522.314748.6838813567.2712522.271552.8183813567.2725522.270152.81313塔尖567.336522.214855.091塔尖567.3375522.213555.087续表15.2问题2模型的建立对于第二个问题，可以分成三个问题，倾角问题，弯曲问题，和扭曲问题。古塔的倾斜角图 1图 2图 3图 4图 5图 6图 7图 8从上述八幅图中，不

15、难看出，古塔已经向着X轴的正方向，Y轴的负方向发生了倾斜，从专业角度来讲倾斜是指基础两端点倾斜方向的沉降差与其距离的比值，所以我们先取第一层的中心和塔尖进行倾斜角的测量。首先测出两点间的空间距离，两点间的空间距离就是古塔的塔身长度S，计算公式：然后，求塔身与垂直方向的夹角的余弦值，公式为：进而求得倾斜角的值：把四年的测量数据分别带入公式可以得到古塔四次测量时的倾斜角，见下表年份1986199620092011倾斜角0.80180.81210.83770.8397古塔的弯曲Matlab curvefitting软件对古塔进行X-Z坐标系中心坐标曲线拟合，Y-Z坐标系中心坐标曲线拟合，结果如下图。

16、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)拟合公式：然后代入曲率公式Figure 4General model Gauss1: f(x) = a1*exp(-(x-b1)/c1)2)Coefficients (with 95% confidence bounds): a1 = 55

17、.34 (53.18, 57.49) b1 = 522.3 (522.3, 522.3) c1 = 0.2665 (0.2421, 0.291)拟合公式：然后代入曲率公式古塔的扭曲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*