Поиск частоты сигнала?

Question

[Славка]

telegra.ph/Chart-01-10 добрый день, по ссылке находися изображение графика, который получается при считывании данных с ацп, это затуающие колебания, правда график обрезан с одной стороны( проблема аппаратная), необходимо найти частоту колебаний. Уже месяц бьюсь над проблемой, получалось найти частоту для согнала, колебающегося относительно постоянного цента, но тут сигнал не симметричный. Может есть готовые методы решения подобных задач ?

Answer 1

[Андрей]

Пробовали ли частотный анализ (преобразование Фурье) ?

Comments

[Славка]

БПФ, не пробовал, все отговаривали, говоря что сигнал может не удовлетворять требованиям для БПФ, сам я в этом мало что понимаю

[pixik]

Славка: БПФ!

[Славка]

pixik: стоит попробовать? Может тогда подскажите реализацию для шарпа

Answer 2

[pixik]

Реализации алгоритмов/Быстрое преобразование Фурье

Comments

[Славка]

А что из себя представляют входные данные? Как мне из моих отсчетов получить комплексные числа?

[pixik]

Славка: честно говоря, не помню. я это делал давно. Попробуйте действительную часть установить значением амплитуды, а мнимую нулём.

Answer 3

[x67]

Серьезно? А вручную по графику частоту и период находить не учили?
Кстати говоря, сигнал у вас судя по всему совсем не шумный, так что если его легонечко фильтрануть, легко убрать постоянную составляющую можно с помощью производной - частота сохраняется, только фаза сдвигается. Но если есть малейшие шумы, то получается мусор. Лучше конечно по хорошему сделать - полосно-пропускающий фильтр

Comments

[Славка]

это уже пропущенные данные через фильтр калмана, но будет реализован фильтр низких частот, нужная мне частота находится в пределах 11-13 ГЦ

[x67]

Славка: ну тогда ведь даже постоянную составляющую не нужно отсеивать, ищете два локальных максимума в пределах 1/5 -- 1/20 секунды и находите временной период между ними, тем самым получая период. А 1/период=частота. Максимумы советую искать следующим образом, берете отрезок размером 1/5 сек., находите минимум и по обе стороны от минимума находите максимумы, причем поиск максимумов нужно продолжать, пока или производная не изменит свой знак(при том сама она должна изменяться немного инерционно, для устойчивости к ошибкам) или пока вы не дойдете до точки, всего на 0,2*(max-min) отличающейся от первого минимума - это для случая, когда не нужно избавляться от постоянной составляющей. Очень грубо с точки зрения алгоритмов, но никаких Фурье-преобразований, очень просто и достаточно надежно.

Answer 4

[Rsa97]

Если есть сами исходные данные графика, то в чём проблема? Найти локальные минимумы и максимумы и посчитать частоту.

Comments

[Славка]

Шумы

[Rsa97]

Славка: Настолько сильные? Тогда предварительно сгладить. Или искать минимумы/максимумы первой производной.
Вам ведь абсолютные значения не нужны, только количество локальных минимумов или максимумов.
Ну и в начале графика вроде сигнал чистый, так что частоту можно по первой трети определить.

Answer 5

[d-stream]

Корректный вариант: предобразование Фурье, естественно перед этим избавиться от "обрезания", далее оттуда получаем основную гармонику.

Малокорректный вариант (в свое время видел в практическом использовании у сейсмологов): находить некую "ось" ("ноль") и считать полупериодом соседние пересечения оси.

"ось"\"ноль" - видимо аппроксимировать полиномом первого порядка из среднеквадратичного.