PMX算法MATLAB实现

function [newVanewVb]=PMX（VaVb） 
    startXorPoint=mod（ceil（rand（1）*10）length（Va） ）; 
    if startXorPoint==0 
        startXorPoint=startXorPoint+1; 
    end    
    xorLength=mod（floor（rand（1）*10）length（Va））; 
    endXorPoint=startXorPoint+xorLength; 
    while（endXorPoint>length（Vb） ） 
        xorLength=mod（floor（rand（1）*10）length（Va））; 
        endXorPoint=startXorPoint+xorLength; 
    end    
    tempa=Va（startXorPoint:endXorPoint）; 
    tempb=Vb（startXorPoint:endXorPoint）; 
    %找出映射关系
    for ix=1:length（tempa） 
        aMappingb（ix1:2）=[Va（startXorPoint+ix-1）Vb（startXorPoint+ix-1）];
    end
    for ix=1:length（tempa） 
        bMappinga（ix1:2）=[Vb（startXorPoint+ix-1）Va（startXorPoint+ix-1）];
    end
    %交换子串
    Va（startXorPoint:endXorPoint）=tempb; 
    Vb（startXorPoint:endXorPoint）=tempa; 
    
    %对Va进行冲突检测
    indexVa = 1;
    while（indexVa        singleGene =  Va（indexVa）;
        while（ ismember（singleGenetempb） == 1 ）
            conflictIndex = find（tempb==singleGene）;
            singleGene = bMappinga（conflictIndex2）;
        end
        V