试题.试卷—--古塔的变形高教杯数学建模c题论文答案.doc

古塔的变形_高教杯数学建模c题论文2013高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： C 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：题目：古塔的变形摘要 古塔是展现中国古代悠久历史的文化载体，是我国古代具有代表性的建筑，它不仅蕴含了丰富的历史信息还见证了古代建筑师的建筑之精妙。本文针对古塔长时间承受自重和外力作用引起的组合变形问题，采用数据处理、几何分析及曲线拟合的方法，利用matlab等数学软件编程、计算，给出了确定古塔各层中心位置的通用方法，根据所给的各次测量数据计算出古塔各层的中心坐标，分析了该塔的倾斜、弯曲、扭曲等变形情况对变形趋势做了预测。 针对问题1，我们首先要根据这几年来四次古塔的数据变化情况，用建模软件MATLAB制作成模型图，用数学建模中拟合的方法来画出塔的基本形状，再确定古塔每层的中心点，建立中心点拟合线方程模型，观察是否有倾斜、扭曲、变形等情况。 针对问题2，古塔倾斜的原因主要与日光照射、地基活动有关。首先朝向阳面的地基水分较少，阴面的地基水分较多，于是1万多吨的塔，开始向水分较多，地基松软的方向倾斜。另外大量的地下水开采影响了地基的稳固，而古塔附近的铁路运输，也会造成震动。而且，受到地基的不均匀沉降、地震、大风等影响，都会有可能倾斜等变形情况。 针对问题3，根据管理部门委托的测绘公司的数据表来看，古塔每年都以很小的角度在偏移，由于各种人为或者自然原因，使得古塔慢慢的倾斜为斜塔，斜塔并不一定都会倒塌，只要塔的重心线（通过重心点所引的垂直线）还在塔的底面积范围内，塔就是安全的。因此纠偏要根据每座塔的具体情况而定。且一般来说，有些塔在倾斜的过程中，原本松软的地基会被渐渐压实，然后与倾斜角度构成新的平衡，便就此稳定下来。关键词:古塔 组合变形 倾斜 弯曲 扭曲 中心位置 变形趋势一 问题复述 古塔是有着特定的形式和风格的东方传统建筑，是中国五千年文明史的载体之一，被誉为中国古代杰出的高层建筑。但由于存在时间久远，受各方面影响，古塔的可能会发生变形。文物部门为了更好地保护古塔，必须对其进行适时的观测，确定各种变形量，根据变形量，预测古塔的变形趋势，最后制定必要的保护措施。请根据附件1提供的4次观测数据，讨论以下问题：1、 给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。2、分析该塔倾斜、弯曲、扭曲等变形情况。3、分析该塔的变形趋势。二 问题的分析 本文研究的是古塔的变形问题。在给出了的一个古塔实例以及相应数据（附件1，该实例古塔的4次观测数据）的条件下，我们需要用建立该古塔的各层中心位置的通用方法，且列表给出各次测量的古塔各层的中心位置。并进一步分析该古塔倾斜、弯曲、扭曲等变形情况，最后分析该古塔的变形趋势。 在问题1中，我们需要确定古塔各层中心位置，并列表给出各次测量的古塔各层中心位置。 首先，我们假设各测量点都是选取得科学合理的位置，都是围绕中心点的，并且同一层测量点大致在同一平面上，由已知数据（附件1）也可以看出它们是大致在同一平面上。那么，我们由已知条件知道每层给出的各测量点的数据，我们通过画三维图形可以看出，那近似于一个的多边形，所以我们就可以把问题转变为求多边形的中心位置，然后记录各个中心位置坐标数据。在问题2中，我们需要分析该塔倾斜、弯曲、扭曲等变形情况。根据问题1中求得的各层的中心位置坐标，然后将各层的中心位置坐标，接连起来，观察它。理论上正常古塔的中心位置坐标连线应该大致是一条垂直于X轴和Y轴平行Z轴的直线。如果中心位置坐标连线，还是直线但不平行于Z轴了，说明该古塔发生了倾斜，如果中心位置坐标连线不是一条直线，而是一条有一个弧度的曲线，那么说明该古塔发生了弯曲，如果中心位置坐标连线不是一条直线，而是一条有多个弧度的曲线，那么说明该古塔发生了扭曲。在问题3中，我们分析该塔的变形趋势时，通过分析每一个中心点的变化趋势，来判断整个塔的变形趋势，与预测古塔以后可能会发生的变形情况。 三 模型假设1、假设倾斜只受地基的沉降影响，忽略其他因素。2、假设各测量点都是选取的都是科学合理的位置。3、假设每层的测量点都是围绕着这一层的中心点。4、假设同一层测量点都大致在同一平面上。5、假设各层测量点构成的几何图形的中心是与这一层的中心位置相重合的。6、假设测量点的位置都是古塔上固定的位置。四 符号说明符号 符号说明 第层中心；层数； 测量点标记； 轴坐标； 轴坐标；轴坐标； 拟合线轴坐标； 拟合线轴坐标； 拟合线轴坐标； 第层中心的轴坐标； 第层中心的轴坐标； 第层中心的轴坐标； 测量点数 时间 时间差 坐标差 坐标差 坐标差 真实值与计算的近似值之差 真实值与计算的近似值之差 方程中的参数 第层的方程参数 五 建立数学模型1、确定古塔各层中心位置的通用方法，并列表给出古塔各层中心坐标根据几何中心计算方法，我们可以得出计算各层中心位置的通用方法。 即：我们把每层的8个测量值，分别求x，y，z的均值，得到的坐标就是每层的中心位置。表1：1986年各层中心位置楼层中心位置坐标x/my/mz/m1566.6648522.71051.78742566.7196522.66847.32023566.7735522.627312.75524566.8161522.594417.07835566.8621522.559121.72056566.9084522.524426.23517566.9467522.508129.83698566.9843522.492433.35099567.0218522.476436.854910567.0569522.423040.172111567.1045522.423044.440912567.1518522.383648.711913567.0850522.740352.8343塔顶567.2473522.243755.1232表1数据分析：我们可以从表中发现，中心位置的x轴坐标，在逐渐增加，而y轴坐标在逐渐减少。而正常的塔，理论而言，x轴和y轴坐标因该是不变的。结论：说明该古塔发生了倾斜变形。表2：1996年各层中心位置层中心位置坐标x/my/mz/m1566.6650522.71021.78302566.7205522.66747.31463566.7751522.625612.75084566.8183522.592217.07515566.8649522.556321.71606566.9118522.521026.22957566.9506522.504229.83228566.9884522.488133.34549567.0265522.471436.848210567.0620522.457240.167611567.1102522.417344.435412567.1578522.377548.707413567.0912522.734052.8300塔顶567.2543522.236655.1198表2数据分析：我们可以从表中发现，中心位置的x轴坐标，在逐渐增加，而y轴坐标在逐渐减少。而正常的塔，理论而言，x轴和y轴坐标因该是不变的。结论：说明该古塔发生了倾斜变形。表3：2009年各层中心位置层中心位置坐标x/my/mz/m1566.7268522.70151.76452566.7640522.66937.30903566.8001522.638412.73234566.8293522.613217.06975566.8603522.586621.70946566.9471522.534226.21107566.9792522.512329.82468567.0305522.479733.33989567.0816522.446636.843810567.1370522.393740.161111567.1799522.354744.432612567.2225522.316048.699813567.2712522.271552.8183塔顶567.336522.214855.091表3数据分析：我们可以从表中发现，中心位置的x轴坐标，在逐渐增加，而y轴坐标在逐渐减少。而正常的塔，理论而言，x轴和y轴坐标因该是不变的。结论：说明该古塔发生了倾斜变形。表4：2011年古塔各层中心位置坐标层中心位置坐标x/my/mz/m1566.7270522.70141.76322566.7462522.66907.29053566.8004522.638712.72694566.8297522.642717.05205566.8610522.586021.70396566.9478522.533526.20457566.9800522.511529.81708567.0313522.478833.33669567.0825522.445736.822210567.1381522.392640.144111567.1810522.353544.424812567.2238522.314748.683813567.2725522.270152.8131塔顶567.3375522.213555.087表4数据分析：我们可以从表中发现，中心位置的x轴坐标，在逐渐增加，而y轴坐标在逐渐减少。而正常的塔，理论而言，x轴和y轴坐标因该是不变的。 结论：说明该古塔发生了倾斜变形。2分析该塔倾斜、弯曲、扭曲等变形情况我们由已知的空间直线的标准点向式方程 整理得直线射影式方程其中；。这样直线可以看作是用这2个方程表示的平面相交的直线，所以可以分别对2个方程进行数据拟合。设表示按拟合方程求得的近似值。一般地，它不同于实测值，两者之差同理可得当Q取最小值时a，b，c，d的值即为方程的系数，即满足下列方程时Q值最小。有方程组（6）可写成根据n组数据点解方程组就可以得到a，b，c，d的值，也就可以得到拟合线方程。然后我们就可以通过中心位置点拟合线和拟合线方程来判断，古塔是否已经发生倾斜，弯曲或扭曲。 六 模型的求解（1）通过前的算法模型，我们建立出1986年中心位置拟合线方程为：图1 1986年古塔数据综合分析分析1986年古塔变形情况：通过观察图1的第4个子图，我们可以看出中心位置点基本上都是在拟合线上或附近的，说明我们的建立的模型正确的，是可靠性较高的模型。 我们通过分析表1中心点的位置坐标，中心点x坐标有逐层增加的趋势，y坐标有减小的趋势，然后结合拟合线方程和图4中的第1个子图，可知中心拟合线不与z轴平行，有向x轴数值增加方向倾斜。我们再分析图1中的第3个子图，这是一个府视图，我们可以观察到古塔的测量点曲线图的曲线密集程度分布是不均匀的，也说明古塔发生变形。结论：综上分析我们可以得出古塔发生了倾斜变形，往x轴数值增加方向与y轴数值减少方向倾斜。（2）通过前的算法模型，我们可以建立起1996年中心位置拟合线方程为：图2 1996年古塔数据综合分析分析1996年古塔变形情况：通过观察图2的第4个子图，我们可以看出中心位置点基本上都是在拟合线上或附近的，说明我们的建立的模型正确的，是可靠性较高的模型。 我们通过分析表2中心点的位置坐标，中心点x坐标有逐层增加的趋势，y坐标有减小的趋势，然后结合拟合线方程和图4中的第1个子图，可知中心拟合线不与z轴平行，有向x轴数值增加方向倾斜。我们再分析图2中的第3个子图，这是一个府视图，我们可以观察到古塔的测量点曲线图的曲线密集程度分布是不均匀的，也说明古塔发生变形。结论：综上分析我们可以得出古塔发生了倾斜变形，往x轴数值增加方向与y轴数值减少方向倾斜。（3）通过前的算法模型，我们可以建立起2009年中心位置拟合线方程为：图3 2009年古塔数据综合分析分析2009年古塔变形情况：通过观察图3的第4个子图，我们可以看出中心位置点基本上都是在拟合线上或附近的，说明我们的建立的模型正确的，是可靠性较高的模型。 我们通过分析表3中心点的位置坐标，中心点x坐标有逐层增加的趋势，y坐标有减小的趋势，然后结合拟合线方程和图4中的第1个子图，可知中心拟合线不与z轴平行，有向x轴数值增加方向倾斜。我们再分析图3中的第3个子图，这是一个府视图，我们可以观察到古塔的测量点曲线图的曲线密集程度分布是不均匀的，也说明古塔发生变形。结论：综上分析我们可以得出古塔发生了倾斜变形，往x轴数值增加方向与y轴数值减少方向倾斜。（4）通过前的算法模型，我们可以建立2011年中心位置拟合线方程为：图4 2011年古塔数据综合分析分析2011年古塔变形情况： 通过观察图4的第4个子图，我们可以看出中心位置点基本上都是在拟合线上或附近的，说明我们的建立的模型正确的，是可靠性较高的模型。 我们通过分析表4中心点的位置坐标，中心点x坐标有逐层增加的趋势，y坐标有减小的趋势，然后结合拟合线方程和图4中的第1个子图，可知中心拟合线不与z轴平行，有向x轴数值增加方向倾斜。我们再分析图4中的第3个子图，这是一个府视图，我们可以观察到古塔的测量点曲线图的曲线密集程度分布是不均匀的，也说明古塔发生变形。结论：综上分析我们可以得出古塔发生了倾斜变形，往x轴数值增加方向与y轴数值减少方向倾斜。表5：各层两次测量期间中心位置变化情况1996-19862009-19962011-2009变化量10.0002-0.0003-0.00440.0618-0.0087-0.01850.0002-0.0001-0.001320.0009-0.0010-0.00560.04350.0019-0.0056-0.0178-0.0003-0.018530.0016-0.0017-0.00440.02500.0128-0.01850.00030.0003-0.005540.0022-0.0022-0.00320.01100.0210-0.00540.00040.0295-0.017750.0028-0.0028-0.0045-0.00460.0303-0.00660.0007-0.0006-0.005660.0034-0.0034-0.00560.03530.0132-0.01850.0007-0.0007-0.006570.0039-0.0039-0.00470.02860.0081-0.00760.0008-0.0008-0.007680.0041-0.0043-0.00550.0421-0.0084-0.00560.0008-0.0009-0.003290.0047-0.0050-0.00670.0551-0.0248-0.00440.0009-0.0009-0.0216100.00510.0342-0.00450.0750-0.0635-0.00650.0011-0.0011-0.0170110.0057-0.0057-0.00550.0697-0.0626-0.00280.0011-0.0012-0.0078120.0060-0.0061-0.00450.0647-0.0615-0.00760.0013-0.0013-0.0160130.0062-0.0063-0.00430.1800-0.4625-0.01170.0013-0.0014-0.0052塔顶0.0070-0.0071-0.00340.0817-0.0218-0.02880.0015-0.0013-0.0040表5数据分析：首先我们分析z轴上的数据变化，我们从z轴数值变化可知，z轴数值一直都是在减少，说明该古塔的各层中心高度一直在降低，但各层间的减少量不一样，说明不是地基下沉引起的变化，说明该古塔一定有倾斜变形。我们再分析下x轴上的数据变化，我们从x轴数值变化可知，该古塔的中心位置坐标的x轴坐标几乎都在增加，有明显的趋势，说明古塔有往x轴正方向倾斜。我们最后在分析下y轴上的数据变化，y轴坐标在2009至1996年间1到7层中心位置坐标有增加而8到塔顶有减少，说明该时期该古塔在y轴方向上可能有弯曲，但看表5中其它两个期间的数据y轴数值都是在减少说明这两个时间，该古塔又有向y轴数值减少方向倾斜。从表5中我们看不到古塔有在某个方向上有，坐标数值增加了，再减少，再增加，或减少了，再增加，再减少，所以我们可以初步认定该古塔没有发生扭曲。结论：综上分析我们可以大致知道此古塔的倾斜趋势是往x轴正方向即x轴坐标增加方向倾斜，而在y轴方向上有弯曲的迹象，并且有向y轴坐标减少的方向倾斜，没有发生扭曲。表6：中心位置的变化情况1996-19862009-19862011-1986变化量10.0002-0.0003-0.00440.0620-0.0090-0.02290.0622-0.0091-0.024220.0009-0.0010-0.00560.04440.0009-0.01120.02660.0006-0.029730.0016-0.0017-0.00440.02660.0111-0.02290.02690.0114-0.028440.0022-0.0022-0.00320.01320.0188-0.00860.01360.0483-0.026350.0028-0.0028-0.0045-0.00180.0275-0.0111-0.00110.0269-0.016760.0034-0.0034-0.00560.03870.0098-0.02410.03940.0091-0.030670.0039-0.0039-0.00470.03250.0042-0.01230.03330.0034-0.019980.0041-0.0043-0.00550.0462-0.0127-0.01110.0470-0.0136-0.014390.0047-0.0050-0.00670.0598-0.0298-0.01110.0607-0.0307-0.0327100.00510.0342-0.00450.0801-0.0293-0.01100.0812-0.0304-0.0280110.0057-0.0057-0.00550.0754-0.0683-0.00830.0765-0.0695-0.0161120.0060-0.0061-0.00450.0707-0.0676-0.01210.0720-0.0689-0.0281130.0062-0.0063-0.00430.1862-0.4688-0.01600.1875-0.4702-0.0212塔顶0.0070-0.0071-0.00340.0887-0.0289-0.03220.0902-0.0302-0.0362表6数据分析：首先我们分析z轴上的数据变化，我们从z轴数值变化可知，z轴数值一直都是在减少，说明该古塔的各层中心高度一直在降低，但各层间的减少量不一样，说明不是地基下沉引起的变化，说明该古塔一定有倾斜变形。我们再分析下x轴上的数据变化，我们从x轴数值变化可知，该古塔的中心位置坐标的x轴坐标几乎都在增加，有明显的趋势，说明古塔有往x轴正方向倾斜。我们最后在分析下y轴上的数据变化，y轴坐标在这4次测量间1到7层中心位置坐标有增加而8到塔顶有减少，说明该时期该古塔在y轴方向上可能有弯曲。 从表6中我们没有发现古塔有在某个方向上有，坐标数值有增加趋势了，再减少，再增加，或有减少趋势了，再增加，再减少，所以我们可以初步认定该古塔没有发生扭曲。结论：综上分析我们可以大致知道此古塔的倾斜趋势是往x轴正方向即x轴坐标增加方向倾斜，而在y轴方向上有弯曲的迹象，并且有向y轴坐标减少的方向倾斜，没有发生扭曲。综上6次数据分析可得结论：该古塔发生了倾斜和弯曲，古塔的倾斜趋势是往x轴正方向即x轴坐标增加方向和y轴坐标的负方向即坐标减少方向倾斜，而在y轴方向上有弯曲，目前该古塔还没有发生扭曲变形。3 分析该塔的变形趋势通过4次测量点的数据，我们建立了一个拟合线趋势分析模型，来预测该古塔的变形趋势，因为从上面的分析我们可知，该古塔没有地基下沉现象，所以我们可知道塔的高度坐标是一个因变量，所以不在考虑：模型正确定性与可用性检验：图5第1层到第6层的中心位置移动趋势拟合图图6 第7层到塔顶的中心位置移动趋势拟合图第1、2、3、6、8、9、10、11、12、13和塔顶的拟合线模型都是与中心位置点非常吻合，说明这些层的模型是基本符合要求的，具有可靠性的，说明此模型是可用的。根据此模型（模型方程11），我们通过分析计算得到表7中每一层的参数数据：表7：各层中心位置变形趋势模型参数参数层abcd10.0044-0.0428-0.00060.005820.0023-0.02050.0001-0.002230.0018-0.01590.0009-0.010740.0008-0.00560.0027-0.03135-0.00030.00560.0021-0.023560.0025-0.02140.0009-0.012270.0020-0.01630.0005-0.009180.0030-0.0255-0.00060.002090.0039-0.0340-0.00180.0126100.0053-0.0475-0.00450.0788110.0050-0.0432-0.00450.0383120.0383-0.0395-0.00440.0372130.0127-0.1193-0.03260.3151塔顶0.0058-0.0504-0.00160.0087表7数据分析：我们把表7中的数据带入对应的参数，就可以得到每一层的中心位置的x和y的坐标随时间的变化情况的模型。分析参数a，a几乎都是正数，所以随着时间的增加，中心位置的x坐标将越来越大，这符合了我们问题2中分析的古塔变形情况，古塔会向x轴数值增大的方向倾斜。分析参数c，c在1到7层为正数，而在8到塔顶，为负数，也验证问题2中分析的古塔变形情况，古塔会向y轴数值减少方向倾斜，和古塔会发生弯曲，弯曲点在7和8层之间。我们也可以根据此模型方程结合模型方程（2）可以计算出z轴的坐标值。结论：该模型是可用的，具有可靠性，和扩展性。通过该模型（模型方程11）我们可以预测以后的时间，该古塔的中心位置坐标是怎么变化的，会变化多少，我们都可以通过此模型计算得到。通过此模型分析该古塔我们得到该古塔会向x轴数值增大的方向逐渐倾斜，同时会向y轴数值减少方向逐渐倾斜，和古塔会逐渐发生弯曲，弯曲点在7和8层之间。七 模型的评价1 、模型的优点（1）运用的模型简单易懂，有很好的实际指导意义。（2）运用表格和图像相结合，对于结果的分析更加清晰。（3）数学软件Matlab和Excel软件的运用提高了结果的可行度，数据更加精确。2、 模型的缺点（1）本题对数据依赖性比较大，只是对题中所给数据做了一个理想化的模型可能与实际不相吻合。（2）我们建立的拟合模型也只是一个预测模型，在实际生活中会有各方面的影响，因此通过我们所建立的模型只能大致体现一个发展趋势而无法精确地描述其发展趋势。 八 模型的改进1、尽量降低对模型的依赖性，提高所做的模型与实际的吻合程度。2、在建立拟合模型的时候要尽多的考虑到一些实际生活的各方面因素对其产生的影响。3、只能把所做的模型当做参考，而不能过多的依赖模型所得的结果。九 参考文献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年。十 附录附件2%1986年古塔第一层，中心位置坐标 x=565.454 562.058 561.39 563.782 567.941 571.255 571.938 569.5; y=528.012 525.544 521.447 518.108 517.407 519.857 523.953 527.356; z=1.792 1.818 1.783 1.769 1.772 1.77 1.794 1.801; xmean=mean(x); ymean=mean(y); zmean=mean(z); xmean,ymean,zmeanxmean = 566.6648ymean = 522.7105zmean = 1.7874%1986年古塔第二层，中心位置坐标x=565.48 562.238 561.663 564.001 567.995 571.165 571.801 569.414;y=527.764 525.364 521.42 518.226 517.563 519.961 523.908 527.141;z=7.326 7.351 7.314 7.301 7.306 7.304 7.324 7.336; xmean=mean(x); ymean=mean(y); zmean=mean(z); xmean,ymean,zmeanxmean = 566.7196ymean = 522.6684zmean = 7.3202%1986年古塔第三层，中心位置坐标x=565.506 562.415 561.931 564.216 568.048 571.076 571.666 569.33;y=527.52 525.188 521.394 518.343 517.716 520.063 523.864 526.93;z=12.761 12.786 12.749 12.736 12.741 12.74 12.758 12.771; xmean=mean(x); ymean=mean(y); zmean=mean(z); xmean,ymean,zmeanxmean = 566.7735ymean = 522.6273zmean = 12.7552%1986年古塔第四层，中心位置坐标x=565.526 562.555 562.144 564.387 568.091 571.005 571.558 569.263;y=527.327 525.047 521.373 518.435 517.838 520.144 523.829 526.762;z=17.084 17.109 17.072 17.059 17.064 17.063 17.081 17.094; xmean=mean(x); ymean=mean(y); zmean=mean(z); xmean,ymean,zmeanxmean = 566.8161ymean = 522.5944zmean = 17.0783%1986年古塔第五层，中心位置坐标x=565.548 562.706 562.373 564.571 568.136 570.929 571.443 569.191;y=527.119 524.896 521.351 518.534 517.969 520.232 523.791 526.581;z=21.726 21.751 21.714 21.701 21.705 21.708 21.723 21.736; xmean=mean(x); ymean=mean(y); zmean=mean(z); xmean,ymean,zmeanxmean = 566.8621ymean = 522.5591