pma3.asv

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

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