Отправлено ВН 09 декабря 2002 г. 21:51 В ответ на: Ответ: отправлено ВН 09 декабря 2002 г. 21:17

На примере квадратного двумерного FFT, M=N.
После прямого двумерного FFT 0 частоты по обоим измерениям находится в левом верхнем углу. Часто это неудобно. Делают перестановку результирующей матрицы так, чтобы 0 частоты был в центре. Делается так. Результирующая матрица на 4 квадранта разбивается. Первый квадрант: l и k от 0 до 0.5*N-1. 2-й: l от 0.5*N до N-1, k от 0 до 0.5*N-1. 3-й: l от 0 до 0.5*N-1, k от 0.5*N до N-1.
4-й:l и k от 0.5*N до *N-1.
Далее меняются местами 4 и 1, 3 и 2.
В результате 0 частоты в центре, вниз от него положительные "столбцовые" частоты, вверх - отрицательные "столбцовые". Вправо - положительные "строчные", влево -отрицательные.
Перед обратным Фурье - обратная перестановка.
Тот же результат для прямого можно получить, если умножить исходные данные s(i,j) на exp(-j*pi*i)*exp(-j*pi*j). Индексы неудачно выбрал, первая j - мнимая единица.

• Спасибо, попробую еще раз... — BEAVIS   (10.12.2002 11:40, пустое)