(论文)关于警力的分布问题

题目：关于警力分布的问题摘要本文针对某区域多个学校周围警力分布的问题，在执勤点限定范围不同的情况下分别建立了数学模型，求得所需最少警员数及相应的执勤点布设方案。问题一将执勤点的布设限定在95个标志点上。首先，我们利用图论中的相关理论分析引入最短路径长矩阵，并由Floyd算法求得最短路径长矩阵；然后，在假设同一警员的辖区内接连发生两起险情间的时间间隔足够长的前提下，我们以配置最少人数的警员为目标，以学校的安全保障和执勤点布设位置的限定为约束条件，通过相关0-1变量的引入，最终建立模型一，并直接利用相应的软件求得所需最少警员数为20名，同时，求解结果还给出了一种最少警员配置下的执勤点布设方案。问题二依然将执勤点的布设限定在95个标志点上，但要求给出具体的执勤点布设方案并进行评价；我们通过分类讨论和枚举，求得最少警员配置下的执勤点布设共有774144种可行方案。为此，首先我们引入优化指标f（所有的执勤点到其所辖学校的最短路径长之和），并认为f越小时方案越优，然后以配置最少人数的警员和求f最小值为优化目标，以学校的安全保障和执勤点布设位置的限定为约束条件，建立起一个双目标规划模型，最后运用分支定界的思想分类讨论、枚举，即获得一种最优的执勤点布设方案。在方案的评价中，我们引入评价指标缺警数，并利用蒙特卡罗模拟仿真进行多次模拟，求得各警员辖区内可能同时出现多起险情（或多起险情发生间隔很短）时该分配方案的平均缺警数，其中在某一时段第一类学校发生险情的概率为0.5，第二类学校发生险情的概率为0.7时，平均缺警数为3，由此可见虽然假设了同一警员的辖区内接连发生两起险情间的时间间隔足够长，但求得的最优分配方案能够较好的应对同一警员辖区内可能同时发生多起险情的情况。问题三将执勤点的布设限定在了道路上，而不再是有限个标志点，这虽然增大了执勤点布设的灵活性，也利于减少配置警员人数，但也给模型的建立和求解带来巨大困难。我们在满足执勤点布设位置的精度要求下增设标志点，认为直线道路可由一系列的标志表示，将问题3转化为与问题1、2类似的问题，最终建立于模型一类似的模型。关键词：警员配置执勤点布设图论最短路Floyd算法蒙特卡罗模拟标志点1一问题重述今年3月23日早晨，福建省南平市实验小学多名无辜学生在校门口被犯罪分子砍杀。该起重大恶性伤害事件引起了某市市委、市政府领导的高度重视，立即召集市公安局、教育局、行政执法局等有关部门和单位，召开加强校园周边特殊时段安全防范工作紧急会议，研究确定了加强校园周边安全防护工作的若干意见。根据要求，公安部门要将学校安保工作纳入综合控制体系，加强社会嫌疑人员监控与防范。继续做好和落实公安部推出的维护校园及周边治安秩序“八条措施”。要在上下学高峰时段统筹派遣警力值勤护卫，加强校园周边巡逻与保卫工作。在学生、幼儿上下学的重点时段，各所中小学、幼儿园附近道路上安排警员执勤点。要做好应急处置工作，对学校险情进行快速反应，及时处置。现有某区域内学校分布如图，设各标志点之间的道路为直线段。假设警员的执勤点布置在标志点，在接警后能以200米/分的速度赶往现场，根据学校人数的规模分类，各类学校要求尽可能在1分钟之内到达，第2类学校要求尽可能在2分钟之内能有第二名警员到达。1．至少需要多少警员？2．选择合理的执勤点位置，给出方案的评价。3．若执勤点布置不限定在标志点，而是限定在道路上，重新讨论上述问题。二基本假设1.假设在同一警员的辖区内接连发生两起险情间的时间间隔足够长。2.假设警员接警后都能即刻以预定速度赶到事发现场，途中无交通堵塞等意外情况出现。3.假设警员到达事发现场后能迅速控制局面，并能顺利完成每一次任务。4.假设险情发生现场可以用一个点的坐标表示，且险情只在学校所在标志点发生。2三符号说明iv或jv：图G中编号为i或j的的顶点，其中,1,2,...,95ijijd：iv到jv的距离，,1,,95ij,即有22()()ijijijdxxyyijw：iv到jv的直达距离，若iv不能直达jv，则令ijw=，,1,,95ijija：iv到jv的最短路径长，1,,95,1,,95ijin：在标志点iv处安排的警员数四问题分析4.1对问题的总体分析题目要求在一个学校分布比较集中的区域内布置警员执勤点，达到两个目标：①能及时处理各类学校险情，其中，对于各类学校都要求警员接警后尽可能在1分钟之内到达，而对于第2类学校还要求尽可能在2分钟之内能有第二名警员到达；②优化警力资源的配置，在问题1中即是要使配置的警员数量达到最少，在问题二中，即要在警员人数配置最少的前提下，确定一种较优的执勤点布设方案，能方便警员执勤及警员间的工作协调。布置警员执勤点的首要目标就是能及时处理各类学校险情，保障安全，而优化警力资源的配置是以保障安全为前提的；因此，我们可以以优化警力资源的配置作为目标，以保障安全和执勤点布设位置的限定作为约束条件进行数学模型的建立(如下图所示)。模型建立的思路4.2对问题1的分析问题1将警员执勤点的布设位置限定在95个标志点上，分别要求求得所需的最少警员数。目标：优化警力资源的配置约束条件:①保障安全；②执勤点布设位置的限定建立模型3分析题图及所给各标志点的坐标数据，由图论相关理论可知各标志点及其间的直线道路可构成一个图，记为((),())GVGEG,其中()VVG是由所有标志点组成的顶点集，()EEG是以各标志点间的直线道路为元素的边集；由于((),())GVGEG既无环又无平行边，因此G是一个简单图。若以iv到jv的直达距离ijw对图((),())GVGEG的边赋权，就得到一个无向赋权图(,,)GVEW，而由该无向赋权图的邻接矩阵W我们可以求得每对顶点之间的最短路（求每对顶点间的最短路可用Floyd算法实现求解），最终建立最短路径长矩阵ijA。为达到优化目标，我们可根据最短路径长矩阵ijA分别搜索出距离第一类学校不超过200米的标志点和距离第二类学校不超过400米的标志点，从而在满足约束条件的前提下同时确定最少配置的警员数。4.3对问题2的分析问题2仍然是将警员执勤点的布设位置限定在95个标志点上，要求我们选择合理的执勤点位置，并给出方案的评价。在问题1中，我们已经分析了如何求得最少警员数的思路，但并未具体分析最少人数警员的具体执勤点位置安排，而这样的可行方案将可能有很多种，因此，我们的任务就是从多种配置最少警员的可行方案中寻求一种最优方案或较优方案；运用分支定界的思想，我们将能获得一种最优或多种较优执勤点布设方案。因为求得的分配方案是在假设1（假设在同一警员的辖区内接连发生两起险情间的时间间隔足够长）的条件下求得的，而实际中这一假设并不容易完全满足。于是，对于方案的评价，我们可以用计算机模拟仿真定量求出在学校险情的发生不满足假设1时，该分配方案缺少的警员数，从而来评价方案的优劣。4.4对问题3的分析问题3要求将执勤点限定在道路上，重新确定所需的最少警员数和警员执勤点。尽管将执勤点的布设限定在道路上，大大增加了执勤点布设的灵活性，也利于求得更少人数的警员配置，但给布设方案的确定带来了很大的困难。容易知道，每一条直线道路是由无数个连续的点组成的集合，但在满足一定精度要求的前提下，我们可以将一条直线道路看作按一定距离沿该直线道路分布的有限多个点组成。由此，我们可以在直线道路上增设若干标记点（足够多以满足求解精度要求），再将执勤点的布设限定在所有标记点上，这样就可把问题3转化为了问题1、2,只是问题3中标志点的数量要多，于是，模型3的模型建立也就迎刃而解了。五模型的建立与求解5.1问题1的模型建立与求解45.1.1数据的预处理⑴根据两点间的距离公式22()()ijijijdxxyy（式5.1）我们可以将标志点的坐标转化为各标志点间的距离，建立矩阵9595ijDd。⑵将所有的标志点从B～Z,1A～1Z,2A～2Z,3A～3R依次记为1v,2v,…,kv,…,95v⑶由题图可得学校分类并标记如下表所示:5.1.2模型一的建立⑴无向赋权图(,,)GVEW的邻接矩阵依据4.2节中问题分析的思路，我们首先由各标志点及标志点间的直线道路建立一个无向赋权图(,,)GVEW，其权值取为iv到jv的直达距离ijw，且有ijijdw(,),(,)ijijijijijijvvevvEvvevvE可直达,即不可直达即（式5.2）则可得无向赋权图G的邻接矩阵9595()ijWw：1111(9595)nmmn⑵每对顶点间的最短路径长设ija为iv到jv的最短路径长，则由ija构成最短路径长矩阵(9595)()ijAa：11119595nmmnaaAaa⑶目标函数的确定设在标志点iv配置警员in名，则目标函数为951miniin（式5.3）⑷约束条件的确定①设集合kAS由所有与学校ksc(1,2,...,19)k间最短路径长小于200米的标记点组成，则引入0-1变量kib：10kibikikvVSvVS时时（式5.4）5由于各类学校都要求警员接警后尽可能在1分钟之内到达，则kAS中至少有一个标志点被配置了警员，即有约束条件一：9511kiiibn(1,2,...,19)k（式5.5）②设集合kPS由所有与第二类学校ksc(14,15,...,19)k间最短路径长小于400米的标记点组成，则引入0-1变量kic：10kicikikvPSvPS时时（式5.6）由于第二类学校还要求尽可能在2分钟之内能有第二名警员到达，则kPS中至少有一个标志点被配置了警员，即有约束条件二：9512kiiicn(14,15,...,19)k（式5.7）③标志点iv处配置警员数in必须大于或等于0，即有约束条件三:0in(1,2,...,95)i（式5.8）④约束条件四:式5.4。⑤约束条件五：式5.6。⑸最终模型的确定综上所述，我们建立问题一的数学模型如下：obj：951miniin9519511,(1,2,...,19)2,(14,15,...,19)..0,(1,2,...,95)0,10,1kiiikiiiikikibnkcnkstnibc5.1.3模型一的求解⑴由式5.1所示两点间距离公式，利用MATLAB进行简单的编程即可求得iv到6jv的距离ijd，建立矩阵9595ijDd。⑵由式5.2可求得赋权邻接矩阵9595()ijWw中所有元素的值。⑶由Floyd算法求解iv到jv的最短路径长ija。对于无向图，0A是对称矩阵，即ijjiaa。Floyd算法又称插点法，其基本思想是：递推产生一个矩阵序列01,,...,,...,knAAAA，该序列中(,)kAij表示表示从顶点iv到jv的路径上所经过的顶点序号不大于k的最短路径长度；计算时用迭代公式111(,)min((,),(,)(,))kkkkAijAijAikAkj，k是迭代次数，,,1,2,...,ijkn；最后，当kn时，nA即是各顶点之间的最短通路值。针对本问题，Floyd算法求解iv到jv最短路的具体步骤如下：step1:给矩阵(9595)()ijAa和9595()ijRr赋初值，令ijijaw，ijrj，1k；其中9595()ijRr用来存放每对顶点之间最短路径上所经过的顶点的序号，k为迭代次数。step2:更新ijijar和，对所有ij和，若ikkjijaaa，则ijikkjaaa，ijrk。step3:若95k，停止迭代，此时(9595)()ijAa即为所求最短路径长矩阵；否则1kk，转step2。⑷根据最短路径长矩阵，利用MATLAB搜索出所有集合kAS(1,2,...,19)k和kPS(14,15,...,19)k包含的元素，列成表格如下：⑸模型的最终求解在完成以上求解后，根据5.1.2中最终建立