while gen < MAXGEN;   
      for a=1:1:NIND  
          X      = phen1(a,:);
          %約束設(shè)置
          W_     = reshape(X,[M,M]);
          trw    = trace(W_);%約束1
          [V,I]  = eig(W_);
          lmax   = max(max(I));
          y1     = func_Fk(X,WW1,rr,lmax,bs{1},As{1});%約束3
          y2     = func_Fk(X,WW1,rr,lmax,bs{2},As{2});%約束3
          ta1    = mean(X.*ga{1,1});
          tb1    = mean(X.*gb{1,1});
          tc1    = mean(X.*gc{1,1});
          td1    = mean(X.*gd{1,1});
          ta2    = mean(X.*ga{2,2});
          tb2    = mean(X.*gb{2,2});
          tc2    = mean(X.*gc{2,2});
          td2    = mean(X.*gd{2,2});
          
          if abs(trw-1)<=0.1 & real(y1)<=0 & real(y2)<=0
              %計(jì)算對(duì)應(yīng)的目標(biāo)值
              epls        = func_obj(X,h__);
              E           = abs(epls);
              JJ(a,1)     = E;
              WW(a,1:16)  = X;
          else
              %計(jì)算對(duì)應(yīng)的目標(biāo)值
              epls        = func_obj(X,h__);
              E           = abs(epls);
              JJ(a,1)     = E/1e3;
              WW(a,1:16)  = X;
          end
      end 
      
      Objvsel=(JJ);    
      [Chrom,Objv]=reins(Chrom,Selch,1,1,Objv,Objvsel);   
      gen=gen+1; 
 
      idx = find(JJ==0);
 
      
      Error(gen) = 1/mean(JJ)-1;
      [VV,II]=min(JJ);
end 
%輸出w
WW(idx,:)= [];
[VV,II]  = min(JJ);
Wopt     = WW(II,:);