CUMCM 2011B 期中 作业

1、数学建模 期中实验电子信息工程全鹏举 20131060282杨杰波20131060178杨 济20131060181交巡警服务平台设置的优化模型摘要城区交巡警服务平台的合理配置，对有效提高交巡警的执法效率有着重要意义，根据该市道路交巡警服务平台分布所面临的实际问题，本文就如何合理地分配各平台的管辖范围、设置交巡警服务平台、调度警务资源等实际课题展开研究。具体方法与结果如下：对于问题一，首先运用Floyd算法，通过MATLAB计算了A区路口节点任意两点之间的最短距离，得到距离矩阵D。然后通过比较，筛选出每个路口对应的最近的交巡警平台，整合得到每个平台负责的区域范围（见表5-1）。发生重大突发事件

2、时，通过分析附件2中的表格数据,可知除去交巡警服务平台所在的路口节点12、14、16之外，只剩下17个平台和10个要封锁的路口节点，问题即可转化为找17个交巡警服务平台分别到10个要封锁的路口节点的最短距离。统计它们在附录二中的最短距离，算出到每个路口的最短时间分别为0.35、10.615、2.506、0.5831、14.214、1.434、8.244、4.575、2.708、4.7518min，基本符合实际情况，该方案即为交巡警服务平台合理分配警力的调度方案。针对A区交巡警平台设置的不合理性，通过比较交巡警平台的工作量和发案率，分别在编号为1、2、5、7、20的交巡警平台处增加平台，可以使得

3、平台工作量的不均衡性和部分地方出警时间过长的情况都得到显著改善。对于问题二，运用Floyd算法划分出了每个交巡警平台负责的区域，结合各个路口节点的发案率，计算出每个交巡警平台的工作量(包括负责节点的数目、发案率之和、发案率×距离之积三种情况)。运用MATLAB编程分别得到交巡警平台与负责节点的数目、发案率之和、发案率×距离之和的条状图，经比较发现，部分平台的工作量显著地大于其他平台，因此该现有设置方案不合理。通过在编号为94、96、167、171、178、179、180、379、382、383、475、476、477、478的路口节点处设置交巡警平台，可以显著地改善这一情况

4、。若P点发生紧急重要事件，可以以事发地点P为中心，采取由内到外，有密到疏的搜捕包围圈的方法，得到了一个合理的调配方案。关键词：最短路问题 Floyd算法 MATLAB编程 交警服务平台设置一、问题简述为了更有效地贯彻实施警察的执法、治安等职能，需要在市区的一些交通要道和重要部位设置交巡警服务平台，每个交巡警服务平台的职能和警力配备基本相同。由于警务资源是有限的，如何根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题。现在就某市设置交巡警服务平台的相关情况，需要我们建立数学模型分析研究下面的问题：（1）附件1和附件2中，分别给出了该市

5、中心城区A的交通网络和现有的20个交巡警服务平台的设置情况示意图以及相关的数据信息。现要求为各交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地。对于重大突发事件，需要调度全区20个交巡警服务平台的警力资源，对进出该区的13条交通要道实现快速全封锁。实际中一个平台的警力最多封锁一个路口，请给出该区交巡警服务平台警力合理的调度方案。根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的情况，拟在该区内再增加2至5个平台，需要确定增加平台的具体个数和合适位置（2）针对全市（主城六区A，B，C，D，E，F）的具体情况，

6、按照设置交巡警服务平台的原则和任务，分析研究该市现有交巡警服务平台设置方案的合理性。如果有明显不合理，请给出解决方案。如果该市地点P（第32个节点）处发生了重大刑事案件，在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警服务平台警力资源的最佳围堵方案。二、模型的假设（1） 题中所提供的数据为真实数据。 （2） 简化问题，只考虑交巡警平台对路口节点的管辖，排除路口之间的关系。（3） 警车行驶正常，时速固定（60km/h）。（4）警车出警迅速，无停滞，道路保持畅通。三、符号的说明符号含义单位L从交巡警平台到达出事地块所行使的最大距离kmt从交巡警平台到达出事地块

7、的出警时间minv恒定的车速Km/h任意两点间的最短距离mm路线起终点之间的距离mm全市所有路口节点到这80个交巡警服务平台的最小距离mm一个282行2列的矩阵无四、问题的分析4.1 问题一的分析本问题要求根据附件2的数据为各交巡警服务平台分配管辖范围，进而给出该区交巡警服务平台警力合理的调度方案；然后结合现实条件，分析如何在A区增加服务平台，使交巡警服务平台工作量和出警时间更合理化，并找出要增加交巡警服务平台的区域，确定其具体位置。首先结合附件2中全市交通路口的节点数据和路线表格，利用MATLAB编程算出A区路口节点的起终点的路线距离，运用Floyd算法，计算得到任意两节点间的最短距离（92

8、*92）矩阵，A区92个节点的散点图见图41（程序见附录一）图41（1） 对第一段话的分析 本段所要解决的是20个交巡警服务平台管辖范围的划分问题。划分原则是交巡警服务平台到达所管辖路口节点的时间最少，即距离最短。利用任意两节点间的最短距离，找出2192路口节点分别到120交巡警服务平台的距离的最小值，这样，路口节点就归该距离最小的交巡警服务平台管辖；然后算出每个交警平台所管辖范围的总发案率，绘制成表格。（2）对第二段话的分析 发生重大突发事件时，需要调度全区20个交巡警服务平台的警力资源，对进出该区的13条交通要道实现快速全封锁。该题同样使用MATLAB解决最短路的编程问题。除去交巡警服务平

9、台所在的路口节点12、14、16，只剩下17个交巡警服务平台和10个要封锁的路口节点，即是找17个交巡警服务平台分别到10个要封锁的路口节点最短距离。根据实际条件：一个平台的警力最多封锁一个路口，因此若有一个交巡警服务平台同时达到两个以上的路口节点距离最短，就找附近的最短距离代替，通过统计最短距离，算出到各个路口的时间，从而确定交巡警服务平台警力合理的调度方案。（3）对第三段话的分析 由于现有交巡警服务平台工作量不均衡和部分路口节点出警时间过长，导致了交巡警平台调度的不合理。此问中我们通过在合适的路口节点增加适当数目交巡警服务平台，使交巡警服务平台的设置方案趋于合理化。这里我们给出的方案为，在

10、工作量大和出警时间长的交巡警服务平台上增加警点，这需要运用第一问的结果，交巡警平台的管辖范围，在管辖范围内，找出工作量大，出警时间长的警点（不超过五个），增加平台即可。4.2 问题二的分析研究该市现有交巡警服务平台设置方案是否合理可以从两个部分着手：工作量是否均衡；有些地方出警时间是否过长（或距离是否过大）。首先根据附件二中的数据，利用MATLAB编程算出A、B、C、D、E、F区路线起终点的距离，运用Floyd算法，计算得到全市任意两节点间的最短距离（ 582*582矩阵），全市582个节点的散点图见图42（程序见附录二）图42用MATLAB画图分析该市80个交巡警平台的工作量是否均衡。 定义

11、工作量为以下三种情况：工作量=负责节点的数目；工作量=发案率之和；工作量=发案率×距离之和。编写程序（程序见附录五）分别可以得到工作量与负责节点的数目、发案率之和、发案率×距离之和的条状图。如果图中纵轴方向高度波动较小 ，则此时工作量相对均衡，即现有交巡警平台设置方案合理，否则即为不合理，所以需要我们给出解决方案。P点发生重大刑事案件，案发三分钟后报警，交巡警平台按案发5分钟后，犯罪嫌疑人可能到达的所有路口节点部署第一道防线，根据警力，再部署案发10分钟，15分钟或者更短的时间间隔里部署更密集的第二道、第三道防线，由此得出调度全市交巡警服务平台警力资源的最佳围堵方案五、模型

12、的建立与求解5.1 问题一的模型建立与求解由题和附图1知，A区设有92个交叉口的节点，20个交警服务平台设置点，13个入城区的路口节点。则平常的路口节点（除去交警服务平台设置点）共有72个。如图3-1，其中120表示20个交警服务平台设置点，（程序见附录三）图5-1首先计算出A区起终点的路线距离1、根据附件2中所给的起点和终点的路线，结合各个节点的坐标位置，用MATLAB计算出两点之间的直线距离，得到92*92距离矩阵：2、根据附件2中的路线，我们可以得到各节点的邻接矩阵：即如果两个点相邻，则邻接矩阵中相对应的元素的值为1，否则为0；例如：8和9这两个点相邻，那么=1。3、根据Floyd算法，

13、为了求出任意两节点之间的距离，需要得到相邻两个节点的直线距离。可以利用距离矩阵的元素与的点乘积得到个各路线间的距离矩阵：4、我们可以将C中不相邻点间距离0改为无穷大（Inf）从而得到标志点与标志点间的权值矩阵： ，即如果8和9之间不相邻，也即不能直接到达，那么C中的=0和=0都将变成和等于无穷大（Inf），否则则等于C中相应元素的数据。其次计算得到任意两节点间的最短距离。有了上面的理论，去除大于92的数据，即不是A区的路段，结合各个节点的坐标位置，运用Floyd算法、MATLAB编程很容易就可求得任意两节点间的最短距离。5、运用Floyd算法求出任意两点间最短距离，得到最短距离矩阵：（1） 为

14、各交巡警服务平台分配管辖范围在MATLAB中编写程序（见附录四）即可筛选出附录一中120列到2192行的最小的最短距离，并记下此刻相对应的位置。分析整合数据就可得到表5-1表51交警平台编号交巡警平台位置标号管辖范围的节点个数管辖范围两点间的距离所管辖的此节点的发案率总发案率A11106716.1940.810.36812.0710.96951.17111.4031.17310.2960.9746.2651.1759.30050.87612.8361.1786.40310.8101.7A2273936.8221.49.74019.1441.74381.7449.48681.1708.60230

15、.97216.0620.8202.1A3355422.7090.95.65512.65916515.240.76618.4020.8302.2A4465718.6820.85.86017.3920.7623.51.26310.3081.46419.3630.8401.7A5594951.29.7508.48531.15112.2930.85216.5940.65311.7081.45620.8370.55823.0191.15915.2090.9502.1A661612.52.5A776305.8312.19.63211.4021.54712.8061.64812.9021.46141.902

16、0.6702.4A883338.27651.45469.30051.2802.4A9953120.5571.68.2345.02491.7354.24261.44510.9511.4902.1A101011001.61.6A111132691.22.02716.4330.81102.62.6A121222517.8892.42.4A131352127.0831.48.5229.05541.42352.42423.8541.11302.2A141411402.52.5A151532847.5181.34.82957.0051.41502.1A16164366.08281.15.03711.182

17、0.13834.0591.21602.6A17173418.51.45.3429.84891.41702.5A18185808.06230.86.1816.70821.48210.7931.1835.38520.91801.9A19193779.84890.83.4794.47210.81901.8A2020108411.752111.5854.47211.2863.60561.48714.6511.18812.9460.9899.48681.49013.0220.99115.9880.99236.0130.82001.9（2）给出该区交巡警服务平台警力合理的调度方案重大突发事件发生时，需调度

18、全区20个交巡警服务平台的警力对进出该区的13条交通要道实现快速地全封锁。由附件2知，位置编号为12、14、16的巡警服务平台恰好设置在13个路口中的其中3个，此时问题即可转化为17个交巡警服务平台的警力对该区10个路口的封锁的调度问题。 即位置编号为1、2、3、4、5、6、7、8、9、10、11、13、15、17、18、19、20的巡警服务台位对标号为21、22、23、24、28、29、30、38、48、62的10个路口节点进行快速全封锁。此时根据附录二的任意两点间最短路的数据结论，我们统计出关于交巡警服务平台到10个路口的距离的表格，其中列A为110个路口的距离；列B为210个路口的距离；

19、列S为2010 个路口的距离。对标号为21-62行的数据分别找出最小值，标出灰色。此时我们可以看出，离所要讨论的10个交通要道较近的交巡警服务平台集中在D,E,F,G,H,I,J,K,L,N列处，而其他相距较远的列A,B,C,M,O,P,Q,R,S数据的时间较长、与最短路距离相差较大，可以不予考虑。所以现只需统计D,E,F,G,H,I,J,K,L,N的数据记录即可，见表52表5245678910111315DEFGHIJK LN21182.73162.35162.65141.66126.99115.3995.11 50.72 27.08 165.6322200.76177.5177.8150.

20、36142.14131.3277.08 32.70 9.06 171.5123214.82191.55191.86164.42156.19145.3891.14 46.75 5.00 185.5624226.54182.85183.16155.72147.5136.6882.44 38.05 23.85 176.8728162.27113.07113.3785.702102.2897.75141.95 186.33 228.08 47.51829155.35106.15106.4680.155104.93107.24151.44 195.82 237.57 57.0053081.0331.82

21、932.1355.83130.6034.9279.11 123.50 165.25 44.0153885.7257.7858.0830.41425.9314.3468.58 112.97 138.33 64.584873.9524.75825.0612.90230.9941.9986.19 130.57 172.32 51.086623.552.55153.3779.91786.7793.36147.61 191.99 213.32 118.1已知实际中一个平台的警力最多封锁一个路口，但图中G、L、N 三列有重叠的数据，即一个平台的警力需要负责多个路口，不满足实际条件，因此需要在同一行中的最小

22、值附近进行适当的调整调度，使其合理。 调整结果如表53的阴影部分，为警力平台封锁对应路口的最短距离。表5345678910111315DEFGHIJK LN21182.73162.35162.65141.66126.99115.3995.11 50.72 27.08 165.6322200.76177.5177.8150.36142.14131.3277.08 32.70 9.06 171.5123214.82191.55191.86164.42156.19145.3891.14 46.75 5.00 185.5624226.54182.85183.16155.72147.5136.6882.

23、44 38.05 23.85 176.8728162.27113.07113.3785.702102.2897.75141.95 186.33 228.08 47.51829155.35106.15106.4680.155104.93107.24151.44 195.82 237.57 57.0053081.0331.82932.1355.83130.6034.9279.11 123.50 165.25 44.0153785.7257.7858.0830.41425.9314.3468.58 112.97 138.33 64.584873.9524.75825.0612.90230.9941.

24、9986.19 130.57 172.32 51.086213.552.55153.3779.91786.7793.36147.61 191.99 213.32 118.1各个路口调度的最短时间表54各个路口的调度信息45678910111315DEFGHIJKLNS(100m)3.5106.1525.065.831142.1414.3482.4445.7527.0847.518t(min)0.35106152.5060.583114.2141.4348.2444.5752.7084.7518由表（54）知在紧急情况下，交巡警到达各个路口的时间，其中最长时间14.214min。（3）确定需要增

25、加平台的具体个数和位置根据表51我们得到120交巡警平台标号，及管辖的路口节点数目，及发案率总和的表格如下表5-5交巡警平台位置标号管辖范围的节点个数总发案率交巡警平台位置标号管辖范围的节点个数总发案率11010.31132.0279.71222.6355.61352.4465.81418.5599.71534.8612.51645.0769.61736.18351855.3958.21933.41011.6201011.5从表格中的数据可知，根据方案一，要增加五个交巡警服务平台，它们的位置分别在标号为1、2、5、7、20的交巡警平台处。5.2 问题二的模型建立与求解 （1）研究该市现有交巡警

26、服务平台设置方案（参见附件）的合理性1） 首先分析定义工作量为平台管辖的数目时，交巡警平台与负责节点的数目的图像关系。根据Floyd算法，得到全市任意两个路口节点之间的最小距离，为582*582的对称矩阵。依据附件2通过MATLAB编程，得到所有路口节点（包括交巡警服务平台）到这80个交巡警服务平台的最小距离（582*80），并找到矩阵中每行的最小值，也即计算出每个路口节点到80个交巡警服务平台的最短距离，该路口节点就归此交巡警平台管辖。将这些数据存放到中，第一列表示第i个路口节点到交巡警平台的最小距离，直到记录完所有到交巡警平台的编号，由中第二行交巡警平台的编号，可以知道每个交巡警平台管辖的

27、路口节点的个数，通过编程，可以得到交巡警平台与管辖的路口数目的图像，横坐标180是交巡警平台的序列号，如下图53所示图52容易知道，如果一个交巡警服务平台所管辖的路口节点越多，它的工作量就会越大，此图在一定程度上反映现有交巡警服务平台设置方案合理与否，但考虑不完全。由图可知交巡警平台序列号为22、24、26、34、41、42、43、47、65、66、70、71、72、73工作量比较大（所管辖路口节点数目超过12个）对应的交巡警平台的位置标号是94、96、98、171、178、179、180、321、382、383、475、476、477、478。2） 然后分析定义工作量为发案率之和时，交巡警平

28、台与发案率之和的图像关系。由附件2中全市交通路口节点数据可知发案率的平均次数越大，所需要的警力就会越大，即交巡警服务平台的工作量也就越大，反之，工作量越小。根据附件2给出的发案率，来计算每个交巡警平台管辖范围的总发案率。运用MATLAB编程，得到交巡警平台和发案率总和的关系图，如图54。由图可知交巡警平台序列号为22、30、34、36、42、43、47、62、65、66、70、71、72、73工作量比较大（平均发案次数超过13.0）对应的交巡警平台的位置标号是94、167、171、173、179、180、321、379、382、383、475、477、478、479.图533） 最后分析定义工

29、作量为时发案率×距离之积，交巡警平台与发案率×距离之积的图像关系。警车的速度恒为60km/h，重大事故发生点离所属的交巡警平台越近，出警时间就越短，工作量就越少，反之，就越远。考虑到最短距离和发芽率，用MATLAB编程得到交巡警平台和发案率与最短距离之积的关系图，如图55。由图可知交巡警平台序列号为28、30、34、36、38、41、42、43、46、62、63、65、68、70、72、73、74工作量比较大（发案率和最短距离之积大于200）对应的交巡警平台的位置标号是100、167、171、173、175、178、179、180、320、379、380、382、385、4

30、75、477、478、479.根据图形的分析，条状图的纵坐标的差值还是很大，比较大的点已经在上述中给出，可知交巡警服务平台设置方案是很不合理的。此时我们提出两种解决方案，这两种方案并用解决问题。第一种是增加交巡警服务平台，第二种是修路。根据上面的分析，可知序列号为22、24、30、34、41、42、43、47、62、65、66、70、71、72、73交巡警平台，对应的交巡警平台的位置标号是94、96、167、171、178、179、180、379、382、383、475、476、477、478.满足三种工作量中的任何两个都较大情况.则在这些交巡警平台上再建立一个交巡警平台，用MATLAB编程，

31、画出增加过交巡警平台后，三种工作量与交巡警平台的条状图。见图56、56、57其中图示第一列均为坐标参考值（程序见附录六）图54图55图56此时的纵轴分布较为集中，即工作量相对集中，此方案合理。这就说明，改进的方案可行。通过观察交巡警平台与管辖的路口数目的图像，我们可以看到，有一些交巡警平台管辖的路口数目为1，也就是只管辖自己本身，这种情况的产生可以认为是由于交巡警平台所处地理位置不优越的原因，所以采用修路解决。（2）为了快速搜捕嫌疑犯，给出调度全市交巡警服务平台警力资源的最佳围堵方案P点发生重大刑事案件，案发三分钟后报警，交巡警平台按案发5分钟后，犯罪嫌疑人可能到达的所有路口节点部署第一道防线

32、，根据警力，再部署案发10分钟，15分钟或者更短的时间间隔里部署更密集的第二道、第三道防线。根据由floyd算法求出的任意两路口节点之间的最小距离，求出在相应时间内犯罪分子可能到达的路口节点，然后对其封锁，若此防线搜捕失败，警力立即调配到下一个防线进行搜捕。由此得出调度全市交巡警服务平台警力资源的最佳围堵方案六、模型的优缺点6.1模型的优点该模型应用附件2，全市交通路口的路线表格中的数据，通过Floyd算法，算出任意两个路口节点之间最短距离，具有真实性，而不是附件1中A图两点间的直线距离，这更加符合实际情况，生活中，任两个路口节点由于地理环境的原因，并不一定都有之间相连的直线通道，这为以后求解

33、问题的可靠性做了保证。整个模型的建立在实际的地图和图表数据，具有很好的现实检验意义。在模型建立的过程中遵循科学的数学逻辑和实际的生活逻辑。该模型通过解决最短路问题建立出了在理想假设下的交巡警平台的最优设置，使出警时间最短化，并且对实际生活中的交巡警服务台的调度具有一定的实际参考价值，对提高交巡警的执法效率有着不可或缺的影响。6.2模型的缺点交巡警服务平台设置的优化模型在实际情况下具有一定的局限性，模型分析中我们以时间不超过3min、警车以速度60km/h匀速行驶作为模型求解的数值条件，但现实中道路拥堵等情况的发生势必会影响警车速度，从而会使实际时间偏大。为了更贴近实际，则应考虑地区人口密度对道

34、路畅通性的影响，即对速度的影响。另外，现实中事故发生点不一定在路口节点处，同时也忽略了实际生活中存在的诸多不定因素，这些都是模型的局限性所在。七、参考文献1 楼世博 金晓龙 李鸿祥等 图论及其应用 M.北京：人民邮电出版社，1982年2 姜启源 谢金星 叶俊 编 数学模型（第三版） M.北京：高等教育出版社，20043 林雪松 林德新等 MATLAB7.0应用集锦 M.北京：机械工业出版社，2006年4 建模论文，警务资源；警务资源合理配置问题研究 2011-9-115杨启帆 康旭升. 编 数学建模M.北京：高等教育出版社，20056 刘胜 范玉顺 尹朝万 基于工作流模型的资源配置优化方法，计

35、算机集成制造系统，第11卷第9期，2006年9月 八、附录附录一clc,clear allp=xlsread('E:/cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口节点数据','A2:A93')x=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口节点数据','B2:B93')y=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口节

36、点数据','C2:C93')plot(x,y,'r*')B=x ylen=length(x)for i=1:len for j=1:len S(i,j)=(B(i,1)-B(j,1)2+(B(i,2)-B(j,2)2)0.5; endendSm=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口的路线','A2:A144')n=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通

37、路口的路线','B2:B144')M=m nf=find(n<=92)lenf=length(f)for i=1:lenf %去除大于92的数据，即不是A区的路段 mm(i)=m(f(i);nn(i)=n(f(i);endm=mm;n=nn;A=sparse(len,len)for t=1:lenf i=m(t);j=n(t); A(i,j)=1;A(j,i)=1endAfor i=1:len for j=1:len C(i,j)=A(i,j)*S(i,j); endendCfor i=1:len for j=1:len if C(i,j)=0 D(i,j)=10

38、00; else D(i,j)=C(i,j) end endendfor i=1:len D(i,i)=0endfor k=1:len for i=1:len for j=1:len if D(i,k)+D(k,j)<D(i,j) D(i,j)=D(i,k)+D(k,j); D(i,j)=D(i,k); end end endendDD的数据如下（插入的excel表格）附录二clc,clear allp=xlsread('E:/cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口节点数据','A2:A583')x

39、=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口节点数据','B2:B583')y=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口节点数据','C2:C583')plot(x,y,'r*')B=x ylen=length(x);for i=1:len; for j=1:len; S(i,j)=(B(i,1)-B(j,1)2+(B(i,2)-B(j,2)2)0.5;

40、endendS;m=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口的路线','A2:A929')n=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口的路线','B2:B929')M=m n;lenf=length(m);A=sparse(len,len);for t=1:lenf; i=m(t);j=n(t); A(i,j)=1;A(j,i)=1;endA;for i=1:len;

41、for j=1:len; C(i,j)=A(i,j)*S(i,j); endendC;B=ones(len)-A;D=C+B*106;for k=1:len; for i=1:len; for j=1:len; if D(i,k)+D(k,j)<D(i,j); D(i,j)=D(i,k)+D(k,j); end end endendD=D-diag(diag(D);D由于D为582*582对称方阵，数据量特别大，现只给出一小部分附录三clc,clear allm=1 1 2 3 3 4 4 5 5 6 7 7 8 8 9 10 11 11 12 12 14 15 15 16 16 17

42、17 1718 18 19 20 21 22 22 23 23 24 24 25 26 26 27 28 28 29 30 30 31 31 3233 33 34 35 36 36 36 36 37 38 38 39 40 41 41 42 43 43 44 45 46 46 4747 47 48 49 49 50 51 51 52 53 53 54 54 55 56 57 57 57 58 60 61 62 62 6364 64 65 66 66 67 67 68 68 69 69 69 70 70 71 71 72 73 73 74 74 75 7677 77 78 79 80 81 82

43、 82 83 84 85 86 86 87 87 88 88 89 89 89 90 91;%路径数据n=75 78 44n45 65 39 63 49 50 59n32 47 9 47 35 34 22 26 25 471 21 7 3114 38 40 42 81 81 83 79 86 22 372 13 13 383 13 25 11 27 10 12 29 15 307 48 32 34 33 34 8 9 45 35 37 16 39 7 39 41 40 2 17 92 43 2 72 3 46 855 48 6 5 61 50 53 51 52 59 56 52 54 55 6

44、3 3 57 58 60 4 59 62 60 4 8564 65 76 66 67 76 44 68 69 75 70 71 1 2 43 72 74 73 74 18 1 80 76 7778 19 79 80 18 82 83 90 84 85 20 87 88 88 92 89 91 20 84 90 91 92;for t=1:143 i=m(t);j=n(t); A(i,j)=1;A(j,i)=1endAB=1 413 359;2 403 343;3 383.5 351;4 381 377.5;5 339 376;6 335 383; 7 317 362;8 334.5 353.5

45、;9 333 342;10 282 325;11 247 301;12 219 316;13 225 270;14 280 292;15 290 335;16 337 328;17 415 335;18 432 371;19 418 374;20 444 394;21 251 277;22 234 271;23 225 265;24 212 290;25 227 300;26 256 301;27 250.5 306;28 243 328;29 246 337;30 314 367;31 315 351;32 326 355;33 327 350;34 328 342.5;35 336 339

46、;36 336 334;37 331 335;38 371 330;39 371 333;40 388.5 330.5;41 411 327.5;42 419 344;43 411 343;44 394 346;45 342 342;46 342 348;47 325 372;48 315 374;49 342 372;50 345 382;51 348.5 380.5;52 351 377;53 348 369;54 370 363;55 371 353;56 354 374;57 363 382.5;58 357 387;59 351 382;60 369 388;61 335 395;6

47、2 381 381;63 391 375;64 392 366;65 395 361;66 398 362;67 401 359;68 405 360;69 410 355;70 408 350;71 415 351;72 418 347;73 422 354;74 418.5 356;75 405.5 364.5;76 405 368;77 409 370;78 417 364;79 420 370;80 424 372;81 438 368;82 438.5 373;83 434 376;84 438 385;85 440 392;86 447 392;87 448 381;88 444.

48、5 383;89 441 385;90 440.5 381.5;91 445 380;92 444 360;%坐标数据 for i=1:92 for j=1:92 S(i,j)=(B(i,2)-B(j,2)2+(B(i,3)-B(j,3)2)0.5; endendSA for i=1:92 for j=1:92 C(i,j)=A(i,j)*S(i,j); endendC附录四for i=1:lenD(i,i)=0;endd=D(21:92,1:20);mind=min(d');num=;for i=1:72aa=find(d(i,:)=mind(i);num=num aa;end附录五

49、clc,clear allp=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口节点数据','A2:A583')x=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口节点数据','B2:B583')y=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口节点数据','C2:C583')pl

50、ot(x,y,'r.')B=x ylen=length(x)for i=1:len for j=1:len S(i,j)=(B(i,1)-B(j,1)2+(B(i,2)-B(j,2)2)0.5; endendSm=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口的路线','A2:A929')n=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交通路口的路线','B2:B929'

51、)M=m nlenf=length(m);A=sparse(len,len)for t=1:lenf i=m(t);j=n(t); A(i,j)=1;A(j,i)=1endAfor i=1:len for j=1:len C(i,j)=A(i,j)*S(i,j); endendCB=ones(len)-AD=C+B*106for k=1:len for i=1:len for j=1:len if D(i,k)+D(k,j)<D(i,j) D(i,j)=D(i,k)+D(k,j); end end endendD=D-diag(diag(D);Dg=xlsread('E: /cumcm2011B附件2_全市六区交通网路和平台设置的数据表','全市交巡警平台','B2:B81')leng=length(g);E=zeros(len,leng)for i=1:leng E(:,i)=D(:,g(i);endS=ze