选址matlab程序

function WaitBarUseage
h = waitbar（0‘Please wait...‘）;
for i=1:100
    waitbar（i/100h[num2str（i）‘%‘]）
    pause（.1）
end
close（h）


n=49;  %总共49个点
A=zeros（nn）;

for i=1:n
    for j=1:n
       if（i==j） A（ij）=0;
       else A（ij）=100000;
       end
    end
end  
A（12）=120;A（13）=270; A（15）=540;A（16）=799; A（115）=420;A（140）=844;
A（23）=210;A（215）=360;A（34）=210;A（315）=311;  A（316）=440; A（45）=530;
A（416）=430;A（427）=630;A（430）=760;A（530）=720; A（540）=1521; A（547）=186;
A（67）=330; A（639）=387; A（640）=727; A（78）=230; A（740）=429;
A（741）=347;A（842）=819; A（910）=280;  A（911）=190; A（915）=840;
A（1011）=279; A（1012）=160; A（1014）=660; A（1015）=680; A（1038）=598;
A（1043）=325; A（1113）=880; A（1114）=640; A（1137）=153; A（1214）=610;
A（1216）=650;A（1217）=540;A（1243）=435;A（1312）=680;A（1319）=1020;
A（1332）=490;A（1336）=266;A（1337）=592;A（1417）=270;A（1418）=640;
A（1419）=860;A（1516）=430;A（1538）=361;A（1543）=349;A（1617）=540;
A（1627）=550;A（1643）=473;A（1644）=285;A（1718）=380;A（1744）=406;
A（1745）=362;A（1718）=380;A（1744）=406;A（1745）=362;A（1819）=780;
A（1824）=1010;A（1845）=508;A（1848）=664;A（1920）=710;A（1921）=580;
A（1934）=130;A（1935）=127;A（1936）=688;A（2021）=560;A（2024）=650;
A（2025）=820;A（2048）=305;A（2149）=270;A（2223）=340;A（2223）=340;
A（2224）=490;A（2225）=1090;A（2227）=910;A（2245）=795;A（2325）=990;
A（2326）=2170;A（2327）=920;A（2425）=650;A（2448）=560;A（2526）=2320;
A（2629）=1940;A（2631）=2627;A（2728）=700;A（2730）=640;A（2744）=637;
A（2746）=304;A（2829）=230;A（2830）=500;A（2831）=1980;A（3047）=554;
A（3335）=36;A（3843）=368;A（4041）=304;A（4042）=929;A（4142）=669;
A（4445）=466;A（4647）=541;
 
for j=1:n
    for i=1:j-1
       A（ji）=A（ij）; %使对称
    end
end  



[mn]=size（A）;

B=zeros（mn）;
B=A;


%各区需求量
%road4.m
a=[1232;974;965;358;223;715;753;989;1391;624;677;487;636;495;603;721;834;642;786;693;156;3257;1126;364;531;51
774;323;194;151;234;246;701;55;233;174;568;761;583;317;204;272;948;1150;401;224;217;366;55];
b=[1123584;1000000000;733400;272080;169480;1000000000;457824;1000000000;1000000000;663936;411616;370120;580032;526680
    733248;876736;760608;585504;955776;526680;1000000000;2475320;684608;276640;403560;31008;941184;196384;1000000000
    114760;1000000000;1000000000;1000000000;1000000000;1000000000;1000000000;1000000000;1000000000;1000000000;289104
    1000000000;1000000000;720480;699200;243808;1000000000;1000000000;1000000000;1000000000];


%利用Flod算法计算最短距离矩阵
for k=1:n    
  for i=1 :n
     for j=1:n        
          t=B（ik）+B（kj）;
        if t     end 
  end
end

 %输出距离矩阵
fid=fopen（‘distance.txt‘‘w‘）;
  for i=1:n
      for j=1:n
          fprintf（fid‘%4d ‘B（ij））;
      end
      fprintf（fid‘\n‘）;
  end
  fclose（fid）;


%问题1计算
d=zeros（n1）;

d1=zeros（n1）;
d2=zeros（n1）;


  
      %计算各点k到把（ij）点作为选定点时的距离
      mins1=1000000000000000;
    for i=1:n-1
       for j=i+1:n