我想在方形格子上模拟一个自我避免的二维随机游走并绘制路径。
到目前为止,我已经写出了这个问题的代码:
n <- 100
x <- 0
y <- 0
randomwalkx <- 0
randomwalky <- 0
for(i in 1:n){
random <- sample(0:3, 1)
if(random == 0){x <- x+1}
if(random == 1){x <- x-1}
if(random == 2){y <- y+1}
if(random == 3){y <- y-1}
xcheck <- x %in% randomwalkx
ycheck <- y %in% randomwalky
if(ycheck == "TRUE" && xcheck == "TRUE"){
if(random==0){x<-x-1}
if(random==1){x<-x+1}
if(random==2){y<-y-1}
if(random==3){y<-y+1}
}else{
randomwalkx <- c(randomwalkx, x)
randomwalky <- c(randomwalky, y)
}
}
plot(randomwalkx,randomwalky,xlab='x',ylab='y')
lines(randomwalkx,randomwalky,type='b')
然而,该路径仅沿一个对角线方向行进。谁能看到我犯的错误或解决这个问题?