fftshift / ifftshift с точки зрения кругового смещения
и я пытаюсь связать fftshift / ifftshift с круговым сдвигом.
N = 5
Y = 0:N-1
X = [0 1 2 3 4]
Когда я жду (X), я получаю
[3 4 0 1 2]
Когда я сдвигаю (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, существует только разница между fftshift и ifftshift, если N нечетно.
Где код кругового сдвига:
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 ++ для поддержки левого (-ve) переключения. )
РЕДАКТИРОВАТЬ: С тех пор я нашел лучший способ сделать это: https://kerpanic.wordpress.com/2016/04/08/more-efficient-ifftshift-fftshift-in-c/
Другие решения
Других решений пока нет …