MFO算法MATLAB代码,【优化求解】飞蛾扑火算法（MFO）【Matlab 192期】

以下是MATLAB实现的MFO算法代码：

function [bestFit, bestInd] = MFO(func, dim, lb, ub, maxIter)

% 参数说明：
% func：优化函数
% dim：优化问题的维度
% lb：每个维度的下界
% ub：每个维度的上界
% maxIter：最大迭代次数

% MFO算法参数设置
N = 20;              % 飞蛾数量
a = 0.2;             % 吸引度系数
b = 1;               % 距离衰减系数
tmax = maxIter;      % 最大迭代次数

% 初始化飞蛾位置和适应度值
X = zeros(N,dim);
F = zeros(N,1);
for i = 1:N
    X(i,:) = lb + (ub-lb).*rand(1,dim); % 随机初始化位置
    F(i) = func(X(i,:));                % 计算适应度值
end

% 记录最佳适应度值和最佳位置
[bestFit, bestInd] = min(F);
bestX = X(bestInd,:);

% 迭代搜索
for t = 1:tmax
    % 计算飞蛾之间的距离
    D = pdist2(X,X);
    D(D==0) = Inf;
    
    % 计算每个飞蛾的吸引度
    A = zeros(N,1);
    for i = 1:N
        for j = 1:N
            A(i) = A(i) + (F(j)<F(i))*exp(-b*D(i,j));
        end
    end
    A = a*A/sum(A);
    
    % 更新飞蛾位置
    for i = 1:N
        % 计算移动方向
        dir = zeros(1,dim);
        for j = 1:N
            if j ~= i
                dir = dir + A(j)*(X(j,:)-X(i,:))/D(i,j);
            end
        end
        
        % 更新位置
        X(i,:) = X(i,:) + dir;
        
        % 边界处理
        X(i,X(i,:)<lb) = lb(X(i,:)<lb);
        X(i,X(i,:)>ub) = ub(X(i,:)>ub);
        
        % 计算适应度值
        F(i) = func(X(i,:));
        
        % 更新最佳位置和最佳适应度值
        if F(i) < bestFit
            bestFit = F(i);
            bestX = X(i,:);
        end
    end
    
    % 显示迭代信息
    disp(['Iteration ' num2str(t) ': Best Fit = ' num2str(bestFit)]);
end

% 返回最佳适应度值和最佳位置
bestFit = -bestFit; % 将最小值转换为最大值
bestInd = -1;
end

使用时，只需要传入优化函数、维度、下界、上界和最大迭代次数等参数即可，如下所示：

% 优化函数
func = @(x) sum(x.^2);

% 优化问题的维度
dim = 10;

% 每个维度的下界和上界
lb = -10*ones(1,dim);
ub = 10*ones(1,dim);

% 最大迭代次数
maxIter = 100;

% 调用MFO函数进行优化
[bestFit, bestInd] = MFO(func, dim, lb, ub, maxIter);

% 显示最佳适应度值和最佳位置
disp(['Best Fit = ' num2str(bestFit)]);
disp(['Best Ind = ' num2str(bestInd)]);

注意，这里的优化函数必须是一个能够计算出某个位置的适应度值的函数。在这里，我使用了一个简单的函数 $f(x)=\sum_{i=1}^n x_i^2$ 作为优化函数进行测试。