pmx.m

patten regnization source从1-14章能运行

function [h1, h2] = pmx(p1, p2, varargin)
% [h1,h2] = pmx(p1, p2, crosspoint1, crosspoint2)
% [h1,h2] = pmx(p1, p2)
%
% p1客 p2狼 何盒 荤惑 背硅(PMX) 备泅  
% 背硅痢捞 林绢瘤瘤 臼栏搁 烙狼肺 急琶
% 背硅痢 crosspoint1绰 背硅痢 crosspoint2焊促 累酒具 窃.

if length(varargin) >= 2
  % if given, accept crosspoints
  c1 = varargin{1};
  c2 = varargin{2};
else
  % if not given, generate crosspoints
  c1 = floor(rand*(length(p1)-1))+1;
  c2 = floor(rand*(length(p2)-1));
  if c2 >= c1
    c2 = c2 + 1;
  else
    tmp = c1;
    c1 = c2;
    c2 = tmp;
  end
end

if c2 <= c1
  error('crossover points are incorrect (c2<=c1)')
end

h1 = p1;  % child1
h2 = p2;  % child2

for i=c1+1:c2
  j = find(h2==p1(i));%h2俊辑 p1(i)狼 蔼捞 乐绰 困摹阑 茫酒辑
  tmp = h2(j);%弊 镑狼 蔼阑 tmp俊 历厘窍咯
  h2(j) = h2(i);%困摹甫 官槽促.
  h2(i) = tmp;
  j = find(h1==p2(i));%h1俊辑 p2(i)狼 蔼捞 乐绰 困摹甫 茫酒辑
  tmp = h1(j);%弊 镑狼 蔼阑 tmp俊 历厘窍咯
  h1(j) = h1(i);%困摹甫 官槽促.
  h1(i) = tmp;
end