【Matlab】智能优化算法_麻雀搜索算法SSA

2023-07-13

【Matlab】智能优化算法_麻雀搜索算法SSA

1.背景介绍
2.数学模型
3.文件结构
4.伪代码
5.详细代码及注释
6.运行结果
7.参考文献

1.背景介绍

麻雀通常是群居的鸟类，有很多种类。它们分布在世界的大部分地区，喜欢生活在人类生活的地方。此外，它们是杂食性鸟类，主要以谷物或杂草的种子为食。众所周知，麻雀是常见的留鸟。与其他许多小鸟相比，麻雀有很强的智慧，有很强的记忆力。请注意，有两种不同类型的圈养家雀，都是生产者和寻食者。生产者积极寻找食物来源，而蹭食者则通过生产者获得食物。此外，证据显示，鸟类通常灵活地使用行为策略，在生产和蹭食之间切换。也可以说，为了寻找他们的食物，麻雀通常使用生产者和搜寻者的策略。

研究表明，个体监视群体中其他人的行为。同时，鸟群中的攻击者，想提高自己的捕食率，就会用高摄入量的同伴来竞争食物资源。此外，当麻雀选择不同的觅食策略时，个体的能量储备可能起到重要的作用，能量储备低的麻雀会更多地蹭吃蹭喝。值得一提的是，位于种群外围的鸟类，更有可能受到捕食者的攻击，并不断尝试获得更好的位置。请注意，位于中心的动物，可能会向它们的邻居靠近，以减少它们的危险领域。我们还知道，所有的麻雀都表现出对一切事物好奇的自然本能，同时它们也总是保持着警惕。例如，当一只鸟确实检测到捕食者时，一个或多个个体发出一声鸣叫，整个群体就会飞走。

2.数学模型

麻雀搜索算法。为了简单起见，我们将麻雀的以下行为理想化，并制定了相应的规则。

生产者通常有高水平的能量储备，并为所有蹭吃蹭喝的人提供觅食区或方向。它负责识别可以找到丰富食物来源的区域。能量储备的水平取决于对个体的健身价值的评估。
一旦麻雀发现捕食者，个体就开始鸣叫作为报警信号。当报警值大于安全阈值时，生产者需要带领所有的蹭食者到安全区域。
每只麻雀只要寻找更好的食物来源，就能成为生产者，但整个种群中生产者和蹭食者的比例是不变的。
具有较高能量的麻雀会被当作生产者。几个饥饿的蹭食者更有可能飞到其他地方寻找食物，以获得更多的能量。
蹭食者跟随能提供最佳食物的生产者去寻找食物。同时，一些蹭食者可能会不断地监视生产者并争夺食物，以增加自己的捕食率。
当意识到有危险时，处于群体边缘的麻雀迅速向安全区域移动，以获得更好的位置，而处于群体中间的麻雀则随机行走，以接近其他人。

在模拟实验中，我们需要使用虚拟麻雀来寻找食物。麻雀的位置可以用以下矩阵表示：

在这里插入图片描述

其中n是麻雀的数量，d表示要优化的变量的维度。然后，所有麻雀的适应度值可以用以下向量表示：

在这里插入图片描述

其中n表示麻雀的数量，FX中每行的值表示个体的适应度值。在SSA中，具有更好适应度值的生产者在搜索过程中优先获得食物。此外，因为生产者负责搜索食物和引导整个人口的流动。所以，生产者可以在比拾荒者更广泛的地方寻找食物。根据规则（1）和（2），在每次迭代过程中，生产者的位置更新如下：

在这里插入图片描述

其中t表示当前迭代，j＝1，2。。。，d.Xt i，j表示迭代时第i个麻雀的第j维的值。itermax是迭代次数最多的常数。α∈（0,1]是一个随机数。R2（R2∈[0,1]）和ST（ST∈[0.5,1.0]）分别表示报警值和安全阈值。Q是服从正态分布的随机数。L表示amatrixof1×d，其中内部的每个元素都是1。

当R2<ST时，这意味着周围没有捕食者，生产者进入广泛搜索模式。如果R2≥ST，则意味着一些麻雀已经发现了捕食者，所有麻雀都需要迅速飞往其他安全区域。

至于拾荒者，他们需要执行规则（4）和（5）。如上所述，一些拾荒者更频繁地监视生产者。一旦他们发现生产者找到了好的食物，他们就会立即离开目前的位置去争夺食物。如果他们赢了，他们可以立即得到生产者的食物，否则他们将继续执行规则（5）。搜寻者的位置更新公式如下所述：

在这里插入图片描述

其中XP是生产者占据的最佳位置。Xworst表示当前全球最差位置。A代表一个1×d的矩阵，里面的每个元素都被随机分配为1或-1，A+=AT（AAT ）-1。当i>n/2时，表明健身值较差的第1个蹭饭的人最有可能挨饿。

在模拟实验中，我们假设这些意识到危险的麻雀占总种群的10%到20%。这些麻雀的初始位置是在种群中随机产生的。根据规则（6），数学模型可以表达如下：

在这里插入图片描述

其中Xbest是当前的全局最优位置。β,作为步长控制参数，是一个均值为0，方差为1的正态分布的随机数。 K∈[-1, 1]是一个随机数。这里fi是当前麻雀的适配值。fg和fw分别是当前全球最佳和最差的适配值。ε是最小的常数，以避免零除法错误。

为了简单起见，当fi > fg时，表示麻雀处于群体的边缘。Xbest代表种群中心的位置，在其周围是安全的。fi = fg表明处于种群中间的麻雀意识到了危险，需要向其他麻雀靠近。K表示麻雀移动的方向，也是步长的控制系数。

3.文件结构

在这里插入图片描述

Get_Functions_details.m			% 基准的全部信息和实现
main.m							% 主函数
SSA.m							% 麻雀搜索算法

4.伪代码

在这里插入图片描述

5.详细代码及注释

5.1 Get_Functions_details.m

function [lb,ub,dim,fobj] = Get_Functions_details(F)
switch F
    case 'F1'
        fobj = @F1;
        lb=-100;
        ub=100;
        dim=30;
        
    case 'F2'
        fobj = @F2;
        lb=-10;
        ub=10;
        dim=30;
        
    case 'F3'
        fobj = @F3;
        lb=-100;
        ub=100;
        dim=30;
        
    case 'F4'
        fobj = @F4;
        lb=-100;
        ub=100;
        dim=30;
        
    case 'F5'
        fobj = @F5;
        lb=-30;
        ub=30;
        dim=30;
        
    case 'F6'
        fobj = @F6;
        lb=-100;
        ub=100;
        dim=30;
        
    case 'F7'
        fobj = @F7;
        lb=-1.28;
        ub=1.28;
        dim=30;
        
    case 'F8'
        fobj = @F8;
        lb=-500;
        ub=500;
        dim=30;
        
    case 'F9'
        fobj = @F9;
        lb=-5.12;
        ub=5.12;
        dim=30;
        
    case 'F10'
        fobj = @F10;
        lb=-32;
        ub=32;
        dim=30;
        
    case 'F11'
        fobj = @F11;
        lb=-600;
        ub=600;
        dim=30;
        
    case 'F12'
        fobj = @F12;
        lb=-50;
        ub=50;
        dim=30;
        
    case 'F13'
        fobj = @F13;
        lb=-50;
        ub=50;
        dim=30;
        
    case 'F14'
        fobj = @F14;
        lb=-65.536;
        ub=65.536;
        dim=2;
        
    case 'F15'
        fobj = @F15;
        lb=-5;
        ub=5;
        dim=4;
        
    case 'F16'
        fobj = @F16;
        lb=-5;
        ub=5;
        dim=2;
        
    case 'F17'
        fobj = @F17;
        lb=[-5,0];
        ub=[10,15];
        dim=2;
        
    case 'F18'
        fobj = @F18;
        lb=-2;
        ub=2;
        dim=2;
        
    case 'F19'
        fobj = @F19;
        lb=0;
        ub=1;
        dim=3;
        
    case 'F20'
        fobj = @F20;
        lb=0;
        ub=1;
        dim=6;     
        
    case 'F21'
        fobj = @F21;
        lb=0;
        ub=10;
        dim=4;    
        
    case 'F22'
        fobj = @F22;
        lb=0;
        ub=10;
        dim=4;    
        
    case 'F23'
        fobj = @F23;
        lb=0;
        ub=10;
        dim=4;            
end

end

% F1

function o = F1(x)
o=sum(x.^2);
end

% F2

function o = F2(x)
o=sum(abs(x))+prod(abs(x));
end

% F3

function o = F3(x)
dim=size(x,2);
o=0;
for i=1:dim
    o=o+sum(x(1:i))^2;
end
end

% F4

function o = F4(x)
o=max(abs(x));
end

% F5

function o = F5(x)
dim=size(x,2);
o=sum(100*(x(2:dim)-(x(1:dim-1).^2)).^2+(x(1:dim-1)-1).^2);
end

% F6

function o = F6(x)
o=sum(abs((x+.5)).^2);
end

% F7

function o = F7(x)
dim=size(x,2);
o=sum([1:dim].*(x.^4))+rand;
end

% F8

function o = F8(x)
o=sum(-x.*sin(sqrt(abs(x))));
end

% F9

function o = F9(x)
dim=size(x,2);
o=sum(x.^2-10*cos(2*pi.*x))+10*dim;
end

% F10

function o = F10(x)
dim=size(x,2);
o=-20*exp(-.2*sqrt(sum(x.^2)/dim))-exp(sum(cos(2*pi.*x))/dim)+20+exp(1);
end

% F11

function o = F11(x)
dim=size(x,2);
o=sum(x.^2)/4000-prod(cos(x./sqrt([1:dim])))+1;
end

% F12

function o = F12(x)
dim=size(x,2);
o=(pi/dim)*(10*((sin(pi*(1+(x(1)+1)/4)))^2)+sum((((x(1:dim-1)+1)./4).^2).*...
(1+10.*((sin(pi.*(1+(x(2:dim)+1)./4)))).^2))+((x(dim)+1)/4)^2)+sum(Ufun(x,10,100,4));
end

% F13

function o = F13(x)
dim=size(x,2);
o=.1*((sin(3*pi*x(1)))^2+sum((x(1:dim-1)-1).^2.*(1+(sin(3.*pi.*x(2:dim))).^2))+...
((x(dim)-1)^2)*(1+(sin(2*pi*x(dim)))^2))+sum(Ufun(x,5,100,4));
end

% F14

function o = F14(x)
aS=[-32 -16 0 16 32 -32 -16 0 16 32 -32 -16 0 16 32 -32 -16 0 16 32 -32 -16 0 16 32;,...
-32 -32 -32 -32 -32 -16 -16 -16 -16 -16 0 0 0 0 0 16 16 16 16 16 32 32 32 32 32];

for j=1:25
    bS(j)=sum((x'-aS(:,j)).^6);
end
o=(1/500+sum(1./([1:25]+bS))).^(-1);
end

% F15

function o = F15(x)
aK=[.1957 .1947 .1735 .16 .0844 .0627 .0456 .0342 .0323 .0235 .0246];
bK=[.25 .5 1 2 4 6 8 10 12 14 16];bK=1./bK;
o=sum((aK-((x(1).*(bK.^2+x(2).*bK))./(bK.^2+x(3).*bK+x(4)))).^2);
end

% F16

function o = F16(x)
o=4*(x(1)^2)-2.1*(x(1)^4)+(x(1)^6)/3+x(1)*x(2)-4*(x(2)^2)+4*(x(2)^4);
end

% F17

function o = F17(x)
o=(x(2)-(x(1)^2)*5.1/(4*(pi^2))+5/pi*x(1)-6)^2+10*(1-1/(8*pi))*cos(x(1))+10;
end

% F18

function o = F18(x)
o=(1+(x(1)+x(2)+1)^2*(19-14*x(1)+3*(x(1)^2)-14*x(2)+6*x(1)*x(2)+3*x(2)^2))*...
    (30+(2*x(1)-3*x(2))^2*(18-32*





查看全文








原文地址：https://blog.csdn.net/Gyangxixi/article/details/131631005
本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：https://www.msipo.com/article-1164.html 如若内容造成侵权/违法违规/事实不符，请联系MSIPO邮箱：3448751423@qq.com进行投诉反馈，一经查实，立即删除！







上一篇：C#实现将小数值四舍五入为整数下一篇：【单片机】MSP430单片机，msp430f5529，DHT11 温湿度检测仪，上限报警，单击双击判定，OLED




相关阅读

• 使用Yfinance和Plotly分析金融数据
• Redis 实现库存扣减
• C#实现将小数值四舍五入为整数
• gitee注册以及使用的简单教程
• css 3d 常用函数
• vue init webpack project 报错处理
• 十、SkyWalking链路追踪
• 二.ethercat 开发总结之ESI XML查看以及logo提取
• pycharm 配置github
• ASP.NET的页面和控件
• 第一百天学习记录：线性代数：行列式（宋浩板书）
• 算法学习 day23
• 机器学习技术（五）——特征工程与模型评估
• IDEA中 application.yaml文件没有绿色的叶子
• 如何在 iframe 中获取接口数据？
• Linux+Docker+Gitlab+Jenkins+虚拟内存
• vue3 修改data后，取消恢复原值。 JSON.parse(JSON.stringify())深拷贝的使用
• 【算法与数据结构】239、LeetCode滑动窗口最大值
• VSCODE 设置同步network Error 问题
• Java实现图片与Base64编码互转























热门文章




2023年安全员-A证证模拟考试题库及安全员-A证理论考试试题
2023-11-30




1、底层世界单片机
2023-12-01




前端实现留言板
2023-12-01




构建私有Registry
2023-12-01




再探Java集合系列—LinkedList
2023-11-29




LeetCode Hot100 42.接雨水
2023-11-29




[ Vue3 ] 组合式API + Setup语法糖 如何获取组件实例？
2023-12-01




C语言--每日选择题--Day29
2023-12-01




centos7 yum安装mysql5.7
2023-12-05




【Visual Studio 新手入门指导】包括项目创建、常用快捷键、美化、项目启动、添加文件等多种基础操作，图文详细，准确无误
2023-12-01



























手机版|MSIPO技术圈
 皖ICP备19022944号-2


Copyright © 2023, msipo.com 




返回顶部