2013高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、我网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： C我们的参赛报名号为（如果赛区设置报名号的话）：sm063所属学校（请填写完整的全名）：山东现代职业学院参赛队员(打印并签名) ：1. 马昱轩2. 李倩3. 梁帅健指导教师或指导教师组负责人(打印并签名)：宋祖芳日期：2013年 9月 16 日赛区评阅编号（由赛区组委会评阅前进行编号）：sm063032013高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：古塔的变形模型摘要某古塔是我国重点保护文物,已有上千年历史。
由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生诸如倾斜、弯曲、扭曲等各种变形。
为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。
管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。
基于附件1提供的4次观测数据：对于问题1，1986年和1996年的观测数据中都缺少13层一个点的数据（因MATLA B程序中用的是循环语句，所以计算时赋予0值），其它各层均给出8个点的观测值，为使所得数据更具真实性，确定古塔各层中心位置的方法更具适用性，本文将每层所给8点构成的图形看做不规则八边形，用中垂线求交点法求得古塔各层中心坐标。
对于问题2，结合问题1的分析，采用垂直投影法[1]求古塔的倾斜度，根据所得数据，分析古塔的倾斜程度（因1986年和1996年13层赋予值后所得数据偏差较大，为使所得数据更具真实性，所以本问题起1986年和1996年13层数据予以舍弃）；弯曲是建立在二维平面上的一条曲线，通过截取古塔过x轴、z轴的界面，求出每层古塔的倾斜度，从而分析得出古塔塔身在x轴、z轴的界面的弯曲程度。
同理，也可分析y轴、z轴的界面的弯曲程度；扭曲同样是采用垂直投影法[1]求古塔每层的倾斜度，建立三维立体空间，根据所得数据，分析古塔塔身的扭曲程度。
对于问题3，利用问题1、2所得数据，进行合理的分析与猜想，进而分析出古塔塔身的倾斜、弯曲、扭曲等变化趋势。
本文的模型解决了题目给出的问题，计算过程中充分尊重观测数据，给出更符合实际的结果。
本文所得结果大部分由图表给出，结合图像，较为直观地表现出古塔变形情况。
结果表明，采用MATLAB数学软件可得出可靠结论。
关键词：古塔变形监测MATLAB垂直投影法倾斜度一、问题的重述由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。
为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。
某古塔已有上千年历史，是我国重点保护文物。
管理部门委托测绘公司先后于1 986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。
请你们根据附件1提供的4次观测数据，讨论以下问题：1. 给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。
2. 分析该塔倾斜、弯曲、扭曲等变形情况。
3. 分析该塔的变形趋势。
二、问题的分析2.1问题的背景某古塔是我国重点保护文物。
已有上千年历史，由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生诸如倾斜、弯曲、扭曲等各种变形。
对现存历史文物古塔的保护，掌握古塔的形变，显得尤为重要，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。
管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。
2.2问题的分析基于附件1提供的4次观测数据：2.2.1问题一的分析本题建立模型是利用中垂线（即中垂线相交于一点为中点)法，求解古塔各层中心的坐标，以确保古塔中心坐标的准确性。
2.2.2问题二的分析根据层与层之间的中心点，塔体倾斜变化采用垂直投影法[1](即垂线法，利用相对位移量进行监测)；弯曲是建立在二维平面上的一条曲线，通过截取古塔过x轴、z 轴的界面，求出每层古塔的倾斜度，从而分析得出古塔塔身在x 轴、z 轴的界面的弯曲程度。
同理，也可分析y 轴、z 轴的界面的弯曲程度；扭曲同样是采用垂直投影法求古塔每层的倾斜度，建立三维立体空间，根据所得数据，分析古塔塔身的扭曲程度。
2.2.3问题三的分析我们在对古塔进行安全性监测方面要满足两个方面：①塔体的倾斜变化；②塔身各段的变化[1]。
针对上述几个方面就采取对应的监测方法。
三、问题的假设1.假设求解古塔各层中心坐标时纵坐标先不予以考虑；2.假设古塔的地基土抗压能力强、韧性大、不易收缩，对古塔的变形不产生任何影响；3.假设对古塔观测数据的缺失可忽略不计；四、符号及文字说明x 表示建立古塔坐标轴的横坐标i x 表示古塔各层的第i (i =1,2,...,8）个点的横坐标y 表示建立古塔坐标轴的纵坐标i y 表示古塔各层的第i （i =1,2，...,8)个点的纵坐标z 表示建立古塔坐标轴竖坐标i z 表示古塔各层的第i (i =1,2,...,8)个点的竖坐标 i k 表示古塔的第i l 条直线的斜率'x 表示两个x 点之间的中点'y 表示两个y 点之间的中点H 表示古塔中心点到水平面的距离 D 表示古塔中心点的位移量Q表示古塔的倾斜度j表示古塔的第j(j=1,2,...,13)层数（MATLAB程序中）i表示古塔各层的第i(i=1,2,...,8）个点（MATLAB程序中）13⨯矩阵（MATLAB程序中）X表示关于横坐标的813⨯矩阵（MATLAB程序中）Y表示关于纵坐标的8K表示各层每两点之间中垂线的斜率（MATLAB程序中）xx表示每层8条中垂线两两相交所得交点的关于横坐标的813⨯矩阵（MATLAB程序中）13⨯矩阵（MATLAB程序yy表示每层8条中垂线两两相交所得交点的关于纵坐标的8中）五、模型的建立与求解（一）问题一的模型建立与求解1.通过中垂线的方法对古塔各层中心坐标进行求解，根据附件1表格数据，1986年和1996年的观测数据中都缺少13层一个点的数据计算时赋予0值。
然后观察数据，古塔的每i层八个点的高度变化不大，假设把古塔的每一层的点看作是在一个平面上，求解古塔各层中心坐标时纵坐标先不予以考虑，之后用平均值的方法计算出。
分析过程如下图1-1：（用MATLAB软件作图（附件1，程序1））1-1 根据上图，做出任意两边的中垂线，交与一点O，如下图所示1-2（在上图上用画图工具作图）：1-2由图可知，'A 点坐标为()','y x ，'B 点坐标)","(y x 。
根据公式2'21x x x +=，得:2'21x x x +=，2'21y y y +=；2"32x x x +=，2"32y y y += ① 根据直线的斜率公式，求出直线AB l 的斜率A k 。
同理，求出直线BC l 的斜率B k 。
1212x x y y k A --=2323x x y y k B --=再根据中垂线定理1'-=kk ，求出中垂线1l ，2l 的斜率1k ，2k ； 122111y y x x k k A --=-=233221y y x x k k B --=-= ② 设中垂线的直线方程1l 为：)'('1x x k y y -=- ③ 设中垂线的直线方程2l 为：)"("2x x k y y -=- ④ 通过方程③、④建立方程组:⎩⎨⎧-=--=-)"(")'('21x x k y y x x k y y解得： 2121"''"k k x k x k y y x -+--=； ''21y x k x k y +-= ⑤再把①、②式中所求的数据代入⑤中，即可解出两中垂线交点O 坐标。
同理，通过数学软件MATLAB 编程求解（附录1，程序2）可得出每层八条中垂线两两相交的八个交点坐标，得出数据并对这八个中垂线的交点求平均值，因1986年和1996年的观测数据中都缺少13层一个点的数据计算时赋予0值得出的数据偏差太大，此时予以舍弃。
塔顶数据是塔顶的四个点求平均值而得。
各层中心坐标如下表：而每层古塔都有一定的高度差，因此，假设每层古塔的中心高度为每层八点的平均高度，即：),,(i i i i z y x O ，因此，古塔每层的中心坐标是：（二）问题二的模型建立与求解1，根据问题一求出的每层的中心坐标数据，通过分析古塔的倾斜、弯曲、扭曲及层与层之间的中心点，确定了塔体倾斜变化采用垂直投影法[1](即垂线法，利用相对位移量进行监测);弯曲、扭曲即为塔身各段的变化采用求倾斜度的方法。
用MATLAB数学软件编程（附件1，程序3）和画图工具画图，如下图2-1：2-1 通过对上图的观察与分析，每层的中心大致在一条倾斜的直线上，采用垂直投影的方法，对古塔的倾斜角进行求解，如下图通过MATLAB数学软件编程和画图工具并用作图（附件1，程序4），建立的数学模型，如下图2-2；2-2通过对图形的观察，()o o o z y x A ,,点为古塔塔顶的中心坐标，()111,,z y x O 点为古塔第一层的中心坐标，()1,,z y x B o o 点是A 点在O 点所在水平面上的投影， 根据倾斜度公式DH Q =,[3]计算出Q 的值：2112121)()(z z y y x x D o o -+-+-）（1z z H o - 21121211)()()(z z y y x x z z Q o o o -+-+--=比较数据：2011200919961986Q Q Q Q >>>分析数据，Q 的值越小，倾斜程度越严重，古塔随着年代的增加而越来越倾斜。