2011年全国大学生数学建模B题论文材料.doc

-2011高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题.我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出.我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性.如有违反竞赛规则的行为,我们将受到严肃处理.我们参赛选择的题号是(从A/B/C/D中选择一项填写): B 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名): 参赛队员 (打印并签名) :1. 2. 3. 指导教师或指导教师组负责人 (打印并签名): 日期: 2011 年 9 月 11 日赛区评阅编号(由赛区组委会评阅前进行编号):2011高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评阅人评分备注全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):交巡警服务平台的设置与调度摘 要由于警务资源是有限的,所以根据城市的实际情况与需求,合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是有关部门面临的一个实际课题.本文着力于通过所给资料,寻找最优化的交巡台设置与调度方案. 按照设置交巡警服务平台的原则和任务,我们首先对问题1用Floyd算法,提出最佳的交巡警服务平台管辖区域划分方案,缩短了出警时间,平衡了工作量,然后采用回溯法,给出了应对突发事件的警力比较合理调度方案；对于问题2,我们将其归结为全局的配置问题,首先用优化后的Floyd算法对该市现有六城区的交巡警服务平台设置进行改进,其次以时间最短、围堵区域最小为原则,提出了应对重大刑事案件的最佳围堵方案.对于问题1,本文将最短时间问题转化为单向最短路径问题.我们没有运用经典的求最短距的Dijkstra算法,采取时间复杂度更简便的Floyd算法,应用Matlab编程,以出警时间最短为原则,将72个交通节点分配给20个交巡警服务平台；对于出现突发事件,本文采用回溯法,以最节省警力、实现全区封锁联动时间（即封锁路口最长时间）最短为目标,成功的实现了应对突发事件时警力的合理调度；对于某些交巡警服务平台工作量大、出警时间过长等问题,本文利用Mathematica对附表2中的数据进行分析,整理分析A区各节点事故发生率后,利用图论的相关知识,提出应增设4个服务平台,基本实现警力的最优配置.最后,借助于Matlab和Mathematica软件,对附件中所提供的数据进行了筛选,去除异常数据,对残缺数据进行适当补充,并从中随机抽取了3组数据(每组8个采样)对理论结果进行了数据模拟,结果显示,理论结果与数据模拟结果吻合良好.而对于问题2,我们对附件中所提供的A,B,C,D,E,F六城区的数据进行了整合与分析,并做出了直观的图表.遵循警情主导警务原则、快速出警原则、方便与安全原则,并结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况,在充分考虑现有警力和财力并确保安全的条件下,科学分析现有平台的数量和具体位置的合理性.数据显示C区和F区的事故发生率较高、交巡警服务平台工作量高于全市平均水平、交巡警服务平台平均每天出警时间过长,针对以上问题我们再次利用均衡二分法,并考虑区域边界处的设点拥挤问题,提出了在C区增设5个交巡平台、F区增设1个交巡平台.对于该市地点P(第32个节点)处发生了重大刑事案件的围堵问题,本文将其归结为资源调配问题.本文合理假设了犯罪嫌疑人的车行驶速度(分三种情况考虑:等于警车速度,警车速度的二倍,警车速度的一半),确定三分钟后犯罪嫌疑人逃逸的可能覆盖范围,从而利用回溯法的思想采用Matlab编程确定犯罪嫌疑人的车的所有可能位置.以时间最短、围堵区域最小为原则,采用改进的穷举算法,快速地形成围堵区域,并实现了围堵区域最小的目的.实现了资源调配问题的优化决策.考虑到该城市未来发展规划,只需对本文所建模型进行适当改进即可,在此不进行详细解答.关键词 最短路径 Floyd算法 回溯法 穷举法 优化决策目 录交巡警服务平台的设置与调度1摘 要11.问题重述12.问题分析12.1对于问题一的分析12.2对问题二的分析13.模型假设24.定义与符号说明25.模型的建立与求解25.1 问题一的模型25.1.1 模型建立25.1.2 模型求解35.2 问题二的模型85.2.1 模型建立85.2.2 模型求解97.模型的评价与推广108. 附件10附件1:用Floyd算法分配个服务平台管辖区域10附件2:邻接矩阵的matlab实现程序22附件3:围堵方案的java实现程序29附件4:全区的交巡警平台有效覆盖范围(有效代表三分钟内可以到达)30附件5:用Mathmatica求数据均值与方差30附件6:输入任意两点的坐标,输出两点间距离30附件7:A区各线路距离311.问题重述“有困难找警察”,是家喻户晓的一句流行语.警察肩负着刑事执法、治安管理、交通管理、服务群众四大职责.为了更有效地贯彻实施这些职能,需要在市区的一些交通要道、人员密集区和重要部位设置交巡警服务平台.每个交巡警服务平台的职能和警力配备基本相同.由于警务资源的有限性,根据城市的实际情况与需求,合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题.本文着力于寻找最优化的设置与调度方案.问题1要求合理分配交巡警服务平台的管辖范围,使其在所管辖的范围内出现突发事件时,尽量能在3分钟内有交巡警(警车的时速为60km/h)到达事发地；对于重大突发事件,给出该区交巡警服务平台警力合理的调度方案,尽快封锁道路；拟在该区内再增加2至5个平台,以减少出警时间、平均工作量,确定需要增加平台的具体个数和位置.问题2要求分析研究该市现有交巡警服务平台设置方案的合理性并给出解决方案；如果该市地点P(第32个节点)处发生了重大刑事案件,在案发3分钟后接到报警,犯罪嫌疑人已驾车逃跑.为了快速搜捕嫌疑犯,给出调度全市交巡警服务平台警力资源的最佳围堵方案.2.问题分析本题所要解决的是A区以及全市的安巡警服务平台设置与调度问题,根据现实生活状况,我们首先要考虑的是警力资源的限制,即要使得所布置的警力尽可能的少.其次是在交巡台数量最少的情况下,力求警员到达现场的时间在3分钟以内,解决突发状况.2.1对于问题一的分析该市中心城区A的交通网络有92个节点和20个交巡警服务平台,要求当突发事件发生时,尽量能在3分钟内有交巡警到达事发地,已知警车的时速为V=60km/h,我们将最短时间转化为最短路问题,应用Floyd算法,求解出A区距离每一节点最近的交巡台,即将该节点分配给该交巡台.对于重大突发事件,要实现对进出该区的13条交通要道进行快速封锁,即需调度交巡台尽快到达13个节点,重复Floyd算法,找出最近交巡台,即可找出调配方案.但需注意的是,有的出入口本来就有交巡台,但为了达最优化,需进行重新分配,故应用回溯法,找到调度方案.现有交巡台工作量不均衡和有些地方出警时间过长,统计A区各个交巡台案发率,计算均值与方差,在案发率较高地带增设交巡台,平衡工作量,尽量缩短出警时间.2.2对问题二的分析对于问题二,是对问题一的进一步改进与推广,在遵循警情主导警务原则,快速出警原则与方便与安全原则,结合辖区地域特征、人口分布和治安状况等实际情况,充分考虑现有警力和财力并确保安全的条件下,设置交巡平台,重复上一问的做法,评估交巡平台的合理性.对于改进方案,应考虑城区内部工作量,城区之间的联系以及城市边界的警力调度.对于突发状况的围堵方案,应在最短时间内对可能逃逸区域进行合围,最小范围内缩小包围圈.3.模型假设1.假设题中所给数据均真实可靠.2.出警时道路恒畅通(无交通事故、交通堵塞等发生),警车行驶正常,警车及肇事车辆行驶时均以60km/h匀速行驶,转弯处不需要花费时间.3.事故均发生在路口节点,两节点连线上认为没有事故发生.4.每条线路行驶都是双向的.5.考虑肇事车辆在P点向各个方向逃逸的概率相等.6.在整个行驶中,车辆只在主要干道行驶.7.发生事故时,忽略反应调度时间.4.定义与符号说明 任意两个标志点与之间的距离标志点间的距离组成的距离矩阵标志点的邻接矩阵邻接矩阵的元素相邻标志点间的距离矩阵相邻标志点与间的距离标志点的权值矩阵 标志点间的最短距离矩阵标志点与之间的最短距离肇事车辆逃逸速度 5.模型的建立与求解5.1 问题一的模型5.1.1 模型建立此问是关于最短路径的模型分析及MATLAB的实现A区道路状况及交巡台的设置如图1所示.本文应用Floyd算法,通过构造距离矩阵,依次找出距离每一节点最近的交巡台,使得有事故发生时,交巡警在最短时间内到达事故现场,以此为依据分配管辖区域.如果道路不通时,认为两端节点的距离为无穷. 图1 A区各节点及服务平台示意图当有重大突发事件时,要对进出该区的13条交通要道进行快速封锁,固定13个出入口,应用回溯法,找到距离节点最近的交巡平台.封锁时间决定于最后到达节点的时间,由于一个平台的警力最多封锁一个路口,至少需调动13个平台的警力.为达到工作量的均衡和出警时间尽可能的短,需进行优化决策.考虑每一节点案发率的不同,在A区增设2到5个平台,使得每一平台的工作量均衡,平均出警时间大体相同.5.1.2 模型求解首先我们可以根据题中所给的各个标志点的坐标,用matlab计算出任意两点之间的直线距离,得到92*92的距离矩阵:根据题中的分布图,我们可以得到各标志点的邻接矩即如果两个点相邻,则邻接矩阵中相对应的元素的值为1,否则为0；例如:3和44这两个点相邻,那么. 根据Floyd算法,我们是要求出任意两节点之间的距离,所以我们需要得到相邻两个结点的直线距离.我们可以利用距离矩阵的元素与的点乘积得到相邻标志点间的距离矩阵:对于D中不相邻点间距离0改为无穷大(Inf)从而得到节点与节点间的权值矩阵:即如果15和10之间不相邻,也即不能直接到达,那么D中的和都将变成和等于无穷大(Inf),否则则等于D中相应元素的数据.运用Floyd算法求出任意两点间最短距离,得到最短距离矩阵:由Floyd算法,运行MATLAB程序,可统计出距离每一节点最近的交巡台的位置,MATLAB运行结果如表1所示.带括号的节点为发生事故时任意交巡台都不能在三分钟内赶到节点.交巡台节点距离交巡台节点距离132127.0831 45718.681513229.055465823.841413235.0000 65916.0312132423.8537 46017.9240122517.8885 4(61)52.105511269.0000 4623.5000112716.4330 46310.308715(28)47.5184 4649.363215(29)57.005236515.23987305.8310 36618.401293120.557216714.915873211.40181756810.79278338.27651695.00009345.0249 2708.60239354.2426 1747111.265016366.0828 27216.4031163711.1818 187319.723116(38)34.0588 1746.26502(39)36.82191756.265024019.1442 1769.800517418.5000 19779.848917429.84891786.40312438.000019794.47212449.846818808.062394510.9508 18816.70828469.3005188210.793574712.806218835.385274812.9021208411.75225495.000020854.47215508.485320863.605055112.8932208714.651155217.1944208812.946455311.7082208914.752235422.7089189019.525635512.6590209116.006055621.437020(92)36.0060表1 该市A区指定节点到交巡警服务平台最短距离由上表可初步确定A区20个交巡台的管辖范围,如表2所示.带括号的节点为发生事故时任意交巡台都不能在三分钟内赶到节点.交巡台序号辖区内节点辖区内案发率交巡台序号辖区内节点辖区内案发率167 68 69 71 74 75 76 78 9.42 40 43 44 70 72 399.7354 55 65 665.6457 60 62 63 646.6549 50 51 52 53 567.7658 594.5730 32 47 48 619833 465931 34 35 458.210 1.61126 274.612 2541321 22 23 248.514 2.515(28) (29)4.81636 37 (38) 51741 425.31873 80 81 82 8371977 793.42084 85 86 87 88 89 91 90 (92)11.5表2 该市A区交巡警服务平台所管辖交叉路口清单图2 A区各交巡台管辖区域示意图 需要说明的是,同一条路整体归一个交巡台管理.当有重大突发事件时,固定13个进出A区的节点,运用回溯法,结合上表,找到距离节点最近的交巡台,以此来达到总体时间的最短,我们一共可以得到四个方案,在这个过程中可以发现,有些交巡台要避免去最近的节点封锁而去较远的节点,以此来节省警力.具体封锁方案如表3、表4所示.最短调度时间均为8.0155.方案一:交巡台过程出入口节点号2路径403938时间3.9822min4路径62时间0.3500min6路径474830时间3.1829min7路径3029时间8.0154min8路径4748时间3.0995min9路径353616时间1.5083min10路径262712时间7.5863min11路径22时间3.2696min12路径2524时间3.5916min13路径23时间0.5000min14路径21时间3.2649min15路径28时间4.7518min16路径14时间6.7417min表3 A区突发事件封锁方案一方案二方案三方案四路口标号平台号路口标号平台号路口标号平台号1213121312101416142314161661691662114211121142210221022122311231423132412241224112815281528152972972973083063093819381738148548948862206220622表4 A区突发事件封锁方案二、三、四在对交巡台均衡工作量,加快出警时间方面,综合各节点的案发率、交巡台到其辖区内任一节点的路程进行综合评估,做出优化决策.在案发率较高地带增设交巡台,以缓解周围交巡台的工作压力,为达均衡工作量的目的,将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所示.交巡台序号负责区域内的节点管辖区域内的案发率平均每天出警时间169 71 74 75 786.63.5069240 43 70 72 7.26.6736354 55 445.24.9314457 60 61 62 636.47.7015549 50 51 52 53 567.75.9456658 59 476.16.3949730 485.93.5258833 46 5.0 2.2748935 454.92.1271101.61126 274.62.394612254.02.86221323 245.73.8239142.515(28) (29)4.814.15801636 37 (38) (39)6.410.05621741 425.32.56891873 80 81 835.93.84381977 793.41.14572085 86 87 (92)6.45.535421222.82.52393132 344.94.39626664 65 67 68 76 5.12.66559082 84 88 89 916.23.2171表5 优化后的A区交巡台管辖区域示意图图3 A区增设平台示意图 图中方块所示节点即为增设平台处.5.2 问题二的模型5.2.1 模型建立对于问题2,对附件中所提供的A,B,C,D,E,F六城区的数据进行整合,做出直观的图表.遵循警情主导警务原则、快速出警原则、方便与安全原则,结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况,充分考虑现有警力和财力并确保安全,科学分析现有平台的数量和具体位置的合理性.对于该市地点P(第32个节点)处发生了重大刑事案件的围堵问题,本文将其归结为资源调配问题.本文合理假设了犯罪嫌疑人的车行驶速度(分三种情况考虑:等于警车速度,警车速度的二倍,警车速度的一半),并确定三分钟后犯罪嫌疑人的车行驶的最远距离,从而利用回溯法的思想采用Matlab编程确定犯罪嫌疑人的车的所有可能位置.以时间最短、围堵区域最小为原则,采用改进的双层Floyd算法,快速地形成围堵区域,并使围堵区域尽可能的小.5.2.2 模型求解全市整体状况如表5所示,数据显示C区和F区的事故发生率较高、交巡警服务平台工作量高于全市平均水平且交巡警服务平台平均每天出警时间过长,针对以上问题本文再次利用问题1的Floyd算法,并考虑区域边界处的设点拥挤问题,本文提出了在C区增加5个服务平台、在F区增加1个服务平台.全市六个城区城区面积城区人口平台数平均人口全区案发率各区平台案发率均值A2260202.727124.56.625B1032180.20466.48.3C22149170.223187.211.012D3837390.19167.87.533E43276150.176119.47.96F27453110.193109.29.927均值53.333313.33330.619112.41678.5595表6 全市整体状况图4 全市增设交巡台位置示意图(方块所示区域)对于P点发生重大刑事案件,动用全市警力进行围堵,我们希望使得包围圈尽可能的小,由于犯罪嫌疑人的车速度未知,我们分以下三种情况进行考虑:1)当犯罪嫌疑人的车速与警车速度同,即.运用穷举法,对肇事车辆可能的逃逸路线进行分析,以3分钟路程为半径,找到肇事车辆逃逸的覆盖范围,如图5所示,其中实线表示可能路径,在此范围内有8、9、10、15号共4个交巡平台,保证这4个平台警力不动,组成第一范围包围圈.图5 肇事车逃逸3分钟内覆盖区域示意图 进一步分析可能的逃逸路线,调度16号交巡台到36号节点,2号交巡台到3号节点,3号交巡台到55号节点,6号交巡台到47号节点,组成第二组半包围,保证对A区的封锁.若肇事车辆经36号节点逃往16号节点,则会与16号交巡台在途中相遇.对于从32号节点经7号节点逃逸到30号和47号节点,存在从A区逃往其他城区的可能,需调动其他城区交巡台的支援.将C区119号交巡台调度到237号节点,将D区320号交巡台调度到371号节点,321号交巡台经368号、369号节点到370号节点,至此,在全市范围内实现全面封锁.2)当犯罪嫌疑人的车速比警车车速小,即,我们令 方法同1),寻找分钟逃逸范围内所覆盖的全部交巡台,经过整合分析,保持7号、8号、9号、15号共4个交巡台原地封锁,10号交巡台到34号节点封锁,6号交巡台到47号节点封锁,16号交巡台到36号节点封锁,3号交巡台经55号节点到46号节点进行封锁,2号交巡台经3号节点到45号节点封锁,在此过程中,10号、2号和3号交巡台会在途中与肇事车辆相遇.3)当犯罪嫌疑人的车速比警车车速大,即,我们令由于肇事车辆逃逸速度较快,可能会逃逸到C区和F区,故需调动C区和F区警力进行围堵. 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区交巡台进行如下调配,实现全市封锁: C区:240号交巡台调到239号节点,170号交巡台调到225节点,167号交巡台调到259节点. D区:320号交巡台调度到371号节点,321号交巡台经368号、369号节点到370号节点. F区:477号交巡台调度到501号节点,518号交巡台调到521号节点,478号节点调到527号节点,484号节点到571号节点.7.模型的评价与推广本文避免了时间复杂度较复杂的Dijkstra算法,选用Floyd算法,在求最短路径上提高了效率,代码编写简单.模型的建立思路清晰,遵循可操作性、科学性、可比性原则,该模型建立出了在较理想状态下交巡警平台的最优设置,减少出警时间,均衡工作量,提高工作效率,在遇突发事件时,可尽快实现道路封锁,给生活中交巡警平台的设立予参考,具有一定的实际应用价值,也可以应用于其他适用区域.模型的运算由矩阵、向量的运算组成,易于用数学软件求解和验证.本模型较好的解决了交巡警平台的最优选址问题,当事故发生时,交巡警可以第一时间到达事发地点,有效的改善了交巡警在执行任务中的效率,在经济迅猛发展的今天,城市加速扩张,人口迅速增长,交巡警平台的设置是平安城市的最好保障.该模型也可运用到其他最优选址问题中去,比如关于消防救援工作最优路径问题、重大生产安全事故应急救援问题、公共交通的最优路径问题等. 同时也可利用该模型算法拓展模型在其他领域的适用范围.该模型也有一定的局限性,如现实中不能时刻都保证道路的畅通性.既不能保证出警的时间总是维持在3分钟之内.忽略了实际地形对于车速的影响以及实际生活中存在的不定因素.参考文献1 徐孝凯,王凤禄,数据结构简明教程第二版,北京:清华大学出版社,2005年4月1日2 李建中,骆吉洲,华章数学译丛第二版,北京:机械工业出版社,2002年6月3 陈庆华等,组合最优化技术及其应用第1版,北京:国防科技大学出版社,1989年8月4 W.T.Tutte,Graph Theory,英国:Cambridge University Press,2001年3月1日8. 附件附件1:用Floyd算法分配个服务平台管辖区域area_a _x,area_a _y=find (location_all _daolu<=92); road_index _a=area_a _x,area_a _y; road_index _aa1=find (road_index _a (:,2)=1);a2=find (road_index _a (:,2)=2);A=road_index _a (a1,1);B=road_index _a (a2,1);c,ia,ib = intersect (A,B);csize (c) for i=1:140 tt=c (i); uu=location_all _daolu (tt,:); uu1=uu (1);uu2=uu (2); vv1=location_a _zuobiao (uu1,:); vv2=location_a _zuobiao (uu2,:); ww1=vv1 (1),vv2 (1); ww2=vv1 (2),vv2 (2); line (ww1,ww2)end for i=1:140 tt=c (i); uu=location_all _daolu (tt,:); uu1=uu (1);uu2=uu (2); vv1=location_a _zuobiao (uu1,:); vv2=location_a _zuobiao (uu2,:); ww1=vv1 (1),vv2 (1); ww2=vv1 (2),vv2 (2); line (ww1,ww2,k)end% ? Error using => line% String argument is an unknown option.for i=1:140 tt=c (i); uu=location_all _daolu (tt,:); uu1=uu (1);uu2=uu (2); vv1=location_a _zuobiao (uu1,:); vv2=location_a _zuobiao (uu2,:); ww1=vv1 (1),vv2 (1); ww2=vv1 (2),vv2 (2); line (ww1,ww2,Color,.8 .8 .8)end save data_b _problem% A区节点间的邻接矩阵load data_b _problem;matric_lingjie=zeros (92,92);xx,yy=find (matric_lingjie=0);matric_lingjie (xx,yy)=inf;for i=1:92 matric_lingjie (i,i)=0;end for i=1:140 tt=c (i); uu=location_all _daolu (tt,:); uu1=uu (1);uu2=uu (2); % 端点序号 vv1=location_a _zuobiao (uu1,:); % 第一个端点坐标 vv2=location_a _zuobiao (uu2,:); % 第二个端点坐标 % 计算端点间距离 distance=sqrt (vv1 (1)-vv2 (1)2+(vv1 (2)-vv2 (2)2); matric_lingjie (uu1,uu2)=distance; matric_lingjie (uu2,uu1)=distance; % 赋值给邻接矩阵endD,R=floyd (matric_lingjie);matric_fenkuai=D (1:20,:);for i=1:20 for j=1:92 if matric_fenkuai (i,j)>30 matric_fenkuai (i,j)=0; end endendti=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);ti (3)=text (location_a _zuobiao (3,1),location_a _zuobiao (3,2)+1.5,3);ti (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);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);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);ti (12)=text (location_a _zuobiao (12,1),location_a _zuobiao (12,2)+1.5,12);ti (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