交巡警服务平台的设置与调度

1、交巡警服务平台的设置与调度【摘要】警察是现代社会中不可或缺的社会角色，肩负着执法、治安与服务群众等重要 职能。为了更好地履行这些职能，交巡警服务平台要合理地分布在城市的各个地区，这 样不仅可以及时响应出警到达案发现场，在遇到突发事件时也可以通过联合调度高效地 行动起来。该论文就交巡警服务平台的设置与调度等实际问题，针对所提出的 5 个问题分别给 出具体的解决方案并给出结果：对于问题 1 要给 A 区的每个服务平台分配管辖范围，即分配其管辖的节点。我们根 据“就近原则”来分配管辖的节点，保证尽量在 3 分钟内有交巡警到达事发地。对此， 借助MATLA编程采用“ Floyd最短路径算法”确定距离每

2、个节点最近的服务平台，从而 得到每个服务平台的管辖范围。对于问题 2 的合理的调度方案的确定，我们在“快速封锁”的原则下，通过调度警 力使得A区在最短时间内被全封锁。20个服务平台对13个路口进行全封锁，而且每个服 务平台最多封锁一个路口，这可划归于一个 0-1规划问题，因此可用LINGO编程求得各 种可选调度方案中 13 个路口封锁时间的最大值取值最小时的调度情况。对于问题 3 增加平台的个数与位置的确定，我们的目的是使各个服务平台的工作量 达到均衡状态而且出警时间过长的问题得到有效解决。为此，我们在出警时间过长的节 点或附近尝试增加新的服务平台，然后计算方差来衡量工作量的均衡程度，比较增加

3、2至 5 个服务平台时的方差，以此确定方差最小的情况为最后的可选方案。这个过程仍然 借助MATLABS序来完成，采用“模拟退火法”来确定工作量达到均衡时新增平台的个数 与位置。对于问题 4对全市服务平台设置方案的合理性的讨论，我们借助问题 1 和问题 3 的 解决方法来确定各区服务平台的管辖范围与新增服务平台的个数与位置。同时对模型进 行优化，考虑到有些服务平台的工作量过少的情况，撤消一些现有的服务平台。借助 MATLA程序，可以给出一个较合理的解决方案，即给出各个分区的服务平台的调整方案。对于问题 5 围堵方案的确定，可将全市的交通网看作一张图，各个节点看作顶点。 同时根据必要的假设：嫌疑犯

4、一直朝远离事发点 P点的方向逃跑，而且不走回路。这时， 将 P 点看作树根，嫌疑犯的可能的逃跑路线便成为一个树，有可能经过的节点便是枝和 叶。这样，就能根据图论的知识，通过 MATLA与LINGOS序，利用“追捕算法”来对各 个分支道路进行有序的封锁排查，进而求得最佳的围堵方案。关键词 ： Floyd 最短路径算法、 0-1 规划、模拟退火法、平台的设置与调度、图论、追捕 1 、问题背景现代社会中，要保证社会的良性运行，人民警察扮演着重要的角色。他们肩负着刑 事执法、治安管理、交通管理、服务群众的重要职能，在保障人民财产安全面前起着举 足轻重的作用。为了更有效地实施这些职能，需要在市区的一些交

5、通要道和重要部位设 置交巡警服务平台。这些交巡警服务平台的职能和警力配备基本相同，共同维护着社会 的治安，保证城市的安全。在交巡警服务平台的设置中我们要充分考虑到警察所担负的社会职责：刑事执法、 治安管理、交通管理、服务群众。在保证人民群众人身与财产安全的前提下还要尽量把 服务工作做好做到位，让人民群众满意，因此，交巡警服务平台的设置是关键还要考虑 到调度上的方面与快捷。而现实中由于警务资源的限制，同时考虑到社会资源的高效分配，不可能在每个交 通要到和重要部位都设置这些交巡警服务平台。因而，如何根据城市的实际情况和需求 合理的配置有限的交巡警服务平台、分配平台的管辖范围、调度警务资源是警务部门

6、面 临的一个实际而重要的问题。2、问题重述就某市设置交巡警服务平台的相关情况，建立数学模型分析研究下面的问题： 附件1中的附图1给出了该市中心城区A的交通网络和现有的20个交巡警服务平台 的设置情况示意图，相关的数据信息见附件2。全市（主城六区A，B，C,D,E，F）的交巡警服务平台的设置的具体情况与其他这些城市的情况见附件2。问题1：请为A区各交巡警服务平台分配管辖范围， 使其在所管辖的范围内出现突发 事件时，尽量能在3分钟内有交巡警（警车的时速为 60km/h）到达事发地。问题2:对于重大突发事件，需要调度 A区20个交巡警服务平台的警力资源，对进 出该区的13条交通要道实现快速全封锁。实

7、际中一个平台的警力最多封锁一个路口。请 给出遇到重大突发事件时的警力合理的调度方案。问题3:根据A区现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实 际情况，拟在该区内再增加2至5个平台，请确定需要增加平台的具体个数和位置。问题4：针对全市的具体情况，按照设置交巡警服务平台的原则和任务， 分析研究该 市现有交巡警服务平台设置方案的合理性。如果有明显不合理，请给出解决方案。问题5:如果该市地点P （第32个节点）处发生了重大刑事案件，在案发 3分钟后 接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警服务 平台警力资源的最佳围堵方案。3、基本假设现针对问题做必要的

8、假设：1、发生突发事件后，警察在接到报警时能立刻出警，时间上的消耗主要考虑到在路 上所花的时间。2、 警察出警时，在去案发现场的路上道路通畅（不发生堵车等现象），警车行驶正 常。3、交巡警服务平台的工作量主要由于出警所引起，而在现场或服务平台处理案件的 工作量忽略。4、任何事故都发生在市区道路的节点上。5、警察出警到达事发地点走最短距离，而且警车始终保持匀速行驶。&犯罪嫌疑人驾驶车辆出逃时，只能通过联通市区与外界的节点逃出市区。7、交巡警服务平台一般不能跨区管理，除非是一些重大事故。8、一般地，驾车逃跑的嫌疑犯不会走回路，而且想尽量离出事地点越远越好。9、嫌疑犯逃跑的车速是恒定的。4、符号说明

9、表1符号数学解释路口（节点）的发案率p, j两节点之间的距离i节点处交巡警服务平台的工作量i节点处的出警时间标准差警车的行驶速度嫌疑犯逃跑车速5、问题分析问题1的分析分配交巡警服务平台管辖范围的目的是当其所管辖的范围内发生突发事件时，尽量 能在3分钟之内有交巡警到达事发地。根据实际生活中的经验，人们当然也希望突发事 件发生时交巡警能够能尽快到达事发地处理事件。因而，在确定每个服务平台的管辖范 围时，应根据“就近原则”一一节点距离某个服务平台最近就属于那个服务平台管辖。 可以考虑借助“ Floyd-Warshall最短路径算法”解决这个问题。问题2的分析当发生重大突发事件时要调度全区 20个交巡

10、警服务平台的警力资源对 13条进出该 区的交通要道实现全面快速的封锁。其中时间是关键因素，我们希望封锁完成得越快越 好，而全圭寸锁的最终完成以最后一个交通要道的圭寸锁实现为标志。从而只须要求圭寸锁最 慢的那个交通要道的封锁时间最短即可。问题3的分析之所以要增加2至5个平台，是为了解决现有交巡警服务平台的工作量不均衡和有 些地方出警时间过长的实际问题。由于管辖范围的划分是依据出警时间来确定的，因此 在这里我们首先应该考虑工作量均衡问题，然后再以出警时间作为辅助来检验增加的服 务平台的合理性。而一般说来，出警时间与工作量具有一定的正相关性，因此在出警时 间较长的节点或其附近增加交巡警服务平台可以在

11、一定程度上协调解决这两个问题。由于工作量既与出警路程有关，又与案件发生频率有关，所以可用案件发生频率与 路程的积来表征工作量的大小，于是对每一个服务平台都有一个相应的工作量；同时又 可用方差来表征工作均衡的程度。因此，通过添加服务平台，可以调整管辖范围，从而 不但可以确保出警时间，也能调整工作量大小，以方差大小来衡量其均衡程度，最终可 得一个添加方案。整个过程主要采用模拟退火算法，由于此题距离较远的节点出警时间不一定能保证， 且工作量相对较大，故可从距离入手来添加服务平台。依据该算法原理，采用逐步寻优 的手段，以工作量方差为检验标准，温度下降时工作量以较大概率趋于均衡，最终将得 到一个较为合理

12、的平台设置方案。问题4的分析根据设置交巡警服务平台的原则与任务一一刑事执法、治安管理、交通管理与服务 群众，同时结合题目的要求可知平台设置的效果：第一是保证接到报警后出警时间不能 过长；第二是各个服务平台的工作量大致均衡。接下来主要依据这两个方面对服务平台 的设置进行优化。对于第一点，依据题意，不管在哪个区，出现突发事件时，尽量在相对短的时间内 有巡警到达，基于这一点，我们可以对各个区内巡警服务平台的管辖范围进行划分。这 就是说，范围的划分是为了保证出警时间不是太长。但考虑到现实中交巡警一般不能跨 区管理，并且每个区的面积、人口密度、服务平台数目也不相同，因此要对不同区分开 考虑，故不同区的出

13、警最长时间的设定标准是不同的。对于第二条，由上面划分的管辖范围，依据问题3中的方法算出各个服务平台的工作量，通过比较，我们可以直观的发现其中的不合理的地方一一一些服务平台之间的工 作量相差太大，有些节点的出警时间的相对较长，可以分如下情况进行优化：(1) 某个服务平台的工作量很大，若将其中某些节点划给其它服务平台管理则可达 相对均衡状态，但划出的这些节点归为其他服务平台则会导致出警时间过长，也即该服 务平台附近没有其他工作量较小的服务平台，两种矛盾无法调和，此时应该考虑添加一个服务平台。（2）某几个服务平台较为集中，且各自工作量相对较小，而将这些平台管辖区域外 的节点调给它时又会导致出警时间过

14、长，此时可以考虑撤消（或转移）某一个或几个服 务平台来平衡工作量。通过以上两方面的调整，采用模拟退火算法思想，通过添加和撤消服务平台，最终 工作量的方差达到最小且基本不变，而且每个节点的出警时间均不会有太大差异，这样 就得到一个相对较为合理的设置方案。问题5的分析当发生重大案件需要全市的交巡警联合围堵时，交巡警需要对嫌疑犯有可能经过的 道路提前进行围堵。将全市的交通网看作一个图，把各个节点看作顶点。我们根据假设， 嫌疑犯希望里出事地点 P节点越来越远，而且嫌疑犯所走的路线不会出线回路。这时， 我们将P点当作根，嫌疑犯可能的逃跑路线构成一个树。交巡警便可以根据逃跑路线所 构成的树对嫌疑犯进行围堵

15、。根据图论的知识，借助MATLAB! LINGO编程便可以得到一个围堵方案。 6模型的建立与分析问题1的求解6.1.1模型的建立与求解根据“就近原则”的要求只需要确定距离每个节点最近的交巡警服务平台是哪个， 从而就可以确定每个交巡警服务平台的管辖范围。这里我们用到“Floyd-Warshall最短路径算法”。这个算法通过不停地加入中间点进行矩阵的迭代，从而解决有向图中每对顶点的最短路径问题，当我们把每个节点看作顶点时，我们便可以求得每个服务平台的 管辖范围一一管辖的节点。由MATLABS程（见附录1），利用Floyd-Warshall最短路径算法可以得到 A区20 个服务平台的管辖范围，即每个

16、服务平台所管辖的节点编号。A区服务平台的工作量见表 2,结果见表3:表2A区交巡警平台编号A区交巡警平台的位置 编号A区交巡警平台的工作量A11A22A33A1818A1919A2020表3A区交巡警平台 编号A区交巡警平台位置编号每个交巡警平台所管辖节点 编号A111、67、68、69、71、73、74、75、 76、 78A222、39、40、43、44、70、72A333、 54、 55、 65、 66A444、57、60、62、63、64A555、49、50、51、52、53、56、5& 59A666A777、30、32、47、48、61A888、 33、 46A999、 31、 34

17、、 35、 45A101010A111111、26、27A121212、25A131313、21、22、23、24A141414A151515、28、29A161616、36、37、38A171717、41、42A181818、80、81、82、83A191919、77、79A202020、84、85、86、87、88、89、90、 91、 92A区各个交巡警平台的相应的管辖范围如图1所示（MATLA程序见附录2）图1 （不同管辖范围用不同的颜色标识）6.1.2模型的分析从结果中可以很直观地看到不同交巡警平台所管辖的周围的节点个数是不同的，从而管辖范围不同，他们的工作量当然也不相同。例如平台A

18、20要管辖10个节点，而平台A6仅管辖1个节点，这显然是不合理的一一利用最短路径的思想划分的管辖范围会使得 服务平台的工作量的不均衡。问题2的求解6.2.1模型的建立与求解当遇到突发事件需要对 A区进行全封锁时，合理的调度方案应该是使出警时间最长 的节点处的出警时间达到最小值，其余的节点处出警时间尽量减小。又由于警车匀速行 驶，所以路程的长短反映出出警时间的长短，故可用距离来代替时间考虑。由于20个服务平台的交巡警封锁13个路口，每个路口要有一个平台的交巡警负责，每个交巡警平台最多只能封锁一个路口，可以考虑用0-1规划来解决。目标函数为13个路口到各个服务平台的距离的集合，求解目标函数就是对这

19、些距离进行搜索，找出达到 13个一一对应时的最小距离，也就是说对所有可选的调度方案，找出每个方案的最大距 离，然后从这些最大距离中取出最小值。对于取出的方案，由于前12个对应距离均比这个距离小，所以这个对应就表示完成封锁所需时间的最小值。借助计算机编程可得结果。其LINGO程序见附录3,运行结果如下：以A7服务平台封锁29号出口节点即可找到符合题意的一种封锁方案，并且这个方 案中的最大距离是所有可行方案中最小的，所以说这个结果是最佳的。而其余的封锁时 间均比完成该封锁用时短，故只用保证第 29号出口节点由A7服务平台负责封锁，其余 的节点处可有不止一种选择，这里我们给出一种相对较佳的方案作参考

20、，见表4：表4路 口 号12141621222324282930384862平台号A10A16A9A14A11A13A12A15A7A6A2A5A4622模型分析LINGO程序的最终结果显示只需要将 29号路口由A7来封锁，其余路口在保证选择与 他们距离短于A7-29之间的距离的平台的情况下，尽量选较短距离。即便如此，调度方 案也不唯一，这里仅给出一个相对较佳的可行方案。问题3的求解6.3.1模型建立与求解为解决服务平台工作量的不均衡和有些地方的出警时间过长的实际问题，可以考虑 添加新的服务平台来分摊工作量较大的服务平台的工作量，又由于工作量大小和出警时 间在一定程度上有正相关性，故一般也可同

21、时减少出警时间过长节点的出警时间。每个服务平台都有一定的工作量，我们定义工作量函数其表示标号为i的服务平台的工作量，nj为其管辖的节点处的发案率。节点i处的出警时间为由于v不变，ti仅与dj成正相关关系，距离长的地方出警时间相应会长。dj为节点与管辖其的服务平台之间的距离。各个服务平台的工作量的均衡程度由各个服务平台的工作量之间的方差s来决定。当我们计划增加2至5个服务平台时，我们最终要使得方差s达到最小值时的工作量 达到最均衡的状态。这时便可以确定所增加的服务平台的个数与位置。对此，我们运用“模拟退火算法2 ”对其进行求解。该算法是寻求各个状态点之间的 一个均衡的状态。而我们正是要通过增加服

22、务平台的个数使得更个服务平台的工作量达 到基本均衡，同时使得出警时间过长的节点的出警时间得到有效缩短。MATLAB?序见附录4，运行结果表5:表5新增服务平台个数新增服务平台所在节点标号521、 24、 28、 38、 92521、 24、 28、 39、 92新增平台如图2:图2 （其中青色加号表示新增平台）6.3.2模型分析根据程序运行的结果，可知在 38或39号节点处增加新的服务平台所得结果差异不 大，因此都可以当作供选方案。由于退火模型存在普遍的缺陷，根据所选的“退火系数” t的不同，得到的结果可以有一定的差异，但对问题没有本质上的影响。如t取不同值时, 92这个节点也可不添加服务平台

23、，但也符合时间限制，对工作量均衡也无显着影响，只 是上面的添加方案可以更加有效，工作量也相对较小，即这是一个较佳方案，故选之。问题4的求解6.4.1模型的建立与求解对现有交巡警服务平台设置的合理性的讨论主要有两个指标来确定一一各服务平台 的工作量的差异与各个节点的出警时间的长短。现实生活中不同区的警力管辖范围仅限于该区，只有发生特殊事件时才有跨区合作 的机会，因此我们将全市六区分开考虑，每个区的服务平台的管辖范围仅限于该区，那 么问题的解决思路与问题3有很多类似的地方。问题的解决过程如下：（1）计算各区的服务平台的工作量，及它们的方差大小。计算各个节点处的出警时 间。（2）根据问题三的方法对每

24、个区的服务平台进行优化。（3）对于服务平台集中且各个服务平台的工作量都不大的情况考虑撤销一个服务平 台节约警力资源。根据求解以上问题的思路，我们把全市划分为6各区，各个区由Floyd-Warshall最短路径算法求得各个服务平台的管辖范围，先通过直观地观察撤销一些明显的工作量 过少的的平台，然后由模拟退火算法通过增加服务平台的个数来均衡各个平台的工作量。由MATLA编程运行得到各个服务平台的工作量如表 6,得到结果如表7：表6交巡警平台编号交巡警平台的位置标号交巡警平台的工作量A11A22A33F9483F10484F11485表7区号新增平台所在节点的标号撤销平台所在节点的标号A21、 24

25、、 28、 38 或 39、 92B：120、14797、99C205、240、261、281、301、313166、169、179、178D333、362322、 325E：391、409、417、452372、 376、 377F493、515484增加与撤销平台如图3所示：（其中蓝色表示原服务平台，青色加号表添加新平台，红色叉号表撤消平台）图3642模型的分析原方案的设计从时间和工作量上来讲绝大部分是合理的，只是有些地方设置不太科 学。事实上，考虑时间因素和工作量因素占不同比重时优化方案也会不同。此处先考虑 的是时间因素，对各个区域的服务平台管辖范围进行划分，然后采用问题3中方法给出工作

26、量并得出撤消和添加的服务平台，然后再以时间限制进行优化检验。当然，采取不 同的均衡标准与时间限制也可以有不同的优化方案。问题5的求解6.5.1模型的建立与求解当发生重大案件需要全市的交巡警联合围堵时，交巡警需要对嫌疑犯有可能经过的 道路提前进行围堵。将全市的交通网看作一个图，把各个节点看作顶点。根据假设，嫌 疑犯希望离出事地点 P节点越远越好，而且嫌疑犯所走的路线不会出现回路。因此，与 同一个顶点相连的任意两个顶点不会连通（即无圈），这样可以将P点当作树根，而嫌疑 犯可能的逃跑路线则构成一个树，这个树采用广度优先搜索BFS的方法得到的（附录5）下面我们用图论的知识对这个围堵问题进行求解:假设嫌

27、疑犯可能的逃跑路线所构成的树为 T,我们的算法如下：（1）把T的全部叶删除，得到子树 ;若Ti是一个顶点组成的，算法终止； 否则转（2）。（2）用T替换T ,转（1）。（3）反复执行（1）和（2）。对于示意图如下图3的树，我们先从v2、v3、v4、v6、v7、v10、v12、v13点出发，分 别向 W、v5、v8、v9 围堵，然后删除 v2、v3、v4、v6、v7、v10、v11 v12、v13 点。接着从 W、*、*、v9点出发，分别向5、U2、U3围堵，依据这样的思路进行下去便可以对嫌疑犯 进行围堵。图3由于嫌疑犯逃跑路径具有一定的随机性，而且其车速又是恒定的。从该假设出发， 利用上述算法

28、，不妨定义度数相同的顶点构成的集合为圈，这样可分不同的圈对其进行 围堵。根据嫌疑犯车速来确定在哪个圈进行围堵， 然后可用上述算法不断去掉“叶子”（外 围顶点）最后把嫌疑犯包围。确定在给定圈内围堵疑犯可否成功方法如下：用问题 2的 模型计算出警察最快封锁该圈的时间（t1 ），比较嫌疑犯在恒定速度下到达该圈的时间（t2 ），以此来判断：当t1 t2时，围堵成功；t1 t2时，围堵失败。由MATLAB程（见附录6）和LINGQ编程（见附录7）得到第三圈下的一个结果： MATLAB!行结果（疑犯逃跑路程）: arealock（D,a,b）ans =LINGO运行结果（围堵过程中的最大距离）：t1Loc

29、al optimal soluti on found. Objective value: Objective bound: In feasibilities:Exte nded solver steps: Total solver iteratio ns: 21.634350.3606 ；60t2时，选择该围堵方案； 模型的分析t242545718.46233当t16.5.2t1ut2时，不选择此方案。由于车速的不确定性，并且嫌疑犯逃跑路径也有一定的随机性，因此围堵方案也有 不确定性。而且在假设中忽略了一些联通的道路，这也对结果造成了一定程度的影响。7、模型评价基于本问题有大量数据的前提，而且

30、问题比较复杂，我们主要借助于计算机编程的 思想，运用各种算法（如Floyd-Warshall最短路径算法、模拟退火算法）解决各个实际 问题，最后得到一些合理的答案。当寻求遇到突发事件时各区警力资源的合理调度问题时，我们要求尽量快速地实现 封锁与围堵，这时根据图论的知识，采用BFS算法并结合假设可得到围堵方案设计的一些准则，依据这些准则可选出围堵方案。但由于不确定的因素太多，无法保证结果的唯一性与准确性，因此只给出一种结果，它还是能够在一定程度上说明问题的。8、参考文献1Thomas Charles , Ronald Clifford Stein, Introduction to Algorit

31、hms M, American:China Machine Press,2009,386-3882 Cabbage,Dan,Zfj3000 等. 模拟退火算法 EB/OL.3 王树禾,图论M,北京：科学出版社，2004, ,45-479、附录附录1Floyed-Warshall 最短路径算法MATLAB?序Floyed-Warshall 算法程序代码如下：% 这个算法采用的是 floyd-warshall 最短路径方法求解% W表示初始时每个点与点之间联通的距离的矩阵化% i 表示起点， j 表示终点%k=0 时% D(i,j) 表示中间点集合为 1,.,k 时的 i 点到 j 点的最短路径，

32、其中 D(i,j)=W(i,j)% path 表示最短路的路径矩阵%function D=floydwarshall(W)n=size(W,1);D=W;path=zeros(n,n);a b=ind2sub(n,n,find(W=inf);for k=1:npath(a(k),b(k)=b(k);endfor k=1:nfor i=1:nfor j=1:nif D(i,j)D(i,k)+D(k,j)D(i,j)=D(i,k)+D(k,j);path(i,j)=path(i,k);endendendend 生成最短路径的函数程序代码如下：function D=setshortestpathA=

33、xlsread(fujian,1);B=xlsread(fujian,2);W=inf(92,92);l=triu(W,1);u=tril(W,-1);tmp=zeros(92,92);W=l+u+tmp; nB=size(B,1);for i=1:nBif all( B(i,1)=92, B(i,2)30% i% position% endhold ona b c=ind2sub(3 3 3,position+6);plot( A(i,2),A(i,3),*,color, x(a,b,c) y(a,b,c) z(a,b,c) );endn=size(B,1);for i=1:nif all(

34、 B(i,1)=92, B(i,2)=92)hold online( A(B(i,1),2),A(B(i,2),2), A(B(i,1),3),A(B(i,2),3),color,k);endendxlabel(A 区交通网络与平台设置的分配管辖范围初步示意图 ); 附录3求解问题2的程序代码LINGO 程序代码如下：model:sets:cordon/c1.c13/;police/p1.p20/;forme/1.13/:f;links(police,cordon):distance,volume;endsets min=smax(f(1),f(2),f(3),f(4),f(5),f(6),f

35、(7),f(8),f(9),f(10),f(11),f(12),f(1 3);for(cordon(J):sum(police(I):volume(I,J)=1;); for(police(I):sum(cordon(J):volume(I,J)=0;);for(links:bin(volume);); for(cordon(J):f(J)=sum(police(I):volume(I,J)*distance(I,J);); data: enddata end附录 4模拟退火算法MATLAB?序代码如下：function sa2addpoliceA=xlsread(fujian);rand(s

36、eed,sum(clock);D=setshortestpath;pposition=1:20;for m=1:5% 计算管辖区域的归属dist=;p2p=;p1p=;i=21;while i=92if isempty(find(i=pposition)tmp1,tmp2=min( D(i,pposition) );dist=dist tmp1;p2p=p2p tmp2;p1p=p1p i;i=i+1;elsei=i+1;continue;endenddist,tmp3=sort(dist);p2p=p2p(tmp3);p1p=p1p(tmp3);% 计算初始方差worksum=zeros(1

37、,20+m-1);i=1;while i=92if isempty(find(i=pposition)worksum( find( (p2p( find(i=p1p) ) pposition)=1) )= .worksum( find( (p2p( find(i=p1p) ) = pposition)=1) + dist(find(i=p1p)*A(i,5);i=i+1;elsei=i+1; continue;end end worksum variance=var(worksum);%SA 模型t=1; fo=variance;j=1; loop=0;%SA 循环开始 while j=72 if loop=20 break;end%