数学建模2：评价类模型

clc;clear;
data=xlsread('day_20.xls');%第一列是借车站号；第二列是还车站号；第三列是用车时长
%P为任意两点间的最小时间
min_time=ones(181)*inf;
for i=1:size(data,1)-1
    if min_time(data(i,1),data(i,2))==inf
        min_time(data(i,1),data(i,2))=data(i,3);
    end
end
min_time;
for i=1:181
    for j=i:181
        min_time(i,j)=min(min_time(i,j),min_time(j,i));
        min_time(j,i)=min(min_time(i,j),min_time(j,i));
    end
end
P=floyd(min_time)
xlswrite('day_20_P.xls',P)

function D=floyd(w)
D=w;
n=length(w);
path=zeros(n);
for i=1:n
    for j=1:n
        if D(i,j)~=inf
            path(i,j)=j;
        end
    end
end
for k=1:n
    for i=1:n
        for j=1:n
            if D(i,k)+D(k,j)<D(i,j)
                D(i,j)=D(i,k)+D(k,j);
                path(i,j)=path(i,k);
            end
        end
    end
end

clc;clear;
x1=xlsread('每时间段各站点的借车频次.xls','Sheet1','B3:BI182');%各时段各站点借车频次
x2=xlsread('每时间段各站点的还车频次.xls','Sheet1','B3:BI182');%各时段各站点还车频次
y1=duiqi(x1);
gaofeng1=gaofengqi(y1);
y2=duiqi(x2);
gaofeng2=gaofengqi(y2);
%第一列为站号，第二列为最大值，第三列为最大值对应时段，
%第四列为次大值，第五列为次大值对应时段，
%各时段对应值为：1、6-6:30；2、6：30-7点；3、7-7：30点；。。。。依次类推。。。30、20-21点；16：21-22点
xlswrite('day_20_gaofeng_jie',gaofeng1);
xlswrite('day_20_gaofeng_huan',gaofeng2);

function y=duiqi(x)
m=[1:181]';%m的第一列为站号,以后各列为各情况下各站统计量
for i=1:size(x,1)
    for j=1:size(x,2)/2
        if x(i,j*2-1)~=0
        m(x(i,j*2-1),j+1)=x(i,j*2);
        end
    end
end
y=m;

function gaofeng=gaofengqi(y)
x=zeros(size(y,1),size(y,1)-1);
x=y(:,2:end);
gaofeng=[];
for i=1:size(x,1)
    [max1,index1]=max(x(i,:));
    [max2,index2]=max([setdiff(x(i,:),x(i,index1)),0]);
    gaofeng=[gaofeng;i,max1,index1,max2,index2];
end

function output=grayrela
x0=xlsread('参考因子与比较因子.xls','Sheet1','B2:F182')'
%计算灰色关联系数
%参考因子与比较因子共同存储在一个矩阵x0中,参考因子位于第一列
%斜率序列
for i=2:length(x0(:,1))
    x1(i,:)=x0(i,:)-x0(i-1,:);
end
%标准化
m=length(x1(1,:));
for i=1:m
x2(:,i)=x1(:,i)/std(x1(:,i));
end
%排序
[y,pos]=sort(x2(:,1));
x2_sorted=x2(pos,:);
% 判定关联性质
n=length(x1(:,1));
k=[1:n]';
for j=1:m
sig_j(j)=qiuhe(k.*x2_sorted(:,j))-qiuhe(x2_sorted(:,j))*qiuhe(k)/n;
end
%caculation of distantion
for j=2:m
    dist_0i(:,j)=abs(sign(sig_j(:,j)./sig_j(:,1)).*x2_sorted(:,j)-x2_sorted(:,1));
end
%计算关联系数
for i=1:n
    for j=1:m
        coef_rela(i,j)=(min(dist_0i)+0.5*max(dist_0i))/(dist_0i(i,j)+0.5*max(dist_0i));
    end
end
for j=1:m
    output(j)=qiuhe(coef_rela(:,j))/n;
end
xlswrite('day_20_Coef',output);
% function output=qiuhe(input)
% output=0;
% for i=1:length(input)
%     output=output+input(i);
% end

function output=qiuhe(input)
output=0;
for i=1:length(input)
    output=output+input(i);
end