通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵
: C, l1 ~' g+ c# ~9 f/ b - result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值& ^8 @6 C! Z; g$ O$ q
- p=1;tb=2:length(a);
* h$ d9 A\" m3 e$ H - while size(result,2)~=length(a)-1
5 O7 k: |; ^% }3 \+ i - temp=a(p,tb);temp=temp(:);\" e# [( Y3 |0 n- Q! V
- d=min(temp); W\" _- A$ E; M$ Q2 g1 R, \
- [jb,kb]=find(a(p,tb) == d,1);
3 R3 Y! ]/ P5 z - j=p(jb);k=tb(kb);
0 U) D: a- Y; D+ j - result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];4 k% O5 G4 x/ \ J+ \
- end
\" B4 e! Z0 i5 K$ r
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示
2 m8 k" G! b# @9 E, Iresult第一行表示起点,第二行表示终点,第三行表示权值
/ ]- |6 A5 y8 a1 S; P例如# v! l) g4 B. ~! N) ?
result = - 1 2 5 4 4 7- G1 G9 L. k* ^. Y
- 2 5 4 6 7 3
- ]: e/ _7 [/ o6 Z - 50 40 10 30 42 45# P% Q F. r8 E& h( @. Z' |7 B
- 起点就是1,终点是2距离50
7 g+ S! K5 l/ K - 然后从2,到5距离40* u% e9 V+ b! J
- 从5到4距离10& s: C( i& X- R) @0 ^- C
- 从4到6距离30+ P5 D& U* @0 Q/ t
- 从4到7距离42% y% z7 u# F0 N' v\" R
- 从7到3距离45! {0 A! L: e' K! r& _. T9 C/ q; K
- 最后生成完一整棵树
6 g3 }, X0 `0 v, J/ C5 l7 S
复制代码
! L' h0 h2 v% D% g) m! [' S( Y |