Про "сшивку", суммирование с перекрытием и "фильтрацию по уму" в случае FFT. Сначала выясняете длину импульсной характеристики M. Если имеете только частотную характеристику, делаете обратное ДПФ от нее и выясняете, при каком M отсчеты становятся меньше требуемой точности вычислений. Потом определеяте потребную длину FFT N. Входной поток разбиваете на блоки по N-M отсчетов, дополняете M нулями и делаете FFT. При необходимости фильтруете и делаете обратное FFT. Таким образом, входные блоки имеют длину N-M, выходные N. Непрерывный выходной поток получаете, объединяя выходные блоки с перекрытием M отсчетов. Внутри зоны перекрытия суммируете одновременные отсчеты из разных блоков. От более раннего блока это затухающие отсчеты, от более позднего - начальные отсчеты переходного процесса. В сумме они обеспечивают математически корректную сшивку блоков. Все, непрерывный выходной поток не имеет переходных процессов внутри себя, только в начале. Хотя, может, ваша задача и не требует обратного преобразования и суммирования с перекрытием.
Сообщение отредактировал V_G - Mar 24 2011, 04:38
|