2013年数学建模全国二等奖获得论文.doc

古塔变形情况的分析与改进 【摘要】“盛世修古建”随着我国经济实力的不断发展，古建的保护和改善也成为了国家所关注的事。因为自然灾害所带来的影响，使古建发生了不同的形变。本案例研究的是古塔变形的问题，要求是在自然的影响下对古塔的变形进行假设和分析得出倾斜，弯曲，扭曲的情况。 通过对问题的假设及分析求解中，我们所建立的模型简单且改进措施方便，并且能推广到更多的古塔保护问题上，具有很大的优势。 最后我们组员结合对本次数学建模的学习，实践，写出了我们的感想。我们的理解阐述了数学建模的概念，步骤以及我们在此过程中遇到的问题。 【关键词】数据分析 投影 线性规划 函数 对比影响 措施 Ⅰ、问题重述 由于长时间承受自重，气温，风力等各种作用，偶然还要受地震，飓风的影响，古塔会产生各种变形，比如倾斜，弯曲，扭曲等。为保护古塔，文物部门需要适合时对古塔进行观测，了解各种变形量，以制定必要的保护措施。 古塔已有上千年历史，是我国重点保护文物。管理部门委托测绘公司先后于1986年7月，1996年8月，2009年3 月和2011年3月对该塔进行了4次观测。 根据附件1提供的4次观测数据，讨论以下问题； 给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。 分析该塔倾斜，弯曲，扭曲等变形情况。 分析该塔的变形趋势。 Ⅱ、符号说明 ——职工工资总额，元； ——2000-2024年职工工资总额，元； Ⅲ、模型假设 假设第一问中所得的图为八边形 2.假设塔的变形过程中倾斜，弯曲，扭曲是互不影响的 3.假设在倾斜过程中没有弯曲与扭曲，在扭曲过程中没有倾斜与弯曲，在弯曲过程中没有扭曲与倾斜 Ⅳ、问题分析 对于问题是一个数据平均处理问题，它涉及到许多变量及假设。假设附件1给出的每层8个数据均选 在塔的方位所测，并且这8个点恰好能够成一个平面，若要求塔的中心相当于求这几个面的中心，再将其连接所得就是塔的中心。 对于问题二，可认为是三个小问，即分别研究古塔的倾斜，弯曲，扭曲三种情况。我们先采取从整体到部分的情况，倾斜是指基础两端点倾斜方向的沉降差与距离的比值，而弯曲即可看成是将塔中心分别投影在面XOZ与面YOZ上，最后将中心投影在面XOY上即可得到扭曲变形情况。 最后应用线性回归的方法看倾斜，弯曲，扭曲三种变形情况中哪种变形更严重或者任二种或三种都存在，进而我们建立一个通俗易懂的模型来给予解决及改进。 Ⅴ、问题求解 问题一 1.由附件（1）可得此塔的大致形状如下图，其中（1986年用蓝色绘制；1996年用绿色绘制；2009年用黄色绘制；2011年用红色绘制）并将四个图画在同一个空间中，可以更清晰的看出此塔的每年变形情况。从而更有利于我们分析此塔的变形情况。（用MATLAB画图程序见附录（一）） 2.对第一层坐标分析得此图。（用MATLAB画图程序见附录（二） 由上图可知该图形并非为平面八边形结构，而是一个空间不规则的图形，而现在我们假设它为规则规则几何体来求解，由所学知识可知，规则的几何图形他的重心即为中心，由几何重心公式 可得出1986年的中心坐标数据（见如下表格），而用拟合的方法得到的数据（见如下表格）。（用MATLAB程序见附录（三））1986（算数中心） 1986（拟合中心） x y z x y z 566.66475 522.7105 1.787375 566.665 522.709 1.787375 566.71963 522.66838 7.32025 566.7225 522.671 7.32025 566.7735 522.62725 12.75525 566.7786 522.635 12.75525 566.81613 522.59438 17.07825 566.8226 522.606 17.07825 566.86213 522.55913 21.7205 566.8698 522.575 21.7205 566.90838 522.52438 26.23513 566.9185 522.544 26.23513 566.94675 522.50813 29.83688 566.9521 522.527 29.83688 566.98425 522.49238 33.35088 566.9846 522.51 33.35088 567.02175 522.47638 36.85488 567.0176 522.493 36.85488 567.05688 522.46238 40.17213 567.0476 522.48 40.17213 567.1045 522.423 44.44088 567.1013 522.439 44.44088 567.15