0

由于从 arctan 函数返回的角度上升到正值，我在尝试解决我的 omega 与时间图时遇到了很多问题。我尝试添加 2*pi 使其返回到第一象限，但无济于事（从答案到使用 arctan / arctan2 绘制从 0 到 2π 的 a）。然后我尝试将数组中的所有项目设为负数以便它继续，但我不确定问题是什么： arctan plot 返回不连续跳跃的值 我也尝试使用 arctan2、math.atan、math.atan2

h = h5py.File("Waveforms/BBH0001/rhOverM_Asymptotic_GeometricUnits.h5", "r")
data = h['Extrapolated_N2.dir/Y_l2_m1.dat']

由于在尝试查找数组的反正切时遇到错误，我不得不定义一个函数，然后将 iy 向量化

def f(z):
    return np.int(math.degrees(np.arctan(z)))

f2 = np.vectorize(f)
z = data[:,2]/data[:,1]
a = f2(z)

这部分代码对你们中的一些人来说可能看起来很痛苦。我尝试使用很多不同的方法来解决问题

i = 1
for i in xrange(len(z)):
    if a[i] > (a[i-1] + 10):
        b = [-x for x in a]

plt.plot(data[:,0], b, label='w')
4

0 回答 0