并试图将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/