实际问题引入3 l& R' F' n( x* O- f
/ J9 L# x# J( {5 ?1 l/ f( v1 N- P# g实这道题的答案,其实就是找到这个图的最小生成树。( p8 W( x4 r- ]6 G# Z% H+ B
6 D5 {8 R8 R& X$ n; U" @ v
Kruskal算法
6 z( v& N2 E4 v. F7 R+ n& t, w此算法可以称为“加边法”,初始最小生成树的边数为 0,每迭代一次就选择一条满足条件的最小代价的边,加入到最小生成树边的集合里面。
+ k% j9 K' M* ?% ?! O2 ^其实核心思想就是贪心思想:通过局部最优达到整体最优
2 ?: V# |$ S. W6 y: x( D3 |. U2 p# O0 k( a8 d8 p' X
将所有的边权进行排序
" ~9 i' O/ V' p8 {2 j不断迭代选择权最小的边,直到所有的点被连起来(边数=节点数-1)。
5 X+ |! ^4 I( N7 v, t在迭代期间,如果边构成了环,就要丢弃该边,因为树中是不存在环的!
4 {$ h& Z3 s0 K: m( C( i整体代码展示0 H% P( H3 \ D7 i
在matlab中,最小生成树的生成直接用minspantree()函数就行。- s=[1,1,1,1,2,2,3,3,4,4,5,5,6];
4 H! C7 v* H\" W& |; w* C, @$ k - t=[2,3,4,5,3,6,5,7,5,6,6,7,7];$ Q& f3 k8 D% f8 ?& k+ R2 Y& H
- w=[35,24,10,25,25,20,15,11,12,30,15,25,18];! i- _1 X! Z! d! Z4 O6 W. A4 ?+ v+ \. o
- names={'1','2','3','4','5','6','7'};7 a9 j6 R- n# y
- G=graph(s,t,w,names);
) F/ o+ l2 q$ a7 \' ~6 F - p=plot(G,"EdgeLabel",G.Edges.Weight);/ T$ m; Q# k1 I/ n4 a# B9 k3 t9 c
- % 求解最小生成树9 a, c) `6 X* J; P: O
- T=minspantree(G,"Method","sparse");
, ^! B! ~: T8 n9 E - % sparse代表的是Kruskal算法 R; o2 e; u+ e9 U# z
- % dense代表的是Prim算法6 U5 m) x/ v1 z! L8 {6 _
- 0 q8 q5 p2 |2 c* W7 p4 D t. {4 S
- % sparse:Kruskal算法
# q\" j; e$ p, q3 S - % 算法按权重对所有的边排序,然后将不构成循环的边添加到树中
: |2 e/ f\" H, r: y1 Z\" a - p=plot(G,"EdgeLabel",G.Edges.Weight);
# \, U) l1 X. N\" ~. S - highlight(p,T,"NodeColor","red","EdgeColor","red");
% Y$ J: B4 A6 h8 g, Q' X - % 将最小生成树的边设置为红色!
' e5 d6 f$ j7 P; ]8 V1 a y/ }
复制代码
5 Z% I3 o. e, V( V! n: {" A: ~- S
/ l# ?2 C- E; B( A) h; _ J+ z, Y% s生成的最小生成树:
$ g3 e) l, u! P1 x
( I/ J' \/ a5 w" W# F" j
我们也可以把最小生成树的边和节点打印出来,也可以把整段路的权加起来看看:
7 f) X. ^4 O# J" g! N
尾声看到这里,相信我们已经学会Kruskal算法寻找最小生成树的过程了,当然,这离数学建模的要求,离我们的目标还非常遥远,博主在不断学习的过程中,也希望可以通过分享学习日记的方式带动大家! , }1 E' c7 r/ E/ o1 [ x% a7 v
3 w; N7 h) e/ J3 p4 v6 k5 A. n" p$ R |