公交换乘

数学建模公交换乘问题的matlab程序共享乘公交看奥运，下面的程序是找出最佳乘公交的路径，直达的线路，一次转乘或两次转乘的线路题目：乘公交，看亚运（中找最佳路线问题）2010年11月12日第16届亚运会在广州举行，为了让全体市民更好看亚运会，广州市政府决定在亚运期间放假3天、以及全体市民可在亚运及亚残运会期间免费坐公交、地铁30个工作日等惠民政策，这一政策的施行在很大程度上加剧了广州市交通出行的困难。

附件一，附件二是广州市公交（包括地铁）的一些信息，为了方便游客看亚运会，请你用数学建模的方法，为广州市设计一个公交线路查询系统，满足查询者的各种不同需求。

请你们解决如下问题：1、在亚运会开幕前，仅考虑公汽线路，给出任意两公汽站点之间线路选择问题的一般数学模型与算法。

并根据附录数据，利用你们的模型与算法，求出以下3对起始站→终到站之间的最佳路线（要有清晰的评价说明）。

(1)、华穗路→交通大厦(2)、越秀桥→山村(3)、江南大道北→策边村程序中用到的lx0.xls数据下载/view/bc9efa4de45c3b3567ec8bfd.html原始数据（原始线路信息、站点编号表、公交说明）附录二：线路.xls 下载地址/view/ff5a660b844769eae009ed94.html有问题QQ869896537尽量解决如果你的题目的数据不是这样的，即数据的量不同，你需要在程序中修改相应数据，例如循环次数的修改等数据Matlab 源程序代码：clear all;clc[a b c]=xlsread('lx0.xls');biao=unique(b);for i=1:length(biao)b(strcmp(b,biao(i)))={i-1};endbiao=biao(2:end);%公交站点编号对照表zhan=cell2mat(b);%公交数字站点信息for i=1:1108fl=zhan(i,:);fl(find(fl==0))=[];if length(fl)>47;zhan(i,:)=0;%去掉地铁线路endif length(fl)<6;qian(i)=2.5;elseif length(fl)<10&&length(fl)>=6;qian(i)=2.0;elseif length(fl)>=10;qian(i)=1.5;endendcx1=inputdlg('起站点：');cx2=inputdlg('到站点：'); for i=1:length(biao)if ismember(cx1,biao(i));shuzizhandian1=i; %将站名换成数字站点endif ismember(cx2,biao(i));shuzizhandian2=i;endendnum1=0;num2=0;for i=1:1:1108for j=1:55y=zhan(i,j);if shuzizhandian1==y;num1=num1+1;gj(num1)=i;%经过起站点的公交号zd(num1)=j;endif shuzizhandian2==y;num2=num2+1;gjend(num2)=i;%经过到起站的公交号zdend(num2)=j;endendendstart=[gj;zd]'; %起始站点可以选择的公交及在公交路线中的位置des=[gjend;zdend]'; %目的地可供选择的公交信息zhi=intersect(gj,gjend)' ; %可选的直达车信息cy=length(zhi); %cy为直达车的个数%******************……下面按照转乘车次数安排线路……***********************%%%clc;if cy~=0;for i=1:cyz1=[];z2=[];z1=find(gj==zhi(i));z2=find(gjend==zhi(i));% if zdend(z2(1,1))-zd(z1(1,1))<0;% zhi(i)=[];% endzs(i)=zdend(z2(1,1))-zd(z1(1,1));endfprintf(' 《《《《公交直达选择信息》》》》\n')xuanzhe=[zhi,zs'];you=sortrows(xuanzhe,2);%优先选择排序for i=1:cyww=zhan(you(i,1),:); %直达线路fprintf('%4d==\t<%3d站>\t',you(i,1),you(i,2));if find(ww==shuzizhandian1)-find(ww==shuzizhandian2)>0;forj=find(ww==shuzizhandian2):find(ww==shuzizhandian1)fprintf('<-- %s',biao{ww(j)});endfprintf('\n');elseforj=find(ww==shuzizhandian1):find(ww==shuzizhandian2)fprintf('--> %s',biao{ww(j)});endfprintf('\n');endendelse%*********************……下面求转乘一次的线路……*************************%fprintf('\n============================================= ==================================\n');fprintf(' 可惜！没有直达车！！\n 亲！可以选择换乘哦\n');ji=0 ;%第ji个可以换乘的公交length(gj);length(gjend);zc=[0];for i=1:length(gj)m=gj(i);for j=1:length(gjend)n=gjend(j);p1=zhan(m,:);p1(find(p1==0))=[];p2=zhan(n,:);p2(find(p2==0))=[];if intersect(p1,p2);ji=ji+1;hc(ji,:)=[m,n]; %转乘的两公交t=intersect(p1,p2); %转乘交点站zc(ji,1:length(t))=t;%转乘信息收集zc(ji,(length(t)+1):5)=0; %扩充来使数据集中放置endendendif max(max(zc))~=0;fprintf('\n============================================= ==================================\n');fprintf(' 可以选择换乘一次\n');patch=[hc ,zc];s=0;for i=1:jizc1=zc(i,:);zc1(find(zc1==0))=[];for j=1:length(zc1)s=s+1;patch1(s,:)=[hc(i,:),zc1(j)];endendfprintf(' \n\t 首坐车\t 转乘车\t 中途下车站点\n');disp(patch1) %for i=1:see1=zhan(patch1(i,1),:); ee2=zhan(patch1(i,2),:);iffind(ee1==shuzizhandian1)-find(ee1==patch1(i,3))<0;long=abs(find(ee1==shuzizhandian1)-find(ee1==patch1(i,3)));fprintf('%4d==\t<%3d站>\t',patch1(i,1),long);t1=long*3;forj=find(ee1==shuzizhandian1):find(ee1==patch1(i,3))fprintf('--> %s',biao{ee1(j)});endfprintf('\n');elselong=abs(find(ee1==shuzizhandian1)-find(ee1==patch1(i,3)));fprintf('%4d==\t<%3d站>\t',patch1(i,1),long);t1=long*3;forj=find(ee1==shuzizhandian1):-1:find(ee1==patch1(i,3))fprintf('<-- %s',biao{ee1(j)});endfprintf('\n');endiffind(ee2==shuzizhandian2)-find(ee2==patch1(i,3))>0;long=abs(find(ee2==shuzizhandian2)-find(ee2==patch1(i,3)));t2=long*3;fprintf('%4d==\t<%3d站>\t',patch1(i,2),long);forj=find(ee2==patch1(i,3)):find(ee2==shuzizhandian2)fprintf('--> %s',biao{ee2(j)});endfprintf('\n');elselong=abs(find(ee2==shuzizhandian2)-find(ee2==patch1(i,3)));t2=long*3;fprintf('%4d==\t<%3d站>\t',patch1(i,2),long);forj=find(ee2==patch1(i,3)):-1:find(ee2==shuzizhandian2)fprintf('<-- %s',biao{ee2(j)});endfprintf('\n');endlufei=qian(patch1(i,1))+qian(patch1(i,2));t(i)=t1+t2+5;T=t1+t2+5;fprintf('!选择（%2d）===========',i);fprintf('!路费=%3d元',lufei);fprintf('!\t\ttime=%3d分钟\n\n',T);disp('-------------------------------------------------------------------------------------------------------------------------------');endzuikuai=find(t==min(t));fprintf('!最快的选择是————》选择（%3d）\n',zuikuai);disp('-------------------------------------------------------------------------------------------------------------------------');else%*********************……下面求转乘两次的线路……*************************%fprintf('\n============================================= ==================================\n');fprintf(' 换乘一次无法找到\n');buyao=union(gj,gjend);%经过起点或终点的公交quan=[1:1108];disan=setdiff(quan,buyao);%可供选择的转乘车ds2=[];ds3=[];for j=1:length(gj)ds2=union(zhan(gj(j),:),ds2);%经过起点的车可以覆盖的地点ds2(find(ds2==0))=[];endfor k=1:length(gjend)ds3=union(zhan(gjend(k),:),ds3);%经过终点的车可以覆盖的地点ds3(find(ds3==0))=[];endl=0;%中间车计数量zj1=0;zuihou=0;for i=1:length(disan)%转乘车数量ds1=[];ds1=zhan(disan(i),:);ds1(find(ds1==0))=[];if intersect(ds1,ds3);%寻找中间车if intersect(ds1,ds2);l=l+1;zjc(l)=disan(i);%可考虑的中间转乘车车号for j=1:length(gj)jj=zhan(gj(j),:);jj(find(jj==0))=[];if intersect(ds1,jj);j1=intersect(ds1,jj);%第1次下车站点% zj1=zj1+1;%xx(zj1,:)=[gj(j),disan(i),j1(1)] ;%第一二辆车及转点for k=1:length(gjend)kk=zhan(gjend(k),:);kk(find(kk==0))=[];if intersect(ds1,kk);j2=intersect(ds1,kk);%第2次下车站点zj1=zj1+1;xx(zj1,:)=[gj(j),j1(1),disan(i),j2(1),gjend(k)];zuihou=1;endendendendendendendif zuihou==0;fprintf('\n============================================= ==================================\n');fprintf(' sorry没有两次以内的转乘安排\n');elsefprintf('\n============================================= ==================================\n');fprintf(' 为你安排的两次转乘路线：\n');disp(' 起点车一下车站点中间车二下车站点到终点车');for i=1:zj1ss1=zhan(xx(i,1),:); ss2=zhan(xx(i,3),:); ss3=zhan(xx(i,5),:);iffind(ss1==xx(i,2))-find(ss1==shuzizhandian1)<0||find(ss2==xx(i,4))-find (ss2==xx(i,2))<0||find(ss3==shuzizhandian2)-find(ss3==xx(i,4))<0;xx(i,:)=0;%去除逆行线路endendxx(all(xx==0,2),:)=[];disp(xx);if length(xx(:,1))==0;fprintf('没有找到正向线');elsefor i=1:length(xx(:,1))ss1=zhan(xx(i,1),:); ss2=zhan(xx(i,3),:); ss3=zhan(xx(i,5),:);long1=find(ss1==xx(i,2))-find(ss1==shuzizhandian1);fprintf('%4d==\t<%3d站>\t',xx(i,1),long1);forj=find(ss1==shuzizhandian1):find(ss1==xx(i,2))fprintf('--> %s',biao{ss1(j)});endfprintf('\n');long2=find(ss2==xx(i,4))-find(ss2==xx(i,2));fprintf('%4d==\t<%3d站>\t',xx(i,1),long2);for j=find(ss2==xx(i,2)):find(ss2==xx(i,4))fprintf('--> %s',biao{ss2(j)});endfprintf('\n');long3=find(ss3==shuzizhandian2)-find(ss3==xx(i,4));fprintf('%4d==\t<%3d站>\t',xx(i,1),long3);forj=find(ss3==xx(i,4)):find(ss3==shuzizhandian2)fprintf('--> %s',biao{ss3(j)});endlufei2=qian(xx(i,1))+qian(xx(i,3))+qian(xx(i,5));% timetwo(i)=3*(long1+long2+long3)+10;T2=3*(long1+long2+long3)+10;fprintf('\n线路《%3d》========!路费=%d元',i,lufei2);fprintf('!\t\ttime=%3d分钟\n\n',T2);fprintf('\n\n');endzuikuai2=find(timetwo==min(timetwo));fprintf('!最快的选择是————》选择线路《%3d》\n',zuikuai2);fprintf('!最短时间：《%3d分钟》\n',min(timetwo(zuikuai2)));disp('-------------------------------------------------------------------------------------------------------------------------');endendendend。