#Browni liikumine
#Ühel teljel
#500 m/s
#68 nm vaba tee keskmine pikkus
#7.3*10^2 põrget sekundis

liikumised=rnorm(20, 0, 40) #mõlemas suunas sd 40 nm
algkoht=500 #nm
filter(liikumised, 1, method="recursive", init=algkoht)

filter(rnorm(20, 0, 40), 1, method="recursive", init=500)

plot(filter(rnorm(20, 0, 40), 1, method="recursive", init=500))

veerge=4
ridu=3
par(mfrow=c(ridu, veerge))
for(i in 1:(veerge*ridu)){
  plot(filter(rnorm(20, 0, 40), 1, method="recursive", init=500),
    ylim=c(0, 1000), ylab="nm")
}
par(mfrow=c(1, 1))


#filter(
#  c(rnorm(20, 0, 40), rnorm(20, 0, 40)),
#  c(1, 1),
#  method="recursive", 
#  init=c(500, 500))

sapply(1:2, function(d) filter(rnorm(20, 0, 40), 1, method="recursive", init=500))

#Tehke tsükliga samale lehele mitu kahemootmelise Browni liikumise joonist
png("c:/jaagup/2016/1/brown.png")
plot(sapply(1:2, function(d) filter(rnorm(200, 0, 40), 1, method="recursive", init=500)), 
     type="l", xlab="x", ylab="y")
dev.off()

par(mfrow=c(2, 2))
for(i in 1:4){
 plot(sapply(1:2, function(d) filter(rnorm(200, 0, 40), 1, method="recursive", init=500)), 
     type="l", xlab="x", ylab="y")
}


sapply(nurgad, sin)

samm=5
funktsioonid=c(sin, cos)
nurgad=rnorm(10, 0, 0.3)
muutused=sapply(1:2, function(nr) sapply(nurgad, funktsioonid[[nr]]))
muutused
muutused[ ,2]
cumsum(muutused[, 1])
filter(muutused[, 1], 1, method="recursive", init=500) #x-koordinaadid
filter(muutused[, 2], 1, method="recursive", init=500) #y-koordinaadid


algkoht=c(500, 500)
sapply(1:2, function(nr) filter(muutused[,nr], 1, method="recursive", init=algkoht[nr]))
plot(sapply(1:2, 
  function(nr) filter(muutused[,nr], 1, method="recursive", init=algkoht[nr]))
   , type="l")


funktsioonid=c(cos, sin)
nurgamuutused=rnorm(1000, 0, 0.3) #igakordne keeramine
nurgamuutused
nurgad=filter(nurgamuutused, 1, method="recursive", init=0)
nurgad
muutused=sapply(1:2, function(nr) sapply(nurgad, funktsioonid[[nr]])) # x ja y-koordinaadi muutused
muutused
algkoht=c(500, 500)
plot(sapply(1:2, function(nr) 
   filter(muutused[,nr], 1, method="recursive", init=algkoht[nr])), type="l")
 

#veebikeskkonna install
#looge kaust moodulite hoidmiseks (nt. c:/python34/rmoodulid)
.libPaths("c:/python34/rmoodulid")
install.packages("shiny")
library(shiny)