Python:一维数组循环卷积

Question

我想知道 numpy/scipy 中是否有用于一维数组循环卷积的函数。 scipy.signal.convolve() function仅提供“模式”但不提供“边界”，而 signal.convolve2d() 函数需要二维数组作为输入。

作为时间序列作业的一部分，我需要这样做来比较开放卷积和循环卷积。

Answer 1

根据卷积定理，可以用傅立叶变换得到圆卷积。

import numpy as np
def conv_circ( signal, ker ):
    '''
        signal: real 1D array
        ker: real 1D array
        signal and ker must have same shape
    '''
    return np.real(np.fft.ifft( np.fft.fft(signal)*np.fft.fft(ker) ))