列车运行图MATLAB程序

1、%第一个%设6:42时间为0时刻，转换如下。a1=6777777889999910101010101011111212121213131414151515161616161717171718181818192021;b1=420102233405215301217222740510152845503343713344351825436404501835573193237212364222476;c1=7*ones(1,length(a1);d1=60*(a1-c1)+b1+18;e1=7778888899910101010101111111112121212131313141515151616

2、161617171818181818191919202121;f1=23415831528375611535810152146513926311424485415244615624472126485923380713184301730173547;g=7*ones(1,length(a1);h1=60*(e1-g)+f1+18;%输入经历时间，求众数，得到一个差值矩阵f。x1=414148414248454141414148484141414841414141414141414141574141414141484148415748414141484148554841'v1=mode(x

3、1)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)-1fori=1:d1(end)-d1(1)+1ifd1(k+1)-y11(i)+m1(k+1)>=5&&y11(i)-d1(k)-n1(k)>=5;z1(i)

4、=y11(i);endendend%将满足条件的列车再进行剔除，求数量最大且满足条件的列车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;endifz1(i)>z1(i+1)r=floor(k-1)/5)+1;foru=1:rg=i-k+1+5*(u-1);ifg>d1(end)-d1(1)+1g=d1(end)-d1(1)+1;endl1(g)=z1(g);endk=1;endendl1(find(l1=

5、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);holdonendy11=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);holdonend%第二个%设6:42时间为0时刻，转换如下。a2=88810101011131415161920;b2=

6、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=18201818181818182018181818'v2=mode(x2)j2=x2-v2;disp(j2)%矩阵剔除，将f中大于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵。

7、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)-1fori=1:d2(end)-d2(1)+1ifd2(k+1)-y21(i)+m2(k+1)>=5&&y21(i)-d2(k)-n2(k)>=5;z2(i)=y21(i);endendend%将满足条件的列车再进行剔除，求数量最大且满足条件的列车r=0;g=0;k=1;p=1;l2=zeros(

8、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;endifz2(i)>z2(i+1)r=floor(k-1)/5)+1;foru=1:rg=i-k+1+5*(u-1);ifg>d2(end)-d2(1)+1g=d2(end)-d2(1)+1;endl2(g)=z2(g);endk=1;endendl2(find(l2=0)=;%找到l中0的位置，并令其为空，即删除w2=l2+v2;y23=41*ones(1,length(l2);y24=59*ones(1

9、,length(w2);fori=1:length(l2)plot(l2(i),w2(i),y23(i),y24(i);holdonendy21=41*ones(1,length(d2);y22=59*ones(1,length(h1);fori=1:length(d2)plot(d2(i),h1(i),y21(i),y22(i);holdonend%第三个%设6:42时间为0时刻，转换如下。a3=789910101111121314151616171819191920;b3=502661820501336565252512848261010395839;c3=7*ones(1,length(

10、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=2323252323222323232324302322232323232323'v3=mode(x3)j3=x3-v3;disp(j3)%矩阵剔除，将f中大于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵。n3=(

11、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)-1fori=1:d3(end)-d3(1)+1ifd3(k+1)-y31(i)+m3(k+1)>=5&&y31(i)-d3(k)-n3(k)>=5;z3(i)=y31(i);endendend%将满足条件的列车再进行剔除，求数量最大且满足条件的列车r=0;g=0;k=1;p=1;l3=zeros(1,le

12、ngth(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;endifz3(i)>z3(i+1)r=floor(k-1)/5)+1;foru=1:rg=i-k+1+5*(u-1);ifg>d3(end)-d3(1)+1g=d3(end)-d3(1)+1;endl3(g)=z3(g);endk=1;endendl3(find(l3=0)=;%找到l中0的位置，并令其为空，即删除w3=l3+v3;y33=59*ones(1,length(l3);y34=82*ones(1,len

13、gth(w3);fori=1:length(l3)plot(l3(i),w3(i),y33(i),y34(i);holdonendy31=59*ones(1,length(d3);y32=82*ones(1,length(h3);fori=1:length(d3)plot(d3(i),h3(i),y31(i),y32(i);holdonend%第四个%设6:42时间为0时刻，转换如下。a4=78888899910101111111313131315151616161818192020;b4=0021284857915255159163957213242567475245443642939;c4

14、=7*ones(1,length(a4);d4=60*(a4-c4)+b4+18;e4=788899991011111112121314141415161616171819202021;f4=2727495616323742019274415254901024351533592931410377;g=7*ones(1,length(a4);h4=60*(e4-g)+f4+18;%输入经历时间，求众数，得到一个差值矩阵f。x4=27272828283528273528282836282828282828282835352728282828'v4=mode(x4)j4=x4-v4;disp

15、(j4)%矩阵剔除，将f中大于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵。n4=(j4+abs(j4)/2;m4=j4;m4(find(j4>0)=0;%选取可行列车z4=zeros(1,d4(end)-d4(1)+1);%(或zeros(1，d(end)-d(1)y41=d4(1):1:d4(end);fork=1:length(x4)-1fori=1:d4(end)-d4(1)+1ifd4(k+1)-y41(i)+m4(k+1)>=5&&y41(i)-d4(k)-n4(k)>=5;z4(i)=y41(i);endende

16、nd%将满足条件的列车再进行剔除，求数量最大且满足条件的列车r=0;g=0;k=1;p=1;l4=zeros(1,length(d4)+50);t4=zeros(1,d4(end)-d4(1)+1);fori=1:d4(end)-d4(1)ifz4(i)+1=z4(i+1)k=k+1;endifz4(i)>z4(i+1)r=floor(k-1)/5)+1;foru=1:rg=i-k+1+5*(u-1);ifg>d4(end)-d4(1)+1g=d4(end)-d4(1)+1;endl4(g)=z4(g);endk=1;endendl4(find(l4=0)=;%找到l中0的位置，并

17、令其为空，即删除w4=l4+v4;y43=82*ones(1,length(l4);y44=110*ones(1,length(w4);fori=1:length(l4)plot(l4(i),w4(i),y43(i),y44(i);holdonendy41=82*ones(1,length(d4);y42=110*ones(1,length(h4);fori=1:length(d4)plot(d4(i),h4(i),y41(i),y42(i);holdonend%第五个%设6:42时间为0时刻，转换如下。a5=78899910111214141516161718;b5=292951193944

18、2211819343917354533;c5=7*ones(1,length(a5);d5=60*(a5-c5)+b5+18;e5=789991010111214141516161818;f5=45457355501837343550553351149;g=7*ones(1,length(a5);h5=60*(e5-g)+f5+18;%输入经历时间，求众数，得到一个差值矩阵f。x5=16161616161616161616161616161616'v5=mode(x5)j5=x5-v5;disp(j5)%矩阵剔除，将f中大于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都

19、是1*20的矩阵。n5=(j5+abs(j5)/2;m5=j5;m5(find(j5>0)=0;%选取可行列车z5=zeros(1,d5(end)-d5(1)+1);%(或zeros(1，d(end)-d(1)y51=d5(1):1:d5(end);fork=1:length(x5)-1fori=1:d5(end)-d5(1)+1ifd5(k+1)-y51(i)+m5(k+1)>=5&&y51(i)-d5(k)-n5(k)>=5;z5(i)=y51(i);endendend%将满足条件的列车再进行剔除，求数量最大且满足条件的列车r=0;g=0;k=1;p=1;

20、l5=zeros(1,length(d5)+50);t5=zeros(1,d5(end)-d5(1)+1);fori=1:d5(end)-d5(1)ifz5(i)+1=z5(i+1)k=k+1;endifz5(i)>z5(i+1)r=floor(k-1)/5)+1;foru=1:rg=i-k+1+5*(u-1);ifg>d5(end)-d5(1)+1g=d5(end)-d5(1)+1;endl5(g)=z5(g);endk=1;endendl5(find(l5=0)=;%找到l中0的位置，并令其为空，即删除w5=l5+v5;y53=110*ones(1,length(l5);y54

21、=126*ones(1,length(w5);fori=1:length(l5)plot(l5(i),w5(i),y53(i),y54(i);holdonendy51=110*ones(1,length(d5);y52=126*ones(1,length(h5);fori=1:length(d5)plot(d5(i),h5(i),y51(i),y52(i);holdonend%第六个%设6:42时间为0时刻，转换如下。a6=891010111313141515161617181818191919;b6=58162463234753165841533031134304056;c6=7*ones(

22、1,length(a6);d6=60*(a6-c6)+b6+18;e6=991011111314151516161717181818191920;f6=153320320404103315581047202851475713;g=7*ones(1,length(a6);h6=60*(e6-g)+f6+18;%输入经历时间，求众数，得到一个差值矩阵f。x6=17171817171717171717171717171717171717'v6=mode(x6)j6=x6-v6;disp(j6)%矩阵剔除，将f中大于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵

23、。n6=(j6+abs(j6)/2;m6=j6;m6(find(j6>0)=0;%选取可行列车z6=zeros(1,d6(end)-d6(1)+1);%(或zeros(1，d(end)-d(1)y61=d6(1):1:d6(end);fork=1:length(x6)-1fori=1:d6(end)-d6(1)+1ifd6(k+1)-y61(i)+m6(k+1)>=5&&y61(i)-d6(k)-n6(k)>=5;z6(i)=y61(i);endendend%将满足条件的列车再进行剔除，求数量最大且满足条件的列车r=0;g=0;k=1;p=1;l6=zeros

24、(1,length(d6)+50);t6=zeros(1,d6(end)-d6(1)+1);fori=1:d6(end)-d6(1)ifz6(i)+1=z6(i+1)k=k+1;endifz6(i)>z6(i+1)r=floor(k-1)/5)+1;foru=1:rg=i-k+1+5*(u-1);ifg>d6(end)-d6(1)+1g=d6(end)-d6(1)+1;endl6(g)=z6(g);endk=1;endendl6(find(l6=0)=;%找到l中0的位置，并令其为空，即删除w6=l6+v6;y63=126*ones(1,length(l6);y64=143*one

25、s(1,length(w6);fori=1:length(l6)plot(l6(i),w6(i),y63(i),y64(i);holdonendy61=126*ones(1,length(d6);y62=143*ones(1,length(h6);fori=1:length(d6)plot(d6(i),h6(i),y61(i),y62(i);holdonend%第七个%设6:42时间为0时刻，转换如下。a7=911121314141517181921;b7=17516426211312223640;c7=7*ones(1,length(a7);d7=60*(a7-c7)+b7+18;e7=91

26、1121314141517181921;f7=3220315721362727375155;g=7*ones(1,length(a7);h7=60*(e7-g)+f7+18;%输入经历时间，求众数，得到一个差值矩阵f。x7=1515151515151515151515'v7=mode(x7)j7=x7-v7;disp(j7)%矩阵剔除，将f中大于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵。n7=(j7+abs(j7)/2;m7=j7;m7(find(j7>0)=0;%选取可行列车z7=zeros(1,d7(end)-d7(1)+1);%(或ze

27、ros(1，d(end)-d(1)y71=d7(1):1:d7(end);fork=1:length(x7)-1fori=1:d7(end)-d7(1)+1ifd7(k+1)-y71(i)+m7(k+1)>=5&&y71(i)-d7(k)-n7(k)>=5;z7(i)=y71(i);endendend%将满足条件的列车再进行剔除，求数量最大且满足条件的列车r=0;g=0;k=1;p=1;l7=zeros(1,length(d7)+50);t7=zeros(1,d7(end)-d7(1)+1);fori=1:d7(end)-d7(1)ifz7(i)+1=z7(i+1)

28、k=k+1;endifz7(i)>z7(i+1)r=floor(k-1)/5)+1;foru=1:rg=i-k+1+5*(u-1);ifg>d7(end)-d7(1)+1g=d7(end)-d7(1)+1;endl7(g)=z7(g);endk=1;endendl7(find(l7=0)=;%找到l中0的位置，并令其为空，即删除w7=l7+v7;y73=143*ones(1,length(l7);y74=158*ones(1,length(w7);fori=1:length(l7)plot(l7(i),w7(i),y73(i),y74(i);holdonendy71=143*one

29、s(1,length(d7);y72=158*ones(1,length(h7);fori=1:length(d7)plot(d7(i),h7(i),y71(i),y72(i);holdonend%第八个%设6:42时间为0时刻，转换如下。a8=910111214141819192021;b8=164961723385040535221;c8=7*ones(1,length(a8);d8=60*(a8-c8)+b8+18;e8=911111214141919202121;f8=325223339546569837;g=7*ones(1,length(a8);h8=60*(e8-g)+f8+18

30、;%输入经历时间，求众数，得到一个差值矩阵f。x8=1616161616161616161616'v8=mode(x8)j8=x8-v8;disp(j8)%矩阵剔除，将f中大于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵。n8=(j8+abs(j8)/2;m8=j8;m8(find(j8>0)=0;%选取可行列车z8=zeros(1,d8(end)-d8(1)+1);%(或zeros(1，d(end)-d(1)y81=d8(1):1:d8(end);fork=1:length(x8)-1fori=1:d8(end)-d8(1)+1ifd8(k+1

31、)-y81(i)+m8(k+1)>=5&&y81(i)-d8(k)-n8(k)>=5;z8(i)=y81(i);endendend%将满足条件的列车再进行剔除，求数量最大且满足条件的列车r=0;g=0;k=1;p=1;l8=zeros(1,length(d8)+50);t8=zeros(1,d8(end)-d8(1)+1);fori=1:d8(end)-d8(1)ifz8(i)+1=z8(i+1)k=k+1;endifz8(i)>z8(i+1)r=floor(k-1)/5)+1;foru=1:rg=i-k+1+5*(u-1);ifg>d8(end)-d8

32、(1)+1g=d8(end)-d8(1)+1;endl8(g)=z8(g);endk=1;endendl8(find(l8=0)=;%找到l中0的位置，并令其为空，即删除w8=l8+v8;y83=158*ones(1,length(l8);y84=174*ones(1,length(w8);fori=1:length(l8)plot(l8(i),w8(i),y83(i),y84(i);holdonendy81=158*ones(1,length(d8);y82=174*ones(1,length(h8);fori=1:length(d8)plot(d8(i),h8(i),y81(i),y82(

33、i);holdonend%第九个%设6:42时间为0时刻，转换如下。a9=910111111121214151617171818191919192020202021212122;b9=34491724523035412748285917583821581117274310213919;c9=7*ones(1,length(a9);d9=60*(a9-c9)+b9+18;e9=911111212121215151717181819191919202020202121212223;f9=5511387155156248949203819242943193238485314201;g=7*ones(1,length(a9);h9=60*(e9-g)+f9+18;%输入经历时间，求众数，得到一个差值矩阵f。x9=2122214323212121212121212