2013高社杯全国大学生数学建模竞赛C题

2013高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：C我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名)：1.2.3.指导教师或指导教师组负责人(打印并签名)：日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：1古塔的变形摘要：本文研究的古塔的变形问题，通过对问题背景及附件资料进行深入地分析，采用数据拟合、求平均值等方法整理出具有科学性的分析数据。通过对建筑物位移监测数据处理方法的研究,采用自回归模型对位移监测数据进行处理,根据建立的模型对具体建筑物的监测点的位移变化量进行预报。经过计算分析,根据位移量之间变化的关系而建立的自回归预测模型具备较高的拟合及预测精度,运用三维坐标系和数学软件将古塔的模型以空间模型的形式表现出来，直观且科学，对于研究古塔的变形具有较高的科学性和说服性。再通过三维坐标之间的回归和三维坐标与时间的回归而分析出古塔的倾斜，弯曲，扭曲等变形状况，通过数学软件的计算及列表列图的方法将结果直观体现，通过大量的计算与分析，运用几何和代数方法将古塔的变形量以数学的方式说明。对于分析古塔变形趋势中，运用了位移差和位移残差平方公式等量及与时间的关系来说明其变形趋势。对于问题一，通过对监测数据的分析，得出此塔为八边形的塔，并通过平均值法求出古塔各层的中心坐标，具体见表（一）。对于问题二，通过问题一对变形监测数据的研究和处理，我们组运用了自回归模型的方法，利用Z和X,Y之间的回归关系，我们运用数学软件计算求出a1和a2，并通过代数及其几何关系，求出每年监测出的古塔倾斜角度。具体见图（2）和表（2）。对于问题中的弯曲问题，我们用古塔的中点高度发生变化的多少来表示弯曲程度，由于弯曲程度主要是随着时间的变化而变化，所以我们用时间和变形的监测数据进行回归拟合得出△Z（表示古塔的弯曲程度），经过数学软件的大量计算，用列表的方式将每次监测所得每层的弯曲程度表示出来，具体见表（3）。对于问题中的扭曲问题，我们用第一层作为基层，即不扭动层，其他层相对于一层扭动了多少度来说明古塔的扭曲程度。为了这一说明，我们取每一层第一点和第五点作直线。并将每一层的直线与第一层的直线做对比，求出两直线之间的夹角，并用此夹角来说明除第一层以外的每一层相对于第一层扭动了多少度。以此角度大小来表现古塔扭曲程度。具体见表（4）对于问题三，是分析古塔的变形趋势，根据监测的变形数据和位移与时间及波动稳定性的关系，列出了一个时期对上一时期的位移差，看出古塔变形的趋势是向哪个方向，具体见图（5）。再通过位移平方差公式，对古塔整体变形的趋势波动进行分析，并结合时间等因素推算出古塔变形趋势在增加。具体见表（5）。关键词：变形数据拟合平均值自回归模型位移变化三维坐标系数学软件几何代数方法位移差位移平方差2一、问题重述1.1基本情况：某古塔在我国已有上千年的历史，是我国的重点保护文物。但由于长时间受到自重、气温、风力、地震、飓风的影响，古塔出现了各种变形，如倾斜、弯曲、扭曲等。为保护古塔，文物部门需要适时对古塔进行观测，了解各种变形量，以制定必要的保护措施，因此管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该古塔进行了4次观测。1.2需解决的问题:问题一：根据附件1提供的4次观测数据，给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。问题二：利用数学模型和所得数据分析该塔倾斜、弯曲、扭曲等变形情况。问题三：分析该塔的变形趋势。二、问题分析2.1问题一的分析：对于问题一，通过观察所给的数据，每一层的8个观测点都位于古塔的每层的八个角落。据此，将古塔的每一层类似看做一个正八边形，可用求正八边形的中点的方法来确定古塔每一层的中点坐标。通过所给每层的观测数据，分别求X，Y，Z坐标的平均值，近似作为古塔各层中心坐标。2.2问题二的分析：对于问题二，主要是研究古塔的倾斜，弯曲，扭曲等变形量，据此分三点来分析问题二。对于倾斜，主要是三维坐标中X，Y轴对于Z轴的倾斜角度，即是对中心轴的倾斜角度，而第一题已经对每年的监测数据进行了总结，得出每次监测各层3的中心坐标，故运用X，Y与Z的自回归模型，再通过代数和几何关系，求出每年古塔X，Y轴的倾斜角度，并加以说明，具体结果见表（2）。对于弯曲，主要是要找出能表示弯曲程度的量，经过查找变形量具体分析的资料，知道弯曲主要是中点的高度（Z）发生了变化，列出时间和每层的中心坐标的回归拟合，并且通过数学软件计算，得出弯曲程度△Z，并列出表直观说明，易于比较。具体结果见表（3）。对于扭曲变量，通过matlab软件做出的古塔空间图形，及查找了扭曲问题的分析的资料，我们最终采用以第一层为基层，其他层相对于古塔扭曲的度数，主要是通过每层的两个点算出每层的斜率，再通过斜率公式算出角度，以此说明古塔每次扭曲程度的变化及大小。具体见表（4）2.3问题三的分析：对于该塔的变形趋势，经过第一题和第二题，已将古塔的变形位置的中心和变形量进行了具体分析。本题要分析变形趋势，具体是从古塔各中心点在各个时期的位移变化量，通过列位移残差公式，得出该塔的变形趋势。见图（5）。再通过整体的位移残差平方公式，得到古塔每层在某个时期的变形趋势波动，具体见表（5）。三、问题假设１．假设该塔为正八边形的塔２．假设该塔的底层不扭曲变形３．假设监测该塔时后两次改变了监测位置４．假设地质的变动对古塔无影响四、建立模型及求解（一）符号说明：1∠A1为古塔与X轴构成的角度42∠A2为古塔与Y轴构成的角度3Z为中点的高度4△Z来表示古塔的弯曲程度5Bi为表示两直线的夹角6△X为X轴的位置偏移量7△Y为Y轴的位置偏移量8△Z1为Z轴的位置偏移量9△S2为古塔整体的位置偏移平方差10i为古塔层数11t为从1886年开始所经过的时间（二）模型建立及求解3.1对问题1的求解根据数据中给出的各层各个点坐标的值，可把古塔的每一层类似的看成一个正八边形。根据数学逻辑和模型，正八边形的中心坐标可用平均值法求出，故塔的中心坐标为（ZYX，，）。附件一是每次测量出的古塔的各层的各个点的坐标，通过平均值法整理数据，得出以下结果：表（一）i19861996200920111X566.8377566.665566.7268566.727Y522.7105522.7102522.7015522.7014Z1.7873751.7831.76451.763252X566.7196566.7205566.764566.7642Y522.6684522.6674522.6693522.669Z7.320257.3146257.3097.290553X566.7735566.7751566.8001566.8004Y522.6273522.6256522.6384522.6387Z12.7552512.7507512.7322512.726884X566.8161566.8183566.8293566.8297Y522.5944522.5922522.6132522.6127Z17.0782517.0751317.0697517.0525X566.8621566.8649566.8604566.861Y522.5591522.5563522.5866522.586Z21.720521.71621.7093821.703886X566.9084566.9118566.9471566.9478Y522.5244522.521522.5342522.5335Z26.2351326.229526.21126.20457X566.9468566.9506566.9792566.98Y522.5081522.5042522.5123522.5115Z29.8368829.8322529.8246329.8178X566.9843566.9884567.0305567.0313Y522.4924522.4881522.4797522.4788Z33.3508833.3453833.3398833.336639X567.0218567.0265567.0816567.0825Y522.4764522.4714522.4466522.4457Z36.8548836.8482536.8437536.8222510X567.0569567.062567.137567.1381Y522.4624522.4572522.3937522.39266Z40.1721340.1676340.1611340.1441311X567.1045567.1102567.1799567.181Y522.423522.4173522.3547522.3535Z44.4408844.4353844.4326344.4248812X567.1518567.1578567.2225567.2238Y522.3836522.3775522.316522.3147Z48.7118848.7073848.6997548.6838813X567.085567.0912567.2712567.2725Y522.7403522.734522.2715522.2701Z52.8342952.8352.8183852.81313塔尖X567.2473567.2544567.336567.3375Y522.2438522.2367522.2148522.2135Z55.1232555.1197555.09155.087其中的X=(x1+x2+……+x8)/8得来，Y，Z值同理。通过matlab做出古塔的模型图片：下图（1）7图（1）3.2对问题二的求解：对于古塔的倾斜，弯曲、扭曲等变形情况的分析，第一题通过科学分析和数学逻辑思维，用平均值法对建筑物位置的监测数据进行总结，得出各个时期各个层的中心坐标值，通过这些数据，运用过对建筑物位移监测数据处理方法的研究,采用自回归模型对位移监测数据进行处理,根据建立的模型对具体建筑物的监测点的位移变化量进行预报。经过计算分析,对古塔的变形情况进行说明。3.2.1关于古塔的倾斜问题，我们建立了空间直角坐标（如上图1）对此进行了讨论，发现古塔与X轴构成的角度为∠A1,与Y轴构成的角为∠A2。因此建立了Zi,t与Xi,t、Yi,t的自回归模型如下式○1,且∠A1=arctana1，∠A2=arctana2。8Zi,t=a0+a1*Xi,t+a2*Yi,t○1从问题一中所求得的平均值可得出下图2：1986(1)1996（11）XYZXYZ566.8377522.71051.787375566.665522.71021.783566.7196522.66847.32025566.7205522.66747.314625566.7735522.627312.75525566.7751522.625612.75075566.8161522.594417.07825566.8183522.592217.07513566.8621522.559121.7205566.8649522.556321.716566.908452