并试图将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 ++中的循环移位来实现这两个功能,这是我已经完成的功能。

非常感谢。

最佳答案

看完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/

关于c++ - fftshift/ifftshift的换档,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/34918807/

10-11 21:05