列车运行图MATLAB程序.docx

1、列车运行图MATLAB程序%第一个%设6:42时间为0时刻，转换如下。a1=6 7 7 7 7 7 7 8 8 9 9 9 9 9 10 10 10 10 10 10 11 11 12 12 12 12 13 13 14 14 15 15 15 16 16 16 16 17 17 17 17 18 18 18 18 19 20 21;b1=42 0 10 22 33 40 52 15 30 12 17 22 27 40 5 10 15 28 45 50 33 43 7 13 34 43 5 18 25 43 6 40 45 0 18 35 57 3 19 32 37 2 12 36 42 22

2、 47 6;c1=7*ones(1,length(a1);d1=60*(a1-c1)+b1+18;e1=7 7 7 8 8 8 8 8 9 9 9 10 10 10 10 10 11 11 11 11 12 12 12 12 13 13 13 14 15 15 15 16 16 16 16 17 17 18 18 18 18 18 19 19 19 20 21 21;f1=23 41 58 3 15 28 37 56 11 53 58 10 15 21 46 51 3 9 26 31 14 24 48 54 15 24 46 15 6 24 47 21 26 48 59 23 38 0 7 1

3、3 18 43 0 17 30 17 35 47;g=7*ones(1,length(a1);h1=60*(e1-g)+f1+18;%输入经历时间，求众数，得到一个差值矩阵f。x1=414148414248454141414148484141414841414141414141414141574141414141484148415748414141484148554841;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(j10)

4、=0;%选取可行列车z1=zeros(1,d1(end)-d1(1)+1);%或zeros1，dend-d(1)y11=d1(1):1:d1(end);for k=1:length(x1)-1for i=1:d1(end)-d1(1)+1if d1(k+1)-y11(i)+m1(k+1)=5 & y11(i)-d1(k)-n1(k)=5; z1(i)=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);for i=1:d

5、1(end)-d1(1) if z1(i)+1=z1(i+1) k=k+1;endif z1(i)z1(i+1) r=floor(k-1)/5)+1;for u=1:r g=i-k+1+5*(u-1);if gd1(end)-d1(1)+1 g=d1(end)-d1(1)+1;end l1(g)=z1(g);end k=1;endendl1(find(l1=0)=;%找到l中0的位置，并令其为空，即删除w1=l1+v1;y13=0.1*ones(1,length(l1);y14=41*ones(1,length(w1);for i=1:length(l1) plot(l1(i),w1(i),y

6、13(i),y14(i);hold onendy11=0.1*ones(1,length(d1);y12=41*ones(1,length(h1);for i=1:length(d1)plot(d1(i),h1(i),y11(i),y12(i);hold onend%第二个%设6:42时间为0时刻，转换如下。a2=8 8 8 10 10 10 11 13 14 15 16 19 20;b2=6 44 58 0 24 53 16 32 23 58 28 19 19;c2=7*ones(1,length(a2);d2=60*(a2-c2)+b2+18;e2=8 9 9 10 10 11 11 13

7、 14 16 16 19 20;f2=24 4 16 18 42 11 34 50 43 16 46 37 37;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;%选取可行列车z2=zeros(1,d2(end)-d2(1)+1);

8、%或zeros1，dend-d(1)y21=d2(1):1:d2(end);for k=1:length(x2)-1for i=1:d2(end)-d2(1)+1if d2(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);for i=1:d2(end)-d2(1) if z2(i)+1=z2(i+1) k=k+1;

9、endif z2(i)z2(i+1) r=floor(k-1)/5)+1;for u=1:r g=i-k+1+5*(u-1);if gd2(end)-d2(1)+1 g=d2(end)-d2(1)+1;end l2(g)=z2(g);end k=1;endendl2(find(l2=0)=;%找到l中0的位置，并令其为空，即删除w2=l2+v2;y23=41*ones(1,length(l2);y24=59*ones(1,length(w2);for i=1:length(l2) plot(l2(i),w2(i),y23(i),y24(i);hold onendy21=41*ones(1,le

10、ngth(d2);y22=59*ones(1,length(h1);for i=1:length(d2)plot(d2(i),h1(i),y21(i),y22(i);hold onend%第三个%设6:42时间为0时刻，转换如下。a3=7 8 9 9 10 10 11 11 12 13 14 15 16 16 17 18 19 19 19 20;b3=50 26 6 18 20 50 13 36 56 52 52 51 28 48 26 10 10 39 58 39;c3=7*ones(1,length(a3);d3=60*(a3-c3)+b3+18;e3=8 8 9 9 10 11 11 1

11、1 13 14 15 16 16 17 17 18 19 20 20 21;f3=13 49 31 41 43 13 36 59 19 15 16 21 51 11 49 33 33 2 21 2;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;

12、m3=j3;m3(find(j30)=0;%选取可行列车z3=zeros(1,d3(end)-d3(1)+1);%或zeros1，dend-d(1)y31=d3(1):1:d3(end);for k=1:length(x3)-1for i=1:d3(end)-d3(1)+1if d3(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)-

13、d3(1)+1);for i=1:d3(end)-d3(1) if z3(i)+1=z3(i+1) k=k+1;endif z3(i)z3(i+1) r=floor(k-1)/5)+1;for u=1:r g=i-k+1+5*(u-1);if gd3(end)-d3(1)+1 g=d3(end)-d3(1)+1;end l3(g)=z3(g);end k=1;endendl3(find(l3=0)=;%找到l中0的位置，并令其为空，即删除w3=l3+v3;y33=59*ones(1,length(l3);y34=82*ones(1,length(w3);for i=1:length(l3) p

14、lot(l3(i),w3(i),y33(i),y34(i);hold onendy31=59*ones(1,length(d3);y32=82*ones(1,length(h3);for i=1:length(d3)plot(d3(i),h3(i),y31(i),y32(i);hold onend%第四个%设6:42时间为0时刻，转换如下。a4=7 8 8 8 8 8 9 9 9 10 10 11 11 11 13 13 13 13 15 15 16 16 16 18 18 19 20 20;b4=0 0 21 28 48 57 9 15 25 51 59 16 39 57 21 32 42 56 7 47 5 24 54 4 36 42 9 39;c4=7*ones(1,length(a4);d4=60*(a4-c4)+b4+18;e4=7 8 8 8 9