Перестраиваемый полосовой фильтр Опции

[РРоман]

Имеется следующая задача:
создать полосовой фильтр с заданной полосой пропускания (~100Гц), при этом необходимо иметь возможность "перемещать" его полосу пропускания заданном диапазоне частот (1Гц-3КГц).
Подскажите пожалуйста дилетанту можно ли вообще создать такой фильтр в цифровой форме? Как?
Можно ли для этого использовать ДПФ? Если да то как (в общих чертах).

Спасибо!

[Guest_TSerg_*]

создать полосовой фильтр с заданной полосой пропускания (~100Гц), при этом необходимо иметь возможность "перемещать" его полосу пропускания заданном диапазоне частот (1Гц-3КГц).
Подскажите пожалуйста дилетанту можно ли вообще создать такой фильтр в цифровой форме? Как?
Можно ли для этого использовать ДПФ? Если да то как (в общих чертах).

Почему бы не взять книжку по цифровым фильтрам, да и не почитать ?
Как я понимаю, у Вас учебная задача

Если Вы научитесь считать ФНЧ любым понятным Вам способом, то фильтры остальных типов легко приводятся к ФНЧ.
Например ПФ это два ФНЧ с требуемыми частотами среза и два сумматора.

Если хочется что-то поинтереснее, то см:
http://www.tsi.lv/RSR/vol04/RSR2000v04-07ae.pdf

[РРоман]

Почему бы не взять книжку по цифровым фильтрам, да и не почитать ?
Как я понимаю, у Вас учебная задача

Если Вы научитесь считать ФНЧ любым понятным Вам способом, то фильтры остальных типов легко приводятся к ФНЧ.
Например ПФ это два ФНЧ с требуемыми частотами среза и два сумматора.

Если хочется что-то поинтереснее, то см:
http://www.tsi.lv/RSR/vol04/RSR2000v04-07ae.pdf

Кое-что по фильтрам читал, но на решение поставленной задачи меня не хватает... Задача-то не учебная, а вполне рабочая.

За ссылочку большое спасибо! Попробуем... Только по-моему алгоритм Герцеля попроще считается, за то описанный в работе вариант, наверное, более гибкий.

[Guest_TSerg_*]

Задача-то не учебная, а вполне рабочая.
Только по-моему алгоритм Герцеля попроще считается, за то описанный в работе вариант, наверное, более гибкий.

Тогда надо точнее описывать задачу и точнее будет ответ.

Вообще-то, алгоритм Герцеля, как разновидность DFT, занимается преобразованием свойств сигнала из временной области в частотную. Понятно также, что алгоритм блочный.
Заниматься его временными откликами - странная затея.

[РРоман]

Вообще-то, алгоритм Герцеля, как разновидность DFT, занимается преобразованием свойств сигнала из временной области в частотную. Понятно также, что алгоритм блочный.
Заниматься его временными откликами - странная затея.

Так, вроде, все по книжке (Сергиенко А.П. "ЦОС"): "ДПФ можно представить набором фильтров, а дискретную фильтрацию можно организовать с помощью ДПФ", там еще и АЧХ двух частотных каналов ДПФ приведены...

[Guest_TSerg_*]

Так, вроде, все по книжке (Сергиенко А.П. "ЦОС"): "ДПФ можно представить набором фильтров, а дискретную фильтрацию можно организовать с помощью ДПФ", там еще и АЧХ двух частотных каналов ДПФ приведены...

"Фильтр" Герцеля даст магнитуду по конкретной анализируемой частоте после обработки блока данных.
Более того, при совпадении частоты настройки "фильтра" и частоты входного сигнала "фильтр" становится неустойчивым, что мы и видели на Вашей картинке, а также см:

Если же входная частота ушла, то "фильтр" имеет затухающий временной отклик.

Но дело в том, что именно временной отклик в этом алгоритме мало кого интересует.

[РРоман]

"Фильтр" Герцеля даст магнитуду по конкретной анализируемой частоте после обработки блока данных.
Более того, при совпадении частоты настройки "фильтра" и частоты входного сигнала "фильтр" становится неустойчивым, что мы и видели на Вашей картинке

Я реализовал алгоритм не полностью, а только ту его часть, которая вычисляется на каждом отсчете:

Код

s(n) = x(n) + 2cos(2πω)s(n − 1) − s(n − 2)

и решил, что это и будет мой фильтр :). Мощность сигнала (magnitude?) я вообще не считал. Правильно я поступал или нет? Можно ли вообще реализовать ВЕРМЕННУЮ фильтрацию с помощью ДПФ или Герцеля???

А какая неравномерность в полосе пропускания, и какое подавление за ее пределами?

Конкретные цифры пока не известны. Фильтр нужен для звука, параметры буду подбирать экспериментально.

Ну кажется в маткаде тоже самое записывается как i=1..длина_массива_h.

тогда не понимаю что означает запись:
hf = h.*sin(2*pi*f/Fs*(1:length(h))
нужно рассчитать для всех элементов массива? Тогда, какого массива?

[alex_os]

тогда не понимаю что означает запись:
hf = h.*sin(2*pi*f/Fs*(1:length(h))
нужно рассчитать для всех элементов массива? Тогда, какого массива?

Все просто делается FIR ФНЧ, его коэф - массив h, размером length(h). каждый элемент массива поэлеметно (.*) перемножаете на синусоиду нужной вам частоты. в результате получается полосовой FIR с коэф. hf, центральная частота полосового фильтра f. ФНЧ h проектируете любым доступным способом исходя из нужных Вам полосы пропускания, затухания и пр. Получается табличка h, sin произвольной частоты легко посчитать. Все это, в большинстве случаев работает.
ps и забудьте про Герцеля .