使用python绘制引力波

经过不懈的努力人们终拿到了黑洞的照片，但是是否还记得之前人们人们发现的引力波吗？想不想自己手绘一个引力波的图形呢？

接下来就手把手教你绘制一个引力波的图像。

绘制引力波最重要的是要有数据源，可以从下面的网址下载我们所需要的数据

数据源

从这个链接，把下面的这些文件下载下来。

需要使用的python库

numpy
matplotlib
scipu.io

接下来就是代码的编写了，新建一个文件，命名为wave.py，然后使用下面的代码

    import numpy as npimport matplotlib.pyplot as pltfrom scipy.io import wavfilerate_h, hstrain= wavfile.read(r"H1_Strain.wav","rb")rate_l, lstrain= wavfile.read(r"L1_Strain.wav","rb")#reftime, ref_H1 = np.genfromtxt('GW150914_4_NR_waveform_template.txt').transpose()reftime, ref_H1 = np.genfromtxt('wf_template.txt').transpose() #使用python123.io下载文件htime_interval = 1/rate_hltime_interval = 1/rate_lfig = plt.figure(figsize=(12, 6))# 丢失信号起始点htime_len = hstrain.shape[0]/rate_hhtime = np.arange(-htime_len/2, htime_len/2 , htime_interval)plth = fig.add_subplot(221)plth.plot(htime, hstrain, 'y')plth.set_xlabel('Time (seconds)')plth.set_ylabel('H1 Strain')plth.set_title('H1 Strain')ltime_len = lstrain.shape[0]/rate_lltime = np.arange(-ltime_len/2, ltime_len/2 , ltime_interval)pltl = fig.add_subplot(222)pltl.plot(ltime, lstrain, 'g')pltl.set_xlabel('Time (seconds)')pltl.set_ylabel('L1 Strain')pltl.set_title('L1 Strain')pltref = fig.add_subplot(212)pltref.plot(reftime, ref_H1)pltref.set_xlabel('Time (seconds)')pltref.set_ylabel('Template Strain')pltref.set_title('Template')fig.tight_layout()plt.savefig("Gravitational_Waves_Original.png")plt.show()plt.close(fig)

欢迎关注公众号，获取更多有趣信息

python绘制引力波相关推荐

怎么用python找论文_如何利用Python绘制学术论文图表方法
论文中图表形式多样,常用的处理工具有excel.MATLAB以及Python等,excel自处理的方法有两个缺陷: 1.当数据较多时,容易出现excel"翻白眼"的现象: 2.需要 ...
python绘制条形图用什么函数_Python绘制正余弦函数图像完整代码
通过python绘制正弦和余弦函数,从默认的设置开始,一步一步地调整改进,让它变得好看,变成我们初高中学习过的图象那样.通过这个过程来学习如何进行对图表的一些元素的进行调整. 01. 简单绘图 mat ...
python绘制所有特征的密度图（density plot）
python绘制所有特征的密度图(density plot) # python绘制所有特征的密度图: n = len(input_feats) #密度图kde,查查分布 cols = 4 rows = ...
python绘制影像组学训练集、测试集对应的ROC曲线以及瀑布图（rad-score 瀑布图）
python绘制影像组学训练集.测试集对应的ROC曲线以及瀑布图(rad-score 瀑布图) # 所有数据的瀑布图 tagets_all = df['label'] tagets_all = df. ...
什么是折线图？怎样用Python绘制？怎么用？终于有人讲明白了（附代码）
来源:大数据DT 本文约5400字,建议阅读10分钟本文为你介绍数据分析时经常用到的折线图,可以用来呈现哪些数据关系?在数据分析过程中可以解决哪些问题?怎样用Python绘制折线图? 01 概述折 ...
柱状图、堆叠柱状图、瀑布图有什么区别？怎样用Python绘制？（附代码）
来源:大数据DT(ID:hzdashuju) 作者:屈希峰,资深Python工程师,知乎多个专栏作者本文约8000字,建议阅读20分钟柱状图是当前应用最广泛的图表之一,你几乎每天都可以在电子产品上 ...
层次聚类python_用python绘制层次聚类图
层次聚类(Hierarchical clustering)代表着一类的聚类算法,这种类别的算法通过不断的合并或者分割内置聚类来构建最终聚类. 聚类的层次可以被表示成树(或者树形图(dendrogram ...
利用Python绘制 3D 体素色温图
简介: 测试了Matplotlib中的体素绘制函数 voxels() ,但是无论是在 AI Studio中执行,还是直接在本地执行,都会出了错误.也就是无法通过 gca(projection='3 ...
python绘制直方图显示数字_Python实现绘制双柱状图并显示数值功能示例
本文实例讲述了Python实现绘制双柱状图并显示数值功能.分享给大家供大家参考,具体如下: # -*- coding:utf-8 -*- #! python3 import matplotlib.py ...

python绘制引力波

使用python绘制引力波

欢迎关注公众号，获取更多有趣信息

python绘制引力波相关推荐

最新文章

热门文章