人工水母搜索算法—matlab代码
Posted 编码雪人
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了人工水母搜索算法—matlab代码相关的知识,希望对你有一定的参考价值。
clc clear foj = @ Sphere; Lb = -100; % 搜索空间下界 Ub = 100; % 搜索空间上界 N_iter = 1000; % 最大迭代次数 n_pop = 50; % 种群个数 d = 10; % 种群维度 beta = 3; gamma = 0.1; Z = zeros(n_pop, d); % 随机生成一个d维向量 Z(1, :) = rand(1, d); % 利用logistic生成n_pop个向量 for i=2:n_pop Z(i,:) = 4.0*Z(i-1,:).*(1-Z(i-1,:)); end % 将z的各个分量载波到对应变量的取值区间 pop = zeros(n_pop, d); fitness = zeros(n_pop, 1); for i=1:n_pop pop(i,:) = Lb + (Ub - Lb)*Z(i,:); fitness(i) = foj(pop(i,:)); end mu = mean(pop, 1); figure scatter(pop(:,1), pop(:,2), \'r*\'); hold on scatter(mu(:,1), mu(:,2),\'ko\'); box on [bestScore, loc] = min(fitness); bestPop = pop(loc, :); % 当前种群中食物数目最多的位置 S = pop; % S是更新后的种群 fmin = zeros(N_iter,1); for t=1:N_iter for i=1:n_pop c = abs((1-t/N_iter)*(2*rand-1)); if c>=0.5 % 洋流位置更新公式 trend = bestPop - beta*rand*mu; % 洋流的方向,即论文中的公式.9 S(i,:) = pop(i,:) + rand(1,d).*trend; % 水母位置更新,即论文中的公式.11 else % 种群内部运动 if rand>(1-c) S(i,:) = pop(i,:) + gamma*rand(1,d).*(Ub - Lb); % 被动运动,即论文公式.12 else JK = randperm(n_pop); if foj(pop(JK(1)))>=foj(pop(JK(2))) Direction = pop(JK(2),:) - pop(JK(1),:); else Direction = pop(JK(1),:) - pop(JK(2),:); end Step = rand(1, d).*Direction; S(i,:) = pop(i,:) + Step; % 主动运动,论文公式.16 end end % 检查边界 S(i,:) = simpleBound(S(i,:), Lb, Ub); Fnew = foj(S(i,:)); % 判断是否更新相应位置和适应度值 if Fnew<fitness(i) fitness(i) = Fnew; pop(i,:) = S(i,:); end % 判断是否更新最优值 if Fnew<bestScore bestScore = Fnew; bestPop = S(i,:); end end fmin(t) = bestScore; % 命令窗口输出 disp([\'Iteration \' num2str(t) \': Best Cost = \' num2str(bestScore)]); % disp([\'Bestpop \' num2str(bestPop)]); end figure semilogy(fmin, \'r-.\')
function x=simpleBound(x,Lb,Ub) % 函数名称:越界处理函数 % param x:水母 % param Lb:变量下界 % param Ub:变量上界 d = length(x); for i=1:d if x(d) > Ub x(d) = (x(d) - Ub) + Lb; elseif x(d) < Lb x(d) = (x(d) - Lb) + Ub; end end end
function [y] = Sphere(xx) d = length(xx); sum = 0; for ii = 1:d xi = xx(ii); sum = sum + xi^2; end y = sum; end
以上是关于人工水母搜索算法—matlab代码的主要内容,如果未能解决你的问题,请参考以下文章
优化算法人工蜂鸟优化算法(AHA)含Matlab源码 1470期