整数空间粒子群算法MATLAB程序

整数空间粒子群算法clear allclcticT =1200; %最大迭代次数PopSize=40; %种群大小c1 =2; %学习因子1c2 = 2; %学习因子2w=0.9; %惯性因子NDim = 60; %空间维数 for i=1:PopSize for j=1:NDim %初始化各粒子初始位置，在有效范围内随机选数 x(i,j)=round(rand*20) + 20; endendvmax=zeros(PopSize,NDim);for i=1:PopSize for j=1:NDim %初始化各粒子最大速度，使粒子不能越出边界 vmax(i,j)=10; end endv = vmax.*rand( PopSize,NDim); %计算初始各粒子的适应度值for i = 1:PopSize fvalue(i) = youhualianxuduo(x(i,:); endpbest = x; %记录各粒子的个体极值点位置fpbest = fvalue; %记录初始个体最佳适应度值fbestval,index = min(fvalue); % 找出全局极值和相应的序号 gbest1=fbestval;g1=x(index,:);%浸入主循环，按照公式依次迭代直到满足 迭代次数for i=1:T for j=1:PopSize if(youhualianxuduo(x(j,:)fpbest(j) pbest(j,:)=x(j,:); fpbest(j)=youhualianxuduo(x(j,:); end if (fpbest(j)x(j,d) a1=0; else a1=c1*(pbest(j,d)-x(j,d); end if pbest(j,d)x(j,d) b1=c1*(pbest(j,d)-x(j,d); else b1=0; end if g1(d)x(j,d) a2=0; else a2=c2*(g1(d)-x(j,d); end if g1(d)x(j,d) b2=c2*(g1(d)-x(j,d); else b2=0; end R1=randi(a1,b1,1,1); R2=randi(a2,b2,1,1); v(j,d)=round(w*v(j,d)+R1+R2; end%更新速度V for n=1:PopSize %速度v的限定，最大值为10 for k=1:NDim if v(n,k)10 v(n,k)=10; end if v(n,k)-10 v(n,k)=-10; end end end x(j,:)=x(j,:)+v(j,:); %更新位置x % 逸出标志 for l=1:PopSize % 阻止逸出 for m=1:NDim if x(l,m)40 x(l,m)=round(rand*20+20); end end end end w=0.9-i*0.5/T gb1(i)=gbest1;endfigure(7)plot(gb1);TempStr=sprin