DFT и IDFT программно для звукового сигнала, преобразовать туда и обратно Опции

[sergiator]

Здравствуйте!
Проблема: с помощью преобразования Фурье преобразую аудио сигнал (речь, файл in.wav в архиве) с временной в частотную форму а потом в обратно во временную.
В результирующем сигнале речь можно разобрать, но присутствует много посторонних звуков (файл out.wav в архиве). Понимаю, что перед каждым преобразованием (или после) нужно применять какие-то фильтры (по-моему, даже называются они оконные).
Объясните популярно, пожалуйста: какие фильтры (название) и в каком месте применять (до или после).

ЗЫ: Скорость преобразования не имеет значения - нужно просто разобраться что к чему. Поэтому использую не FFT, а DFT, код переделал от сюда (обратное) и от сюда (прямое)

Сообщение отредактировал sergiator - Sep 20 2013, 11:05

Прикрепленные файлы

 example.zip ( 208.97 килобайт ) Кол-во скачиваний: 12

 

[sergiator]

Здравствуйте!
Проблема: с помощью преобразования Фурье преобразую аудио сигнал (речь, файл in.wav в архиве) с временной в частотную форму а потом в обратно во временную.
В результирующем сигнале речь можно разобрать, но присутствует много посторонних звуков (файл out.wav в архиве). Понимаю, что перед каждым преобразованием (или после) нужно применять какие-то фильтры (по-моему, даже называются они оконные).
Объясните популярно, пожалуйста: какие фильтры (название) и в каком месте применять (до или после).

ЗЫ: Скорость преобразования не имеет значения - нужно просто разобраться что к чему. Поэтому использую не FFT, а DFT, код переделал от сюда (обратное) и от сюда (прямое)

Так, я натупил, но проблема осталась:
1) натупил: неправильно преобразовал данные из byte pcm (формат входящего файла) в double (необходимого для FFT) и обратно
2) проблема: все равно есть посторонние звуки (файл out2.wav)

Прикрепленные файлы

 out2.wav ( 149.67 килобайт ) Кол-во скачиваний: 6

 

[ukpyr]

попробуйте умножить вх.сигнал на какое-нибудь окно:
http://ru.wikipedia.org/wiki/%D0%9E%D0%BA%...%80%D1%8C%D0%B5

[andyp]

Вам бы про свертку в частотной области (convolution in frequency domain) почитать.

[Corner]

Для начала надо понять, что Фурье в лоб для звука не пригодно. Надо преобразовывать с перекрытием. Если берете, для примера, 8192 точки то перекрытие желательно 6144 точки (сдвиг на 2048). И при обратном преобразовании результат надо усреднять по числу перекрывающихся областей (в данном случае по 4-ем).

Сообщение отредактировал Corner - Sep 20 2013, 19:11