搜索
2011高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：B我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名)：1.2.3.指导教师或指导教师组负责人(打印并签名)：日期：2011年9月12日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：1交巡警服务平台的设置和调度摘要“有困难找警察”，是家喻户晓的一句流行语。警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能。为了更有效地贯彻实施这些职能，需要在市区的一些交通要道和重要部位设置交巡警服务平台。本文通过定性与定量分析、建立优化模型，为交巡警服务平台的设置和调度提供参考。在第一个问题中，选择Dijkstra最短路径算法，利用Matlab软件，先根据城区A交通路口的路线，求出表示各节点之间是否直接相连的0-1矩阵，然后根据城区A各节点坐标求出城区A各节点距离的权值矩阵（若两节点内无路则权值为无穷大），接着把权值矩阵化为最短距离矩阵。根据需要变化最短距离矩阵，建立0-1规划模型，目标是使得出警时间最短（转化为出警距离最短计算），列出最优化方程，最后利用Lingo软件进行求解，得出服务平台管辖路口节点以及堵截路口的最合理方案。综合考虑交巡警服务平台的发案率和出警时间，采用动态加权平均的方法算出各个交巡警服务平台的忙碌值。然后进行排名。取大于平均值的前九名，在城区A增加2~5个服务平台时，综合这些节点周围交通节点的密集程度，决定在A区增加三个服务平台，分别为A20附近的节点90（440.5，381.5），A1、A2和A3区域内的节点67（401，359），A4和A5区域内的节点56（354，374）。在第二个问题中，首先对各城区现有平台设置的合理性进行评估。引入负荷距离法、方差分析法，求得方差、偏差距离、单位平台处理案件数等参数，得出结论：城区C、F服务平台的负担太大，而且警力配置不均匀；城区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）。最后通过比较调度前后的该城区的偏差距离、方差、单位平台处理按键数的变化，评估解决方案的合理性。在围堵犯罪嫌疑人的时候，采用画树状图的方法，以三分钟为一个层次，结合概率知识。无论他选择从哪条路出城，得出的围堵方案都能在报警后六分钟之内抓住犯罪嫌疑人。具体方案为：第一个三分钟出动服务平台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进行围堵。如果第一个三分钟时已经围堵到了犯罪嫌疑人，那就不用出动第二个三分钟的四个平台，可以节省警力，而且能确保抓住犯罪嫌疑人。关键字：Dijkstra最短路径算法0-1规划负荷距离法方差树状图21.问题的重述“有困难找警察”，是家喻户晓的一句流行语。警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能。为了更有效地贯彻实施这些职能，需要在市区的一些交通要道和重要部位设置交巡警服务平台。每个交巡警服务平台的职能和警力配备基本相同。由于警务资源是有限的，如何根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题。以给出的条件为例，一共有五个问题需要解决。问题一：为各交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地。问题二：对于重大突发事件，需要调度全区20个交巡警服务平台的警力资源，对进出该区的13条交通要道实现快速全封锁。实际中一个平台的警力最多封锁一个路口，请给出该区交巡警服务平台警力合理的调度方案。问题三：根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实际情况，拟在该区内再增加2至5个平台，请确定需要增加平台的具体个数和位置。问题四：针对全市（主城六区A，B，C，D，E，F）的具体情况，按照设置交巡警服务平台的原则和任务，分析研究该市现有交巡警服务平台设置方案（参见附件）的合理性。如果有明显不合理，请给出解决方案。问题五：如果该市地点P（第32个节点）处发生了重大刑事案件，在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警服务平台警力资源的最佳围堵方案。2.模型的假设（1）交通路口之间的路线都是直线；（2）每条路都是双向的；（3）每个交巡警服务平台的警力相当；（4）每个交巡警服务平台出警速度均为60km/h，不存在堵车等现象；（5）每个交巡警服务平台最多只能封锁一个路口；（6）犯罪嫌疑人对逃跑路线的选择是随机的；（7）犯罪嫌疑人逃跑时所行路线不会重复；（8）犯罪嫌疑人最终目的为逃出该市，不会在城区内躲藏；（9）犯罪嫌疑人行车速度与警车速度一致，同为60km/h；（10）警方通讯时间忽略不计。3.模型的符号说明符号意义iB节点i的忙碌值xC案发地理重心的横坐标yC案发地理重心的纵坐标d偏差距离ijd从节点i到节点j的权值()nD城区A每两个节点的最短距离矩阵3iD参考附件2全市交通路口节点数据的第i个节点ixD节点iD的横坐标iyD节点iD的纵坐标E城区A出入口节点与20个交巡警服务平台的最短距离矩阵ijf节点iD到节点jD的距离xG平台地理重心的横坐标yG平台地理重心的纵坐标ipt节点i的出警时间2S城区各个服务平台案件处理数的方差iV节点i的案发率iW城区A各节点距离的权值矩阵ijx平台i是否管辖节点j（是则为1，否则为0）ijy平台i是否负责堵截节点j（是则为1，否则为0）4.对问题的分析问题一中，我们首先明确最后得出的结果是城区A内20个交巡警服务平台与交通路口的路线的对应关系。路线即两个路口节点之间的部分，所以我们将问题转化为求各交巡警服务平台与路口节点的对应关系，则服务平台管理从平台到该节点的最短路径。对Matlab对原始数据进行处理，用传统的Dijkstra最短路径算法求出城区A每两个路口节点的最短路径的距离。然后运用线性0-1规划模型，列出优化方程，用Lingo软件进行求解。问题二中，可以参考对问题一的分析，最后得出结果应为城区A内20个交巡警服务平台与13个出入市区的路口的对应关系。由于一个平台的警力最多封锁一个路口，也就是每个路口必须指派一个平台的警力，属于指派问题，可以用匈牙利算法解决。匈牙利算法是0-1整数规划的特殊形式，所以可以列出优化方程用Lingo软件求解。问题三中，经过统计各个交警服务平台的工作量和出警时间，可以看出有的服务平台的工作量比较大，而有的服务平台出警时间过长。为了有一个统一的评价，我们对工作量和出警时间做了动态加权平均，建立出综合评价模型。然后对服务平台排名，排在前几名的就是工作量和出警时间相对较大和较长的。所以我们就应该在这些服务平台周围增设服务平台。问题四中，从三个角度评估平台分布的合理性。第一个角度考虑到警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能，希望服务平台尽量与相应案发地点的距离达到最小。这里我们参考负荷距离法，引入“平台地理重心”和“案发地理重心”两个自定义概念，两者分别代表平台分布位置与案发地点分布位置，以这两个重心之间的距离（称为“偏差距离”）为评价标准，距离越小，则该地的平台分布越合理。分别求出六城区的偏差距离，并结合各城区的面积与人口，进行讨论、评估。第二个角度通过求解每一个城区各个服务平台案件处理数的方差，然后比较各城区之间的方差，方差越大表示该区警力分配越悬殊，以此评估服务平台分布的合理性。第三个角度简单地比较各城区单位平台处理案件数，然后判断那些城区警力资源4紧张，需要增设服务平台。最后综合以上三个角度，作出相应的对策。5.模型的建立与求解求解的过程分为三部分。第一部分是对原始数据进行处理，作出求解时可以直接使用的数据表格。第二部分是针对问题的分析建立模型，得出优化方程。第三部分利用软件对模型进行求解，得出最终结果。5.1问题一5.1.1数据的处理首先从附件2提取三部分的数据：1、城区A各路口节点（共92个）横纵坐标；2、涉及城区A路口节点的路线；3、20个交巡警服务平台对应的路口节点标号。根据以上数据，根据92个路口节点的横纵坐标，制作每两个节点之间距离的矩阵。如果两节点之间没有路可以贯通则用表示，得出矩阵W。ijd为矩阵W中的元素，ijf为节点i到节点j的距离，则有：,ijij0ijijijfd若与直接相连，若与不直接相连，若将该赋权图的权值矩阵W输入，按照Dijkstra方法，反复使用迭代公式：()(1)(1)(1),,{1,2,...,}d{d,d+d}minkkkkijijikkjijkn，就可以得到最终结果()nD。()nD即为最短距离矩阵，每个数值都代表两点之间最短路径的距离。这个过程由Matlab软件实现，具体算法参考附录1.然后得出服务平台与城区A各路线节点的对应最短距离矩阵D。5.1.2模型的建立目标为让交巡警从服务平台到每一个路口节点的总距离最小，然后还要满足每个结点只需要一个交巡警服务平台负责。运用线性0-1规划模型，如果服务平台位于i的交巡警负责j点，则记ijx为1，否则为0.再结合5.1.1中得到的最短距离矩阵D，可以得出如下数学表达式：209211minijijijDx(1)s.t.52011,1,2,...,;ijixjn(2)01,1,2,...;ijxijn或，(3)根据题意，交巡警服务平台所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地，那么当事发地点利服务平台的距离在：360360千米以内的情况下才可以满足条件。又知地图距离和实际距离的比例是1:100000,即1毫米对应100米，则有：(,)(,)30Dijxij,1,2,...;ijn(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，