消防站选址 建模论文Word文档下载推荐.docx
《消防站选址 建模论文Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《消防站选址 建模论文Word文档下载推荐.docx（12页珍藏版）》请在冰豆网上搜索。
模型一即重心法求取模型给出了一个最佳选址点，模型二即中心法求取模型给出了两个选址点。
关键词：
重心法迭代法中心选址样本位置点C语言编程
一、
问题重述
随着滨州市城市规划的实施和城市建设的发展,特别是西区、东区开发和旧城改造的加快，滨州市原有的消防站点已经不能满足城市日益发展的要求。
2009年8月8日，滨州格林豪泰大酒店（原帝都宾馆），建筑面积11000平方米，8月8日上午9时38分左右，正在装修中二楼突然起火。
二楼三楼正在装修的工人发现着火后虽然也进行了扑救，但火势过猛扑救无望后工人全部撤出现场。
10点18分滨州市消防支队接到报警，出动25辆消防车、120名消防官兵进行了扑救，大火于中午12时30分扑灭，二楼1500平方米被大火烧毁，造成了重大的经济损失。
请根据滨州市区的有关情况，收集诸如各居民小区的位置，规模，各大工厂的位置和规模等相关数据，通过数学建模的方法，再选出1-2个地址供消防站选择。
最后，根据你们建模分析的结果，给有关部门如消防支队写一份报告，提出选址具体建议。
二、模型假设
1、从市区内的各小区、工厂区、学校、商业区根据规模大小、发生火灾可能性大小抽取各样本位置点，并在坐标轴上标出；
2、假设各位置点发生火灾的可能性根据生活常识可以得出；
3、假设消防车在到达火灾地点的途中没有障碍；
4、假设道路曲折程度在可接受范围之内。
三、符号说明
模型一：
di：
样本位置点i到消防站的距离；
vi：
样本位置点i发生火灾的可能性；
xi：
样本位置点i的横坐标；
yi：
样本位置点i的纵坐标；
minz：
各样本位置点i到消防站最短距离之和；
x0：
消防站的横坐标；
y0：
消防站的纵坐标。
模型二;
Vi:
选取的样本位置点i；
D（Vi,Vj）表示各各点的距离;
l（Vi）：
表示行中最大值；
L（Vi）：
表示各行的和
四、建立模型
问题的模型分析及结果
重心法求取模型
将市区比例缩小建立坐标轴，抽取的样本位置点在上面标出，各样本位置点坐标见附录。
由公式：
x0=（∑xivi）/（∑vi）
y0=（∑yivi）/（∑vi）
x0、y0消防站站的地址坐标，即这些样本位置点的重心位置，到它们距离之和最短，有
minz=∑vidii遍历所有样本位置点
…（0）
在坐标平面中设待定消防站的位置为（x0，y0）
x0=
...
（1）
y0=
…
（2）
其中di可以坐标间公式
di=
…（3）
该方法求解过程：
1、确定各目的地点的坐标；
2、不考虑距离因素，用重心公式估算初始选址点：
…（4）
…（5）
（4）、（5）式中含有未知数，不能一次求解，故使用迭代法。
具体步骤为：
1．由
（1）式，先将区域内所有应急点的几何重心点坐标（x0，y0）作为初始地点；
2．将（x0，y0）代入式
（2）中计算出与（x0，y0）相对应的总距离z0；
3．将（x0，y0）代入式（3）中计算出改善的服务点（x1，y1）；
4．利用式
（1）计算出与（x1，y1）相对应的总距离z1；
5．比较z0与z1，若z1≥z0，则（x0，y0）为最优解，计算完毕，若z1<
z0，则返回第3步计算，将（x*，y*）代人式（3）中计算出另一改善的服务点（x2，y2），反复迭代，直至zk+1≥zk，求出最优解为止。
经过C语言编程对这一模型的分析，得出
，
，各个样本位置点与这一点最短距离之和
，得出消防站大致位置在新立河东路滨州学院附近，如下图黑点所示。
模型二：
中心选址法求取模型
在市区东西两侧各选取一个点作为消防站位置。
1、确定各抽象点，如下图所示
注意：
Vi后面的数字或汉字仅代表在那附近选取的点。
其比例尺为1厘米＝500米
2、分析一个有对称中心的图形其中心，自然就取其对称中心。
若一个不规则图形怎样定义其中心呢？
在此，我们用圆去覆盖该图形，其中半径最小的圆的圆心既定义为该图形的中心.
对滨州选取的参考区域如下图一：
（一）在市区西侧选址
V1：
左侧起点V2：
6.32米V3：
3.2公里V4：
4.1公里V5：
4.5公里
V6：
中间起点V7：
5.2公里V8：
总长6.3公里处V9：
总长2.0公里处
V10：
总长3.4公里处
例如：
V3：
3.2公里就是3.2公里处的那个点是V3
下图反应对上面地图的抽象，表示从一点出发可以到达各个点
V1V2V3V4V5
V6V7V8V9V10
表格符号：
l（Vi）表示行中最大值；
L（Vi）表示各行的和。
D（Vi,Vj）
V1
V2
V3
V4
V5
V6
V7
V8
V9
V10
l（Vi）
L（Vi）
0.6
2.7
2.8
3.4
3.2
4.1
6.8
7.4
2.6
2.9
4.2
6.9
7.3
34
0.9
1.3
1.4
2.0
2.1
5.0
4.9
0.44
0.76
1.2
1.5
4.4
4.7
0.62
0.7
4.0
4.6
0.86
0.8
3.6
1.1
3.5
4.5
3.3
l（Vi）=min{（Vi,V1）,d（Vi,V2）……,d（Vi,V10）}
实际上l（Vi）为以Vi为圆心，l（Vi）为半径覆盖所有顶点的圆的半径。
根据定义，取min{l（Vi）∣i=1……10}=l（V5）=4.0为半径，V5为圆心的圆是覆盖所有定点的最小的圆，所以应选为定点v5为最优地址。
L（Vi）=d（Vi,V1）+d（Vi,V2）+……+d（Vi,V10）为以出发到其它各点最短路距离之和。
min{L（v1）,L（V2）,……,L（V10）}=L（V5）=18.66以为最优选址，到各点距离之和最小为18.66。
西侧建在阳光贵苑附近，即上图黑点所示。
（二）在市区东侧选址
右侧起点V2：
2.6公里V3：
3.4公里V4：
4.3公里V5：
5.2公里
6.1公里V7：
1.2公里V8：
8.7公里处V9：
总长9.8公里处
总长6.3公里处V11:
中间起点
3.2公里就是3.2公里处的那个点是V3。
表中符号：
表示各行的和
V11
l（Vi
2.5
2.3
3.0
3.1
5.6
6.3
1.6
1.7
2.2
3.9
1.0
0.94
0.88
1.8
2.4
4.3
0.98
3.7
3.8
26
14.2
15.0
0.79
然后计处：
l（Vi）=min{（Vi,V1）,d（Vi,V2）……,d（Vi,V11）}
根据定义，取min{l（Vi）∣i=1……11}=l（V6）=３.７为半径，V6为圆心的圆是覆盖所有定点的最小的圆，所以应选为定点V6为最优地址。
L（Vi）=d（Vi,V1）+d（Vi,V2）+……+d（Vi,V11）为以出发到其它各点最短路距离之和。
min{L（v1）,L（V2）,……,L（V11）}=L（V5）=19.03以为最优选址，到各点距离之和最小为19.03。
L（V6）=20.26与L（V5）=19.03进行比较，可得出结论选择V6处较好。
由此得出东侧建在滨城区劳动局附近，即上图黄点所示。
五、模型的评价与推广
模型一采用重心法求得了距市区中各小区、工厂、学校、政府机关及商业区距离和最小的且能有效出警的消防站地址，用此法也可选出两个最佳地址，即将东西两侧的样本点划分，中间有所相互覆盖，将各点带入C程序即可求得。
模型二运用中心选址法求得到各样本位置点距离之和最小的两个最佳地址，在这两位置建立消防站可以有效的出警，保障广大市民的人身财产安全。
两个模型方法简单明晰，可以对消防站的选址进行估计。
六、参考文献
[1]姜启源，《数学模型》，北京：
高等教育出版社，2005.
[2]重心法在应急系统选址中的应用，
NewsID=2009623162411，2009-6-23
[3]谭浩强，《C语言程序设计教程》（第三版），北京：
高等教育出版社，2007
[4]城市消防站的布局与选址，
[5]冯杰黄力伟等，《数学建模原理与案例》，北京：
科学出版社，2007
附录
附录一：
模型一各样本位置点坐标及火灾发生系数
位置点
坐标
发生系数
1
（0.46，2.7）
0.2
8
（5.6，2.9）
2
（0.62，205）
9
（7.5，2.2）
0.1
3
（0.61，3.5）
10
（8.9，2.4）
4
（3.0，3.4）
11
（8.0，3.9）
5
（3.9，4.2）
0.08
12
（5.8，6.0）
0.12
6
（4.2，3.2）
13
（8.8，4.9）
7
（4.1，1.9）
14
（10.1，2.4）
0.15
附录二：
C语言程序
#include<
stdio.h>
math.h>
structxuanzhi
{
doublex,y,v;
}xz[14]={{0.46,207,0.2},{0.62,2.5,0.2},
{0.61,3.5,0.1},{3.0,3.4,0.1},{3.9,4.2,0.08},
{4.2,3.2,0.1},{4.1,1.9,0.1},{5.6,2.9,0.12},
{7.5,2.2,0.1},{8.9,2.4,0.2},{8.0,3.9,0.1},
{5.8,6.0,0.12},{8.8,4.9,0.08},{10.1,2.4,0.15}}
;
main（）
structxuanzhi;
inti;
doubled[14];
doublex0=0.0,y0=0.0,min=0.0,z=0.0,t=0.0,m=0.0,w=0.0,
x1=0.0,y1=0.0,m1=0.0,m2=0.0,w2=0.0,t1=0.0,m3=0.0,w3=0.0;
for（i=0;
i<
14;
i++）
{
m+=xz[i].v*xz[i].x;
w+=xz[i].v;
m1+=xz[i].v*xz[i].y;
}
x0=m/w;
y0=m1/w;
d[i]=sqrt（（xz[i].x-x0）*（xz[i].x-x0）+（xz[i].y-y0）*（xz[i].y-y0））;
m3+=xz[i].v*xz[i].x/d[i];
t1+=xz[i].v*xz[i].y/d[i];
w3+=xz[i].v/d[i];
x1=m3/w3;
y1=t1/w3;
do
{
x0=x1;
y0=y1;
m2+=xz[i].v*xz[i].x/d[i];
t+=xz[i].v*xz[i].y/d[i];
w2+=xz[i].v/d[i];
x1=m2/w2;
y1=t/w2;
m2=0.0;
t=0.0;
w2=0.0;
while（fabs（x0-x1）-0.00000001>
0&
&
fabs（y0-y1）-0.00000001>
0）;
printf（"
x0=%12lf,y0=%12lf\n"
x0,y0）;
min=xz[i].v*d[i];
z+=min;
minz=%12lf\n"
z）;
}
升级会员 