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

这个报错可能是因为输入的 `vpv` 和 `ipv` 数组长度不匹配，导致在计算 `p(i)` 时索引超出了数组维度。你可以检查一下输入的 `vpv` 和 `ipv` 数组长度是否一致，或者在使用这两个数组时进行长度的判断。例如：

for i = 1:num
    if i > length(vpv) || i > length(ipv)
        break;
    end
    p(i) = vpv(i) * ipv(i) / dc(i);
    % ...
end

这样就可以避免索引超出数组维度的问题。