function?Z=floyd（ady）
??????????????n=size（a1）;???????????????%矩陣維數
??????????????Z=zeros（n1）;??????????????%節點重要度向量???????????????????????????
??????????????%由鄰接矩陣a得到直接距離矩陣H
??????????????H=zeros（size（a））;
??????????????for?i=1:n
??????????????????for?j=1:n
??????????????????????if?j==i
??????????????????????????H（ij）=0;
??????????????????????elseif?a（ij）==1
??????????????????????????H（ij）=1;
??????????????????????else
??????????????????????????H（ij）=inf;
??????????????????????end
??????????????????end
??????????????end
??????????????
??????????????%用Floyd法計算節點收縮前的最短距離矩陣D
??????????????D=H;???????????????????????????
??????????????for?k=1:n
??????????????????for?i=1:n
??????????????????????for?j=1:n
??????????????????????????if?D（ik）+D（kj）??????????????????????????????D（ij）=D（ik）+D（kj）;?????????????????????????????
??????????????????????????end?
??????????????????????end?
??????????????????end????????????????????????????????????
??????????????end
??????????????
??????????????%計算節點重要度
??????????????D2=zeros（size（D））;??????????????
??????????????for?i=1:n??????????????????
??????????????????%得到與節點i鄰接的節點向量I
??????????????????I=zeros（10）;
??????????????????t=0;
??????????????????for?j=1:n
??????????????????????if?a（ij）==1
??????????????????????????t=t+1;
??????????????????????????I=[Ij];
??????????????????????end
??????????????????end
??????????????????
????????????????%計算收縮前的最短距離矩陣D2
????????????????for?p=1:n
????????????????????for?q=1:n
????????????????????????if?p~=i&q~=i
????????????????????????????if?D（pi）+D（iq）==D（pq）
????????????????????????????????D2（pq）=D（pq）-2;
????????????????????????????elseif?D（pi）