数学建模社区-数学中国

标题: 智能优化算法:麻雀搜索算法-附代码 [打印本页]

作者: 杨利霞    时间: 2021-8-6 16:56
标题: 智能优化算法:麻雀搜索算法-附代码
; B! q6 U% X. I& X1 v5 H3 B; {
智能优化算法:麻雀搜索算法-附代码
) c/ K! p( m$ v0 t3 m3 o2020智能优化算法:麻雀搜索算法7 U7 M1 q9 R' n$ v' ]7 f5 H( E
文章目录0 x' J' W9 T. q9 g' {* Z$ M
2020智能优化算法:麻雀搜索算法* w; i4 z  T% y) l# Y
1.算法原理
2 I" j: k3 l2 P% ^$ t3 A7 {0 R2.算法结果
; j. S4 S' R; d3.参考文献
; j1 D: R& D* ~. }3 L4.Matlab代码
) J4 M, U1 S0 `) g* Y0 {9 R4 _% \5.Python代码
1 \# X8 C! n; q6 U+ C: w8 ]% B3 A6 ~6 K) C4 g

! a2 C: ]/ ~4 ~+ [( [+ k摘要:麻雀搜索算法(Sparrow Search Algorithm, SSA)是于2020年提出的。SSA 主要是受麻雀的觅食行为和反捕食行为的启发而提出的。该算法比较新颖,具有寻优能力强,收敛速度快的优点2 M! t- v/ j8 F5 M  s7 ~8 M
1.算法原理. Y. V; ]3 p. q
建立麻雀搜索算法的数学模型,主要规则如下所述:* ]1 X- V9 @9 g/ w( m' c0 o6 l
% O9 @! Y, h$ K- ]6 E: Z8 u* T
0 M+ O; e) [6 a0 _/ ]5 S# T
发现者通常拥有较高的能源储备并且在整个种群中负责搜索到具有丰富食物的区域,为所有的加入者提供觅食的区域和方向。在模型建立中能量储备的高低取决于麻雀个体所对应的适应度值(Fitness Value)的好坏。
4 H! N2 |/ H+ }3 E. g7 R4 m  b一旦麻雀发现了捕食者,个体开始发出鸣叫作为报警信号。当报警值大于安全值时,发现者会将加入者带到其它安全区域进行觅食。
1 J5 F/ C* P$ T+ R- L' L发现者和加入者的身份是动态变化的。只要能够寻找到更好的食物来源,每只麻雀都可以成为发现者,但是发现者和加入者所占整个种群数量的比重是不变的。也就是说,有一只麻雀变成发现者必然有另一只麻雀变成加入者。9 m2 D! Y/ c7 o9 T
加入者的能量越低,它们在整个种群中所处的觅食位置就越差。一些饥肠辘辘的加入者更有可能飞往其它地方觅食,以获得更多的能量。% M5 v4 e" o: p% a5 @( e* F
在觅食过程中,加入者总是能够搜索到提供最好食物的发现者,然后从最好的食物中获取食物或者在该发现者周围觅食。与此同时,一些加入者为了增加自己的捕食率可能会不断地监控发现者进而去争夺食物资源。5 z5 ?5 ?. [- H0 t& @
当意识到危险时,群体边缘的麻雀会迅速向安全区域移动,以获得更好的位置,位于种群中间的麻雀则会随机走动,以靠近其它麻雀。
" ^( e1 j9 l8 c# q$ B6 M在模拟实验中,我们需要使用虚拟麻雀进行食物的寻找,由n只麻雀组成的种群可表示为如下形式:! x' `9 a1 y  \7 J- d- q" g+ {
: a2 Y8 x8 F& u7 @, u4 `
2.png ​        ! E! ]7 o/ _$ p% A) D0 T
(1)( K% D: r! o, p/ t# A5 G
7 C1 m" ], t* p" r% C3 M
  f- V) O1 w4 `7 d( a3 V& x
其中,d dd 表示待优化问题变量的维数,n nn 则是麻雀的数量。那么,所有麻雀的适应度值可以表示为如下形式:
: z/ c6 u0 x$ v9 b# V6 x​        3.png 0 U' f- Z: G0 Z" ^, F; p
(2)0 T3 j8 P% M3 J6 N% }
% k# p; m7 o% N3 i2 D( G6 T8 \

* ?6 J' w1 O& K7 H其中,f 表示适应度值。
# y7 @$ A3 Q+ Q% |: A" X
* l% _- J5 Q/ Y

& E. p& b% M& W+ M在 SSA 中,具有较好适应度值的发现者在搜索过程中会优先获取食物。此外,因为发现者负责为整个麻雀种群寻找食物并为所有加入者提供觅食的方向。因此,发现者可以获得比加入者更大的觅食搜索范围。根据规则(1)和规则(2),在每次迭代的过程中,发现者的位置更新描述如下:
2 E! o* L: {) j. H8 X (3) 4.png 8 ~9 O1 a8 r* m7 g8 Q# z

: j! l& @1 [- |' C1 ^8 K2 x
" n% A7 d9 @& g8 t2 Y
其中,t tt 代表当前迭代数,j = 1 , 2 , 3 , . . . , d j =1, 2, 3, . . . , dj=1,2,3,...,d。i t e m m a x item_{max}item
6 C& x: R/ }8 `4 ~0 O4 Kmax
' c0 |3 F7 X& Q( [  [5 ?7 e​       
* t9 G% y9 ~2 o' W7 e& ~
. J. J9 @! B" o& J是一个常数,表示最大的迭代次数。X i j X_{ij}X
& D: ~( C: u, Z- I" M$ M* yij6 ?1 X2 j& O3 E- d- m0 p
​       
: {, [5 [6 V2 h$ k9 S0 A) X 表示第 i ii 个麻雀在第j jj 维中的位置信息。α ∈ ( 0 , 1 ] α∈(0, 1]α∈(0,1]是一个随机数。R 2 ( R 2 ∈ [ 0 , 1 ] ) R_2(R_2∈[0,1])R
* w+ [( d, r: E1 [2 A& v20 W4 @$ _) }/ f* d9 {0 n
​        % s3 F! s% C3 W: p& q7 l& L
(R , a* R# X. s& C% a3 A
20 q5 f. m& |# u# P
​        ' d3 q# `/ Q. c
∈[0,1])和 S T ( S T ∈ [ 0.5 , 1 ] ) ST(ST∈[0.5,1])ST(ST∈[0.5,1])分别表示预警值和安全值。Q QQ 是服从正态分布的随机数。L LL 表示一个 1 × d 1×d1×d 的矩阵,其中该矩阵内每个元素全部为 1。
9 \3 i  L0 e" _! M/ z  o! ?0 r2 T" }7 r; _/ l/ B

3 y+ w5 g" M8 M1 n, E! G  G当 R 2 < S T R2< STR2<ST 时,这意味着此时的觅食环境周围没有捕食者,发现者可以执行广泛的搜索操作。如果 R 2 ≥ S T R2≥ STR2≥ST,这表示种群中的一些麻雀已经发现了捕食者,并向种群中其它麻雀发出了警报,此时所有麻雀都需要迅速飞到其它安全的地方进行觅食。
0 f1 m% y& r7 Z( I9 T+ n% @% A/ ?2 }5 N; A, {+ E2 x

% T% [/ K6 P, F对于加入者,它们需要执行规则(3)和规则(4)。如前面所描述,在觅食过程中,一些加入者会时刻监视着发现者。一旦它们察觉到发现者已经找到了更好的食物,它们会立即离开现在的位置去争夺食物。如果它们赢了,它们可以立即获得该发现者的食物,否则需要继续执行规则(4)。加入者的位置更新描述如下:
3 l0 A: _' C) P' r9 Q$ o& j. w3 p" g 5.png
9 i6 d9 W2 [1 H5 B$ [​        * S4 G' Z" ~, K" W" n* s1 I( D" X
则表示当前全局最差的位置。A AA表示一个 1 × d 1×d1×d 的矩阵,其中每个元素随机赋值为 1 或-1,并且 A + = A T ( A A T ) − 1 A^+=A^T(AA^T)^{-1}A 9 S" k0 k  A$ X% h. g; }
+
8 [. V  \" _' K2 s( C =A " B+ ?* R+ N. A
T
3 o8 P1 F, w( N0 Y  v/ m3 u' V! T (AA " B1 U8 V$ |! t7 `  [2 w. _4 z
T
. s+ B, a8 p2 \! x" b! V ) 7 Z4 z# G9 v' ]. Z
−13 e8 ?  ?/ q* r3 y) f
。当i >n/2 时,这表明,适应度值较低的第 i 个加入者没有获得食物,处于十分饥饿的状态,此时需要飞往其它地方觅食,以获得更多的能量。1 [4 z. z& e5 I0 c* ^! h8 z0 U

9 X: n# g3 b+ p% E3 I3 [
( S( x( ~) Q+ W$ x' Z1 P
在模拟实验中,我们假设这些意识到危险的麻雀占总数量的 10% 到 20%。这些麻雀的初始位置是在种群中随机产生的。根据规则(5),其数学表达式可以表示为如下形式:
5 m/ x/ A* b% W4 ^; c4 Q 6.png
. @% |  b0 A% g" X (5)
3 I& I* b7 |2 b, y% q9 H; F4 @8 K/ O. k+ v, \8 v8 o
$ c  }1 C& V* A. g) ^* }! n
其中,其中 X b e s t X_{best}X , Z" ^' x, x9 G, D; o0 v0 {# h
best
4 `/ S3 n1 o, m! a2 Z/ X2 }+ Q​       
9 D8 _0 P& o5 J) c. a( h  E 是当前的全局最优位置。β ββ 作为步长控制参数,是服从均值为 0,方差为 1 的正态分布的随机数。K ∈ [ − 1 , 1 ] K∈[-1,1]K∈[−1,1]是一个随机数,fi则是当前麻雀个体的适应度值。f g f_gf 1 O4 S5 p0 T+ P' `2 z4 P
g
: B7 ?; O, R/ u0 l9 D" Z​       
2 J1 A4 `; Z4 h9 U' v( @ 和 f w f_wf 3 i/ W$ P" N: p$ M; c
w
8 k. V- S. Z' m$ R# F5 t​        ! r5 G# ~; W7 ^% s  \& f
分别是当前全局最佳和最差的适应度值。ε \varepsilonε 的常数,以避免分母出现零。9 Y$ E; s  O9 @) ^$ ~3 S. l; G. D
, [: @9 T3 _8 V
( e2 Y* w; E8 C; d
为简单起见,当 f i > f g f_i >f_gf ; v& ]3 D' A  g# k$ a
i4 s8 r' Y% i" u
​        ; v# G0 u" B, o* u
>f ! O7 ?" v. z, h( H
g' R$ Z7 Q  S4 d6 {
​       
# W; D& r& E' i+ b: Q& K0 i 表示此时的麻雀正处于种群的边缘,极其容易受到捕食者的攻击。X b e s t X_{best}X 2 j3 j8 S+ e0 r! ^
best
; G- d: f4 W. I9 c) r2 Z! t​        # _! D+ O. K3 ?3 I" |# T4 n: h
表示这个位置的麻雀是种群中最好的位置也是十分安全的。f i = f g f_i = f_gf . W* w1 o5 i  f7 z& S0 g+ ]
i
9 J, p8 R! y7 `( U​       
3 W# y6 D" Q( F =f
# K% y! |) {$ J) a: s; ng
$ h/ T. {5 \  _3 J/ ]​       
2 s9 n' z. [  q- r& x 时,这表明处于种群中间的麻雀意识到了危险,需要靠近其它的麻雀以此尽量减少它们被捕食的风险。K KK 表示麻雀移动的方向同时也是步长控制参数。
3 g9 z5 h+ s& L! _! P+ x3 d7 ]  Y9 d
" X% ~% a% q+ u5 w
算法流程% u& ~7 g) k) |# e/ b
0 h( t4 i% L, X6 f) D

. w. r" m. Z) M  g) h& w) xStep1: 初始化种群,迭代次数,初始化捕食者和加入者比列。6 i$ v9 ]  G- C
0 V+ \! T0 g) {

) u# b4 Y( _! Z1 M% lStep2:计算适应度值,并排序。0 C& F( {9 n& K: {/ Q
3 h% ^7 f; O  g+ Y

) ^6 l8 n" ~5 F/ U6 IStep3:利用式(3)更新捕食者位置。
9 ^) c( j8 q  y" @& A
" `: L# F. f( i! S2 n5 S! L

6 v* s/ O( H4 H4 {# k) E/ UStep4:利用式(4)更新加入者位置。* I8 ^0 R5 \8 w1 M8 y
; v8 R1 |: x& k3 j5 Q( ]
: A& \/ \- K6 O$ ]( S# x+ X8 ]
Step5:利用式(5)更新警戒者位置。
; q) m) C! ?1 O( j: j+ U: y) S' B! _& P. Z
" ~6 i5 ?& n. x; b- {  S
Step6:计算适应度值并更新麻雀位置。
0 J6 e% j9 Z% ?& Z) L# a% r4 O* W$ ~, g! o
8 O% m* F2 @$ {3 o+ i6 l
Step7:是否满足停止条件,满足则退出,输出结果,否则,重复执行Step2-6;) @1 K+ x  e4 U5 o7 v( v
* R' F7 `) K+ t, {& q
. A5 Z( ]6 B( J% V! [+ F
2.算法结果+ s$ S8 B. }$ d' E5 \
1.png
" u: ^9 L) o2 `7 G5 \
) Z- ^& {' k! r# k5 v# u
' }( c7 k: V/ w6 w3 _9 P. o- p% [3 ?
3.参考文献
3 T. T6 R. E) ?/ I1 d% m& `1 P: b" `[1] Xue J , Shen B . A novel swarm intelligence optimization approach: sparrow search algorithm[J]. Systems ence & Control Engineering An Open Access Journal, 2020, 8(1):22-34.2 p3 ?* d4 D1 \! l9 f/ ?) c

- {1 u! |4 i2 }
) N2 }, P" }) @% e! Y
4.Matlab代码% x3 I# A$ ]9 U' i+ Y0 R
麻雀搜索算法" v3 c5 R1 k+ P& k0 g* ?4 @
改进算法:. ^3 l, n% b3 M
1.基于反向策略的麻雀搜索算法) _) r. {! e, J* c
2.基于Tent混沌映射的麻雀搜索算法
. o( `5 ?: z; U% w3.基于Logistic混沌映射的麻雀搜索算法
& _% k' y7 F( N4.基于Circle混沌映射的麻雀搜索算法
' V  N" Y9 ^1 D; o( q5.基于Piecewise混沌映射的麻雀搜索算法# f: E1 {0 Y9 ~4 r' H, B2 [! m
6.基于Chebyshev混沌映射的麻雀搜索算法
* `7 Q- A- h/ Q* s' A7.基于Sine混沌映射的麻雀搜索算法: e7 b# G% z/ w- t; X! w
8.基于Singer混沌映射的麻雀搜索算法1 G6 V. ]# t3 Z8 M  F
9.基于迭代混沌映射的麻雀搜索算法. L$ ^" d# C7 L& I
10.基于Sinusoidal混沌映射的麻雀搜索算法
* R& D% ]8 C! l3 @4 `& i11.基于随机游走改进的麻雀搜索算法
! d! B  |" E/ d2 c7 s& W12.基于萤火虫改进的麻雀搜索算法$ \( {) g4 X: x& y: m
13.基于精英反向策略的麻雀搜索算法1
- S2 F. O6 C: V14.基于levy飞行改进的麻雀搜索算法% N) A9 {3 U1 o, u4 h
15.基于自适应t分布的麻雀算法1 x2 T2 [5 R' d! m6 H+ ?
改进麻雀文献复现代码:
+ o& D: K/ V' N) y9 Z1.混沌麻雀。
5 O2 N3 f. O5 g+ y' ~参考文献:[1]吕鑫,慕晓冬,张钧,王震.混沌麻雀搜索优化算法[J/OL].北京航空航天大学学报:1-10[2020-11-16].https://doi.org/10.13700/j.bh.1001-5965.2020.0298.
0 r: v2 h7 ^; S( t
  p2 Y# x8 A" G: \* i
: k1 d* c; p0 a/ `
2.融合柯西变异和反向学习的改进麻雀算法
/ x# }, ~5 y% t4 L8 r3 Y/ p8 v[1]毛清华,张强.融合柯西变异和反向学习的改进麻雀算法[J/OL].计算机科学与探索:1-12[2020-12-16].http://kns.cnki.net/kcms/detail/11.5602.tp.20201203.1601.006.html.
5 V; T; z2 _0 @- V& Y- j: _/ I2 q
1 ^, m+ q0 v# o7 W+ J
* K  K2 m; i, Z
[3. 混合正弦余弦算法和Lévy飞行的麻雀算法(ISSA)]* r) Y, l* f0 T3 g6 D1 i! |
[1]毛清华,张强,毛承成,柏嘉旋.混合正弦余弦算法和Lévy飞行的麻雀算法[J/OL].山西大学学报(自然科学版):1-6[2021-04-09].https://doi.org/10.13451/j.sxu.ns.2020135.( X3 f! f" B1 p, R

3 a$ b* J% X+ z+ N) W, f% G

/ J4 \* H7 q5 ~: O7 X3 I[4.基于 Sobol 序列和纵横交叉策略的麻雀搜索算法(SSASC)]
" A9 h1 a) a+ j8 n[1]段玉先,刘昌云.基于 Sobol 序列和纵横交叉策略的麻雀搜索算法[J/OL].计算机应用. https://kns.cnki.net/kcms/detail/51.1307.TP.20210525.1453.002.html) x# G0 [5 {6 n
+ Q: {( i6 N6 u5 f' P
2 S. A$ c% z4 C& P/ d7 x( [5 h* b
5.Python代码
* j6 l" S% u: r0 B* x麻雀搜索算法; x9 A: c( `% q6 D, }; `' U
改进算法:/ M3 A- Y+ D+ K9 @( T. N
基于Sinusoidal混沌映射的麻雀搜索算法 python 代码# X0 a* X$ O; ~2 q- f7 S
基于迭代混沌映射的麻雀搜索算法 python 代码
" T4 Q/ |" {% z" {0 [基于Singer混沌映射的麻雀搜索算法 python 代码
& [' ~) r+ n4 s- k. `5 \7 ?) m基于Sine混沌映射的麻雀搜索算法 python代码
9 y' N$ ^# h. e+ R8 l( d: V基于Piecewise混沌映射的麻雀搜索算法 python代码
8 ~. F4 n) n! _! s基于Logistic混沌映射的麻雀搜索算法 python 代码
7 l8 b# C# B3 u/ F2 E; u/ `) n- @基于Circle混沌映射的麻雀搜索算法 python 代码7 i% a  d$ a; [( J  P2 T
基于Chebyshev混沌映射的麻雀搜索算法 python代码: U, W2 d; N; S; ?- u
基于Tent混沌映射的麻雀搜索算法 python代码7 k! Q1 Z4 B' L8 z; b
基于反向策略的麻雀搜索算法 python代码
4 o( X8 n# m5 x& A  c" X1 I6 N基于精英反向策略的麻雀搜索算法1 python代码0 h7 `! y9 p2 v+ f: \6 @
基于精英反向策略的麻雀搜索算法2 python 代码
$ e0 {5 \5 u4 p) F2 f基于萤火虫改进的麻雀搜索算法 python代码; L; x; h/ o/ c  i7 [
基于levy飞行改进的麻雀搜索算法 python 代码- A/ }. {6 U6 p2 E5 P! F$ S( Z
基于随机游走改进的麻雀搜索算法+ q( I) o& v: e9 z0 O8 p% A# c
基于自适应t分布的麻雀搜索算法
" o4 ^' x4 ]* `& W. n改进麻雀文献复现代码:
& \0 G& e2 ?7 d3 C' m. {1.混沌麻雀4 l; Y$ ~& B) Q" K- v, x' N9 x7 Y4 D
参考文献:[1]吕鑫,慕晓冬,张钧,王震.混沌麻雀搜索优化算法[J/OL].北京航空航天大学学报:1-10[2020-11-16].https://doi.org/10.13700/j.bh.1001-5965.2020.0298.4 O* m8 a+ W* n

: v4 [' c& U* y; g" m

: P& h( T) f2 T( `2.混合正弦余弦算法和Lévy飞行的麻雀算法(ISSA)& ]" X( m. Z' _- A1 U
[1]毛清华,张强,毛承成,柏嘉旋.混合正弦余弦算法和Lévy飞行的麻雀算法[J/OL].山西大学学报(自然科学版):1-6[2021-04-09].https://doi.org/10.13451/j.sxu.ns.2020135.; v. e/ _" Y, q2 a$ n2 Y" r
————————————————
& J0 W7 p2 l* v+ U; H, ~版权声明:本文为CSDN博主「Jack旭」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。( \1 Q; E- m, y. A
原文链接:https://blog.csdn.net/u011835903/article/details/108830958
9 r7 ^1 {; w  r" ~3 i! @( Q: X2 a  P) Q: R2 ?

8 ^  ?) i6 ~" c1 N) b$ {/ i
作者: 1051373629    时间: 2021-8-17 17:13
谢谢!
, ]* p% Z: v6 R: H




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5