我想先计算出我使用以下代码的加速度数据的hpf,然后进行FFT。
# Make Function
def filter(s):
# Signal after hpf
fc = 0.1
b = 0.08
N = int(np.ceil((4 / b)))
if not N % 2: N += 1
n = np.arange(N)
sinc_func = np.sinc(2 * fc * (n - (N - 1) / 2.))
window = np.blackman(N)
sinc_func = sinc_func * window
sinc_func = sinc_func / np.sum(sinc_func)
# reverse function
sinc_func = -sinc_func
sinc_func[int((N - 1) / 2)] += 1
#s = list(df[' acc.y'])
new_signal = np.convolve(s,sinc_func)
a=np.fft.fft(new_signal)
return a
问题是这现在返回一个数组,我想计算此Numpy的单值幅度,即我需要FFT我计算出的信号的幅度