function [pbest,Aim_best]=PSO(Aim,Num,M,D,C1,C2,K)
% function [pbest,Aim_best]=PSO(Aim,Num,M,C1,C2,D,K,Vmin,Vmax,Xmin,Xmax)
%粒子群优化算法
%pbest为全局最优解
%Aim_best为最优解时的函数值
%Aim是目标函数
%Num是粒子个数
%M为最大迭代次数
%D自变量个数
%C1、C2为学习因子
%K是收敛因子
%Vmax为最大速度
%Vmin
%Xmax
%Xmin
D=6;
%初始化种群
for i=1:1:Num
for j=1:1:D
x(i,j)=randn; %随机初始化位置
v(i,j)=randn; %随机初始化速度
end
end
%计算各粒子适应度，初始化pi和pg
for i=1:1:Num
p(i)=Aim(x(i,:));%目标函数最优值
y(i,:)=x(i,:);%最优位置pbest
end
pg=x(Num,:);%全局最优位置
for i=1:(Num-1)
if Aim(x(i,:))<Aim(pg)
pg=x(i,:);
end
end
%主循环
for i=1:1:M
for j=1:1:Num
v(j,:)=K*(v(j,:)+C1*rand*(y(j,:)-x(j,:))+C2*rand*(pg-x(j,:)));%更新速度
x(j,:)=x(j,:)+v(j,:);%更新位置
if(Aim(x(j,:))<p(j))
p(j)=Aim(x(j,:));
y(j,:)=x(j,:);
end
if(Aim(x(j,:))>p(j))
pg=y(j,:);
end
end
end
pbest=pg';
Aim_best=Aim(pg);
