быстрый FIR для комплексного сигнала Опции

[Relayer]

Приветствую!

в качестве предисловия - занимаюсь опенсорсным SDR софтом. возник следующий проблемс. имеем квадратурные каналы I/Q. применяем к ним комплексное ффт. хотим через fast convolution сделать band pass filter. а фигушки коэффициенты FIR"а рассчитанные через sync+окно вещественные - результат свертки симметричен относительно 0Hz. в принципе я "выкрутился" - построил low pass коэфф и комплексно домножил на сдвиг sin-cos. работает, фильтрует, но чует мое сердце что где-то есть общая метода рассчета комплексных коэфф FIR-фильтров для быстрой свертки. подскажите кто чем может

[GetSmart]

Пост №4

Только зачем это Вам? Вы же фильтрацию в частотной области делаете, или я не понял чего-то? Для этого достаточно домножить результат комплексного ДПФ на вещественную функцию "АЧХ" фильтра, и сделать ОДПФ.

И что получим? Симметрию относительно 0. Где там польза для автора?

Пост №14

Ему бы ещё теорему о циклической свёртке изучить не мешало.

Свёртке чего (в данном контексте) с чем? Вы её знаете? И что же такое "полезное" эта теорема говорит?

Сообщение отредактировал GetSmart - Apr 16 2008, 06:51

[Relayer]

И что получим? Симметрию относительно 0. Где там польза для автора?

вобщем проблему более менее разрешил. точнее сказать убедился в том, что выбранный мной изначально метод - правильный и рабочий. ачх, точнее коэфф на которые надо домножать в частотной области у меня получились комплексные. вчера дописал код и вчерновую запустил - артефактов нет.
всетаки остановился на bandpass/stop. low/highpass тоже работает, но с некоторыми оговорками - об этом позже. итак у нас есть Flo, FHi. строим windowed-sinc filter. lowpass если требуется bandpass ачх и highpass - если bandstop. частота среза этого фильтра = (Fhi-Flo)/2. получаем вещественные коэфф фильтра. их N/2 штук, при этом FFT работает с окном размером N. старшую половину блока зануляем. причина связана с тем что используется техника overlap&add. теперь генерируем комплексную экспоненту с частотой (Flo+Fhi)/2 и поэлементно домножаем (комплексно) на нее сгенерированные коэфф фильтра. от полученного комплексного массива выполняем прямое FFT и получаем коэфф на которые будет домножаться (комплексно) сигнал в частотной области.
теперь касательно low/highpass в комплексной плоскости. их тоже можно сгенерировать таким образом. к примеру если нам необходим lowpass с частотой среза Flo можно сгенерировать по рассмотренной выше методике bandpass c полосой -Fs/2..Flo. но вылезла мелкая бяка - на результирующей ачх есть небольшой пик в районе F/s (хотя там должно быть подавление). мое предположение - это происходит из за свойств цикличности.