telegra.ph/Chart-01-10 добрый день, по ссылке находися изображение графика, который получается при считывании данных с ацп, это затуающие колебания, правда график обрезан с одной стороны( проблема аппаратная), необходимо найти частоту колебаний. Уже месяц бьюсь над проблемой, получалось найти частоту для согнала, колебающегося относительно постоянного цента, но тут сигнал не симметричный. Может есть готовые методы решения подобных задач ?
Серьезно? А вручную по графику частоту и период находить не учили?
Кстати говоря, сигнал у вас судя по всему совсем не шумный, так что если его легонечко фильтрануть, легко убрать постоянную составляющую можно с помощью производной - частота сохраняется, только фаза сдвигается. Но если есть малейшие шумы, то получается мусор. Лучше конечно по хорошему сделать - полосно-пропускающий фильтр
Славка: ну тогда ведь даже постоянную составляющую не нужно отсеивать, ищете два локальных максимума в пределах 1/5 -- 1/20 секунды и находите временной период между ними, тем самым получая период. А 1/период=частота. Максимумы советую искать следующим образом, берете отрезок размером 1/5 сек., находите минимум и по обе стороны от минимума находите максимумы, причем поиск максимумов нужно продолжать, пока или производная не изменит свой знак(при том сама она должна изменяться немного инерционно, для устойчивости к ошибкам) или пока вы не дойдете до точки, всего на 0,2*(max-min) отличающейся от первого минимума - это для случая, когда не нужно избавляться от постоянной составляющей. Очень грубо с точки зрения алгоритмов, но никаких Фурье-преобразований, очень просто и достаточно надежно.
Славка: Настолько сильные? Тогда предварительно сгладить. Или искать минимумы/максимумы первой производной.
Вам ведь абсолютные значения не нужны, только количество локальных минимумов или максимумов.
Ну и в начале графика вроде сигнал чистый, так что частоту можно по первой трети определить.
Корректный вариант: предобразование Фурье, естественно перед этим избавиться от "обрезания", далее оттуда получаем основную гармонику.
Малокорректный вариант (в свое время видел в практическом использовании у сейсмологов): находить некую "ось" ("ноль") и считать полупериодом соседние пересечения оси.
"ось"\"ноль" - видимо аппроксимировать полиномом первого порядка из среднеквадратичного.
