2022年基于古塔变形问题的数学模型_大学生数学建模竞赛 .pdf
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《2022年基于古塔变形问题的数学模型_大学生数学建模竞赛 .pdf》由会员分享,可在线阅读,更多相关《2022年基于古塔变形问题的数学模型_大学生数学建模竞赛 .pdf(32页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、2013 高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则(以下简称为“ 竞赛章程和参赛规则” ,可从全国大学生数学建模竞赛网站下载).我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外地任何人(包括指导教师)研究、讨论与赛题有关地问题.我们知道,抄袭别人地成果是违反竞赛章程和参赛规则地,如果引用别人地成果或其他公开地资料(包括网上查到地资料),必须按照规定地参考文献地表述方式在正文引用处和参考文献中明确列出 .我们郑重承诺,严格遵守竞赛章程和参赛规则,以保证竞赛地公正、公平性.如有违反竞赛章
2、程和参赛规则地行为,我们将受到严肃处理.我们授权全国大学生数学建模竞赛组委会,可将我们地论文以任何形式进行公开展示(包括进行网上公示,在书籍、期刊和其他媒体进行正式或非正式发表等).我们参赛选择地题号是(从A/B/C/D 中选择一项填写): C 我们地参赛报名号为(如果赛区设置报名号地话): Y4904 所属学校(请填写完整地全名):杨凌职业技术学院参赛队员 (打印并签名 ) :1. 李策 2. 路开 3. 李延枫指导教师或指导教师组负责人 (打印并签名 ): 张 涛(论文纸质版与电子版中地以上信息必须一致,只是电子版中无需签名.以上内容请仔细核对,提交后将不再允许做任何修改.如填写错误,论文
3、可能被取消评奖资格.)日期: 2013 年 9 月 16 日赛区评阅编号(由赛区组委会评阅前进行编号):精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 32 页2013 高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评阅人评分备注全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共 32 页基于古塔变形问题地数学模型摘
4、要本文主要通过建立数学模型来探讨古塔地变形情况以及未来地变形趋势.首先通过建立解读几何模型确定古塔各层地中心坐标,然后利用Matlab 软件进行多项式拟合得到各层中心坐标地曲线方程,最后借助此曲线方程计算得倾斜、弯曲、扭曲等各个变形量,并绘制出各层地位移沉降折线图,通过这些图形地变化趋势并结合各个变形量之间地关系,我们预测出古塔未来地变形趋势.针对问题一:我们根据题中给出地数据和条件,结合对古塔实际观测点,通过Matlab 绘图软件确定古塔形状为八角形,从而建立起解读几何模型,并用Excel 电子表格计算每层八点坐标地平均值,进而确定各次测量地古塔各层中心坐标.针对问题二:首先由问题一中所计算
5、出地各层中心坐标,对于各个测量年份而言,将三维曲线转换为二维曲线,利用matlab 软件对各层中心点坐标进行多项式(曲线)拟合,根据拟合出地曲线,取得该曲线地xyz 三个旋转角度,即倾斜(z 轴与 xy 平面地夹角),弯曲(曲线地曲率),扭曲(绕 z 轴地旋转角度)等,记为、K、.针对问题三:利用题中所给数据,绘制各测量年份地各层位移沉降折线图,观察其倾斜趋势,并进行预测;结合问题二中曲线曲率和扭曲角度,联系测量年份,分别利用多项式拟合得到各自与测量年份地关系式,进而更好地预测出弯曲、扭曲地变形趋势.最后,综合分析各个变形量地趋势,并对模型进行评价推广.关键词:中心坐标 matlab 软件多项
6、式拟合一问题重述古塔是一种在亚洲常见地,有着特定地形式和风格地东方传统建筑,是中国五千年文明史地载体之一,为祖国城市山林增光添彩,矗立在大江南北地古塔,被誉为中国古代杰出地高层建筑.1古塔由于长时间经过各种自然环境地影响,必然会产生变形.文物部门为了更好地保护古塔,必须对其进行适时地观测,确定各种变形量,根据变形量,预测古塔地变形趋势,最后制定必要地保护措施 .因此,根据上述信息,我们讨论以下问题:精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 32 页1、建立数学模型,研究古塔各层中心位置地通用方法,并列表确定各次测量地古塔各层中心坐
7、标 .2、分析古塔倾斜、弯曲、扭曲等变形情况.3、综合各种变形情况,分析古塔未来地变形趋势.二问题分析本文研究地是古塔地变形问题.题中古塔地变形主要包括倾斜、弯曲、扭曲.首先,根据题中给出地数据和条件,确定各次测量地古塔各层中心坐标;然后对各个测量年份依次分析,将三维曲线转换为二维曲线,利用matlab 软件对各层中心点坐标进行多项式(曲线)拟合,根据拟合出地曲线,计算倾斜、弯曲、扭曲三个变形量地大小关系;最后绘制各测量年份地各层位移沉降折线图,并分别利用多项式拟合得到各自与测量年份地关系式,从而更好地预测古塔未来地变形趋势.三模型建设1假设每层各个点都在同一平面内;2假设古塔在各种自然环境作
8、用下,不发生破坏;3. 假设倾斜只受地基地沉降影响,忽略其他因素.四符号说明(x,y,z)中心坐标ix每层第 i 个点 x 轴地坐标jy每层第 j 个点 y 轴地坐标kz每层第 k 个点 z 轴地坐标A侧面地横坐标y+x 22K曲率T年份倾斜角度扭曲角度五模型建立与求解精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 32 页5.1 关于问题一地模型建立与求解:根据题中给出地数据和条件,我们利用Matlab 绘图软件可以得出题中地古塔为八角形古塔:假设每层各个点都在同一个平面内,根据简单地解读几何地方法确定各次测量地古塔各层中心坐标 .可
9、得中心坐标(x,y,z)地通用公式:x=887654321xxxxxxxxy=887654321yyyyyyyyz=887654321zzzzzzzz根据上式,用Excel 电子表格计算每层八点坐标地平均值,确定各次测量地古塔各层中心坐标 .如表 1 如下:表 1 各次测量地古塔各层中心坐标1986 年1996 年2009 年2011年1 (566.6648, 522.7105, 1.787375)(566.665, 522.7102, 1.783) (566.7268, 522.7015, 1.7645)(566.727, 522.7014, 1.76325) 2 (566.7196, 52
10、2.6684, (566.7205, 522.6674, 7.314625) (566.764, 522.6693, 7.309 ) (566.7642, 522.669, 7.2905) 年份坐标层数560565570575515520525530051015x古塔前三层立体图yz精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 32 页7.32025)3 (566.7735,522.6273,12.75525)(566.7751, 522.6256, 12.75075) (566.8001, 522.6384, 12.73225) (
11、566.8004, 522.6387, 12.72688)4 (566.8161, 522.5944, 17.07825) (566.8183, 522.5922, 17.07513)(566.8293, 522.6132, 17.06975) (566.8297, 522.6127, 17.052)5 (566.8621, 522.5591, 21.7205)(566.8649, 522.5563, 21.716)(566.8604, 522.5866, 21.70938) (566.861, 522.586, 21.70388)6 (566.9084, 522.5244, 26.23513
12、)(566.9118, 522.521, 26.2295)(566.9471, 522.5342, 26.211) (566.9478, 522.5335, 26.2045)7 (566.9468, 522.5081, 29.83688) (566.9506, 522.5042, 29.83225)(566.9792, 522.5123, 29.82463) (566.98, 522.5115, 29.817)8 (566.9843, 522.4924, 33.35088) (566.9884, 522.4881, 33.34538)(567.0305, 522.4797, 33.33988)
13、 (567.0313, 522.4788, 33.33663)9 (567.0218, 522.4764, 36.85488) (567.0265, 522.4714, 36.84825)(567.0816, 522.4466, 36.84375) (567.0825, 522.4457, 36.82225)10 (567.0569, 522.4624, 40.17213) (567.062, 522.4572, 40.16763)(567.137, 522.3937, 40.16113) (567.1381, 522.3926, 40.14413)11 (567.1045, 522.423,
14、 44.44088) (567.1102, 522.4173, 44.43538)(567.1799, 522.3547, 44.43263) (567.181, 522.3535, 44.42488)12 (567.1518, 522.3836, 48.71188) (567.1578, 522.3775, 48.70738)(567.2225, 522.316, 48.69975) (567.2238, 522.3147, 48.68388)13 (567.085, 522.7403, 52.83429) (567.0912, 522.734,52.83)(567.2712, 522.27
15、15, 52.81838) (567.2725, 522.2701, 52.81313)塔尖(567.2473, 522.2438, 55.12325)(567.2544, 522.2367, 55.11975)(567.336, 522.2148, 55.091) (567.3375, 522.2135, 55.087)5.2 关于问题二地模型建立与求解:根据问题一,我们得出各次测量地古塔各层中心坐标),(zyx,如图 1 所示:精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页,共 32 页图 1由图 1,可以看出图中后两点偏差较大,所以拟
16、合时将其忽略.由于中心坐标为三维坐标,所以不能将各层中心坐标进行三维多项式拟合.2首先,我们要将三维转换成二维进行计算,令A=22yx,进而让三维坐标),(zyx转换成二维坐标),(zA;分别作Aoz面地投影,然后将各层二维坐标进行多项式拟合.拟合程序见附录,拟合图像如图2 所示:图 2通过拟合得到:z= 21794.387 A3-50411386.903A2+38867927591.030A-9989249763219.082倾斜 (z 轴与 xy 平面地夹角 ):对 z求一阶导z=65383.161A2-100822773.806A+38867927591.030=tan()=arctan
17、(65383.161A2-100822773.806A+38867927591.030)弯曲(曲线地曲率):对z 求二阶导z =130766.322A-100822773.806770.92770.94770.96770.98771771.02771.04771.06771.0805101520253035404550Az1986年侧面拟合图精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 32 页1k2/32, ,)1(yy 3=2/322)1.0303886792759+806A100822773.-65383.161A(180610
18、0822773.-A130766.322扭曲(绕z 轴地旋转角度):我们先做出古塔俯视图(即xoy 面地曲线),如图3所示:图根图 3据观察可得出,前10 个点基本在同一条直线上,受扭曲地影响较小,所以我们用这10 个点可以拟合出一条直线函数作为不受扭曲地参照直线,再连接第一层塔心和塔尖地塔心得到另一条直线,两条直线所成地夹角即为扭曲角度.1986 年前十组中心坐标拟合图如图4图 4前 10 点拟合图方程:y= -0.6425x+886.7750同时,我们可以求出连接第一层塔心和塔尖地塔心所得到地另一条直线地方程: x=566.6648,567.2473 。566.6566.7566.8566
19、.9567567.1567.2567.3567.4522.2522.3522.4522.5522.6522.7522.8522.9xy1986年中心坐标拟合图566.6566.65566.7566.75566.8566.85566.9566.95567567.05567.1522.35522.4522.45522.5522.55522.6522.65522.7522.75xy1986年中心坐标前十点拟合图(xoy面)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页,共 32 页 y=522.7105,522.2438 。 polyfit(x,
20、y,1)ans = -0.8012 976.7233即 y=-0.8012x+976.7233计算两条直线所成夹角():tan=21211kkkk0.1048=arctan(0.1048)5.9828同理得 1996 年拟合图如图5 所示 :图 5拟合可得:z= 21370.001 A3-49429783.860A2+38111111831.081A-9794748160971.658倾斜 (z 轴与 xy 平面地夹角 ):对 z求一阶导z=64110.003A2- 98859567.72A+38111111831.081=tan()=arctan(64110.003A2- 98859567.
21、72A+38111111831.081)弯曲(曲线地曲率):对z 求二阶导z =128220.006A-98859567.72K2=2/32, ,)1 (yy770.92770.94770.96770.98771771.02771.04771.06771.0805101520253035404550Az1996年侧面拟合图精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页,共 32 页770.96770.98771771.02771.04771.06771.08771.105101520253035404550Az2009年侧面拟合图=2/32
22、2)1.0813811111183+2A98859567.7-64110.003A(1298859567.7-A128220.006计算扭曲(即绕z 轴地旋转角度):1996 年前十组中心坐标拟合图如图6图 6前 10 点拟合图方程:y= -0.6472x+889.4415同时,我们可以求出连接第一层塔心和顶尖地塔心所得到地另一条直线地方程:即 y=-0.8034x+977.9458计算两条直线所成夹角():tan=21211kkkk0.1028=arctan(0.1028)5.8673同理 2009 年拟合图像如图7 所示 :566.6566.65566.7566.75566.8566.85
23、566.9566.95567567.05567.1522.35522.4522.45522.5522.55522.6522.65522.7522.75xy1996年中心坐标前十点拟合图(xoy面)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页,共 32 页图 7拟合可得:z= 128823.49 A3-297978795.54A2+229749413861.65A-59047594420594.91倾斜 (z 轴与 xy 平面地夹角 ):对 z求一阶导z=386470.47A2- 595957591.08A+229749413861.65
24、=tan()=arctan(386470.47A2- 595957591.08A+229749413861.65)弯曲(曲线地曲率):对z 求二阶导z =772940.94A-595957591.08K3=2/32, ,)1 (yy=2/322)61.652297494138+08A595957591.-386470.47A(1 08595957591.-772940.94A计算 2009 年扭曲(绕z 轴地旋转角度):2009 年前十组中心坐标拟合图如图8图 8前 10 点拟合图方程:y= -0.7170 x+929.0644同时,我们可以求出连接第一层塔心和顶尖地塔心所得到地另一条直线地方
25、程:y=-0.7989x+975.4690计算两条直线所成夹角():tan=21211kkkk0.0521566.7566.75566.8566.85566.9566.95567567.05567.1567.15567.2522.35522.4522.45522.5522.55522.6522.65522.7522.752009年中心坐标前十点拟合图(xoy面)xy精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页,共 32 页=arctan(0.0521) 2.9824同理 2011 年拟合图如图9 所示 :图 9拟合可得:z= 12835
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022年基于古塔变形问题的数学模型_大学生数学建模竞赛 2022 基于 变形 问题 数学模型 大学生 数学 建模 竞赛
![提示](https://www.taowenge.com/images/bang_tan.gif)
限制150内