%形成初始的最小生成樹
%生成初始網絡節點框架
xx=[1.9761.0560.4801.9280.1963.6400.5242.8760.1841.0080.6643.3600.5480.9163.4242.8562.4883.2722.8763.1122.3482.1283.3003.4402.3041.1722.3882.9443.616];
yy=[1.0901.0261.3041.7981.0760.4740.9141.8081.6021.5861.8220.9040.4300.1821.1920.1820.2721.7381.5601.3940.1120.3340.4741.4901.5560.3540.5061.1960.718];
plot（xxyy‘ko‘）
%對節點進行編號
for?i=1:29
????temp=[‘??‘int2str（i）];%**********************int2str（i）為將整形轉換為字符串型；‘??‘為在數字前加一個空格
????text（xx（i）yy（i）temp）%進行編號***********在圖形中指定的位置（xy）上顯示字符串temp
end
hold?on
%計算所有節點之間的距離
A=zeros（2929）%********************29X29的0矩陣
for?i=1:29
????for?j=1:29
????A（ij）=sqrt（（xx（i）-xx（j））^2+（yy（i）-yy（j））^2）%生成一個29*29的支路長度的對稱陣
????end
end
A（A==0）=inf;%把a=0的變成無窮值，為了下面的權值篩選即對角線元素為Inf
result=[];%存放的邊即最小生成樹的邊，形成樹矩陣