Week 4.1 - Genetic Algorithm

clear all 
close all
clc
%Defining our search space
x= linspace (0,0.6,150);
y= linspace (0,0.6,150);

%Creating a 2Dimensional mesh
[xx yy] =meshgrid(x,y);

%Evaluating the stalagamite function
for i=1:length(xx)
    for j= 1: length(yy)
        input_vector(1)=xx(i,j);
        input_vector(2)=yy(i,j);
        f(i,j)=stalagmite(input_vector);
    end
end
surf (xx,yy,f);
shading interp;
[x,fval]= ga(@stalagmite,2)

function [f] = stalagmite(input_vector)
x= input_vector(1);
y= input_vector(2);


f1_x = (sin(5.1*pi*x + 0.5)).^6;

f1_y = (sin(5.1*pi*y + 0.5)).^6;

t1 = ((x-0.0667).^2)/(0.64);

f2_x = exp(-4.*log(2)*t1);

t2 = ((y-0.0667).^2)/(0.64);

f2_y = exp(-4.*log(2)*t2);

f = -(f1_x.*f2_x.*f1_y.*f2_y);

end