六行python代码的爱心曲线_6行python代码的爱心曲线
polar-heart.png
前些日子在做绩效体系的时候,遇到了一件囧事,居然忘记怎样在Excel上拟合正态分布了,尽管在第二天重新拾起了Excel中那几个常见的函数和图像的做法,还是十分的惭愧。实际上,当时有效偏颇了,忽略了问题的本质,解决数据分析和可视化问题,其实也是Python的拿手好戏。
例如,画出指定区间的一个多项式函数:
Python 代码如下:
import numpy as np
import matplotlib.pyplot as plt
X = np.linspace(-4, 4, 1024)
Y = .25 * (X + 4.) * (X + 1.) * (X - 2.)
plt.title('$f(x)=\\frac{1}{4}(x+4)(x+1)(x-2)$')
plt.plot(X, Y, c = 'g')
plt.show()
通过numpy的linspace方法来确定横坐标x的取值范围,列出方程,然后调用matplotlib的pyplot画出函数曲线即可。numpy 是一个用python实现的科学计算包,包括一个强大的N维数组对象Array和成熟的函数库,有用于整合C/C++和Fortran代码的工具包,提供了实用的线性代数、傅里叶变换和随机数生成函数等工具,可以理解成Matlab。
记得中学的时候,我问老师三角函数到底有啥用?(无知者无畏)老师反问我,“如果给了你一块洋铁,怎样才能剪出煤炉烟囱的拐弯呢?”,现在仍然记得老师的这个例子,哪些看似抽象的数学公式,实际上是自己不知道她们的应用场景而已。
sin_cos.png
Python代码如下:
import numpy as np
import matplotlib.pyplot as plt
X = np.linspace(0, 2 * np.pi, 100)
YSinValues = np.sin(X)
YCosValues = np.cos(X)
plt.plot(X, YSinValues)
plt.plot(X, YCosValues)
plt.show()
Matplotlib 是一个 Python 的 2D绘图库,甚至可以生成出版质量级别的图形。
对于那些正态分布而言,Python 画起来也就相当简单了:
import numpy as np
import matplotlib.pyplot as plt
def pdf(X, mu, sigma):
a = 1. / (sigma * np.sqrt(2. * np.pi))
b = -1. / (2. * sigma ** 2)
return a * np.exp(b * (X - mu) ** 2)
X = np.linspace(-6, 6, 1000)
for i in range(3):
samples = np.random.standard_normal(10)
mu, sigma = np.mean(samples), np.std(samples)
plt.plot(X, pdf(X, mu, sigma), color = '.66')
plt.plot(X, pdf(X, 0., 1.), color = 'b')
plt.show()
为了不显得单调,这里多画了几条曲线。只要算出方差和均值,从excel中读出哪些数值就可以拟合正态分布了。
normal.png
回归到主题,关于爱心线,有这样一个凄美的爱情故事。
迩来流浪于吴越,一片闲云空皎洁。
300多年前,斯德哥尔摩的街头,落魄的笛卡尔过着乞讨的生活,全部的财产破破烂烂的衣服和随身所带的几本数学书籍。清高的笛卡尔没有开口请求路人施舍,只是默默地低头在纸上写写画画,潜心于他的数学世界。 一个宁静的午后,笛卡尔照例坐在街头的阳光中研究数学问题,身边过往的人群,喧闹的车马队伍,都无法对他造成干扰。
有美一人,婉如清扬。邂逅相遇,与子偕臧。
“你在干什么呢?”扭过头,笛卡尔看到一张年轻秀丽的睑庞,一双清澈的眼睛如湛蓝的湖水,楚楚动人,她就是瑞典的小公主,国王最宠爱的女儿克里斯汀。 她蹲下身,拿过笛卡尔的数学书和草稿纸,和他交谈起来。言谈中,他发现,这个小女孩思维敏捷,对数学有着浓厚的兴趣。
几天后,笛卡尔意外地接到通知,国王聘请他做小公主的数学老师。满心疑惑的笛卡尔跟随侍卫一起来到皇宫,他听到了从远处传来的银铃般的笑声。他看到了那天在街头偶遇的女孩子,从此,他当上了公主的数学老师。
情不知所起,一往而深
公主的数学在笛卡尔的悉心指导下突飞猛进,他们之间也开始变得亲密起来。笛卡尔向她介绍了直角坐标系,代数与几何可以结合起来,也就是日后笛卡尔创立的解析几何学雏形。 在笛卡尔的带领下,克里斯汀走进了奇妙的坐标世界,她对曲线着了迷。每天的形影不离也使他们彼此产生了爱慕之心。
在瑞典这个浪漫的国度里,一段纯粹而美好的爱情悄然萌发。
念去去,千里烟波,暮霭沉沉楚天阔。
他们的恋情传到了国王的耳中,国王大怒,下令将笛卡尔处死。在克里斯汀的苦苦哀求下,国王将他放逐回国,公主被软禁在宫中。
当时,欧洲大陆正在流行黑死病。笛卡尔回到法国后不久,便染上重病。在生命进入倒计时的那段日子,他日夜思念的还是街头偶遇的那张温暖笑脸。他每天坚持给她写信,盼望着她的回音。然而,这些信都被国王拦截下来,公主一直没有收到他的任何消息。
欲知心里事,看取腹中书.
在笛卡尔给克里斯汀寄出第十三封信后,他永远地离开了这个世界。此时,被软禁在宫中的小公主依然徘徊在皇宫的走廊里,思念着远方的情人。
这最后一封信上没有写一句话,只有一个方程:r=a(1-sinθ)。
国王以为这个方程里隐藏着两个人的秘密,便把全城的数学家召集到皇宫,但是没有人能解开这个函数式。他不忍看着心爱的女儿每天闷闷不乐,便把这封信给了她。拿到信的克里斯汀欣喜若狂,她立即明白了恋人的意图,找来纸和笔,把图形画了出来,一颗心形图案出现在眼前,克里斯汀泪流满面,这条曲线就是著名的“心形线”。
国王去世后,克里斯汀继承王位,便立刻派人去法国寻找心上人的下落,收到的却是笛卡尔去世的消息,留下了一个永远的遗憾…… 这封享誉世界的另类情书,据说至今还保存在欧洲笛卡尔的纪念馆里。
这个故事的出处无从知道,网络上流传着各种各样的版本,甚至在百度百科也有着这个故事。后来,有人考证了真实性,认为这是一个美丽的谎言,但并不妨碍人们对爱心线喜爱。
在直角坐标系中,爱心线的方程的python 表达为:x** 2+ y** 2 + a * x= a * sqrt(x** 2+y** 2) 和 x** 2+ y** 2 - a * x= a * sqrt(x** 2+y** 2)
通过x 来求对应的y值很麻烦,就像软件设计中的“万能层”那样,可以采用参数方程来表示:
x=a*(2*cos(t)-cos(2*t))
y=a*(2*sin(t)-sin(2*t))
具体的python代码如下:
import numpy as np
import matplotlib.pyplot as plt
a = 1
t = np.linspace(0 , 2 * np.pi, 1024)
X = a*(2*np.cos(t)-np.cos(2*t))
Y = a*(2*np.sin(t)-np.sin(2*t))
plt.plot(Y, X,color='r')
plt.show()
代表爱心的心形线来了:
xy-heart.png
但这不是六行代码呀?也不是r=a(1-sinθ)呀? 的确如此,那是极坐标系,python 的matplotlib同样支持极坐标系的,爱心线的六行pyton代码如下:
import numpy as np
import matplotlib.pyplot as plt
T = np.linspace(0 , 2 * np.pi, 1024)
plt.axes(polar = True)
plt.plot(T, 1. - np.sin(T),color="r")
plt.show()
这样,得到的就是封面中的图像了:
polar-heart.png
心形线确实是爱心满满,如果融入了忧伤会是怎样呢?
import numpy as np
import matplotlib.pyplot as plt
x = np.linspace(-8 , 8, 1024)
y1 = 0.618*np.abs(x) - 0.8* np.sqrt(64-x**2)
y2 = 0.618*np.abs(x) + 0.8* np.sqrt(64-x**2)
plt.plot(x, y1, color = 'r')
plt.plot(x, y2, color = 'r')
plt.show()
这样就得到了另一个爱心线:
narrow_heart.png
网络上还有关于爱心线的各种漂亮实现,也充满了各种各样的情绪,但对于每一种,基本上都可以用python 相对简洁的实现。
实际上,绘图很简单,难的是那些曲线方程的表达以及实际的应用场景需求,比如螺旋线。
rope.png
进一步,还可以画出各种的3维视图,例如:
3d.png
作为数据分析乃至大数据处理的最后一个环节,就是所谓洞见的可视化,python 可以说是其中的一个简单实用的工具。
参考阅读
http://matplotlib.org/
《Python 数据可视化编程实战》
《Python 数据分析实战》
六行python代码的爱心曲线_6行python代码的爱心曲线相关推荐
- 六行python代码的爱心曲线_6行python代码的爱心线
前些日子在做绩效体系的时候,遇到了一件囧事,居然忘记怎样在Excel上拟合正态分布了,尽管在第二天重新拾起了Excel中那几个常见的函数和图像的做法,还是十分的惭愧.实际上,当时有效偏颇了,忽略了问题 ...
- python程序员的浪漫-千行生日祝福代码
python程序员的浪漫-千行生日祝福代码
- python 代码_6行Python代码的爱心曲线
前些日子在做绩效体系的时候,遇到了一件囧事,居然忘记怎样在Excel上拟合正态分布了,尽管在第二天重新拾起了Excel中那几个常见的函数和图像的做法,还是十分的惭愧.实际上,当时有效偏颇了,忽略了问题 ...
- python9行代码_如何用9行Python代码编写一个简易神经网络
原标题:如何用9行Python代码编写一个简易神经网络 Python部落(python.freelycode.com)组织翻译,禁止转载,欢迎转发. 学习人工智能时,我给自己定了一个目标--用Pyth ...
- python自动生成鸡汤文_20行python代码实现鸡汤智能生成器
"Don't think of the overwhelming majority of the impossible." "不要去想不可能之事" " ...
- python红色的颜色表达式_50行Python代码实现视频中物体颜色识别和跟踪(必须以红色为例)...
目前计算机视觉(CV)与自然语言处理(NLP)及语音识别并列为人工智能三大热点方向,而计算机视觉中的对象检测(objectdetection)应用非常广泛,比如自动驾驶.视频监控.工业质检.医疗诊断等 ...
- python语音分割_用7行Python代码构建自己的有声读物
点击关注我哦 欢迎关注 "小白玩转Python",发现更多 "有趣" 有声读物是我们可以通过音频听取一本书或者其他作品的内容,是现下一种很受欢迎的阅读方式.类似 ...
- python 知乎关系图谱_5000行python代码+可视化60W数据,告诉你知乎用户不为人知的事...
一次完整的python分析+可视化展示,是什么样的? 比如我想知道知乎用户的学历,是否都是985呢?我还想知道知乎最受关注的话题都是些什么?高端人士都喜欢看什么书呢?"人在XX,刚下飞机?& ...
- python必背100源代码-100行Python代码实现一款高精度免费OCR工具
近期Github开源了一款基于Python开发.名为Textshot的截图工具,刚开源不到半个月已经500+Star. 这两天抽空看了一下Textshot的源码,的确是一个值得介绍的项目. 相对于大多 ...
最新文章
- eclipse下classes文件夹无法发布到tomcat的问题--tomcat发布慢的问题
- 【ZOJ - 2972】Hurdles of 110m (dp)
- 关于c#静态构造函数
- criteria函数_干货铺 | 二级MS office考试中一些常考函数(2)
- Font Manager :字体管理及批量安装工具
- 【java】调用System.gc()会发生什么?
- A加载B,B发生的变化
- MATLAB保存数据为dat格式,[转载] Matlab中的数据以.txt或.dat格式保存
- 劳荣枝潜逃 23 年落网,多亏了它!
- 外媒:巴基斯坦将成为南亚地区首个测试5G通讯的国家
- mysql alisql_对比!!mysql5.7 mysql8.0 alisql5.6同样插入400万数据的对比
- 帆软填报报表的使用教程
- 兆骑科创创新创业赛事活动承办,科创企业,企业孵化器
- 计算机应用基础网课作业答案,知到网课答案计算机应用基础(新)全部答案
- 《全民学乒乓》学习笔记
- Python写的Web spider(网络爬虫)
- 人工智能中神经网络与自然语言处理共生关系_CodingPark编程公园
- 浅析正向代理、反向代理、透明代理
- iOS_异常堆栈报告分析
- WeUI:专为开发微信HTML5应用的开源Web UI组件库
热门文章
- 信号与系统第一章--基本知识
- JAVA日记之mybatis-2代理开发/动态sql/sql抽取/核心配置标签 ----喝最烈的酒.
- 《挑战程序设计竞赛》--初级篇习题POJ部分【动态规划】
- c语言写桌面程序unity,Unity和iOS原生界面交互示例
- jdk9、jdk10、jdk11、jdk12、jdk13新特性
- window11安装JDK1.8【jdk-8u121】
- leetcode数组相关简单习题,玉米迪迪的刷题之旅(*╹▽╹*)
- 什么是ICTI认证?ICTI认证有什么好处?
- 响铃:国际竞争变标准之争,长袖善舞的海尔这次又走在了前面
- itunes在计算机中丢失,iTunes用户莫名丢失音乐文件 苹果将发布更新