Отправлено SM 25 февраля 2004 г. 12:33 В ответ на: Только ifft(abs(fft(signal))^2) - модуль забыл. отправлено SM 24 февраля 2004 г. 16:06

это я писал про автокорреляцию. А кросс-корреляция через FFT выглядит вот так:

out = ifft(fft(X) .* conj(fft(Y)));

Где X и Y могут быть как добиты нулями, так и нет, в зависимости, интересует ли циркулярность при корреляции. Ну ссответственно так как Y есть образцовый сигнал, то conj(fft(Y)) считается заранее и хранится в памяти.

• Спасибо, перечитал Оппенгейма и понял что overlap-add корреляция через FFT это то что надо, хочу коррелировать по тысяче точек, а в тайм домене больше 120-180 не получается — AntZ   (26.02.2004 03:19, пустое)
  • Еще (+) — SМ   (26.02.2004 09:36, 278 байт)