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;  endendif c2 <= c1  error('crossover points are incorrect (c2<=c1)')endh1 = p1;  % child1h2 = p2;  % child2for 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