《列车运行图MATLAB程序》由会员分享，可在线阅读，更多相关《列车运行图MATLAB程序（19页珍藏版）》请在金锄头文库上搜索。
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(x1)j1
3、=x1-v1;disp(j1)%矩阵剔除，将f中大于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵。n1=(j1+abs(j1)/2;m1=j1;m1(find(j10)=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)=y11(i);endendend%将满足条件的列
4、车再进行剔除，求数量最大且满足条件的列车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);ifgd1(end)-d1(1)+1g=d1(end)-d1(1)+1;endl1(g)=z1(g);endk=1;endendl1(find(l1=0)=;%找到l中0的位置，并令其为空，即删除w1=l1+v1;y
5、13=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=644580245316322358281919;c2=7*one
6、s(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的矩阵。n2=(j2+abs(j2)/2;m2=j2;m2(find(j20)=0
7、;%选取可行列车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(1,length(d2)+50);t2=zeros(1,d2(end)-d2(1)+1);fori=1:d2(end
8、)-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);ifgd2(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,length(w2);fori=1:length(l2)plot(l2(i),w2(i),y23(i),y24(i);holdon
9、endy21=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(a3);d3=60*(a3-c3)+b3+18;e3=889910111111131415161617171819202021;f3
10、=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=(j3+abs(j3)/2;m3=j3;m3(find(j30)=0;%选取可行列车z3=zeros(1,d3(end)-d3(1)+1);%
11、(或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,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;endifz3
12、(i)z3(i+1)r=floor(k-1)/5)+1;foru=1:rg=i-k+1+5*(u-1);ifgd3(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,length(w3);fori=1:length(l3)plot(l3(i),w3(i),y33(i),y34(i);holdonendy31=59*ones(1,length(d3);y32=82*on
13、es(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=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(j4)%矩阵剔除，将f中大于0的用0代替得以矩阵m，将f小于0的数代替得一矩阵n,n,m都是1*20的矩阵。n4=(j4+abs(j4)/2;m4=j4;m4(find(j40)=0;%选取可行列车
