这段代码中 p(i) = vpv(i) * ipv(i) / dc(i);这行代码出现报错，报错内容为索引超出数组维度。索引值 2 超出数组 'vpv' 的有效范围 [1-1]。请修改并给我一段正确的代码

根据错误提示，可以看出 `vpv` 只有一个元素，而你在第 2 个元素上进行了索引，导致出现了错误。所以在这里需要进行检查 `vpv` 和 `ipv` 数组的维度，确保它们的维度是一致的。

下面给出修改后的代码：

function [duty,iterations] = PSOMPPT(vpv,ipv)
% 初始化
persistent p dc dbest counter iteration iter_max num gbest gbest_idx particles velocity;

if isempty(num)
    num=10;
endif
if isempty(p)
    p=zeros(1,num);
    dbest=0;
    counter=0;
    iteration=0;
    iter_max=15;
    gbest=0;
    gbest_idx=0;
endif
if isempty(dc)
    dc=linspace(0,0.7,num);
endif
if isempty(particles)
    particles = rand(num, 3);
    particles(:, 1) = particles(:, 1) .* (1 - 0) + 0; % 初始位置随机生成
    particles(:, 2) = zeros(num, 1); % 初始速度为 0
    particles(:, 3) = inf(num, 1); % 初始个体最佳适应度设为正无穷
endif
if isempty(velocity)
    velocity = rand(num, 1);
    velocity = velocity .* (1 - 0) + 0; % 初始速度随机生成
endif

% 计算适应度值
for i = 1:num
    p(i) = vpv(i) * ipv(i) / dc(i);
    if p(i) > gbest % 更新全局最优适应度
        gbest = p(i);
        gbest_idx = i;
    end
    if p(i) < particles(i, 3) % 更新个体最优适应度
        particles(i, 3) = p(i);
        particles(i, 1) = dc(i);
    end
end

% 更新速度和位置
for i = 1:num
    velocity(i) = velocity(i) + rand() * (particles(i, 1) - dc(i)) + rand() * (dc(gbest_idx) - dc(i));
    dc(i) = dc(i) + velocity(i);
    if dc(i) > 1 % 约束条件
        dc(i) = 1;
    end
    if dc(i) < 0 % 约束条件
        dc(i) = 0;
    end
end

% 输出结果
duty = dc(1);
counter = counter + 1;
if counter > 100 % 重置控制周期
    counter = 1;
    iteration = iteration + 1;
end
iterations = iteration;
if iteration >= iter_max % 达到迭代次数上限
    duty = dc(gbest_idx);
    return;
end

需要注意的是，这段代码需要传入两个数组 `vpv` 和 `ipv`，并且这两个数组的维度应该相同，否则也会出现类似的错误。