См. посты №4 и №14.

Пост №4

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

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

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

вобщем проблему более менее разрешил. точнее сказать убедился в том, что выбранный мной изначально метод - правильный и рабочий. ачх, точнее коэфф на которые надо домножать в частотной области у меня получились комплексные. вчера дописал код и вчерновую запустил - артефактов нет.
всетаки остановился на 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 (хотя там должно быть подавление). мое предположение - это происходит из за свойств цикличности.

Вообще-то, я хотел упомянуть о "теореме запаздывания", но уж больно автор темы, гм... того. Теорема о свёртке же говорит о том, как надо фильтровать сигнал в частотной области.

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

вот именно - наиспражнялись вы здесь вволю. ну да ладно - вас ждет тема по lossless audio compression. тема TTA-кодека вами совсем не раскрыта. пора бы с ней покончить - за три дня что вы отмалчиваетесь вполне можно было разобраться в алгоритме. или вы язык ЦЭ в школе не проходили?

Relayer, я в терминах немного не бум-бум. У фильтра (bandpass/stop) получилась крутизна с обоих сторон одинаковая, а симметрия вместо нуля сдвинута на другую частоту?

Преобразование Гильберта не помогло што ли? Мне вот интересно, если восстановить Im часть коэффициентов фильтра по Re части, потом ещё умножить на комплексную экспоненту (или на что другое), то можно убрать из фильтрации скажем все отрицательные частоты?

Сообщение отредактировал GetSmart - Apr 17 2008, 12:34

вот как выглядит ачх - на входе белый шум (комплексный)
bandpass 3-6kHz:

bandstop 3-6kHz:

lowpass 3kHz с теми артефактами справа о которых я писал ранее


крутизна по бокам одинаковая. с гильбертом что-то у меня не пошло. толи не то крутил, толи не туда

Если последний рисунок "свернуть в трубочку", то вполне идеальный фильтр без артефактов. ИМХО единственный способ сделать идеальный low/highpass фильтр в положительных и отрицательных частотах - это "условная" свёртка, о которой я писал раньше.

естественно. артефакт это левый скат ачх который "завернулся" вправо. в зависимости от выбора окна будет более или менее выражен.



гм ... давайте подробнее об этом. "теоретики" поразбежались - есть шанс разобраться в этом вопросе

Тока сначала мне надо понять что имеется на входе. Как я где-то по ссылке прочитал - это комплексное представление ПЧ радиотракта. Непонятно только чем отфильтрованное и какой с полосой. Если в этой полосе будут только частоты -Fs/2..0..Fs/2, то без проблем на выход пропускать все отрицательные, а положительные фильтровать через FIR. Или наоборот (в любых комбинациях).

Хотя принципиальной выгоды в этом может и не быть. Т.к. в комплексном представлении при оцифровке все частоты в входном сигнале как бы сворачиваются "в трубочку" и вылезают в полосе -Fs/2..0..Fs/2. Например -Fs/2 - 1Hz вылезет в полосе Fs/2 - 1 Hz. ИМХО чтобы таких артефактов не было с "реальными" неидеальными аналоговыми фильтрами перед АЦП наверно нужно оцифровывать полосу сигнала заметно шире, а потом использовать только часть (серединку).