2

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

4

0 回答 0