【数学建模集训系列】公交查询系统的matlab实现-只含公交的查询

clear
startP=input('请输入起点：');
endP=input('请输入终点：');
tic
load SSdata_con;
%load SSdata_tim;
%load SSdata_L;
%load SSdata_mon;
load SLdata_sl;
disp('------------------------------直达--------------------------------------')
direct=0;
result_mat=[];
%-----------
for i=1:520
if SLmat1(startP,i,1)>0&&SLmat1(endP,i,1)>0&&(SLmat1(endP,i,1)-SLmat1(startP,i,1)>0)
distance=SLmat1(endP,i,1)-SLmat1(startP,i,1);
money=calc_money(i,distance);
time=distance*3;
%direct=strvcat(direct,num2str(i));
direct=1;
result=[i distance 0 money time];
result_mat=[result_mat(:,:);result(1,:)];
%disp(['可通过L' num2str(i) '经过' num2str(distance) '个站点上行直达----价格：' num2str(money) '----时间：' num2str(time)])
end
if SLmat1(startP,i,2)>0&&SLmat1(endP,i,2)>0&&(SLmat1(endP,i,2)-SLmat1(startP,i,2)>0)
distance=SLmat1(endP,i,2)-SLmat1(startP,i,2);
money=calc_money(i,distance);
time=distance*3;
%direct=strvcat(direct,num2str(i));
direct=1;
result=[i distance 1 money time];
result_mat=[result_mat(:,:);result(1,:)];
%disp(['可通过L' num2str(i) '经过' num2str(distance) '个站点下行直达----价格：' num2str(money) '----时间：' num2str(time)])
end
end
if size(result_mat,1)~=0
money_sort=sortrows(result_mat,4);%按钱排序
time_sort=sortrows(result_mat,5);%按时间排序
disp('按价格排序的最佳路线')
if size(money_sort,1)<3 hang=size(money_sort,1);
else hang=3;
end
for i=1:hang
if money_sort(i,3)==0
disp(['可通过L' num2str(money_sort(i,1)) '经过' num2str(money_sort(i,2)) '个站点上行直达----价格：' num2str(money_sort(i,4)) '----时间：' num2str(money_sort(i,5))])
elseif money_sort(i,3)==1
disp(['可通过L' num2str(money_sort(i,1)) '经过' num2str(money_sort(i,2)) '个站点下行直达----价格：' num2str(money_sort(i,4)) '----时间：' num2str(money_sort(i,5))])
end
end
disp('按时间排序的最佳路线')
for i=1:hang
if time_sort(i,3)==0
disp(['可通过L' num2str(time_sort(i,1)) '经过' num2str(time_sort(i,2)) '个站点上行直达----价格：' num2str(time_sort(i,4)) '----时间：' num2str(time_sort(i,5))])
elseif time_sort(i,3)==1
disp(['可通过L' num2str(time_sort(i,1)) '经过' num2str(time_sort(i,2)) '个站点下行直达----价格：' num2str(time_sort(i,4)) '----时间：' num2str(time_sort(i,5))])
end
end
end
if direct==0 disp('没有直达的车'); end
disp('-----------------------------换乘一次-----------------------------------')
result_mat1=[];result_mat2=[];
result2=[];result11=[];direct1=0;
for i=1:3957
if SS_matrix_con(startP,i)==1&&SS_matrix_con(i,endP)==1
direct1=1;
%disp('In')
for j1=1:520
for k=1:2
if SLmat1(startP,j1,k)>0&&SLmat1(i,j1,k)>0&&SLmat1(i,j1,k)-SLmat1(startP,j1,k)>0
%disp(num2str(j1))
distance1=SLmat1(i,j1,k)-SLmat1(startP,j1,k);
money1=calc_money(j1,distance1);
time1=distance1*3;
result10=[j1 distance1 k money1 time1 i];
result_mat1=[result_mat1(:,:);result10];
%disp('A')
end
end
end
%size( result_mat1)
%disp('q')
%result_mat1
for j2=1:520
for k=1:2
%disp('B')
if SLmat1(i,j2,k)>0&&SLmat1(endP,j2,k)>0&&(SLmat1(endP,j2,k)-SLmat1(i,j2,k)>0)
%disp('B')
distance2=SLmat1(endP,j2,k)-SLmat1(i,j2,k);
money2=calc_money(j2,distance2);
time2=distance2*3;
result11=[j2 distance2 k money2 time2];
result_mat2=[result_mat2(:,:);result11];
%result_mat2
end
end
end
%size( result_mat2)
%disp_str=(['先通过L' num2str(j) '经过' num2str(distance) '个站点上行到S' num2str(i) ','] );
for m=1:size(result_mat1,1)
for n=1:size(result_mat2,1)
result12=[result_mat1(m,:) result_mat2(n,:)];
result2=[result2;result12];
end
end
%result2
end
result_mat1=[];result_mat2=[];%重要的一句
end
if direct1~=0
result_disp=[result2(:,:),result2(:,4)+result2(:,10),result2(:,5)+result2(:,11)+5];
%----------------------------------------
if size(result_disp,1)~=0
money_sort=sortrows(result_disp,12);%按钱排序
time_sort=sortrows(result_disp,13);%按时间排序
if size(money_sort,1)<5
hang=size(money_sort,1);
else hang=5;
end
disp('按价格排序的最佳路线')
for i=1:hang
if money_sort(i,3)==1 w1=['上行'];
else w1=['下行'];
end
if money_sort(i,9)==1
w2=['上行'];
else w2=['下行'];
end
disp(['先通过L' num2str(money_sort(i,1)) '经过' num2str(money_sort(i,2)) '个站点' w1 '到S' num2str(money_sort(i,6)) '，再通过L' num2str(money_sort(i,7)) '经过' num2str(money_sort(i,8)) '个站点' w2 '到终点----费用：' num2str(money_sort(i,12)) '----时间：' num2str(money_sort(i,13))])
end
disp('按时间排序的最佳路线')
for i=1:hang
if time_sort(i,3)==1
w1=['上行'];
else w1=['下行'];
end
if time_sort(i,9)==1
w2=['上行'];
else w2=['下行'];
end
disp(['先通过L' num2str(time_sort(i,1)) '经过' num2str(time_sort(i,2)) '个站点' w1 '到S' num2str(time_sort(i,6)) '，再通过L' num2str(time_sort(i,7)) '经过' num2str(time_sort(i,8)) '个站点' w2 '到终点----费用：' num2str(time_sort(i,12)) '----时间：' num2str(time_sort(i,13))])
end
end
else
disp('不能换乘一次')
end
disp('-----------------------------换乘两次-----------------------------------')
result2_mat1=[];result2_mat2=[];result2_mat3=[];
result3=[];result31=[];
for i=1:3957
for j=1:3957
if SS_matrix_con(startP,i)==1&&SS_matrix_con(i,j)&&SS_matrix_con(j,endP)==1
%第一个换乘点
for k=1:520
for m=1:2
if SLmat1(startP,k,m)>0&&SLmat1(i,k,m)>0&&(SLmat1(i,k,m)-SLmat1(startP,k,m)>0)
distance1=SLmat1(i,k,m)-SLmat1(startP,k,m);
money1=calc_money(k,distance1);
time1=distance1*3;
result20=[k distance1 m money1 time1 i];
result2_mat1=[result2_mat1(:,:);result20];
%disp_str=(['先通过L' num2str(k) '经过' num2str(distance) '个站点上行到S' num2str(i) ','] );
end
end
end
%第二个换乘点
for k=1:520
for m=1:2
if SLmat1(i,k,m)>0&&SLmat1(j,k,m)>0&&(SLmat1(j,k,m)-SLmat1(i,k,m)>0)
distance2=SLmat1(j,k,m)-SLmat1(i,k,m);
money2=calc_money(k,distance2);
time2=distance2*3;
result21=[k distance2 m money2 time2 j];
result2_mat2=[result2_mat2(:,:);result21];
%disp_str=([disp_str '再在S' num2str(i) '通过L' num2str(k) '经过' num2str(distance) '个站点上行到S' num2str(j) ',']) ;
end
end
end
%第三个换乘点
for k=1:520
for m=1:2
if SLmat1(endP,k,m)>0&&SLmat1(j,k,m)>0&&(SLmat1(endP,k,m)-SLmat1(j,k,m)>0)
distance3=SLmat1(endP,k,m)-SLmat1(j,k,m);
money3=calc_money(k,distance3);
time3=distance3*3;
result22=[k distance3 m money3 time3];
result2_mat3=[result2_mat3(:,:);result22];
%disp([disp_str '再在S' num2str(j) '通过L' num2str(k) '经过' num2str(distance) '个站点上行到终点' '----价格：' num2str(money) '----时间：' num2str(time)])
end
end
end
%拼接矩阵
if size(result2_mat1,1)~=0&&size(result2_mat2,1)~=0&&size(result2_mat3,1)~=0
for m1=1:size(result2_mat1,1)
for n1=1:size(result2_mat2,1)
for k1=1:size(result2_mat3,1)
result22=[result2_mat2(n1,:) result2_mat3(k1,:)];
result31=[result2_mat1(m1,:) result22];
result3=[result3(:,:);result31];
end
end
end
end
end
result2_mat2=[];result2_mat3=[]; result2_mat1=[];
end
end
result2_disp=[result3(:,:),result3(:,4)+result3(:,10)+result3(:,16),result3(:,5)+result3(:,11)+result3(:,17)+10];
if size(result2_disp,1)~=0
money_sort=sortrows(result2_disp,18);%按钱排序
time_sort=sortrows(result2_disp,19);%按时间排序
if size(money_sort,1)<5
hang=size(money_sort,1);
else hang=5;
end
disp('按价格排序的最佳路线')
for i=1:hang
if money_sort(i,3)==1 w1=['上行'];
else w1=['下行'];
end
if money_sort(i,9)==1
w2=['上行'];
else w2=['下行'];
end
if money_sort(i,15)==1
w3=['上行'];
else w3=['下行'];
end
disp(['先通过L' num2str(money_sort(i,1)) '经过' num2str(money_sort(i,2)) '个站点' w1 '到S' num2str(money_sort(i,6)) '，再通过L' num2str(money_sort(i,7)) '经过' num2str(money_sort(i,8)) '个站点' w2 '到S' num2str(money_sort(i,12)) ',再通过L' num2str(money_sort(i,13)) '经过' num2str(money_sort(i,14)) '个站点' w3 '到终点----费用：' num2str(money_sort(i,18)) '----时间：' num2str(money_sort(i,19))])
end
disp('按时间排序的最佳路线')
for i=1:hang
if time_sort(i,3)==1 w1=['上行'];
else w1=['下行'];
end
if time_sort(i,9)==1
w2=['上行'];
else w2=['下行'];
end
if time_sort(i,15)==1
w3=['上行'];
else w3=['下行'];
end
disp(['先通过L' num2str(time_sort(i,1)) '经过' num2str(time_sort(i,2)) '个站点' w1 '到S' num2str(time_sort(i,6)) '，再通过L' num2str(time_sort(i,7)) '经过' num2str(time_sort(i,8)) '个站点' w2 '到S' num2str(time_sort(i,12)) ',再通过L' num2str(time_sort(i,13)) '经过' num2str(time_sort(i,14)) '个站点' w3 '到终点----费用：' num2str(time_sort(i,18)) '----时间：' num2str(time_sort(i,19))])
end
end

toc