matplotlib绘制引力波

详细注释

'''
引力波绘制
'''
import numpy as np
import matplotlib.pyplot as plt
from scipy.io import wavfile#读取波形的库
'''
从配置文档中读取时间相关数据
r表示转义，当出现特殊字符时加一个r表示原来的字符
H1_Strain.wav、L1_Strain.wav:音频文件
wf_template.txt：引力波的理论模型
np.genfromtxt():执行两个循环：将文件的每一行转化为字符串序列和将每一个字符串序列转化为相应的数据类型；即是一个两行的矩阵
transpose():矩阵转置
'''
rate_h,hstrain=wavfile.read(r"H1_Strain.wav","rb") #rate:速率 strain：数据矩阵
rate_l,lstrain=wavfile.read(r"L1_Strain.wav","rb")
reftime,ref_H1=np.genfromtxt(r"wf_template.txt").transpose()#reftime:时间序列 ref_H1:信号的数据'''
求导数，得到波形的时间间隔
'''
htime_interval=1/rate_h
ltime_interval=1/rate_h'''
丢失信号起始点
'''
htime_len=hstrain.shape[0]/rate_h
#hstrain.shape[0]表示读取矩阵的第一维的长度即数据点的个数，除以坐标轴rate就可以得到函数在坐标轴上的总长度
htime=np.arange(-htime_len/2,htime_len/2,htime_interval)
#绘制出关于原点对称的图像，以-1/2为起点1/2，interval为时间间隔 创造时间序列htime
ltime_len=lstrain.shape[0]/rate_l
ltime=np.arange(-ltime_len/2,ltime_len/2,ltime_interval)'''
绘制图形y
绘制以时间为x轴，应变数据为Y轴的图像，并设置标题和坐标轴的标签
'''
fig=plt.figure(figsize=(12,6))#创绘图空间12*6plth=fig.add_subplot(221)
plth.plot(htime,hstrain,'y')
plth.set_xlabel('Tiem (seconds)')
plth.set_ylabel('H1 Strain')
plth.set_title('H1 Strain')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("Wave.png")
plt.show()
plt.close(fig)

matplotlib绘制引力波相关推荐

嵩天老师Python面向对象-53 Matplotlib绘制引力波
python绘制引力波
使用python绘制引力波经过不懈的努力人们终拿到了黑洞的照片,但是是否还记得之前人们人们发现的引力波吗?想不想自己手绘一个引力波的图形呢? 接下来就手把手教你绘制一个引力波的图像. 绘制引力波最重 ...
matplotlib绘制图表，设置刻度标签、最大最小刻度、字体大小，label位置、刻度轴箭头等
matplotlib绘制图表,设置刻度标签.最大最小刻度.字体大小,label位置.刻度轴箭头等 1. 效果图 2. 源码 2.1 仅使用普通轴ax + fontdict 源码 2.2 使用mpl设置 ...
使用KMeanCluster对多个区域进行聚类，并结合Matplotlib绘制中心点、最大最小距离点
使用KMeanCluster对多个区域进行聚类,并结合Matplotlib绘制中心点.最大最小距离点 1. 效果图 2. 源码 2.1 原始数据--xq.txt 2.2 源码参考这篇博客将演示如何 ...
Matplotlib绘制象限图——支持中文、箭头、自定义轴线交点
Matplotlib绘制象限图--支持中文.箭头.自定义轴线交点 1. 效果图 2. 原理 2.1 绘制象限图步骤 2.1 添加文字到图表 3. 源码参考这篇博客将介绍如何使用matplotlib ...
python绘制三维轨迹_Python学习（一） —— matplotlib绘制三维轨迹图
在研究SLAM时常常需要对其输出的位姿进行复现以检测算法效果,在ubuntu系统中使用Python可以很好的完成相关的工作. 一. Ubuntu下Python的使用在Ubuntu下使用Python有 ...
python使用matplotlib绘制一条正弦曲线（plot函数可视化sine plot）
python使用matplotlib绘制一条正弦曲线(plot函数可视化sine plot) 目录 python使用matplotlib绘制一条正弦曲线(plot函数可视化sine plot) #导入 ...
Python使用matplotlib绘制透明背景的可视化图像并保存透明背景的可视化结果（transparent background）
Python使用matplotlib绘制透明背景的可视化图像并保存透明背景的可视化结果(transparent background) 目录
Python使用matplotlib绘制分组对比柱状图（bar plot）可视化时汉语（中文）标签显示成了框框□□、什么情况、我们有解决方案
Python使用matplotlib绘制分组对比柱状图可视化时(bar plot)汉语(中文)标签显示成了框框□□.什么情况.我们有解决方案目录

matplotlib绘制引力波

matplotlib绘制引力波相关推荐

最新文章

热门文章