GeneticAlgorithm.rar_Hill-climbing_Matlab Code Hill _genetic_hil

function GeneticAlgorithm clear all; clc; %-----initialize------ chess = zeros(8); %put Queens Randomly for i=1:8 for j=1:8 chess(i,j)= randint(1,1,[1,8]); end end crossoverprob = .8; mutationprob = .01; %-----calculate current heuristic------ disp(chess); for i=1:8 current_heuristic(i) = computeHeuristic(chess(i,:)); fittness(i) = 28 - current_heuristic(i); end step_count = 0; while(1) step_count = step_count +1; [IsFound,k] = IsOneHeuristic(current_heuristic); if (IsFound) == true chess(k,:) break; end %SELECTION------------------------------------------------------- % find each state probability for i=1:8 prob(i) = fittness(i)/(sum(fittness)); end for i=1:8 point(i) =sum(prob(1:i)); end %Choose 8 number randomly and set their Domains and set new generation newchess = [,]; for i=1:8 random(i) = rand(1); Domain(i) = GetDomain(random(i),point); newchess(i,:) = chess(Domain(i),:); end %CROSS-OVER------------------------------------------------------ for i=1:crossoverprob*8 randdigit = randint(1,1,[1,8]); rand1 = randint(1,1,[1,8]); rand2 = randint(1,1,[1,8]); while(rand1 == rand2) rand2 = randint(1,1,[1,8]); end newchess(rand1,:) = cat(2,newchess(rand1,1:randdigit),newchess(rand2,randdigit+1:8)); newchess(rand2,:) = cat(2,newchess(rand2,1:randdigit),newchess(rand1,randdigit+1:8)); end %MUTATION-------------------------------------------------------- mute = rand(8); for i=1:8 for j=1:8 if(mute(i,j))<(mutationprob) muterow(i,j) = 1; newchess(i,j) = randint(1,1,[1,8]); else muterow(i,j) = 0; end end end chess = newchess; for i=1:8 current_heuristic(i) = computeHeuristic(chess(i,:)); fittness(i) = 28 - current_heuristic(i); end end step_count %Function Defenition---------------------------------------------------- function heuristic = computeHeuristic(chess) heuristic = 0; for i=1:7 rowi = chess(i); for j=i+1:8 rowj = chess(j); if(rowi-rowj == 0) || (abs(rowi-rowj)==j-i) heuristic = heuristic+1; end end end function [IsFound,k]=IsOneHeuristic(current_heuristic) k = 0; IsFound = false; for i=1:8 if(current_heuristic(i)) == 0 IsFound = true; k = i; break; end end function Domain = GetDomain(r,point) if (r) <= point(1) Domain = 1; else for i=2:8 if (point(i-1) < r) && (r <= point(i)) Domain = i; break; end end end