Отправлено ВН 09 декабря 2002 г. 21:17 В ответ на: 2-х мерное преобразование Фурье или Уолша отправлено BEAVIS 09 декабря 2002 г. 19:25

Прямое Фурье:
S(l,k)=(1/(M*N))*SUMMA(SUMMA(s(i,j)*exp(-j*2*pi((l*i/M)+(k*j/N)))))
Внутренняя сумма по j, от 0 до N-1.
Внешняя по i, от 0 до M-1.
Или наоборот.
Обратное:
s(i,j)=SUMMA(SUMMA(S(l,k)*exp(j*2*pi((l*i/M)+(k*j/N)))))
Внутренняя сумма по k, от 0 до N-1.
Внешняя по l, от 0 до M-1.
Или наоборот.
Преобразования обратимы. Скорее всего у Вас ошибка в реализации.
Простой вариант реализации.
1. Добиваете 0 по каждому из измерений так, чтобы N и M стали ближайшими степенями 2. Для определенности: i-строки, j-столбцы.
2. Каждую строку пропускаете через одномерное прямое FFT размером N.
3. Над результатом пункта 2 - каждый столбец через одномерное прямое FFT размером M. Получаете спектр.
Обратное, чтобы не путаться, в обратном порядке, т.е. сначала п.3, затем 2. Одномерные FFT, естественно, обратные.
Итого для прямого: M FFT размером N + N FFT размером M.
Столько же для обратного.
Для упрощения входные данные, s(i,j), если они действительные, лучше считать комплексными, с 0 мнимой частью.

• Ответ: Забыл добавить. — ВН   (09.12.2002 21:51, 901 байт)
  • Спасибо, попробую еще раз... — BEAVIS   (10.12.2002 11:40, пустое)