古塔变形.docVIP

古塔变形 摘要 根据题目所给出的坐标数据，由于部分数据缺失，我们根据所给的其他数据确定等差数列,更可能的接近实际数据,并复值缺失的数据，将数据全部整理出来。我们用最小二乘法的方法确定中心坐标。由于古塔长时间承受外界的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲，根据中心坐标解释古塔的变化，及分析它的形变趋势。用合理的数学定义来解决问题。 在问题一中，我们先用图形拟合让八个点尽可能在一个平面上，再用最小二乘法拟合找到中心点。用拟合、最小二乘法及MATLAB编程解决问题 在问题二中，根据每个阶层的中心点，链接起来,塔尖与底层中心的水平距离与塔高的比值既塔的倾斜度α;根据弯曲形变,我们定义弯曲率为K，利用参数方程和二次拟合出的空间曲线的曲率，来描述古塔的弯曲率。 总的来说我们用拟合、最小二乘法、空间曲线拟合、MATLAB编程等常用的方法来解决实际问题，具有较好的推广性。 关键词 拟合、最小二乘法、空间曲线拟合 、空间曲线参数方程 1问题重述 由于长时间承受自重、气温、风力等各种作用，偶然还要受地震、飓风的影响，古塔会产生各种变形，诸如倾斜、弯曲、扭曲等。为保护古塔，文物部门需适时对古塔进行观测，了解各种变形量，以制定必要的保护措施。 某古塔已有上千年历史，是我国重点保护文物。管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该塔进行了4次观测。 请你们根据附件1提供的4次观测数据，讨论以下问题： 1. 让我们讨论出古塔各层中心位置的通用方法，根据给出数据，并列表给出各次测量的古塔各层中心坐标。 2. 分析该塔倾斜、弯曲、扭曲等变形情况。 3. 分析该塔的变形趋势。 2模型假设 由于中国古代建筑追求对称、规则，我们假设古塔为八边形。 假设每次测量数据的点都是固定的。 假设古塔只受倾斜、弯曲、扭曲这三个因素影响，不受长时间承受自重、气温、风力、地震、飓风等的影响。 3变量说明 i为观测层数，j为观测点，为中心点坐标 j=(1,2,3,﹒﹒﹒8); i=(1,2,3,﹒﹒﹒13) K为观测次数（1,2,3,4） j=(1,2,3,﹒﹒﹒8); i=(1,2,3,﹒﹒﹒13) ；第k次测量古塔各层中心拟合曲线。 4模型准备 根据题目所给出的坐标数据，由于部分数据缺失，我们根据后俩年所给的其他数据确定前俩年数据，用等差数列,更可能的接近实际数据,并复值缺失的数据，将数据全部整理出来。下面是补充的数据： 13 1 566.308 525.092 52.866 13 1 566.3142 525.0857 52.864 2 564.716 523.616 52.878 2 564.7222 523.6098 52.877 3 564.481 521.521 52.897 3 564.4872 521.5148 52.889 4 565.91 519.893 52.88 4 565.9161 519.8868 52.878 5 569.8678 523.172 52.79 5 569.8691 523.1158 52.793 6 569.701 521.05 52.703 6 569.7072 521.0438 52.696 7 569.897 523.188 52.794 7 569.9032 523.1817 52.789 8 568.582 524.822 52.822 8 568.588 524.8157 52.817 5模型建立 根据1986年所给的数据，用MATLAB编程我们做出这个古塔的空间形状，其他年份的古塔也是这个样子，古塔是一个八边形，如下： 模型分析 拟我们用到数据拟合让八个点尽可能在一个平面上。用最小二乘法合图形确定平面。 平面方程的一般表达式： 因此拟合平面方程为 最小二乘法思想，建立模型 6模型求解 对各层观测点最接近的平面方程取得极小值 求最小值的偏导数为0 利用MATLAB编程，解得每次测量各层的拟合平面系数如图一； 1969年平面拟合系数 1 -0.0008 0.0034 0.472 2 -0.0008 0.0036 5.8872 3 -0.1035 -0.1706 160.1073 4 -0.0844 -0.1391 13