第五题A 十队 庞威 王仁和 冯奇.docx

1、第五题A 十队 庞威 王仁和 冯奇2009高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A

2、 我们的参赛报名号为（如果赛区设置报名号的话）： A甲1909 所属学校（请填写完整的全名）： 海军航空工程学院 参赛队员 (打印并签名) ：1. 庞威 2. 王仁和 3. 冯奇 指导教师或指导教师组负责人 (打印并签名)： 孙玺菁 日期： 2010 年 8月 19 日赛区评阅编号（由赛区组委会评阅前进行编号）：2009高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：免费自行车交通系统服务网点布局规划摘要：本文是在

3、针对网点设置和自行车的数量分布问题进行分析的，建立了一个包括四个评价指标的数学模型。对问题一，在已知网点与车辆分布状况之下，我们首先提出了四个评价指标，分别为超市、地铁站点、十字路口和居民区对网点的综合影响，辐射区域的人口密度，人均占有率，距离指标。通过计算这四个评价要素对于网点的影响来来衡量现有网点与车辆分布状况。对问题二，通过在地图中标定120个网点(包括原来的17网点)，然后通过问题一定义的评价指标，筛选选出指标最好的100个网点。通过每个网点所覆盖的总的人数与所有的覆盖区域的总人数的比例来确定每个网点的自行车数。对问题三，在第二问120点的基础上，假设自行车总数从1600到3600变化

4、，借助0-1变量进行线性规划。求出每种自行车数情况的最优解，在定义的新评价指标下找出最优的自行车总数，相应求得网点数及分布情况。关键字：加权法，图像处理，像素法，0-1整数规划1 问题重述城区现有人口15 万，地域面积约22.9平方公里。已知规划中的地铁站有5个，图上AE点，预计高峰时间人流量在4000-5000人/站，其余时间10002000人/站。大型社区有两个，社区C有1.4万人，社区C有2.8万人，其余地区，除山地、湖泊和河流区域外，可以认为人口是均衡分布的。大型超市有三个，预计高峰时间人流量在3000人/座，其余时间1000人/座。最大可能方便居民使用，应优先考虑交通枢纽和地点人流量

5、，根据现实中调查可以推断：早晨在社区周边的网点车辆数较多，下午下班时在地铁站和超市附近网点的车辆数较多。十字路口的人流量一般较大。网点之间的距离一般控制在300米1000米之间。目前该地区现有17个网点，600辆免费自行车。现在我们来解决以下问题1. 设定一个评价标准来衡量现有网点与车辆分布状况。2. 在规划中要在图中增加到100个网点和3600辆车，决定网点位置跟每个网点的车辆数，使在我们的评价指标下达到最优。3.目前只有110万元左右，建设一个网点需5000元，投入一辆自行车的成本约300元，为了尽可能实现主要居民区网点平均间距500米的公共交通体系，并最大程度服务居民，则需要在此地区建立

6、多少个，如何分布网点并确定每个网点的车辆数。本文中所有没有表头的表格均按列读取。2 问题分析与模型建立2.1问题1的分析与模型建立 对网点的设定的影响因素有多种，为了计算得方面，我们在此仅考虑了一下四个主要影响因素，并将其作为4个指标，这4个指标分别为：（1）网点之间的距离越大，距离免费自行车的人的满意度越低，但如果网点之间的距离很小的，会造成资源的浪费，因此自行车网点之间的距离应该合适才行即为题目要求的300到1000米，则设定以某个网点为对象距离他最近的网点的距离减去650米的绝对值的倒数为指标，越接近650则越大，越满足距离要求。（2）考虑超市、地铁站点、十字路口和居民区对于设定自行车网

7、点的综合影响，设为综合影响系数，其中越大，表明对于设定自行车网点的影响越大，越有必要在此点设定自行车网点，以不同权重分别乘以网点到超市、地铁站点、十字路口和居民区距离的倒数的和为指标来表征影响系数。（3）由于自行车网点之间的距离应设定在300到1000米之间，我们在本文中设定为1000米，并以各个网点为圆心，以1000米为半径画圆来求网点辐射区域中的人口密度和人均占有率两个指标。但是由于本题中的含有山地河流，在计算时需要排除这些面积，非常繁琐，因此我们用与圆内接正方形和圆外切正方形面积相同的正方形来代替此圆，此正方形的两边与地图的两边界平行，见图1。图1 以网点为圆心1000米为半径的辐射区域

8、具体的计算公式如(1)所示： 为圆的半径。 （1）然后我们定义辐射区域的人口密度为在每一个圆中人的总数除以圆的面积（除去河流和山川的面积）。要计算圆中的的人口总数，需要像素人口密度，定义像素人口密度为每个不同像素对应的不同的人数。因为地图是有每个像素组成的，我们运用Matlab中的Imread命令读出地图的像素，由于不同颜色图的像素不一样，因此我们利用总的像素数减去山地河流湖泊的像素数，得到的便是有人居民居住地方的像素数，然后我们利用每个地区总的人数除以这个地区的像素数得到的便是像素人口密度。为了计算每个自行车网点所覆盖的人数，只需要将这个网点所覆盖的像素人口密度求和即可。用自行车网点覆盖的人

9、口总数除以该网点覆盖的面积，得到一个数，即为辐射区域的人口密度，辐射区域的人口密度越大，则在此地点设定自行车网点越有必要。（4）人均占有率定义为以网点为圆心，1000米为半径网点的上午与下午的自行车数的平均值除以该圆中总的人数。人均占有率越大，则越应该在此设定自行车网点且在这个网点的自行车数就越多。2.1问题二的分析与模型建立要在图中增加到100个自行车网点，由于网点数目比较多，而区域相对较小，因此我们分析可能设定为自行车网点的103点，从题目中可以看出在十字路口和超市、居民区、地铁站附近设为自行车网点，然后从这些可能设定为自行车网点的点中筛选出对于我们评价标准最好的83个自行车网点。1.对距

10、离指标、人口密度和居民区、超市、地铁站点和十字路口分别对于设定自行车网点的影响较易求得；2.计算人均占有率时按网点辐射区域的人口数占所有网点辐射人口总数比例来计算每个网点的自行车数目。按网点距离社区，超市地铁，十字路口距离不同将所有网点分为三类，根据第一问的统计表格求得每一类网点的上下午自行车数变化平均数，进而得到120个网点的上下午自行车数目的变化这样就可以按第一问的指标计算人均占有率。2.3问题三的分析与模型建立 比较第三问与前两问发现要求的变化有三点：1.距离要求是主要居民区网点平均距离为500米；2.限制可用资金为110万元；3.要求尽可能的服务居民。根据不同的要求，在第二问预选的12

11、0个网点的基础上，进行分析。假设自行车总数变化范围为1600到3600，在此基础上设定01变量进行线性规划。因为当自行车数为3600时，建立网点的可用费用只剩下20000元，也就是说只能建4个网点，很明显这是不合实际的，因此车辆数不可能超过3600。 为了尽可能的服务居民，对各个指标加权重后求和，根据对于居民的相关程度不同分布权重，已体现尽可能服务于居民。 3 模型的假设（1）假设居民在选择自行车网点的时候，总是选择离自己最近的自行车网点且如果网点距离居民的距离大于1000，人们则拒绝去此网点使用自行车；（2）假设在山地、湖泊和河流等区域内无居民居住；（3）假设此地区的交通便利，图中只标出了主

12、要的道路，每个网点之间可直线到达；（4）第二问中丁字路口作为十字路口考虑；（5）自行车不会出现意外情况丢失或被盗，损坏后能马上被修好并投入使用。4 符号说明：居民区、超市、地铁站点和十字路口分别对于设定自行车网点的影响系数；：居民区、超市、地铁站点和十字路口分别对于自行车网点的影响权重；：将进行标准化后的数据；：以各个网点为圆心，1000米为半径的圆的近似面积；：分别为自行车网点、地铁站、超市的横坐标；：分别为自行车网点、地铁站、超市的纵坐标；:第个网点到最近居民区、十字路口、超市和地铁站点的距离；：各网点和与之距离最近的网点之间的距离；：对进行标准化后的数据；：山地河流湖泊、居民区及其他地区

13、的像素人口密度；：以各个网点为圆心，以1000米为半径，在这个圆以内的人口总数；：每个网点的辐射区域人口密度；：标准化后的人口密度；：表示每像素点上的人口；：为第站点自行车的个数；：网点辐射区域占有率；：将标准化后的数据；：总的评价指标。5 模型求解5.1问题一的解答5.1.1首先根据图可以读出各个自行车网点和地铁站点以及超市的坐标，如表1、2、3所示。表1 自行车各个网点的坐标（单位：像素）116010532319784368278446451008561010586934297764931489384869906771109346121192382512545702133759914756

14、250154041991618520617636771表2 地铁站点坐标（单位：像素）a429195B510627C510876D4691027E899640表3 超市坐标（单位：像素）510773604254791884（1）由图易知各个网点和与之距离最近的网点，求出这两网点之间的距离如表4所示。表4 各网点和与之距离最近的网点的距离网点网点距离（米）12133321297131721545266542666147567145578105169112431085161192441217497131543514755715134351615906173215由表5可得只有1，2网点之间的距离大

15、于1千米，3,4,5,9,11,17的最短距离小于300，均不满足距离要求。根据分析中的距离指标的定义得 对表4中数据进行标准化处理即 （2）将数据代入（2）式得表5。表5 网点之间距离标准化数据0.16130221410.80830.80830.28440.3860.41670.88480.4167088110.430.49430.3860.49430.2371（2）影响安排自行车网点的因素根据其距居民区、超市、地铁站点以及十字路口不同而不同，设定其权值为0.1，0.2，0.3或0.4，数越大，对相应网点的影响越小，并记为，网点到居民区、超市、地铁站点以及十字路口距离越近，所赋的权重越小，所

16、得乘积的倒数求和作为指标。由此可得影响每个网点的因素的权值见表6。表6 影响每个网点的因素网点居民区十字路口超市地铁站点10.10.20.40.320.10.20.30.430.30.10.20.440.20.10.40.350.30.10.40.260.20.10.30.470.10.20.30.480.40.10.30.290.40.10.30.2100.40.10.302110.40.10.30.2120.20.10.30.4130.20.10.40.3140.20.20.30.3150.30.20.30.2160.30.10.40.2170.30.10.30.3然后分别求出各个网点到最

17、近居民区、十字路口、超市和地铁站点的距离，我们定义影响系数为影响设定自行车网点的程度的量，可有 （3）其中越大，则在此设定自行车网点越合适。并对数据进行标准化处理即 （4）由附录中的程序得表7表7 标准影响系数0.8060.6950.3230.3190.2140.9340.8850.5060.3700.3470.1380.6640.12810.2670.3980.243（3）设山地河流湖泊的像素人口密度为=0，居民区CI的像素人口密度为，居民区CII的像素人口密度为，其他地区为，以各个网点为圆心，以1000米为半径，在这个圆以内的人口总数为记为，每像素点上的人口为，所以辐射区域的人口密度，则

18、（5） （6）由附录中的程序计算得到人均人口密度，对人均人口密度进行标准化处理即 （7）得表8。表8 标准辐射区域的人口密度0.66831.00000.52880.76620.89050.53970.78010.48270.49030.47960.47660.72340.48120.89550.48990.47210.6501（4）定义人均占有率为 （8）为第站点自行车的个数。这里取上午与下午的车辆数平均值。由附录中的程序得到人均占有率，对人均占有率数据进行标准化处理即 （9）得表9。表9 标准化后的人均占有率1.00000.56530.49890.21930.21770.27930.4590

19、0.31230.76860.31430.39540.57300.23500.46290.23080.19960.5217将各个网点对应的评价指标求和即 （10）代入数据得表10。表10 总的评价指标2.6362.482.112.132.032.5111.722.511.551.892.391.342.741.481.481.302.42由表10得总的评价指标值越大，则在该点设定为自行车网点就越好，因此由表10得设定为自行车网点的点有好到差的顺序见表11。表11 设定为自行车网点的点由好到差的顺序13168217114351079141512165.2问题二的解答5.2.1 标定可能设为自行车网

20、点的点如图2所示，图中圆圈表示可能设定为自行车网点的点。我们在十字路口、超市、地铁车站以及居民区附近都设定为自行车网点，但这些网点总和不到83个，从图中可以看出有些网点之间的距离过大，因此我们在这些点之间再设定一些网点，由此得出来的网点总数和原来已经设定好的网点总数为120个，我们根据第一问中设定的标准对这103个网点进行筛选，从中选出83个网点。图2 可能的自行车网点分别计算出图2中各个可能网点及已有网点的坐标，然后分别计算出各个网点距离超市、地铁站、距其最近的居民区及距其最近的十字路口路（包含丁字路口）的距离。 根据第一问的方法再计算出各个网点的人均人口密度、人均占有率、以及网点之间距离，

21、求出总的指标，从中筛选出较好的83个网点。5.2.2 网点各个指标的计算分别测量出可能设定为自行车网点的坐标，根据这些网点的坐标计算出距其最近的十字路口、超市、居民区地铁站点的距离，然后根据距离的大小设定相应的权重，即距离最近的为主要因素0.1，距离较远的为0.2，次之的为0.3，最远的为0.4，权重越大，代表因素对自行车网点的影响越小。由附录得这120个网点的影响系数，见表12表12 120网点的影响系数0.36530.37680.71040.27700.18450.19590.49860.45220.72050.25250.71600.75730.24630.28470.18370.357

22、20.17920.68900.12620.17530.70680.72300.70610.12500.69220.12850.68560.68620.69030.13850.72970.73100.11250.68650.70720.82470.25760.21530.11270.10240.68660.69180.70760.75320.98790.70620.69730.68810.69860.09340.68570.69000.69930.71380.73550.10420.08170.68710.69240.69890.70860.71440.70980.69670.69160.123

23、50.17520.74420.70830.71660.70880.70300.69680.22880.72990.27980.85820.73430.71150.71950.76320.79850.71380.70320.71700.70240.69510.76360.69250.11760.75460.73380.77001.00000.72810.72610.71130.70150.71690.75650.27610.09360.21520.80260.86950.24200.75110.72530.70470.10810.69300.14600.71540.70410.10950.098

24、60.71610.78940.23820.6885由附录中的网点之间的坐标计算得到这120个网点之间的距离指标标准化数据，见表13。表13 120个网点之间的距离指标的标准化数据0.61560.44230.31450.30390.30110.22210.22830.26190.27360.23020.1920.30380.36740.25750.24290.41950.23020.23460.23280.40230.29960.29110.22430.28710.25750.27390.34970.28140.24070.25860.22430.41030.25750.2560.33990.4

25、0230.33350.25130.29920.4220.30380.31660.27530.39530.25810.25180.23910.61560.35450.1920.26510.42020.25270.24880.24490.36220.3560.23280.3090.39120.25750.28090.34410.23910.34050.36020.6230.25860.25180.37240.26180.24880.34050.43520.35450.25130.27340.33390.26180.24490.29250.50490.64380.23650.24230.26190.

26、26660.23750.29250.42020.24070.23650.25760.85390.36220.2560.34970.28140.39530.25760.23520.34410.50650.32220.23460.28720.33350.25810.73380.34410.23350.27530.27390.26510.30190.22210.32980.43250.2335下面计算120各网点的人均占有率，首先规定每个网点的辐射区域的车辆数与此区域的人口数目占所有辐射区域的车辆总数的比例成正比即 （11）按每个网点距离社区，超市地铁，十字路口的距离不同将网点分为三类，及距离社区的距离最近则分为社区类。根据题目提供的统计数据，计算出每一类网点的上午下午车辆数目平均变化率(=1,2,3)。网点的人均占有率为 （12）推导后不难得到 （13）由附录中的程序计算得到120网点之间的人均占有率，见表14。表14 120个网点的人均占有率1.00001.00000.77840.77840.77840.77841.00000.77840.77840.75680.77840.77840.77841.00000.75680.77840.77840.77840.77840.77840.77