function [Swarm, PBEST] = initiate(psoOptions,Obj)
Swarm = struct( 'individual',struct(...
  'dim' ,struct(...
    'number',[], ...%個體位置
    'speed',[], ...%個體速度
    'best',[]),'number',[],'speed',[],'best',[]), ...%個體最優
   'P',[], ... 
    'X',[],'Y',[]);
%X是粒子群活動的范圍的初始值。
B=ones(1,Obj.Dim,psoOptions.population);
for j = 1 : psoOptions.population

    Swarm(j).X = ( ones(psoOptions.Vars.PopSize, 1) * (Obj.ub - Obj.lb) ) .* rand(psoOptions.Vars.PopSize, Obj.Dim) + ones(psoOptions.Vars.PopSize, 1) * Obj.lb;
    for i = 1 : psoOptions.Vars.PopSize 
      for d=1:Obj.Dim 
        Swarm(j).individual(i).dim(d).number = Swarm(j).X(i,d);            
        Swarm(j).individual(i).dim(d).speed = rand * psoOptions.SParams.Vmax;             
        Swarm(j).individual(i).dim(d).best = Swarm(j).X(i,d);    
        
     B(i,d,j)=Swarm(j).individual(i).dim(d).best;
      
    
     %  PBEST(j,d)=PBEST;   
       
      Swarm(j).individual(i).number=[Swarm(j).individual(i).dim.number];
      Swarm(j).individual(i).best=[Swarm(j).individual(i).dim.best];
    % Swarm(j).individual(i).P=[Swarm(j).individual(i).dim(d).best];
      end
                       Swarm(j).P =[Swarm(j).individual(i).dim.best];
     Swarm(j).individual(i).number=[Swarm(j).individual(i).dim.number];
     Swarm(j).individual(i).best  =[Swarm(j).individual(i).dim.best];
    end
     
    %Swarm(j).Y= feval('LJcluster', Swarm(j).X(i,:));
    %[PBEST, index] = min(Swarm(j).Y);%pbest是一個1。10列的向量。index是最小值的索引
   %是一個j群的全局最優的實數
end
B=min(B);
PBEST=feval('LJcluster', B);