www.gusucode.com > 几种多目标优化算法集合,包括MOEAD,MOPSO,NNIA,NSGA2等 > MOEAD/realmutate.m
function ind = realmutate(ind, domains, rate) %REALMUTATE Summary of this function goes here % Detailed explanation goes here % double rnd, delta1, delta2, mut_pow, deltaq; % double y, yl, yu, val, xy; % double eta_m = id_mu; eta_m=20; numVariables = size(domains,1); if (isstruct(ind)) a = ind.parameter; else a = ind; end for j = 1:numVariables if (rand() <= rate) y = a(j); yl = domains(j,1); yu = domains(j,2); delta1 = (y - yl) / (yu - yl); delta2 = (yu - y) / (yu - yl); rnd = rand(); mut_pow = 1.0 / (eta_m + 1.0); if (rnd <= 0.5) xy = 1.0 - delta1; val = 2.0 * rnd + (1.0 - 2.0 * rnd) * (xy^(eta_m + 1.0)); deltaq = (val^mut_pow) - 1.0; else xy = 1.0 - delta2; val = 2.0 * (1.0 - rnd) + 2.0 * (rnd - 0.5) * (xy^ (eta_m + 1.0)); deltaq = 1.0 - (val^mut_pow); end y = y + deltaq * (yu - yl); if (y < yl) y = yl; end if (y > yu) y = yu; end a(j) = y; end end if isstruct(ind) ind.parameter = a; else ind = a; end end