【正文】 end end for i=1:92 for j=1:92 if Y(i,j)~=0 Y(i,j)=1。 end if i==39|i==48|i==58|i==67|i==87|i=20 Y(i,i)=1。 end end end for j=1:92 Y(:,j)=r(j)*Y(:,j)。 end s=sum(Y,2) s=nonzeros(s) d=var(s)。 for j=1:92 [h(j),H(j)]=max(Y(:,j))。 end for i=1:92 z(i)=x(i,H(i)) end 8．2．1求解問題二第一問的MATLAB程序clc,clearx=xlsread(39。39。,39。sheet139。,39。A1:BU7339。)r=xlsread(39。39。,39。全市交通路口節(jié)點數(shù)據(jù)39。,39。E94:E16639。)n=length(x)。g=8。h=zeros()。 Y=zeros(n,n)。 D=zeros(n,n)。for i=1:n %挑出1g 行的數(shù)據(jù) if i=g for j=1:n D(i,j)=x(i,j)。 end end end for i=1:n for j=1:n if D(i,j)==0 D(i,j)=10000。 end end endfor j=1:n%找出里j點最近的平臺作為管理平臺。 [y(j),I(j)]=min(D(:,j))。end for j=1:n Y(I(j),j)=y(j)。 end for j=1:n if j=g %剔除平臺相互管理管理 Y(:,j)=0 end end for i=1:n %把Y矩陣變成01矩陣，用于發(fā)案率求和 for j=1:n if Y(i,j)~=0 Y(i,j)=1。 end if i=g Y(i,i)=1。 end end end for j=1:n Y(:,j)=r(j)*Y(:,j)。 end s=sum(Y,2)%將i點所管理的所有j點發(fā)案率累加求和 s=nonzeros(s)%去除0點 d=var(s)。%對工作量（發(fā)案率）求方差 for j=1:n [h(j),H(j)]=max(Y(:,j))。%輸出ij對應點和其j發(fā)案率 end for i=1:n z(i)=x(i,H(i))%對出警時間 end（改進之前求解方差的程序）clc,clearx=xlsread(39。39。,39。sheet139。,39。A1:BU7339。)r=xlsread(39。39。,39。全市交通路口節(jié)點數(shù)據(jù)39。,39。E94:E16639。)n=length(x)。g=8。h=zeros()。 Y=zeros(n,n)。 D=zeros(n,n)。for i=1:n %挑出1g，m,n 行的數(shù)據(jù) if i=g for j=1:n D(i,j)=x(i,j)。 end end if i==12392|i==14592|i==15592 for j=1:n D(i,j)=x(i,j)。 end endend for i=1:n for j=1:n if D(i,j)==0 D(i,j)=10000。 end end endfor j=1:n [y(j),I(j)]=min(D(:,j))。end for j=1:n Y(I(j),j)=y(j)。 end for j=1:n if j=g|j==12392|j==14592|j==15592 %剔除平臺相互管理管理 Y(:,j)=0。 end end for i=1:n for j=1:n if Y(i,j)~=0 Y(i,j)=1。 end if i==12392|i==12492|i==15592|i=g Y(i,i)=1。 end end end for j=1:n Y(:,j)=r(j)*Y(:,j)。 end s=sum(Y,2)。 s=nonzeros(s)。 d=var(s)。 for j=1:n [h(j),H(j)]=max(Y(:,j))。 end for i=1:n z(i)=x(i,H(i))。 end（改進之后求解方差的程序） 27