Быстрая свертка! Как?, помогите, пожалуйста, с алгоритмом Опции

[coolibin]

Нужен алгоритм быстрой свертки(желательно подробный). Если кто то опишет или даст линк, я буду только благодарен

[DRUID3]

Нужен алгоритм быстрой свертки(желательно подробный). Если кто то опишет или даст линк, я буду только благодарен

Свертка во-временнОй области соответствует умножению спектров базиса Фурье в частотной. Спектры быстро вычисляюЦЦо через FFT. Вот Вам и метод.

Собственно и "гугл" рулит: http://alglib.sources.ru/fft/fastconvolution.php

[rloc]

Нужен алгоритм быстрой свертки(желательно подробный). Если кто то опишет или даст линк, я буду только благодарен

Свертка во-временнОй области соответствует умножению спектров базиса Фурье в частотной. Спектры быстро вычисляюЦЦо через FFT. Вот Вам и метод.

Позвольте поинтересоваться, а какой длительности сигнал нужно сворачивать? Может со сжатием в частотной области поспешили? Точно знаю, что фирма Intel в своих библиотеках MKL перед сверткой определяет длительность и в зависимости от этого уже сворачивает. Исходные данные даны не полностью, а задачку уже решаете.

[shasik]

В Matlab'e приведенный выше пример полностью прокатывает, поэтому еще раз обращаю внимание на вычисление прямого и обратного БПФ.
А еще смущает строка:

Код

for(i = 0; i != uSizeA; i++)

Если будете делать фильтрацию в лоб и с помощью БПФ, а затем посмотрите графики выходных отфильтрованных сигналов, то у Вас может наблюдаться небольшой сдвиг (в зависимости от того, как Вы поступили с начальными условиями - считали, что до этого сигнал был равен 0 либо подождали, пока не прибудет необходимое количество отсчетов). И если для сравнения результатов Вы вычисляли СКО сигналов - результаты будут не совпадать