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

[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

[Pavia]

Внису немного смуты в головы сдешних обитателей.

Sergey Semenov
Как-то странно, что вы боитесь фильтра 8 порядка? Так как это очень низкий порядок для ЦОС. Или вы его в аналоге хотите делать?
А во вторых я не понядл слово "проще" ?

Напоминаю что процесс передескритизации с понижении частоты состоит из 2-х пунктов.
1. Вначале фильтрация.
2. Потом децимация

Если пропустить фильтрацию, то появяться ложные частоты. А.Б.Сергиенко_Цифровая обработка сигналов.djvu стр 231

Отсюда вывод от передескритизаци вы выиграете только при наличии дальнейшей обработки.

amaora
Алгоритм Герцеля имеет недостатки.
1. Узкая полоса. Что годиться только для идеальной синусойды. А в реалии она может плавать +-1Гц. А вот для реальной лучше делать фильтр с полосой в несколько Гц.
2. Для выхода на режим требует несколько периудов скажем так 10, а то точность будет маленькая.

Sergey Semenov
Что касаеться оптимизации, то надо смотреть разные фильтры:
линейный фильтр с КИХ характеристикой, Баттерворта, RC, рекуретный,линейный фильтр с БИХ характеристикой, эллептический.

ViKo
То что у вас преведено в 9 сообщении в ЦОС подобное можно сделать фильтром 2 порядка. Вернее есть фильтр скользящего среднего который реализуеться просто при наличии RAM. Просто запоминаем в суматоре.

r_2=1/(2*r); // константа
цикл по i
S[i]=S*r_2;
S=S+A[i+r]-A[i-r];
i=i+1;
конец цикла
Где r зависит от частоты среза.

А вообще я тоже только начинающий в ЦОС.