引力波（matplotlib绘制）

照着老师的代码敲了一遍，但出了一些错误，直接复制也是同样的错误，找了一会儿才知道在第9，10行的1应该改为1.这样算出来的值才是浮点数，而不是0了。

import numpy as np
import matplotlib.pyplot as plt
from scipy.io import wavfile

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("wf_template.txt").transpose()

htime_interval = 1./rate_h
ltime_interval = 1./rate_l
fig = plt.figure(figsize=(12,6))

#丢失信号起始点
htime_len = hstrain.shape[0]/rate_h
htime = 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_l
ltime = 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)