通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵# x# p; l! y) }; x7 ?' } R3 A, P
- result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值
5 A$ ]& c# u/ T% a$ k Q/ r- w3 k: l1 W- [ - p=1;tb=2:length(a);4 q+ t' U8 z4 @, g5 Y) ~
- while size(result,2)~=length(a)-1 g1 J6 H2 {0 C2 n6 Q1 N
- temp=a(p,tb);temp=temp(:);
7 N4 ~0 V4 h$ x1 X- | - d=min(temp);
' k7 e* I# s7 M9 C, u- G1 { - [jb,kb]=find(a(p,tb) == d,1);
4 U3 {: A# t7 P! |* k, D - j=p(jb);k=tb(kb);2 c0 Y. ?4 e( R( s/ ~; o
- result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];
4 B& \$ c\" s) ^7 S5 W - end8 K- H& z& E% U! h
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示
6 j0 U6 r* b5 R' v0 `5 [2 nresult第一行表示起点,第二行表示终点,第三行表示权值* P! E1 H* E2 p* x, Q* f
例如' I$ x' X! v! A3 k+ j1 N: C
result = - 1 2 5 4 4 73 k2 H\" q8 _% ?& [$ A
- 2 5 4 6 7 3
- W* r! Z\" R1 l5 z% j - 50 40 10 30 42 45 o: c; H\" o# F- S: j
- 起点就是1,终点是2距离50
( Y; Z+ T+ b$ m4 U. y\" H% z9 C - 然后从2,到5距离40
T2 u; z5 L\" j r- L3 f1 Y2 ^* E - 从5到4距离10
, N\" k# `9 \; c7 V: w( e+ W - 从4到6距离30
7 M7 O/ K3 m+ [0 y2 S& Q1 [% K% w - 从4到7距离42
. B a\" `( R8 b: [ - 从7到3距离45( I% h9 Z& T+ b+ `, D% p! T6 S: i# n
- 最后生成完一整棵树
+ z! [: _, [# J5 E+ K9 N# Q/ l& s
复制代码
; K3 o. ?+ n% P5 Q) ?% q. k; M! } |