【正文】 ast(x(k+1,:)))]。 l=length(xx(1,:))。 for i=1:(l1) j=1。 while i+j=l if w(xx(1,i),xx(1,i+j))1 w(xx(1,i),xx(1,i+j))=1。 end j=j+1。 end end elseif x(k,3)==2 %下行 elseif x(k,3)==3 %環(huán)線 xx=[x(k,4:last(x(k,:))1) x(k,4:last(x(k,:))1)]。 l=length(xx(1,:))/2。 for i=1:l j=1。 while i+j=length(xx(1,:)) if w(xx(1,i),xx(1,i+j))1 w(xx(1,i),xx(1,i+j))=1。 end j=j+1。 end end end endendfor i=1:length(w(:,1)) w(i,i)=0。end w1_2=w。save(39。C:\Documents and Settings\wly\桌面\新建文件夾\1\b\B2007data\w1 2\39。,39。w1_239。)程序3：計(jì)算問題一模型三時(shí)間費(fèi)用率矩陣qijclc。clear。x=load(39。C:\Documents and Settings\wly\桌面\新建文件夾\1\b\B2007data\39。)。w=inf(3957,3957)。for k=1:929 clear xx。 if x(k,3)==0 %往返 xx=[x(k,4:last(x(k,:))) fliplr(x(k,4:last(x(k,:))1))]。 l=length(xx(1,:))。 for i=1:(l1) for j=(i+1):l q=(ji)*3。 if w(xx(1,i),xx(1,j))q w(xx(1,i),xx(1,j))=q。 end end end elseif x(k,3)==1 %上行 xx=[x(k,4:last(x(k,:))) x(k+1,4+1:last(x(k+1,:)))]。 l=length(xx(1,:))。 for i=1:(l1) for j=(i+1):l q=(ji)*3。 if w(xx(1,i),xx(1,j))q w(xx(1,i),xx(1,j))=q。 end end end elseif x(k,3)==2 %下行 elseif x(k,3)==3 %環(huán)線 xx=[x(k,4:last(x(k,:))1) x(k,4:last(x(k,:))1)]。 l=length(xx(1,:))/2。 for i=1:l j=1。 while i+j=length(xx(1,:)) q=j*3。 if w(xx(1,i),xx(1,i+j))q w(xx(1,i),xx(1,i+j))=q。 end j=j+1。 end end end kendw=w+5。for i=1:length(w(:,1)) w(i,i)=0。endw1_3=w。save(39。C:\Documents and Settings\wly\桌面\新建文件夾\1\b\B2007data\w1_3_時(shí)間\39。,39。w1_339。)程序4：改進(jìn)的Dijkstra程序clc。clear。%%設(shè)定虛擬節(jié)點(diǎn)load(39。C:\Documents and Settings\wly\桌面\新建文件夾\1\b\B2007data\w1_3_時(shí)間\39。)。load(39。D:\39。)。w=w1_3。clear w1_3。startp=3359。endp=1828。l1=length(w(1,:))+1。l2=length(w(1,:))+2。for i=l1:l2 w(i,i)=0。 for j=1:l2 w(i,j)=inf。 w(j,i)=inf。 endendw(l1,startp)=0。w(endp,l2)=0。%%改進(jìn)的Dijkstra算法m=length(wdao)。n=size(w,1)。w1=w(l1,:)。l(1,1:n)=0。ll=l。z(1,1:n)=0。for i=1:n l(i)=w1(i)。 z(i)=l1。endclear w1。s(1,1:n)=0。s(1)=l1。u=s(1)。k=1zhuan=zeros(1,l2)1。%轉(zhuǎn)車次數(shù) 1為最大2次 0為最大1次 +1為直達(dá)while kn if zhuan(1,u)=1 amp。amp。 l(1,u)=l(1,endp) %只更改一次轉(zhuǎn)車之內(nèi)能到達(dá)的節(jié)點(diǎn) for i=1:m for j=1:k tem=wdao(i)。 if tem~=s(j) amp。amp。 l(tem)l(u)+w(u,tem) l(tem)=l(u)+w(u,tem)。 zhuan(1,tem)=zhuan(1,u)+1。 z(tem)=u。 end end end ll=l。 for i=1:m for j=1:k if wdao(i)==s(j) ll(wdao(i))=inf。 end end end else ll(u)=inf。 end lv=inf。 v=startp。 %只探測(cè)一次轉(zhuǎn)車之內(nèi)能到達(dá)的節(jié)點(diǎn) for i=1:m tem=wdao(i)。 if ll(tem)lv lv=ll(tem)。 v=tem。 end end s(k+1)=v。 if v==endp break。 end k=k+1 u=s(k)。endl(1,end)save(39。d:\39。)。程序5：上行與下行起點(diǎn)終點(diǎn)是否相接clc。clear。x=load(39。C:\Documents and Settings\wly\桌面\新建文件夾\1\b\B2007data\39。)。c=1。for k=2:929 last1=4。 last2=4。 if x(k,3)==2 while x(k1,last1+1)~=0 last1=last1+1。 end while x(k,last2+1)~=0 last2=last2+1。 end test(c,1)=(x(k,4)==x(k1,last1))。 test(c,2)=(x(k1,4)==x(k,last2))。 c=c+1。 endendtestt=sum(test)。程序6：用計(jì)算實(shí)例的結(jié)果計(jì)算轉(zhuǎn)車路徑（以w1 m3為例）clear。load(39。D:\39。)s=length(z(1,:))1。t=length(z(1,:))。ii=1。while z(t)~=s lj(ii)=t。 ii=ii+1。 t=z(t)。endlj(ii)=t。lj(ii+1)=z(t)。lj=fliplr(lj)。lj1_3_1=lj。save(39。d:\39。,39。lj1_3_139。)%、終點(diǎn)的轉(zhuǎn)車路徑。其它程序要么與已附程序類似，要么重要性不大，限于篇幅在此不全部列出。34