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

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

Ну все правильно Вы сделали. Взяли ФНЧ и перенесли на нужную комплексную частоту. Чего тут еще придумаешь?

Синфазная и квадратурная составляющие ничем по частотным свойствам друг от друга не отличаются. Поэтому их можно пропускать через одинаковые фильтры. Или, как рекомендуется в литературе, нужно разработать фильтр с семплом Fs, а I и Q пропускать через него поочередно с семплом Fs/2 (они же сдвинуты на 90 град).
Если нужна теория, почитай на сайте Ксайлинкса даташиты:
http://www.xilinx.com/ipcenter/catalog/log...docs/da_fir.pdf
http://www.xilinx.com/support/documentatio...piler_ds534.pdf
http://www.xilinx.com/ipcenter/catalog/log...ocs/mac_fir.pdf

Здесь много теоретической информации (даже если с плисинами не работаешь)

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

Это совершенно из другой оперы.
Здесь речь идёт о фильтрации через ДПФ.

Сообщение отредактировал Stanislav - Apr 11 2008, 07:24

ну ... эээ ... как бы хотелось теории и вских вумных книжек на эту тему с bandpass и bandstop это все получается. а вот low/high pass да в комплексной таким макаром уже не получится. в принципе оно и не особо нужно, но так сказать "для общего развития" не помешало бы


вы еще предложите значения для ачх фильтра брать только из "+1" и "0"

Честно говоря не очень понятно о чём Вы рассказывете как на Привозе, приходится догадываться
Если я правильно догадался, то вот есть такая фигня - синтез всяких фильтров посредством переноса спектра. Здесь есть Application Note
Multi-rate filter design

у них была ещё программа для синтеза таких крутых фильтров, называется mrf. Так это AN для той программы
Как водится программка есть в сети и крякнута

за ссылку спасибо - вечерком почитаю. а рассказываю я вот о чем - к примеру хочу я спроектировать фильтр который пропускает все отрицательные частоты и положительные предположим до F. т.е. lowpass но в комплексной области через быструю свертку.

Что Вы называете быстрой свёрткой? Быстрые свёртки бывают разные. Посредством преобразование Фурье?


непонятна постановка задачи. Если у Вас фильтрация идёт через ДПФ и известна частотная характеристика - то причём здесь FIR. Умножил спектры -и всё. Или Вам нужно получить заданную импульсную характеристику?

Ладно, забудьте :-) В том материале, что я привёл строятся эффективные полосовые фильтры в исходном домене по следующей методике 1) перенос спектра в 0 2) фильтрация посредством фильтра НЧ
3) возвращают спект на место на заданную частоту

Вычислительный эффект получают передискретизацией перед фильтрацией НЧ в соответствии с полосой фильтра

ну естественно. через него родимого - я же написал об этом в начале темы



йя-йя. танцуем от windowed sinc-filter



ну я собственно все это и сделал ... в частотной области. lowpass/highpass сдвинутый на нужное смещение по частоте. как результат для комплексного сигнала такая схема реализует bandpass/bandstop фильтр для произвольных Flo и Fhi, в том числе и отрицательны

Не понял юмора.
Вы это о чём, уважаемый? Прошу пояснить.


По-моему, он просто не знает, как нужно фильтровать сигнал в частотной области.
А также как умножить вещественное число на комплексное.

я это к эффекту Гиббса если вы не сразу соображаете



хочу вам напомнить уважаемый вы наш всезнайка и писатель (интересно - найдется ли такая тема где вы не продемострировали свою эрудицию) - для личной переписки есть личные сообщения. впредь прошу вас воздержаться от писанины в этой теме. а особенно не по существу. если вам не терпится выяснить со мной отношения - мы можем их продолжить выяснять в той самой ветке где вы публично продемострировали ваше полное незнание самых элементарных принципов компресии

Полегче на поворотах, уважаемый.
Вероятно, для Вас будет открытием, что явление Гиббса присутствует в любом КИХ-фильтре, включая и тот, что Вы пытаетесь реализовать. Однако, это не имеет отношения к теме, а мой вопрос заключается в следующем: отчего это Вы решили, что я буду предлагать Вам сделать к-ты АЧХ фильтра, равными только 0 и 1? Вот на него мне и хочется прочитать ответ.

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

Тогда ещё вопрос: где Вы у меня видели "писанину", и что написанное в теме Вашим покорным слугой не правильно?

Форум открыт для общения; не Вам запрещать мне высказывать собственное мнение, которое я, к тому же, могу подкрепить строгим доказательством.
А вот Вам не мешало бы почитать собственные посты, и попытаться со стороны расшифровать, о чём же хотел поведать миру их автор?


Это я-то?
Обязательно продолжим.

В общем-то Вы правильно использовали одно из свойств (частотного сдвига) преобразования Фурье.

http://ru.wikipedia.org/wiki/%D0%9F%D1%80%....80.D1.8C.D0.B5

Умные книжки - любой учебник по функциональному анализу или по ЦОС (Рабинер-Голд, Оппенгейм-Шафер, отечественные авторы Гольденберг-Матюшкин и др.)

можно здесь посмотреть:
http://dsp-book.narod.ru/books.html


Успехов.

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

ну в этом то я не сомневался



к сожалению по большей части везде описывается dsp-обработка вещественных сигналов. а у меня комплексный - в этом то собственно и проблема.