这段代码中 gbest = pbest(cur_bestid, :)显示报错。请修改后再发送一遍

好的，以下是修改后的代码：

function [duty, iterations] = PSOMPPT(vpv, ipv)
    persistent p u dc dbest counter iteration iter_max num;
    if isempty(num)
        num = 10;
    end
    if isempty(p)
        p = zeros(1, num);
        dbest = 0;
        counter = 0;
        u = 1;
        iteration = 0;
        iter_max = 15;
    end
    if isempty(dc)
        dc = linspace(0, 0.7, num);
    end
    iterations = iteration;
    if iterations <= iter_max
        if (counter >= 1 && counter <= 100)
            duty = dc(u);
            counter = counter + 1;
            return;
        end
        if (u >= 1 && u <= num)
            p(u) = vpv * ipv;
        end
        u = u + 1;
        if (u < num + 1)
            duty = dc(u);
            counter = 1;
            return;
        end
        u = 1;
        counter = 1;
        iteration = iteration + 1;
        w = 0.729;
        c1 = 1.494;
        c2 = 1.494;
        dim = num;
        swarm_size = 50;
        max_iter = 100;
        min_bound = zeros(1, dim);
        max_bound = ones(1, dim);
        x = repmat(min_bound, swarm_size, 1) + rand(swarm_size, dim) .* (repmat(max_bound, swarm_size, 1) - repmat(min_bound, swarm_size, 1));
        v = rand(swarm_size, dim);
        pbest = x;
        for i = 1:swarm_size
            if p(i) > pbest(i)
                pbest(i) = p(i);
            end
        end
        [gbestval, gbestid] = max(pbest);
        gbest = repmat(min_bound, 1, dim) + rand(1, dim) .* (repmat(max_bound, 1, dim) - repmat(min_bound, 1, dim));
        for iter = 1:max_iter
            for i = 1:swarm_size
                v(i, :) = w * v(i, :) + c1 * rand(1, dim) .* (pbest(i, :) - x(i, :)) + c2 * rand(1, dim) .* (gbest - x(i, :));
                x(i, :) = x(i, :) + v(i, :);
                for j = 1:dim
                    if x(i, j) > max_bound(j)
                        x(i, j) = max_bound(j);
                    elseif x(i, j) < min_bound(j)
                        x(i, j) = min_bound(j);
                    end
                end
                p(i) = vpv * ipv * x(i, u);
                if p(i) > pbest(i)
                    pbest(i, :) = x(i, :);
                end
            end
            [cur_bestval, cur_bestid] = max(pbest);
            if cur_bestval > gbestval
                gbestval = cur_bestval;
                gbest = pbest(cur_bestid, :);
            end
        end
        dbest = gbest(u);
        dc1 = EPOUpdateDuty(dbest, dc, iteration, iter_max, num);
        dc = dc1;
        duty = dc(u);
        return;
    else
        duty = dbest;
        return;
    end
end

function D = EPOUpdateDuty(dbest, d, iter, iter_max, num)
    D = zeros(1, num);
    dup = zeros(1, num);
    persistent s;
    if isempty(s)
        s = 0.5;
    end
    res = 0.01;
    if iter > iter_max
        iter = iter_max;
    end
    eta = (res / s) ^ (iter / iter_max);
    s = s * eta;
    for i = 1:num
        deltaD = s * (2 * rand() - 1);
        if d(i) == dbest
            dup(i) = dbest;
        else
            dup(i) = dbest + deltaD;
        end
        if dup(i) > 1
            dup(i) = 1;
        end
        if dup(i) < 0
            dup(i) = 0;
        end
    end
    D = dup;
end

请注意，修改的地方只有一行，即将 `gbest = pbest(cur_bestid, :)` 修改为 `gbest = pbest(cur_bestid, :);`，其他部分没有变化。