古塔变形
566.7268566.764566.8001566.8293566.8604566.9471566.9792567.0305567.0816567.137522.7015522.6693522.6384522.6132522.5866522.5342522.5123522.4797522.4466522.39371.76457.30912.7322517.0697521.7093826.21129.8246333.3398836.8437540.16113 2011年：
566.727566.7642566.8004566.8297566.861566.9478566.98567.0313567.0825567.1381522.7014522.669522.6387522.6127522.586522.5335522.5115522.4788522.4457522.39261.763257.290512.7268817.05221.7038826.204529.81733.3366336.8222540.14413 散点连线图：
6
古塔变形
将四年的放在一起看：
7
古塔变形
现在大致可以知道1996到2009年这个时间段古塔变形较大，下面进行进一步量化分析： （一）、倾斜程度：
对中心点进行线性拟合，求出拟合直线与竖坐标轴的夹角，夹角的余弦值可作为倾斜程度的标识：
8
古塔变形
即：
nlcosa??|n||l|
xlxl2?yl2?zl2
对各个年份的中心点进行拟合，采用线性插值法得：
9
古塔变形
放大后：
解得的数据为： 年份 1986 1996 2009 2011
有表格可以看出，塔身与地面的夹角从1.560142变化到1.559192，逐渐靠近地面，下面利用Matlab对其拟合分析：
余弦值 0.010654 0.010775 0.011579 0.011604 角度值 正切值 余弦值 0.010654 0.010775 0.011579 0.011604 1.560142 93.85223 1.560021 92.80484 1.559217 86.35388 1.559192 86.17069
10
2013高教社杯全国大学生数学建模竞赛C题论文-冷磊
