c++ - 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/