数学建模 古塔的变形

2013高教社杯全国大学生数学建模竞赛承诺我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：C 我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名)：1.2.3.指导教师或指导教师组负责人(打印并签名)： （论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。）日期：2013年9月13日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）： 评阅人 备注 全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：对古塔变形问题的数学建模摘要 中国古语有云，“救人一命胜造七级浮屠”，所谓浮屠也就是大众口中的“塔”。在中国辽阔的大地上，古塔的踪影随处可见。它们造型精美、结构巧妙，成为可多得的独特景观。早起的古塔，主要是阁楼式的建筑，从唐朝经过两宋至辽、金，是我国古塔发展的高峰时期，特别是唐和两宋，古塔的建造达到了空前繁荣度，总量较以前大增，材料也更为丰富，除了木材和砖、石以外，还使用了铜、铁、琉璃等、材料上有木塔为主转为以石塔为主，平面则由四方形逐渐演变为六角和八角形。 由于古塔建造年代久远，碍于当时的科技技术，或多或少总存在一些问题，最常见的是地基不均匀沉降，从而导致塔体倾斜。由于古塔建造的年代不同，地质情况千差万别，建造后所经历的人为破坏和自然力破坏也不尽相同，因此对古塔倾斜后处理与加固技术多种多样，归纳起来不外乎两类，即复位纠偏和加补强。根据管理部门调查表，分别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符号说明-------------第层中心；--------------------层数；---------------------测量点标记；-------------轴坐标；-------------------轴坐标；------------------轴坐标；------------------拟合线轴坐标；------------------拟合线轴坐标；-------------------拟合线轴坐标；---------------------第层中心的轴坐标；----------------------第层中心的轴坐标；----------------------第层中心的轴坐标；----------------------测量点数----------------------时间----------------------时间差--------------------坐标差--------------------坐标差--------------------坐标差---------------------真实值与计算的近似值之差---------------------真实值与计算的近似值之差--------------------方程中的参数-------------------------第层的方程参数5模型建立与求解确定古塔各层中心位置的通用方法，并列表给出古塔各层中心坐标 根据几何中心计算方法，我们可以得出计算各层中心位置的通用方法。即：我们把每层的8个测量值，分别求x，y，z的均值，得到的坐标就是每层的中心位置。 表1：1986年各层中心位置 楼层 中心位置坐标 x/m y/m z/m 1 2 3 4 5 6 7 8 9 10 11 12 13 塔顶 表1数据分析： 我们可以从表中发现，中心位置的x轴坐标，在逐渐增加，而y轴坐标在逐渐减少。而正常的塔，理论而言，x轴和y轴坐标因该是不变的。 注：表1数据的计算代码见附件5。 结论：说明该古塔发生了倾斜变形。 表2：1996年各层中心位置 层 中心位置坐标 x/m y/m z/m 1 2 3 4 5 6 7 8 9 10 11 12 13 塔顶 表2数据分析： 我们可以从表中发现，中心位置的x轴坐标，在逐渐增加，而y轴坐标在逐渐减少。而正常的塔，理论而言，x轴和y轴坐标因该是不变的。 注：表2数据的计算代码见附件5。 结论：说明该古塔发生了倾斜变形。 表3：2009年各层中心位置 层 中心位置坐标 x/m y/m z/m 1 2 3 4 5 6 7 8 9 10 11 12 13 塔顶 表3数据分析： 我们可以从表中发现，中心位置的x轴坐标，在逐渐增加，而y轴坐标在逐渐减少。而正常的塔，理论而言，x轴和y轴坐标因该是不变的。 注：表3数据的计算代码见附件5。 结论：说明该古塔发生了倾斜变形。 表4：2011年古塔各层中心位置坐标 层 中心位置坐标 x/m y/m z/m 1 2 3 4 5 6 7 8 9 10 11 12 13 塔顶 表4数据分析： 我们可以从表中发现，中心位置的x轴坐标，在逐渐增加，而y轴坐标在逐渐减少。而正常的塔，理论而言，x轴和y轴坐标因该是不变的。 注：表4数据的计算代码见附件5。 结论：说明该古塔发生了倾斜变形。分析该塔倾斜、弯曲、扭曲等变形情况 我们由已知的空间直线的点向式方程 整理得直线射影式方程 其中；；；。这样直线可以看作是用这2个方程表示的平面相交的直线，所以可以分别对2个方程进行数据拟合。设表示按拟合方程求得的近似值。一般地，它不同于实测值，两者之差同理可得当Q取最小值时a，b，c，d的值即为方程的系数，即满足下列方程时Q值最小。有令方程组（6）可写成其中： 根据n组数据点解方程组就可以得到a，b，c，d的值，也就可以得到拟合线方程。然后我们就可以通过中心位置点拟合线和拟合线方程来判断，古塔是否已经发生倾斜，弯曲或扭曲。（1）通过前的算法模型，我们建立出1986年中心位置拟合线方程为： 注：a，b，c，d参数的计算代码见附件6。图11986年古塔数据综合分析 分析1986年古塔变形情况： 通过观察图1的第4个子图，我们可以看出中心位置点基本上都是在拟合线上或附近的，说明我们的建立的模型正确的，是可靠性较高的模型。 我们通过分析表1中心点的位置坐标，中心点x坐标有逐层增加的趋势，y坐标有减小的趋势，然后结合拟合线方程和图4中的第1个子图，可知中心拟合线不与z轴平行，有向x轴数值增加方向倾斜。 我们再分析图1中的第3个子图，这是一个府视图，我们可以观察到古塔的测量点曲线图的曲线密集程度分布是不均匀的，也说明古塔发生变形。 注：图1生成的代码见附件7。 结论：综上分析我们可以得出古塔发生了倾斜变形，往x轴数值增加方向与y轴数值减少方向倾斜。（2）通过前的算法模型，我们可以建立起1996年中心位置拟合线方程为： 注：a，b，c，d参数的计算代码见附件6。图21996年古塔数据综合分析 分析1996年古塔变形情况： 通过观察图2的第4个子图，我们可以看出中心位置点基本上都是在拟合线上或附近的，说明我们的建立的模型正确的，是可靠性较高的模型。 我们通过分析表2中心点的位置坐标，中心点x坐标有逐层增加的趋势，y坐标有减小的趋势，然后结合拟合线方程和图4中的第1个子图，可知中心拟合线不与z轴平行，有向x轴数值增加方向倾斜。 我们再分析图2中的第3个子图，这是一个府视图，我们可以观察到古塔的测量点曲线图的曲线密集程度分布是不均匀的，也说明古塔发生变形。 注：图2生成的代码见附件7。 结论：综上分析我们可以得出古塔发生了倾斜变形，往x轴数值增加方向与y轴数值减少方向倾斜。（3）通过前的算法模型，我们可以建立起2009年中心位置拟合线方程为： 注：a，b，c，d参数的计算代码见附件6。图32009年古塔数据综合分析 分析2009年古塔变形情况： 通过观察图3的第4个子图，我们可以看出中心位置点基本上都是在拟合线上或附近的，说明我们的建立的模型正确的，是可靠性较高的模型。 我们通过分析表3中心点的位置坐标，中心点x坐标有逐层增加的趋势，y坐标有减小的趋势，然后结合拟合线方程和图4中的第1个子图，可知中心拟合线不与z轴平行，有向x轴数值增加方向倾斜。 我们再分析图3中的第3个子图，这是一个府视图，我们可以观察到古塔的测量点曲线图的曲线密集程度分布是不均匀的，也说明古塔发生变形。 注：图3生成的代码见附件7。 结论：综上分析我们可以得出古塔发生了倾斜变形，往x轴数值增加方向与y轴数值减少方向倾斜。（4）通过前的算法模型，我们可以建立2011年中心位置拟合线方程为： 注：a，b，c，d参数的计算代码见附件6。图42011年古塔数据综合分析 分析2011年古塔变形情况： 通过观察图4的第4个子图，我们可以看出中心位置点基本上都是在拟合线上或附近的，说明我们的建立的模型正确的，是可靠性较高的模型。 我们通过分析表4中心点的位置坐标，中心点x坐标有逐层增加的趋势，y坐标有减小的趋势，然后结合拟合线方程和图4中的第1个子图，可知中心拟合线不与z轴平行，有向x轴数值增加方向倾斜。 我们再分析图4中的第3个子图，这是一个府视图，我们可以观察到古塔的测量点曲线图的曲线密集程度分布是不均匀的，也说明古塔发生变形。 注：图4生成的代码见附件7。 结论：综上分析我们可以得出古塔发生了倾斜变形，往x轴数值增加方向与y轴数值减少方向倾斜。 表5：各层两次测量期间中心位置变化情况 1996-1986 2009-1996 2011-2009 变化量 1 2 3 4 5 6 7 8 9 10 11 12 13 塔顶 表5数据分析： 首先我们分析z轴上的数据变化，我们从z轴数值变化可知，z轴数值一直都是在减少，说明该古塔的各层中心高度一直在降低，但各层间的减少量不一样，说明不是地基下沉引起的变化，说明该古塔一定有倾斜变形。 我们再分析下x轴上的数据变化，我们从x轴数值变化可知，该古塔的中心位置坐标的x轴坐标几乎都在增加，有明显的趋势，说明古塔有往x轴正方向倾斜。 我们最后在分析下y轴上的数据变化，y轴坐标在2009至1996年间1到7层中心位置坐标有增加而8到塔顶有减少，说明该时期该古塔在y轴方向上可能有弯曲，但看表5中其它两个期间的数据y轴数值都是在减少说明这两个时间，该古塔又有向y轴数值减少方向倾斜。 从表5中我们看不到古塔有在某个方向上有，坐标数值增加了，再减少，再增加，或减少了，再增加，再减少，所以我们可以初步认定该古塔没有发生扭曲。 注：表5数据的计算代码见附件7。 结论：综上分析我们可以大致知道此古塔的倾斜趋势是往x轴正方向即x轴坐标增加方向倾斜，而在y轴方向上有弯曲的迹象，并且有向y轴坐标减少的方向倾斜，没有发生扭曲。表6：中心位置的变化情况 1996-1986 2009-1986 2011-1986 变化量 1 2 3 4 5 6 7 8 9 10 11 12 13 塔顶 表6数据分析： 首先我们分析z轴上的数据变化，我们从z轴数值变化可知，z轴数值一直都是在减少，说明该古塔的各层中心高度一直在降低，但各层间的减少量不一样，说明不是地基下沉引起的变化，说明该古塔一定有倾斜变形。 我们再分析下x轴上的数据变化，我们从x轴数值变化可知，该古塔的中心位置坐标的x轴坐标几乎都在增加，有明显的趋势，说明古塔有往x轴正方向倾斜。 我们最后在分析下y轴上的数据变化，y轴坐标在这4次测量间1到7层中心位置坐标有增加而8到塔顶有减少，说明该时期该古塔在y轴方向上可能有弯曲。 从表6中我们没有发现古塔有在某个方向上有，坐标数值有增加趋势了，再减少，再增加，或有减少趋势了，再增加，再减少，所以我们可以初步认定该古塔没有发生扭曲。 注：表6数据的计算代码见附件7。 结论：综上分析我们可以大致知道此古塔的倾斜趋势是往x轴正方向即x轴坐标增加方向倾斜，而在y轴方向上有弯曲的迹象，并且有向y轴坐标减少的方向倾斜，没有发生扭曲。 综上6次数据分析可得结论：该古塔发生了倾斜和弯曲，古塔的倾斜趋势是往x轴正方向即x轴坐标增加方向和y轴坐标的负方向即坐标减少方向倾斜，而在y轴方向上有弯曲，目前该古塔还没有发生扭曲变形。分析该塔的变形趋势 通过4次测量点的数据，我们建立了一个拟合线趋势分析模型，来预测该古塔的变形趋势，因为从上面的分析我们可知，该古塔没有地基下沉现象，所以我们可知道塔的高度坐标是一个因变量，所以不在考虑：模型正确定性与可用性检验：图5第1层到第6层的中心位置移动趋势拟合图图6第7层到塔顶的中心位置移动趋势拟合图 分析图5和图6： 第1、2、3、6、8、9、10、11、12、13和塔顶的拟合线模型都是与中心位置点非常吻合，说明这些层的模型是基本符合要求的，具有可靠性的，说明此模型是可用的。 注：图5和图6的生成代码见附件8。根据此模型（模型方程11），我们通过分析计算得到表7中每一层的参数数据：表7：各层中心位置变形趋势模型参数 参数层 a b c d 1 2 3 4 5 6 7 8 9 10 11 12 13 塔顶 表7数据分析：我们把表7中的数据带入对应的参数，就可以得到每一层的中心位置的x和y的坐标随时间的变化情况的模型。 分析参数a，a几乎都是正数，所以随着时间的增加，中心位置的x坐标将越来越大，这符合了我们问题2中分析的古塔变形情况，古塔会向x轴数值增大的方向倾斜。 分析参数c，c在1到7层为正数，而在8到塔顶，为负数，也验证问题2中分析的古塔变形情况，古塔会向y轴数值减少方向倾斜，和古塔会发生弯曲，弯曲点在7和8层之间。 我们也可以根据此模型方程结合模型方程（2）可以计算出z轴的坐标值。 注：表7数据的计算代码见附件8。 结论：该模型是可用的，具有可靠性，和扩展性。通过该模型（模型方程11）我们可以预测以后的时间，该古塔的中心位置坐标是怎么变化的，会变化多少，我们都可以通过此模型计算得到。 通过此模型分析该古塔我们得到该古塔会向x轴数值增大的方向逐渐倾斜，同时会向y轴数值减少方向逐渐倾斜，和古塔会逐渐发生弯曲，弯曲点在7和8层之间。6模型分析与评价从模型的假设来分析 在模型的假设中，我们所作的假设都非常符合于现实情况和实事的需求，如每层的测量点都是围绕在这层的中心点，如果不围绕，那说明测量点只是取的某个局部区域的测量点，不具有代表性。从模型的真实性来分析模型是在现有的数据条件下，总结出的合理的规律建立起来的模型，而根据模型求解出来的又如实的反映现实中的情况。如在运用空间直线拟合分析中心位置模型的可靠性与古塔倾斜情况。从模型的求解方法来分析 根据实际问题，分别对问题进行分类取量分析、数据拟合分析、线性规划分析等多种方法来进行建模求解，使模型更接近于真实性，准确性，和可靠性。 综上分析：该模型的具有合理性，真实性等优点，但它只是针对于该古塔的具体实例，缺乏模型的扩展性，对其中一些细节问题考虑得过于粗糙，也有其不足之处。如地域性，有的地方多大风，有的地方多地震。7模型推广 问题1中所建立的模型可以用来求古塔的中心位置，只要有合理的，符合我们假设的测量点，我们就可以计算出古塔的中心位置。 问题2中有用到中心位置点拟合的模型，可以用来分析古塔的变形情况，在得到中心位置数据后，我们就可以通过分析中心位置点的拟合线情况判断古塔发生什么变形。 问题3我们可以古塔的预测模型可以预测将来的某年的各层中心位置变化情况。近而我们可以在根据问题2中的分析模型得知，古塔会发生什么样的变形，是倾斜，还是扭曲还是弯曲，如果有弯曲，或扭曲，弯曲点在哪，扭曲点又在哪。8对古塔的维护建议 由于古塔建造的年代不同，地质情况千差万别，建造后所经历的人为破坏和自然力破坏也不尽相同，因此，对古塔倾斜后处理与加固技术多种多样，归纳起来不外乎两类，即复位纠偏和加固补强。 另外应该根据当地的地理环境，气候环境制定合理的保护与维护措施，如古塔处是什么地质，常年刮什么风。对以经发生了变形的古塔，一定要及时修护，同时分析变形的原因，针对变形原因，制定预防。9参考文献[1]梁海奎，古塔变形测量方法探讨，城市勘测，2011年03期：113-114页，2011。[2]陈淑贞，空间直线方程的解题探讨，海南师范大学学报(自然科学版)，2011年03期：348-351页，2011。[3]陈兆斗高瑞，高等数学（工本），北京海淀区成府路205号：北京大学出版社，2006。[4]尚继红程晓伟，古塔加固工程实例，甘肃科技，2008年第7期：133-135页，2008。[5]胡志晓，古塔倾斜观测和数据分析，江苏建筑，2011年第6期：34-35页，2011。[6]陈东佐康玉庆，中国古塔的维修与保护，太原大学学报，2006年第4期：75-78页，2006。[7]桂占吉陈修焕杨亚辉，基于MATLAB高等数学实验，武汉：华中科技大学出版社，2010年。10附录附件1 1986年观测数据 层 点 坐标 x/m y/m z/m 1 1 2 3 4 5 6 7 8 2 1 2 3 4 5 6 7 8 3 1 2 3 4 5 6 7 8 4 1 2 3 4 5 6 7 8 5 1 2 3 4 5 6 7 8 6 1 2 3 4 5 6 7 8 7 1 2 3 4 5 6 7 8 8 1 2 3 4 5 6 7 8 9 1 2 3 4 5 6 7 8 10 1 2 3 4 5 6 7 8 11 1 2 3 4 5 6 7 8 12 1 2 3 4 5 6 7 8 13 1 2 3 4 5 　 　 　 6 7 8 塔尖 1 2 3 4 1996年观测数据 层 点 坐标 x/m y/m z/m 1 1 2 3 4 5 6 7 8 2 1 2 3 4 5 6 7 8 3 1 2 3 4 5 6 7 8 4 1 2 3 4 5 6 7 8 5 1 2 3 4 5 6 7 8 6 1 2 3 4 5 6 7 8 7 1 2 3 4 5 6 7 8 8 1 2 3 4 5 6 7 8 9 1 2 3 4 5 6 7 8 10 1 2 3 4 5 6 7 8 11 1 2 3 4 5 6 7 8 12 1 2 3 4 5 6 7 8 13 1 2 3 4 5 　 　 　 6 7 8 塔尖 1 2 3 4 2009年观测数据 层 点 坐标 x/m y/m z/m 1 1 2 3 4 5 6 7 8 2 1 2 3 4 5 6 7 8 3 1 2 3 4 5 6 7 8 4 1 2 3 4 5 6 7 8 5 1 2 3 4 5 6 7 8 6 1 2 3 4 5 6 7 8 7 1 2 3 4 5 6 7 8 8 1 2 3 4 5 6 7 8 9 1 2 3 4 5 6 7 8 10 1 2 3 4 5 6 7 8 11 1 2 3 4 5 6 7 8 12 1 2 3 4 5 6 7 8 13 1 2 3 4 5 6 7 8 塔尖 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 2011年观测数据 层 点 坐标 x/m y/m z/m 1 1 2 3 4 5 6 7 8 2 1 2 3 4 5 6 7 8 3 1 2 3 4 5 6 7 8 4 1 2 3 4 5 6 7 8 5 1 2 3 4 5 6 7 8 6 1 2 3 4 5 6 7 8 7 1 2 3 4 5 6 7 8 8 1 2 3 4 5 6 7 8 9 1 2 3 4 5 6 7 8 10 1 2 3 4 5 6 7 8 11 1 2 3 4 5 6 7 8 12 1 2 3 4 5 6 7 8 13 1 2 3 4 5 6 7 8 塔尖 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　附件2%1986年古塔第一层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第二层，中心位置坐标x=[]; y=[]; z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第三层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第四层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第五层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第六层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第七层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第八层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第九层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第十层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第十一层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第十二层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔第十三层，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1986年古塔塔顶，中心位置坐标x=[];y=[];z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmean xmean=ymean=zmean= 附件3%1996年古塔第一层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第二层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第三层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第四层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第五层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第六层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第七层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第八层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第九层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第十层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第十一层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第十二层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔第十三层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%1996年古塔塔顶，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=附件4%2009年古塔第一层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第二层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第三层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第四层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第五层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第六层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第七层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第八层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第九层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第十层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第十一层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第十二层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔第十三层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2009年古塔塔顶，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=附件5%2011年古塔第一层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2011年古塔第二层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zmeanxmean=ymean=zmean=%2011年古塔第三层，中心位置坐标>>x=[];>>y=[];>>z=[];>>xmean=mean(x);>>ymean=mean(y);>>zmean=mean(z);>>xmean,ymean,zme