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

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

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

我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。

如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。

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

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： C我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名) ：1.2.3.指导教师或指导教师组负责人(打印并签名)：（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。

以上内容请仔细核对，提交后将不再允许做任何修改。

如填写错误，论文可能被取消评奖资格。

）日期： 2013 年 9 月 13 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：对古塔变形问题的数学建模摘要中国古语有云，“救人一命胜造七级浮屠”，所谓浮屠也就是大众口中的“塔”。

在中国辽阔的大地上，古塔的踪影随处可见。

它们造型精美、结构巧妙，成为可多得的独特景观。

早起的古塔，主要是阁楼式的建筑，从唐朝经过两宋至辽、金，是我国古塔发展的高峰时期，特别是唐和两宋，古塔的建造达到了空前繁荣度，总量较以前大增，材料也更为丰富，除了木材和砖、石以外，还使用了铜、铁、琉璃等、材料上有木塔为主转为以石塔为主，平面则由四方形逐渐演变为六角和八角形。

由于古塔建造年代久远，碍于当时的科技技术，或多或少总存在一些问题，最常见的是地基不均匀沉降，从而导致塔体倾斜。

由于古塔建造的年代不同，地质情况千差万别，建造后所经历的人为破坏和自然力破坏也不尽相同，因此对古塔倾斜后处理与加固技术多种多样，归纳起来不外乎两类，即复位纠偏和加补强。

根据管理部门调查表，分别4次观测一个上千年的古塔在长时间承受自重、气温、风力，地震、飓风的影响下，古塔会产生各种变形，如倾斜、弯曲、扭曲等，为保护古塔变形情况，文物部门需要适时对古塔进行观测，了解各种变形情况，以制定必要的保护措施。

我们根据4次测量的数据，建立了古塔倾斜变形的预测模型。

根据这个模型，我们可以比较准确的预测出每一年的变形情况，然后根据情况文物部对其制定合理的保护措施。

针对问题1，我们首先要根据这几年来四次古塔的数据变化情况，用建模软件MATLAB制作成模型图，用数学建模中拟合的方法来画出塔的基本形状，再确定古塔每层的中心点，建立中心点拟合线方程模型，观察是否有倾斜、扭曲、变形等情况。

针对问题2，古塔倾斜的原因主要与日光照射、地基活动有关。

首先朝向阳面的地基水分较少，阴面的地基水分较多，于是1万多吨的塔，开始向水分较多，地基松软的方向倾斜。

另外大量的地下水开采影响了地基的稳固，而古塔附近的铁路运输，也会造成震动。

而且，受到地基的不均匀沉降、地震、大风等影响，都会有可能倾斜等变形情况。

针对问题3，根据管理部门委托的测绘公司的数据表来看，古塔每年都以很小的角度在偏移，由于各种人为或者自然原因，使得古塔慢慢的倾斜为斜塔，斜塔并不一定都会倒塌，只要塔的重心线（通过重心点所引的垂直线）还在塔的底面积范围内，塔就是安全的。

因此纠偏要根据每座塔的具体情况而定。

且一般来说，有些塔在倾斜的过程中，原本松软的地基会被渐渐压实，然后与倾斜角度构成新的平衡，便就此稳定下来。

关键词：古塔变形变形趋势预测倾斜弯曲扭曲维护文物中心位置0 问题背景在中国辽阔美丽的大地上，随处都可以看到古塔的踪影。

这些千姿百态的古塔，其造型之美，结构之巧，雕刻、装饰之华丽，均堪与我国其它种类的古代建筑相比。

然而，在我国早期的古代建筑物中有楼有阁，有台有榭，有廊有庑，有民居有桥梁有陵墓，唯独没有塔。

原来塔这种建筑并不是我国的固有类型，而是外国的一种建筑。

在传入我国以后，塔又和我国原有的建筑形式相结合，形成了一种具有中国民族传统特色的新的建筑类型。

由于古塔距今历史久远，期间经过各种风吹雨打，自然以及人为损坏，因此需要相关管理部门对古塔做测量，并且针对不同的变形情况制定不同的维护方法，以保证古塔的保留，防止历史建筑物的坍塌导致文化遗产的流逝与人员伤害。

1 问题重述古塔，是中国五千年文明史的载体之一，古塔为祖国城市山林增光添彩，塔被佛教界人士尊为佛塔。

矗立在大江南北的古塔，被誉为中国古代杰出的高层建筑。

但由于存在时间久远，受各方面影响，古塔的可能会发生变形。

如果不能及时发生古塔的变形，去及时维护，古塔就有可能遭受无法挽回的损失。

古塔的变形，由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。

为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。

某古塔已有上千年历史，是我国重点保护文物。

管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。

请你们根据附件1提供的4次观测数据，讨论以下问题：问题1：给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。

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

问题3：分析该塔的变形趋势。

2 问题分析由给出了的一个古塔实例以及相应数据（附件1，该实例古塔的4次观测数据）的条件下，要求建立该古塔的各层中心位置的通用方法，且列表给出各次测量的古塔各层的中心位置。

并进一步分析该古塔倾斜、弯曲、扭曲等变形情况，最后分析该古塔的变形趋势。

问题1中要我们确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心位置。

首先，我们假设各测量点都是选取得科学合理的位置，都是围绕中心点的，并且同一层测量点大致在同一平面上，由已知数据（附件1）也可以看出它们是大致在同一平面上。

那么，我们由已知条件知道每层给出的各测量点的数据，我们通过画三维图形可以看出，那近似于一个的多边形，所以我们就可以把问题转变为求多边形的中心位置，然后记录各个中心位置坐标数据。

问题2中要我们分析该塔倾斜、弯曲、扭曲等变形情况。

我们可以根据问题1中求得的各层的中心位置坐标，然后将各层的中心位置坐标，接连起来，观察它。

理论上正常古塔的中心位置坐标连线应该大致是一条垂直于X 轴和Y 轴平行Z 轴的直线。

如果中心位置坐标连线，还是直线但不平行于Z 轴了，说明该古塔发生了倾斜，如果中心位置坐标连线不是一条直线，而是一条有一个弧度的曲线，那么说明该古塔发生了弯曲，如果中心位置坐标连线不是一条直线，而是一条有多个弧度的曲线，那么说明该古塔发生了扭曲。

问题3要我们分析该塔的变形趋势时，我们通过分析每一个中心点的变化趋势，来判断整个塔的变形趋势，与预测古塔以后可能会发生的变形情况。

3 模型假设假设1：各测量点都是选取的都是科学合理的位置。

假设2：每层的测量点都是围绕着这一层的中心点。

假设3：同一层测量点都大致在同一平面上。

假设4：各层测量点构成的几何图形的中心是与这一层的中心位置相重合的。

假设5：测量点的位置都是古塔上固定的位置。

4 符号说明j O -------------第j 层中心；j --------------------层数(123...13)j ∈、、；i ---------------------测量点标记(123...)i ∈、、； x -------------x 轴坐标；y -------------------y 轴坐标；z ------------------z 轴坐标；X ------------------拟合线x 轴坐标；Y ------------------拟合线y 轴坐标；Z -------------------拟合线z 轴坐标；0jx ---------------------第j 层中心的x 轴坐标； 0j y ----------------------第j 层中心的y 轴坐标； 0jz ----------------------第j 层中心的z 轴坐标； n ----------------------测量点数t ----------------------时间t ∆----------------------时间差x ∆--------------------x 坐标差y ∆--------------------y 坐标差z ∆--------------------z 坐标差xQ ---------------------x 真实值与计算的近似值之差 y Q ---------------------y 真实值与计算的近似值之差,,,,,,M N P a b c d --------------------方程中的参数,,,j j j ja b c d -------------------------第j 层的方程参数5 模型建立与求解5.1确定古塔各层中心位置的通用方法，并列表给出古塔各层中心坐标根据几何中心计算方法，我们可以得出计算各层中心位置的通用方法。

101000010(x ,y ,z )n i j i n i j j j j i j n i j i x x n y o y n z z n ===⎧⎪⎪=⎪⎪⎪⎪=⎨⎪⎪⎪⎪=⎪⎪⎩∑∑∑即：我们把每层的8个测量值，分别求x ，y ，z 的均值，得到的坐标就是每层的中心位置。

楼层 中心位置坐标 x/m y/m z/m1 566.6648 522.7105 1.78742 566.7196 522.6684 7.32023 566.7735 522.6273 12.75524 566.8161 522.5944 17.07835 566.8621 522.5591 21.72056 566.9084 522.5244 26.23517 566.9467 522.5081 29.83698 566.9843 522.4924 33.35099 567.0218 522.4764 36.854910 567.0569 522.4230 40.172111 567.1045 522.4230 44.440912 567.1518 522.3836 48.711913 567.0850 522.7403 52.8343塔顶 567.2473 522.2437 55.1232表1数据分析：我们可以从表中发现，中心位置的x 轴坐标，在逐渐增加，而y 轴坐标在逐渐减少。