列车运行图MATLAB程序.doc

%第一个
%设６：４2时间为0时刻，转换如下．
a1=［6　7 ７ 7 7 7 7 8 8 9　9 9 ９　９ 10 10 10　10　10　1０ １1 11 12 1２ １2 １２　13 13 14　14　1５ 15 １5　16 １6　16 16 17　17 17 １7 １8 18　18 １８ 19 ２0 ２1］；
ｂ1=[42 0 1０　２2　33 40　５2 15　30 12 17 2２ 27 4０　5　1０ 15 2８　4５ 50　33 43　7 13 3４ 4３ 5 18 ２5　43　６　40 45 0 18　３5　５７ 3 1９ 3２　37 ２　１2 ３６ 42　22　47　6］；
ｃ1=7*ｏｎes(1,lｅｎgth(a1）)；
ｄ１=６0*（a1—c1)+b1＋１8；
e1=[7 7 7　8 8　8 8　8 9 9 9 10 10　10 １０ 1０ 11 1１ 11　11 1２　12　12　12 13　1３ 13　1４ 1５ 15 15 16 16 １6 16 17　17 18 １8　18　18 1８ 1９ １9 1９ 20 21 21］;
f1=［２3 41　58　3 １5 ２8 37 5６ １1　53 ５8 10　15 21　46 ５１　3 9 2６ 31 １４ 2４ 48　５４ 15　24 46 15 6 24 ４7　2１ 26 48 ５9 23 38 0 7　1３ 18 43 0　1７　30　17 ３５　４７］;
g=7＊onｅs(1,length（a1）)；
ｈ1=60*（e1—g）+f1＋18；
%输入经历时间，求众数,得到一个差值矩阵f。
x1＝［４1
41
48
41
42
4８
４5
41
41
4１
41
４８
48
41
41
41
48
41
4１
41
4１
41
41
４1
４1
４1
４１
57
4１
41
4１
４1
4１
48
4１
４8
41
５7
4８
41
41
４１
48
41
４8
5５
48
41］’;
v１=mode(x1)
j１=x1-v1;
disｐ(j1）
%矩阵剔除，将f中大于0得用0代替得以矩阵m，将ｆ小于0得数代替得一矩阵ｎ，ｎ，ｍ都就是１*２０得矩阵。
n1=（j1+abｓ(j1）)/2；
m１=j1；
ｍ1（fｉnd(j1＞0)）＝0;
%选取可行列车
z1＝ｚeros(1，ｄ１（enｄ)-d1(１）+１）;％(或zeros（１,d（end）—ｄ(1））)
y11=d1(1）：１：d1(enｄ）；
foｒ k=１:ｌength(x1)—1
fｏr i=1：ｄ1(end）—ｄ１(1)+1
iｆ d１(ｋ+1）—ｙ1１(i)+ｍ1（k＋１）〉＝５ ＆＆ y11(i)—d１(ｋ）-n1（k）〉=5;
　 z1(i)＝y１1(ｉ)；
　　ｅｎd
enｄ
end
%将满足条件得列车再进行剔除,求数量最大且满足条件得列车
r=０;
ｇ=0;
k＝1;
p=1;
l１=zeros(1,lengtｈ(ｄ１)+5０）;
t1=zeros(1，d１(end）—d1(１)＋1）；
foｒ i=1：d1(eｎd)—d１（1)
if z1(i)+1==ｚ1（i+1)
　 ｋ=ｋ+1;
ｅnd
iｆ z１(i）〉ｚ1（ｉ+1）
　　 r＝fｌｏor（（k—1）／5)+1;
　 fｏｒ u=1:ｒ
　 g=i—k+１+5＊（ｕ-1);
　 　 　 　 if　g〉d1（end)-d1（１)+1
　 　　ｇ=d１(end)－d1(1）+1;
　 　　 end
　　 　 l１(g)=z1（g);
eｎd
ｋ=1；
eｎｄ
end
l1（fｉｎｄ(l１＝=0））=［]；%找到l中0得位置，并令其为空，即删除
w1=l1+ｖ１;
y１3=0、1*ones（1,lｅngｔh(l1));
y1４=41*oｎes(１，lｅngth(w1));
for ｉ＝1:lengｔh(l1）
　 ｐlｏt（［ｌ1（ｉ),w1(ｉ)]，[ｙ13（ｉ）,y14（i)]）；
holｄ ｏn
eｎｄ
y11=0、1＊onｅｓ(１,ｌeｎgth（d1）)；
y12＝4１*oｎes(１，lenｇtｈ（h1)）；
fｏｒ i=1：lｅngth(d1)
pｌoｔ([d1（i),ｈ１(i)］，[y11(i),y12(i）]);
hold oｎ
ｅnd
%第二个
%设6:42时间为0时刻,转换如下。
ａ2＝[8 8　８ １０ 10 1０ １１ １３ 1４　15　16 19 20］;
b2＝[6　4４ ５8 0 24 53 16 32　23 58　28 19 １9];
c2=7＊ｏnｅs(1,leｎgth(a2））;
d2=60*(ａ2—c２)+b2＋18;
e2=［8　9　9 10 １0 11　11　１３ １4 1６ 16 19 ２０］;
f2=[2４　４　16 18 42 11　３4 50 4３　16　４6　37　37］；
g=7*oｎeｓ（１，ｌeｎｇｔh（a２）);
h１=６0＊(e２－g)+f２+18;
%输入经历时间，求众数，得到一个差值矩阵f。
x2=[１8
20
18
18
18
1８
１８
18
2０
１8
18
1８
１8
]';
v2=modｅ(x２)
j２=ｘ2-v２;
disp(j2）
%矩阵剔除，将f中大于0得用0代替得以矩阵m,将f小于0得数代替得一矩阵n,n,ｍ都就是1＊20得矩阵．
n2=(j２＋ａbs(j２））/２;
m2=j2;
ｍ2(ｆind（ｊ2〉0)）＝0;
%选取可行列车
z2＝zeｒoｓ(1，d2(enｄ）-d2(1)+１);％（或zeros(1,ｄ(end)—d(1)))
y21=d2(1）:1:d2(eｎd);
foｒ k=1:lengtｈ（x2）-１
ｆor　ｉ=1:d2（end)-ｄ2（１)+1
　if d2（ｋ+1)—y2１（i）＋m2（k+１)>=5　＆& y２1(i）—d2（k)－n２(k)>=５;
　　　 z2(i)=ｙ21(i）;
　 　enｄ
　　ｅｎd
eｎｄ
%将满足条件得列车再进行剔除,求数量最大且满足条件得列车
r=0;
g=0；
k=1;
ｐ＝1;
l2＝ｚｅros（1,lenｇth（ｄ2）+5０）；
t2=zｅｒos（1，ｄ２(eｎd)—d2（1)+1);
fｏr i=1：d２(enｄ)—d２（１)　
if　ｚ2（i)+１=＝z２（i+1）
k=ｋ+1;
ｅnd
if　z2(i）>z2(ｉ+1)
　 r=ｆｌoｏr((ｋ—1)／5)+1;
fｏr u=1：ｒ
　 　　ｇ=i-k+１＋5＊（u－１)；
if g〉d２（enｄ)—d2(１）+1
　　　 　 　 g=d2（ｅｎd)—d2(1）＋1;
　 end
　 　 ｌ2(ｇ）＝z2(g);
　 end
k=1；
ｅnｄ
end
ｌ2（find（l2==0)）＝［］;％找到l中0得位置,并令其为空,即删除
w2=l2+v2；
y２3＝41＊ones(1，leｎｇth(ｌ2)）;
y24=59*ones(1,ｌength(w2）);
fｏr ｉ=1：lｅnｇth(l2）
　 　pｌot([l2（i)，w2(i)],[y２３(i）,ｙ2４(i)］）；
hｏld　ｏｎ
end
y２1=41*ones(1，leｎgth(d2））；
ｙ22=59*ones（1，length（h1)）;
ｆor i＝１：ｌeｎgth(d2)
plot(［d2(i),h1(ｉ)],[y2１（ｉ)，y2２（ｉ）］);
hold oｎ
eｎd
％第三个
%设６:42时间为0时刻，转换如下．
a3=[7 8 9 ９ 10　10 11 11　12　1３ 14 1５ １6 16 17 18 １9 19　１９ ２0］;
b3=［５０ ２６ 6　１8 ２0 50 13 3６ ５6　5２　5２ 51　28 48 ２6 1０ 10　39 5８ 3９］;
c3=7*oneｓ（1,length（ａ３））;
d3=60*（a3—ｃ3)+b3+1８；
e3=［8　８ ９ 9 1０ １1 11 11　１３ 14 15 １6 16 1７ 1７　18　19 20 20 21];
ｆ3=[１３ 4９　31 41　43 １３ 36 59 19 15 16 21 ５1 １1 ４９　33 ３3 2　21 2］；
g=7*ｏnes（1，length(a3))；
ｈ３=60＊（e3—g)+ｆ3＋18;
％输入经历时间，求众数，得到一个差值矩阵f。
ｘ3=［23
２3
2５
23
23
２2
23
23
23
23
２4
30
23
22
23
23
2３
２3
2３
23]’;
v3=mｏｄe(x3)
j３=ｘ3-ｖ3；
disp（j３）
%矩阵剔除，将f中大于０得用0代替得以矩阵m,将f小于0得数代替得一矩阵ｎ,n,m都就是１*20得矩阵.
n３=（j３+abｓ(j3)）/2;
m３=j３;
m3(ｆind（ｊ3＞0))=0;
%选取可行列车
z3=zerｏs(1，d３(end）-d3(1)＋1);％(或zeroｓ（1,d（eｎd)—d(1)）)
y31=ｄ3（１):1：d３(end);
ｆｏr ｋ＝１：lｅngｔｈ（x3)-1
for i=1:d3(enｄ)-d3(1)+1
ｉｆ ｄ3(k+1）—y３1（ｉ)＋m3(k+1)〉＝５ ＆& y31(i)-d3(k)—ｎ3(k)〉=5;
　ｚ３（i)＝y31（ｉ)；
enｄ
　ｅｎd
ｅnｄ
％将满足条件得列车再进行剔除,求数量最大且满足条件得列车
r=0;
g=0;
ｋ=1；
p＝1;
l３=ｚｅros(１,ｌengtｈ(ｄ3）＋5０）;
ｔ3=zeｒｏｓ(1,d3(eｎｄ)—d3(1)+1)；
ｆor　　i=1:d3（ｅnｄ)—d3(１）
ｉf z3(i)+1=＝z3（i＋１)
　 k=k+１;
enｄ
if z3(i)〉z3（i+１）
　 r=fｌoor((ｋ—1)/5)+1；
　for　ｕ=1:ｒ
　 　g=i—k+１＋5*(u—1)；
　 　　 　if ｇ〉ｄ3(eｎd)-ｄ３(１)+1
　 　 　　g＝ｄ3(eｎd）-ｄ3(1）+1；
　 　 end
　 　　l3(g)＝ｚ３(ｇ);
　ｅnｄ
　 k=1；
enｄ
end
l3(fiｎd(l3==0）)=[］;％找到l中0得位置,并令其为空,即删除
w3＝ｌ3+v3;
y33=59*ones(1，ｌength（l3))；
ｙ34=８2*oｎes（1,ｌenｇth（w3）);
for i=1：length（l３)
　　 plot([ｌ3(i),w3(i)]，［ｙ3３(i),y３４(i)]);
hｏｌｄ oｎ
end
y31=59*onｅs(1,length(d3));
y3２＝８2*oneｓ(1，leｎｇth(ｈ３))；
for i=1:length（d3)
pｌot([d3(i),ｈ３(ｉ）］,[y３1（i)，y３2（i)］）;
hold on
end
%第四个
%设6:42时间为0时刻,转换如下。
ａ4＝［7　8 ８ ８ ８ ８ 9 9 9　10 10　11 1１ １１ 13　13　13　1３　1５　１5 16 16 １6 18 １8 19 2０ 20］；
b4=［0 0 21 ２8 48 57 9 15　２5　51 ５9 16　39 5７ 21 32 4２　5６ 7 47 ５ 24　５４　４ 36　42　9 39];
c4=7＊ones（1,length（ａ4)）;
d4=６0*(a4—c４）+b4+１8；
ｅ4=[7 8　8 8 9 9　9 9 10 1１ 1１　１１ 12 12 1３　１4 14 14 15　1６ 1６　16 1７ 18 19 20 ２０ 21];
f4＝[27　2７ 49 56　16 ３2 37　42　0 19 27 ４４ 1５ ２5 ４9　0 10 24 ３５　1５　３3 5９　2９ 3１ 4　１0　3７ 7]；
g=7*oneｓ(1,length(a４));
h4＝6０＊(e4-g）+ｆ４+18;
%输入经历时间,求众数,得到一个差值矩阵ｆ.
ｘ4=[２7
27
２8
２8
２8
35
28
27
35
2８
２8
28
36
２8
2８
２8
28
28
２8
28
２8
3５
35
27
２8
２8
２８
28]’;
ｖ4＝mｏｄe(x４)
j4＝ｘ４—v4;
dｉsｐ(j４)
％矩阵剔除,将f中大于０得用０代替得以矩阵m，将ｆ小于０得数代替得一矩阵n，ｎ，m都就是1＊20得矩阵。
ｎ４=（ｊ４+abs（ｊ4）)／2;
m4=j4；
ｍ4(fｉnd(j4＞０))=0；
%选取可行列车
z4=zｅroｓ(1,d4（ｅnｄ）—ｄ4（1)+1);%（或zｅrｏｓ（1，d（enｄ)-ｄ(1)))
y41=ｄ4（1）：１:d4(enｄ);
for ｋ＝１:lenｇth（x４)-１
for i=1:d4(ｅnd)-ｄ４(１）＋１
if d4(k+1）-y4１(i)+m4(k+１)>=5 && y41（i)-d４(k）-ｎ4(k）>=5;
　　 z4(ｉ）=ｙ41(i）;
　 ｅnd
　　　 end
ｅnｄ
%将满足条件得列车再进行剔除,求数量最大且满足条件得列车
r=０;
ｇ＝0;
k=1;
ｐ=1;
l4=ｚeｒｏｓ（１，ｌengｔｈ(d４）+50);
t4=zerｏs(１，ｄ4（enｄ)—d4(1)+1)；
for i=1：d4（end)—d4(1) 　
if z4（i）+1==z4（i＋1)
k=k＋1；
end　
if z4(i）〉ｚ4(i+1)
　r=floｏｒ(（k-1)／5)+1；
for u=１:r
　 　 　g=i-k+1＋5＊(u-1);
　 　if g>ｄ4(ｅnｄ)—d４(１)+１
　 　 ｇ=d4(ｅｎｄ)-d4(１）+1;
　 end
　　 　 　 l4(g)=z4(g）;
　 eｎｄ
　 　k=1；
enｄ
eｎd
ｌ4(find（l4==0））＝［];%找到l中0得位置,并令其为空,即删除
ｗ4=l4+ｖ4；
y43=８2*onｅｓ（１,ｌenｇth(l4）);
y44=1１０*ones(1,lengtｈ(w4）);
ｆoｒ　i=1：lengｔｈ（l4）
plot([l4(i),w4(i)］，［y４3(i）,y４4（i)］);
holｄ　ｏｎ
ｅｎd
y41=8２*oｎes（1,ｌenｇｔh(d4））;
y４２＝110*oneｓ(1,length（h4)）；
ｆor i=1：lｅngth（d4）
ｐlｏt([d4(ｉ）,h4(i)],［ｙ41(i）,y42(i)])；
ｈｏld　oｎ
end
％第五个
%设6:４２时间为0时刻，转换如下。
ａ5=［7 8 8　9 9 9　10 １1　１2　1４ 14　１5　１6　1６ 17 1８］;
b５=[29　29 ５1　19 3９ ４4 2 ２1 18 1９ 34　39 １7 3５　45 33]；
ｃ5=7*onｅs(１,lengtｈ(ａ５)）;
d5＝６０*（ａ５—c5)+ｂ5+１8;
e5=[７ 8 9 ９ 9 10　10 1１ 12　１4 14 15 16　16 １8 1８］;
f5=[4５ 45　7 ３５ ５5 0 18 3７ 34 ３5 5０ 5５ ３3 51　1 ４9];
g=7*ones(1，leｎgth(a5）);
h5=６0*(ｅ5-g）＋f5+18；
％输入经历时间，求众数,得到一个差值矩阵f。
x５＝[１6
１６
１６
16
16
１6
１6
16
１6
16
16
１6
16
16
16
1６]＇;
v5=mode(ｘ５)
j5=x5—v5;
dｉsｐ(j5）
％矩阵剔除，将ｆ中大于0得用0代替得以矩阵m,将f小于0得数代替得一矩阵n,n,ｍ都就是1＊２０得矩阵。
ｎ5=(j5+aｂs(ｊ5)）/2;
ｍ5＝ｊ5;
m５(fiｎd(j5〉0))=0;
％选取可行列车
z５＝ｚｅros（1,d５(ｅnd)-ｄ5(1)+1);%(或zｅroｓ(1,ｄ（eｎd）－d（1））)
y51=d5（1):１：d５(end）;
for k=1:length(x5)—1
foｒ　i=1:d5(eｎｄ）—d５（1)+1
　iｆ ｄ5（ｋ+１）-y5１(i)＋m5(ｋ+1）〉＝5 ＆＆ y5１(i)—d5(k)—n５(k)〉=5；
z５（ｉ)=ｙ51(ｉ）;
　 　 enｄ
　 enｄ
enｄ
％将满足条件得列车再进行剔除,求数量最大且满足条件得列车
r＝0;
ｇ=0；
ｋ=1;
p=1;
l5＝ｚeroｓ（1,lengtｈ(d５)+５0);
t５=zｅｒos（１,d5(end)—d5（1)+１）;
for i=1:ｄ5（end）-ｄ5(1)
if z５(i)+1＝＝z5(i+1）
　 　 k=k＋1;
ｅｎd
if ｚ5（i)>z5(i+1)
　 r＝flooｒ（(ｋ—１)/5）+1;
　 　fｏr　ｕ＝1：r
　　 　g=ｉ—k＋１+5*(u—１)；
　 　 if　ｇ>d５(ｅnd)—ｄ5(1）+1
　　 g=d5(end）－d5(1）+1；
　 　 enｄ
　 l5(g）=z5(g);
enｄ
　k=１；
end
ｅnd
l5(find(l５==0))=［];%找到l中０得位置，并令其为空,即删除
w５=l５+v5;
ｙ53=１10＊onｅs（1，lenｇｔh（l5)）;
ｙ5４=１2６*oｎes(１,length(w5))；
for ｉ=１：ｌｅnｇｔh（ｌ5）
　 pｌot（［l５(i）,ｗ5（i)]，[ｙ5３(i),y54(i）]);
ｈoｌd on
end
y５１=1１0*ｏnes（1,leｎgｔh(d5))；
y52＝1２6*oｎes(1,lengtｈ(ｈ5));
ｆｏr i=１:ｌength（ｄ5)
plｏt（[d5（i),h５（i)],[y５1(ｉ),y52(i)]）;
ｈｏｌd on
end
%第六个
%设6:4２时间为0时刻,转换如下。
a6=［8 9 10 １0　1１　13 13 １４ 15 15 16 16 １7 18 18 １8 19　１９ 19]；
b６=[58 1６ ２ 46 3 23　４7　53 16　５8　41　５3 30 ３ 11 34 3０ 40　56］；
c6＝7*onｅs(1，lｅｎgｔｈ(a６）);
ｄ6＝6０*(a6—c6）＋ｂ6+1８；
e6=［9　9　10 11 １1　１３ 14　15 1５ 16 1６ 17 17 1８ 18 18　19 19 ２０］；
ｆ6=［15 33 20 3 2０ ４0　4 10 3３ 15　５8 １０ ４7　2０ 28 51 47 5７ １３]；
g=7＊oneｓ（１，lｅngｔｈ(a6))；
ｈ6=60＊(ｅ６-g)+ｆ6+18;
％输入经历时间,求众数,得到一个差值矩阵f。
x6=[17
17
18
17
17
１7
1７
17
1７
17
17
1７
１7
17
1７
17
1７
17
17]＇;
ｖ6=ｍode（ｘ6)
j6=x6－v6；
disｐ(ｊ6)
％矩阵剔除，将f中大于0得用０代替得以矩阵ｍ，将f小于0得数代替得一矩阵ｎ,ｎ,m都就是1*2０得矩阵.
ｎ6＝（j6+abｓ(j6)）/2;
m6=j6；
m6（ｆind(ｊ６>０))=0;
%选取可行列车
z6=zeros(1,d６(eｎｄ)－d6（1）＋1）;%(或ｚｅｒos(1,d(enｄ）—ｄ(１)））
ｙ61=ｄ６(1):１:d６(ｅnd)；
for k＝1：length（x6）—1
for　i＝１:ｄ6（ｅnd）-d6(1）＋1
ｉf　d6(k+1)—y61（ｉ）＋m6(k+1）〉=5 &&　y6１（i)—d6（k）-n6(k)>＝５;
　 　z6(i)=y６1(ｉ）;
　 ｅnｄ
　　 　end
end
%将满足条件得列车再进行剔除,求数量最大且满足条件得列车
r＝0;
g＝0；
k=1;
p=１；
l6=zeros（1,leｎｇth(d６）+5０)；
t6=zeｒｏs（1,d6（eｎd）—d6(1)＋1);
fｏr　 i=１:d６(eｎd）—ｄ6(1）
ｉf z6(ｉ)+1==z6（i+1)
k=k+１;
eｎd　
if z6(ｉ)>z6(ｉ+1)
　 r＝ｆloor((k-１）／5)+1;
　 for u=1:r
　　 g=i—k+1＋5＊（u-１）;
　 　ｉf g〉d6（end)—d６(１)+1
　 g=d６（eｎd）—d6(1）+1;
　　 ｅnd
　　　　 l6(g)＝z６（g);
　　　 eｎｄ
　 ｋ=1;
end
end
l6(fｉnd(l6＝=0))=［]；％找到l中０得位置，并令其为空,即删除
w6=ｌ6＋v6;
y６3＝126＊oneｓ（1,lenｇｔh（l6))；
y64=143＊onｅs(1，lｅnｇth（w6）);
fｏr i＝１:lengｔh(ｌ6)
　　 plot([ｌ6(i),w6(ｉ）］,[y63(i）,y６4(ｉ)])；
hold oｎ
end
y61=1２6＊ｏnｅs（１,length（d6)）；
y62＝14３*oneｓ(1,leｎｇth(ｈ６));
foｒ　i=1:lenｇｔh（d6)
plot([d６（ｉ),h6(i)]，[y61(i),y62(i）]);
hold oｎ
ｅnd
％第七个
%设6：42时间为0时刻,转换如下．
a７=[9 11 １2 13 14 １４　15 17 18 19　2１］;
ｂ7=[17 5 1６ 42 6　２1 13　12 2２ ３６ ４0]；
c7=7＊ｏneｓ（1，leｎｇth(a7));
d7=60＊(a７-c７)+b７+１8;
ｅ7=[９ 1１ １2 1３ 14 14　15　17　１8 19 21]；
f7=［32　20　３1 ５7 2１ ３６ 27 ２7　37　5１ 5５］;
g=7*oneｓ(１,lengtｈ(a7))；
h７=60*(e７-g)+f7+18;
%输入经历时间，求众数,得到一个差值矩阵ｆ。
ｘ7=[15
15
15
1５
15
15
15
15
15
１5
15]';
ｖ7=mｏde(x7)
j7=x7—v7;
ｄiｓp(j7)
％矩阵剔除，将ｆ中大于0得用0代替得以矩阵m,将ｆ小于0得数代替得一矩阵ｎ，ｎ，m都就是1＊2０得矩阵。
n7=(j7+abs（j７）)/2；
ｍ７=ｊ７；
m7(fｉnｄ（j7>０))＝0;
%选取可行列车
z7＝ｚeｒoｓ（1,d７(end)-d7(1)+1);%（或zeros（1,ｄ(enｄ）—d（1))）
y7１=d７（１):1:ｄ7(eｎｄ);
fｏr　k=1：ｌeｎｇth(ｘ7）—１
ｆor　i=1:ｄ７(end)-d7(1)＋1
　iｆ ｄ7(k+1)—y７1(i）＋m7（k+１)>=5 ＆＆ y71(i)-d7(k）—n７(ｋ)＞＝5;
　 ｚ7(i）=y71（i)；
　　eｎd
　 enｄ
end
％将满足条件得列车再进行剔除，求数量最大且满足条件得列车
r=0；
g＝０；
ｋ=1；
p=1;
ｌ７=zerｏs(1,leｎgth（d７)+５0）;
t７＝zerｏs（1,d7（ｅｎd）－ｄ7（１)+1);
fｏr 　i＝1：d７(ｅnd)-ｄ7(1） 　
iｆ　z7(ｉ)＋１==z7(ｉ+１)
　 ｋ＝k+1;
eｎd　
if z7(ｉ)＞ｚ7(i+1)
　r=fｌooｒ（(k—1)/5)+１;
　 for ｕ＝1：r
　 g=i—k+1+5*（u-1）;
　 iｆ ｇ〉d7(enｄ)-ｄ7（1)+1
g=d7（end)-ｄ7（1）+1;
　 eｎd
ｌ7（ｇ)＝z7(g）;
ｅｎd
　 ｋ＝1；
enｄ
enｄ
ｌ7(find(l7==０))＝[］；％找到l中0得位置,并令其为空,即删除
w7=l７+v7;
ｙ73＝143*ｏnes(１,lenｇth(l７））;
y74=158＊onｅs(１,length(ｗ７）);
ｆor　i=1:leｎｇth(l7）
　 ｐlｏt（[ｌ7(i）,ｗ7（i）],［ｙ７3(ｉ）,y74(i）]）；
hold　on
ｅnd
y７１＝143＊onｅs(1，lｅngth(d7))；
y72＝１5８＊ones(１，lenｇｔh（h7）);
ｆｏr i=１：ｌength(ｄ7)
plot（[ｄ７(i),h7(ｉ)],[ｙ7１(i),ｙ7２（i)])；
hold on
end
％第八个
％设6:4２时间为0时刻,转换如下．
ａ８=［9 10 11 1２　１4 14 1８ 1９ 19 20 21];
b8=[16 49 6　17 23 38 5０　40 53 52　2１];
ｃ８＝7*ｏnｅｓ(1,lenｇth(a８));
d8=６0*(a８—c8)＋b8＋１8；
e8=［９ 11 11 12 14 14 1９　19 20 2１ 21］;
ｆ8=[３2 5 22 3３ 39 54 ６ 56　9 8 3７];
ｇ=7*oｎeｓ(1，length（a8));
ｈ8＝60＊(e8—g)+f8+18;
％输入经历时间，求众数,得到一个差值矩阵f.
x8=[１６
16
16
16
16
１６
16
１6
16
１6
16］’；
v8＝ｍodｅ（x8)
j8=x８—ｖ８;
disｐ(j8)
%矩阵剔除,将f中大于0得用０代替得以矩阵m,将ｆ小于0得数代替得一矩阵n,n,m都就是１＊２0得矩阵。
n８=(j8+abs(ｊ8)）/2；
m８＝j8;
ｍ8(ｆｉnd（j8〉0))=0;
%选取可行列车
ｚ8=zeｒos（1，d8(enｄ）－ｄ8(1)+１）;%(或zｅros(1，ｄ(end)-d(1))）
y81=ｄ8(1)：1:d8(eｎd）;
ｆoｒ k=1:ｌeｎgtｈ(x8)—1
fｏr ｉ=1：ｄ8(eｎd)—d8（１)+1
if d8（k+１）—ｙ8１（i)+m8(k+1）>=5　&&　y8１(i)—d８(k）—n8（k)〉＝5;
　 z8(i）=y８1(i);
　 　　end
　 　ｅnd
ｅｎｄ
%将满足条件得列车再进行剔除，求数量最大且满足条件得列车
r=0;
g=0;
k=1;
p=１；
l８＝zｅros(1,length(d8）＋50);
t８=zｅｒos(1，d８(end)-d８(1)+1)；
for i＝1:d８（eｎｄ)-d8（1) 　
ｉf　ｚ8(i)+1==ｚ８(i＋１）
　　　 k=k+1;
eｎd
ｉｆ ｚ8(i)>z８(ｉ+1）
　 r＝floor(（k—1)/５）+1;
　 　ｆoｒ u=1：ｒ
　 　 　　g＝i—k+1+5＊(ｕ-1);
　 　 　if g〉d8(enｄ)－d８（１）＋1
　 　 　 g=d8(end)—d8（1）＋１;
　 　 eｎd
　 　l8(ｇ）＝z8（ｇ）;
　end
　　 k=1;
end
end
l8(find(l8＝=0）)=[]；%找到ｌ中０得位置，并令其为空,即删除
w８=l8+v8；
y83=１58*ones（1，length（ｌ8)）;
y8４=1７４＊onｅs(1，lｅngtｈ（ｗ8));
foｒ i＝1:ｌength(ｌ8）
　 plot(［l８(ｉ),w8(i)］，[y８3（i)，y84(i)］）；
hold on
enｄ
y81=１58＊ones（1,lｅnｇth（d8));
y８2＝1７4*ones(1，ｌeｎgth(h８）)；
foｒ i=１：ｌengｔｈ（d8)
ｐｌot（［ｄ8（i)，ｈ８(i）]，[y81（i),ｙ８2(i）］);
ｈold on
end
%第九个
％设6:４２时间为0时刻，转换如下.
a9=［９ 10 11　11 １１ 12　１２　14 １5　１6 17 １７ 18 1８ 19　１9 １9 １９ 20 ２0 2０ 20 21　21 21 22];
b9=[３４　４９　17 24 52 30 35 4１ 27 48　28 59 １７ ５8 ３ ８ 21 ５8 １１ 17 27 43 1０　２１ ３9 １９]；
c9=7*oｎes(1，leｎgｔｈ（a9))；
d９=60*(a9-c9）+b９＋１8；
e9=[９　11 11 12 12 12 １2　1５ 1５ １7　17 18 18　1９ １9　１９ １9 20 20 20 20 21 21　21 22　2３］;
f9=[55 11 ３8　７ 15 51　56 2 48 9 ４9　2０ ３8 19 2４ ２９ ４3　１９ 3２　3８　48　５　31 42 ０　1];
g=7*ones（1，length(ａ９)）;
h9＝６０＊（e９-g)+f9+１8；
%输入经历时间，求众数，得到一个差值矩阵f。
x9＝[２1
22
21
43
２3
２１
２1
21
2１
21
２1
21
2１
21
21
２1
22
２1
２１
２1
21
22
21
21
21
42]’;
ｖ9=mｏde(x９)
ｊ９＝x9—v9；
ｄｉｓp（j9)
%矩阵剔除,将f中大于0得用0代替得以矩阵m，将f小于0得数代替得一矩阵n，n,m都就是1*２０得矩阵。
n9=（j9+ａｂs(j９))/2；
m9＝j9；
m9(fｉnd(j9〉0））=0;
％选取可行列车
z9=zeros（1,d９(eｎｄ）-d9(1)＋1)；％(或ｚeroｓ(１,d（ｅnd)－d(1))）
y９１=d9(１）:1:d9(enｄ);
foｒ k=1:ｌengtｈ(ｘ9）-1
for　i＝1:d９（end）—d9（1）＋1
iｆ d9（k+１)—ｙ91(i)+m9（k+1)〉＝5　＆& y９1(i)—ｄ9(ｋ）-n9（ｋ)＞＝5；
　ｚ9（ｉ)=y91（i）；
　　 　 ｅnd
　 　　ｅｎｄ
eｎd
％将满足条件得列车再进行剔除,求数量最大且满足条件得列车
r=0；
g=０;
ｋ=1；
p=１;
ｌ9=zeroｓ(1,lengtｈ(d９)＋50);
t9=zeroｓ（１,d9（ｅnd)-d９（1)+1)；
for i＝1:d９(eｎd）—d9（1） 　
ｉｆ z９(ｉ）＋１==z9(i+1)
　 k＝k+1;
enｄ　
iｆ ｚ９（ｉ)〉z9(i+1)
ｒ=floor((k—1)/５）+1；
　for　ｕ=1:ｒ
　　 g=i-ｋ+1+5*(u-1）；
　 　 　if g〉d9(enｄ)－d9(1）+１
　 　　 ｇ=d9(end)—d9(1）+1;
　 　 ｅnｄ
　　 l9(g）=ｚ9（g)；
　 end
k＝1；
end
enｄ
ｌ9(find(l9=＝0))＝[］;％找到l中0得位置，并令其为空，即删除
w9＝l９＋v9;
ｙ93=174＊ｏnｅs（１,ｌｅｎgｔｈ（l9));
y9４＝195*ones（1,ｌengｔｈ(w9)）;
for ｉ=1:lengｔh(l9)
　　 plｏt(［l9（ｉ),w９(ｉ）］,［y9３（ｉ)，y94(i)］）;
holｄ oｎ
eｎd
ｙ9１=174*oｎes（1,lenｇth(d9）);
y９2=19５*ones(１,leｎgtｈ（h9））;
for i=１:ｌeｎgtｈ(d9)
ｐloｔ(［d9(i）,h9(ｉ）］,[y91(ｉ),y９2(ｉ）］);
hoｌd on
eｎd
ｇtext(’北京西’）；
gtexｔ（’保定东’)；
gtext('定州东’);
ｇtexｔ(＇石家庄北＇）;
ｇteｘt（'邢台东’）;
gteｘｔ('邯郸东’）;
gteｘt（'安阳东’)；
gtext（’鹤壁东'）；
gtｅxt(’新乡东＇)；
gtext('郑州东＇）；