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

[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

[amaora]

Идентифицировать постоянную соствляющую и фазу 1Гц сигнала, в остатке будет phi.

[Sergey Semenov]

Идентифицировать постоянную соствляющую и фазу 1Гц сигнала, в остатке будет phi.

А как? Можете сказать что почитать, куда посмотреть?

Я, к несчастью, очень неопытен в ЦОС.

Как я понимаю, если делать FFT с небольшим числом точек (512) - тогда разрешение будет порядка 40 гц и
этот 1Гц сигнал вылезет на других частотах. А позволить делать FFT с большим числом точек тоже
будет затруднително, я боюсь на финальной платформе не хватит ресурсов