【单目标优化求解】混合正弦余弦算法和Lévy飞行改进麻雀算法求解单目标优化问题【含Matlab源码 1653期】

✅博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，Matlab项目合作可私信。
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭：行百里者，半于九十。

更多Matlab仿真内容点击👇
Matlab图像处理（进阶版）
路径规划（Matlab）
神经网络预测与分类（Matlab）
优化求解（Matlab）
语音处理（Matlab）
信号处理（Matlab）
车间调度（Matlab）

⛄一、麻雀算法简介（具体理论知识见参考文献[3]）

麻雀搜索算法(Sparrow Search Algorithm, SSA)是于2020年提出的。SSA 主要是受麻雀的觅食行为和反捕食行为的启发而提出的。该算法比较新颖，具有寻优能力强，收敛速度快的优点。
1 算法原理
建立麻雀搜索算法的数学模型，主要规则如下所述：
（1）发现者通常拥有较高的能源储备并且在整个种群中负责搜索到具有丰富食物的区域，为所有的加入者提供觅食的区域和方向。在模型建立中能量储备的高低取决于麻雀个体所对应的适应度值(Fitness Value)的好坏。
（2）一旦麻雀发现了捕食者，个体开始发出鸣叫作为报警信号。当报警值大于安全值时，发现者会将加入者带到其它安全区域进行觅食。
（3）发现者和加入者的身份是动态变化的。只要能够寻找到更好的食物来源，每只麻雀都可以成为发现者，但是发现者和加入者所占整个种群数量的比重是不变的。也就是说，有一只麻雀变成发现者必然有另一只麻雀变成加入者。
（4）加入者的能量越低，它们在整个种群中所处的觅食位置就越差。一些饥肠辘辘的加入者更有可能飞往其它地方觅食，以获得更多的能量。
（5）在觅食过程中，加入者总是能够搜索到提供最好食物的发现者，然后从最好的食物中获取食物或者在该发现者周围觅食。与此同时，一些加入者为了增加自己的捕食率可能会不断地监控发现者进而去争夺食物资源。
（6）当意识到危险时，群体边缘的麻雀会迅速向安全区域移动，以获得更好的位置，位于种群中间的麻雀则会随机走动，以靠近其它麻雀。
在模拟实验中，我们需要使用虚拟麻雀进行食物的寻找，由n只麻雀组成的种群可表示为如下形式：

其中，d表示待优化问题变量的维数，n则是麻雀的数量。那么，所有麻雀的适应度值可以表示为如下形式：

其中，f表示适应度值。
在SSA中， 具有较好适应度值的发现者在搜索过程中会优先获取食物。此外， 因为发现者负责为整个麻雀种群寻找食物并为所有加入者提供觅食的方向。因此，发现者可以获得比加入者更大的觅食搜索范围。根据规则(1)和规则(2)，在每次迭代的过程中，发现者的位置更新描述如下：

其中， t代表当前迭代数， j=1， 2， 3， …， d.item maz是一个常数，表示最大的迭代次数。Xy表示第i个麻雀在第j维中的位置信息。xE(0，1]是一个随机数。R2(R2E[0，1])和ST(STe[0.5， 1] ) 分别表示预警值和安全值。Q是服从正态分布的随机数.L表示一个1×d的矩阵， 其中该矩阵内每个元素全部为
1.
当R2<ST时，这意味着此时的觅食环境周围没有捕食者，发现者可以执行广泛的搜索操作。如果R2≥ST，这表示种群中的一些麻雀已经发现了捕食者，井向种群中其它麻雀发出了警报，此时所有麻雀都需要迅速飞到其它安全的地方进行觅食。对于加入者，它们需要执行规则(3)和规则(4)。如前面所描述，在觅食过程中，一些加入者会时刻监视着发现者。一旦它们察觉到发现者已经找到了更好的食物，它们会立即离开现在的位置去争夺食物。如果它们赢了，它们可以立即获得该发现者的食物，否则需要继续执行规则(4)。加入者的位置更新描述如下：

其中， X， 是目前发现者所占据的最优位置， X worst则表示当前全局最差的位置。A表示一个1×d的矩阵， 其中每个元素随机赋值为1或-1，并且A+=A(AA)-.当i>n/2时，这表明，适应度值较低的第i个加入者没有获得食物，处于十分饥饿的状态，此时需要飞往其它地方觅食，以获得更多的能量。在模拟实验中，我们假设这些意识到危险的麻雀占总数量的10%到20%。这些麻雀的初始位置是在种群中随机产生的。根据规则(5)，其数学表达式可以表示为如下形式：

其中， 其中X best是当前的全局最优位置。β作为步长控制参数， 是服从均值为0， 方差为1的正态分布的随机数。KE[-1， 1] 是一个随机数，则是当前麻雀个体的适应度值。f，和fw分别是当前全局最佳和最差的适应度值。e的常数，以避免分母出现零。
为简单起见， 当f：>f， 表示此时的麻雀正处于种群的边缘， 极其容易受到捕食者的攻击。X best表示这个位置的麻雀是种群中最好的位置也是十分安全的。f；=f，时，这表明处于种群中间的麻雀意识到了危险，需要靠近其它的麻雀以此尽量减少它们被捕食的风险。K表示麻雀移动的方向同时也是步长控制参数。

2 算法流程

Step1： 初始化种群，迭代次数，初始化捕食者和加入者比列。
Step2：计算适应度值，并排序。
Step3：利用式（3）更新捕食者位置。
Step4：利用式（4）更新加入者位置。
Step5：利用式（5）更新警戒者位置。
Step6：计算适应度值并更新麻雀位置。
Step7：是否满足停止条件，满足则退出，输出结果，否则，重复执行Step2-6；

⛄二、部分源代码

%% 清除环境变量
clear
clc

%% 参数设置
N = 30; % 种群规模
Function_name = ‘F2’; % 从F1到F23的测试函数的名称（本文中的表1、2、3）
Max_iteration = 100; % 最大迭代次数
cnt_max = 30;
% 加载所选基准函数的详细信息
[lb, ub, dim, fobj] = Get_Functions_details(Function_name);

Curve_GWO = zeros(1, Max_iteration);
Curve_PSO = zeros(1, Max_iteration);
Curve_MFO = zeros(1, Max_iteration);
Curve_SSA = zeros(1, Max_iteration);
Curve_SSA1 = zeros(1, Max_iteration);
Curve_SSA2 = zeros(1, Max_iteration);
Curve_ISSA = zeros(1, Max_iteration);

for cnt = 1:cnt_max
disp([‘第’, num2str(cnt), ‘次迭代’]);
% 初始化种群位置
X = initialization(N, dim, ub, lb);

[GWO_Best_score(cnt), GWO_Best_pos(cnt, :), GWO_Curve] = GWO(X, N, Max_iteration, lb, ub, dim, fobj);
[PSO_Best_score(cnt), PSO_Best_pos(cnt, :), PSO_Curve] = PSO(X, N, Max_iteration, lb, ub, dim, fobj);
[MFO_Best_score(cnt), MFO_Best_pos(cnt, :), MFO_Curve] = MFO(X, N, Max_iteration, lb, ub, dim, fobj);
[SSA_Best_score(cnt), SSA_Best_pos(cnt, :), SSA_Curve] = SSA(X, N, Max_iteration, lb, ub, dim, fobj);
[SSA1_Best_score(cnt), SSA1_Best_pos(cnt, :), SSA1_Curve] = SSA1(X, N, Max_iteration, lb, ub, dim, fobj);
[SSA2_Best_score(cnt), SSA2_Best_pos(cnt, :), SSA2_Curve] = SSA2(X, N, Max_iteration, lb, ub, dim, fobj);
[ISSA_Best_score(cnt), ISSA_Best_pos(cnt, :), ISSA_Curve] = ISSA(X, N, Max_iteration, lb, ub, dim, fobj);

Curve_GWO = Curve_GWO+GWO_Curve;
Curve_PSO = Curve_PSO+PSO_Curve;
Curve_MFO = Curve_MFO+MFO_Curve;
Curve_SSA = Curve_SSA+SSA_Curve;
Curve_SSA1 = Curve_SSA1+SSA1_Curve;
Curve_SSA2 = Curve_SSA2+SSA2_Curve;
Curve_ISSA = Curve_ISSA+ISSA_Curve;

end

Curve_GWO = Curve_GWO/cnt_max;
Curve_PSO = Curve_PSO/cnt_max;
Curve_MFO = Curve_MFO/cnt_max;
Curve_SSA = Curve_SSA/cnt_max;
Curve_SSA1 = Curve_SSA1/cnt_max;
Curve_SSA2 = Curve_SSA2/cnt_max;
Curve_ISSA = Curve_ISSA/cnt_max;

std_GWO = std(GWO_Best_score);
std_PSO = std(PSO_Best_score);
std_MFO = std(MFO_Best_score);
std_SSA = std(SSA_Best_score);
std_SSA1 = std(SSA1_Best_score);
std_SSA2 = std(SSA2_Best_score);
std_ISSA = std(ISSA_Best_score);

best_GWO = max(GWO_Best_score);
best_PSO = max(PSO_Best_score);
best_MFO = max(MFO_Best_score);
best_SSA = max(SSA_Best_score);
best_SSA1 = max(SSA1_Best_score);
best_SSA2 = max(SSA2_Best_score);
best_ISSA = max(ISSA_Best_score);

worst_GWO = min(GWO_Best_score);
worst_PSO = min(PSO_Best_score);
worst_MFO = min(MFO_Best_score);
worst_SSA = min(SSA_Best_score);
worst_SSA1 = min(SSA1_Best_score);
worst_SSA2 = min(SSA2_Best_score);
worst_ISSA = min(ISSA_Best_score);

mean_GWO = mean(GWO_Best_score);
mean_PSO = mean(PSO_Best_score);
mean_MFO = mean(MFO_Best_score);
mean_SSA = mean(SSA_Best_score);
mean_SSA1 = mean(SSA1_Best_score);
mean_SSA2 = mean(SSA2_Best_score);
mean_ISSA = mean(ISSA_Best_score);

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]毛清华,张强,毛承成,柏嘉旋.混合正弦余弦算法和Lévy飞行的麻雀算法[J].山西大学学报(自然科学版)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合