乘公交 看奥运.docx

论文中matlab运行程序并未出结果所以这篇论文中并没有给出结果，只有思路。摘要本文要解决的问题是以08年北京奥运会为背景而提出的。人们为了能现场 观看奥运会，必然会面对出行方式与路线选择的问题。因此如何快速、高效地从 众多可行路线中选出最优路线成为了解决此问题的关键。本文以以下三个影响较 大的因素：第一是换乘次数；第二是乘车时间；第三是乘车费用。针对问题一，在只考虑公汽系统的情况下，我们考虑了直达，一次换乘，和 两次换乘的情况，并分别讨论了在不同换乘情况下的乘车时间以及乘车费用，用 matlab编程求出了任意两站点间的最佳乘车路线以及换车的地点。针对问题二，在问题一的基础上加入了地铁线路，因此增加了选择线路，同 时也增加了通过步行经过地铁转乘的站点，改进问题一中的求解方法，增加线路 和可达点求解。针对问题三，在问题二的基础上又增加了步行这种情况，在适当站点步行， 可以节省交通费用而且不会消耗过多时间。关键词：最佳路线换乘次数乘车时间乘车费用一、问题重述传承华夏五千年的文明，梦圆十三亿华夏儿女的畅想，2008年8月8日这 个不平凡的日子终于离我们越来越近了！在观看奥运的众多方式之中，现场观看 无疑是最激动人心的。为了迎接2008年奥运会，北京公交做了充分的准备，首 都的公交车大都焕然一新，增强了交通的安全性和舒适性，公交线路已达800 条以上，使得公众的出行更加通畅、便利。但同时也面临多条线路的选择问题。 为满足公众查询公交线路的选择问题，某公司准备研制开发一个解决公交线路选 择问题的自主查询计算机系统。这个系统的核心是线路选择的模型与算法，另外还应该从实际情况出发考 虑，满足查询者的各种不同需求。需要解决的问题有：1、仅考虑公汽线路，给出任意两公汽站点之间线路选择问题的一般数学模型与 算法。并根据附录数据，利用模型算法，求出以下6对起始站到终到站最佳路线。 (1)、S3359-S1828 (2)、S1557-S0481 (3)、S0971-S0485(4)、S0008-S0073(5)、S0148-S0485(6)、S0087-S36762、同时考虑公汽与地铁线路，解决以上问题。3、假设又知道所有站点之间的步行时间，请你给出任意两站点之间线路选择问 题的数学模型。二、模型假设1、所有公交线路的开班、收班时间相同。2、公车不会因为堵车等因素延长行驶时间。3、各条线路不会有新的调整与变化。4、环线可以以任意站作为起点站和终点站，并且是双向的。5、除环线以外的线路，到达终点站后，所有的人都必须下车。6、人们对换乘车次数尽量少的偏好程度总是大于对花费时间相对短和花费金钱 相对少的偏好程度。7、同一地铁站对应的任意两个公汽站之间可以通过地铁站换乘，且无需支付地 铁费三、符号定义符号表示意义LAi第i条包含初始站点的线路，i = 1,2, ,mLB j第j条包含目标站点的线路，j = 1,2； ,sLCk第k条中间线路，k =1,2, ,wailLA上的第l个站点，l = 1,2； ,m ib jrLB j上的第r个站点，r = 1,2； , tc kuLCk上的第u个站点，u =1,2,，vxi乘客在第i段线路上乘坐的站数y乘客在一次地铁线路上乘坐的总站数z1公汽换乘公汽的次数z2地铁换乘地铁的次数z3地铁换乘公汽的次数z4公汽换乘地铁的次数4.1问题一4.1.1问题分析仅考虑公汽线路的情况下，首先，需要根据题目给出的公交线路信息数据， 对每条线路进行抽象处理，将分上下行的线路、双向行驶的线路和环行线路抽象 为两条。然后，主要考虑公众最关心的乘车因素，即转乘次数。在最少转乘次数 的基础上考虑共众对其他因素的需求，按照先后顺序考虑行程时间和乘车费用， 给出供公众选用的多种参考方案。并考虑以时间为主要目标的情况下，建立最优 化模型确定任意两站点行程时间最短的方案。4.1.2数据预处理(1) 将上下行的线路分为两条。中下行路线标号需要在原标号的基础上加上520，用以 区分上行线和下行线。例如：L002的上行线为L002:S3748-S2160-S1223-S1404-S2377-S1477-S2017-S2019-S1321-S1381-S138 3-S1691-S3766-S1729-S2654-S3231-S3917-S2303-S1327-S3068-S2833-S1733-S 2113-S2636-S0012-S1968-S0004L002的下行线为L522:S0004-S1968-S0012-S2636-S2113-S2112-S2833-S0618-S1327-S2303-S391 7-S3231-S2654-S1729-S3766-S1691-S1383-S1381-S1321-S2019-S2017-S1477-S 1404-S1223-S2160-S3748(2) 如果下行线是上行线原路返回，那么存储的两行数据中的站点信息刚好顺序 颠倒。例如L001分为L001和L521L001:S0619-S1914-S0388-S0348-S0392-S0429-S0436-S3885-S3612-S0819-S3524-S0820-S3914-S0128-S0710L521：S0710-S0128-S3914-S0820-S3524-S0819-S3612-S3885-S0436-S0429-S0392-S0348-S0388-S1914-S0619(3) 如果是环行线，则将其分为顺时针线和逆时针线。如下图所示：D顺时针线路：A-B-C-D逆时针线路：A-D-C-B以L017环行路公交线路为例L017:S3748-S2160-S0732-S3078-S2808-S2816-S3028-S1123-S3029-S2764-S2543-S2742-S2533-S1839-S2751-S2755-S2937-S1929-S1007-S0940-S1907-S2085-S0609-S0483-S0604-S2650-S3693-S1659-S2962-S0622-S0456-S0427-S0582-S0577-S1895-S3648-S0668-S3081-S3078-S2082-S0683-S2160-S3748L537:S3748-S2160-S683-S2082-S3078-S3081-S668-S3648-S1895-S577-S582-S427-S456-S622-S2962-S1659-S3693-S2650-S604-S483-S609-S2085-S1907-S940-S1007-S1929-S2937-S2755-S2751-S1839-S2533-S2742-S2543-S2764-S3029-S1123-S3028-S2816-S2808-S3078-S732-S2160-S37484.1.3模型建立il0i 1,2, ,m，所有包含目标站点b j0将所有包含初始站点a的线路LA ,LA , ,LA建成一个集合S； 1 < l < n ；12m0的线路 LB/ LB , , LB建成一个集合G；1 < yt; j 1,2, , s。S la, LA2,LA,G MLB、,LB ,LB ,i 1,2, ; m ；LB b > b , , * > b ； j =1,2, , , , s。（1）直达的情况.既如图所示：从A到B的直达情况示意AB-4对于直达的情况，当 SnG 0 时，存在 LA、LB ； 1 < i < m ； 1 < j < s ；使得 LA. = LB ；艮口 LA、LB,为同一线路。此线路既包含初始站点a又包含目标站点bj,。若lo <，那么，此线路为所求直达线路。若i0 >，或者当snG=。时，考虑换乘一次的线路。（2）换乘一次的情况如图所示的情况，其中C为换乘点。当有LA和LB相交时，存在LA、LB , 1 < i < m，1 < j < s，有 a g LA.及-b r，即a ,b r为同一站点。b r g LB , 1 < l < n，1 < r < t。使得 a若l < l < n，1 < r < r，那么，从初始站点a乘坐线路LA，行驶至站点a，00il0iil即在站点b，换乘线路LB至目标站点b。即jrjWo /、a LA） a - b T LB ） bil0iil jrjjr0若不满足l0 < l < n，.1 < r <,，或者，当无任何LA和LB.相交时，考虑换乘 两次的线路。（3）换乘两次的情况如图，是换乘两次的情况，其中C,D是换乘点。-4记 LC , LC , , LC , LC = c r c r r c , k = 1,2, , w，有LC wS , 12wk k 1k 2kvkLCk w G , k = 1,2, , w，且满足LCk与LA、LB都相交时，即线路LC既不包含初始站点a又不包含目标站点b , 1 < l < n , 1 < r < t。 kilojrQ00但是存在c e LC及a e LA，使得c = a，kUk ilikUil存在七e LC*及b e LB，使得匕=b ,即c 、a.为同一站点，且c 、b为同一站点。1 < k < w， 1 < i < m， 1 < j < s， 袂 ilku2jr1 < u < v，1 < u2 < v，1 < l < n，1 < r < t。若l < l < n， 1 < u < u < v， 1 < r < r，那么，从初始站点a乘坐LA线路, 0120il0i行驶至站点a，即在站点c ，换乘LC线路至站点c ，即在站点b，换乘LBilk%kku2jrj线路至目标站点b 。即J、arLA )r a cr(LC )-r c brLB)r bil0iilkUkku2jrjjr0若不满足l < l <n，1 < u < u < v,1< r < r ，或者，当不存在满足条件的LC时, 0120k说明需要换乘三次才能够到达目标站点。（4）最优化模型的建立x = 0;0 <x < 20;20<X <40;，ix > 40设f为乘坐公交线路的费用函数:;f （勺十、3,总时间函数：（1）（2）T = £ x + 5z(0 < z < 2)i=1总费用函数：F = £ f (x )ii=1其中气表示乘客在公交线路七上乘坐的站数；Z1表示公汽换乘公汽的次数。目标：找出任意给定的两站点的乘车线路，使广和尸相对最小。4.1.4算法的解释由于人们的对换乘车次数尽量少的偏好程度总是大于对花费时间和金钱相 对少的偏好程度，我们将优先考虑换乘车次数尽量少，然后再考虑花费时间相对 短、花费金钱相对少，对得出的所有结果中进行筛选。对于整个算法对路线搜索的具体算法步骤为：第一堤输入乘车始点站和终点站，在数据库的站点线路矩阵中查询是否有相 同路线的公汽(1) 有相同路线公汽，则查询站点间有直达线路，输出所有可直达车辆信息， 结束程序；(2) 没有相同路线公汽，则查询站点间需要转乘，进入第二步。第二步；将始点和终点输入公汽站点矩阵，查询是否有相同的站点(1) 有相同站点，则查询站点可以通过一次转乘达到，输出所有一次转乘的车 次和中转站的信息，结束程序；(2) 没有相同站点，则查询站点间需要转乘两次，进入第三步。第三步：找最后一个中转站，在终点站所在线路中沿车驶向的反向找一个站点 当做终点站称为伪终点站，根据伪终点站和始点站输入到公汽站点矩阵，查询 是否有相同的站点(1) 有相同站点，则查询站点可以通过两次转乘达到，输出所有两次转乘的车 次和中转站的信息，结束程序；(2) 没有相同站点，则查询站点间需要转乘三次(此论文不考虑换乘3次及3 次以上的情况) 以上的搜索方法得到公交查询系统的流程，具体流程图如下输入起始点4.1.5模型的求解根据以上算法和前面建立的模型一，用matlab对程序运行（程序见附录） 就可以得出不同目标下的最优路线。4.2问题二4.2.1问题分析针对问题二，将公汽与地铁同时考虑，找出可行路线，然后寻找最优路线。 对于地铁线路，也可以将其作为公交线路，本质上没有什么区别，只不过乘车费 用、时间，换乘时间不一样罢了。因此地铁站可等效为公交站，地铁和公交的转 乘站即可作为两者的交汇点。因此该模型的公交换乘路线模型与模型一中的基本 相同。4.2.2算法分析由于假设同一地铁站对应的任意两个公汽站之间可以通过地铁站换乘且无 需支付地铁费，那么不妨把同一地铁站所对应的几个公汽站合并成一个站。地铁线路T1 = D 01 c D 02 c c D 23，T2 = D24 c D25 c c D39 c D24。1、可以乘坐地铁的线路。（1）若初始站点和目标站点都在地铁线路T1或者T2上，那么，只乘坐地铁T1或 者T 2便可以直达。其中，若都在线路T 2上，就选择经过站数最少的方向。若初始站点和目标站点分别在地铁线路T1和T2上，那么，需要进行一次地 铁换乘地铁才能到达。（2）若只有初始站点或只有目标站点在地铁线路上，则需要换乘公汽才能到达 目标站点。初始站点a0 w Tp，p = 1,2，目标站点b史T1且b史T2，b e LB。当有LB和地铁相交时，即存在LB，有b . eLB，使得b . e Ty，q = 1或2。1 < i < m， 1 < . < s。若1 < r < %，那么，从初始站点a0 （记为Da）乘坐地铁线路，行驶全站点b （记为Db），换乘公汽线路LB至目标站点b。1 < a < 39，1 < b < 39。即0ail0Da T (jp ) T bDb Tjr 1-1Da T (jp)(Iq)(LB ) T b( q = p )jjr0-T b Db T LB.) t b(q 丰 p )其中，q壬p时需要地铁换乘地铁。若不满足1 < r < %,或者当没有这样的LB时，说明在地铁换乘公汽后，还 需要进行公汽换乘公汽。由于这样的情况几乎不存在，故不作考虑。目标站点b gTp，初始站点a任T1且a任T2，p = 1,2jr0il0il0同理可得结论。（3）若初始站点和目标站点都不在地铁线路上，则先乘坐公汽，换乘地铁, 再由地铁换乘公汽。地铁线路既和LA.相交又和LB相交时，即地铁线路既不包含初始站点ai又不包含目标站点b。但是存在LA.、LB ,1 < i < m， 1 < j < s，有a g LA，使得 a g Tp，记 a为 Da，b g LB，使得 b g Tq，记 b 为 Db，p = 1,2，q = 1或2，1 < a < 39，1 < b < 39。< l < n，1 < r < r，那么，从初始站点a乘坐LA线路，行驶至站点a （记0il0iil为Da），换乘地铁线路至站点bjr （记为Db），换乘LB,线路至目标站点b。即bjr0, (LB )j(q = p)bjr0a (LA ) a Da (Tp) b Db LB )il0iiljrjL S)- *0 (单)(Tq) - b D后其中，q丰p时需要地铁换乘地铁。若不满足l0 < l < n，1 < r <%，或者不存在LA、LB都与地铁线路相交，说 明需要在地铁线路前或后进行公汽与公汽的换乘。由于这样的情况几乎不存在， 故不作考虑。2、只乘坐公汽的线路。完全排除地铁线路，与解决问题一的方法相同。4.2.3模型建立设f， g分别为乘坐公交和地铁线路的费用函数:'0,1,f 巴2,。，x = 0;0 Z x < 20;20 < X < 40;x > 40T = 3 x + 2.5y + 5z + 4z + 7z + 6zi1234i=1总费用函数：(0 < zi < 2，X z. < 2)i =1（3）总时间函数：（4）F = X f (x ) + g (y)ii =1其中xi表示乘客在公交线路L上±乘坐的站数；y表示乘客在一次地铁线路上乘坐的总站数；z ,z ,z ,z分别表示公汽换乘公汽，地铁换乘地铁，地铁换乘 1234公汽，公汽换乘地铁的次数。目标：找出任意给定的两站点的乘车线路，使T和F相对最小。4.2.4模型求解首先，我们通过matlab编程（程序见附录）作出了两条铁路的位置关系图， 如下图所示。T1与T2铁路位置关系图同样将地铁线路等效为公交线路得出任意两个站点间的可行线路，再将目标函数分别 用上述模型建立的模型表达式表达，用matlab进行编程（程序见附录）求得出 考虑地铁情况的最优路线。4.3问题三4.3.1问题分析问题三是在前面问题的基础上，加入了步行这一较为自主化的“交通工具”， 使得原本的选择最优线路模型不再使用，于是我们这里建立了一个线路综合评价 模型，通过分类讨论的方式，提供适合各种情况的线路选择方案，从而解决在三 种交通工具并行时的路线选择问题。4.3.2模型建立我们使用综合评价的方法对不同的情况进行分类，分类的标准是：乘坐的交 通工具。由于这里有三种交通工具：步行、地铁、公共汽车，所有我们有 种不 同的搭乘方式，可分为三类：搭乘一种交通工具，搭乘两种交通工具，搭乘三种 交通工具。下面我们来分别讨论。搭乘一种交通工具 搭乘公共汽车这时可以根据问题1的模型来解决选择线路问题。 搭乘地铁这种方式只适用于起始站、属于地铁的相邻公汽站台矩阵，这时可以通过在 地铁直达。 步行步行可以到达任何终点站，但从实际出发，我们认为步行距离不能过长，这里我 们认为当起始站、相隔不到两站时，可以采取步行；起始站相隔超过两站时， 我们认为步行不实际。搭乘两种交通工具 搭乘公共汽车和地铁这时就是等同于运用问题2的模型来选择线路方法。 搭乘公共汽车和步行这时我们可以在问题1的基础上考虑加入步行。这里我们只提供在以下情况下采 取步行较合适：a. 为减少换乘次数：两公汽站台只相隔一站，可采取步行，从而减少了一次的换 乘；b. 为减少路费：因为分段计价的公共汽车票价为：020站：1元；2140站：2元；40站以上：3元。所以在终点站或换乘站在21站处时，我们可提前在20 站出下车，再步行1站到达终点站或换乘站，从而节省1元；c. 当出现两站之间出现如图示情况时，可考虑不行到达终到站。 搭乘地铁和步行a. 当起始站、不属于或不全属于地铁的相邻公汽站台矩阵时，需要通过步行 才能到达终点；b. 由于地铁只有T1、T2两条线路，且两条线路只能在D12和D18处转乘，因此 在需要换乘地铁时，可考虑步行到达地铁站D12或D18。搭乘三种交通工具(地铁、公共汽车、步行)已知所有站点之间的步行时间，假设任意两站点，、j之间的步行时间为，通过转乘从z站到达j站所用时间，转乘用时t，则我们可在上述的基础上加入一个 0判断语句，即t +10与T比较大小，如果t +10较小，则说明通过转乘更省时；反之，步行更省时同时也省钱，则采取步行从，站到达j站。(4) 最优化模型的建立设f,g分别为乘坐公交和地铁线路的费用函数：0,1,2,3x n = 0; J0, y=。； g=k|3, y > 0.0 < X n < 20;20 < X n < 40;x n > 40.根据实际情况，在地铁线路上不考虑步行。我们可以在初始站点、目标站点 或换乘站点的附近考虑步行，即在任意公交线路L , 1 < i < 3上最多下车一次。否则，若在某个L ,1 <i< 3上下车步行两次，则在L上需要多购买车票一次，同时 消耗的时间更多，此做法既违反常理，又不经济实惠。设在线路L, i 1,2，止步行的站数为n, 0 n x，相邻公汽站步行时间为t,那么总时间函数：T 3 3 (x. n)i 1总费用函数：F3-，t n 2.5y 5z 4zi 17z 6z ,(5)f (x n ) g (y),i 1目标：找出任意给定的两站点的乘车线路，使T和F相对最小。五、模型评价及推广本题中所构建的模型简单易懂，操作简单，涵盖了所有路线的选择情兄，并 且相对符合“乘公交，看奥运”的主题，解决了公交线路的选择问题，使公众的 出行更加通畅便利。但是本模型忽略了人流、车流拥挤的状况。因此在实际生活 中有限制。对于若干条从某一初始站点到目标站点的线路，我们可以设计一种带记忆功 能的系统，即乘客选择某路径的次数越多，说明此路径是比较优的路径，为以后 选择路径提供必要的信息。系统使用的时间越长，为乘客提供的信息越全面，越 准确，系统也越智能化。这样就可以为乘客需求量最大的一条增加班次，以满足 更多人的需要。在假设中提到，所有线路的开班、收班时间相同，但事实并非如此。那么可 以在模型的设计中加入线路运行的时间元素，使乘客查询时只显示正在运行的线 路。六、参考文献1大学数学实验姜启源，邢文训，谢金星杨顶辉，北京：清华大学出 版社，20004MATLAB工具箱应用苏金明编七、附录问题一的程序代码(直达的线路)x1=input('please input starting station:');y1=input('please input the terminal :');i1,j1=find(a=x1);i2,j2=find(a=y1);m,n=size(i1);p,q=size(i2);r=0;for i=1:mfor j=1:pif i1(i,n)=i2(j,q)nv=find(x1=a(i1(i,n),:);nu=find(y1=a(i2(j,q),:);if nv<nur=r+1;t(r)=i1(i,n);endendendendif r=0disp(t)elset=0endj1j2%直达的输出说明t是线路j1是起点站在该线路的第几个站j2是终点站在该线路的第几个站问题一的程序代码(换乘一次的线路)x1=input('请输入起点站:');y1=input('请输入终点站:');W=input('输入最多经过站点的个数:');i1,j1=find(a=x1);% 记录行和列i2,j2=find(a=y1);m,n=size(i1);p,q=size(i2);for i=1:mfor j=1:pro=0;if i1(i,n)=i2(j,q)mv=a(i1(i,n),:);mu=a(i2(j,q),:);mo,no=size(mv);po,qo=size(mu);for io=1:nofor jo=1:qoif mv(mo,io)=mu(po,jo)ad=find(a(i1(i,n),:)=x1); %x1 所在的位置 bd=find(a(i2(j,q),:)=y1);%y1 所在的位置ao=find(mv(mo,io)=a(i1(i,n),:); %转站点在 x1 所在列的位置bo=find(mv(mo,io)=a(i2(j,q),:); %转站点在 y1 所在列的位置 if ad<ao&bo<bd&(ao-ad+bd-bo)<Wro=ro+1;to(ro)=mv(mo,io);tka(ro)=ao-1;tji(ro)=bo-1;endendendendif ro=0disp('中转站点')disp(to)disp('中转站点在始发线上的位置)disp(tka)disp('中转站点在抵达线上的位置)disp(tji)vo(1)=i1(i,n);vo(2)=i2(j,q);disp('始发线和抵达线)a(vo,1)disp('起点站位置')ad-1disp('终点站位置')bd-1endendendend问题一的程序代码(换乘两次的线路)x1=input('请输入起点站:');y1=input('请输入终点站:');W=input('输入最多经过站点的个数:');i1,j1=find(a=x1);i2,j2=find(a=y1);m,n=size(i1);p,q=size(i2);vp,vb=size(a);tto=0;%寻找不包含起点和终点的线路for iu=1:vpvc=a(iu,:);rpp=find(x1=vc);rpq=isempty(rpp);tpp=find(y1=vc);tpq=isempty(tpp);if rpq=1&tpq=1tto=tto+1;uu(tto)=iu;endendfor ey=1:size(uu,2)eyy=a(uu(1,ey),:);for ex=1:mexx=a(i1(ex,n),:);for ez=1:pezz=a(i2(ez,q),:);mn=size(exx,2);iq=0;ih=0;%寻找exx和eyy的相同元素，赋值到tqfor i=1:mnfor ig=1:mnif exx(i)=eyy(ig)iq=iq+1;tq(iq)=exx(i);endendend%寻找eyy和ezz的相同元素，赋值到tnfor i=1:mnfor ig=1:mnif ezz(i)=eyy(ig)ih=ih+1;tn(ih)=ezz(i);endendendif iq=0&ih=0 %eyy 与 exx，ezz 都有交点for di=1:iqwq=find(tq(di)=eyy); %exx与eyy的相交元素在eyy中的位置 for dh=1:ihwh=find(tn(dh)=eyy); %ezz与eyy的相交元素在eyy中的位置if wq<whLe=find(x1=a(i1(ex,n),:)-1;Lf=find(eyy(wq)=a(i1(ex,n),:)-1;%起点站和第一个转站点在第一条线路上的位置Lg=find(eyy(wq)=a(uu(1,ey),:)-1;Lh=find(eyy(wh)=a(uu(1,ey),:)-1;%中转站点在第二线路上的位置Lj=find(eyy(wh)=a(i2(ez,q),:)-1;Lk=find(y1=a(i2(ez,q),:)-1;%第二个转站点和终点站在第三线路上的位置if Le<Lf&Lg<Lh&Lj<Lk&(Lf-Le+Lh-Lg+Lk-Lj)<35disp('起点站第一个转站点第二个转站点终点站在所在线的位置) weizhi=Le,Lf,Lg,Lh,Lj,Lkdisp('转两次站经过的路线)；L=a(i1(ex,n),1),a(uu(1,ey),1),a(i2(ez,q),1)disp('中转站点')；zhongzhuan=eyy(wq),eyy(wh)endendendendendendendend问题二的程序代码y1=input('finish');i1,i2=find(x1=a)tt=0;t=size(i1,1);for p=1:tmv=a(i1(p,1),:);n=size(mv,2);for tp=1:nif(mv(tp)>5000&mv(tp)<6000)tt=tt+1;mm(tt)=a(i1(p,1),1);breakendendenddisp('起点经过地铁的线路') disp(mm)ttt=0i1,i2=find(y1=a);t=size(i1,1);for p=1:tmv=a(i1(p,1),:);n=size(mv,2);for tp=1:nif(mv(tp)>5000&mv(tp)<6000) ttt=ttt+1;mmm(ttt)=a(i1(p,1),1);breakendendenddisp('终点经过地铁的线路')disp(mmm)