《中国数学建模竞赛试题-A-奥运会临时超市网点设计优秀论.doc》由会员分享,可在线阅读,更多相关《中国数学建模竞赛试题-A-奥运会临时超市网点设计优秀论.doc(16页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、欢迎光临中国数学建模网奥运会临时超市网点设计的数学模型摘要 本文针对题目要求,分统计问卷,测算人流量分布,网点设计三个步骤逐步进行,在满足三个基本要求的基础上,构建出了以两种超市个数为变量的整数规划模型,基本上解决了北京奥运会临时超市网点设计问题。在步骤一里,我们通过 SPSS分析软件对调查问卷进行了统计分析,得出观众在出行和餐饮方式上的偏好规律如下:公交车 出租私车地铁中餐西餐商场34%19%9%38%22.5%52.5%25%步骤二中,我们通过 floyd算法,并结合步骤一中得到的数据,我们最终测算出 20个商区的人流量分布:C1C2C3C4B1B2B3B4B5B62.32% 2.96%
2、2.33% 5.08% 4.44% 4.00% 6.40% 4.2% 4.57% 9.72%A1A2A3A4A5A6A7A8A9A107.09% 3.84% 4.03% 4.64% 5.47% 11.1% 5.26% 4.43% 4.03% 4.06%步骤三中,我们结合最优产出理论证明了同一商区内大小超市的最优数量比,基于步骤二的结果建立模型,并通过科学假定参数,通过 Lingo软件得到了一组网点设计方案,以供决策者参考:C125C225C325C438B126B226B339B426B527B6414大超市小超市A1310A225A326A427A538A6516A738A826A926A1
3、026大超市小超市最后针对模型的特点,我们阐释了其方法的科学性,结果的现实性和一些改进之处。关键词:floyd算法 最优产出理论 整数规划感谢您对网站建设的支持和参与1欢迎光临中国数学建模网一、 问题重述 在即将到来的 2008 年北京奥运会的主赛场周边地区需要建立 MS(迷你超市)网以满足游客的消费需求。MS在地点、大小类型和总量方面有三个基本要求:满足奥运会期间的购物需求、分布基本均衡和商业上赢利。题中给出了真实地图的简化图,以及在附录中给出的三次调查的数据,现要求做到:1、根据问卷得出观众在出行、用餐和购物等方面所反映的规律。2、测算图2中20个商区的人流量分布(用百分比表示)。3、仅考
4、虑两种规模的MS,给出图2中20个商区内MS网点的设计方案(即每个商区内不同类型MS的个数),以满足上述三个基本要求。4、说明解决方法的科学性,并说明结果是贴近实际的。二、 问题的分析 问题的总括: 通过统计调查问卷,得出运动会观众出行和用餐的需求偏好及购物欲望,并以此为根据测算出北京奥运会体育馆周边各个商区的人流量分布,从而更进一步结合实际,设计出各个商区内MS网点的最优分布。问题的宗旨:1、网点能满足观众的购物欲望(非餐饮方面)。2、两种类型的MS网点分布基本均匀。3、使得期间网点商业净利润为正,并且尽可能大。问题的关键与难点:1、找准问卷反映的规律。2、测量最短路径,从观众的角度选择路径
5、。3、把两种MS规模及商区的面积按实际量化建模,定量设计网点。综上,问题的解决过程要分找规律、求人流量比例,设计网点三个步骤。(即题目顺序)三、 模型的基本假设 1、对步骤二的假设:1)奥运游客对出行方式和餐饮方式的偏好与调查所得的规律一致。2)奥运会期间(指某一天)每位观众平均出行两次,一次为进出场馆,一次为餐饮。3)观众根据手中的地图,选择最短路径,且按原路返回。4)国家体育场(鸟巢)容量为10万人,国家体育馆容量为6万人,国家游泳中心(水立方)容量为4万人。三个场馆的每个看台容量均为感谢您对网站建设的支持和参与2欢迎光临中国数学建模网1万人,出口对准一个商区,各商区面积相同。5)观众从各
6、个交通工具的站点下车后,均采取步行的方式到达体育场馆。2、对步骤三的假设:1)网点的商圈是以网点为中心的圆。2)只考虑两种大小规模的商圈,且大小超市吸收顾客的能力(可视为收益)与其商圈面积的大小成正比。3)不考虑商区内超市网点的竞争,任何网点的商圈彼此相离。4)在简化模型中暂且认为,每天一个商区内的需求量均保持一定。四、 符号说明 二 十 个 商 区 的 集 合 , 即Ai,B j ,Ck,i = 1.4, j = 1.6,k = 1.10中的元素,即 Nn商区中大规模超市的个数商区中小规模超市的个数一个大规模超市的总成本(包括构建成本,经营成本等)一个小规模超市的总成本CbCsRbRsD商区
7、中所有大超市的营业收入之和商区中所有小超市的营业收入之和商区的总需求商区一天的人流量Q 五、 进一步分析与模型的建立 步骤一:统计调查问卷调查问卷是获得客户需求以更好的制定决策的主要方式,对此调查问卷整感谢您对网站建设的支持和参与3欢迎光临中国数学建模网理和统计后所得出的规律将对后来问题的决策有着决定性的影响。对于这三次共 10600份调查问卷,我们通过运用 Excel将部分原数据转化成0-1矩阵并将其转录到 SPSS分析软件上,分以下几项进行统计研究。1、偏好:分别对每一次调查中观众对于出行方式和用餐方式的偏好进行百分比求值:(%)表 1公交车 公交车南北 东西出租 私车地铁东 地铁西 中餐
8、 西餐 商场第一次第二次第三次17.516.81617.117.417.219.4 8.818.6 9.218.8 9.118.418.919.418.819.119.422.422.622.452.552.352.825.125.124.8观察到三次统计结果十分接近,可以认为三次调查问卷的性质完全相同,故综合考虑,结果如下:(%)表 2公交车 出租私车地铁中餐西餐商场341993822.552.5252、性别与消费额:女平均消费额(按问卷上的等级):2.739男平均消费额:2.31看出女子消费者的消费量要高于男子消费者的消费量。3、年龄与消费额:20岁以下:2030:3050:1.992.7
9、92.501.5850岁以上:4、人均消费额:将每档消费额取中点值,结合人数比例求得人均消费额Re=201.7元其中消费额的分布直方图如下:消费额(非餐饮)500040003000200010000123456消费额(非餐饮)感谢您对网站建设的支持和参与4欢迎光临中国数学建模网5、其他:通过 SPSS软件的相关性分析,可认为出行方式与消费额,用餐方式与消费额均没有明显规律和关系。步骤二:测算人流量百分比根据图二,我们描绘出进一步简化的图 1:问题可以考虑为找出每个出口的人根据自己的偏好在最短路径上经过的商区汇总。偏好的分布我们已由步骤一得到,因此问题的关键是找出最短路程。用画图工具测出相关的路
10、程,简化整理成赋权图:(我们先舍弃了与路口不相接的商区)经计算机程序的计算(见附录(1),我们找出了从图中所选的六个商区出发,选择各出行用餐方式的最短路径,进而我们补充汇编了从任一商区到任一种出行和用餐场所途经的商区汇总矩阵(见附录(2)。但从实际角度出发,观众感谢您对网站建设的支持和参与5欢迎光临中国数学建模网在某几种情况下(如从 B1出发去公共汽车车站)很难判断两种路径细微的差别,因此我们在这里把不能通过观察判断的情况找出,假定观众随机的选择两种路径,进行修正。程序和修正后的商区汇总矩阵见附录。经过统计观众进行各种出行和用餐所经过的各商区的数目,我们得到了一个 20 7的矩阵 K1(见附录
11、(3),表示选择每种出行或用餐方式的观众所经过各商区的人流量百分比。按公交车,地铁,出租,私车,中餐,西餐,商场用餐顺序,现定义偏好向量 H和人流向量 J,使得 J=QH,表示了选择各偏好的往返人数。其中,Q=400000,表示一天观众出行或用餐的往返人流量H=0.34 0.19 0.9 0.38 0.225 0.525 0.25,表示各出行和用餐方式在同类中的比例。再令 K2= K1J,则根据 K2各行向量之和可求得一天内任一商区的人流量Q,从而得到商区的人流量百分比。综上,步骤二要建线性模型如下:J=QHK2= K1J具体数值计算留在模型求解中计算。(1)(2)步骤三:最优网点设计这是问题
12、的难点,也是我们问题研究的最终目的。一、难点处理:由于无法掌握四年后的具体数据以进行分析求解,我们只能通过两个途径逐步解决:1)模型的简化。研究对象是一个非常复杂的动态系统,具有不确定性。可以考虑先通过进一步假设简化建模。首先,在某一商区内一天的人流量Q和需求量 D是固定的,其中的 D可以由调查问卷测算的人均消费额 Re和Q确定,有:D =Re Q(3)另外,大小商圈可以看成是两个面积比值为 p(p1)的圆形区域,依据基本假设,大小超市吸收顾客的能力之比也为 p,则有:Rbp N=(4)Rsn再有,顾客的总需求量 D为超市提供了利润的来源,这里认为,在满足需感谢您对网站建设的支持和参与6欢迎光
13、临中国数学建模网求的情况下有:D = Rb + Rs = Re Q(5)2)参数估值。在复杂性之外,有很多必要的参数也是未知的,这就给模型建立后的求解带来了很多不便。可考虑自估参数或查阅经验数据,对模型的可行性进行检验,或者保留参数,待日后确定。二、问题的突破口先分析一个商区内的网点分布,由于大小超市均涉及到数量,成本和收益的权衡问题,因此我们把大小规模的超市看成两种产品,且有:Cb =pCs(6)我们要解决的问题是:在市场(总需求)一定的情况下,两种超市的最优产出是多少。由(3)(4)不难得出:Rb =p Np N D+ n(7)(8)大规模超市的边际收益MRb = RbND n p(p N
14、 + n )=2同理有:小规模超市的边际收益RsD N p(p N + n )MRs = n=(9)2MRb = MCb = Cb根据经济学中的最优产出理论,当MRs = MCs = Cs成立时,利润是最大的,再代入(7)(8)两式,等式两边分别相除得:n= Cbs= p(10)NC(对最优产出理论的简略证明见附录(4)这样,我们得到了两种超市的最优数量配比率,即等于其成本之比的倒数。三、确定模型我们得出的最佳个数安排,仅仅是对于该商区内的超市赢利最大的必要条件,并非充分条件,因此有必要深入的优化组合。在保证商业赢利的同时,使得商区内大小超市的布局满足另外的一些基本要求:感谢您对网站建设的支持
15、和参与7欢迎光临中国数学建模网a)满足购物需求。则有:S D代入(3)即:A(n + p N ) D = ReQ(11)其中,A为一个小商区的供应量。b)保证赢利。Cs n + Cb N R = ReQc)比例接近 p。(12)考虑到 N,n是整数,因此我们不能强要求 n= p,我们转而将第三个N约束条件放宽,即允许其在一定范围内波动,假设其波动的大小为(0 1),这时大小超市个数的最佳比例约束为:n p N(13)在保证这两个前提的情况下,我们将商区内大小超市的净利润之和作为目标函数,把(11)(12)(13)作为约束条件综合考虑得到整数线性规划模型:Max ReQ Cs n Cb NCs
16、n + Cb N ReQA(n + p N ) ReQ目标函数赢利约束需求约束n p N最佳比例约束六、模型求解 1、步骤二的求解经对最短路线的统计,可得:K2=K1J=K1HQ=感谢您对网站建设的支持和参与8欢迎光临中国数学建模网2615.47846.22615.42615.42615.43923.17846.22615.43923.17846.22615.43923.16538.59153.81176926154117699153.86538.53923.13931.032620.693931.0310482.86551.723931.032620.693931.036551.721572
17、4.12620.693931.036551.729172.4111793.126206.911793.19172.416551.723931.031965.55241.41965.51310.31965.53275.97862.13275.91965.51310.3131035896.64586.23275.91965.51310.31965.53275.94586.25896.6692.31385692.31385692.310382077692.31038207769233115242317311038692.310381731242331152046.41364.32046.45457.
18、11364.32728.68185.72728.61364.31364.3136436139.347751859.321239.541859.324958.185577.953718.636197.723718.634958.1812395.41239.541859.323098.864338.415577.9512395.45577.954338.413098.861859.324963.773309.184963.7713236.715718.612409.414891.315718.615718.634746.414891.34963.773309.184963.778272.95182
19、00.56618.363309.183309.186618.363410.72046.41364.32046.43410.747756139.3对于 K2通过横向求和得出一天内各商区的人流量,转化成百分比得到如下结果:C1C2C3C4B1B2B3B4B5B62.32% 2.96% 2.33% 5.08% 4.44% 4.00% 6.40% 4.2% 4.57% 9.72%A1A2A3A4A5A6A7A8A9A107.09% 3.84% 4.03% 4.64% 5.47% 11.1% 5.26% 4.43% 4.03% 4.06%2、步骤三的求解由于Cs,Cb,p,A,均为未知参数,这里我们采取
20、合理的自估,来检验此模型的可行性和有效性。当然,前四个自设数据不能代表将来实际的数据,我们能做的只是尽量贴切实际的估计,但实际分布结果还要以实际的Cs,Cb,p,A为标准在模型中求得。通过查阅相关资料,我们获悉小商亭营业面积为 10平方米左右,国家体感谢您对网站建设的支持和参与9欢迎光临中国数学建模网育场的面积为 7万平方米等信息,随即我们按比例估算了商区的面积,为 43200平方米左右,进一步考虑,我们假设Cs =30000,Cb =90000,p=3,A=15000, =0.5运用 lingo求解得出 20个商区的分布情况为:C125C225C325C438B126B226B339B426
21、B527B6414大超市小超市A1310A225A326A427A538A6516A738A826A926A1026大超市小超市(程序见附录(5)七、模型的检验和修整模型二建立在 D = R = ReQ上,即市场均衡,但我们了解到,现实中由于人心主观的一些不确定性的因素的存在,市场往往达不到完全均衡的状态。我们注意到,当 A(n + p N ) ReQ时,并不一定意味着满足了观众的购买需求。因为当货物所剩无几时,人们的购买欲望会大幅降低,整体来看,实际消费额 L时,已基本满足顾客,再增加超市便会使商业利润降低。我们假定,商区的人流量Q越大,潜在购买力就越大,其对应的值就越大,且两者成正比例关系
22、,我们可以定义饱和系数t = Q Re,描述一个商区的S饱和程度,其中,S指该商区内部的所有商圈面积之和,即 S=SbN + Ssn经上述分析,我们的模型二可以进行适当改进如下:改进的模型Max kReQ Cs n Cb NCs n + Cb N ReQ k目标函数赢利约束ReQA(n + p N ) 需求约束tSsn p N最佳比例约束其中 Ss表示小超市商圈的面积这样,我们的未知参数由原来的五个变为现在的七个,新增了两个参数。其中,k值应由决策者根据市场中人们的消费信心来内部拟定,在保证基本满足需求的前提下,k值要大于 0. 95。感谢您对网站建设的支持和参与11欢迎光临中国数学建模网而
23、t值则是这所有二十个商区共有的数量属性,决策者可以考虑针对其中一个商区进行小规模的商业试验,估计出一个比较可靠的 t值,从而做出更加合理的决策。八、方法和结果的评价一、模型的操作性该模型突出要点,简单有效,只要测出参数值,编程是很简单的,易为很多人所接受;另外,我们得出的超市最佳数量比也很容易由超市的成本算出。二、方法的科学性首先,观众选择路径,构成通过某商区的人流量本是一个不确定的事件,但不确定的事件大量累积后,便会呈现出稳定性的一面,而此问题中的人数达到了上万甚至上十万,可以说提供了稳定性形成的条件,因此此模型假设固定商区的人流量固定是较为合理的。由于大小超市的资金投入,产出,及收益风险均
24、是市场中商品的一般性质,像公司决策生产两种大小不同的产品一样,我们要在配比及总产出上进行权衡,因此,此模型将两种超市看作是两种产品,并赋予它们商品的特征,从而进行边际收益的分析是有根据的。我们在修整的模型里考虑了消费者心理的因素,从更本质的角度考虑表象上的问题,可以使模型的有效性增强。此外,各个商区的人流量基本不同,但各商区的性质是相同的,因此,我们对单一的商区进行研究得到的模型是适用于整个整体的。三、结果的现实性在步骤二中,我们对计算机求最短路径和观众选择最短路径在个别情况下加以区别,认为当两种路径较为接近时,顾客随机的选择一条路径。这样划分是很贴近现实的,因此得到的结果也应该较为准确。在步
25、骤三中,我们在自估参数值时考虑的较为谨慎。通过查阅相关的资料并经过一定的测算,如小商亭营业面积为 10平方米左右,国家体育场的面积为7万平方米等,我们按比例估算了商区的面积,为 43200平方米左右,从而经过模型计算得出结果。从结果的来源看,这是经过我们贴近现实的假设,科学合理的分析得到的;从结果本身意义看,在一个人流量最大的商区内建小迷你超市,建 21个应不足为过,而且,不同商区超市数量之间的关系也和我们的预期较为一致。九、参考文献 1高鸿业,西方经济学(微观部分),北京:中国人民大学出版社,2003.92洪毅等,经济数学模型(第二版),广州:华南理工大学出版社,2002.123沈俊等,数据
26、结构C+实现,北京:科学出版社,2002.7感谢您对网站建设的支持和参与12欢迎光临中国数学建模网八、附录 1、求观众出行最短路径的 Matlab的 m文件如下:functionD,R=floyd(a)n=size(a,1);D=afor i=1:nfor j=1:nR(i,j)=j;endendRfor k=1:nfor i=1:nfor j=1:nif D(i,k)+D(k,j)D(i,j)D(i,j)=D(i,k)+D(k,j);R(i,j)=R(i,k);endendendkDRend2、公交车C1C2地铁出租私车C1C2C3C4B1B2B3B4B5B6A1A2C1C4C1C2C2C1
27、CC2C2CC4C3C4C2C3C2C3C2C2CC4B1B2B3B2B3B3C3CC4C4C4B1 B2B3/B6B2B3B1B6B1 B2B3/B6B2B3B3B2B1B6B3BB6B4B5B6B5B6B3B4 B3/B5B6B5B6B4B3B5B4B3B3BB6A1B4 B3/B5B6B5B6B6B6B6A1AA6A2A3A4A5A6A1AA6A2A3A4A5A6A1A2A1A2A1感谢您对网站建设的支持和参与13欢迎光临中国数学建模网A3A4A5A6A7A8A9A10A3A4A5A6A4A5A6A5A6A3A4A5A6A4A5A6A5A6A3A2A1A3A2A1A4A3A2A1A5A4
28、A3A2A1A6AA1A4A3A2A1A5A4A3A2A1A6AA1A6A6A7A6A7A6A7A8A9A10A1A8A9A10A1A9A10A1A7A8A9A10A1A8A9A10A1A9A10A1A8A7A6A9A8A7A6A10A9A8A7A6A8A7A6A9A8A7A6A10A9A8A7A6A10A1A10A1中餐西餐商场C1C1C4C1C4B3BB6C1C4C2C3C4B1B2B3B4B5B6A1A2A3A4A5A6A7A8A9A10C2CC4C3C4C2CC4B3BB6C3C4B3BB6C4B3BB6B1B6C2CC4C3C4C4C4B1B2B3B2B3B1B6B2B1B6B3BB
29、6B2B1B6B3BB6B4B5B6B5B6B3B4B3B4B5B6B5B6B5B4B3B3BB6A1B6B6A1AA6A1B3BB6A2A1A2A3A4A5A6A3A4A5A6A4A5A6A5A6A2A1B3BB6A3 A2A1B3BB6/A4A5A6A4A5A6A3A2A1A4A3A2A1A5A4A3A2A1A6AA1A7A8A9A10A1A8A9A10A1A9A10A1A10A1A5A6A6A6A7A6A7A6A8A7A6A9A8A7A6A10A9A8A7A6A8A7A6A9A10A1B3BB6A10A1B3BB63、K1=0.01923 0.02586207 0.0259 0.0192
30、3 0.027 0.0206897 0.0236220470.05769 0.01724138 0.069 0.03846 0.018 0.0137931 0.0157480310.01923 0.02586207 0.0259 0.01923 0.027 0.0206897 0.0236220470.01923 0.06896552 0.0172 0.03846 0.071 0.0551724 0.0629921260.01923 0.04310345 0.0259 0.01923 0.018 0.062069 0.074803150.02885 0.02586207 0.0431 0.02
31、885 0.036 0.0413793 0.059055118感谢您对网站建设的支持和参与14欢迎光临中国数学建模网0.05769 0.01724138 0.1034 0.05769 0.107 0.0689655 0.0708661420.01923 0.02586207 0.0431 0.01923 0.036 0.04137930.02885 0.04310345 0.0259 0.02885 0.018 0.05517240.074803150.074803150.05769 0.10344828 0.0172 0.05769 0.018 0.137931 0.1653543310.0
32、1923 0.01724138 0.1724 0.19231 0.179 0.0137931 0.0708661420.02885 0.02586207 0.0776 0.08654 0.08 0.0206897 0.0236220470.04808 0.04310345 0.0603 0.06731 0.063 0.0344828 0.0157480310.06731 0.06034483 0.0431 0.04808 0.045 0.0482759 0.0236220470.08654 0.07758621 0.0259 0.02885 0.027 0.062069 0.039370079
33、0.19231 0.17241379 0.0172 0.01923 0.018 0.137931 0.0866141730.08654 0.07758621 0.0259 0.02885 0.027 0.062069 0.0314960630.06731 0.06034483 0.0431 0.04808 0.045 0.0482759 0.0157480310.04808 0.04310345 0.0603 0.06731 0.063 0.0344828 0.0157480310.02885 0.02586207 0.0776 0.08654 0.08 0.0206897 0.0314960
34、634、证明:设总成本函数为 C = C(q1,q2),销售这两种产品的总收益函数为R = R(q1,q2),则厂商的利润函数为:L = R C = R(q1,q2) C(q1,q2)由极值存在的必要条件: L = R C = 0q q1 q11L= R C = 0q2q2q2现记MR = Rq,MR2 = 为边际收益,并记MC1 = qC ,MC2 = 为边际成RC1q2q211本。于是就有结论:获得最大利润的充分必要条件为:MR1 = MC12MR = MC25、Lingo程序如下(只取其中一个,其他类同)model:max=1622000-30000*a-90000*b;30000*a+90000*b=1622000;感谢您对网站建设的支持和参与15欢迎光临中国数学建模网a/b=2.5;gin(a);gin(b);end感谢您对网站建设的支持和参与16