交巡警服务平台的设置和调度模型

1、PAGE 第 PAGE 2页，共 NUMPAGES 1页2013财经大学数学建模竞赛B题交巡警服务平台的设置与调度模型参赛队员: 参赛队编号：20130032013年8月7日8月11日 交巡警服务平台的设置与调度 摘要本文针对交巡警服务平台的设置与调度等相关问题进行建立模型和研究。首先，运用matlab软件对节点和相应距离组成的三列矩阵进行预处理得到邻接矩阵，基于floyd算法得到任意两点间的最短距离矩阵juli。然后针对问题（一）和问题（二）,依据lingo11.0对问题进行指派优化，得到合理的结论。 问题（1.1） 关于A区的交巡警平台得管辖分配，依据原则一“尽量在3分钟有交巡警到达事发地

2、”的前提下，以各交巡警平台任务量方差最小为目标函数，建立01规划模型，从而得出最佳分配方案。 问题（1.2） 依据木桶盛水原理，封堵A区13个路口方案的好坏取决于最晚到达指定封锁路口的交巡警到达时间的长短。建立以最晚到达时间最短为目标的优化模型，建立0-1规划模型。借助于lingo11.0软件编程，从而得出最佳交巡警平台调度方案。 问题（1.3） 对于个别节点20个中任何一个交巡警平台无论如何在三分钟之赶达事故发生地我们增设4个交巡警平台，分别在节点28、38、61、92 。通过这四个交巡警平台的增设使交巡警尽量在3分钟之赶到事故发生地且警务平台任务量趋于平均。问题（2.1）综合全市各区实际情

3、况，我们把案发率、各区人口、各区管辖面积作为一个相关因子，考虑各区节点到交巡警服务平台是否符合设置原则，一：“尽量在三分钟之赶到案发现场”；原则二“保持警务平衡”作为评判标准，确定各区是否设置合理。对于不合理区域，需要新增平台数建立的平台尽可能少，工作量尽可能的均匀，建立合理模型。问题（2.2）对P点案发三分钟后进行围堵。在模型建立应该遵循围最小原则，出动警员最小为原则利用对P进行围堵。利用lingo11.0软件编程和matlab软件编程从而得到最佳围堵方案。关键词：floyd算法 0-1规划 多目标优化，1问题重述警察肩负着重大职能，为了更好的贯彻职能，更好的为人民服务，需要在市区的重要部位

4、建立服务平台。如何在警察资源有限的情况下，根据城市实际情况和需求设置警务平台和有效的调度警务资源是警务部面临的实际问题。现在就针对某市设置交巡警服务平台的相关情况，建立数学模型分析问题如下：（一）（1.1）依据附件1和附件2中所给的该市区的网络图和数据，对各交巡警的服务平台进行合理分配，使其在所管辖的围出现突发事件时，尽量能在3分钟有交巡警（警车的时速为60km/h）到达事发地。（1.2）给出一个合理的交警服务平台的合理调度方案，在发生重大突发事件时，调度全区20个交巡警服务平台的警力资源，对进出该区的13条交通要道实现快速全封锁。实际中一个平台的警力最多封锁一个路口（1.3）根据现有交巡警服

5、务平台的工作量不均衡和有些地方出警时间过长的实际情况，拟在该区再增加2至5个平台，请确定需要增加平台的具体个数和位置。（二）（2.1）针对全市（主城六区A，B，C，D，E，F）的具体情况，按照设置交巡警服务平台的原则和任务，分析研究该市现有交巡警服务平台设置方案（参见附件）的合理性。如果有明显不合理，请给出解决方案。（2.2）如果该市地点P（第32个节点）处发生了重大刑事案件，在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警服务平台警力资源的最佳围堵方案。2 问题分析此题是研究交巡警的合理分配和调度的数学建模问题。需要我们通过建立合理的数学模型，进行多目标

6、优化，对交巡警进行合理的分配。首先根据附件提供的数据，通过matlab进行数据预处理，利用三列矩阵通过三列转邻接矩阵函数在matlab软件中得到相邻两节点的领接距离矩阵。基于floyd算法，借助于matlab软件进行编程，得到任意两点最短路径的距离矩阵juli.随后针对各个问题进行深入分析和处理得到合理解释。针对问题（一）问题（1.1）：依据交巡警服务平台原则之一：“尽量在三分钟之赶到事故的发生地”，我们通过floyd算法，计算出任意两点之间的最短距离。依据原则一，将各个节点到交巡警服务平台最近的原则，将交巡警管辖节点进行分类。时间对于任何一个交巡警服务平台都无法在3分钟之赶到的节点，我们采取

7、就近原则。问题（1.2）基于预处理得到的所有节点的距离矩阵距离，用matlab提取出20个服务平台分别到13个路口的最短距离 ,用此数据组合成一个2013新的矩阵,通过借助lingo11.0软件进行编程（value=1为表示该服务平台管辖这个节点出口，value=0，则说明该服务平台不接管该节点出口）得到最优调度结果在最短的时间能对A区13个出口进行封锁。问题（1.3）依据问题（1.1）所得出结果即各服务平台所管辖的路口，为了更合理的安排和调度交巡警服务平台，我们增设4个交巡警服务平台，使得28、29、38、39、61、92这几个节点能在事故发生三分钟之有交巡警的赶到，我们通过对（1.1）各节

8、点之间的最短距离调用，找出这六个点的之间任意两点之间的距离，在原则一的前提下将新增平台设在节点28、38、61、92.针对问题（二）问题（2.1）综合全市各区实际情况，我们把案发率、各区人口、各区管辖面积作为一个相关因子，考虑各区节点到交巡警服务平台是否符合设置原则，一：“尽量在三分钟之赶到案发现场”；原则二“保持警务平衡”作为评判标准，确定各区是否设置合理。对于不合理区域，需要新增平台数建立的平台尽可能少，工作量尽可能的均匀，建立合理模型。问题（2.2）对P点案发三分钟后进行围堵。在模型建立应该遵循围最小原则，出动警员最小为原则对P进行围堵。案发三分钟后，警员开始出动，可以计算出警员出动时，

9、嫌犯的最大活动围，然后根据对嫌犯可能所在的区域计算在他跑出路口之前将路口封锁最短时间。3模型假设与符号的说明3.1模型的假设为方便问题的求解与说明，我们对该模型做出如下前提假设：、所有的道路至少都可以双向通车，即车辆是可以相向形式。、车速不受任何外界因素影响，保持恒定。、事故发生和接到服务平台接受信息的时延忽略不计。、所有警车行动路线都选择最短路线。、事故发生，嫌疑犯选择逃跑路线是随机的。3.2符号说明三列矩阵：sanlie邻接矩阵:linjie 最短距离矩阵juli为 第i个节点到第j个节点的最短距离； X 0-1 整数规划矩阵；4模型的求解 首先对数据进行处理，每条路线的节点建立三列矩阵s

10、anlie，运用matlab编程编写函数threetolinjie得到邻接矩阵linjie,采用floyd算法得到任意两点之间的最短距离，作为路线的权重juli。其基本算法步骤为：、输入三列矩阵sanlie（2）、计算邻接矩阵（threetolinjie） linjie （3）、 计算任意两点之间的最短距离juli 数据详见附录1通过floyd算法计算出来的任意两点之间的最短距离为后面解题基础。 41.1要对交巡警服务平台进行合理分配建立数学模型 依据我们的数据预处理得到邻接矩阵，我们用matlab分别提取出A区20个服务平台到92个节点的距离矩阵mzyqJL，建立以交巡警额速度恒定的情况下，

11、服从原则一“尽量在三分钟之赶到事故发生地”下模型，通过函数zhaoshu依据最短距离分配原则，将各个节点分配到服务平台。图纸上1mm为现实中100米，警车速度为60km/h在图纸上为10mm/d。通过floyd算法计算出任意两点之间的距离矩阵juli中利用matlab提取92个点分别到20个警务平台的距离矩阵，利用zhaoshu函数将各个节点进行管辖划分，其中（ 1 i92 , 1 j20 ）建立目标函数：划分为第j个警务平台管辖。S.t 最后求的划分结果如下：平台序号管辖的节点案发率总和平台11、67、68、69、71、73、74、75、76、7810.3平台22、40、43、44、70、7

12、28.3平台33、54、55、65、665.6平台44、57、60、62、63、646.6平台55、49、50、51、52、53、56、58、599.7平台66、2.5平台77、30、32、47、489平台88、33、465平台99、31、34、35、458.2平台1010、1.6平台1111、26、274.6平台1212、254平台1313、21、22、23、248.5平台1414、2.5平台1515、2.1平台1616、36、373.8平台1717、41、425.3平台1818、80、81、82、836.1平台1919、77、793.4平台2020、84、85、86、87、88、89、90

13、、9110.7对辖区进行优化：对于在三分钟，警务人员无法到达的节点，我们按距离最近进行划分。最后得到划分结果如下：平台序号管辖的节点案发率总和平台11、67、68、69、71、73、74、75、76、7810.3平台22、40、43、44、70、72、399.7平台33、54、55、65、665.6平台44、57、60、62、63、646.6平台55、49、50、51、52、53、56、58、599.7平台66、2.5平台77、30、32、47、48、619.6平台88、33、465平台99、31、34、35、458.2平台1010、1.6平台1111、26、274.6平台1212、254平台

14、1313、21、22、23、248.5平台1414、2.5平台1515、28、294.8平台1616、36、37、385平台1717、41、425.3平台1818、80、81、82、836.1平台1919、77、793.4平台2020、84、85、86、87、88、89、90、91、9211.5对辖区再优化：为保持警务量均衡，在此采用lingo11.0 辖区节点按照在原则一下使得平均警务量方差最小为原则，保持警务量平衡进行优化:S.t 优化的结果：平台序号管辖节点案发率总和平台11、67、71、73、74、756.4平台22、3、39、42、727.9平台343、44、54、55、65、68、

15、767.4平台44、57、60、62、63、64、667.4平台55、49、53、58、596.7平台66、50,51、52、565.5平台79、30、32、48、617.7平台834、35、47、465.9平台97、16、33、376.5平台1010、1.6平台1111、26、274.6平台1212、254平台1313、21、22、23、248.5平台1414、2.5平台1515、28、29、316.4平台168、36、38、456.1平台1717、40、41、706.5平台1820、78、79、80、85、84、877.5平台1918、19、69、77、81、828.1平台2083、86、

16、88、89、90、91、927.24.1.2对出入A区的13条交通路口封锁的解决方案及模型木桶原理告诉我们：木桶的盛水容量大小取决于其中最短的那块木块的长度。所以对于某地发生大型案件，须封锁该地区时，封锁效果的好坏取决于最晚到达指定封锁地点的交巡警的到达时间的长短。因此该问题的目标函数为最晚到达封锁路口的交巡警的到达最短。采用0-1整数规划模型，0-1整数规划矩阵为(1i13,1j20).A区实现步骤：（1）、通过matlab提取每个出口到20个平台的距离矩阵juli13(i,j),将juli13（i，j）作为边的权重，运用lingo11.0得到最优的围堵方案。目标函数：结果如下：警务平台24

17、57101112131415161920出口节点16483029222423122128143862个点距离73.880773.958631.829380.154677.079238.052864.7759.7732.649747.518467.416676.392864.4888封锁路口的时间7.388077.395863.182938.015467.707923.805286.4775.9773.264974.751846.741667.639286.44888我们很清楚围堵这13个出口至少需要8.015分钟。4.1.3交巡警服务平台模型的优化改进对于问题1.1，我们发现节点28、29、3

18、8、39、61、92无法在三分钟之有交巡警赶到案发现场。基于原则一“尽量在三分钟之赶到案发现场”和任务均衡原理增设N个警务平台。实现步骤：通过matlab得出这六个点之间的任意两点之间最短距离矩阵juli6，S.t 节点 节点2889383961922809.4868145.014142.014119.8815231.4088299.48680153.6817150.6817110.3947240.076538145.014153.681703100.715386.394839142.014150.68173097.715389.394861119.8815110.3947100.715397

19、.71530145.999892231.4088240.076586.394889.3948145.99980所以在点28、39、61、92点增加服务平台可以使得所有节点在三分钟之可以在三分钟之赶到警务人员。管辖区分区如下：平台序号管辖节点平台11、67、71、73、74、75平台22、3、42、72平台343、44、54、55、65、68、76平台44、57、60、62、63、64、66平台55、49、53、58、59平台66、50,51、52、56平台79、30、32、48平台834、35、47、46平台97、16、33、37平台1010、平台1111、26、27平台1212、25平台13

20、13、21、22、23、24平台1414、平台1515、31平台168、36、45平台1717、40、41、70平台1820、78、79、80、85、84、87平台1918、19、69、77、81、82平台2083、86、88、89、90、91平台2828、29平台3938、39平台6161平台92924.2.1对全市交巡警服务平台的评价及优化：对全市六区的分配和调度合理不合理我们可以根据每个区的警务平台的平均警务量和到达每个节点是时间综合考虑，我们根据原则和任务，将各区的警务平台辖区进行分配：结果如下表显示：A区B区C区D区E区F区平均值不能赶到664712333523总共有的节点数9272

21、1545210310896.83333333平均警务量6.2258.311.0127.5337.969.9278.492833333平台数208179151113.33333333平均每个平台管辖个数，4.699.0585.7786.86679.8187.420116667各区面积22103221383432274239.1666667各区人口60214973765355.33333333本文将各个因素运用主成分分析，将主要因素归结为不能赶到、平均警务量、平均每个平台的管辖个数。在此建立多目标判定函数y=1*x1+2*x2+3*x3(x1不能在三分钟赶到的节点个数，x2个平台平均警务量，x3平

22、均每个平台管辖个数)，基于模糊原理，不妨将各目标变量的权重定位1=2=3=1，得出判定结果如下表：ABCDEF平均值目标函数值15.2795457620.9269678759.8907654222.9639721443.3274254849.4874469635.10179292 由上表可发现，C、E、F区判定函数明显高于全市判定函数值，该区存在明显不合理，依据警务平台原则和任务“警务平台到仅量3分钟到达事故发生地，和保持警务平衡”。本文通过对各区增加警务平台的解决方案，使各区警务平台设置趋于合理。我们选择在C、E、F区增加警务平台，增加节点如下：C平台序号辖区节点平台169333、169、1

23、70平台173337、173、174、175、176、177平台179344、180、182平台185349、165、185平台186346、186、165平台187352、193、194平台189353、189、190、191平台198360、196、197、198平台208364、200、201、202、203、204、206、207、208、209、211平台212205、212平台166166平台167167平台171171平台175175平台178178平台183183平台184184平台188188平台192192E平台序号管辖节点平台385374、382、383、394、385、

24、386、387、388、389、390、391、394、395平台378375、376、377、378、379平台391380、381、392、393、391平台372372、373平台396396平台397397平台398398平台399399平台400400平台401401、402平台403403平台404404平台405405平台406406F平台序号管辖节点平台486484、485、486、487、488、489平台506506、平台507507、平台498494、497、498、499、500、501平台502502、503平台505505、516平台482480、482平台4754

25、75平台476476平台478478平台479479平台4834834.22 在全市围的围堵方案对于P点嫌犯逃跑的围确定。对问题的假设：、嫌犯的逃离速度恒定与警车保持一致，这样就不存在警车从后面追击问 题。、嫌犯逃跑路线随机选择。、嫌犯逃跑只能走互通的道路。P点逃行路线说明：该红线所围区域为嫌疑犯在P点案发三分钟后所能到达的最大围，由于嫌疑犯只能在互通道路上逃行，故只能在如图“黄点”位置。找P点的步骤：（1）、由于案发地32号节点距离出入A区的30号、48号、16号路口的距离分别为17、24、33，所以该嫌犯在案发3分钟后有逃出A区的可能，为了达到快速搜捕嫌犯的目的，围堵的围越小越好。找出3分

26、钟后P逃犯可到达的点，按不精确计算可到达节点：23729151016393555661235有A，C两区的可能。（2）、把上面到达的点，看看所属区域，分别对所属区进行出口封锁。由于警务平台封锁时间差P点根据节点顺次走下去，下次可能到达的点及所属区： 2362382472823937028261438560 4384065445449505960234 173所属区有A,C,F。根据上提计算，所有A区的警务人员最快需要8.01546分钟围堵A好A区的出口，堵住C需要9.021359分钟，堵住F需要7,041784分钟，此时，最长时间是需要9.021359分钟围住A,C,F区，总长时间为12.02

27、1359。而此时逃犯能跑到的节点为：1234567891011141516171926272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980170171172173189190215216217218219224226227228229230231232233234235236237238239240241242243244245246247253482488532548549550 558560561562这些节点仍属于A,C,F区

28、域，逃犯只能跑到如图的区域，跑不出去。最大围也在A,C,F。所以此围堵方案合理。 图P2- P点最后被围堵的最大活动围各区警务平台围堵出口详细情况：A区：出口站点 16 48 30 29 22 24 23 12 21 28 14 38 62警务平台2 4 5 7 10 11 12 13 14 15 16 19 20C区：出口站点 177 202 203 264 317警务平台 175 177 179 166 181F区：出口站点483 541 572 578警务平台483 484 485 479 5 模型的评价与推广、综合运用floyd算法，matlab软件，excel，0-1规划，及ling

29、o11.0对数据进行深入分析依据交巡警平台设置原则和任务对交警平台做出合理设置方案，保证最高效使用警力，最优效服务群众。、整个模型思路清楚，假设合理。在较理想的情况下，对交巡警的调度和安排具有较强的可操作性。该模型的缺点是在考虑问题的时候主要考虑的是尽量在三分钟赶到事故发生地，和保持警务平衡，没有充分的考虑人口和管辖面积因素。、在模型的推广中。该模型具有较强的参照性。在遇到重大突发事件中对路口进行封锁，对警务平台的封锁路线具有一定的参照性。在封锁模型中，该模型的局限在于不能对实际路况进行合理的考虑。、随着经济社会的发展，新中国城镇化的过程中，在交巡警平台最优选址中，该模型也可运用该模型进行最优

30、化选择。6 参考文献1：钱湔.运筹学M.：科学，2000 2：2:薛定宇，.初等运用数学效果的 matlab 求解M.:清华大学出版 社，2004.,8 3：石辛民，郝正清.基于 matlab 的适用数值计算M.：清华大学出版， 交通大学，2006,2 5：4:金星，优化建模与lingo 7：附录附录：函数求三分钟能赶到的节点function mzyqJL=zhaoshu(a) % a为juli for i=1:154 k=0; min=a(i,1); m=1; for j=1:154 if mina(i,j) min=a(i,j); m=j; end if a(i,j)=30 k=k+1;

31、m1=a(i,j); else m1=inf; m2=i; m3=j; h=m1,m3,m2; end mzyqJL(i,j)=m1; end %mzyqJL(i,36)=k; mzyqJL(i,155)=min; mzyqJL(i,156)=m;end函数floyd算法：function d,r=floyd(a) n=size(a,1); d=a; for i=1:n for j=1:n r(i,j)=j; end end r; for k=1:n for i=1:n for j=1:n if d(i,k)+d(k,j)0);for(links(i,j):x(i,j)=1);for(link

32、s(i,j):bin(x(i,j);Lingo 封锁路口 model:sets:comand/A1 .A20/;cross/1.13/;roads(comand,cross):w;tsx(comand,cross):x;endsetsdata:w=file(C:UsersAdministratorDesktopjulichukouC.txt);enddatafor(comand(i):sum(cross(j):x(i,j)=1);for(cross(j):sum(comand(i):x(i,j)=1);for(roads(i,j):w(i,j)*x(i,j)=cyc);min=cyc;for(

33、tsx:bin(x);end 得到结果的部分数据： X( A1, 1) 0.000000 0.000000 X( A1, 2) 0.000000 0.000000 X( A1, 3) 0.000000 0.000000 X( A1, 4) 0.000000 0.000000 X( A1, 5) 0.000000 0.000000 X( A1, 6) 0.000000 0.000000 X( A1, 7) 0.000000 0.000000 X( A1, 8) 0.000000 0.000000 X( A1, 9) 0.000000 0.000000 X( A1, 10) 0.000000 0.

34、000000 X( A1, 11) 0.000000 0.000000 X( A1, 12) 0.000000 0.000000 X( A1, 13) 0.000000 0.000000 X( A2, 1) 0.000000 0.000000 X( A2, 2) 0.000000 0.000000 X( A2, 3) 1.000000 0.000000 X( A2, 4) 0.000000 0.000000 X( A2, 5) 0.000000 0.000000 X( A2, 6) 0.000000 0.000000 X( A2, 7) 0.000000 0.000000 X( A2, 8)

35、0.000000 0.000000 X( A2, 9) 0.000000 0.000000 X( A2, 10) 0.000000 0.000000 X( A2, 11) 0.000000 0.000000 X( A2, 12) 0.000000 0.000000 X( A2, 13) 0.000000 0.000000 X( A3, 1) 0.000000 0.000000 X( A3, 2) 0.000000 0.000000 X( A3, 3) 0.000000 0.000000 X( A3, 4) 0.000000 0.000000 X( A3, 5) 0.000000 0.00000

36、0 X( A3, 6) 0.000000 0.000000 X( A3, 7) 0.000000 0.000000 X( A3, 8) 0.000000 0.000000 X( A3, 9) 0.000000 0.000000 X( A3, 10) 0.000000 0.000000 X( A3, 11) 0.000000 0.000000 X( A3, 12) 0.000000 0.000000 X( A3, 13) 0.000000 0.000000 X( A4, 1) 0.000000 0.000000 X( A4, 2) 0.000000 0.000000 X( A4, 3) 0.00

37、0000 0.000000 X( A4, 4) 0.000000 0.000000 X( A4, 5) 0.000000 0.000000 X( A4, 6) 0.000000 0.000000 X( A4, 7) 0.000000 0.000000 X( A4, 8) 0.000000 0.000000 X( A4, 9) 0.000000 0.000000 X( A4, 10) 0.000000 0.000000 X( A4, 11) 0.000000 0.000000 X( A4, 12) 1.000000 0.000000 X( A4, 13) 0.000000 0.000000 X(

38、 A5, 1) 0.000000 0.000000 X( A5, 2) 0.000000 0.000000 X( A5, 3) 0.000000 0.000000 X( A5, 4) 0.000000 0.000000 X( A5, 5) 0.000000 0.000000 X( A5, 6) 0.000000 0.000000 X( A5, 7) 0.000000 0.000000 X( A5, 8) 0.000000 0.000000 X( A5, 9) 0.000000 0.000000 X( A5, 10) 1.000000 0.000000 X( A5, 11) 0.000000 0

39、.000000 X( A5, 12) 0.000000 0.000000 X( A5, 13) 0.000000 0.000000 X( A6, 1) 0.000000 0.000000 X( A6, 2) 0.000000 0.000000 X( A6, 3) 0.000000 0.000000 X( A6, 4) 0.000000 0.000000 X( A6, 5) 0.000000 0.000000 X( A6, 6) 0.000000 0.000000 X( A6, 7) 0.000000 0.000000 X( A6, 8) 0.000000 0.000000 X( A6, 9)

40、0.000000 0.000000 X( A6, 10) 0.000000 0.000000 X( A6, 11) 0.000000 0.000000 X( A6, 12) 0.000000 0.000000 X( A6, 13) 0.000000 0.000000 X( A7, 1) 0.000000 0.000000 X( A7, 2) 0.000000 0.000000 X( A7, 3) 0.000000 0.000000 X( A7, 4) 0.000000 0.000000 X( A7, 5) 0.000000 0.000000 X( A7, 6) 0.000000 0.00000

41、0 X( A7, 7) 0.000000 0.000000 X( A7, 8) 0.000000 0.000000 X( A7, 9) 1.000000 80.15460 X( A7, 10) 0.000000 0.000000 X( A7, 11) 0.000000 0.000000 X( A7, 12) 0.000000 0.000000 X( A7, 13) 0.000000 0.000000 X( A8, 1) 0.000000 0.000000 X( A8, 2) 0.000000 0.000000 X( A8, 3) 0.000000 0.000000 X( A8, 4) 0.00

42、0000 0.000000 X( A8, 5) 0.000000 0.000000 X( A8, 6) 0.000000 0.000000 X( A8, 7) 0.000000 0.000000 X( A8, 8) 0.000000 0.000000 X( A8, 9) 0.000000 0.000000 X( A8, 10) 0.000000 0.000000 X( A8, 11) 0.000000 0.000000 X( A8, 12) 0.000000 0.000000 X( A8, 13) 0.000000 0.000000 X( A9, 1) 0.000000 0.000000 X(

43、 A9, 2) 0.000000 0.000000 X( A9, 3) 0.000000 0.000000 X( A9, 4) 0.000000 0.000000 X( A9, 5) 0.000000 0.000000 X( A9, 6) 0.000000 0.000000 X( A9, 7) 0.000000 0.000000 X( A9, 8) 0.000000 0.000000 X( A9, 9) 0.000000 0.000000 X( A9, 10) 0.000000 0.000000 X( A9, 11) 0.000000 0.000000 X( A9, 12) 0.000000

44、0.000000 X( A9, 13) 0.000000 0.000000 X( A10, 1) 0.000000 0.000000 X( A10, 2) 0.000000 0.000000 X( A10, 3) 0.000000 0.000000 X( A10, 4) 0.000000 0.000000 X( A10, 5) 1.000000 0.000000 X( A10, 6) 0.000000 0.000000 X( A10, 7) 0.000000 0.000000 X( A10, 8) 0.000000 0.000000 X( A10, 9) 0.000000 0.000000 X

45、( A10, 10) 0.000000 0.000000 X( A10, 11) 0.000000 0.000000 X( A10, 12) 0.000000 0.000000 X( A10, 13) 0.000000 0.000000 X( A11, 1) 0.000000 0.000000 X( A11, 2) 0.000000 0.000000 X( A11, 3) 0.000000 0.000000 X( A11, 4) 0.000000 0.000000 X( A11, 5) 0.000000 0.000000 X( A11, 6) 0.000000 0.000000 X( A11,

46、 7) 1.000000 0.000000 X( A11, 8) 0.000000 0.000000 X( A11, 9) 0.000000 0.000000 X( A11, 10) 0.000000 0.000000 X( A11, 11) 0.000000 0.000000 X( A11, 12) 0.000000 0.000000 X( A11, 13) 0.000000 0.000000 X( A12, 1) 0.000000 0.000000 X( A12, 2) 0.000000 0.000000 X( A12, 3) 0.000000 0.000000 X( A12, 4) 0.

47、000000 0.000000 X( A12, 5) 0.000000 0.000000 X( A12, 6) 1.000000 0.000000 X( A12, 7) 0.000000 0.000000 X( A12, 8) 0.000000 0.000000 X( A12, 9) 0.000000 0.000000 X( A12, 10) 0.000000 0.000000 X( A12, 11) 0.000000 0.000000 X( A12, 12) 0.000000 0.000000 X( A12, 13) 0.000000 0.000000 X( A13, 1) 1.000000

48、 0.000000 X( A13, 2) 0.000000 0.000000 X( A13, 3) 0.000000 0.000000 X( A13, 4) 0.000000 0.000000 X( A13, 5) 0.000000 0.000000 X( A13, 6) 0.000000 0.000000 X( A13, 7) 0.000000 0.000000 X( A13, 8) 0.000000 0.000000 X( A13, 9) 0.000000 0.000000 X( A13, 10) 0.000000 0.000000 X( A13, 11) 0.000000 0.00000

49、0 X( A13, 12) 0.000000 0.000000 X( A13, 13) 0.000000 0.000000 X( A14, 1) 0.000000 0.000000 X( A14, 2) 0.000000 0.000000 X( A14, 3) 0.000000 0.000000 X( A14, 4) 1.000000 0.000000 X( A14, 5) 0.000000 0.000000 X( A14, 6) 0.000000 0.000000 X( A14, 7) 0.000000 0.000000 X( A14, 8) 0.000000 0.000000 X( A14

50、, 9) 0.000000 0.000000 X( A14, 10) 0.000000 0.000000 X( A14, 11) 0.000000 0.000000 X( A14, 12) 0.000000 0.000000 X( A14, 13) 0.000000 0.000000 X( A15, 1) 0.000000 0.000000 X( A15, 2) 0.000000 0.000000 X( A15, 3) 0.000000 0.000000 X( A15, 4) 0.000000 0.000000 X( A15, 5) 0.000000 0.000000 X( A15, 6) 0

51、.000000 0.000000 X( A15, 7) 0.000000 0.000000 X( A15, 8) 1.000000 0.000000 X( A15, 9) 0.000000 0.000000 X( A15, 10) 0.000000 0.000000 X( A15, 11) 0.000000 0.000000 X( A15, 12) 0.000000 0.000000 X( A15, 13) 0.000000 0.000000 X( A16, 1) 0.000000 0.000000 X( A16, 2) 1.000000 0.000000 X( A16, 3) 0.000000 0.000000 X( A16, 4) 0.000000 0.000000 X( A16, 5) 0.000000 0.000000 X( A16, 6) 0.000000 0.000000 X( A16, 7) 0.000000 0.0