FFT

卷积

卷积定理

循环卷积

如果B是一个block-circulant 矩阵，则：
if this convolution is periodic, A(PSF) is then a (block) circulant matrix

fPSF_decon= (fft2(ifftshift(PSF)));

如果滤波器在时间/空间中居中，则需要在进行FFT之前进行ifftshift，以便去除在时域中居中的线性相位2.
见知乎

1

用以上方法计算卷积时，若两个信号长度相差很多，则较短者须补上相当多的零，太不经济。而且在某些情况下，例如较短的 h[n] 是一个 FIR 滤波器而较长的 x[n] 是未知长度的输入（像语音）时，直接用以上方法要等所有的输入都收到后才能开始算输出信号，太不方便。这时可以把 x[n] 分割成许多适当长度的区块（称为 block convolution），然后一段一段的处理。经过滤波后的段落再仔细的连接起来，借由输入或输出的重叠来处理区块连接的部分。这两种做法分别称为重叠-储存之卷积法和重叠-相加之卷积法。

1. https://www.wikiwand.com/zh-cn/%E5%9C%93%E5%91%A8%E6%91%BA%E7%A9%8D
2. https://dsp.stackexchange.com/questions/47249/result-of-inverse-fft-is-sometimes-shifted-in-real-space
3. https://blog.csdn.net/shenxiaolu1984/article/details/50884830