2012西南交通大学数学建模校赛C题-景区灭火.doc

西南交通大学2012年大学生数学建模竞赛题目:B参赛队员1参赛队员2参赛队员3姓名刘童超王枝李若晗学号200857282010470220104693学院数学学院信息学院信息学院专业统计计软计软电话1592894886213548154032****************@******@313033320@西南交通大学教务处西南交通大学实验室及设备管理处西南交通大学数学建模创新实践基地景区灭火的数学模型【摘要】本文采用网格划分的方法,将连续性问题离散化,建立了图论及其相关模型。同时运用MATLAB的图形处理能力进行了三维制图及一维二维插值,运用C++进行了Dijsktra等算法的编程计算,进而合理的解决了问题。第一问中,考虑到等高线的缺失是由于“破损”,我们舍弃了曲线模拟,而采用了一维插值的方法,并用MATLAB给出了插值曲线,并直观的将曲线拟合至原等高线,发现其效果良好。对于插值结果与直观观察的差异,我们给出了误差分析,并解释了原因。第二问中,在已知等高线高度的情况下,我们采用了二维插值的方法,并利用MATLAB软件画出了三维地形图,将景区外貌直观的呈现了出来,在计算地表面积时,我们采用了划分网格、近似求值的方法,利用MATLAB所给出的网格平面与水平面的夹角,估算出了地表面积,,对于其误差,我们也进一步给出了分析。第三问中,我们利用第二问中所求出的高度矩阵,用网格中心点代替此网格,给出了任意两点的空间距离,即任意两点的权重,从而建立了一个图论模型,对于该无向图,我们采用Dijkstra算法利用C++,确定出了最佳路线,并运用MATLAB作图直观的将路线做了出来,并估算出最优路线的空间距离长度约为4567m。第四问中,我们将着火点简化为几个最有可能发生火灾并且救援不方便的点,建立了一个目标规划的模型,然后在一定范围内,对消防点进行了假设,利用第三问的C++程序求出了到着火点的最长时间,移动消防点求出了最优消防站的地址。计算得出结论:在给定的坐标系下,最优消防站点的位置位于点(28,25)。【关键词】:网格离散化大型稀疏阵三次样条插值最短路线MATLABC++一、问题的提出某国家级森林公园的地形等高图如图1所示。由于该风景区植被丰富,拥有大量的国家级重点保护动植物,因此旅游管理部门在图1的A点设置了景区消防站,当景区发生火灾时能及时控制和消灭火情。图1说明:该图水平及竖直方向以10m为单位,山高以50m为单位。请你利用所学数学知识回答以下问题:1、由于人为原因,图1所示的等高图出现了局部破损的情况,请利用数学模型修补好该地图;2、在完成第一问的基础上,结合数学模型建立该景区的三维地形图,并估计该景区的地表面积;3、某天图1所示的B点发生了火灾,于是需要从景区消防站派遣消防员去B点灭火,建立模型确定最佳灭火路线。4、如果需要对景区消防站进行重新选址,请建立模型确定合理的消防站地址。二、(一个像素单位为10m),相邻两等高线的高程差为50m;,故假设消防员步行灭火;;,其等高线之间的山体均匀线性变化;;。三、符号说明1.——第个网格,即其编号方式为先行后列,.——表示整个景区的表面积;3.——表示第个网格的面积;4.——表示第个网格与水平面的夹角;5.——表示点到点的权值。四、问题分析第一问中,需要修补该地形图缺失的等高线,主要有两种方法:一是平面插值,二是曲线拟合。由于曲线拟合的原则是使各点距拟合曲线的距离平方和最小,故拟合曲线基本不过已知点。而考虑到题目所给前提为“等高图破损”,故等高线需过已知点,因此我们选用平面插值方法。插值方法有很多种,比如拉格朗日插值,艾米尔特插值,三次样条插值等,我们可以根据各自优劣来选择其插值方式,进而利用MATLAB软件编程计算,绘制出等高线。第二问中,需要画出景区的三维地形图,经过分析,可以利用MATLAB强大的数据处理功能对问题进行求解。对于未知坐标点的高程,可以用MATLAB进行二维插值,从而绘制出地形图。在需要计算地表面积的时候,可以根据假设5的平面假设,利用分块求和的方法,计算出地表面积。第三问中,需要寻找最佳路线,使消防员由消防站A赶到火灾点B的时间最短,即在假设3的速度均匀的条件下,求一条最短时间的路线。可以将最短时间问题简化为最短路径的问题。而对于每两个已经简化的点,由C++编程可以算出其路径,然后根据Dijkstra算法,可以求出其最短路径,即最佳灭火路线。为了将直观将结果展现出来,我们可以做出它的