全国数学建模大赛B组问题详解.docx

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

2、的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 2011 年 9 月 12 日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：交巡警服务平台的设置和调度摘要“有困难找警察”，是家喻户晓的一句流行语。警察肩负着刑事执法、治安管理、交通管理、

3、服务群众四大职能。为了更有效地贯彻实施这些职能，需要在市区的一些交通要道和重要部位设置交巡警服务平台。本文通过定性与定量分析、建立优化模型，为交巡警服务平台的设置和调度提供参考。在第一个问题中，选择Dijkstra最短路径算法，利用Matlab软件，先根据城区A交通路口的路线，求出表示各节点之间是否直接相连的0-1矩阵，然后根据城区A各节点坐标求出城区A各节点距离的权值矩阵（若两节点内无路则权值为无穷大），接着把权值矩阵化为最短距离矩阵。根据需要变化最短距离矩阵，建立0-1规划模型，目标是使得出警时间最短（转化为出警距离最短计算），列出最优化方程，最后利用Lingo软件进行求解，得出服务平台管

4、辖路口节点以及堵截路口的最合理方案。综合考虑交巡警服务平台的发案率和出警时间，采用动态加权平均的方法算出各个交巡警服务平台的忙碌值。然后进行排名。取大于平均值的前九名，在城区A增加25个服务平台时，综合这些节点周围交通节点的密集程度，决定在A区增加三个服务平台，分别为A20附近的节点90（440.5，381.5），A1、A2和A3区域内的节点67（401，359），A4和A5区域内的节点56（354，374）。在第二个问题中，首先对各城区现有平台设置的合理性进行评估。引入负荷距离法、方差分析法，求得方差、偏差距离、单位平台处理案件数等参数，得出结论：城区C、F服务平台的负担太大，而且警力配置不

5、均匀；城区D、E服务平台的地理分布与发案的地理分布相差较大，不能及时赶到发案地点。再针对各个地区的不同情况（人口、面积、发案率、平台分布疏密程度），经过科学分析，得出方案为：C区增加节点305（200，487）、节点300（206，507）、节点207（333，511）为三个新服务平台，F区增加节点506（358，195）、节点522（371，244）为两个新的服务平台；D区中位于坐标为（70，377）的服务平台D3移动到节点360（76.355），E区中位于坐标为（90，198）的服务平台E15移动到节点422（74，198）。最后通过比较调度前后的该城区的偏差距离、方差、单位平台处理按键数

6、的变化，评估解决方案的合理性。在围堵犯罪嫌疑人的时候，采用画树状图的方法，以三分钟为一个层次，结合概率知识。无论他选择从哪条路出城，得出的围堵方案都能在报警后六分钟之内抓住犯罪嫌疑人。具体方案为：第一个三分钟出动服务平台A5、A6、A10、A15、A16、A2、A3、A4、A17、C8、C6、C4、C7和F1，分别派往节点5、6、10、15、16、3、55、60、41、232、244、240、242、561进行围堵。第二个三分钟出动服务平台C2、C3、D1和 D2，分别派往节点248、168、349、369进行围堵。如果第一个三分钟时已经围堵到了犯罪嫌疑人，那就不用出动第二个三分钟的四个平台，

7、可以节省警力，而且能确保抓住犯罪嫌疑人。关键字：最短路径算法 0-1规划 负荷距离法 方差 树状图1.问题的重述“有困难找警察”，是家喻户晓的一句流行语。警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能。为了更有效地贯彻实施这些职能，需要在市区的一些交通要道和重要部位设置交巡警服务平台。每个交巡警服务平台的职能和警力配备基本相同。由于警务资源是有限的，如何根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题。以给出的条件为例，一共有五个问题需要解决。问题一：为各交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽

8、量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地。问题二：对于重大突发事件，需要调度全区20个交巡警服务平台的警力资源，对进出该区的13条交通要道实现快速全封锁。实际中一个平台的警力最多封锁一个路口，请给出该区交巡警服务平台警力合理的调度方案。问题三：根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实际情况，拟在该区内再增加2至5个平台，请确定需要增加平台的具体个数和位置。问题四：针对全市（主城六区A，B，C，D，E，F）的具体情况，按照设置交巡警服务平台的原则和任务，分析研究该市现有交巡警服务平台设置方案（参见附件）的合理性。如果有明显不合理，请给出解决方案。问题五：

9、如果该市地点P（第32个节点）处发生了重大刑事案件，在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警服务平台警力资源的最佳围堵方案。2.模型的假设（1）交通路口之间的路线都是直线；（2）每条路都是双向的；（3）每个交巡警服务平台的警力相当；（4）每个交巡警服务平台出警速度均为60km/h，不存在堵车等现象；（5）每个交巡警服务平台最多只能封锁一个路口；（6）犯罪嫌疑人对逃跑路线的选择是随机的；（7）犯罪嫌疑人逃跑时所行路线不会重复；（8）犯罪嫌疑人最终目的为逃出该市，不会在城区内躲藏；（9）犯罪嫌疑人行车速度与警车速度一致，同为60km/h；（10）警方通

10、讯时间忽略不计。3.模型的符号说明符号意义节点的忙碌值案发地理重心的横坐标案发地理重心的纵坐标偏差距离从节点到节点的权值城区A每两个节点的最短距离矩阵参考附件2全市交通路口节点数据的第个节点节点的横坐标节点的纵坐标城区A出入口节点与20个交巡警服务平台的最短距离矩阵节点到节点的距离平台地理重心的横坐标平台地理重心的纵坐标节点的出警时间城区各个服务平台案件处理数的方差节点的案发率城区A各节点距离的权值矩阵平台是否管辖节点（是则为1，否则为0）平台是否负责堵截节点（是则为1，否则为0）4.对问题的分析问题一中，我们首先明确最后得出的结果是城区A内20个交巡警服务平台与交通路口的路线的对应关系。路线

11、即两个路口节点之间的部分，所以我们将问题转化为求各交巡警服务平台与路口节点的对应关系，则服务平台管理从平台到该节点的最短路径。对Matlab对原始数据进行处理，用传统的Dijkstra最短路径算法求出城区A每两个路口节点的最短路径的距离。然后运用线性0-1规划模型，列出优化方程，用Lingo软件进行求解。问题二中，可以参考对问题一的分析，最后得出结果应为城区A内20个交巡警服务平台与13个出入市区的路口的对应关系。由于一个平台的警力最多封锁一个路口，也就是每个路口必须指派一个平台的警力，属于指派问题，可以用匈牙利算法解决。匈牙利算法是0-1整数规划的特殊形式，所以可以列出优化方程用Lingo软

12、件求解。问题三中，经过统计各个交警服务平台的工作量和出警时间，可以看出有的服务平台的工作量比较大，而有的服务平台出警时间过长。为了有一个统一的评价，我们对工作量和出警时间做了动态加权平均，建立出综合评价模型。然后对服务平台排名，排在前几名的就是工作量和出警时间相对较大和较长的。所以我们就应该在这些服务平台周围增设服务平台。问题四中，从三个角度评估平台分布的合理性。第一个角度考虑到警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能，希望服务平台尽量与相应案发地点的距离达到最小。这里我们参考负荷距离法，引入“平台地理重心”和“案发地理重心”两个自定义概念，两者分别代表平台分布位置与案发地点分

13、布位置，以这两个重心之间的距离（称为“偏差距离”）为评价标准，距离越小，则该地的平台分布越合理。分别求出六城区的偏差距离，并结合各城区的面积与人口，进行讨论、评估。第二个角度通过求解每一个城区各个服务平台案件处理数的方差，然后比较各城区之间的方差，方差越大表示该区警力分配越悬殊，以此评估服务平台分布的合理性。第三个角度简单地比较各城区单位平台处理案件数，然后判断那些城区警力资源紧张，需要增设服务平台。最后综合以上三个角度，作出相应的对策。5.模型的建立与求解 求解的过程分为三部分。第一部分是对原始数据进行处理，作出求解时可以直接使用的数据表格。第二部分是针对问题的分析建立模型，得出优化方程。第

14、三部分利用软件对模型进行求解，得出最终结果。5.1问题一5.1.1 数据的处理 首先从附件2提取三部分的数据：1、城区A各路口节点（共92个）横纵坐标；2、涉及城区A路口节点的路线；3、20个交巡警服务平台对应的路口节点标号。 根据以上数据，根据92个路口节点的横纵坐标，制作每两个节点之间距离的矩阵。如果两节点之间没有路可以贯通则用表示，得出矩阵W。为矩阵W中的元素，为节点i到节点j的距离，则有： 将该赋权图的权值矩阵W输入，按照Dijkstra方法，反复使用迭代公式：，就可以得到最终结果。即为最短距离矩阵，每个数值都代表两点之间最短路径的距离。这个过程由Matlab软件实现，具体算法参考附录

15、1.然后得出服务平台与城区A各路线节点的对应最短距离矩阵D。5.1.2 模型的建立 目标为让交巡警从服务平台到每一个路口节点的总距离最小，然后还要满足每个结点只需要一个交巡警服务平台负责。运用线性0-1规划模型，如果服务平台位于的交巡警负责点，则记为1，否则为0.再结合5.1.1中得到的最短距离矩阵D，可以得出如下数学表达式： (1)s.t. (2) (3) 根据题意，交巡警服务平台所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地，那么当事发地点利服务平台的距离在： 千米以内的情况下才可以满足条件。又知地图距离和实际距离的比例是1:100000,即1

16、毫米对应100米，则有： (4)5.1.3 模型的求解 将以上数学语言转化为计算机语言，将程序输入Lingo软件求解。具体程序参考附录2，求解结果如下（仅保留非零变量），见表1。表1 城区A各交巡警平台的管辖范围交巡警平台编号管辖范围编号A11，67，68，69，71，73，74，75，76，78A22，39，40，43，44，70，72A33，54，55，65，66A44，57，60，62，63，64A55，49，50，51，52，53，56，58，59A66A77，30，32，47，48，61A88，33，46A99，31，34，35，45A1010A1111，26，27A1212，25A

17、1313，21，22，23，24A1414A1515，28，29A1616，36，37，38A1717，41，42A1818，80，81，82，83A1919，77，79A2020，84，85，86，87，88，89，90，91，925.2 问题二5.2.1 数据的处理 参考附件2全市出入口的位置中出入A的路口标号（一共13个交通节点），以及5.1.1中得到的最短距离矩阵D，取出在A城区内出口节点与20个交巡警服务平台相对应的最短距离矩阵E。5.2.2 模型的建立用0-1规划模型，目标是让交巡警以最短的时间到达封锁点。因为交巡警速度是一定的，所以根据：其中s为距离，v为速度，t为时间。则求得距

18、离最小的情况即为最优解。数学表达式如下： (5)s.t. (6) (7) (8)5.2.3 模型的求解用Lingo软件编程(程序参考附录3)，运行求解，得出一下结果，见表2。表2 A区出口与对应的封锁平台序号出入A区的路口标号负责交巡警平台编号112A12214A9316A16421A14522A10623A13724A11828A15929A81030A71138A21248A51362A4说明：交巡警服务平台A14虽然离交通节点14最近，但为了使全局最优，我们选择了让平台A14去封锁更远一些的交通节点21，服务平台A13到交通节点22和23的距离都比较近，但由于一个服务平台最多只能封锁一个

19、路口所以选择了A13去封锁较远的交通节点23。5.3 问题三 5.3.1 数据的处理 首先根据问题一求出的20个服务平台的分布范围和题目附表中的各个点的发案率，求出每个服务平台每天的工作量，即所管辖的范围内的发案率之和，见表3。 然后再根据每个服务平台跟它所管辖的交通节点的距离之和求出每个平台的出警时间，因为速度一定，可以用距离表示时间的长短，见表3表3 平台的案发率与出警时间平台序号发案率出警时间110.389.7691829.798.1174735.669.0096346.669.2448959.7113.145562.5079.684.842998517.5770198.240.7755

20、8101.60114.625.4330312417.88854138.564.99225142.50154.8104.523716551.32331175.318.34886186.130.94912193.414.320992011.5121.93685.3.2模型的建立和求解 (1) 确定增加平台的个数 首先对发案率和出警时间进行动态加权平均。对于一个服务平台来说工作量大要比出警时间长更重要，更能影响服务平台的工作效率。所以我们将工作量的权重定为0.6，出警时间的权重定为0.4。计算他们的忙碌值。 计算公式为： 计算出他们的综合排名，如表4所示。表4 各服务平台的忙碌值和排名平台A1A2A

21、3A4A5A6A7A8A9A10忙碌值42.0845.0630.9631.6551.071.5039.6910.0321.230.96排名53972186151120平台A11A12A13A14A15A16A17A18A19A20忙碌值12.939.5531.091.5044.6823.5210.5116.037.76855.67排名13168194101412171 我们由表4可以知道，交巡警服务平台A20、A5、A2、A15、A1、A4、A7、A13、A3这些点是排名靠前的，所以拟定在这些点周围增加平台。我们将这些点表示在图1上，我们增加的平台最好能减轻至少一个点的工作量和出警时间，所以我

22、们将这些点分为六个区域。图1 A区的交通节点与平台设置示意图由图1可知，我们可将排名靠钱的点分为几个区域，A20与其所管辖的范围为第一个区域，A1、A2和A3 及其他们管辖范围为第二个区域，A4和A5及其管辖范围为第三个区域，A15为第四个区域，A7为第五个区域，A13为地六个区域。再综合来看A15 虽然排名靠前，但由于它管辖的区域节点数太少，而且发案率比较低，所有在其周围增加服务平台成本高。A7、和A13的排名比较靠后，管辖的范围也不是很大，所以决定也不再在这两个区域增加平台。所以为了尽量减少成本，而且使忙碌值尽量小，所以决定只在A区增加三个平台最好。（2）确定增加平台的位置要确定增加的平台

23、的位置，我们采用重心法。重心法是一种选择销售中心位置，从而使销售成本降低的方法。它把销售成本看成运输距离和运输数量的线形函数。此种方法利用地图确定各点的位置，并将一坐标重叠在地图上确定各点的位置。这里采用这种方法给交巡警服务平台选址。这里用A20管辖的第一个区域为例，将20，84，85，86，87，88，89，90，91，92的节点的横纵坐标以及案发率（），使用公式： （9） （10） 可以得出重点坐标为（443.287，385.1261），同理可得另外两个区域的重点为：（400.0391，351.3574）、（360.6503，377.5613）。然后我们再在途中寻找离这三个重心较近且能分担

24、工作量和出警时间的交通节点，即为我们要求的增加的平台的位置A20附近的节点90（440.5，381.5），A1、A2和A3区域内的节点67（401，359），A4和A5区域内的节点56（354，374）。5.4 问题四5.4.1 负荷距离法分析5.4.1.1 模型的引入 首先讲述负荷距离法的基本思想。单一设施选址中要用到多种分析方法：定性与定量分析方法，以即将定量与定性分析相结合的选址度量法等方法。负荷距离法就是一种单一设施选址的方法。 负荷距离法（load-distance method）的目标是在若干个候选方案中，选定一个目标方案，他可以使总负荷（货物、人或其他）移动的距离最小。我们首先定

25、义“平台地理重心”，它是一个城区所有交巡警服务平台所在节点的横、纵坐标分别取其平均值得到的一个坐标，它代表该城区服务平台的平均位置。 然后我们参考负荷距离法以及加权平均法，定义“案发地理重心”。在一个城区中，以各个路线节点的案发率为权数，对该城区所有路线节点的横、纵坐标进行加权平均，得到一个坐标，它代表该城区发案的平均位置。最后我们求得城区内平台地理重心与案发地理重心的绝对距离，以该距离的大小评价现有各城区交巡警服务平台设置方案的合理性。5.4.1.2 模型的建立与求解 对于任意一个城区内任意一路线节点，它的重点坐标为。直接引用公式（9）、（10）。其中为节点的案发率，所求得的坐标为案发地理重

26、心。对于一个共有n个服务平台的城区内，平台节点（i=1.n）,它的重点坐标为。 （11） （12）所求得的坐标为平台地理重心。各区的偏差距离.根据上述计算方法，可以得出如下表格：表5 各城区平台地理重心与案发地理重心分布情况城区平台地理重心案发地理重心偏差距离A341.40344.40352.01348.4311.35B150.3196.44154.31101.366.35C269.47438.82257.51442.2912.45D61.00353.3372.36338.4218.75E187.57206.33174.54207.3513.06F388.59272.91379.88267.1

27、910.425.4.1.3 模型的评价根据表格，理论上时希望所有的服务平台都向特定的方向（从平台地理重心到案发地理重心的方向）都移动相应的偏差距离。但是由于实际再该处不一定有节点，而且落实到局部不一定适用。该模型给移动平台的方向与距离提供一定的参考。5.4.2 城区内各服务平台案件处理数方差的分析5.4.2.1 模型的建立在概率论和统计学中，一个随机变量的方差描述的是它的离散程度，也就是该变量离其期望值的距离。在概率论和统计学中，一个随机变量的方差描述的是它的离散程度，也就是该变量离其期望值的距离。首先明确目标是求每一个城区各个服务平台案件处理数的方差，并把它们的方差进行比较，对现有交巡警服务

28、平台设置方案的合理性进行评价。这里以城区A为例：参考城区A各交巡警平台的管辖范围（表1），将每个服务平台（A1-A20）平均每日处理案件总数目作出统计，结果如表6所示。表6 城区A各交巡警服务平台每日处理案件数量交巡警平台序号平均每日处理案件数量A110.3A29.1A35.6A46.6A59.7A62.5A79.6A85A98.2A101.6A114.6A124A138.5A142.5A154.8A165A175.3A186.1A193.4A2011.5然后求出所有服务平台平均每日处理案件总数目的方差，结果为：7.8185。用同样的方法得出城区BF各个服务平台案件处理数的方差，得出如下表格：表7 各城区服务平台案件处理数的方差城区ABCDEF7.818513.532530.328811.346716.886424.28025.4.2.1 模型的评价这里所求的城区各个服务平台案件处理数的方差，它代表一个城区内所有服务平台处理案件数目的离散程度，这个方差值越大，则代表这个城区的警力分配不均衡，应该进行调配。如图我们看到城区C、F的值很大，说明这两个城区有的服务平台过于繁忙，有的