m0_62740426 2022-08-10 15:35 采纳率: 33.3%
浏览 28
已结题

位置 1 中的索引超过数组边界(不得超过 1报错

一个简单的粒子群算法,在第32行标注的地方一直报错
我检查了这些数组的大小都是一样的啊,为啥会报错 ,看了我一个小时555

img

img

clc
clear

c1 = 1.4955;
c2 = 1.4955;
w = 1;

maxgen = 100;
sizepop = 50;
nvar = 2;
vmax = 5 ;
vmin = -5;
popmax = 100;
popmin = -100;
pbest = [];
for i = 1 : sizepop
pop(i,:) = (popmax-popmin) * rand(1,nvar) + popmin;
V(i,:) = (vmax-vmin) * rands(1,nvar) + vmin;
funtion(i) = fun(pop(i,:));
fitness(i) = 1/funtion(i);
end
[bestfitness,bestindex] = max(fitness);
bestfuntion = funtion(bestindex);
bestpop = pop(bestindex,:);
pbest = pop;
pbestfuntion = funtion;
pbestfitness = fitness;

for i = 1 : maxgen

for j = 1 : sizepop
V(j,:) = w.V(j,:) + c1rand*(pbest(j,:)-pop(j,:))+c2rand(bestpop-pop(j,:));

    V(j, find(V(j,:)) > vmax) = vmax;
    V(j, find(V(j,:)) < vmin) = vmin;
    
    pop(j,:) = pop(j,:) + V(j,:);
  
    pop(j,find(pop(j,:) > popmax)) = popmax;
    pop(j,find(pop(j,:) < popmin)) = popmin;     
    funtion(j) = fun(pop(j,:));
    fitness(j) = 1/funtion(j);
 

end
for j = 1 : sizepop
if fitness(j) > pbestfitness(j)
pbestfitness(j) = fitness(j);
pbestfuntion(j) = funtion(j);
pbest = pop(j,:);

   end
   if   fitness(j) > bestfitness
       bestfitness = fitness(j);
       bestfuntion = funtion(j);
       bestpop = pop(j,:);
   end

end

end

  • 写回答

1条回答 默认 最新

  • falwat 2022-08-10 15:49
    关注

    不好排查问题出在哪儿。建议你自己在出错的地方设置一个断点,然后在命令窗口里依次测试表达式的每个运算

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 8月18日
  • 已采纳回答 8月10日
  • 创建了问题 8月10日

悬赏问题

  • ¥15 MATLAB动图问题
  • ¥15 有段代码不知道怎么理解,const isToken = (config.headers || {}).isToken === false
  • ¥100 采用栈和深度优先算法对我的代码进行修改显示路径
  • ¥15 pointnet2包安装
  • ¥20 射频功率问题,解答者有酬谢!
  • ¥80 构建降雨和积水的预测模型
  • ¥15 #Qt Transform setTransform()在鼠标拖动移动视角是一致在原地不动,无法变换视角(细微观察似乎视图有在原地抖动),无法变换视角(细微观察似乎视图有在原地抖动)
  • ¥50 如何利用无人机拍摄的数码照片测量鸟卵的长短径
  • ¥100 github贡献者给与奖励
  • ¥15 使用DS18B20+ESP8266获取温度数据返回-127.00