@Владислав56

Т.З. Записан сердечный ритм, файл в wave формате, нужно оставить сигнал в полосе частот от 2 до 50 Гц, использовать программные фильтры
Мое решение в лоб, читаем один аудиокадр, т.к. звук стерео, то делим на левый и правый канал, затем применяем БПФ к каждому из них, получаем частоты, собираем это все в массив, фильтруем до заданной полосы частот.
Подскажите, правильно ли я рассуждаю и как потом из полученных частот собрать новый wave файл?

0

@VTsaregorodtsev

Если бы на кадры Вы сигнал не делили - то обратным преобразованием Фурье. Но при покадровой обработке такой подход приведёт к негладкой сшивке на границах кадров.
Поэтому, видимо, лучше взять обычный цифровой КИХ-фильтр. Но там частоты рядом с границей полосы пропускания всё равно будут пролезать - ибо крутизна спада фильтра будет 10-20-30-... (в зависимости от порядка фильтра) децибел на октаву.
Но можно попробовать делать Фурье на отрезках длиной в 3 кадра (добавлять к текущему кадру предыдущий и следующий), обратное фурье - тоже на трёх этих кадрах, но сшивать затем только СРЕДНИЕ части=кадры. Негладкости на границах должны сильно уменьшиться.


Как вариант: поскольку до программирования пока ещё не дошли - можно попросить тутошнего модератора перекинуть тему в раздел https://www.cyberforum.ru/digi... rocessing/
Там спецы должны быть получше, чем я (я в DSP и/или электронике - чисто по верхам).
А затем с уточнённым решением вернётесь сюда, если возникнут вопросы именно по программированию=реализации данного там совета.


Ну и мой вопрос к Вашему вопросу: зачем отфильтрованный сигнал собирать обратно в файл - его кто-то будет ушами слушать?
Т.е., может, задача у Вас не совсем удачно (или не совсем окончательно) поставлена?
Может, для принятия решения (для некоторого кадра = для некоторого момента времени) проге таки хватит анализа энергии сигнала в выделенной полосе частот - а человек-контролёр пусть при необходимости слушает исходный=нефильтрованный файл?

0

Рыжий Лис

Глупый вопрос, но что мешает открыть файл в любом аудио-редакторе и отсечь ненужные частоты?

0

Рыжий Лис

0

@Владислав56

Вы правы, по краям будут сильные всплески, но можно сделать первоначальную обработку это кадра (блэкмен), в начале уменьшить громкость и постепенно плавно увеличивать, а потом ближе к концу кадра - уменьшать, но все равно получается "шум" в спектре из-за постоянного увеличения и уменьшения громкости.
А если вообще не делить на кадры и полностью закинуть data из wave под fft, при это разделить на каналы, потом сделать конкатенацию при этом ещё как-нибудь нормализовать, отфильтровать получившийся массив, те частоты, которые не проходят - занулить, потом сделать обратное преобразование Фурье, так сработает?
Можно дальше не собирать в аудиофайл, это я уже решил так сделать, чтобы услышать результат фильтрации.

0

Рыжий Лис

Возьмите фильтр чебушева: https://docs.scipy.org/doc/sci... freqs.html

0