函数或变量‘pbest_fit’无法识别

非常抱歉，这是我的失误。代码中使用了变量 `pbest_fit`，但是没有进行初始化。请在代码开头添加以下语句：

pbest_fit = inf(N,1); % 个体最优适应度矩阵

这样就可以解决 `pbest_fit` 无法识别的问题。

完整代码如下：

% 果蝇优化投影寻踪算法
% 作者：AI助手

% 初始化参数
N = 30; % 果蝇数量
Max_iter = 100; % 最大迭代次数
wmax = 0.9; % 最大惯性权重
wmin = 0.4; % 最小惯性权重
c1 = 2; % 加速系数1
c2 = 2; % 加速系数2
Lb = -50; % 变量下界
Ub = 50; % 变量上界
D = 3; % 变量维度
X = zeros(N,D); % 果蝇位置矩阵
V = zeros(N,D); % 果蝇速度矩阵
fitness = zeros(N,1); % 果蝇适应度矩阵
pbest = X; % 个体最优位置矩阵
pbest_fit = inf(N,1); % 个体最优适应度矩阵
gbest = zeros(1,D); % 全局最优位置向量
gbest_fit = inf; % 全局最优适应度
iter = 0; % 迭代次数

% 初始化果蝇位置和速度
for i = 1:N
    X(i,:) = Lb + (Ub-Lb).*rand(1,D);
    V(i,:) = rand(1,D);
end

% 迭代寻优
while iter < Max_iter
    % 计算果蝇适应度
    for i = 1:N
        fitness(i) = 1/(1 + sum((X(i,:).^2))); % 目标函数为 Rosenbrock 函数
        if fitness(i) < pbest_fit(i)
            pbest(i,:) = X(i,:);
            pbest_fit(i) = fitness(i);
        end
        if fitness(i) < gbest_fit
            gbest = X(i,:);
            gbest_fit = fitness(i);
        end
    end
    
    % 更新果蝇速度和位置
    w = wmax - iter*(wmax-wmin)/Max_iter; % 更新惯性权重
    for i = 1:N
        r1 = rand(1,D);
        r2 = rand(1,D);
        V(i,:) = w*V(i,:) + c1*r1.*(pbest(i,:) - X(i,:)) + c2*r2.*(gbest - X(i,:));
        X(i,:) = X(i,:) + V(i,:);
        % 边界处理
        X(i,:) = max(X(i,:),Lb);
        X(i,:) = min(X(i,:),Ub);
    end
    
    % 更新迭代次数
    iter = iter + 1;
end

% 输出全局最优解
disp(['Global best fitness: ' num2str(gbest_fit)]);
disp(['Global best solution: ' num2str(gbest)]);

希望这次能够帮到您！