最近技能树重启了每个月一次的赠书活动:数学不难,恋爱才难(文末赠送通关秘籍5本!) 可能是因为提到了恋爱,然后有人留言提到了在R里面绘制一个爱心: 还蛮好玩的,激起了广大粉丝的创造热情,5本书已经送出去了,让我们看看大家的创意吧!
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")
Fig 2x <- seq(-1,1,length = 5000) y1 <- (1-x<sup>2)</sup>(1/2)+(x<sup>2)</sup>(1/3) y2 <- -(1-x<sup>2)</sup>(1/2)+(x<sup>2)</sup>(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+40*t+1200)*sin(pi*t/180) y<- 0.01*(-t^2+40*t+1200)*cos(pi*t/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)
Fig3参考文献: https://blog.csdn.net/kMD8d5R/article/details/79327603 https://www./make-a-valentines-heart-with-r/
然后是龙妈的:值得一提的是这个颜值不错,如下: 据说还有分手函数图。。真的是。。。分手函数:17 x^2 – 16|x|y + 17 y^2 + 150/|5 x + sin(5 y)| < 225http://www./blog/archives/4139 PS: 据说笛卡尔最后还是没能等到 自己等了很久的那个人 但是这不重要,重要的是。。。群里好多有心型图啊 不知道作者们愿不愿意放上来分享一下(哈哈哈哈) 各位端午节快乐,没脱单的早日脱单 以及,下次再有人问你R语言能干什么,告诉他可以拿去表白 和大神们比画图我还是嫩了点(顶锅盖逃走)
|