绘制引力波

最新推荐文章于 2023-11-12 15:45:03 发布

一腔诗意喂了猫

最新推荐文章于 2023-11-12 15:45:03 发布

阅读量505

点赞数

CC 4.0 BY-SA版权

分类专栏：数据分析文章标签： Python plot

本文链接：https://blog.youkuaiyun.com/qq_36142336/article/details/78807288

数据分析专栏收录该内容

9 篇文章

订阅专栏

本文通过使用Python中的NumPy和Matplotlib等库，读取并绘制了来自两个不同探测器（H1和L1）的引力波数据。同时，还展示了预定义的引力波模板信号。通过对这些信号进行可视化，可以直观地比较不同来源引力波的特点。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

（本文是北理工嵩天老师MOOK课程笔记之一）

# -*- coding:utf-8 -*-
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()
hi=1/rate_h
li=1/rate_l

htime_len=hstrain.shape[0]/rate_h
htime=np.arange(-htime_len/2,htime_len/2,hi)
ltime_len=lstrain.shape[0]/rate_l
ltime=np.arange(-ltime_len/2,ltime_len/2,li)
fig=plt.figure(figsize=(12,6))

plth=fig.add_subplot(221)
plth.plot(htime,hstrain,'y')

pltl=fig.add_subplot(222)
pltl.plot(ltime,lstrain,'g')

pltt=fig.add_subplot(212)
pltt.plot(reftime,ref_H1)

fig.tight_layout()
plt.show()
plt.close(fig)