交巡警服务平台的设置与调度数模论文

交巡警服务平台的设置与调度数模论文PAGEPAGE28交巡警服务平台的设置与调度摘要本文讨论了交巡警服务平台的设置和调度问题。交巡警在日常道路管理中发挥着不可替代的作用，但警务资源有限，为了使警务资源发挥作用最大必须对有限的交巡警服务平台、平台管辖范围、警务资源等进行合理配置。在问题一第1问中要对交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地。我们将其转化为最短路径问题应用Floyd算法求出A区每个道路结点到交巡警平台的最短路径，通过任意节点归离它最近的服务平台管理的原则将其分配给该交巡警平台管辖。在问题一第2问中我们提出了两个中不同的模型，分别针对总体出警时间和个体出警时间，第一种方案最快完成封锁需，第二种最快完成封锁的时间为8分钟，因此选后者最为最佳解决模型。具体调配方案见表1右侧。在问题一第3问中，我们以方差作为衡量交巡警服务平台的工作量的标准，以出警时间尽量小于三分钟作为控制出警时间的标准，建立优化模型，在MATLAB的计算下得出：再增加5个交巡警服务平台,分别在节点39、48、58、68、87能使工作量均衡量最小，方差从8.43314降到2.89。在问题二第1问中，对各区交巡警服务平台设置是否合理分为两个方向，其一为交巡警服务平台设置数目是否合理，其二为交巡警服务平台设置的地理位置是否合理。对前者采用主成分分析法，对后者沿用问题一中的模型，经分析可知其平台设置不合理，需改进，我们在满足出警时间尽量在3分钟内、保证各交巡警服务平台工作量均衡的原则下提出了解决方案，即增加市区服务平台数，最后结果显示，我们的解决发案是各区平台方差明显下降，超出3分钟出警时间点明显减少。在问题二第2问中，该问题实质为指派问题,使用0-1规划解决交巡警服务平台是否对嫌疑犯进行搜捕的问题,实现三分钟后交巡警可以在较短的时间内搜捕到嫌疑犯。关键字最短路径，Floyd算法，优化，主成分分析法，0-1规划1．问题重述1．1问题背景“有困难找警察”，是家喻户晓的一句流行语。警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能。为了更有效地贯彻实施这些职能，需要在市区的一些交通要道和重要部位设置交巡警服务平台。每个交巡警服务平台的职能和警力配备基本相同。由于警务资源是有限的，如何根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题。1．2问题内容现需建立数学模型解决以下问题：（1）图一给出了该市中心城区A的交通网络和现有的20个交巡警服务平台的设置情况示意图，相关的数据信息见附件2。请为各交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地。对于重大突发事件，需要调度全区20个交巡警服务平台的警力资源，对进出该区的13条交通要道实现快速全封锁。实际中一个平台的警力最多封锁一个路口，请给出该区交巡警服务平台警力合理的调度方案。根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实际情况，拟在该区内再增加2至5个平台，请确定需要增加平台的具体个数和位置。（2）针对全市（主城六区A，B，C，D，E，F）的具体情况，按照设置交巡警服务平台的原则和任务，分析研究该市现有交巡警服务平台设置方案（参见附件）的合理性。如果有明显不合理，请给出解决方案。如果该市地点P（第32个节点）处发生了重大刑事案件，在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警服务平台警力资源的最佳围堵方案。图1：A区的交通网络与平台设置的示意图图2：全市六区交通网络与平台设置的示意图说明：（1）图中实线表示市区道路；红色线表示连接两个区之间的道路；（2）实圆点“·”表示交叉路口的节点，没有实圆点的交叉线为道路立体相交；（3）星号“*”表示出入城区的路口节点；（4）圆圈“○”表示现有交巡警服务平台的设置点；（5）圆圈加星号“”表示在出入城区的路口处设置了交巡警服务平台；（6）附图2中的不同颜色表示不同的区。2．模型假设（1）假设该市每条道路路况相同；（2）假设警车在每条路况中都保持60km/h（即1000m/min）的速度；（3）假设交巡警到达了某条路段的交叉路口节点即控制了该路段；（4）假设交巡警工作地理范围仅为交通网络图中给出的道路范围；（5）假设交巡警处理案件的时间为定值；（6）假设计算交巡警工作量时不考虑行驶时间；（7）假设不同的区互相不管辖，同一个区仅管辖本区内的路口节点；（8）假设嫌疑犯以的速度逃逸。3．符号规定区编号的交巡警平台的坐标，区编号的节点的坐标，编号的交警平台到达编号的节点的距离，单位：m编号的交巡警平台到达编号的节点的时间,单位：min警车速度，单位：1000m/min编号的交巡警平台管辖的道路结点的集合（建立模型时不包括点本身，最终结果中包含点）0-1变量各交巡警服务平台处理的案件次数的方差编号的交巡警服务平台管辖范围内编号点的案发率各节点平均案发率区内交巡警服务平台数4．问题分析该题为典型的图论与资源合理配置的问题，我们运用经典的Floyd算法和动态规划解决分配交巡警管辖范围、警务调动、服务平台设置合理性分析等问题。4．1问题一的分析4．1．1问题一（Ⅰ）的分析问题一（Ⅰ）要为各交巡警服务平台分配管辖范围，我们将其转化为求最短路径问题，若使警车尽可能在三分钟之内到达事发地，即求出距离交警平台3000m之内的道路节点有哪些，若有道路节点距所有交警平台均大于3000m，则采取就近将其分配给距离最近的交警平台。4．1．2问题一（Ⅱ）的分析问题一（Ⅱ）要实现对13个路口的封锁，而实际中一个平台的警力最多封锁一个路口，为实现快速全封锁必须使封锁时间最短，现有两种模型：eq\o\ac(○,1)面向全体交巡警服务平台，实现总体调动时间最短,节省警务资源；eq\o\ac(○,2)面向一个交巡警服务平台，实现一个交巡警服务平台的调动时间最短。在结果中进行比较，选出较优方案。4．1．3问题一（Ⅲ）的分析问题一（Ⅲ）要解决交巡警服务平台的工作量不均衡和出警时间过长的问题，我们以每个交巡警服务平台处理案件次数的方差作为衡量工作量的标准，以交巡警到达管辖范围的最大时间小于作为衡量出警时间的标准，优化求解增加的交巡警平台数以及增加的具体位置。4．2问题二的分析4．2．1问题二（Ⅰ）的分析由上题可知交巡警服务平台设置要求即保证出警时间较短、保证各交巡警服务平台工作量均衡。使用主成分分析法对各区交巡警服务平台设置数目是否合理进行分析；沿用问题一第三问模型对各区交巡警服务平台地理位置的设置是否合理进行分析，并对不合理的设置给出具体的解决方案。4．2．2问题二（Ⅱ）的分析要实现对嫌疑犯的快速搜捕，应确定三分钟之内嫌疑犯可能到达的路口节点，再根据这些节点确定搜捕方案。5．模型的建立与求解5．1问题一的模型建立与求解Ⅰ：分配各交巡警服务平台的管辖范围5．1．Ⅰ．1数据的基本处理①在MATLAB中引用Excel表格将每两个相邻节点（车辆通过时不经过第三个节点）的距离计算出来；②在算法中引入，对于相邻的两个节点记为，对于不直接相连的两个节点记为；假设该图邻接矩阵为,存放各节点之间的距离，其中；；；又该图为无向图，为对称矩阵，。5．1．Ⅰ.2模型的建立要保证警车尽量会在3分钟内到达突发事件现场，则有：eq\o\ac(○,1)找出相应的点对应点：（一）当对应唯一一个点时，则；（二）当对应多个点时，；（三）当没有对应的点时，采取就近原则，。求解结果：；；；；；；；；；；；；；；；；；；；。具体管辖散点图见下图：图3Ⅱ.合理的调度方案5．1．Ⅱ.1数据处理在该题的解答过程中引入0-1变量，对没有经过的相邻节点记为0，对被经过的两个相邻节点记为1，其中：编号为的交巡警服务平台，；：第个出入区的路口标号，，其对应的路口标号为5．1．Ⅱ.2两种模型的建立与比较建立：模型一：面向全体交巡警服务平台，实现总体调动时间最短实质为优化问题，其模型为：模型二：面向一个交巡警服务平台，实现最后完成调度的交巡警服务平台调动时间最短，即最终完成调动时间最短。实质为优化问题，其模型为：经LINGO和MATLAB计算后可得：模型一模型二出入区路口编号时间/min出入区路口编号时间/min211383.9288113624.8413620.35211383.928512482.475833166.02710300.5831512482.47892910.4930610303.293161.532579298.0155105227.7079105227.7079117243.8053117243.8053136230.5126236.477142140131125.977121120144213.265158284.7518158284.75181642110.0060162146.7117表1模型一时间总计：，出警最短时间为；模型二时间总计为，出警最短时间为。比较：针对总体考虑时，模型一时间较短，但其出警时间较长；针对个体考虑时，模型二出警时间较短，但其总体调度时间较长。按照“尽快出警”原则，当发生重大突发事件时，应在最短的时间内封锁路口，模型二比模型一更符合题意，因此我们选择模型二，具体调配方案见表1右侧。Ⅲ确定需要增加平台的具体个数和位置5．1．Ⅲ.1模型的建立若使各交巡警服务平台工作量均衡，则有方差取最小值；若使出警时间较短，则应满足交巡警服务平台到其管辖的节点最长时间尽量小于。方差计算方法：优化模型：对于：①当对应唯一一个点时，则；②当对应多个点时，；③当没有对应的点时，采取就近原则，。求解结果：新增交巡警服务台编号：39、48、58、68、87原始方差：8.4314新增结点后方差：2.89新增交巡警服务台后其管辖节点及时间分别为：编号案发率123456789101112131415，16171819203948586787表2由上表可知，增加新交巡警服务平台之后仅有编号15的交巡警服务平台到达其管辖节点时间较长，大于三分钟。5．2问题二的模型建立与求解Ⅰ合理性分析5．2．Ⅱ.1模型的建立按照要求，交巡警平台的设置主要遵循以下原则[1]:

1、警情主导警务原则：根据管区道路交通流量、拥堵状况、治安复杂情况、发案量高低，科学确定平台管控区域；

2、快速处警原则：城区接警后确保快速到达现场；

3、方便与安全原则：按照醒目、规范，方便群众和确保安全的原则，科学设置平台。

平台设置在遵循上述三大原则的基础上，应当结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况，在充分考虑现有警力和财力并确保安全的条件下，科学确定平台的数量和具体位置。先对题目中给出的各区的交巡警服务平台分配管辖范围，模型沿用问题一第一问。对交巡警平台设置数目的合理性分析（主成分分析）：交巡警服务平台的设置数目主要与该区面积、人口、工作量（案发率）、工作量方差等有关。区域编号指标ABCDEF地域面积/22.00103.00221.00383.00432.00274.00人口/万人60.0021.0049.0073.0076.0053.00区域每天总发案次数117.7566.40187.2067.80151.1796.92工作量均值a/每天发案率5.898.3011.017.5310.088.81工作量方差8.4315.4632.1412.7715.2030.00出警时间大于3的节点数6.004.0043.0012.0029.0034.00表3用表示 区，各区的地域面积、人口数、区域每天总的发案次数、工作量均值、工作量方差、出警时间大于3的节点数的取值分别记作，构造矩阵（1）数据标准化处理：对原始数据标准化处理，将各指标值转化为标准化指标，有：其中：即为第个指标的样本均值和样本标准差。对应地称……标准化指标变量。（2）计算相关系数矩阵:相关系数矩阵，有其中：是第个指标与第个指标的相关系数。（3）计算特征值与特征向量计算相关系数矩阵的特征值及对应的标准化特征向量，其中，由特征向量组成3个新的指标变量其中：是第1主成分，是第2主成分，…,是第6主成分。（4）选择个主成分，计算综合评价值。eq\o\ac(○,1)计算特征值的信息贡献率和累积贡献率。称为主成分的信息贡献率；并且称为主成分的累积贡献率。当接近于1时，选择前个指标变量作为个主成分，代替原来的指标变量，从而对个主成分进行综合评判。eq\o\ac(○,2)计算综合得分：其中为第个主成分的信息贡献率，根据综合得分值就可进行评价。经MATLAB计算后，以四个主成分的贡献率为权重，构建主成分综合评价模型：，把各区的四个主成分值代入上式，可以得到各区的排名和综合评价结果，如下表所示：不合理度排名CEFDAB综合评价值1.19311.09320.3389-0.1788-0.9783-1.4681表4由上表可知，C、E两区交巡警平台设置数目较不合理。对交巡警平台设置地理位置的合理性分析：本模型沿用问题一第三问的模型，分别对六个区进行合理性评价：计算各交巡警服务平台工作量方差及出警时间：其中：：区交巡警服务平台工作量方差；：区交巡警服务平台编号；：区非交巡警服务平台编号；：编号的交巡警服务平台所管辖节点的案发率；：区每个节点（路口）平均案发率；：区交巡警服务平台个数。对于：①当对应唯一一个点时，则；②当对应多个点时，；③当没有对应的点时，采取就近原则，。经MATLAB计算后得：区号评价指标ABCDEF方差8.431415.456732.1398512.76515.198330.0042工作量均值5.88758.311.011767.53333310.077788.811111出警时间>3的点数6443122934表5由上表可知C、E、F、B四区工作量较大，C、E、F三个区有较多的路口节点无法满足出警时间尽可能小于三分钟的条件。5．2．Ⅱ.2解决方案针对工作量较大、出警时间较长、交巡警服务平台设置较少等问题，参照问题一第二问的解决方案，在某些路口节点增设交巡警服务平台，使得各区每个交巡警服务平台工作量均值减小、方差减小，同时保证出警时间减短。沿用问题一第二问的模型，现得出解决方案如下表：各区编号

项目ABCDEF增加点数5312365改进前工作量方差8.4315.4632.1412.7715.2030.00改进后工作量方差2.877.435411.1998.114.73337.16改进前工作量均值5.898.3011.017.5310.088.81改进后工作量均值4.985.5083334.4666675.656.2166675.783333改进前大于3点数6.004.0043.0012.0029.0034.00改进后大于3点数2212597表6新增加的点对应原图中的编号为：改进后前后对比折线图如下：其中横坐标分别代表六个区图4图5图6Ⅱ围堵方案5．2．Ⅱ.1模型的建立沿用第一问的程序，求出嫌疑犯在三分钟之后仍在区内。设编号为的交巡警平台到达围捕节点的时间为，凶犯到达点的时间为，建立0-1变量，，当调动点交巡警时，反之为0；当嫌疑犯到达点时，反之为0可建立模型如下：经……（暂未求出）6．模型的评价与推广评价：评价交巡警平台设置是否合理主要有“质”“量”两个指标。“质”又分为两个小指标，其一为出警时间，其二为工作量均衡；“量”即交巡警服务平台的数量是否与该区面积人口相适宜。我们在建立模型的过程中全面参考了题目中给出的数据，以最短路径Floyd算法解决出警时间问题，以方差来衡量各区工作量是否均衡的问题，用主层次分析交巡警服务平台数目设置是否合理。针对出警时间问题，由图中给出的路口节点数据可以得知总有一些节点距离其他点较远，无法实现出警时间小于三分钟的要求；使用方差衡量工作量不均衡时，每当增加一个交巡警服务平台就会使方差减小一定数值，在计算过程中我们没有排除数量增加导致的方差减小的干扰因素；在对整个市区交巡警服务平台设置是否合理进行评价时，我们参考了改进后的区工作量均值，而没有给出实际应该参考的理论值。推广：本题解决了城市规划中关于站点设置的问题，在实际生活中，该模型也可以用于其他规划问题。7．参考文献[1]/question/317403473.html[2]数学建模算法与应用，司守奎、孙玺菁，国防工业出版社8.附录8．1附录清单8．1．1求解问题一第一问的MATLAB程序；8．1．2求解问题一第二问的MATLAB程序；8．1．2求解问题一第二问的MATLAB程序；8．1．3求解问题二第一问的MATLAB程序；8．1．4求解问题二第二问的MATLAB程序；8．2附录正文8．2．1.1求解问题一第一问的MATLAB程序；clc,clearx=[413.0000403.0000383.5000381.0000339.0000335.0000317.0000334.5000333.0000...282.0000247.0000219.0000225.0000280.0000290.0000337.0000415.0000432.0000...418.0000444.0000251.0000234.0000225.0000212.0000227.0000256.0000250.5000...243.0000246.0000314.0000315.0000326.0000327.0000328.0000336.0000336.0000...331.0000371.0000371.0000388.5000411.0000419.0000411.0000394.0000342.0000...342.0000325.0000315.0000342.0000345.0000348.5000351.0000348.0000370.0000...371.0000354.0000363.0000357.0000351.0000369.0000335.0000381.0000391.0000...392.0000395.0000398.0000401.0000405.0000410.0000408.0000415.0000418.0000...422.0000418.5000405.5000405.0000409.0000417.0000420.0000424.0000438.0000...438.5000434.0000438.0000440.0000447.0000448.0000444.5000441.0000440.5000...445.0000444.0000];y=[359.0000343.0000351.0000377.5000376.0000383.0000362.0000353.5000342.0000...325.0000301.0000316.0000270.0000292.0000335.0000328.0000335.0000371.0000...374.0000394.0000277.0000271.0000265.0000290.0000300.0000301.0000306.0000...328.0000337.0000367.0000351.0000355.0000350.0000342.5000339.0000334.0000...335.0000330.0000333.0000330.5000327.5000344.0000343.0000346.0000342.0000...348.0000372.0000374.0000372.0000382.0000380.5000377.0000369.0000363.0000...353.0000374.0000382.5000387.0000382.0000388.0000395.0000381.0000375.0000...366.0000361.0000362.0000359.0000360.0000355.0000350.0000351.0000347.0000...354.0000356.0000364.5000368.0000370.0000364.0000370.0000372.0000368.0000...373.0000376.0000385.0000392.0000392.0000381.0000383.0000385.0000381.5000...380.0000360.0000];n=[1.7 2.1 2.2 1.7 2.1 2.5 2.4 2.4 2.1 1.6 2.6 2.4 2.2 2.5 2.1 2.6 2.5 1.9 1.8 1.9 1.4 1.4 2.4 1.1 1.6 1.2 0.8 1.3 1.4 2.1 1.6 1.5 1.4 1.7 1.4 1.1 0.1 1.2 1.4 1.7 1.4 1.4 1.7 1.1 1.4 1.2 1.6 1.4 1.2 1.1 0.8 0.6 1.4 0.9 1 0.5 0.8 1.1 0.9 0.7 0.6 1.2 1.4 0.8 0.7 0.8 0.8 0.9 1.1 0.9 1.1 0.8 0.9 1.1 0.8 1.1 0.8 0.8 0.8 0.8 1.4 1.1 0.9 1 1.2 1.4 1.1 0.9 1.4 0.9 0.9 0.8];fori=1:92forj=1:92d(i,j)=sqrt((x(j)-x(i))^2+(y(j)-y(i))^2);endendn=xlsread('B题附件2_全市六区交通网路和平台设置的数据表.xls','全市交通路口的路线','A2:A144');m=xlsread('B题附件2_全市六区交通网路和平台设置的数据表.xls','全市交通路口的路线','B2:B144');s=zeros(92,92);fori=1:143ifm(i)<=92s(n(i),m(i))=1;endendn=92;a=s.*d;a=a+a';fori=1:nforj=1:nifi~=j&a(i,j)==0a(i,j)=inf;endendendfork=1:nfori=1:nforj=1:nifa(i,j)>a(i,k)+a(k,j)a(i,j)=a(i,k)+a(k,j);endendendendt=a/10;forj=21:92[c(j-20),I(j-20)]=min(t(1:20,j));%I显示的管理点，c是相应的时间；end8．2．1.2求解问题一第二问的MATLAB、LINGO程序MATLAB:clc,clearx=[413.0000403.0000383.5000381.0000339.0000335.0000317.0000334.5000333.0000...282.0000247.0000219.0000225.0000280.0000290.0000337.0000415.0000432.0000...418.0000444.0000251.0000234.0000225.0000212.0000227.0000256.0000250.5000...243.0000246.0000314.0000315.0000326.0000327.0000328.0000336.0000336.0000...331.0000371.0000371.0000388.5000411.0000419.0000411.0000394.0000342.0000...342.0000325.0000315.0000342.0000345.0000348.5000351.0000348.0000370.0000...371.0000354.0000363.0000357.0000351.0000369.0000335.0000381.0000391.0000...392.0000395.0000398.0000401.0000405.0000410.0000408.0000415.0000418.0000...422.0000418.5000405.5000405.0000409.0000417.0000420.0000424.0000438.0000...438.5000434.0000438.0000440.0000447.0000448.0000444.5000441.0000440.5000...445.0000444.0000];y=[359.0000343.0000351.0000377.5000376.0000383.0000362.0000353.5000342.0000...325.0000301.0000316.0000270.0000292.0000335.0000328.0000335.0000371.0000...374.0000394.0000277.0000271.0000265.0000290.0000300.0000301.0000306.0000...328.0000337.0000367.0000351.0000355.0000350.0000342.5000339.0000334.0000...335.0000330.0000333.0000330.5000327.5000344.0000343.0000346.0000342.0000...348.0000372.0000374.0000372.0000382.0000380.5000377.0000369.0000363.0000...353.0000374.0000382.5000387.0000382.0000388.0000395.0000381.0000375.0000...366.0000361.0000362.0000359.0000360.0000355.0000350.0000351.0000347.0000...354.0000356.0000364.5000368.0000370.0000364.0000370.0000372.0000368.0000...373.0000376.0000385.0000392.0000392.0000381.0000383.0000385.0000381.5000...380.0000360.0000];fori=1:92forj=1:92d(i,j)=sqrt((x(j)-x(i))^2+(y(j)-y(i))^2);endendn=xlsread('B题附件2_全市六区交通网路和平台设置的数据表.xls','全市交通路口的路线','A2:A144');m=xlsread('B题附件2_全市六区交通网路和平台设置的数据表.xls','全市交通路口的路线','B2:B144');s=zeros(92,92);fori=1:143ifm(i)<=92s(n(i),m(i))=1;endendn=92;a=s.*d;a=a+a';fori=1:nforj=1:nifi~=j&a(i,j)==0a(i,j)=inf;endendendfork=1:nfori=1:nforj=1:nifa(i,j)>a(i,k)+a(k,j)a(i,j)=a(i,k)+a(k,j);endendendendt=a/10;zz=[12141621222324282930384862];fori=1:20forj=1:13z(i,j)=t(i,zz(j))endendzzLINGO：模型一：model:title封锁路径最快问题;sets:plot/1..13/:a;plot2/1..20/:b;links(plot2,plot):c,d;endsetsdata:c=22.236216.02859.286819.293421.096222.501822.893219.001219.515812.08345.880911.85014.885220.463914.12977.388117.394719.197520.603021.121017.228917.743610.31123.982210.30956.035118.352312.76726.025616.032217.835019.240519.009315.117315.63198.19966.09388.19794.393421.997415.00858.266918.273520.076321.481822.654416.226915.53538.10304.86107.39590.350017.628212.96966.228016.234617.749519.155118.285211.306910.61533.18299.42112.47585.255117.658813.00026.258616.265217.780119.185618.315811.337510.64593.21359.45172.50645.337314.914910.90124.159614.166215.036316.441815.57208.57028.01550.58317.35271.29027.991714.09259.43392.692312.698914.213815.619414.749610.228010.49323.06085.88543.09958.677313.01078.27421.532511.539213.132014.537613.66789.775710.72443.49234.72574.19949.33677.586612.77576.95679.51077.70799.11358.243614.194915.14357.911410.14988.618614.76083.79148.337311.39505.07233.26964.67513.805318.633219.581912.349814.588213.056919.1992011.950314.54338.68536.88256.47703.591621.781422.730115.498017.736416.205122.34745.97705.973312.71492.70830.90550.50002.385422.808323.757016.524916.120817.232021.331811.950306.74173.26505.06776.47338.358718.049918.916711.484310.147512.191415.358517.029613.29816.556416.563017.150918.556517.68674.75185.70054.40159.74965.108611.810114.54336.7417010.006611.809413.214915.100311.308312.17504.74273.40595.44988.616921.892114.90328.161618.168219.971021.376522.549218.657119.523912.09154.755712.79867.820524.247218.514511.772821.779523.582224.987824.904221.012221.526814.09458.366913.69936.734422.546516.961510.219820.226422.029223.434823.203619.311519.826212.39387.639311.99865.033726.945821.213114.471424.478126.280927.686427.602923.010822.319214.886911.065614.17986.4489;enddata@for(links(i,j):@bin(d(i,j)));!对d（i,j）进行0，1分布;min=@sum(links(i,j):c(i,j)*d(i,j));@for(links(i,j):@bin(d(i,j)));!对d（i,j）进行0，1分布;@sum(links(i,j):d(i,j))=13;!必须有十三个队派出去;@for(plot(j):@sum(plot2(i):d(i,j))>=1);!每个地方的警队必须大于等于1;@for(plot2(i):@sum(plot(j):d(i,j))<=1);!每个警队的警力不足只能派出到一个地方;模型二：model:title封锁路径最快问题;sets:plot/1..13/:a;plot2/1..20/:b;links(plot2,plot):c,d;endsetsdata:c=22.236216.02859.286819.293421.096222.501822.893219.001219.515812.08345.880911.85014.885220.463914.12977.388117.394719.197520.603021.121017.228917.743610.31123.982210.30956.035118.352312.76726.025616.032217.835019.240519.009315.117315.63198.19966.09388.19794.393421.997415.00858.266918.273520.076321.481822.654416.226915.53538.10304.86107.39590.350017.628212.96966.228016.234617.749519.155118.285211.306910.61533.18299.42112.47585.255117.658813.00026.258616.265217.780119.185618.315811.337510.64593.21359.45172.50645.337314.914910.90124.159614.166215.036316.441815.57208.57028.01550.58317.35271.29027.991714.09259.43392.692312.698914.213815.619414.749610.228010.49323.06085.88543.09958.677313.01078.27421.532511.539213.132014.537613.66789.775710.72443.49234.72574.19949.33677.586612.77576.95679.51077.70799.11358.243614.194915.14357.911410.14988.618614.76083.79148.337311.39505.07233.26964.67513.805318.633219.581912.349814.588213.056919.1992011.950314.54338.68536.88256.47703.591621.781422.730115.498017.736416.205122.34745.97705.973312.71492.70830.90550.50002.385422.808323.757016.524916.120817.232021.331811.950306.74173.26505.06776.47338.358718.049918.916711.484310.147512.191415.358517.029613.29816.556416.563017.150918.556517.68674.75185.70054.40159.74965.108611.810114.54336.7417010.006611.809413.214915.100311.308312.17504.74273.40595.44988.616921.892114.90328.161618.168219.971021.376522.549218.657119.523912.09154.755712.79867.820524.247218.514511.772821.779523.582224.987824.904221.012221.526814.09458.366913.69936.734422.546516.961510.219820.226422.029223.434823.203619.311519.826212.39387.639311.99865.033726.945821.213114.471424.478126.280927.686427.602923.010822.319214.886911.065614.17986.4489;enddata@for(links(i,j):@bin(d(i,j)));!对d（i,j）进行0，1分布;min=v;@for(links(i,j):v>=c(i,j)*d(i,j));@for(links(i,j):@bin(d(i,j)));!对d（i,j）进行0，1分布;@sum(links(i,j):d(i,j))=13;!必须有十三个队派出去;@for(plot(j):@sum(plot2(i):d(i,j))>=1);!每个地方的警队必须大于等于1;@for(plot2(i):@sum(plot(j):d(i,j))<=1);!每个警队的警力不足只能派出到一个地方;8．2．1.3求解问题一第三问的MATLAB程序clc,clearx=xlsread('mintime.xls','sheet1','A1:CN92');r=xlsread('B题附件2_全市六区交通网路和平台设置的数据表.xls','全市交通路口节点数据','E2:E93');Z=zeros(92,92);h=zeros();Y=zeros(92,92);D=zeros(92,92);fori=1:92%挑出1-20，m,n行的数据ifi<=20forj=1:92D(i,j)=x(i,j);endendifi==39|i==48|i==58|i==67|i==87forj=1:92D(i,j)=x(i,j);endendendfori=1:92forj=1:92ifD(i,j)==0D(i,j)=10000;endendendforj=1:92[y(j),I(j)]=min(D(:,j));endforj=1:92Y(I(j),j)=y(j);endforj=1:92ifj<=20|j==39|j==48|j==58|j==67|j==87%剔除平台相互管理管理Y(:,j)=0endendfori=1:92forj=1:92ifY(i,j)~=0Y(i,j)=1;endifi==39|i==48|i==58|i==67|i==87|i<=20Y(i,i)=1;endendendforj=1:92Y(:,j)=r(j)*Y(:,j);ends=sum(Y,2)s=nonzeros(s)d=var(s);forj=1:92[h(j),H(j)]=max(Y(:,j));endfori=1:92z(i)=x(i,H(i))end8．2．1求解问题二第一问的MATLAB程序clc,clearx=xlsread('mintimeB.xls','sheet1','A1:BU73')r=xlsread('B题附件2_全市六区交通网路和平台设置的数据表.xls','全市交通路口节点数据','E94:E166')n=length(x);g=8;h=zeros();Y=zeros(n,n);D=zeros(n,n);fori=1:n%挑出1-g行的数据ifi<=gforj=1:nD(i,j)=x(i,j);endendendfori=1:nforj=1:nifD(i,j)==0D(i,j)=10000;endendendforj=1:n%找出里j点最近的平台作为管理平台。[y(j),I(j)]=min(D(:,j));endforj=1:nY(I(j),j)=y(j);endforj=1:nifj<=g%剔除平台相互管理管理Y(:,j)=0endendfori=1:n%把Y矩阵变成0-1矩阵，用于发案率求和forj=1:nifY(i,j)~=0Y(i,j)=1;endifi<=gY(i,i)=1;endendendforj=1:nY(:,j)=r(j)*Y(:,j);ends=sum(Y,2)%将i点所管理的所有j点发案率累加求和s=nonzeros(s)%去除0点d=var(s);%对工作量（发案率）求方差forj=1:n[h(j),H(j)]=max(Y(:,j));%输出i-j对应点和其j发案率endfori=1:nz(i)=x(i,H(i))%对出警时间end（改进之前求解方差的程序）clc,clearx=xlsread('mintimeB.xls','sheet1','A1:BU73')r=xlsread('B题附件2_全市六区交通网路和平台设置的数据表.xls','全市交通路口节点数据','E94:E166')n=length(x);g=8;h=zeros();Y=zeros(n,n);D=zeros(n,n);fori=1:n%挑出1-g，m,n行的数据ifi<=gforj=1:nD(i,j)=x(i,j);endendifi==123-92|i==145-92|i==155-92forj=1:nD(i,j)=x(i,j);endendendfori=1:nforj=1:nifD(i,j)==0D(i,j)=10000;endendendforj=1:n[y(j),I(j)]=min(D(:,j));endforj=1:nY(I(j),j)=y(j);endforj=1:nifj<=g|j==123-92|j==145-92|j==155-92%剔除平台相互管理管理Y(:,j)=0;endendfori=1:nforj=1:nifY(i,j)~=0Y(i,j)=1;endifi==123-92|i==124-92|i==155-92|i<=gY(i,i)=1;endendendforj=1:nY(:,j)=r(j)*Y(:,j);ends=sum(Y,2);s=nonzeros(s);d=var(s);forj=1:n[h(j),H(j)]=max(Y(:,j));endfori=1:nz(i)=x(i,H(i));end（改进之后求解方差的程序）基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造\t"_bla