Есть устройство на tms320f2810. там реализовано БПФ для вещественного сигнала.
Возникла необходимость в реализации преобразования гильберта там же.
из [url="http://www.dsplib.ru/content/hilbert/hilbert.html"] ясно, что для выполнения сего для начала делается прямое БПФ, затем обнуляется отрицательная часть спектра и удваивается положительная, после чего выполняется обратное БПФ и на выходе получаем синфазную и квадратурную составляющую.

Проблема в следующем: данный алгоритм работает для комплексного сигнала. То есть если имеем 1024 отсчета сигнала, то на вход необходимо подать 2048 отсчетов, где для комплексной составляющей пишем нули. Хотелось бы сэкономить память и попытаться использовать вещественный БПФ для преобразования гильберта, когда на вход подавались бы те же 1024 отсчета.

Преобразователь Гильберта делается как КИХ фильтр, и только потом уже КИХ фильтр вычисляется с помощью быстрой свёртки на основе FFT, если нужно.

Если стадии FFT используют входной буффер (что, называется in-place), то сэкономить не удастся. Уже после первой стадии бабочек получите комплексный сигнал.

Разумеется, можно сэономить написав свое FFT c хитрой первой стадией и не in-place

PS Вы статью по Вашей ссылке до конца дочитали? Там в конце приводится схема быстрых вычислений с FFT половинного размера.

Сообщение отредактировал andyp - Apr 15 2013, 08:13

Дочитать до конца дочитал но там на ускорение алгоритма... Оказывается еще нужно захватывать полный период сигнала чтоб минимизировать ошибку преобразования. Применение окна не позволило добиться приемлемых результатов. подумываю о ких фильтрах

С этого надо было начинать, если частоты не слишком близко к нулю и половине частоты дискретизации, то возможно достаточно будет десятка умножителей в этом фильтре.