论坛君：尽管笛卡尔和瑞典公主的故事已被证实只是杜撰，但因这个故事出名的心形函数被广为流传。今天又是一个虐单身狗的日子，面对各种毫无新意的表白方式，让我们来看看理工科式的表白~

---

一切都是从一个故事开始的：

1650年，斯德哥尔摩街头，一个宁静的午后，笛卡尔邂逅了18岁的瑞典公主克里斯汀。

机遇巧合，一段纯粹、美好的爱情悄然萌发。

然而，没过多久，他们的恋情传到了国王的耳朵里。国王大怒，将笛卡尔放逐，而公主被软禁在宫中。 身体孱弱的笛卡尔不久便染上重病，在生命进入倒计时的那段日子，他日夜思念的还是街头偶遇的那张温暖的笑脸。他每天坚持给她写信，盼望着她的回音。然而，这些信都被国王拦截下来。在笛卡尔给克里斯汀寄出第十三封信后，他永远地离开了这个世界。此时，被软禁在宫中的小公主依然徘徊在皇宫的走廊里，思念着远方的情人。

这最后一封信上没有写一句话，只有一个方程：r=a(1-sinθ) ——传说中著名的心形函数。

不过，事实告诉我们，除非你的目标妹子也是一只Geeker（至少会用Mathematica或者MATLAB等软件），否则像笛卡尔这样单给一个函数的结局大概就是你推公式别人推妹子了……表白什么的还是选择更浅显易懂的方法吧。

接下来看看用各种数据分析软件做出的心形模型，保证亮瞎你的表白对象！

Stata

在这个特殊的日子里，特意搜集了利用Stata绘制心形线的程序，代码如下：


以下为输出结果：

以上程序基于stata12运行，如果有stata14的话，还可以输入中文！祝大家好运！

R

接下来是R绘制心形线，代码如下：


结果如下：

注：操作平台为Rstudio，由于手头电脑没有装R，所以上述程序没有验证，就交给大家自己实践了！祝好运！

以上程序来自：
http://mp.weixin.qq.com/s?__biz=MzA4NTAyMjQ2Mg==&mid=206788501&idx=1&sn=53e07723fdb395a787e57ea53bffd332&scene=21#wechat_redirect

Matlab

利用matlab抱得美人归：


效果如斯：

SAS

1）准备心形图的图形数据


2）准备图形上的文字

这个最重要，字一定要大，而且千万记得署上自己的名字。过年的时候收到几条没有署名的祝福短信，郁闷的我也不知道该怎么问人家名字。


3）画图


用SAS画曲线有多种方式：

1）用ScatterPlot画图是最自然的方式，不过需要10000个点心形图看上去才比较舒服。因为曲线的边有些锯齿，所以我不太喜欢。

2）用Series画图，这个需要点技巧。对数据不做特殊处理的话，画出来的是这样，不过也很美很特别。

Heart

3）我建议大家动手试试其他的语句，比如pbsplineplot，会有各种惊喜。

为了能出下面这个心形图，我特意加了side变量，然后通过lineattrs=(color=red)强行将两组数据用红色画线，而不是默认的一蓝一红。

Do you love SAS? Yes, I DO.（行动才是爱！）

Python

爱心的python表示：

1.第一种

2.第二种

3.第三种

4.第四种

最后，感谢那些做出巨大贡献的程序猿，没有他们的代码，表白怎能如此风骚！！！


（关注微信+10个论坛币）
链接地址：https://bbs.pinggu.org/thread-2192087-1-1.html