Как удалить ненужную прямую линию на участке?

Question

Данные в массиве freq ниже не отсортированы.

Есть ли удобный способ удалить ненужную прямую линию слева направо на моем участке?

import pylab as pl 
import numpy as np 
timepoints=np.loadtxt("timesequence",usecols=(1,),unpack=True,delimiter=",") 
t=np.histogram(timepoints,bins=500)[0] 
sp = np.fft.fft(t) 
freq = np.fft.fftfreq(t.shape[-1],d=0.0005) 
print freq 
pl.plot(freq*2*np.pi, np.sqrt(sp.real**2+sp.imag**2)) 
pl.show() 

Answer 1

Как вы его, сюжет начинается в точке пика ноль, затем работает свой путь направо, затем переходит в крайнее левое положение и работает свой путь назад к центру. Это НЕ простой вариант времени (справа налево), если это так, как вы думаете.

Одним из способов является построение положительных точек в «частоте» отдельно от отрицательных точек в «freq». Замените строку кода pl.plot следующим образом:

mask = freq>=0 
pl.plot(freq[mask]*2*np.pi, np.sqrt(sp[mask].real**2+sp[mask].imag**2)) 
pl.plot(freq[~mask]*2*np.pi, np.sqrt(sp[~mask].real**2+sp[~mask].imag**2)) 

ps, ​​Вам необходимо установить цвета, чтобы они соответствовали друг другу.