tests/RastriginAux.R

In elthonf/EMFGeneticos: Algoritmos Geneticos para o curso da USP-EACH

x<-seq(-5.12,5.12,length=100)

y<-x

f<-function(x,y) {
    20+(x^2-10*cos(2*3.14*x))+
       (y^2-10*cos(2*3.14*y)) }

z<-outer(x,y,f)

z[is.na(z)]<-1

persp(x,y,z,theta=30,phi=30,expand=0.5,col="red",
         ltheta=90,shade=0.50,ticktype="detailed",d=5,r=1)



#x %o% y
?persp



z <- 2 * volcano        # Exaggerate the relief
x <- 10 * (1:nrow(z))   # 10 meter spacing (S to N)
y <- 10 * (1:ncol(z))   # 10 meter spacing (E to W)
## Don't draw the grid lines :  border = NA
par(bg = "slategray")
persp(x, y, z, theta = 135, phi = 30, col = "green3", scale = FALSE,
      ltheta = -120, shade = 0.75, border = NA, box = FALSE)




par(bg = "white")
x <- seq(-1.95, 1.95, length = 30)
y <- seq(-1.95, 1.95, length = 35)
z <- outer(x, y, function(a, b) a*b^2)
nrz <- nrow(z)
ncz <- ncol(z)
# Create a function interpolating colors in the range of specified colors
jet.colors <- colorRampPalette( c("blue", "green") )
# Generate the desired number of colors from this palette
nbcol <- 100
color <- jet.colors(nbcol)
# Compute the z-value at the facet centres
zfacet <- z[-1, -1] + z[-1, -ncz] + z[-nrz, -1] + z[-nrz, -ncz]
# Recode facet z-values into color indices
facetcol <- cut(zfacet, nbcol)
persp(x, y, z, col = color[facetcol], phi = 30, theta = -30)

par(op)


volcano