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