circshift 方面的 fftshift / ifftshift

Question

并试图将 fftshift/ifftshift 与循环移位联系起来。

N = 5 
Y = 0:N-1

X = [0 1 2 3 4]

当我 fftshift(X) 时，我得到

[3 4 0 1 2]

当我 ifftshift(X) 时，我得到

[2 3 4 0 1]

如何将 fftshift/ifftshift 与循环移位联系起来？它只是简单地向不同方向移动 X 中的数字吗？

我需要知道这一点，因为我正在尝试根据 C++ 中的循环移位实现这两个函数，这是我已经完成的函数。

非常感谢。

Answer 1

看了Matlab代码，没有直接使用循环移位，而是使用了Matlab语法。

说 N = 否。元素数

要实施 fftshift，

circularShiftRightBy = floor(N/2)

要实现 ifftshift，

circularShiftRightBy = ceil(N/2)

是N/2，如果N是奇数，fftshift和ifftshift之间只有区别。

其中循环移位码为：

template<typename ty>
void circshift(ty *out, const ty *in, int xdim, int ydim, int xshift, int yshift)
{
 for (int i =0; i < xdim; i++) {
   int ii = (i + xshift) % xdim;
   if (ii<0) ii = xdim + ii;
   for (int j = 0; j < ydim; j++) {
     int jj = (j + yshift) % ydim;
     if (jj<0) jj = ydim + jj;
     out[ii * ydim + jj] = in[i * ydim + j];
   }
 }
}

（从 fftshift/ifftshift C/C++ source code 修改为也支持左移 (-ve)。）

编辑：后来我找到了更好的方法：https://kerpanic.wordpress.com/2016/04/08/more-efficient-ifftshift-fftshift-in-c/

circshift 方面的 fftshift / ifftshift

fftshift / ifftshift in terms of circshift

c++

matlab

fft

ifft