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

[Relayer]

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

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

[GetSmart]

Т.е., буквально, фильтруем отдельно действительную часть сигнала, затем комплексную, складываем и получаем комплексный результат.

По одной части сигнала кажется невозможно отличить частоты +F и -F (симметрия от 0) т.к. они совпадают с противоположной, сдвинутой на 180 град. Так что фильтр должен работать сразу с комплексным представлением. Либо предворительно сдвинуть спектр в части, подлежащей фильтрации.

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

Из чистого любопытства, для чего понадобилась такая необычная (ИМХО) фильтрация?

Сообщение отредактировал GetSmart - Apr 14 2008, 09:45

[Relayer]

По одной части сигнала кажется невозможно отличить частоты +F и -F (симметрия от 0) т.к. они совпадают с противоположной, сдвинутой на 180 град. Так что фильтр должен работать сразу с комплексным представлением. Либо предворительно сдвинуть спектр в части, подлежащей фильтрации.

о чем я и говорил

Из чистого любопытства, для чего понадобилась такая необычная (ИМХО) фильтрация?

вот для этого: http://code.google.com/p/sdrlab

Не знаю, поможет или нет. Но мы для расчета комплексного фильтра (чтобы не было симметрии относительно 0) делали так.
1. Вычисляли обычный полосовой фильтр необходимой конфигурации в матлабе. При этом получался фильтр с действительными коэффициентами, который симметричен относительно 0.
2. Делали преобразование Гильберта над фильтром. При этом одна половинка исчезает. .

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

PS 2 fontp и blackfin - господа, вы наверное не совсем поняли - я тут никому ничего доказывать не собираюсь. высшее математическое есть, но по цос у меня самообразование. и по работе я специализируюсь в областях далеких от электроники и цос. а ветку открыл именно с целью чтобы профессионалы не зубоскалили и выпендривались друг перед другом, а помогли реальным советом