И ещё вопрос про фильтрацию сигнала, Помогите отфильтровать сигнал новичку Опции

[Sergey Semenov]

Добрый день,

Есть сигнал :



Частота дискретизации: Fs= 40000

Сигнал представляет из себя сумму 1Гц тона и полезного акустического сигнала phi(t), что-то типа такого

F(t) = A + B * t + sin(2 * pi * 1 * t) + phi(t)

А и B - константы

Вопрос: как можно отфильтровать исходный сигнал, чтобы выделить акустический сигнал phi(t), который имеет частоту от 20Гц?
В дальнейшем этот сигнал хочется отфильтровать, дицимировать и проанализировать в диапазоне от 20Гц до 2 Кгц.

High-pass-filter c cut-off frequency = 5Hz получился у меня со слишком высоким порядком и будет сложен для реализации. Notch filter получился неудачно, пока не понял почему. В каком направлении вообще нужно идти?

Спасибо!

В приложении сам сигнал длительность 2.8 сек
 matlab.mat.zip ( 404.32 килобайт ) Кол-во скачиваний: 132

[ViKo]

Вот что получил ВЧ-фильтром Баттерворта 8 порядка на 200 Гц. Когда пробую понизить частоту, увеличить порядок, фильтр ломается, видимо, из-за переполнений.

Эскизы прикрепленных изображений

 

[Sergey Semenov]

Вот что получил ВЧ-фильтром Баттерворта 8 порядка на 200 Гц. Когда пробую понизить частоту, увеличить порядок, фильтр ломается, видимо, из-за переполнений.

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

Спасибо всем за ответы! попробую пойти по пути, предложенному thermit.