古塔的变(胡、罗、钟)形精品.docVIP

承 诺 书 我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛 赛区评阅编号（由赛区组委会评阅前进行编号）： 2013高教社杯全国大学生数学建模竞赛 编 号 专 用 页 赛区评阅编号（由赛区组委会评阅前进行编号）： 赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注 全国统一编号（由赛区组委会送交全国前编号）： 全国评阅编号（由全国组委会评阅前进行编号）： 古塔的变形 摘要 本文研究了古塔变形的问题。由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。 为解决问题1，我们先讨论研究了中古塔各层中心位置的通用方法。为此我们建立了模型一和优化模型二。通过这四年的数据拟合得出未变形前的古塔的中心轴方程（空间直线），利用中心轴求出层高H为4.254491m，进而可求出未变形前的古塔各层中心坐标（见正文表8），，然后再用类似的方法，得用层高可求出测量的这四年各次测量的古塔各层中心坐标（见正文表3）。 为解决问题2，我们首先分析了该塔的倾斜情况并建立了模型三。用未变形前的古塔中心轴与这四年各年的中心轴（有变形）之间的夹角来表示各年的塔身倾斜情况。利用三角函数知识，可求出四年的空间倾斜角度的正切值分别为0.0119471，0.01211，0.014978，0.015015。 接着我们分析了该塔的弯曲情况并建立了模型四以及优化模型五。通过数据拟合以及曲率公式曲率k= ，得出4个年份的曲率依次为0.000424，0.00046，0.000143，0.000145，可以看到古塔弯曲程度先增加，接着有所减少，然后再增加； 对于古塔的扭曲，通过三次多项式拟合，得出测量的这四年挠率依次为0.0000069221，0.00000000000002988，0.000282427666666667，0.000360233333333333。古塔程度；问题3利用塔的倾斜度变化来比较得出塔呈现出递增的倾斜趋势，倾斜度分别由tant1986=0.11947上升至tant19966=0.012114?再上升到tant2009=0.014978接着升到tant2011=0.015015；接着把各个年份的曲率的进行比较、得知塔的弯曲程度先由K1986=0.000424上升至K1996=0.00046下降至K2009=0.000143再上升为K2011=0.000145因此??最后把各年平均挠率进行比较，分析出古塔的扭曲程度先由i1986=0.0000069221减少到i1996=0.00000000000002988再下降为i2009=0.000282427666666667再上升为i2011=0.000360233333333333,因此古塔呈现出倾斜度，弯曲程度和扭曲程度关键词：????数据拟合?斜率??曲率??挠率??MATLAB软件?古塔，是中国五千年文明史的载体之一，古塔为祖国城市山林增光添彩，塔被佛教界人士尊为佛塔。矗立在大江南北的古塔，被誉为中国古代杰出的高层建筑，但由于长时间受各种因素的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解塔变形的情况，以制定必要的保护措施。现管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。经过4个次的测量，为了解古塔的变形情况，首先需要测量出古塔各个塔层的中心坐标来判断塔中心轴线（各个塔层中心坐标的连线）是否发生偏移并且为了方便于以后对古塔的测量，故需要寻求出一个测量的通用方法，以便于以后能够更好的测出。在通过4个年份4次分别对古塔每个塔层的多点观测数据，分析出该古塔在1986年到2011年间倾斜、弯曲、扭曲的变化情况和各个年之间的趋势分析。 ：表示三维空间中的点的坐标 ：表示第层的中心坐标，=1，2，……，13。 h：表示塔的高度 H：表示古塔各层的层高 N：表示塔的层数 tan t：表示古塔的倾斜程度（斜率） i：表示古塔的扭曲程度（挠率） k：表示古塔的弯曲程度（曲率） 问题分析 4.1对于问题一： 先把塔层定义为一层，然后把每年这一层的八个观测进行算术平均处理，接着再对剩下的塔层进行相同处理。随后再对每一层4个年份进行算术平均，得出13个塔层4个年份算术平均的X Y Z数据，然后对每个项的13组数据进行数据拟合，最终拟合出一个相似于塔的中心线的方程组，然后把两次算术平均处理得出的Z代入进入方程组得出X拟合，再把X拟合代入方程