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]

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

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

非常感谢。

Answer 1

在查看Matlab代码后，它不直接使用循环移位，而是使用Matlab语法。

说N =不。元素

实施fftshift，

circularShiftRightBy = floor(N/2)

实施ifftshift，

circularShiftRightBy = ceil(N/2)

如果N为奇数，则为N / 2，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/

在转移方面的fftshift / ifftshift

1 个答案: