列车运行图MATLAB程序.docx
《列车运行图MATLAB程序.docx》由会员分享，可在线阅读，更多相关《列车运行图MATLAB程序.docx（27页珍藏版）》请在冰豆网上搜索。
列车运行图MATLAB程序
%第一个
%设6:
42时间为0时刻，转换如下。
a1=[6777777889999910101010101011111212121213131414151515161616161717171718181818192021];
b1=[420102233405215301217222740510152845503343713344351825436404501835573193237212364222476];
c1=7*ones（1,length（a1））;
d1=60*（a1-c1）+b1+18;
e1=[7778888899910101010101111111112121212131313141515151616161617171818181818191919202121];
f1=[23415831528375611535810152146513926311424485415244615624472126485923380713184301730173547];
g=7*ones（1,length（a1））;
h1=60*（e1-g）+f1+18;
%输入经历时间，求众数，得到一个差值矩阵f。
x1=[41
41
48
41
42
48
45
41
41
41
41
48
48
41
41
41
48
41
41
41
41
41
41
41
41
41
41
57
41
41
41
41
41
48
41
48
41
57
48
41
41
41
48
41
48
55
48
41]';
v1=mode（x1）
j1=x1-v1;
disp（j1）
%矩阵剔除，将f于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵。
n1=（j1+abs（j1））/2;
m1=j1;
m1（find（j1>0））=0;
%选取可行列车
z1=zeros（1,d1（end）-d1
（1）+1）;%〔或zeros〔1，d〔end〕-d
（1）〕〕
y11=d1
（1）:
1:
d1（end）;
fork=1:
length（x1）-1
fori=1:
d1（end）-d1
（1）+1
ifd1（k+1）-y11（i）+m1（k+1）>=5&&y11（i）-d1（k）-n1（k）>=5;
z1（i）=y11（i）;
end
end
end
%将满足条件的列车再进展剔除，求数量最大且满足条件的列车
r=0;
g=0;
k=1;
p=1;
l1=zeros（1,length（d1）+50）;
t1=zeros（1,d1（end）-d1
（1）+1）;
fori=1:
d1（end）-d1
（1）
ifz1（i）+1==z1（i+1）
k=k+1;
end
ifz1（i）>z1（i+1）
r=floor（（k-1）/5）+1;
foru=1:
r
g=i-k+1+5*（u-1）;
ifg>d1（end）-d1
（1）+1
g=d1（end）-d1
（1）+1;
end
l1（g）=z1（g）;
end
k=1;
end
end
l1（find（l1==0））=[];%找到l中0的位置，并令其为空，即删除
w1=l1+v1;
y13=0.1*ones（1,length（l1））;
y14=41*ones（1,length（w1））;
fori=1:
length（l1）
plot（[l1（i）,w1（i）],[y13（i）,y14（i）]）;
holdon
end
y11=0.1*ones（1,length（d1））;
y12=41*ones（1,length（h1））;
fori=1:
length（d1）
plot（[d1（i）,h1（i）],[y11（i）,y12（i）]）;
holdon
end
%第二个
%设6:
42时间为0时刻，转换如下。
a2=[88810101011131415161920];
b2=[644580245316322358281919];
c2=7*ones（1,length（a2））;
d2=60*（a2-c2）+b2+18;
e2=[89910101111131416161920];
f2=[2441618421134504316463737];
g=7*ones（1,length（a2））;
h1=60*（e2-g）+f2+18;
%输入经历时间，求众数，得到一个差值矩阵f。
x2=[18
20
18
18
18
18
18
18
20
18
18
18
18
]';
v2=mode（x2）
j2=x2-v2;
disp（j2）
%矩阵剔除，将f于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵。
n2=（j2+abs（j2））/2;
m2=j2;
m2（find（j2>0））=0;
%选取可行列车
z2=zeros（1,d2（end）-d2
（1）+1）;%〔或zeros〔1，d〔end〕-d
（1）〕〕
y21=d2
（1）:
1:
d2（end）;
fork=1:
length（x2）-1
fori=1:
d2（end）-d2
（1）+1
ifd2（k+1）-y21（i）+m2（k+1）>=5&&y21（i）-d2（k）-n2（k）>=5;
z2（i）=y21（i）;
end
end
end
%将满足条件的列车再进展剔除，求数量最大且满足条件的列车
r=0;
g=0;
k=1;
p=1;
l2=zeros（1,length（d2）+50）;
t2=zeros（1,d2（end）-d2
（1）+1）;
fori=1:
d2（end）-d2
（1）
ifz2（i）+1==z2（i+1）
k=k+1;
end
ifz2（i）>z2（i+1）
r=floor（（k-1）/5）+1;
foru=1:
r
g=i-k+1+5*（u-1）;
ifg>d2（end）-d2
（1）+1
g=d2（end）-d2
（1）+1;
end
l2（g）=z2（g）;
end
k=1;
end
end
l2（find（l2==0））=[];%找到l中0的位置，并令其为空，即删除
w2=l2+v2;
y23=41*ones（1,length（l2））;
y24=59*ones（1,length（w2））;
fori=1:
length（l2）
plot（[l2（i）,w2（i）],[y23（i）,y24（i）]）;
holdon
end
y21=41*ones（1,length（d2））;
y22=59*ones（1,length（h1））;
fori=1:
length（d2）
plot（[d2（i）,h1（i）],[y21（i）,y22（i）]）;
holdon
end
%第三个
%设6:
42时间为0时刻，转换如下。
a3=[789910101111121314151616171819191920];
b3=[502661820501336565252512848261010395839];
c3=7*ones（1,length（a3））;
d3=60*（a3-c3）+b3+18;
e3=[889910111111131415161617171819202021];
f3=[13493141431336591915162151114933332212];
g=7*ones（1,length（a3））;
h3=60*（e3-g）+f3+18;
%输入经历时间，求众数，得到一个差值矩阵f。
x3=[23
23
25
23
23
22
23
23
23
23
24
30
23
22
23
23
23
23
23
23]';
v3=mode（x3）
j3=x3-v3;
disp（j3）
%矩阵剔除，将f于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵。
n3=（j3+abs（j3））/2;
m3=j3;
m3（find（j3>0））=0;
%选取可行列车
z3=zeros（1,d3（end）-d3
（1）+1）;%〔或zeros〔1，d〔end〕-d
（1）〕〕
y31=d3
（1）:
1:
d3（end）;
fork=1:
length（x3）-1
fori=1:
d3（end）-d3
（1）+1
ifd3（k+1）-y31（i）+m3（k+1）>=5&&y31（i）-d3（k）-n3（k）>=5;
z3（i）=y31（i）;
end
end
end
%将满足条件的列车再进展剔除，求数量最大且满足条件的列车
r=0;
g=0;
k=1;
p=1;
l3=zeros（1,length（d3）+50）;
t3=zeros（1,d3（end）-d3
（1）+1）;
fori=1:
d3（end）-d3
（1）
ifz3（i）+1==z3（i+1）
k=k+1;
end
ifz3（i）>z3（i+1）
r=floor（（k-1）/5）+1;
foru=1:
r
g=i-k+1+5*（u-1）;
ifg>d3（end）-d3
（1）+1
g=d3（end）-d3
（1）+1;
end
l3（g）=z3（g）;
end
k=1;
end
end
l3（find（l3==0））=[];%找到l中0的位置，并令其为空，即删除
w3=l3+v3;
y33=59*ones（1,length（l3））;
y34=82*ones（1,length（w3））;
fori=1:
length（l3）
plot（[l3（i）,w3（i）],[y33（i）,y34（i）]）;
holdon
end
y31=59*ones（1,length（d3））;
y32=82*ones（1,length（h3））;
fori=1:
length（d3）
plot（[d3（i）,h3（i）],[y31（i）,y32（i）]）;
holdon
end
%第四个
%设6:
42时间为0时刻，转换如下。
a4=[78888899910101111111313131315151616161818192020];
b4=[0021284857915255159163957213242567475245443642939];
c4=7*ones（1,length（a4））;
d4=60*（a4-c4）+b4+18;
e4=[78889
升级会员 