2011年数学建模B题国家一等奖

小时

2011高教社杯全国大学生数学建模竞赛
, {& G( }2 L# u  G1 o
+ H8 f* M8 S2 l9 J5 a% v9 I* I承  诺  书
- F8 m3 l' K! C9 |. n
5 J. D  |- t$ z' z% m: W我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
- E# ?9 |2 H, b; a, {2 J. X我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题.
我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出.
5 m9 _7 [, q( k- A/ ?$ |我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性.如有违反竞赛规则的行为,我们将受到严肃处理.
8 a# Q" _6 g2 N8 K+ U7 x
% a3 J4 ?6 f- \  a$ Q4 o! a8 I( T我们参赛选择的题号是(从A/B/C/D中选择一项填写): B                     
        我们的参赛报名号为(如果赛区设置报名号的话):                           
所属学校(请填写完整的全名):            
参赛队员 (打印并签名) :1.                                       
                       2.                                       
                       3.                                         
+ ?5 N. E. i! I8 p8 O% x5 M0 E指导教师或指导教师组负责人  (打印并签名):               

                                              日期:  2011  年  9 月 11 日





赛区评阅编号(由赛区组委会评阅前进行编号):

6 q1 |0 Z  e/ o3 }2011高教社杯全国大学生数学建模竞赛
  Y. }* i( L/ U3 k8 F
! ~) ~$ g/ v8 n编 号 专 用 页
& p( T9 B8 _, C; |- l: p4 w: p

/ N+ k( O( g0 Y8 g; z
% J2 }" P( D9 M9 u2 q赛区评阅编号(由赛区组委会评阅前进行编号):

, C& l2 J" {) F0 S8 `
3 A# l7 r0 ?9 L( Q% m
& C: s. t) `; ?+ c4 L3 ]赛区评阅记录(可供赛区评阅时使用):
8 l! V% z7 B) h( r. i* p
评
阅
; ?) l) ?: \7 d0 o人
                                                                               


5 e( |% C+ d) i评
5 _+ }2 H+ `$ j# e9 e8 G3 h- |5 F. f分
                                                                               
2 _: r* ^3 Z2 M# |
# D; D' ?* m2 `# C& i3 n; U备
注
9 ?6 Q1 T) V$ M6 b4 W' V                                                                               

/ ?; y+ T' @3 h7 N, _: j# D5 ~
- m, O' p) p2 K- \2 y; ^& A
& t% p0 |$ p6 t- \+ F  Z
全国统一编号(由赛区组委会送交全国前编号):


, C! y) @# F3 U" n" z( ]
6 c- q0 \- q, u5 M
  K  U! M- W1 A2 ^
' Z0 o5 l' v- D9 X+ E7 p& @全国评阅编号(由全国组委会评阅前进行编号):
& w5 J* f, i8 t0 X4 u* a3 Q
9 N8 _( Y% ^9 t/ T% r# b- V


0 r* G8 \' g1 Z
交巡警服务平台的设置与调度
5 i' T" A  v( W2 Q+ p2 g) I摘   要
由于警务资源是有限的,所以根据城市的实际情况与需求,合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是有关部门面临的一个实际课题.本文着力于通过所给资料,寻找最优化的交巡台设置与调度方案. 按照设置交巡警服务平台的原则和任务,我们首先对问题1用Floyd算法,提出最佳的交巡警服务平台管辖区域划分方案,缩短了出警时间,平衡了工作量,然后采用回溯法,给出了应对突发事件的警力比较合理调度方案；对于问题2,我们将其归结为全局的配置问题,首先用优化后的Floyd算法对该市现有六城区的交巡警服务平台设置进行改进,其次以时间最短、围堵区域最小为原则,提出了应对重大刑事案件的最佳围堵方案.
% R+ s0 ?) V/ ^' v- p! O对于问题1,本文将最短时间问题转化为单向最短路径问题.我们没有运用经典的求最短距的Dijkstra算法,采取时间复杂度更简便的Floyd算法,应用Matlab编程,以出警时间最短为原则,将72个交通节点分配给20个交巡警服务平台；对于出现突发事件,本文采用回溯法,以最节省警力、实现全区封锁联动时间（即封锁路口最长时间）最短为目标,成功的实现了应对突发事件时警力的合理调度；对于某些交巡警服务平台工作量大、出警时间过长等问题,本文利用Mathematica对附表2中的数据进行分析,整理分析A区各节点事故发生率后,利用图论的相关知识,提出应增设4个服务平台,基本实现警力的最优配置.最后,借助于Matlab和Mathematica软件,对附件中所提供的数据进行了筛选,去除异常数据,对残缺数据进行适当补充,并从中随机抽取了3组数据(每组8个采样)对理论结果进行了数据模拟,结果显示,理论结果与数据模拟结果吻合良好.
3 b& R) Y4 g- F7 }$ v: G; z而对于问题2,我们对附件中所提供的A,B,C,D,E,F六城区的数据进行了整合与分析,并做出了直观的图表.遵循警情主导警务原则、快速出警原则、方便与安全原则,并结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况,在充分考虑现有警力和财力并确保安全的条件下,科学分析现有平台的数量和具体位置的合理性.数据显示C区和F区的事故发生率较高、交巡警服务平台工作量高于全市平均水平、交巡警服务平台平均每天出警时间过长,针对以上问题我们再次利用均衡二分法,并考虑区域边界处的设点拥挤问题,提出了在C区增设5个交巡平台、F区增设1个交巡平台.对于该市地点P(第32个节点)处发生了重大刑事案件的围堵问题,本文将其归结为资源调配问题.本文合理假设了犯罪嫌疑人的车行驶速度(分三种情况考虑:等于警车速度,警车速度的二倍,警车速度的一半),确定三分钟后犯罪嫌疑人逃逸的可能覆盖范围,从而利用回溯法的思想采用Matlab编程确定犯罪嫌疑人的车的所有可能位置.以时间最短、围堵区域最小为原则,采用改进的穷举算法,快速地形成围堵区域,并实现了围堵区域最小的目的.实现了资源调配问题的优化决策.
考虑到该城市未来发展规划,只需对本文所建模型进行适当改进即可,在此不进行详细解答.

7 I  v& v( q# q- u; `) |  Z
" \6 Z6 w4 H2 S) x% Z0 \, n0 [关键词    最短路径  Floyd算法  回溯法  穷举法  优化决策


* g4 A+ N% K- z3 s/ f7 N% u5 C# G
目   录
交巡警服务平台的设置与调度        1
摘   要        1
1 q! B/ |- }, n, e9 N, N" ]( [( r1.问题重述        1
) w; W/ o1 k& ]" r1 q+ o7 D( ~2.问题分析        1
2.1对于问题一的分析        1
2.2对问题二的分析        1
3.模型假设        2
4.定义与符号说明        2
  K) q- p: n. a* n, x- S5.模型的建立与求解        2
5.1 问题一的模型        2
5.1.1 模型建立        2
5.1.2 模型求解        3
: i6 m6 w1 W, p- i2 p5.2 问题二的模型        8
5.2.1 模型建立        8
  M3 ?5 ^0 x! k5.2.2 模型求解        9
* c" \- S7 v; C7 g* c0 g4 x7.模型的评价与推广        10
: c% ~6 N3 T- ]2 ^; e6 o7 k8. 附件        10
( E7 Y7 p) P* F6 t7 I! h! Q附件1:用Floyd算法分配个服务平台管辖区域        10
附件2:邻接矩阵的matlab实现程序        22
附件3:围堵方案的java实现程序        29
附件4:全区的交巡警平台有效覆盖范围(有效代表三分钟内可以到达)        30
附件5:用Mathmatica求数据均值与方差        30
7 v( Q7 y* {3 w5 N  H( h附件6:输入任意两点的坐标,输出两点间距离        30
附件7:A区各线路距离        31
+ I1 b$ V8 |2 v; i
) x  f. s! X% J$ Z


: C+ l- F% G. T4 |, E: g

3 G% u/ V* ~* h) N
1.问题重述
0 d* R; _  l7 x( |7 |/ @1 c“有困难找警察”,是家喻户晓的一句流行语.警察肩负着刑事执法、治安管理、交通管理、服务群众四大职责.为了更有效地贯彻实施这些职能,需要在市区的一些交通要道、人员密集区和重要部位设置交巡警服务平台.每个交巡警服务平台的职能和警力配备基本相同.由于警务资源的有限性,根据城市的实际情况与需求,合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题.本文着力于寻找最优化的设置与调度方案.
1 e. |0 G' N: E' S问题1要求合理分配交巡警服务平台的管辖范围,使其在所管辖的范围内出现突发事件时,尽量能在3分钟内有交巡警(警车的时速为60km/h)到达事发地；对于重大突发事件,给出该区交巡警服务平台警力合理的调度方案,尽快封锁道路；拟在该区内再增加2至5个平台,以减少出警时间、平均工作量,确定需要增加平台的具体个数和位置.
问题2要求分析研究该市现有交巡警服务平台设置方案的合理性并给出解决方案；如果该市地点P(第32个节点)处发生了重大刑事案件,在案发3分钟后接到报警,犯罪嫌疑人已驾车逃跑.为了快速搜捕嫌疑犯,给出调度全市交巡警服务平台警力资源的最佳围堵方案.
2.问题分析
* k6 b7 {9 a4 p" J8 w- S本题所要解决的是A区以及全市的安巡警服务平台设置与调度问题,根据现实生活状况,我们首先要考虑的是警力资源的限制,即要使得所布置的警力尽可能的少.其次是在交巡台数量最少的情况下,力求警员到达现场的时间在3分钟以内,解决突发状况.
2.1对于问题一的分析
$ w) v5 g! Z7 M$ y该市中心城区A的交通网络有92个节点和20个交巡警服务平台,要求当突发事件发生时,尽量能在3分钟内有交巡警到达事发地,已知警车的时速为V=60km/h,我们将最短时间转化为最短路问题,应用Floyd算法,求解出A区距离每一节点最近的交巡台,即将该节点分配给该交巡台.
5 P2 _, J# o! }5 J& b$ E9 B对于重大突发事件,要实现对进出该区的13条交通要道进行快速封锁,即需调度交巡台尽快到达13个节点,重复Floyd算法,找出最近交巡台,即可找出调配方案.但需注意的是,有的出入口本来就有交巡台,但为了达最优化,需进行重新分配,故应用回溯法,找到调度方案.
现有交巡台工作量不均衡和有些地方出警时间过长,统计A区各个交巡台案发率,计算均值与方差,在案发率较高地带增设交巡台,平衡工作量,尽量缩短出警时间.
( x5 H- o  c; a  C/ L6 J$ ]5 N2.2对问题二的分析
! K0 R" b0 p7 N1 p8 a- m# O$ U5 T对于问题二,是对问题一的进一步改进与推广,在遵循警情主导警务原则,快速出警原则与方便与安全原则,结合辖区地域特征、人口分布和治安状况等实际情况,充分考虑现有警力和财力并确保安全的条件下,设置交巡平台,重复上一问的做法,评估交巡平台的合理性.对于改进方案,应考虑城区内部工作量,城区之间的联系以及城市边界的警力调度.
对于突发状况的围堵方案,应在最短时间内对可能逃逸区域进行合围,最小范围内缩小包围圈.
8 U9 U% l' L# A0 p& Y3.模型假设
  A- D% T' T3 F3 u9 d$ S1.假设题中所给数据均真实可靠.
, J6 ~( L  u$ S. l2.出警时道路恒畅通(无交通事故、交通堵塞等发生),警车行驶正常,警车及肇事车辆行驶时均以60km/h匀速行驶,转弯处不需要花费时间.
+ R  n+ @9 i( N0 J3.事故均发生在路口节点,两节点连线上认为没有事故发生.
/ ]' t6 Y, R5 o0 B+ M  d4.每条线路行驶都是双向的.
5.考虑肇事车辆在P点向各个方向逃逸的概率相等.
4 e8 _6 j1 c0 a% e6.在整个行驶中,车辆只在主要干道行驶.
7.发生事故时,忽略反应调度时间.
4.定义与符号说明           

# X7 r* M6 [0 d  Z9 l) c6 p7 L任意两个标志点 与 之间的距离
9 Z" ~8 f( r- t, N$ Q8 j
, d; a8 D. P$ t) W$ Y
) p$ @) Q! F) }标志点间的距离组成的距离矩阵

标志点的邻接矩阵

3 h) I, Q* I+ p! f/ c# N6 a邻接矩阵的元素
0 R+ v! K5 ?- Q7 Y. ?, h/ |
相邻标志点间的距离矩阵

( n& @+ m) o8 M% f$ o' ?相邻标志点 与 间的距离
) T- c1 A# ~- m# @, H; m

标志点的权值矩阵
  
2 I6 ^; j4 `+ }+ ^( [( Y7 g; n& \标志点间的最短距离矩阵
- q4 ]* z3 L  g
4 j* |* [; }3 n3 ^; T2 B# m, q0 A$ w! {标志点 与 之间的最短距离


/ \5 G" V8 x+ c( f  V' z* `9 e肇事车辆逃逸速度
" A+ \) L* {, A1 W+ }" N- u
  
: h+ F: O& C  i: P+ z5.模型的建立与求解
" m7 z% t4 _/ C+ r- N% y* ^8 V7 d5.1 问题一的模型
5.1.1 模型建立
此问是关于最短路径的模型分析及MATLAB的实现A区道路状况及交巡台的设置如图1所示.本文应用Floyd算法,通过构造距离矩阵,依次找出距离每一节点最近的交巡台,使得有事故发生时,交巡警在最短时间内到达事故现场,以此为依据分配管辖区域.如果道路不通时,认为两端节点的距离为无穷.
                   图1  A区各节点及服务平台示意图
当有重大突发事件时,要对进出该区的13条交通要道进行快速封锁,固定13个出入口,应用回溯法,找到距离节点最近的交巡平台.封锁时间决定于最后到达节点的时间,由于一个平台的警力最多封锁一个路口,至少需调动13个平台的警力.
7 B$ ?) }2 f6 e& y& E1 M" i% ?, P7 e为达到工作量的均衡和出警时间尽可能的短,需进行优化决策.考虑每一节点案发率的不同,在A区增设2到5个平台,使得每一平台的工作量均衡,平均出警时间大体相同.
5.1.2 模型求解
首先我们可以根据题中所给的各个标志点的坐标,用matlab计算出任意两点之间的直线距离,得到92*92的距离矩阵:
( [( r7 \: R. h. {( H8 s/ F% `
# ?- b' i/ ^5 a+ y根据题中的分布图,我们可以得到各标志点的邻接矩
+ D$ L& q' P" Z8 u7 _" J
即如果两个点相邻,则邻接矩阵中相对应的元素的值为1,否则为0；例如:3和44这两个点相邻,那么  .
    根据Floyd算法,我们是要求出任意两节点之间的距离,所以我们需要得到相邻两个结点的直线距离.我们可以利用距离矩阵的元素 与 的点乘积得到相邻标志点间的距离矩阵:

/ X7 X2 m) Y1 H+ _0 Q" O' ?对于D中不相邻点间距离0改为无穷大(Inf)从而得到节点与节点间的权值矩阵:

" ^  v6 H: \9 B% s( g, t即如果15和10之间不相邻,也即不能直接到达,那么D中的 和 都将变成 和 等于无穷大(Inf),否则则等于D中相应元素的数据.
' M4 O* C3 {1 n* M1 K运用Floyd算法求出任意两点间最短距离,得到最短距离矩阵 :
! \8 l4 n: o; w
+ e! q" K3 h6 `  z* _  P由Floyd算法,运行MATLAB程序,可统计出距离每一节点最近的交巡台的位置,MATLAB运行结果如表1所示.带括号的节点为发生事故时任意交巡台都不能在三分钟内赶到节点.
% a  R0 o3 X  U% C3 c' f) K2 X- ~交巡台—节点        距离        交巡台—节点        距离
  A+ ~5 D+ u. V13—21        27.0831         4—57        18.6815
13—22        9.0554        6—58        23.8414
3 `3 W! _# G2 K+ u13—23        5.0000         6—59        16.0312
13—24        23.8537         4—60        17.9240
12—25        17.8885         4—(61)        52.1055
3 L8 l, S/ j& |0 R# a2 ^+ f/ `11—26        9.0000         4—62        3.5000
; x( w9 K0 A2 u' x" d  q/ V11—27        16.4330         4—63        10.3087
15—(28)        47.5184         4—64        9.3632
% ~0 X# E7 U, B4 }15—(29)        57.0052        3—65        15.2398
2 Q& `2 ?- }6 [  x" _) S7—30        5.8310         3—66        18.4012
9—31        20.5572        1—67        14.9158
7—32        11.4018        1—75—68        10.7927
8—33        8.2765        1—69        5.0000
) `3 H, w' C- M9—34        5.0249         2—70        8.6023
9—35        4.2426         1—74—71        11.2650
16—36        6.0828         2—72        16.4031
16—37        11.1818         18—73        19.7231
2 F' g7 e+ s) `16—(38)        34.0588         1—74        6.2650
2—(39)        36.8219        1—75        6.2650
2—40        19.1442         1—76        9.8005
1 R: R, i' L3 P2 {17—41        8.5000         19—77        9.8489
17—42        9.8489        1—78        6.4031
2—43        8.0000        19—79        4.4721
2—44        9.8468        18—80        8.0623
/ E$ K' z. N  _, ?9—45        10.9508         18—81        6.7082
8—46        9.3005        18—82        10.7935
0 B. G/ B& m% r0 H7—47        12.8062        18—83        5.3852
7—48        12.9021        20—84        11.7522
/ r: c1 ~  ?- R. L5—49        5.0000        20—85        4.4721
1 V; L' J6 }. a1 c5—50        8.4853        20—86        3.6050
5—51        12.8932        20—87        14.6511
5—52        17.1944        20—88        12.9464
7 q% v! }7 V1 o4 U' i5—53        11.7082        20—89        14.7522
3—54        22.7089        18—90        19.5256
! U" y% [. v5 Q# k* V8 }  H3—55        12.6590        20—91        16.0060
5—56        21.4370        20—(92)        36.0060
表1  该市A区指定节点到交巡警服务平台最短距离
由上表可初步确定A区20个交巡台的管辖范围,如表2所示.带括号的节点为发生事故时任意交巡台都不能在三分钟内赶到节点.
' Q- x) o+ y. C2 E. D交巡台序号        辖区内节点        辖区内案发率        交巡台序号        辖区内节点        辖区内案发率
% x# }8 H) T, [  _5 Y6 f1        67 68 69 71 74 75 76 78         9.4        2         40 43 44 70 72 39        9.7
3        54 55 65 66        5.6        4        57 60 62 63 64        6.6
5        49 50 51 52 53 56        7.7        6        58 59        4.5
7        30 32 47 48 61        9        8        33 46        5
, G! G: x4 ^1 y6 i9        31 34 35 45        8.2        10                 1.6
% l) n2 g: w; h) Z$ `0 H11        26 27        4.6        12         25        4
, k, |* j' S- C+ T+ }: r# i9 `, B13        21 22 23 24        8.5        14                 2.5
15        (28) (29)        4.8        16        36 37 (38)         5
17        41 42        5.3        18        73 80 81 82 83        7
19        77 79        3.4        20        84 85 86 87 88 89 91 90 (92)        11.5
表2  该市A区交巡警服务平台所管辖交叉路口清单

; M* s1 F3 U, i
/ l9 c/ y$ G( A4 m# J图2  A区各交巡台管辖区域示意图
    需要说明的是,同一条路整体归一个交巡台管理.
3 w& N1 X& g6 H7 e( c/ G+ r9 ]当有重大突发事件时,固定13个进出A区的节点,运用回溯法,结合上表,找到距离节点最近的交巡台,以此来达到总体时间的最短,我们一共可以得到四个方案,在这个过程中可以发现,有些交巡台要避免去最近的节点封锁而去较远的节点,以此来节省警力.具体封锁方案如表3、表4所示.最短调度时间均为8.0155.
- }% ?+ r- {2 y' m0 C) f方案一:
0 \* _: r0 u/ U( H  s5 I+ F交巡台        过程        出入口节点号
. }& G4 H/ }3 h, S0 D3 `2        路径        40 39
38
8 r+ L' ~: L  R( S4 U: C# X2 e: [        时间        3.9822min       
4        路径                62
& v' n$ M$ q$ v1 b* _  Q        时间        0.3500min       
3 O: r5 D: q4 i7 ^) o6        路径        47 48        30
/ K/ u. x6 a: J7 a, |) C        时间        3.1829min       
7        路径        30        29
3 a# U4 r9 U9 G  U( m5 }+ Q0 W        时间        8.0154min       
8        路径        47        48
        时间        3.0995min       
9        路径        35 36        16
        时间        1.5083min       
10        路径        26 27        12
' {9 I3 U, u( Y; H3 F        时间        7.5863min       
11        路径                22
        时间        3.2696min       
12        路径        25        24
3 C2 B2 `; M) r7 S# B0 n        时间        3.5916min       
" G- z& v* W( k# W13        路径                23
        时间        0.5000min       
14        路径                21
        时间        3.2649min       
15        路径                28
        时间        4.7518min       
16        路径                14
        时间        6.7417min       
表3  A区突发事件封锁方案一
; @: o4 s0 j8 \! Z9 M$ L6 c$ I: }方案二        方案三        方案四
路口标号        平台号        路口标号        平台号        路口标号        平台号
12        13        12        13        12        10
2 ~$ N# ^. p+ L; D  e& k14        16        14        23        14        16
16        6        16        9        16        6
4 I+ B% D* {' [3 K+ p21        14        21        11        21        14
5 C! w5 ~+ t/ Z% x/ p( H2 q22        10        22        10        22        12
% c; j" Y4 y" D  q0 w- M* ]& J23        11        23        14        23        13
6 d2 E* \8 |: C0 S6 W! M24        12        24        12        24        11
28        15        28        15        28        15
29        7        29        7        29        7
30        8        30        6        30        9
$ G/ ?" t- P4 U" b# P/ K8 B/ b38        19        38        17        38        1
48        5        48        9        48        8
; g; D; e" D4 S1 \8 u, _* U- T+ k62        20        62        20        62        2
表4  A区突发事件封锁方案二、三、四
4 u9 c8 t0 Y/ i+ X在对交巡台均衡工作量,加快出警时间方面,综合各节点的案发率、交巡台到其辖区内任一节点的路程进行综合评估,做出优化决策.
在案发率较高地带增设交巡台,以缓解周围交巡台的工作压力,为达均衡工作量的目的,将32号节点从7号交巡台归到8号交巡台,44号节点从2号交巡台归到3号交巡台,39号节点从2号交巡台归到16号交巡台,47号节点从6号交巡台归到7号交巡台,61号节点从7号交巡台归到4号交巡台.这样,A区每交巡台平均每天处理案件数从6.1950件、方差6.8289降到每天处理5.1917件、方差2.2182,极大的协调了工作量.对于个别节点的重新划分,会增加出警时间,但在总体上平均每天的出警时间大大缩短了.
综上考虑,共增设4个交巡台,重新分配的结果如表4所示.
7 L/ [! n+ m# F+ |5 J交巡台序号        负责区域内的节点        管辖区域内的案发率        平均每天出警时间
+ x+ v1 ~7 Q) R. T. Y3 P1        69 71 74 75 78        6.6        3.5069
' h6 L5 g6 A# d2        40 43 70 72         7.2        6.6736
" M7 i& S* l3 s$ i1 f! ]6 r3        54 55 44        5.2        4.9314
5 x0 ?& ^# M+ X& E4        57 60 61 62 63        6.4        7.7015
) g4 k4 C* n; G' A' X0 x4 b5        49 50 51 52 53 56        7.7        5.9456
/ [) J$ q! J7 [6        58 59 47        6.1        6.3949
9 M$ s/ i, M+ [9 a5 X7        30 48        5.9        3.5258
8        33 46           5.0             2.2748
& n2 f0 ]5 S' r  J3 p% h9        35 45        4.9        2.1271
10                1.6       
4 v0 O! O3 _8 i3 m* g11        26 27        4.6        2.3946
2 x: K" H# W2 X12        25        4.0        2.8622
% J6 ^' F) d- l$ K% j& L% @! T13        23 24        5.7        3.8239
14                2.5       
7 @$ s, x1 }- a: m' y' _+ P$ b15        (28) (29)        4.8        14.1580
7 l. T" i* k, p. A& j) ?0 {16        36 37 (38) (39)        6.4        10.0562
/ g& E3 r' Y  x17        41 42        5.3        2.5689
18        73 80 81 83        5.9        3.8438
9 K- k; C: I  F' V- ~# x9 A' s19        77 79        3.4        1.1457
4 Q) J( {- ]" p8 G* Z20        85 86 87 (92)        6.4        5.5354
- H8 }; G* B8 }/ D; m- C21        22        2.8        2.5239
. n  }, o( W2 h5 T) Q2 f- Z31        32 34        4.9        4.3962
7 b# y% O8 Q. n9 S66        64 65 67 68 76         5.1        2.6655
- y6 g# l* @" H; W( {0 K90        82 84 88 89 91        6.2        3.2171
& M# N: Q$ r0 u8 _: v2 l表5  优化后的A区交巡台管辖区域示意图
9 |, G5 u/ r  Q+ v" F1 A: \# n/ y
图3  A区增设平台示意图
    图中方块所示节点即为增设平台处.
5.2 问题二的模型
* c& O& D, J2 V! N) G& e( x5.2.1 模型建立
( Q3 U% E" D7 L# G& g. x: [对于问题2,对附件中所提供的A,B,C,D,E,F六城区的数据进行整合,做出直观的图表.遵循警情主导警务原则、快速出警原则、方便与安全原则,结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况,充分考虑现有警力和财力并确保安全,科学分析现有平台的数量和具体位置的合理性.
对于该市地点P(第32个节点)处发生了重大刑事案件的围堵问题,本文将其归结为资源调配问题.本文合理假设了犯罪嫌疑人的车行驶速度(分三种情况考虑:等于警车速度,警车速度的二倍,警车速度的一半),并确定三分钟后犯罪嫌疑人的车行驶的最远距离,从而利用回溯法的思想采用Matlab编程确定犯罪嫌疑人的车的所有可能位置.以时间最短、围堵区域最小为原则,采用改进的双层Floyd算法,快速地形成围堵区域,并使围堵区域尽可能的小.
5.2.2 模型求解
全市整体状况如表5所示,数据显示C区和F区的事故发生率较高、交巡警服务平台工作量高于全市平均水平且交巡警服务平台平均每天出警时间过长,针对以上问题本文再次利用问题1的Floyd算法,并考虑区域边界处的设点拥挤问题,本文提出了在C区增加5个服务平台、在F区增加1个服务平台.
% ^* I2 ~1 C& g8 o+ O9 X' q全市六个城区        城区面积        城区人口        平台数        平均人口        全区案发率        各区平台案发率均值
  s& N" w. w2 e9 N& eA        22        60        20        2.727        124.5        6.625
B        103        21        8        0.204        66.4        8.3
# K/ _3 X- Z% {: J/ ~C        221        49        17        0.223        187.2        11.012
7 W/ Y) _, S( z& f5 s. Z$ cD        383        73        9        0.191        67.8        7.533
$ y& V; }% J. ]E        432        76        15        0.176        119.4        7.96
F        274        53        11        0.193        109.2        9.927
均值                53.3333        13.3333        0.619
112.4167
' B9 i6 l9 A2 A" O" Z7 C7 w; \: `* k, D8.5595
表6  全市整体状况
% {% ?" F: b  x1 U" |' N1 v
图4  全市增设交巡台位置示意图(方块所示区域)
3 E: @" |! z: Z& ^. o( j  Y对于P点发生重大刑事案件,动用全市警力进行围堵,我们希望使得包围圈尽可能的小,由于犯罪嫌疑人的车速度未知,我们分以下三种情况进行考虑:
1)当犯罪嫌疑人的车速与警车速度同,即 .
运用穷举法,对肇事车辆可能的逃逸路线进行分析,以3分钟路程为半径,找到肇事车辆逃逸的覆盖范围,如图5所示,其中实线表示可能路径,在此范围内有8、9、10、15号共4个交巡平台,保证这4个平台警力不动,组成第一范围包围圈.

0 O7 l& ]0 L2 v' t- t图5 肇事车逃逸3分钟内覆盖区域示意图
    进一步分析可能的逃逸路线,调度16号交巡台到36号节点,2号交巡台到3号节点,3号交巡台到55号节点,6号交巡台到47号节点,组成第二组半包围,保证对A区的封锁.若肇事车辆经36号节点逃往16号节点,则会与16号交巡台在途中相遇.
$ P8 ]0 F4 B# N8 e2 h3 y1 E对于从32号节点经7号节点逃逸到30号和47号节点,存在从A区逃往其他城区的可能,需调动其他城区交巡台的支援.将C区119号交巡台调度到237号节点,将D区320号交巡台调度到371号节点,321号交巡台经368号、369号节点到370号节点,至此,在全市范围内实现全面封锁.
1 B% n  \  p* w4 A2)当犯罪嫌疑人的车速比警车车速小,即 ,我们令
    方法同1),寻找分钟逃逸范围内所覆盖的全部交巡台,经过整合分析,保持7号、8号、9号、15号共4个交巡台原地封锁,10号交巡台到34号节点封锁,6号交巡台到47号节点封锁,16号交巡台到36号节点封锁,3号交巡台经55号节点到46号节点进行封锁,2号交巡台经3号节点到45号节点封锁,在此过程中,10号、2号和3号交巡台会在途中与肇事车辆相遇.
3)当犯罪嫌疑人的车速比警车车速大,即 ,我们令
) X3 l  A+ w& T6 O% z) f7 R( g由于肇事车辆逃逸速度较快,可能会逃逸到C区和F区,故需调动C区和F区警力进行围堵.
. m" \/ k" F& s4 C   A区将20号交巡台调到62号节点,16号交巡台调到36号节点,2号交巡台经40号到39号节点,17号交巡台调到41号节点,15号、10号、4号、3号、5号、7号、8号、9号交巡台原地封锁,其余交巡台向其邻近的路口节点进行增援.经过分析,肇事车辆可能由28号、48号、30号进入C区及A、D两区的交汇地带,或由16号节点逃逸到F区,在此,对C区、D区、F区交巡台进行如下调配,实现全市封锁:
7 A3 X! B* I* j6 l+ `7 b2 _   C区:240号交巡台调到239号节点,170号交巡台调到225节点,167号交巡台调到259节点.
   D区:320号交巡台调度到371号节点,321号交巡台经368号、369号节点到370号节点.
   F区:477号交巡台调度到501号节点,518号交巡台调到521号节点,478号节点调到527号节点,484号节点到571号节点.
7.模型的评价与推广
; s+ [6 c" U- |; f3 g1 s本文避免了时间复杂度较复杂的Dijkstra算法,选用Floyd算法,在求最短路径上提高了效率,代码编写简单.
模型的建立思路清晰,遵循可操作性、科学性、可比性原则,该模型建立出了在较理想状态下交巡警平台的最优设置,减少出警时间,均衡工作量,提高工作效率,在遇突发事件时,可尽快实现道路封锁,给生活中交巡警平台的设立予参考,具有一定的实际应用价值,也可以应用于其他适用区域.模型的运算由矩阵、向量的运算组成,易于用数学软件求解和验证.
本模型较好的解决了交巡警平台的最优选址问题,当事故发生时,交巡警可以第一时间到达事发地点,有效的改善了交巡警在执行任务中的效率,在经济迅猛发展的今天,城市加速扩张,人口迅速增长,交巡警平台的设置是平安城市的最好保障.该模型也可运用到其他最优选址问题中去,比如关于消防救援工作最优路径问题、重大生产安全事故应急救援问题、公共交通的最优路径问题等. 同时也可利用该模型算法拓展模型在其他领域的适用范围.
1 r) V$ H9 i6 B  F1 n) D/ s该模型也有一定的局限性,如现实中不能时刻都保证道路的畅通性.既不能保证出警的时间总是维持在3分钟之内.忽略了实际地形对于车速的影响以及实际生活中存在的不定因素.

参考文献
[1]         [徐孝凯,王凤禄],《数据结构简明教程》第二版,北京:清华大学出版社,2005年4月1日
% V) S, _6 J+ F, R" D. Q[2]         [李建中,骆吉洲],《华章数学译丛》第二版,北京:机械工业出版社,2002年6月
7 T% @, X0 _) |: P5 v. _[3]         [陈庆华等],《组合最优化技术及其应用》第1版,北京:国防科技大学出版社,1989年8月
; r, m0 A* Z7 j4 B) Y) C[4]         [W.T.Tutte],《Graph Theory》,英国:Cambridge University Press,2001年3月1日
8.附件
' B2 G1 ~" Q: G; r7 v! @附件1:用Floyd算法分配个服务平台管辖区域
4 d6 B* E& X0 o& \  N7 M[area_a _x,area_a _y]=find (location_all _daolu<=92);
, m/ }, D$ H! d$ i road_index _a=[area_a _x,area_a _y];
road_index _a
- z" f8 t$ p* w/ L) C3 ]6 a
a1=find (road_index _a (:,2)==1);
4 Z0 H1 v+ N# d/ Y/ l0 Na2=find (road_index _a (:,2)==2);
! l" U/ l8 |; Y6 T- Q) ^7 I# [7 `A=road_index _a (a1,1);B=road_index _a (a2,1);[c,ia,ib] = intersect (A,B);
c

size (c)
/ M, V2 s* C+ G
- i+ V2 r0 d  C% P4 M( | for i=1:140
    tt=c (i);
0 F) _! n+ L8 M% M% d7 E, ?    uu=location_all _daolu (tt,;
    uu1=uu (1);uu2=uu (2);
    vv1=location_a _zuobiao (uu1,;
+ ]* l* A0 l0 o  ^5 J    vv2=location_a _zuobiao (uu2,;
    ww1=[vv1 (1),vv2 (1)];
    ww2=[vv1 (2),vv2 (2)];
    line (ww1,ww2)
end
6 w, l4 _5 |& H for i=1:140
& l) ?$ z  Y8 A: g4 ~+ Z    tt=c (i);
2 _7 |% t& p& R  Q" ]. E/ i    uu=location_all _daolu (tt,;
    uu1=uu (1);uu2=uu (2);
    vv1=location_a _zuobiao (uu1,;
, Q; Y+ `% k7 p    vv2=location_a _zuobiao (uu2,;
1 \) |: F: h& e8 ?; @9 p- Y    ww1=[vv1 (1),vv2 (1)];
    ww2=[vv1 (2),vv2 (2)];
9 f; \) G, L6 D: p+ O    line (ww1,ww2,'k')
( ]1 l/ c- t: q8 b9 R, x# f1 {1 b! {end
. M5 n* u% d  q2 ~& {  }! B% ??? Error using ==> line
% String argument is an unknown option.

/ \  P6 @& W+ X* f- s& w/ U0 `for i=1:140
    tt=c (i);
" Z" i1 s& m* h4 B- e( K    uu=location_all _daolu (tt,;
0 h" U2 F2 y! a; `. W6 a. [    uu1=uu (1);uu2=uu (2);
+ o/ h) {( c% K2 M) `3 p) o- g/ G; K    vv1=location_a _zuobiao (uu1,;
    vv2=location_a _zuobiao (uu2,;
6 |/ Y3 y, w6 Y. |8 k" ^    ww1=[vv1 (1),vv2 (1)];
    ww2=[vv1 (2),vv2 (2)];
+ a' D; C" Z* i; h  I    line (ww1,ww2,'Color',[.8 .8 .8])
end
save data_b _problem

) T6 k  w; e$ \# e% k* I% A区节点间的邻接矩阵
7 r+ i3 ^% K2 N! N) F" D2 i" vload data_b _problem;
& g% [+ R9 `& Smatric_lingjie=zeros (92,92);
[xx,yy]=find (matric_lingjie==0);matric_lingjie (xx,yy)=inf;
" J0 F# ]- j3 r4 G" sfor i=1:92
        matric_lingjie (i,i)=0;
3 e/ [7 j3 N# f5 a: z7 `/ Kend
# Z; u# q2 [( ]+ G& I

. I8 y5 B  m. H! q" v for i=1:140
4 P8 U) X+ A! m2 @8 H3 j    tt=c (i);
    uu=location_all _daolu (tt,;
6 w7 G/ \3 _( M- h    uu1=uu (1);uu2=uu (2);                       % 端点序号
, \; p* ^5 h0 |, S/ D: O% T- j    vv1=location_a _zuobiao (uu1,;           % 第一个端点坐标
    vv2=location_a _zuobiao (uu2,;           % 第二个端点坐标
4 {. c+ Z/ L$ d6 T4 c    % 计算端点间距离
' f/ i; f; }: O# C2 p    distance=sqrt ((vv1 (1)-vv2 (1))^2+(vv1 (2)-vv2 (2))^2);
$ W7 m% I' H; Z. n0 a    matric_lingjie (uu1,uu2)=distance;
$ ^" g4 n# X, a3 {    matric_lingjie (uu2,uu1)=distance;        % 赋值给邻接矩阵
end

3 _% c1 L2 `$ w- ~[D,R]=floyd (matric_lingjie);
; N( a: M. s0 w! P' {. V
2 \8 d  H' ]7 B/ l8 T' {7 nmatric_fenkuai=D (1:20,;
) ^+ H* l5 {" g5 Y- N
  w2 [( [- j8 S. x( g  _4 Cfor i=1:20
% x" k+ l6 ~! ]* h    for j=1:92
        if matric_fenkuai (i,j)>30
6 ~6 m& h' F* o! @" O' k            matric_fenkuai (i,j)=0;
$ Q1 J; |* O: k- u        end
1 q" B/ x' b5 z    end
end

ti=zeros (1,92);
ti (1)=text (location_a _zuobiao (1,1),location_a _zuobiao (1,2)+1.5,'1');
ti (2)=text (location_a _zuobiao (2,1),location_a _zuobiao (2,2)+1.5,'2');
7 w; E1 b' W% ~; {2 ^0 X) cti (3)=text (location_a _zuobiao (3,1),location_a _zuobiao (3,2)+1.5,'3');
5 o. _* \  L+ d3 y1 r8 Q. Hti (4)=text (location_a _zuobiao (4,1),location_a _zuobiao (4,2)+1.5,'4');
ti (5)=text (location_a _zuobiao (5,1),location_a _zuobiao (5,2)+1.5,'5');
ti (6)=text (location_a _zuobiao (6,1),location_a _zuobiao (6,2)+1.5,'6');
ti (7)=text (location_a _zuobiao (7,1),location_a _zuobiao (7,2)+1.5,'7');
! ^4 S8 u' W5 E# k: C% a3 Q/ F8 Sti (8)=text (location_a _zuobiao (8,1),location_a _zuobiao (8,2)+1.5,'8');
6 v6 s& p0 A6 m: B. l; G9 ?% A0 A" |ti (9)=text (location_a _zuobiao (9,1),location_a _zuobiao (9,2)+1.5,'9');
ti (10)=text (location_a _zuobiao (10,1),location_a _zuobiao (10,2)+1.5,'10');
ti (11)=text (location_a _zuobiao (11,1),location_a _zuobiao (11,2)+1.5,'11');
; u/ a% S! E/ B/ k. D9 s2 i/ }ti (12)=text (location_a _zuobiao (12,1),location_a _zuobiao (12,2)+1.5,'12');
+ J, @% N6 ~: J# T! A2 L  gti (13)=text (location_a _zuobiao (13,1),location_a _zuobiao (13,2)+1.5,'13');
ti (14)=text (location_a _zuobiao (14,1),location_a _zuobiao (14,2)+1.5,'14');
ti (15)=text (location_a _zuobiao (15,1),location_a _zuobiao (15,2)+1.5,'15');
ti (16)=text (location_a _zuobiao (16,1),location_a _zuobiao (16,2)+1.5,'16');
ti (17)=text (location_a _zuobiao (17,1),location_a _zuobiao (17,2)+1.5,'17');
ti (18)=text (location_a _zuobiao (18,1),location_a _zuobiao (18,2)+1.5,'18');
ti (19)=text (location_a _zuobiao (19,1),location_a _zuobiao (19,2)+1.5,'19');
ti (20)=text (location_a _zuobiao (20,1),location_a _zuobiao (20,2)+1.5,'20');
ti (21)=text (location_a _zuobiao (21,1),location_a _zuobiao (21,2)+1.5,'21');
+ A" w, W- p6 Z1 Zti (22)=text (location_a _zuobiao (22,1),location_a _zuobiao (22,2)+1.5,'22');
3 u9 O" ?5 j/ t: }ti (23)=text (location_a _zuobiao (23,1),location_a _zuobiao (23,2)+1.5,'23');
* h: |4 k5 ~7 z/ U6 y! b8 ]8 Dti (24)=text (location_a _zuobiao (24,1),location_a _zuobiao (24,2)+1.5,'24');
ti (25)=text (location_a _zuobiao (25,1),location_a _zuobiao (25,2)+1.5,'25');
9 E2 E* Z% H7 g% C0 C2 Qti (26)=text (location_a _zuobiao (26,1),location_a _zuobiao (26,2)+1.5,'26');
+ s9 l- k# N: U. O- s& A9 @8 B* nti (27)=text (location_a _zuobiao (27,1),location_a _zuobiao (27,2)+1.5,'27');
ti (28)=text (location_a _zuobiao (28,1),location_a _zuobiao (28,2)+1.5,'28');
ti (29)=text (location_a _zuobiao (29,1),location_a _zuobiao (29,2)+1.5,'29');
. j: L0 T) R, E4 }- E, uti (30)=text (location_a _zuobiao (30,1),location_a _zuobiao (30,2)+1.5,'30');
3 @9 {4 f5 w% H$ K6 V5 Oti (31)=text (location_a _zuobiao (31,1),location_a _zuobiao (31,2)+1.5,'31');
ti (32)=text (location_a _zuobiao (32,1),location_a _zuobiao (32,2)+1.5,'32');
& V6 j3 F$ f" A* H- C1 `2 jti (33)=text (location_a _zuobiao (33,1),location_a _zuobiao (33,2)+1.5,'33');
ti (34)=text (location_a _zuobiao (34,1),location_a _zuobiao (34,2)+1.5,'34');
2 T. |" e1 E& w1 Jti (35)=text (location_a _zuobiao (35,1),location_a _zuobiao (35,2)+1.5,'35');
: M* q. a% o" y, O, s! G# qti (36)=text (location_a _zuobiao (36,1),location_a _zuobiao (36,2)+1.5,'36');
ti (37)=text (location_a _zuobiao (37,1),location_a _zuobiao (37,2)+1.5,'37');
+ l$ m0 h2 A7 ^5 tti (38)=text (location_a _zuobiao (38,1),location_a _zuobiao (38,2)+1.5,'38');
$ l, e4 A% ?( P( ]& k- O* Zti (39)=text (location_a _zuobiao (39,1),location_a _zuobiao (39,2)+1.5,'39');
ti (40)=text (location_a _zuobiao (40,1),location_a _zuobiao (40,2)+1.5,'40');
ti (41)=text (location_a _zuobiao (41,1),location_a _zuobiao (41,2)+1.5,'41');
. u7 _; O( f1 Z3 Qti (42)=text (location_a _zuobiao (42,1),location_a _zuobiao (42,2)+1.5,'42');
ti (43)=text (location_a _zuobiao (43,1),location_a _zuobiao (43,2)+1.5,'43');
3 D) U$ B4 ~+ E8 y3 d( ^) H& k8 hti (44)=text (location_a _zuobiao (44,1),location_a _zuobiao (44,2)+1.5,'44');
# g' q7 F; E4 `6 V& L' zti (45)=text (location_a _zuobiao (45,1),location_a _zuobiao (45,2)+1.5,'45');
: `& ~' j4 w7 l- K" n% X4 jti (46)=text (location_a _zuobiao (46,1),location_a _zuobiao (46,2)+1.5,'46');
& d7 o4 v- O. nti (47)=text (location_a _zuobiao (47,1),location_a _zuobiao (47,2)+1.5,'47');
ti (48)=text (location_a _zuobiao (48,1),location_a _zuobiao (48,2)+1.5,'48');
$ v! U" |9 [$ `0 W( b7 R% |1 Mti (49)=text (location_a _zuobiao (49,1),location_a _zuobiao (49,2)+1.5,'49');
$ p) J" z! S0 |3 p* dti (50)=text (location_a _zuobiao (50,1),location_a _zuobiao (50,2)+1.5,'50');
2 ]) w* r( ~( W. H1 i; X+ W8 oti (51)=text (location_a _zuobiao (51,1),location_a _zuobiao (51,2)+1.5,'51');
ti (52)=text (location_a _zuobiao (52,1),location_a _zuobiao (52,2)+1.5,'52');
ti (53)=text (location_a _zuobiao (53,1),location_a _zuobiao (53,2)+1.5,'53');
) q0 @0 h2 f; V- e& Iti (54)=text (location_a _zuobiao (54,1),location_a _zuobiao (54,2)+1.5,'54');
ti (55)=text (location_a _zuobiao (55,1),location_a _zuobiao (55,2)+1.5,'55');
. v" C% k9 p% @( O! ^ti (56)=text (location_a _zuobiao (56,1),location_a _zuobiao (56,2)+1.5,'56');
2 o; S8 Z2 _8 |# Y' t9 C5 ati (57)=text (location_a _zuobiao (57,1),location_a _zuobiao (57,2)+1.5,'57');
+ n8 O- `3 M" T% ]5 [3 q. M$ V% {ti (58)=text (location_a _zuobiao (58,1),location_a _zuobiao (58,2)+1.5,'58');
ti (59)=text (location_a _zuobiao (59,1),location_a _zuobiao (59,2)+1.5,'59');
- w6 e( e7 w  I+ ]+ gti (60)=text (location_a _zuobiao (60,1),location_a _zuobiao (60,2)+1.5,'60');
ti (61)=text (location_a _zuobiao (61,1),location_a _zuobiao (61,2)+1.5,'61');
- u) r0 Z3 i( z: {3 K2 H% ^/ tti (62)=text (location_a _zuobiao (62,1),location_a _zuobiao (62,2)+1.5,'62');
' g4 r. H2 L9 [: `# {9 Rti (63)=text (location_a _zuobiao (63,1),location_a _zuobiao (63,2)+1.5,'63');
ti (64)=text (location_a _zuobiao (64,1),location_a _zuobiao (64,2)+1.5,'64');
ti (65)=text (location_a _zuobiao (65,1),location_a _zuobiao (65,2)+1.5,'65');
ti (66)=text (location_a _zuobiao (66,1),location_a _zuobiao (66,2)+1.5,'66');
ti (67)=text (location_a _zuobiao (67,1),location_a _zuobiao (67,2)+1.5,'67');
ti (68)=text (location_a _zuobiao (68,1),location_a _zuobiao (68,2)+1.5,'68');
. q+ ^4 g! H/ @4 R+ Q3 H- Lti (69)=text (location_a _zuobiao (69,1),location_a _zuobiao (69,2)+1.5,'69');
6 G; _2 s3 X% ]ti (70)=text (location_a _zuobiao (70,1),location_a _zuobiao (70,2)+1.5,'70');
  y6 v, L- O1 k- h) r6 Q' b, vti (71)=text (location_a _zuobiao (71,1),location_a _zuobiao (71,2)+1.5,'71');
% U6 W% D5 ^9 A- @ti (72)=text (location_a _zuobiao (72,1),location_a _zuobiao (72,2)+1.5,'72');
! E9 ?7 E1 S- g( ^: xti (73)=text (location_a _zuobiao (73,1),location_a _zuobiao (73,2)+1.5,'73');
ti (74)=text (location_a _zuobiao (74,1),location_a _zuobiao (74,2)+1.5,'74');
ti (75)=text (location_a _zuobiao (75,1),location_a _zuobiao (75,2)+1.5,'75');
$ Y0 R0 P+ [! x9 jti (76)=text (location_a _zuobiao (76,1),location_a _zuobiao (76,2)+1.5,'76');
1 ?- W( v# a. @% C- D/ d4 Dti (77)=text (location_a _zuobiao (77,1),location_a _zuobiao (77,2)+1.5,'77');
ti (78)=text (location_a _zuobiao (78,1),location_a _zuobiao (78,2)+1.5,'78');
* `7 S4 X; `" y$ @* bti (79)=text (location_a _zuobiao (79,1),location_a _zuobiao (79,2)+1.5,'79');
. K; p/ E, @( l9 P! Cti (80)=text (location_a _zuobiao (80,1),location_a _zuobiao (80,2)+1.5,'80');
) I7 E6 a$ O! I# k% |8 \- Eti (81)=text (location_a _zuobiao (81,1),location_a _zuobiao (81,2)+1.5,'81');
( \4 z- H& ^. c9 Y9 Cti (82)=text (location_a _zuobiao (82,1),location_a _zuobiao (82,2)+1.5,'82');
ti (83)=text (location_a _zuobiao (83,1),location_a _zuobiao (83,2)+1.5,'83');
ti (84)=text (location_a _zuobiao (84,1),location_a _zuobiao (84,2)+1.5,'84');
8 F2 g- F9 w8 ^% U9 ~; uti (85)=text (location_a _zuobiao (85,1),location_a _zuobiao (85,2)+1.5,'85');
+ k3 f) ?1 v% Fti (86)=text (location_a _zuobiao (86,1),location_a _zuobiao (86,2)+1.5,'86');
ti (87)=text (location_a _zuobiao (87,1),location_a _zuobiao (87,2)+1.5,'87');
0 V" w& E3 E( Y% }" s7 Uti (88)=text (location_a _zuobiao (88,1),location_a _zuobiao (88,2)+1.5,'88');
ti (89)=text (location_a _zuobiao (89,1),location_a _zuobiao (89,2)+1.5,'89');
ti (90)=text (location_a _zuobiao (90,1),location_a _zuobiao (90,2)+1.5,'90');
  Q! k1 F: m. |+ [ti (91)=text (location_a _zuobiao (91,1),location_a _zuobiao (91,2)+1.5,'91');
' G1 [5 c  L$ A- [0 q7 o% R6 f6 Xti (92)=text (location_a _zuobiao (92,1),location_a _zuobiao (92,2)+1.5,'92');
8 x( O" f+ A, ?+ K: h$ j
% ?- o* x7 ]! a3 f, B3 m
fenpei=[1        1
1        67
1        68
1        69
, l+ n. P- R7 X9 \1        71
1        73
- m( J7 L3 y4 U& z$ z1        74
1        75
- c$ P- F! Y# f8 D* A1 {1        76
1        78
2        2
2        39
2 B' J' Y( J9 v: E+ I! h, k2        40
9 d6 n. l( c" h; ~% {2 c* m2        43
; u4 w9 H8 W% c# d+ {* j; L2        44
4 d) @4 y7 z& M! Q2        70
: q* a2 M! H2 u) f1 j! {5 Q1 R2        72
! `- _  \0 p& r  [9 x4 n/ h3        3
3        54
3        55
3        65
3        66
4        4
4        57
4        60
4        62
4        63
& R/ R/ k7 j0 X0 V  ^4        64
5        5
5        49
5        50
8 ?, r2 @! \, W2 u" B/ L5        51
5        52
  l3 \2 F3 q" }4 [, x. g2 t$ F5        53
/ J, d  f$ O2 _) B5        56
5        58
5        59
6        6
7        7
7        30
7        32
7        47
7        48
. t4 ~$ R. ~9 F2 l7 A* q  a6 y4 {7        61
8        8
8        33
  N7 D% J4 ^& T& D0 @8 U( O8        46
! ^9 p' a. v% p9        9
9        31
9        34
9        35
) V; j* z0 q" z8 c8 M9        45
0 t( c, Y( l( _: G10        10
3 m, H$ J, z7 m4 r. ~11        11
7 D+ i  x6 a% s) S  K- p11        26
11        27
8 W5 Z( |% ]- v: k# e7 \12        12
12        25
13        13
' @4 z3 B3 D" f  E% h13        21
13        22
: j5 e: J; C0 I) b13        23
13        24
14        14
15        15
- p# [& M2 w! |3 U$ j, T) V15        28
: u/ Y: A+ g( V3 g1 k2 ~# Z6 ~# ?15        29
( S- h  x# L  o. H( P8 ^16        16
16        36
16        37
16        38
17        17
17        41
17        42
18        18
18        80
# U, y: l. L, \% `. u18        81
18        82
18        83
/ o# j5 l  T/ m, a' Y7 i1 j19        19
6 g- ~! F, {; c% E19        77
% B8 E& s' d6 [19        79
20        20
20        84
20        85
# L* S" D  Y* c, s  _- z20        86
( v7 g4 `9 i6 @) f' R20        87
* \. g5 D. y0 [3 I: n  {20        88
20        89
& v8 l/ R5 O8 V. J, G20        90
20        91
20        92];

for i=1:92
( J1 @, N$ ~" y+ a3 s: L    if fenpei (i,1)==1
' _3 H3 H/ a8 P' [7 Y        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b');hold on;
7 y/ w6 Z5 k! _' L8 p; Q    elseif fenpei (i,1)==2
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'g');hold on;
. e6 S% g* h- L' G    elseif fenpei (i,1)==3
# }$ m, z4 r& t        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'r');hold on;
$ a9 a, ?! X" U    elseif fenpei (i,1)==4
* A6 T8 d' o1 l  b* D8 [        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'c');hold on;
    elseif fenpei (i,1)==5
4 ~8 N1 K: T8 w6 |        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'m');hold on;
) T3 ~. {7 |) m4 |# R$ z    elseif fenpei (i,1)==6
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'y');hold on;
    elseif fenpei (i,1)==7
3 R9 }) _  B9 k  Z- d5 X        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'k');hold on;
3 |  ?$ Z% V7 ~' e    elseif fenpei (i,1)==8
$ N3 m) y6 w: q7 K# n$ L        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b+');hold on;
$ K- j. N, M. [* A$ R5 f5 B/ Z    elseif fenpei (i,1)==9
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bo');hold on;
    elseif fenpei (i,1)==10
6 D# q$ N8 v' R, t  g" N        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bs');hold on;
    elseif fenpei (i,1)==11
1 x$ R' Y  U2 g        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'g+');hold on;
, w/ @' W& {7 L    elseif fenpei (i,1)==12
$ c/ A( P) D7 `/ u/ }        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'go');hold on;
    elseif fenpei (i,1)==13
# Y& `7 j8 L  A% m        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'gs');hold on;
    elseif fenpei (i,1)==14
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'r+');hold on;
& n) S3 g2 v3 R; \1 z    elseif fenpei (i,1)==15
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'ro');hold on;
  P9 R, _0 \6 }+ ]. N& v# o+ o1 X- O    elseif fenpei (i,1)==16
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'rs');hold on;
    elseif fenpei (i,1)==17
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'c+');hold on;
    elseif fenpei (i,1)==18
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'co');hold on;
9 I) g0 A/ G: ?  U- O/ H  i& t! N- {    elseif fenpei (i,1)==19
2 g# `, m! J& u        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'cs');hold on;
    else
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'m+');hold on;
    end
- s  E0 l' u5 I3 mend

# k% [8 j, Y" |# C* Xfor i=1:92
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b');hold on;
end
axis ([190,460,250,410]);
2 X. y6 E3 o" `3 l6 x" Hfor i=1:140
    tt=c (i);
    uu=location_all _daolu (tt,;
    uu1=uu (1);uu2=uu (2);
1 S5 ]. l! }. x% L1 }    vv1=location_a _zuobiao (uu1,;
& m& d) i2 x! ]+ {( D& y: b    vv2=location_a _zuobiao (uu2,;
& q5 e: n7 C& W8 `    ww1=[vv1 (1),vv2 (1)];
; L: @# h" z- ^2 N! r6 N% q6 V    ww2=[vv1 (2),vv2 (2)];
    line (ww1,ww2,'Color',[.5 .5 .5])
end
$ H; n; J+ \+ K( y. e
ti=zeros (1,92);
ti (1)=text (location_a _zuobiao (1,1),location_a _zuobiao (1,2)+1.5,'1');
ti (2)=text (location_a _zuobiao (2,1),location_a _zuobiao (2,2)+1.5,'2');
. s+ |# x* w9 s) Wti (3)=text (location_a _zuobiao (3,1),location_a _zuobiao (3,2)+1.5,'3');
$ M+ ~) A, O) t6 y- n0 P" X' _; mti (4)=text (location_a _zuobiao (4,1),location_a _zuobiao (4,2)+1.5,'4');
. @" V5 C* ^  p6 ~6 D3 f& ]5 Uti (5)=text (location_a _zuobiao (5,1),location_a _zuobiao (5,2)+1.5,'5');
; c$ x% ?+ [* z5 k& v2 n+ Eti (6)=text (location_a _zuobiao (6,1),location_a _zuobiao (6,2)+1.5,'6');
: E, T9 M$ v( h6 }3 xti (7)=text (location_a _zuobiao (7,1),location_a _zuobiao (7,2)+1.5,'7');
8 Q% D2 I( X/ e* Q5 ~ti (8)=text (location_a _zuobiao (8,1),location_a _zuobiao (8,2)+1.5,'8');
ti (9)=text (location_a _zuobiao (9,1),location_a _zuobiao (9,2)+1.5,'9');
, ?9 r8 F' E: R2 A. }( Oti (10)=text (location_a _zuobiao (10,1),location_a _zuobiao (10,2)+1.5,'10');
5 M" f! j8 v# S2 x4 B3 _0 n% ?ti (11)=text (location_a _zuobiao (11,1),location_a _zuobiao (11,2)+1.5,'11');
ti (12)=text (location_a _zuobiao (12,1),location_a _zuobiao (12,2)+1.5,'12');
+ u/ S6 K: p' Wti (13)=text (location_a _zuobiao (13,1),location_a _zuobiao (13,2)+1.5,'13');
ti (14)=text (location_a _zuobiao (14,1),location_a _zuobiao (14,2)+1.5,'14');
ti (15)=text (location_a _zuobiao (15,1),location_a _zuobiao (15,2)+1.5,'15');
ti (16)=text (location_a _zuobiao (16,1),location_a _zuobiao (16,2)+1.5,'16');
ti (17)=text (location_a _zuobiao (17,1),location_a _zuobiao (17,2)+1.5,'17');
! {! S9 a' Q/ j0 M& e3 v* Cti (18)=text (location_a _zuobiao (18,1),location_a _zuobiao (18,2)+1.5,'18');
; M  a* [0 _# J7 F5 C4 X' h7 U* [& Uti (19)=text (location_a _zuobiao (19,1),location_a _zuobiao (19,2)+1.5,'19');
ti (20)=text (location_a _zuobiao (20,1),location_a _zuobiao (20,2)+1.5,'20');
ti (21)=text (location_a _zuobiao (21,1),location_a _zuobiao (21,2)+1.5,'21');
1 P9 ]0 T( j* \9 o6 y& d" ~% J; bti (22)=text (location_a _zuobiao (22,1),location_a _zuobiao (22,2)+1.5,'22');
- |0 Z- X. C! Z' X/ Iti (23)=text (location_a _zuobiao (23,1),location_a _zuobiao (23,2)+1.5,'23');
" a% T6 J1 e  R2 M% u* wti (24)=text (location_a _zuobiao (24,1),location_a _zuobiao (24,2)+1.5,'24');
ti (25)=text (location_a _zuobiao (25,1),location_a _zuobiao (25,2)+1.5,'25');
ti (26)=text (location_a _zuobiao (26,1),location_a _zuobiao (26,2)+1.5,'26');
ti (27)=text (location_a _zuobiao (27,1),location_a _zuobiao (27,2)+1.5,'27');
% J4 Q! @& h" a! |8 c& qti (28)=text (location_a _zuobiao (28,1),location_a _zuobiao (28,2)+1.5,'28');
8 B9 e% X0 [! i% \; y2 xti (29)=text (location_a _zuobiao (29,1),location_a _zuobiao (29,2)+1.5,'29');
ti (30)=text (location_a _zuobiao (30,1),location_a _zuobiao (30,2)+1.5,'30');
' L& X" }3 t$ X9 ^1 ?! C# o, gti (31)=text (location_a _zuobiao (31,1),location_a _zuobiao (31,2)+1.5,'31');
ti (32)=text (location_a _zuobiao (32,1),location_a _zuobiao (32,2)+1.5,'32');
% E! S4 \1 k2 eti (33)=text (location_a _zuobiao (33,1),location_a _zuobiao (33,2)+1.5,'33');
ti (34)=text (location_a _zuobiao (34,1),location_a _zuobiao (34,2)+1.5,'34');
ti (35)=text (location_a _zuobiao (35,1),location_a _zuobiao (35,2)+1.5,'35');
ti (36)=text (location_a _zuobiao (36,1),location_a _zuobiao (36,2)+1.5,'36');
8 o  Y2 y+ l) s5 ?ti (37)=text (location_a _zuobiao (37,1),location_a _zuobiao (37,2)+1.5,'37');
% a  K1 p" {/ g$ p6 E: dti (38)=text (location_a _zuobiao (38,1),location_a _zuobiao (38,2)+1.5,'38');
6 N) Y6 X2 `4 Q# i2 dti (39)=text (location_a _zuobiao (39,1),location_a _zuobiao (39,2)+1.5,'39');
+ X3 S7 c# i0 Dti (40)=text (location_a _zuobiao (40,1),location_a _zuobiao (40,2)+1.5,'40');
- z- O) I( v/ v5 V  v4 hti (41)=text (location_a _zuobiao (41,1),location_a _zuobiao (41,2)+1.5,'41');
ti (42)=text (location_a _zuobiao (42,1),location_a _zuobiao (42,2)+1.5,'42');
ti (43)=text (location_a _zuobiao (43,1),location_a _zuobiao (43,2)+1.5,'43');
ti (44)=text (location_a _zuobiao (44,1),location_a _zuobiao (44,2)+1.5,'44');
( j/ P) C8 C5 i# H2 @ti (45)=text (location_a _zuobiao (45,1),location_a _zuobiao (45,2)+1.5,'45');
ti (46)=text (location_a _zuobiao (46,1),location_a _zuobiao (46,2)+1.5,'46');
1 ]6 J# |* o/ ?. b- V+ ]) |  X" Qti (47)=text (location_a _zuobiao (47,1),location_a _zuobiao (47,2)+1.5,'47');
ti (48)=text (location_a _zuobiao (48,1),location_a _zuobiao (48,2)+1.5,'48');
$ H; s) X" ]5 c1 L3 ^ti (49)=text (location_a _zuobiao (49,1),location_a _zuobiao (49,2)+1.5,'49');
- O+ U) o, }& @2 I" m- ^ti (50)=text (location_a _zuobiao (50,1),location_a _zuobiao (50,2)+1.5,'50');
; n. f. Z7 [: j) y" }: iti (51)=text (location_a _zuobiao (51,1),location_a _zuobiao (51,2)+1.5,'51');
6 X: X1 F3 C" ?; B: T: k: Eti (52)=text (location_a _zuobiao (52,1),location_a _zuobiao (52,2)+1.5,'52');
ti (53)=text (location_a _zuobiao (53,1),location_a _zuobiao (53,2)+1.5,'53');
, ^& y# w3 t. O! K9 z6 d/ F) mti (54)=text (location_a _zuobiao (54,1),location_a _zuobiao (54,2)+1.5,'54');
ti (55)=text (location_a _zuobiao (55,1),location_a _zuobiao (55,2)+1.5,'55');
ti (56)=text (location_a _zuobiao (56,1),location_a _zuobiao (56,2)+1.5,'56');
ti (57)=text (location_a _zuobiao (57,1),location_a _zuobiao (57,2)+1.5,'57');
ti (58)=text (location_a _zuobiao (58,1),location_a _zuobiao (58,2)+1.5,'58');
ti (59)=text (location_a _zuobiao (59,1),location_a _zuobiao (59,2)+1.5,'59');
2 w" S4 j/ @/ F+ Y2 {ti (60)=text (location_a _zuobiao (60,1),location_a _zuobiao (60,2)+1.5,'60');
ti (61)=text (location_a _zuobiao (61,1),location_a _zuobiao (61,2)+1.5,'61');
( j& G4 x0 ~; z5 K5 R# gti (62)=text (location_a _zuobiao (62,1),location_a _zuobiao (62,2)+1.5,'62');
ti (63)=text (location_a _zuobiao (63,1),location_a _zuobiao (63,2)+1.5,'63');
% m4 g  U% _2 Rti (64)=text (location_a _zuobiao (64,1),location_a _zuobiao (64,2)+1.5,'64');
: z8 z& \; r4 t6 t0 Q" Hti (65)=text (location_a _zuobiao (65,1),location_a _zuobiao (65,2)+1.5,'65');
ti (66)=text (location_a _zuobiao (66,1),location_a _zuobiao (66,2)+1.5,'66');
ti (67)=text (location_a _zuobiao (67,1),location_a _zuobiao (67,2)+1.5,'67');
4 ^( G: d4 g" h& {3 Vti (68)=text (location_a _zuobiao (68,1),location_a _zuobiao (68,2)+1.5,'68');
ti (69)=text (location_a _zuobiao (69,1),location_a _zuobiao (69,2)+1.5,'69');
ti (70)=text (location_a _zuobiao (70,1),location_a _zuobiao (70,2)+1.5,'70');
ti (71)=text (location_a _zuobiao (71,1),location_a _zuobiao (71,2)+1.5,'71');
ti (72)=text (location_a _zuobiao (72,1),location_a _zuobiao (72,2)+1.5,'72');
6 H/ j5 N$ B, P$ t1 y& F0 |ti (73)=text (location_a _zuobiao (73,1),location_a _zuobiao (73,2)+1.5,'73');
' ^: Q2 A. J# a" Y# Oti (74)=text (location_a _zuobiao (74,1),location_a _zuobiao (74,2)+1.5,'74');
" _) h6 j6 {9 E+ R& m+ Z: Zti (75)=text (location_a _zuobiao (75,1),location_a _zuobiao (75,2)+1.5,'75');
3 L! c8 V+ j, q! P: ^5 q( c$ Iti (76)=text (location_a _zuobiao (76,1),location_a _zuobiao (76,2)+1.5,'76');
& _' j* F. o9 t8 Nti (77)=text (location_a _zuobiao (77,1),location_a _zuobiao (77,2)+1.5,'77');
1 [, O% C) ^/ v- h) b) B' V5 Ati (78)=text (location_a _zuobiao (78,1),location_a _zuobiao (78,2)+1.5,'78');
ti (79)=text (location_a _zuobiao (79,1),location_a _zuobiao (79,2)+1.5,'79');
ti (80)=text (location_a _zuobiao (80,1),location_a _zuobiao (80,2)+1.5,'80');
5 O" C8 s& D0 ~* u0 Lti (81)=text (location_a _zuobiao (81,1),location_a _zuobiao (81,2)+1.5,'81');
0 f# U" F4 q" x0 C; j4 A; L$ Zti (82)=text (location_a _zuobiao (82,1),location_a _zuobiao (82,2)+1.5,'82');
ti (83)=text (location_a _zuobiao (83,1),location_a _zuobiao (83,2)+1.5,'83');
ti (84)=text (location_a _zuobiao (84,1),location_a _zuobiao (84,2)+1.5,'84');
ti (85)=text (location_a _zuobiao (85,1),location_a _zuobiao (85,2)+1.5,'85');
' z6 m/ n  w% x) J0 G) m' R6 iti (86)=text (location_a _zuobiao (86,1),location_a _zuobiao (86,2)+1.5,'86');
ti (87)=text (location_a _zuobiao (87,1),location_a _zuobiao (87,2)+1.5,'87');
ti (88)=text (location_a _zuobiao (88,1),location_a _zuobiao (88,2)+1.5,'88');
. s% ^* [0 h( v  fti (89)=text (location_a _zuobiao (89,1),location_a _zuobiao (89,2)+1.5,'89');
2 E4 J$ D# O( p1 W! jti (90)=text (location_a _zuobiao (90,1),location_a _zuobiao (90,2)+1.5,'90');
4 L8 \8 g/ ]0 |- t. l7 |- _. Xti (91)=text (location_a _zuobiao (91,1),location_a _zuobiao (91,2)+1.5,'91');
ti (92)=text (location_a _zuobiao (92,1),location_a _zuobiao (92,2)+1.5,'92');
   
help voronoi
[vx,vy] = voronoi (location_a _zuobiao _x,location_a _zuobiao _y);

for i=1:92
. W! k9 \. i* G$ X$ ~) q+ y    if fenpei (i,1)==1
; i" {! P3 P, D, q$ C" I# D        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bh');hold on;
    elseif fenpei (i,1)==2
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'gh');hold on;
  c$ j( A. n6 w& b$ s* M5 X# N% A    elseif fenpei (i,1)==3
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'rh');hold on;
    elseif fenpei (i,1)==4
1 O4 y, B$ ?) }# q# V5 c: E! m        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'ch');hold on;
    elseif fenpei (i,1)==5
: E  j- H) Z- y( n" ^: f9 B1 I$ l        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'mh');hold on;
    elseif fenpei (i,1)==6
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'yh');hold on;
: z' U& M* r$ a; `$ s    elseif fenpei (i,1)==7
' Z1 K2 e! m% Y& N9 t+ g        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'kh');hold on;
1 l( h3 r6 @5 {    elseif fenpei (i,1)==8
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b+');hold on;
    elseif fenpei (i,1)==9
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bo');hold on;
& T5 l, v3 _; u# ~    elseif fenpei (i,1)==10
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bs');hold on;
. h$ p4 [# X/ m% C, j    elseif fenpei (i,1)==11
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'g+');hold on;
5 e* i9 h& Z' F; f5 q# W; p) `    elseif fenpei (i,1)==12
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'go');hold on;
& C& A) h8 Y& T* ~# ^    elseif fenpei (i,1)==13
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'gs');hold on;
    elseif fenpei (i,1)==14
/ j+ w6 o/ d( C9 b* K        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'r+');hold on;
    elseif fenpei (i,1)==15
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'ro');hold on;
! v# R8 a9 `/ m( r4 Y  d+ j! B    elseif fenpei (i,1)==16
! {0 Q2 I3 k- }3 j" r: S9 ?        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'rs');hold on;
. H7 _& J7 h# b( a# A, \    elseif fenpei (i,1)==17
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'c+');hold on;
    elseif fenpei (i,1)==18
" H* P& I% L% V  O    plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'co');hold on;
& G- a) M  Z/ n: k& X6 i! V8 b6 l5 E    elseif fenpei (i,1)==19
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'cs');hold on;
    else
        plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'m+');hold on;
    end
end
附件2:邻接矩阵的matlab实现程序
data=[
1   75
, r! I4 P1 d- U9 B' ^2 V5 B6 y, {1   78
2   44
! H/ I' D& o5 ?4 i1 l7 {2 x/ s0 U3   45
3   65
4   39
4   63
5   49
5 f' w7 I$ E8 \  e5   50
0 ]  W/ M0 h. @. [6   59
7   32
7   47
8   9
8   47
9   35
10  34
% c' e! Z( R2 Z+ }* B: K11  22
11  26
12  25
12  471
6 T/ B- D, K: A+ N. q14  21
5 x8 Q& q6 O4 C+ W: ?) m15  7
. W2 y* K$ i3 d# Y$ D( y* K  }' \15  31
16  14
16  38
17  40
2 v, i' m. l4 b) a' ~' F" h2 M3 _17  42
17  81
18  81
18  83
19  79
20  86
21  22
22  372
22  13
7 l. u" }* M5 C% L2 a  j+ |5 z23  13
/ y5 E8 N7 r- \% ?7 r/ F6 e2 e23  383
2 v$ v3 V' m4 z/ v+ L! @0 w24  13
24  25
5 }$ l! M% q% p6 n25  11
26  27
26  10
1 `: p0 G) z6 M+ x1 U27  12
8 S# x; l: {% f. e6 W( {0 e28  29
' v/ I5 ]7 D. w/ J0 y) A28  15
29  30
) e- {, S' R! N4 P9 k30  7
* v! f( ]5 R6 G: o+ Z: Z) v30  48
31  32
5 n) o' k7 H! P31  34
32  33
33  34
33  8
- g/ A0 G& w4 A& a2 j34  9
35  45
36  35
36  37
9 y4 c$ w% v6 L' l  G4 {36  16
. N* U- v( t" x' f: T! Z/ I- E/ N3 A36  39
. g$ ]8 d) t  q37  7
38  39
38  41
0 p+ \2 W0 y  h3 O2 H; I' x. d. g: d; k39  40
6 J& |* Q0 A, W$ e40  2
41  17
41  92
4 f) z: S0 T# q  x42  43
43  2
43  72
" i6 H% [! i8 L/ D/ d" g5 d9 p$ Q44  3
45  46
46  8
9 D$ p, A; s% @- p9 r* e9 B- l46  55
47  48
) {. M7 e3 e: X: X- G- ?47  6
4 E% A! s* ?$ C- n! _47  5
: m5 t1 u8 J" C5 A1 e! F7 d. ^. p48  61
# {* `3 r9 _. y/ H4 z2 ^5 b6 P49  50
49  53
6 ]- J7 z# z; ~( [. q  g/ ~3 t2 P50  51
51  52
51  59
' h3 t( M7 G  A52  56
53  52
" O) ]! A3 k4 v' A2 P53  54
54  55
  v% G: m6 D+ @$ u/ s54  63
55  3
2 S3 N3 ^5 J# L: U* }# _3 {56  57
57  58
7 O+ P/ b* s" ]$ f0 G57  60
/ Z: K# D  \7 K  y+ N57  4
58  59
  k4 i+ V' I' b) t* t, e60  62
61  60
62  4
62  85
63  64
64  65
, W& I* U2 n( i  w* R64  76
, f# q8 W1 b6 R8 y" e65  66
66  67
# O; `# i% _2 r5 r2 V66  76
67  44
5 S' x9 B: j: _9 Q67  68
68  69
9 _; A. E  b% m+ S8 D1 R68  75
69  70
69  71
69  1
70  2
70  43
71  72
9 `! r, ?0 D2 x- R* J  A71  74
7 ~" _) S' H+ D72  73
73  74
$ b5 ~8 s  k, N1 |9 e9 ?- {; j/ Y; Z73  18
5 E; {5 w- U$ b1 @# A2 C74  1
74  80
$ C% w- j* v- X- O75  76
- f$ s1 T. r- \0 W) Y1 y76  77
77  78
77  19
; D1 F. m! l+ r+ G78  79
79  80
80  18
81  82
82  83
82  90
( T7 m. n- N1 B! X/ H0 w% s83  84
. `  s+ R6 f3 _( K84  85
85  20
86  87
86  88
87  88
, t  T* y; F) a9 W2 K7 z87  92
: ]* f- v, v% _0 {88  89
88  91
89  20
89  84
89  90
90  91
91  92
7 P2 t) y! |3 y9 F% b];
x=data(:,1);
% `+ G. r+ l9 X  J4 R% M) vy=data(:,2);
" i. E5 O! c2 i- G- Lz=[6.264982043
6.403124237
9.486832981
42.46469122
15.23975065
45.60975773
10.30776406
5
8.485281374
16.03121954
11.40175425
% a% `7 j( P' [2 a) Q7 z12.80624847
  z# C+ z$ s& m+ ]7 d& ?11.5974135
20.79663434
; _( {4 N3 Y3 G9 a( L4.242640687
+ ~$ C. f, Z, P- C' f2 q9 b. D& s' i49.21635907
32.69556545
9 }/ M1 j# t* y, y# u; l$ ^) T9 a9
% Q+ O& x6 G5 k4 l9 L" n$ k  ~6 N+ ^17.88854382
384.4697647
" m  Y8 A4 y+ q& |& `32.64965543
38.18376618
40
67.41661516
8 l8 L4 W- Y( Z& w& c34.05877273
0 u9 j* E2 O# x3 r! h' }( f$ i26.87936011
9.848857802
40.22437072
) s. }- _1 g6 E6.708203932
5.385164807
* j2 |" v0 U' {0 m1 F4.472135955
( \  x0 H. N& N- Q. L: |# A3.605551275
18.02775638
( y" l" h! w$ Q1 o& x! f" h% }- q358.0460864
9.055385138
. [6 @, g* k/ n& \; b. b2 C5
347.6348659
23.85372088
18.02775638
20.02498439
7.433034374
$ r3 E2 K) G# B& }! k; _, _35.38361203
1 F; |. L3 N) |% P33.04920574
9.486832981
47.51841748
74.3236167
5.830951895
7.071067812
11.70469991
( e" H5 H" O2 Z$ o/ O15.53222457
0 L- f  x$ D) m7 }6 x- t5.099019514
6 U4 A5 V0 N" s$ ~4 Z% q7.566372975
8.276472679
9 S* i8 S: Y+ `: {3 }' Y5.024937811
$ i7 Z+ i5 p4 g' B1 Q5 l! q6.708203932
. _( {+ q* {" v! ~5
0 u! A: V% N. k( d  ^% r) L5.099019514
6.08276253
( L! M  I* ?7 }/ n9 D$ q7 ]1 x35.0142828
30.41381265
, N7 G. }( K  X# y8 h" @! t1 C3
+ N  [9 F2 d  `( m. q) S40.07804885
- X5 ^5 ~6 c( z. Y) H5 Q$ l0 B3 F17.67766953
: Q+ m6 D& o( x3 c" J- E0 o# D4 x; e19.14418972
5 v$ Z) m/ h* g' \( ?8.5
46.31684359
* S, e/ B; J5 Y- |& D8.062257748
8
& T" F9 _; y# I7 N8.062257748
11.62970335
$ A5 f/ ^7 r0 }% a! ?8 }9 m7 s6 r6
' H, p0 Q4 q* I* Z9.300537619
29.42787794
3 P1 s% y' g+ j- c- ?' E' i3 C10.19803903
14.56021978
5 W& a! i4 `$ R3 b6 W" R% O. s56.26944108
29
10.44030651
6.708203932
3.807886553
4.301162634
1 y1 T0 V) V+ G# @2.915475947
4 Y' T$ s0 w/ o; _; a: J/ b/ s. X4.242640687
8.544003745
4 s. _$ ~! s" A% z22.8035085
10.04987562
4 {# w9 C3 _9 [5 {24.18677324
12.6589889
12.3794184
7.5
3 N; Z2 `! ]4 U, b' U1 H9 g8.139410298
18.68154169
7.810249676
13.89244399
% m3 @$ `7 d" t& q9 W8 E- [34.71310992
+ M9 |+ |5 V* U: f! J- F3.5
$ `; s( _, a3 ?# x; \4 Y% Y60.01666435
! }# v2 D1 P; H9.055385138
5.830951895
13.15294644
1 m: J3 D0 g; G5 K3 k3 F0 L3.16227766
4.242640687
+ Y" d, w" s- ]- k9.219544457
: ?& Q2 T: Q% g+ p4 O7 }2 f14.76482306
/ S- Y( ], M$ ]4.123105626
7.071067812
4.527692569
* G5 y* W9 P2 H$ f8 l5.385164807
' h8 w* ^9 t+ T1 W" h6.403124237
5
8.602325267
3 Q$ i6 Z1 E0 R6 k, e3 Z/ d  t" a7.615773106
5
) y/ C$ }  ?5 [8 _' I6.103277808
% r; S) D5 q+ M- |: b- a5 L8.062257748
4.031128874
2 \2 n  K, F, O19.72308292
. o7 f+ P5 Z) g; ^* p0 n; N6.264982043
16.91892432
3.535533906
- P$ f) w+ O6 F4.472135955
10
9.848857802
& z5 `& I0 i* n  ^4 j7 O6.708203932
4.472135955
- [( u. E1 f$ `8.062257748
: r* N% ]& S* b+ R9 |7 N5.024937811
9 n& T3 t0 w- X  S3 g' Z5.408326913
8.732124598
9.848857802
7.280109889
4.472135955
# @( w& l2 I9 S, A3 F2 {11.04536102
9.340770846
4.031128874
( Z2 Z# C8 z# u- A9 q  _9 A' l21.37755833
4.031128874
, l. s! \- P. N) @3 w6 a  z3.041381265
9.486832981
3
' `9 g% K% u2 {4 o/ C, D) {3.535533906
4.74341649
20.02498439
];
xx=zeros(92);
& E0 [9 n- t" `for i=1:143
    xx(x(i),y(i))=z(i);
7 k5 V. K( ^$ n4 r3 `0 u    xx(y(i),x(i))=z(i);
end
* h5 T& J' l: P3 k, M* u3 W. Pfor i=1:92
" s+ ~$ k' x  b+ k    for j=1:92
    if xx(i,j)==0
& j. f0 N: M0 \: Q5 Z        xx(i,j)=10000;
: C. q! X3 S3 F# G# C$ O- X) Q2 _3 m* R: F& a    end
5 ]  O  [) p- [* s    end
4 i7 M- T- g0 _* k9 kend
for i=1:92
8 o& v# s% {2 ?/ S4 l/ G    xx(i,i)=0;
9 a* z' e! E6 G7 b/ Oend

" R$ K& M* X' Z1 [& N2 k4 P/ O0 x附件3:围堵方案的java实现程序
( K4 T- D% R9 q% Z* u1 F4 fpublic void SearchMax(ArrayList<SerachNode>  searchnode ,int[][]quantu,ArrayList<Node> allNode)
{
while(searchnode.size()>0)
4 p& L; K* a4 G% |9 U& N( w{
SerachNode tmp=searchnode.get(0);
searchnode.remove(0);

: p- }9 ?* y: `$ F: Zfor(int j=0;j<quantu.length;j++) //扩展当前结点
6 S. t( o! C' l5 J0 t{
int i=tmp.getBianhao();
: S3 j/ o7 ]* ]7 Sif(quantu[i][j]>0)
4 l) c, h* Y6 c# r8 Q7 B{
: P& X0 X( ^! S
# ~0 J/ ]$ @1 }4 DSerachNode newNode=new SerachNode(allNode.get                                                  (j),tmp.getCurrentQuanzhi()+quantu[i][j]) ;

$ t5 T- }! j! e& \9 _# ~9 k. O) X  Yif(!bounderMAX(newNode))
{
9 Q9 t6 r* e! U5 b  |& s  q$ XMaxresultSet.add(newNode);
' R5 Z* K5 p$ m7 R8 i3 Mcontinue;
1 K) d$ B$ k2 f3 E
}


if(!this.contain(searchnode, newNode))
/ B. O3 y* h. j) Z{
* j1 R7 s+ W9 l4 z+ x; ^. msearchnode.add(newNode);
! y0 t! n: |8 J9 I, G% n}
}}}
' }- {( f8 y5 {1 @
5 z" }3 ^4 M+ |- s3 c附件4:全区的交巡警平台有效覆盖范围(有效代表三分钟内可以到达)
9 k$ e' W4 H8 k! J/ x# J

/ k6 @# Z+ E+ f. m* K7 l附件5:用Mathmatica求数据均值与方差
data={9.1,5.6,7.6,9,8.2,4.6,8.5,4.8,5.3,3.4,8.3,7.2,4.5,5,1.6,4,2.5,6.4,7.9,10.6};
Mean[data]
Variance[data]

附件6:输入任意两点的坐标,输出两点间距离
. p6 ]# t1 U% ^! }# b$ p# G! P程序中以节点10和34间距离为例（Mathematica）
0 C7 i0 V1 o) x2 q& o' ~In[1]:
+ g8 j" ]& ~* D" Q: A0 lClear["Glibal'*"]
  b1 Z1 [% t* i; r6 s9 q: ox1=328;
x2=282;
y1=342.5;
y2=325;
d=Sqrt[(x1-x2)^2+(y1-y2)^2]
; ?; P" V4 t: l7 c3 U- C/ _7 {: NOut[1]:49.2164

8 D1 Q! B( Y0 q0 f: a% N# w附件7:A区各线路距离
路线起点
" f' N5 V2 B8 b/ c标号        路线终点
- B. R/ o* m6 V/ g) W. j标号        起点
/ \+ M! o, c+ K横坐标        起点
纵坐标        终点
横坐标        终点
纵坐标        距离
6 _! k& G8 t1 n) @1        75        413        359        418.5        356        6.2649820431
1        78        413        359        417        364        6.4031242374
2 ?/ X, q8 O$ Y2 j2        44        403        343        394        346        9.4868329805
3        45        383.5        351        342        342        42.464691215
3        65        383.5        351        395        361        15.239750654
$ O4 O! t' a' Y$ z% l9 u, O; w4        39        381        377.5        371        333        45.609757728
4        63        381        377.5        391        375        10.307764064
& ^; B1 ^/ x9 P5 c/ J) J* i3 T; o0 y5        49        339        376        342        372        5
5        50        339        376        345        382        8.4852813742
6        59        335        383        351        382        16.031219542
7        32        317        362        326        355        11.401754251
7        47        317        362        325        372        12.806248475
8        9        334.5        353.5        333        342        11.597413505
8        47        334.5        353.5        325        372        20.796634343
: m# W! G2 U0 |9 b9 M4 H$ |! w5 }9        35        333        342        336        339        4.2426406871
10        34        282        325        328        342.5        49.216359069
11        22        247        301        234        271        32.695565449
& N2 F* L4 R) o! j3 \" w11        26        247        301        256        301        9
/ V& {. w  w  ?% m$ X. o12        25        219        316        227        300        17.88854382
12        471        219        316        155        316        64
0 L! e2 I' ?' H( e14        21        280        292        251        277        32.649655435
8 i9 K9 V' D) r) f15        7        290        335        317        362        38.183766184
15        31        290        335        314        367        40
16        14        337        328        280        292        67.416615163
- D) H( G8 w- O16        38        337        328        371        330        34.058772732
. v5 A# K7 k6 n% i3 q# {9 z17        40        415        335        388.5        330.5        26.879360111
' O$ M# F( n4 I* W% R' q* r17        42        415        335        419        344        9.8488578018
17        81        415        335        438        368        40.224370722
' ]- V, L+ }' G7 J18        81        432        371        438        368        6.7082039325
18        83        432        371        434        376        5.3851648071
19        79        418        374        420        370        4.472135955
20        86        444        394        447        392        3.6055512755
9 E  U; I" ~2 z: S: [21        22        251        277        234        271        18.027756377
; o, X/ \  i. h: ]% j0 D6 K8 F22        372        234        271        232.5        264        7.1589105316
* U/ h  V0 g' G- S6 |9 @' V22        13        234        271        225        270        9.0553851381
23        13        225        265        225        270        5
9 N4 m# Z% ^5 R' y8 y23        383        225        265        192        264        33.015148038
24        13        212        290        225        270        23.853720884
4 v/ a0 ]$ V8 _# A/ j  ]5 d/ V4 z: `: O24        25        212        290        227        300        18.027756377
25        11        227        300        247        301        20.024984395
8 ~$ ?; B0 Y+ z/ O+ m# [# U26        27        256        301        250.5        306        7.4330343737
26        10        256        301        282        325        35.383612026
27        12        250.5        306        219        316        33.049205739
# f- h# X) D+ }+ I: `2 h% M28        29        243        328        246        337        9.4868329805
28        15        243        328        290        335        47.518417482
1 I% W# O! D( D$ [29        30        246        337        314        367        74.323616704
% t: g; W. P" i1 \, Y30        7        314        367        317        362        5.8309518948
2 j1 K3 f5 l" a! g30        48        314        367        315        374        7.0710678119
31        32        315        351        326        355        11.704699911
31        34        315        351        328        342.5        15.532224567
32        33        326        355        327        350        5.0990195136
7 b6 e+ B. y& E2 r33        34        327        350        328        342.5        7.5663729752
; _1 X- ^( B2 B33        8        327        350        334.5        353.5        8.2764726786
( c# y7 Y1 E  H- M34        9        328        342.5        333        342        5.0249378106
35        45        336        339        342        342        6.7082039325
36        35        336        334        336        339        5
36        37        336        334        331        335        5.0990195136
  ~9 a5 O$ Q) [8 z/ ~/ l36        16        336        334        337        328        6.0827625303
36        39        336        334        371        333        35.0142828
- o. C! d  Y8 h* v8 P; n- C37        7        331        335        317        362        30.413812651
/ L6 C) l2 y7 a0 ~8 l, r" a2 X38        39        371        330        371        333        3
38        41        371        330        411        327.5        40.078048855
39        40        371        333        388.5        330.5        17.67766953
' A* u7 g4 D2 v1 b& A2 T- J. G40        2        388.5        330.5        403        343        19.144189719
41        17        411        327.5        415        335        8.5
5 C8 b% }. \4 o6 k" \6 v. h& c41        92        411        327.5        444        360        46.316843588
! B' F5 y8 E7 \8 E6 \42        43        419        344        411        343        8.0622577483
43        2        411        343        403        343        8
43        72        411        343        418        347        8.0622577483
44        3        394        346        383.5        351        11.62970335
45        46        342        342        342        348        6
5 T, \9 D1 d3 `0 Y' u* X' d46        8        342        348        334.5        353.5        9.3005376189
( i, O1 d2 q/ G" X/ |2 Z46        55        342        348        371        353        29.427877939
8 ^; T; w7 a- t47        48        325        372        315        374        10.198039027
5 ^+ y# F3 {: F4 w' ]4 c47        6        325        372        339        376        14.560219779
47        5        325        372        381        377.5        56.269441085
5 o+ k3 M" Q2 g' U$ h2 p  y  g) {, M48        61        315        374        335        395        29
9 j# K: U4 k9 v' m$ M& m49        50        342        372        345        382        10.440306509
: x" `; L% j7 i49        53        342        372        348        369        6.7082039325
% E6 a; Y# \" J8 D  d: C% y50        51        345        382        348.5        380.5        3.8078865529
; x* _6 f1 r# N$ k6 F51        52        348.5        380.5        351        377        4.3011626335
51        59        348.5        380.5        351        382        2.9154759474
52        56        351        377        354        374        4.2426406871
- e% O3 q6 O" [: J! \& r# q# I53        52        348        369        351        377        8.5440037453
- w) c2 b2 |" E. W$ U" p53        54        348        369        370        363        22.803508502
54        55        370        363        371        353        10.049875621
54        63        370        363        391        375        24.186773245
55        3        371        353        383.5        351        12.658988901
5 H' A$ h4 D7 j; X" B. h56        57        354        374        363        382.5        12.379418403
57        58        363        382.5        357        387        7.5
57        60        363        382.5        369        388        8.139410298
57        4        363        382.5        381        377.5        18.681541692
9 B9 x# }( t1 y2 o0 `58        59        357        387        351        382        7.8102496759
: f, `1 i& v1 x% m; ]& B6 B4 o& Q60        62        369        388        381        381        13.892443989
61        60        335        395        369        388        34.713109915
62        4        381        381        381        377.5        3.5
! A: M3 J; M+ h! S6 k# o62        85        381        381        440        392        60.016664352
63        64        391        375        392        366        9.0553851381
4 u% M. Z* Y* h' C6 Y64        65        392        366        395        361        5.8309518948
9 P$ U: W8 b, r$ u' z3 j0 Q64        76        392        366        405        368        13.152946438
65        66        395        361        398        362        3.1622776602
+ M- v: k+ c* T- `66        67        398        362        401        359        4.2426406871
, j- I; E+ |+ W66        76        398        362        405        368        9.2195444573
67        44        401        359        394        346        14.76482306
67        68        401        359        405        360        4.1231056256
68        69        405        360        410        355        7.0710678119
; W4 N1 G# x1 E0 o/ B1 Y68        75        405        360        405.5        364.5        4.5276925691
69        70        410        355        408        350        5.3851648071
69        71        410        355        415        351        6.4031242374
69        1        410        355        413        359        5
70        2        408        350        403        343        8.602325267
2 A  E7 l- f1 g4 o6 x, K70        43        408        350        411        343        7.6157731059
% R- O1 G+ j  v4 w# g71        72        415        351        418        347        5
. T+ \2 c1 r" R& n9 M6 Y' Z4 z71        74        415        351        418.5        356        6.1032778079
9 r, ^! e' c( d' r) `72        73        418        347        422        354        8.0622577483
73        74        422        354        418.5        356        4.0311288741
73        18        422        354        432        371        19.723082923
; G# d* Z4 t  m9 G74        1        418.5        356        413        359        6.2649820431
4 d* o4 e  u  @& Z( z( u+ z74        80        418.5        356        424        372        16.918924316
; Y" u4 F  N  p8 ?75        76        405.5        364.5        405        368        3.5355339059
76        77        405        368        409        370        4.472135955
77        78        409        370        417        364        10
+ d. E% k! o/ b; c7 q77        19        409        370        418        374        9.8488578018
4 ~; M, m8 m4 O+ X8 ?78        79        417        364        420        370        6.7082039325
6 S& c, W6 {' P0 m4 l$ I5 K79        80        420        370        424        372        4.472135955
80        18        424        372        432        371        8.0622577483
( J. I# q) K' o( t1 L5 c81        82        438        368        438.5        373        5.0249378106
82        83        438.5        373        434        376        5.4083269132
) k( z4 l" o) m3 D82        90        438.5        373        440.5        381.5        8.7321245983
83        84        434        376        438        385        9.8488578018
9 N' k# r! h% L+ e84        85        438        385        440        392        7.2801098893
85        20        440        392        444        394        4.472135955
! r4 K  [1 K0 F  g3 U  ^86        87        447        392        448        381        11.045361017
9 j& o& ?6 @( {7 z86        88        447        392        444.5        383        9.3407708461
2 G1 y+ P; E9 C! v4 F' s) m5 x0 c87        88        448        381        444.5        383        4.0311288741
87        92        448        381        444        360        21.377558326
88        89        444.5        383        441        385        4.0311288741
88        91        444.5        383        445        380        3.0413812651
5 H6 o" W: G/ x/ Z" |9 {2 f! S& s+ `0 V89        20        441        385        444        394        9.4868329805
( k" U; l4 G0 [& e# b89        84        441        385        438        385        3
" K  O; O% N8 O4 V89        90        441        385        440.5        381.5        3.5355339059
& O$ ?" z% k: `* @# C90        91        440.5        381.5        445        380        4.7434164903
* [1 N! j, f# m9 Q" D0 k6 u4 |; W91        92        445        380        444        360        20.024984395
* l# R: L9 O. u- i) o
6 x: \1 U  [$ u3 M  P0 k( X

liupeng723911

楼主你太好了.........

hbdkfk2

谢谢楼主！！！！！

467857726de

不错哦  谢谢楼主

qqwhw2012

好东西
2 n7 ~; L# u! w* V* |111

往东就亮

谢谢楼主分享！！！！

霜之咏叹调

haobaniyignla  

李梦龙33

谢谢楼主分享，正在做。。。

莫西莫西

灰常好，就是看着太费劲

Dou1

太详细了，谢谢楼主