Matlab的图论工具箱 5——电缆最短问题

森之张卫东

1. <p>clc, clear
   
2. x=[0     5    16    20    33    23    35    25    10];
   
3. y=[15    20    24    20    25    11     7     0     3];
   
4. xy=[x;y];
   
5. d=mandist(xy); %求xy的两两列向量间的绝对值距离
   
6. d=tril(d); %截取matlab工具箱要求的下三角矩阵
   
7. b=sparse(d) %转化为稀疏矩阵
   
8. [ST,pred]=graphminspantree(b,'Method','Kruskal')  %调用最小生成树的命令
   
9. st=full(ST); %把最小生成树的稀疏矩阵转化成普通矩阵
   
10. TreeLength=sum(sum(st))  %求最小生成树的长度
    
11. view(biograph(ST,[],'ShowArrows','off','ShowWeights','on')) %画出最小生成树</p><p>%{
    
12. h = view(biograph(b,[],'ShowArrows','off','ShowWeights','on'));
    
13. set( h.Nodes(path),'Color',[1 0.4 0.4] )
    
14. fowEdges=getedgesbynodeid( h, get( h.Nodes(path),'ID' ) );
    
15. revEdges=getedgesbynodeid( h, get( h.Nodes(fliplr(path)),'ID' ) );
    
16. edges=[fowEdges;revEdges];
    
17. set(edges,'LineColor',[1 0 0]); %为了将来打印清楚，边画成黑色
    
18. set(edges,'LineWidth',1.5);  %线型宽度设置为1.5
    
19. %}</p><p><p>
    
20. </p>

复制代码