我想计算我使用以下代码的加速度数据的 First 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 的信号的幅度