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

1、高社杯全国大学生数学建模竞赛C题2013高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： C

2、我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：古塔的变形摘要： 本文研究的古塔的变形问题，通过对问题背景及附件资料进行深入地分析，采用数据拟合、求平均值等方法整理出具有科学性

3、的分析数据。通过对建筑物位移监测数据处理方法的研究, 采用自回归模型对位移监测数据进行处理, 根据建立的模型对具体建筑物的监测点的位移变化量进行预报。经过计算分析, 根据位移量之间变化的关系而建立的自回归预测模型具备较高的拟合及预测精度,运用三维坐标系和数学软件将古塔的模型以空间模型的形式表现出来，直观且科学，对于研究古塔的变形具有较高的科学性和说服性。再通过三维坐标之间的回归和三维坐标与时间的回归而分析出古塔的倾斜，弯曲，扭曲等变形状况，通过数学软件的计算及列表列图的方法将结果直观体现，通过大量的计算与分析，运用几何和代数方法将古塔的变形量以数学的方式说明。对于分析古塔变形趋势中，运用了位移

4、差和位移残差平方公式等量及与时间的关系来说明其变形趋势。 对于问题一，通过对监测数据的分析，得出此塔为八边形的塔，并通过平均值法求出古塔各层的中心坐标，具体见表（一）。 对于问题二，通过问题一对变形监测数据的研究和处理，我们组运用了自回归模型的方法，利用Z和X,Y之间的回归关系，我们运用数学软件计算求出a1和a2，并通过代数及其几何关系，求出每年监测出的古塔倾斜角度。具体见图（2）和表（2）。对于问题中的弯曲问题，我们用古塔的中点高度发生变化的多少来表示弯曲程度，由于弯曲程度主要是随着时间的变化而变化，所以我们用时间和变形的监测数据进行回归拟合得出Z（表示古塔的弯曲程度），经过数学软件的大量计

5、算，用列表的方式将每次监测所得每层的弯曲程度表示出来，具体见表（3）。对于问题中的扭曲问题，我们用第一层作为基层，即不扭动层，其他层相对于一层扭动了多少度来说明古塔的扭曲程度。为了这一说明，我们取每一层第一点和第五点作直线。并将每一层的直线与第一层的直线做对比，求出两直线之间的夹角，并用此夹角来说明除第一层以外的每一层相对于第一层扭动了多少度。以此角度大小来表现古塔扭曲程度。具体见表（4）对于问题三，是分析古塔的变形趋势，根据监测的变形数据和位移与时间及波动稳定性的关系，列出了一个时期对上一时期的位移差，看出古塔变形的趋势是向哪个方向，具体见图（5）。再通过位移平方差公式，对古塔整体变形的趋势

6、波动进行分析，并结合时间等因素推算出古塔变形趋势在增加。具体见表（5）。关键词：变形 数据拟合 平均值 自回归模型 位移变化三维坐标系 数学软件 几何 代数方法 位移差 位移平方差一、问题重述1.1基本情况：某古塔在我国已有上千年的历史，是我国的重点保护文物。但由于长时间受到自重、气温、风力、地震、飓风的影响，古塔出现了各种变形，如倾斜、弯曲、扭曲等。为保护古塔，文物部门需要适时对古塔进行观测，了解各种变形量，以制定必要的保护措施，因此管理部门委托测绘公司先后于1986年7月、1996年8月、2009年3月和2011年3月对该古塔进行了4次观测。1.2需解决的问题:问题一：根据附件1提供的4次

7、观测数据，给出确定古塔各层中心位置的通用方法，并列表给出各次测量的古塔各层中心坐标。问题二：利用数学模型和所得数据分析该塔倾斜、弯曲、扭曲等变形情况。问题三：分析该塔的变形趋势。二、问题分析 2.1问题一的分析： 对于问题一，通过观察所给的数据，每一层的8个观测点都位于古塔的每层的八个角落。据此，将古塔的每一层类似看做一个正八边形，可用求正八边形的中点的方法来确定古塔每一层的中点坐标。通过所给每层的观测数据，分别求X，Y，Z坐标的平均值，近似作为古塔各层中心坐标。2.2问题二的分析：对于问题二，主要是研究古塔的倾斜，弯曲，扭曲等变形量，据此分三点来分析问题二。 对于倾斜，主要是三维坐标中X，Y

8、轴对于Z轴的倾斜角度，即是对中心轴的倾斜角度，而第一题已经对每年的监测数据进行了总结，得出每次监测各层的中心坐标，故运用X，Y与Z的自回归模型，再通过代数和几何关系，求出每年古塔X，Y轴的倾斜角度，并加以说明，具体结果见表（2）。 对于弯曲，主要是要找出能表示弯曲程度的量，经过查找变形量具体分析的资料，知道弯曲主要是中点的高度（Z）发生了变化，列出时间和每层的中心坐标的回归拟合，并且通过数学软件计算，得出弯曲程度Z，并列出表直观说明，易于比较。具体结果见表（3）。对于扭曲变量，通过matlab软件做出的古塔空间图形，及查找了扭曲问题的分析的资料，我们最终采用以第一层为基层，其他层相对于古塔扭曲

9、的度数，主要是通过每层的两个点算出每层的斜率，再通过斜率公式算出角度，以此说明古塔每次扭曲程度的变化及大小。具体见表（4）2.3问题三的分析： 对于该塔的变形趋势，经过第一题和第二题，已将古塔的变形位置的中心和变形量进行了具体分析。本题要分析变形趋势，具体是从古塔各中心点在各个时期的位移变化量，通过列位移残差公式，得出该塔的变形趋势。见图（5）。再通过整体的位移残差平方公式，得到古塔每层在某个时期的变形趋势波动，具体见表（5）。三、问题假设1 假设该塔为正八边形的塔2 假设该塔的底层不扭曲变形3 假设监测该塔时后两次改变了监测位置4 假设地质的变动对古塔无影响四、建立模型及求解（一）符号说明：

10、1A为古塔与X轴构成的角度2A为古塔与Y轴构成的角度3 Z为中点的高度4 Z来表示古塔的弯曲程度5 Bi为表示两直线的夹角6X为X轴的位置偏移量7Y为Y轴的位置偏移量8Z1为Z轴的位置偏移量9S2为古塔整体的位置偏移平方差10 i为古塔层数11 t为从1886年开始所经过的时间（二）模型建立及求解3.1 对问题1的求解 根据数据中给出的各层各个点坐标的值，可把古塔的每一层类似的看成一个正八边形。根据数学逻辑和模型，正八边形的中心坐标可用平均值法求出，故塔的中心坐标为（）。 附件一是每次测量出的古塔的各层的各个点的坐标，通过平均值法整理数据，得出以下结果：表（一）i198619962009201

11、11X566.8377566.665566.7268566.727Y522.7105522.7102522.7015522.7014Z1.7873751.7831.76451.763252X566.7196566.7205566.764566.7642Y522.6684522.6674522.6693522.669Z7.320257.3146257.3097.29053X566.7735566.7751566.8001566.8004Y522.6273522.6256522.6384522.6387Z12.7552512.7507512.7322512.726884X566.8161566.8

12、183566.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.97925

13、66.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.13

14、81Y522.4624522.4572522.3937522.3926Z40.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.272

15、5Y522.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）图（1）3.2对问题二的求解：对于古塔的倾斜，弯曲、扭曲等变形情况的分析，第一题通过科学分析和数学逻辑思维，用平均值法对建筑物位置的监测数据进行总结，得出各个时期各个层的中心坐标值，

16、通过这些数据，运用过对建筑物位移监测数据处理方法的研究, 采用自回归模型对位移监测数据进行处理, 根据建立的模型对具体建筑物的监测点的位移变化量进行预报。经过计算分析, 对古塔的变形情况进行说明。3.2.1关于古塔的倾斜问题，我们建立了空间直角坐标（如上图1）对此进行了讨论，发现古塔与X轴构成的角度为A,与Y轴构成的角为A。因此建立了i,t与i,t、i,t的自回归模型如下式,且A=arctan a1 ， A=arctan a2 。i,t=a0+a1*i,t+a2*i,t 从问题一中所求得的平均值可得出下图2：1986(1)1996（11）XYZXYZ566.8377522.71051.7873

17、75566.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.9084522.524426.23513566.9118522.52126.2295566.9468522.508129.83688566.9506522.504229.8

18、3225566.9843522.492433.35088566.9884522.488133.34538567.0218522.476436.85488567.0265522.471436.84825567.0569522.462440.17213567.062522.457240.16763567.1045522.42344.44088567.1102522.417344.43538567.1518522.383648.71188567.1578522.377548.70738567.085522.740352.83429567.0912522.73452.83567.2473522.243

19、855.12325567.2544522.236755.11975a2a1a0a2a1a013.15844112.1343-70421.432.12952115.3736-82169.82009（24）（112011(26)XYZXZY566.7268522.70151.7645566.727522.70141.76325566.764522.66937.309566.7642522.6697.2905566.8001522.638412.73225566.8004522.638712.72688566.8293522.613217.06975566.8297522.612717.052566

20、.8604522.586621.70938566.861522.58621.70388566.9471522.534226.211566.9478522.533526.2045566.9792522.512329.82463566.98522.511529.817567.0305522.479733.33988567.0313522.478833.33663567.0816522.446636.84375567.0825522.445736.82225567.137522.393740.16113567.1381522.392640.14413567.1799522.354744.432635

21、67.181522.353544.42488567.2225522.31648.69975567.2238522.314748.68388567.2712522.271552.81838567.2725522.270152.81313567.336522.214855.091567.3375522.213555.087a2a1a0a2a1a010.6258792.54033-57992.79.61765591.59094-56927.7 图2因此，在所观测的4年里古塔对于X轴和Y轴的倾斜角度如下表（2）1986(1)1996(11)2009(24)2011(26)89.4890589.5034

22、89.38189.3744685.6540589.99984.623684.064表（2）由分析和计算得出每次监测时对中心轴的倾斜角度都有变大的趋势，古塔的变形不容小觑。3.2.2 弯曲关于古塔的弯曲，我们考虑到古塔弯曲后，中点的高度会发生改变，即Z发生了变化，故我们用Z来表示古塔的弯曲程度，即= -Z =+*t具体计算值如下表（3）1986(1)1996(11)层数ZZ11.7873751.7889560.0015811.7831.779195-0.003827.320257.322030.001787.3146257.312763-0.00186312.7552512.757890.002

23、6412.7507512.74644-0.00431417.0782517.08050.0022517.0751317.07242-0.00271521.720521.721170.0006721.71621.71515-0.00085626.2351326.237640.0025126.229526.2255-0.004729.8368829.837940.0010629.8322529.83087-0.00138833.3508833.35086-2E-0533.3453833.345570.00019936.8548836.856590.0017136.8482536.84672-0.0

24、01531040.1721340.174150.0020240.1676340.16526-0.002371144.4408844.440935E-0544.4353844.435770.000391248.7118848.713950.0020748.7073848.70483-0.002551352.8342952.835820.0015352.8352.82763-0.002371455.1232555.127680.0044355.1197555.11221-0.007542009(24)2011(26)层数ZZ11.76451.7657210.0012211.763251.76364

25、80.00039827.3097.300718-0.008287.29057.2988640.008364312.7322512.73154-0.0007112.7268812.729250.00237417.0697517.06191-0.0078417.05217.06030.0083521.7093821.70732-0.0020621.7038821.706120.00224626.21126.20971-0.0012926.204526.207290.00279729.8246329.82168-0.0029529.81729.820270.00327833.3398833.3387

26、-0.0011833.3366333.337640.00101936.8437536.8339-0.0098536.8222536.831930.009681040.1611340.1537-0.0074340.1441340.151920.007791144.4326344.42905-0.0035844.4248844.428020.003141248.6997548.69297-0.0067848.6838848.691140.007261352.8183852.81699-0.0013952.8131352.815350.002221455.09155.09210.001155.08755.0890.002表（3） 表中表示的是古塔每次每层的中心点的变化及弯曲程度，对于1986和1996及2009和2011的比较，前者（即1986和1996）弯曲程度明显加大，这不仅由于时间长达十年的关系，还有各个方面的影响，比如倾斜，风化，地壳运动等对变形位移数据的影响，肯定的是古塔的弯曲程度在很小很小范围内变化，后者（2009与2011）虽没有前者那么明显，但是考虑到时间因素的话，平均每年的弯曲程度的变化还是高于前者的，所以古