巴恩斯利蕨 The Barnsley Fern

巴恩斯利蕨学习链接

下面用R画一下:

npts<-50000

point_mat<-matrix(data=NA,nrow=npts,ncol=2)

list_fun<-list(

    f1=function(x,y){c(0., 0.16*y)}  # R里的匿名函数

    f2=function(x,y){c(0.85*x + 0.04*y, -0.04*x + 0.85*y + 1.6)},

    f3=function(x,y){c(0.2*x - 0.26*y, 0.23*x + 0.22*y + 1.6)},

    f4=function(x,y){c(-0.15*x + 0.28*y, 0.26*x + 0.24*y + 0.44)}

)

x<-0;y<-0

for (i in 1:npts) {

   point_mat[i,]<-c(x,y)

   random<-sample(list_fun,size=1,replace=TRUE,

                  prob=c(0.01, 0.85, 0.07, 0.07))[[1]]

   vec<-random(x,y)

   x<-vec[1];y<-vec[2]

}

# 画图

plot(point_mat[,1],point_mat[,2],

     type='p',pch=19,col='#008000',cex=0.1,

     xlab='',ylab='')

图形如下：

巴特西

巴恩斯利蕨 The Barnsley Fern

最新文章

热门文章