pma.asv

用MATLAB编写的

%用人口迁移算法求函数的最大值
%求目标函数的值
%function y=maxf(x)
%y=(1-2*sin(3*pi*x)^20+sin(20*pi*x)^20)^20;
format long
N=4;%人口规模
l=20;%人口流动次数
k=0.001;%收缩系数
h=1e-6;%人口压力参数
m=20;%迭代次数
a=0;%搜索空间
b=1;%搜索空间
center=[];
zy=[];%记录最优值
zyd=[];%记录最优点
d=[];%计算目标函数值
kk=0;
X=[];
tt=0;
%步骤1开始
while(kk<m)
X=(rand(1,N)*(b-a))+a;%随机产生搜索空间的N个点
c=(b-a)/(2*N);%区域的邻域
%步骤1结束
%步骤2开始
for i=1:1:N
center(i)=X(i);
d(i)=maxf(X(i));%计算各点的函数值
end
%步骤2结束
%步骤3开始
[w,i]=min(d); %记录最优函数值和最优点
zy(numel(zy)+1)=w;%记录最优函数值和最优点
zyd(numel(zyd)+1)=X(i);%记录最优函数值和最优点
%步骤3结束
while(1)
for i=1:1:N
%步骤4开始
X(i)=2*c*rand+(center(i)-c);
if(X(i)>b)
    X(i)=b;
end
if(X(i)<a)
    X(i)=a;
end
%步骤4结束
%步骤5开始
d(i)=maxf(X(i));
%步骤5结束
end
%步骤6开始
[w,i]=min(d); %记录最优函数值和最优点
zy(numel(zy)+1)=w;%记录最优函数值和最优点
zyd(numel(zyd)+1)=X(i);%记录最优函数值和最优点
%步骤6结束
%步骤7开始
tt=tt+1;
if(tt>=l)
    break
    
end
end
%步骤7结束
%步骤8开始
[w,i]=min(zy);%读取最优值和最优点
yyd=zyd(i);
a1=yyd-c;%优惠区域
if(a1>=b)
    a1=b-c;
end
if(a1<a)
    a1=a;
end
b1=yyd+c;%优惠区域
if(b1>b)
    b1=b;
end
if(b1<=a1)
    b1=a1+c;
end
X=(rand(1,N)*(b1-a1))+a1;%随机产生新搜索空间的N个点，替换原来的点
%步骤8结束
%步骤9开始
for i=1:1:N
center(i)=X(i);%第i个区域的中心
d(i)=maxf(X(i)); %计算各点的函数值
end
%步骤9结束
%步骤10开始
[w,i]=min(d); %记录最优函数值和最优点
zy(numel(zy)+1)=w; %记录最优函数值和最优点
zyd(numel(zyd)+1)=X(i); %记录最优函数值和最优点
%步骤10结束
%步骤11开始
while(1)
c=(1-k)*c;%收缩优惠区域
%步骤11结束
%步骤12开始
a2=zyd(numel(zyd))-c;
if(a2>=b)
    a2=b-c;
end
if(a2<a)
    a2=a;
end
b2=zyd(numel(zyd))+c;%优惠区域
if(b2>b)
    b2=b;
end
if(b2<=a2)
    b2=a2+c;
end
X=(rand(1,N)*(b2-a2))+a2;%随机产生新搜索空间的N个点，替换原来的点
%步骤12结束
%步骤13开始
for i=1:1:N
center(i)=X(i);
d(i)=maxf(X(i)); %计算各点的函数值
end
%步骤13结束
%步骤14开始
[w,i]=min(d); %记录最优函数值和最优点
zy(numel(zy)+1)=w ;  %记录最优函数值和最优点
zyd(numel(zyd)+1)=X(i) ; %记录最优函数值和最优点
%步骤14结束
%步骤15开始
if(~(c>h))
    break
end
end
kk=kk+1;
end
%步骤15结束
%步骤16开始
%步骤16结束
[w,i]=min(d);
zuiyouzhi=w
zuiyoudian=X(i)