绘制引力波
代码
import numpy as np#科学计算
import matplotlib.pyplot as plt # 绘图用的库
from scipy.io import wavfile #读取下载好的音频文件
#数据源:
#http://python123.io/dv/H1_Strain.wav
#http://python123.io/dv/L1_Strain.wav
#http://python123.io/dv/wf_template.txt
#.wav为声音文件
#此处的写法默认为与代码文件保存到同一个地方了
#wavefile.read用于读取下载好的音频文件
#读取的数据分别赋值给read_h和hstrain,字符串前面加r表示原始字符串
rate_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()