第3个回答 2010-08-20
奥运会临时超市网点设计模型
(小三黑体,题目直接用竞赛试题题目,不必另起)
摘要 (一级标题,4号黑体,居中)
(论文其他内容小4号宋体字,单倍行距,左侧装订)
本文根据题目附录中提供的问卷调查数据,利用关系数据库查询语言,从不同侧面进行了准确统计,找出了运动会期间观众在出行方式、餐饮方式以及消费额(非餐饮)三方面所反映的规律:大部分(约72%)的观众坐公交和地铁出行;过半数(约52%)的观众选择西餐作为餐饮方式;绝大部分(约88%)的观众消费额在300以下,其中200到300之间人数约占44%。
根据观众在出行方式、餐饮方式以及消费额(非餐饮)三方面所反映的规律,对不同消费档次(非餐饮)的观众进行统计,分别测算出题目(图2)中20个商区的人流量分布:A1:6.83% A2:5.09% A3:5.63% A4:6.19% A5:6.72% A6:11.73% A7:5.04% A8:4.49% A9:3.95% A10:3.40%
B1:2.81% B2:2.26% B3:4.55% B4:3.95% B5:4.49% B6:7.27% C1:1.69% C2:2.60% C3:5.39% C4:5.84%
在解决了问题1、2的基础上,对不同消费档次的观众赋予不同消费档次指数,然后,通过对综合购买力的分析以及对各消费档次观众的消费水平进行全面、综合考查,并以此为依据对问题3建立了线性优化模型,运用数学软件MATLAB编程对模型进行二维搜索,得到了模型最优解,设计出了各商区两种类型迷你超市网MS的分布方案:
商区
网类型 A1 A2 A3 A4 A5 A6 A7 A8 A9 A10
小MS个数 5 4 4 4 5 8 4 3 3 2
大MS个数 5 4 4 5 5 9 4 4 3 3
商区
网类型 B1 B2 B3 B4 B5 B6 C1 C2 C3 C4
小MS个数 2 3 4 3 4 6 2 2 4 3
大MS个数 2 1 3 3 3 5 1 2 4 5
最后,通过综合分析,我们建立的模型能够准确描述各商区消费水平,得出两种不同类型MS个数分布基本均衡,既满足了奥运会期间的购物需求,又考虑了商业赢利。
关键词(一级标题,四号黑体,居中)
人流量;二维搜索;消费档次指数;线性优化模型;综合购买力(3-5个)
(第一页只有摘要和关键词,而且论文从这一页开始编页号,页码居中)
一. 问题的提出(一级标题,四号黑体,居中)
2008年北京奥运会的建设工作已经进入全面设计和实施阶段。奥运会期间,在比赛主场馆的周边地区需要建设由小型商亭构建的临时商业网点,称为迷你超市(Mini Supermarket, 以下记做MS)网,以满足观众、游客、工作人员等在奥运会期间的购物需求,主要经营食品、奥运纪念品、旅游用品、文体用品和小日用品等。在比赛主场馆周边地区设置的这种MS,在地点、大小类型和总量方面有三个基本要求:满足奥运会期间的购物需求、分布基本均衡和商业上赢利。
图1给出了比赛主场馆的规划图。作为真实地图的简化,在图2中仅保留了与本问题有关的地区及相关部分:道路(白色为人行道)、公交车站、地铁站、出租车站、私车停车场、餐饮部门等,其中标有A1-A10、B1-B6、C1-C4的黄色区域是规定的设计MS网点的20个商区。
为了得到人流量的规律,一个可供选择的方法,是在已经建设好的某运动场(图3)通过对预演的运动会的问卷调查,了解观众(购物主体)的出行和用餐的需求方式和购物欲望。假设我们在某运动场举办了三次运动会,并通过对观众的问卷调查采集了相关数据,参照采集的数据,请你按以下步骤对图2的20个商区设计MS网点:
1. 根据附录中给出的问卷调查数据,找出观众在出行、用餐和购物等方面所反映的规律。
2. 假定奥运会期间(指某一天)每位观众平均出行两次,一次为进出场馆,一次为餐饮,并且出行均采取最短路径。依据1的结果,测算图2中20个商区的人流量分布(用百分比表示)。
3. 如果有两种大小不同规模的MS类型供选择,给出图2中20个商区内MS网点的设计方案(即每个商区内不同类型MS的个数),以满足上述三个基本要求。
4. 阐明你的方法的科学性,并说明你的结果是贴近实际的。
(图2,图3请见附录2)。
二. 问题假设(一级标题,四号黑体,居中)
1.奥运会期间(指某一天)每位观众平均出行两次,一次为进出场馆,一次为餐饮,并且出行均采取最短路径。
2.观众在一天内的行程如下:
进场馆——>出场餐饮——>餐饮完回场馆——>出场馆
且进场馆和出场馆路径相同,出场餐饮和餐饮完回场路径相同。
3.出场餐饮与餐饮完回场馆时不考虑出行方式,只按餐饮方式采取最短路径。
4.各场馆内进出口与看台一一对应(即进场时一个进口只能到达唯一确定看台,出场时一个出口对应唯一看台,看台之间不能相互跨越)。
5.每位观众通过出行或餐饮路径上所有商区(包括看台出口所对的商区)。
6.三个场馆人数固定(A区为10万人,B区为6万人,C区为4万人),每个看台人数固定,均为1万人(即商区A1、A2、A3、A4、A5、A6、A7、A8、A9、A10、B1、B2、B3、B4、B5、B6、C1、C2、C3、C4对应的二十个看台每个均为一万人)。
7.观众在奥运期间的出行方式、餐饮方式、消费额档次均不变,且服从问卷调查所得规律。
三. 假设合理性分析及说明(一级标题,四号黑体,居中)
根据最短路径原则,观众从各车站或停车场到场馆往返路径相同;同理,餐饮往返路径也相同。因此只须考虑观众看完比赛从场馆到车站或停车场的路径(下称第一类路径)以及观众出场馆到达餐饮地点的路径(下称第二类路径)即可。即对各商区人流量只须计算这两类路径的人流量,各商区总人流量为观众走这两类路径人流量的2倍。为方便计算,本模型中人流量仅为第一类和第二类路径人流量之和。从图2可以看出,各场馆到餐饮地点或者无车可乘或者相距很近无须乘车,故在观众出场馆餐饮时只根据餐饮方式采取最短路径,忽略出行方式。
四. 符号约定(一级标题,四号黑体,居中)
W: 出行方式为公交(东西);
N: 出行方式为公交(南北);
E: 出行方式为地铁东;
R: 出行方式为地铁西;
P: 出行方式为私车;
T: 出行方式为出租;
C: 餐饮方式为中餐;
F: 餐饮方式为西餐;
B: 餐饮方式为商场;
五. 模型建立与求解(一级标题,四号黑体,居中)
1. 问题1求解
根据附录中给出的问卷调查数据,我们利用数据库编程(Visual Basic +SQL关系数据查询语言)首先统计得出了三次问卷调查中按年龄、出行方式、餐饮方式、消费水平分档的各类人数,如表1所示。
……………………………………………
………………………………………………
………………………………………………..
为了能清楚看出观众在出行、用餐和购物等方面反映的情况,用百分比表示各出行方式、餐饮方式、消费额档次人群的分布情况,如表2所示:(略)
………………………………………………….
…………………………………………………
……………………………………………………
2.问题2求解
商区人流量与平均购物欲望是影响商区选址的主要因素。各商区人流量与观众出行方式、餐饮方式有关。商区人流量的消费档次水平分布,体现了该商区人流的平均购物欲望。因此,以消费档次水平为划分标准,分别按出行方式及餐饮方式对人群进行统计,不同消费档次水平人数及百分比表示如表3所示:
……………………………………………..
……………………………………………...
……………………………………………
3.问题3求解
…………………………………………..
………………………………………….
商区Z的综合购买力(百万元)H =商区Z各个消费档次购买力之和。
各个消费档次购买力为:该消费档次人流量╳消费档次指数
根据以上标准可以建立以总出售能力最小作为目标函数的模型:
Min f=m1╳( + + )+m2╳( + + )
约束条件为:
╳m1+ ╳m2>= (i=1,2……10)
╳m1+ ╳m2>= (j=1,2……6)
╳m1+ ╳m2>= (k=1,2,3,4)
, , , , , >=1且为整数
m1<m2
满足以上条件的 , , , , , (i=1,2……10;j=1,2……6 ;i=1,2……4)的值即为每个商区两类不同规模迷你超市的个数。利用MATLAB线性优化并同时对其进行编程,对m1,m2的值进行二维搜索,并通过逐步缩短m1,m2的搜索步长以提高精确度,确定最优解。
步长为0.2时的最优搜索结果如下(程序见附录1):
………………………………..
………………………………….
………………………………….
模型最优解分析:
(1)由上图所示的设计方案可以看出,在每个商区内大小规模的MS网分布比较均衡。计算所有场馆中大型MS网个数:75,小型MS网个数:75。即从整体范围大小规模的MS网分布比较均衡,即“分布基本均衡”得到满足。
(2)模型中把“满足购物需求”作为建模的条件,把“商业上赢利”作为目标函数。故求得的最优解必满足以上两要求。
综合以上分析,由此模型所得的20个商区内MS网的设计方案为符合要求的可行方案。
六.模型合理性分析及问题4解答(一级标题,四号黑体,居中)
模型利用数据库编程(Visual Basic +SQL关系数据查询语言)统计三次调查问卷,分别得出三次问卷调查中观众在出行、用餐和购物等方面所反映的规律,三次所得规律大致相同,故采取三次调查的平均值以更符合一般实际。
模型以预演运动会调查问卷所得统计数据作为建模所用数据,即把实际情况所反映的规律运用到所建模型中,并根据实际情况作出合理假设,使建模条件和实际情况更加接近。
在对问题2进行求解时,由于问题中没有给出路径距离的有关数据,模型根据图示合理制定相应的最短路径的行程路线,提供了一个最短行程的可行性方案。
模型用销售能力来衡量商区网点的规模,使设计问题合理地转化为最优化问题,以使所得结果在满足题目基本要求的前提下达到最优。从MS设计方案可以看出,各商区内大小两种不同类型的MS分布基本均衡,与各商区的综合购买力相匹配,在解模过程中运用数学软件MATLAB处理了线性最优化问题,大大提高解题速率,并使模型精确度更高。
参考文献(一级标题,四号黑体,居中)
[1] 王沫然 MATLAB 6.0与科学计算 北京 电子工业出版社 2001年9月
[2] 钱颂迪等 运筹学 北京 清华大学出版社 1990年1月
[3] 孙洪祥等 概率论与随机过程 北京 北京邮电大学出版社 2001年2月
[4] 姜启源 数学模型 北京 高等教育出版社 1993年8月
[5] 王能超 计算方法简明教程 北京 高等教育出版社 2004年1月
附录1:
利用MATLAB进行二维搜索的线性优化程序
min_value=100000;
t=1:1:40;
j=1:1:40;
m1=1;
m2=1;
while m1>=1 && m1<=4
m2=m1+0.2;
while m2<=7
s1=0;
vlb=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1];
vub=[];
a=[-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0,0;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2,0;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-m2];
%b=[-20.1020,-13.4069,-15.3143,-16.7830,-18.4440,-35.2131,-16.7830,-15.0859,-13.4069,-11.7279,-10.0663,-8.3783,-13.4027,-10.0663,-11.7363,-21.7978,-5.9456,-7.7796,-14.9323,-17.3817];
b=[-20.1027,-15.0825,-16.7618,-18.4591,-20.1203,-35.2161,-15.1094,-13.4121,-11.7328,-10.0535,-8.3921,-6.7038,-13.4033,-11.7418,-13.4121,-21.7996,-5.0290,-7.7802,-15.8503,-17.3827];
c=[m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m1,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2,m2];
[x,lam]=lp(c,a,b,vlb,vub)
for i=1:20
s1=s1+x(i)*m1+x(20+i)*m2;
end
if min_value>s1
min_value=s1;
t=x;
p=m1;
q=m2;
end
m2=m2+0.2;
end
m1=m1+0.2;
end
plot(j,x);
附录2:
图二
图三