程序员的浪漫--R语言画小心心

2019-06-05 本文已影响18人寥廓江天万里霜

Fig 1.

n=50000
r=0.6
r2=(1-r*r)^0.5
X=rnorm(n)
Y=X*r+r2*rnorm(n)
Y=ifelse(X>0,Y,-Y)
plot(X,Y,col="red")

1.JPG

Fig 2

x <- seq(-1,1,length = 5000)
y1 <- (1-x²⁾(1/2)+(x²⁾(1/3)
y2 <- -(1-x²⁾(1/2)+(x²⁾(1/3)
plot(c(rev(x),x),c(rev(y2),y1), type = 'l' ,xlim = c(-1,1), xlab = "",ylab = "", axes = F,col=2)

好像有点单调对吧？？？

那我们来参考一下加州伯克利大佬们的浪漫

library("RColorBrewer")
col<-brewer.pal(3,'Set1')[1]

RColorBrewer这个包可以玩好久

t<-seq(0,60,len=100)
x<--0.01(-t^2+40t+1200)sin(pit/180)
y<- 0.01(-t^2+40t+1200)cos(pit/180)
plot(c(-12,14),c(0,20),type = 'n',axes = T,xlab = '',ylab = '')

一定要先画箭头

arrows(-11,8,14,8,lwd=8,col='gold')
arrows(14,8,-11,8,lwd=8,col='gold',angle = 120)
polygon(x,y,col=col,border = NA)
polygon(-x,y,col=col,border = NA)
lines(x,y,lwd=1)
lines(-x,y,lwd=1)