matlab 公交网络矩阵程序,【数学建模集训系列】公交查询系统的matlab实现-站点和站点之间直达矩阵...

%==============================================================

% 功能: 输出两个矩阵，分别是站点之间的联通矩阵、站点经过的线路统计

% 日期: 8/8/2011

%==============================================================

clear,clc

tic %计时开始

fid=fopen('Bus.txt','r');%打开数据

if fid>0 disp('数据文件打开成功！')

else disp('打开失败'),return

end

a=[];

b=[];

maxS=0;%初始化最大站点数

%计算最大站点数

while ~feof(fid)

a=fgetl(fid);

if isempty(a)

%disp('原路返回')

else

n=find(a=='S');

for i=1:length(n)

if str2num(a(n(i)+1:n(i)+4))>maxS

maxS=str2num(a(n(i)+1:n(i)+4));

end

end

end

end

disp(['站点数量为: ' num2str(maxS)])

clear a b maxS

%------------------------------------------------------------

%fclose(fid);%关闭文件

%fid=fopen('Bus.txt','r');%重新打开文件

fseek(fid,0,-1);%修改文件指针，跳到文件头

%fgetl(fid)

%---------------------SS_matrix--------------------------

%SS_matrix_con表示是否直接连通，SS_matrix_tim表示时间

%SS_matrix_mon表示价格,SS_matrix_L,表示连通两站点间的公交路线

%SS_matrix_con=zeros(maxS);%站点间联通矩阵，初始化全零

a=[];

temp_L=[];

while ~feof(fid)

a=fgetl(fid);

%记录当前L

n_L=find(a=='L');

if ~isempty(n_L)

temp_L=str2num(a(n_L+1:n_L+3));

end

%temp_L

n=find(a=='S');

if ~isempty(n)

for i=1:length(n)

for j=i:length(n)

a1=str2num(a(n(i)+1:n(i)+4));

a2=str2num(a(n(j)+1:n(j)+4));

SS_matrix_con(a1,a2)=1; %SS连通矩阵

%SS_matrix_L(a1,a2)=temp_L;%SS的线路矩阵

if a1~=a2

SS_matrix_tim(a1,a2)=(j-i)*3;

end

end

end

end

end

%------------处理原路返回的路线----------------

temp_L=1;

n=[];

fseek(fid,0,-1);%修改文件指针，跳到文件头

for i=3:(520*4-1)

temp=fgetl(fid);

%记录当前L

n_L=find(a=='L');

if ~isempty(n_L)

temp_L=str2num(a(n_L+1:n_L+3));

end

if ~isempty(temp)

if temp(1)=='S'

n=find(temp=='S');

for i=length(n):-1:1

for j=i:-1:1

a1=str2num(temp(n(i)+1:n(i)+4));

a2=str2num(temp(n(j)+1:n(j)+4));

SS_matrix_con(a1,a2)=1;

SS_matrix_tim(a1,a2)=(i-j)*3;

%SS_matrix_L(a1,a2)=temp_L;%SS的线路矩阵

end

end

end

end

end

%--------------------------------------------

%disp('站点之间互相连通的稀疏矩阵为')

%SS_xishu=sparse(SS_matrix);

clo=fclose(fid);%关闭数据文件

if clo==0 disp('数据文件成功关闭！')

else disp('数据文件关闭失败！'),return

end

save D:\matlab_files\SSdata_con SS_matrix_con

disp('站点S间的连通矩阵已经保存为D:\matlab_files\SSdata_con.mat')

save D:\matlab_files\SSdata_tim SS_matrix_tim

disp('站点S间的时间矩阵已经保存为D:\matlab_files\SSdata_tim.mat')

%save D:\matlab_files\SSdata_L SS_matrix_L

%disp('站点S间的线路矩阵已经保存为D:\matlab_files\SSdata_L.mat')

toc %结束计时

clear