Enhanced Jaya algorithm(一种增强的Jaya算法)

2024-01-07 15:20

本文主要是介绍Enhanced Jaya algorithm(一种增强的Jaya算法),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一种增强的Jaya算法:

参考文献:`Enhanced Jaya algorithm: A simple but efficient optimization method for constrained engineering design problems.


Jaya算法:

Jaya算法是由印度学者Venkata Rao 于2016年首次提出。该算法的特点在于无算法相关的参数且只有进化策略。算法首先进行种群初始化,其数学表达式如下:
x i = l + ( u − l ) × λ x_{i} =l+\left ( u-l \right ) \times \lambda xi=l+(ul)×λ
式中, l 和 u 表示变量的取值下界和上界, λ 表示服从均匀分布的随机数。 式中,l和u表示变量的取值下界和上界,\lambda 表示服从均匀分布的随机数。 式中,lu表示变量的取值下界和上界,λ表示服从均匀分布的随机数。
进化策略公式如下:
v i = x i + λ 1 × ( x B e s t − ∣ x i ∣ ) − λ 2 × ( x W o r s t − ∣ x i ∣ ) , i = 1 , 2 , ⋯ , N v_{i} =x_{i} +\lambda _{1} \times \left ( x_{Best}-\left | x_{i} \right | \right ) -\lambda _{2}\times \left ( x_{Worst}-\left | x_{i} \right | \right ),i=1,2,\cdots ,N vi=xi+λ1×(xBestxi)λ2×(xWorstxi),i=1,2,,N
式中, λ 1 , λ 2 表示两个随机数,其取值范围为 [ 0 , 1 ] ; x B e s t 表示当前种群个体中的最优个体; x W o r s t 表示当前种群个体中的最差个体; N 表示种群个体数目 式中,\lambda {1} ,\lambda {2}表示两个随机数,其取值范围为[0,1]; x_{Best}表示当前种群个体中的最优个体;x_{Worst}表示当前种群个体中的最差个体;N表示种群个体数目 式中,λ1,λ2表示两个随机数,其取值范围为[0,1]xBest表示当前种群个体中的最优个体;xWorst表示当前种群个体中的最差个体;N表示种群个体数目

接着,利用贪婪策略,选择进入下一次迭代的种群个体,其数学表达式如下:
x i = { v i , i f f ( v i ) ≤ f ( x i ) x i i f f ( v i ) > f ( x i ) x_{i} =\left\{\begin{matrix}v_{i},if f\left ( v_{i} \right ) \le f\left ( x_{i} \right ) & \\x_{i} if f\left ( v_{i} \right ) > f\left ( x_{i} \right ) & \end{matrix}\right. xi={vi,iff(vi)f(xi)xiiff(vi)>f(xi)


增强的Jaya算法(EJaya):

EJaya算法细化Jaya算法的进化策略,EJaya算法包含一个全局搜索策略和局部搜索策略。
(1)局部搜索策略:
通过引入种群个体的平均解,结合种群个体中的最优个体和最差个体分别计算局部上吸引点和局部下吸引点,并以此为基础提出一种新的种群局部进化策略。
种群个体平均解: M = 1 N ∑ i = 1 N x i , i = 1 , 2 , ⋯ , N 种群个体平均解:M=\frac{1}{N} \sum_{i=1}^{N} x_{i} ,i=1,2,\cdots ,N 种群个体平均解:M=N1i=1Nxi,i=1,2,,N
局部上吸引点: P u = λ 3 × x B e s t + ( 1 − λ 3 ) × M 局部上吸引点:P_{u} =\lambda _{3}\times x_{Best} +\left ( 1-\lambda _{3} \right ) \times M 局部上吸引点:Pu=λ3×xBest+(1λ3)×M
局部下吸引点: P l = λ 4 × x W o r s t + ( 1 − λ 4 ) × M 局部下吸引点:P_{l} =\lambda _{4}\times x_{Worst} +\left ( 1-\lambda _{4} \right ) \times M 局部下吸引点:Pl=λ4×xWorst+(1λ4)×M
新的局部进化策略: v i = x i + λ 5 × ( P u − x i ) − λ 6 × ( P l − x i ) , i = 1 , 2 , ⋯ , N 新的局部进化策略:v_{i} =x_{i} +\lambda _{5} \times \left ( P_{u}-x_{i} \right ) -\lambda _{6}\times \left ( P_{l}-x_{i} \right ),i=1,2,\cdots ,N 新的局部进化策略:vi=xi+λ5×(Puxi)λ6×(Plxi),i=1,2,,N
(2)全局进化策略:
引入历史种群个体来扩大搜索空间,提高算法的全局搜索性能。历史个体生成的数学表示如下:
X o l d = { X , i f P s w i t h ≤ 0.5 X o l d , i f P s w i t h > 0.5 X_{old} =\left\{\begin{matrix}X,ifP_{swith} \le 0.5 & \\X_{old} ,ifP_{swith} > 0.5 & \end{matrix}\right. Xold={X,ifPswith0.5Xold,ifPswith>0.5
X o l d = p e r m u t i n g ( X o l d ) X_{old}=permuting\left ( X_{old}\right ) Xold=permuting(Xold)
式中, P s w i t h 表示一个 [ 0 , 1 ] 内的随机数 式中,P_{swith}表示一个[0,1]内的随机数 式中,Pswith表示一个[0,1]内的随机数

EJaya算法的流程图

在这里插入图片描述

function [BestValue,XTarget,BestCost]=EJAYA(fobj,nPop,nVar,VarMin,VarMax,MaxIt)
% 参考文献:Enhanced Jaya algorithm: A simple but efficient optimization method for constrained engineering design problems
%%输入参数
%%fhd----------------目标方程
%%nPop---------------种群个体数目 
%%nVar---------------变量数目
%%VarMin-------------变量取值的下界
%%VarMax-------------变量取值的上界
%%MaxIt--------------算法的最大迭代次数
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%输出参数
%%BestCost-----------每次迭代的最优适应度值(用于绘制迭代收敛曲线)
%%BestValue----------最优的适应度值
%%XTarget------------最优解(个体位置)
for i=1:nPopX(i,:)=VarMin+rand(1,nVar).*(VarMax-VarMin); %种群初始化f(i) = fobj(X(i,:));
end
gen=1;% 初始化算法迭代次数
[BestCost(1),ind]=min(f);
XTarget=X(ind,:);
old=X;% 初始化历史种群
% 主要迭代步骤
while(gen+1 <= MaxIt) [row,col]=size(X);[~,tindex]=min(f);Best=X(tindex,:); [~,windex]=max(f);worst=X(windex,:);xnew=zeros(row,col);fi1=rand;go1=1-fi1;fi2=rand;go2=1-fi2;ULP=(go1*Best+fi1*mean(X)-(X(i,:))); %Eq.(4)DLP=(go2*worst+fi2*mean(X)-(X(i,:))); %Eq(6)gl=rand;if gl<0.5old=X; %Eq.(8)elseold=old(randperm(nPop),:); %Eq.(9) endfor i=1:row fi=rand;if fi<0.5xnew(i,:)=(X(i,:))+rand(1,nVar).*ULP-rand(1,nVar).*DLP;%Eq.(7)elsexnew(i,:)=(X(i,:))+randn.*(old(i,:)-(X(i,:))) ; %Eq.(10)endendfor i=1:rowxnew(i,:) = max(min(xnew(i,:),VarMax),VarMin);% 边界检查fnew(i,:) = fobj(xnew(i,:));endfor i=1:nPopif(fnew(i)<f(i))X(i,:) = xnew(i,:);f(i) = fnew(i);endendgen = gen+1;[BestCost(gen),ind]=min(f);XTarget=X(ind);
end
BestValue=min(f);
end

这篇关于Enhanced Jaya algorithm(一种增强的Jaya算法)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/580354

相关文章

不懂推荐算法也能设计推荐系统

本文以商业化应用推荐为例,告诉我们不懂推荐算法的产品,也能从产品侧出发, 设计出一款不错的推荐系统。 相信很多新手产品,看到算法二字,多是懵圈的。 什么排序算法、最短路径等都是相对传统的算法(注:传统是指科班出身的产品都会接触过)。但对于推荐算法,多数产品对着网上搜到的资源,都会无从下手。特别当某些推荐算法 和 “AI”扯上关系后,更是加大了理解的难度。 但,不了解推荐算法,就无法做推荐系

康拓展开(hash算法中会用到)

康拓展开是一个全排列到一个自然数的双射(也就是某个全排列与某个自然数一一对应) 公式: X=a[n]*(n-1)!+a[n-1]*(n-2)!+...+a[i]*(i-1)!+...+a[1]*0! 其中,a[i]为整数,并且0<=a[i]<i,1<=i<=n。(a[i]在不同应用中的含义不同); 典型应用: 计算当前排列在所有由小到大全排列中的顺序,也就是说求当前排列是第

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

综合安防管理平台LntonAIServer视频监控汇聚抖动检测算法优势

LntonAIServer视频质量诊断功能中的抖动检测是一个专门针对视频稳定性进行分析的功能。抖动通常是指视频帧之间的不必要运动,这种运动可能是由于摄像机的移动、传输中的错误或编解码问题导致的。抖动检测对于确保视频内容的平滑性和观看体验至关重要。 优势 1. 提高图像质量 - 清晰度提升:减少抖动,提高图像的清晰度和细节表现力,使得监控画面更加真实可信。 - 细节增强:在低光条件下,抖

【数据结构】——原来排序算法搞懂这些就行,轻松拿捏

前言:快速排序的实现最重要的是找基准值,下面让我们来了解如何实现找基准值 基准值的注释:在快排的过程中,每一次我们要取一个元素作为枢纽值,以这个数字来将序列划分为两部分。 在此我们采用三数取中法,也就是取左端、中间、右端三个数,然后进行排序,将中间数作为枢纽值。 快速排序实现主框架: //快速排序 void QuickSort(int* arr, int left, int rig

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

秋招最新大模型算法面试,熬夜都要肝完它

💥大家在面试大模型LLM这个板块的时候,不知道面试完会不会复盘、总结,做笔记的习惯,这份大模型算法岗面试八股笔记也帮助不少人拿到过offer ✨对于面试大模型算法工程师会有一定的帮助,都附有完整答案,熬夜也要看完,祝大家一臂之力 这份《大模型算法工程师面试题》已经上传CSDN,还有完整版的大模型 AI 学习资料,朋友们如果需要可以微信扫描下方CSDN官方认证二维码免费领取【保证100%免费

dp算法练习题【8】

不同二叉搜索树 96. 不同的二叉搜索树 给你一个整数 n ,求恰由 n 个节点组成且节点值从 1 到 n 互不相同的 二叉搜索树 有多少种?返回满足题意的二叉搜索树的种数。 示例 1: 输入:n = 3输出:5 示例 2: 输入:n = 1输出:1 class Solution {public int numTrees(int n) {int[] dp = new int

Codeforces Round #240 (Div. 2) E分治算法探究1

Codeforces Round #240 (Div. 2) E  http://codeforces.com/contest/415/problem/E 2^n个数,每次操作将其分成2^q份,对于每一份内部的数进行翻转(逆序),每次操作完后输出操作后新序列的逆序对数。 图一:  划分子问题。 图二: 分而治之,=>  合并 。 图三: 回溯:

最大公因数:欧几里得算法

简述         求两个数字 m和n 的最大公因数,假设r是m%n的余数,只要n不等于0,就一直执行 m=n,n=r 举例 以18和12为例 m n r18 % 12 = 612 % 6 = 06 0所以最大公因数为:6 代码实现 #include<iostream>using namespace std;/