使用粒子群优化算法计算Shubert函数的最优解function f = Shubertx x为一个2维向量 计算Shubert函数的值s1 = 0;s2 = 0;for i = 15 s1 = s1 + icosi+1x1+i; s2 = s2 + icosi+1x2+i;endf = prodsinrepmat1521repmatx1x215+repmat1;215202 s1

这个错误提示是因为在计算适应度时，需要把两个变量x和y传入到Shubert函数中，但是这个函数只接受一个2维向量作为输入，因此需要修改适应度计算的部分代码。

可以将适应度计算的部分代码修改为：

% 计算适应度 for i = 1:popSize fit(i) = Shubert(pop(i,:)); end

% 更新最优解和适应度 [currentBestFit, currentBestIdx] = min(fit); if currentBestFit < bestFit bestFit = currentBestFit; bestPos = pop(currentBestIdx,:); end

% 记录历史最优解和适应度 bestFitHistory(iter) = bestFit; bestPosHistory(iter,:) = bestPos;

其中，Shubert函数的输入仍然是一个2维向量，适应度计算只需要调用Shubert函数即可，不需要再传入x和y。