pma.m

自己编写的改进人口迁移算法求解函数最优值

%用人口迁移算法求函数的最大值
%求目标函数的值
format long
N=3;%人口规模
l=10;%人口流动次数
k=0.01;%收缩系数
h=1e-8;%人口压力参数
m=10;%迭代次数
a1=3;%搜索空间
b1=6;%搜索空间
a2=3;%搜索空间
b2=5;%搜索空间
center1=[];
center2=[];
zy=[];%记录最优值
zyd=[];%记录最优点
zydd=[];%记录最优点
d=[];%计算目标函数值
kk=0;
X1=[];
X2=[];
tt=0;
%步骤1开始
while(kk<m)
X1=(rand(1,N)*(b1-a1))+a1;%随机产生搜索空间的N个点
X2=(rand(1,N)*(b2-a2))+a2;%随机产生搜索空间的N个点
c1=(b1-a1)/(2*N);%区域的邻域
c2=(b2-a2)/(2*N);%区域的邻域
%步骤1结束
%步骤2开始
for i=1:1:N
center1(i)=X1(i);
center2(i)=X2(i);
d(i)=maxf(X1(i),X2(i));%计算各点的函数值
end
%步骤2结束
%步骤3开始
[w,i]=min(d); %记录最优函数值和最优点
zy(numel(zy)+1)=w;%记录最优函数值和最优点
zyd(numel(zyd)+1)=X1(i);%记录最优函数值和最优点
zydd(numel(zydd)+1)=X2(i);
%步骤3结束
while(1)
for i=1:1:N
%步骤4开始
X1(i)=2*c1*rand+(center1(i)-c1);
X2(i)=2*c2*rand+(center2(i)-c2);
if(X1(i)>b1)
    X1(i)=b1;
end
if(X1(i)<a1)
    X1(i)=a1;
end
if(X2(i)>b2)
    X2(i)=b2;
end
if(X2(i)<a2)
    X2(i)=a2;
end
%步骤4结束
%步骤5开始
d(i)=maxf(X1(i),X2(i));
%步骤5结束
end
%步骤6开始
[w,i]=min(d); %记录最优函数值和最优点
zy(numel(zy)+1)=w;%记录最优函数值和最优点
zyd(numel(zyd)+1)=X1(i);%记录最优函数值和最优点
zydd(numel(zydd)+1)=X2(i);
%步骤6结束
%步骤7开始
tt=tt+1;
if(tt>=l)
    break
    
end
end
%步骤7结束
%步骤8开始
[w,i]=min(zy);%读取最优值和最优点
yyd=zyd(i);
yydd=zydd(i);
a11=yyd-c1;%优惠区域
a21=yydd-c2;
if(a11>=b1)
    a11=b1-c1;
end
if(a11<a1)
    a11=a1;
end
if(a21>=b2)
    a21=b2-c2;
end
if(a21<a2)
    a21=a2;
end
b11=yyd+c1;%优惠区域
b21=yydd+c2;
if(b11>b1)
    b11=b1;
end
if(b11<=a11)
    b11=a11+c1;
end
if(b21>b2)
    b21=b2;
end
if(b21<=a21)
    b21=a21+c2;
end
X1=(rand(1,N)*(b11-a11))+a11;%随机产生新搜索空间的N个点，替换原来的点
X2=(rand(1,N)*(b21-a21))+a21;%随机产生新搜索空间的N个点，替换原来的点
%步骤8结束
%步骤9开始
for i=1:1:N
center1(i)=X1(i);%第i个区域的中心
center2(i)=X2(i);%第i个区域的中心
d(i)=maxf(X1(i),X2(i)); %计算各点的函数值
end
%步骤9结束
%步骤10开始
[w,i]=min(d); %记录最优函数值和最优点
zy(numel(zy)+1)=w; %记录最优函数值和最优点
zyd(numel(zyd)+1)=X1(i);%记录最优函数值和最优点
zydd(numel(zydd)+1)=X2(i);
%步骤10结束
%步骤11开始
while(1)
c1=(1-k)*c1;%收缩优惠区域
c2=(1-k)*c2;%收缩优惠区域
%步骤11结束
%步骤12开始
a12=zyd(numel(zyd))-c1;
a22=zydd(numel(zydd))-c2;
if(a12>=b1)
    a12=b1-c1;
end
if(a12<a1)
    a12=a1;
end
if(a22>=b2)
    a22=b2-c2;
end
if(a22<a2)
    a22=a2;
end
b12=zyd(numel(zyd))+c1;%优惠区域
b22=zydd(numel(zydd))+c2;%优惠区域
if(b12>b1)
    b12=b1;
end
if(b12<=a12)
    b12=a12+c1;
end
if(b22>b2)
    b22=b2;
end
if(b22<=a22)
    b22=a22+c2;
end
X1=(rand(1,N)*(b12-a12))+a12;%随机产生新搜索空间的N个点，替换原来的点
X2=(rand(1,N)*(b22-a22))+a22;%随机产生新搜索空间的N个点，替换原来的点
%步骤12结束
%步骤13开始
for i=1:1:N
center1(i)=X1(i);
center2(i)=X2(i);
d(i)=maxf(X1(i),X2(i)); %计算各点的函数值
end
%步骤13结束
%步骤14开始
[w,i]=min(d); %记录最优函数值和最优点
zy(numel(zy)+1)=w ;  %记录最优函数值和最优点
zyd(numel(zyd)+1)=X1(i);%记录最优函数值和最优点
zydd(numel(zydd)+1)=X2(i);
%步骤14结束
%步骤15开始
if(~(max(c1,c2)>h))
    break
end
end
kk=kk+1;
end
%步骤15结束
%步骤16开始
%步骤16结束
[w,i]=min(d);
zuiyouzhi=w
zuiyoudian1=X1(i)
zuiyoudian2=X2(i)