纠错：nPop = 30; % 初始粒子群微粒个数 dim = 2; % 空间维数 maxIter1 = 300; % 最大迭代次数 w = 0.9; % 惯性因子 c1 = 2; % 个体学习因子 c2 = 2; % 社会学习因子 limit = [75,150;-40,40]; % 位置参数限制 vlimit = [-5, 5;-5, 5]; % 速度参数限制 for i = 1:dim x(:,i) = limit(i, 1) + (limit(i, 2) - limit(i, 1)) * rand(nPop, 1); %初始种群的位置 end v = rand(nPop, dim); % 初始化种群速度 xm = x;% 每个个体的历史最佳位置 ym = zeros(1, dim);% 种群的历史最佳位置 fxm = zeros(nPop, 1);% 每个个体的历史最佳适应度 fym = - inf;% 种群历史最佳适应度 for ii = 1:length(xm(:, 1)) x0 = ceil(xm(ii, 1)); y0 = ceil(xm(ii, 2)); xx(ii, 1) = x0_1(x0,y0); yy(ii, 1) = x0_2(x0,y0); zz(ii, 1) = f(x0,y0); end plot3(xx,yy,zz, 'r*');

修改后的代码如下：

nPop = 30; % 初始粒子群微粒个数 
dim = 2; % 空间维数 
maxIter1 = 300; % 最大迭代次数 
w = 0.9; % 惯性因子 
c1 = 2; % 个体学习因子 
c2 = 2; % 社会学习因子 
limit = [75, 150; -40, 40]; % 位置参数限制 
vlimit = [-5, 5; -5, 5]; % 速度参数限制  

% 生成初始种群的位置
for i = 1:dim     
    x(:,i) = limit(i, 1) + (limit(i, 2) - limit(i, 1)) * rand(nPop, 1); 
end 

v = rand(nPop, dim); % 初始化种群速度 
xm = x; % 每个个体的历史最佳位置 
ym = zeros(1, dim); % 种群的历史最佳位置 
fxm = zeros(nPop, 1); % 每个个体的历史最佳适应度 
fym = -inf; % 种群历史最佳适应度 

% 计算种群初始位置的适应度，并绘制散点图
for ii = 1:length(xm(:, 1))     
    x0 = ceil(xm(ii, 1));     
    y0 = ceil(xm(ii, 2));     
    xx(ii, 1) = x0;     
    yy(ii, 1) = y0;     
    zz(ii, 1) = f(x0, y0); 
end  
plot3(xx, yy, zz, 'r*');

修改内容：

1. 对于注释，将百分号后面的空格去掉。
2. 将初始种群位置的计算放到一个循环内，避免重复计算。
3. 修改了变量名 x0_1 和 x0_2，改为正确的 x0 和 y0。