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

Оценивать мат. ожидание процесса.

А постоянная составляющая по какому-то закону изменяется (в зависимости от какого-то внешнего параметра, температуры, например) или непредсказуемо?

Закон не известен. Можно с некоторыми оговорками предположить, что постоянная составляющая не изменяется. Но это грубовато.


Это понятно. Как это сделать, если на входе только кусочки синусоиды, которые постепенно формируют период. А оценивать постоянную составляющую нужно сразу.

Вообще-то я делал так:
- усреднял полученный буфер и рассматривал его как первое приближение.
- усреднял следующий буфер - и в качестве оценки постоянной составляющей брал среднее от двух предыдущих оценок - следующее приближение и т.д.
Такая оценка в итоге сама превращалась в синусоиду...

По идее, для получения постоянной составляющей надо устремить время усреднения к бесконечности (в разумных пределах, конечно). Я бы использовал достаточно большое плывущее окно, по которому бы и усреднял.
Если девайсу Вашему памяти нехватает для большого плывущего окна - пусть компьютерная программа этим занимается.

Можно попробовать сделать цифровой НЧ фильтр, со срезом на низкой частоте (заведомо меньшей частоты Вашего полезного сигнала)

Мат.ожидание (среднее в частном случае) это далеко не "постоянная" составляющая процесса. Достаточно усреднить синусоиду в районе pi/2 чтобы это увидеть. Если оценку постоянной составляющей надо делать быстро, по каждому буферу, то можно попробовать метод МНК. Или ФНЧ как предлагали выше.

Большое окно - это хорошо. Расчетом занимается и так комп, просто ему и так задач хватает.
Буду думать как сделать расчет по окну более эффективным. А можно сделать размер окна адаптивным?

Фильтр НЧ с постоянной времени заведомо во много периодов синусоиды поможет.

Простейший рекурсивный y(n) = (1-e)* y(n-1) + e*x(n) , если всё равно какой в смысле времени реакции, е очень мало

Но если это слишком медленно, то нужно сшивать буфера, оценивать период синусоиды и суммировать по периоду

Зделать FFT и 0 отщёт будет = пост состовляющей.

1

Сообщение отредактировал ivan219 - Nov 21 2008, 19:58

Круто!
Средство на все случаи жизни.


Глубокая мысль. А как это сделать практически, не подскажете?


Некошерно. В лоб решение, без мысли...

А как сшивать, ежели постоянка?


М-да, лекарство не лучше болезни... Условиям задачи не удовлетворяет. А, между тем, решение существует.

Для вычисления параметров синусоиды и определения постоянной составляющей теоретически достаточно 4-х отсчётов сигнала.
Решением системы 4-х уравнений будут амплитуда, частота, фаза синусоиды и постоянка.
В Вашей постановке с учётом дополнения задача имеет точное решение по содержимому только одного буфера.
Возможность практической реализации ограничивается необходимой точностью при заданных параметрах выборки и заданном уровне помех.
Оценки ПС в окнах выборки можно усреднять (с "забыванием", если подставка медленно меняется).

Криво, уж простите... Ничего хорошего так не выйдет.

А чем не подходит это
http://electronix.ru/forum/index.php?act=S...f=2&t=54351

Всем спасибо!
Лично я остановился на варианте уважаемого fontp, как наиболее приемлемом алгоритмически для моей задачи.

Сообщение отредактировал vladimir_m - Nov 22 2008, 12:32

А каким должен быть коэффициент K Вашего чудо-алгоритма?

Ну, можно, конечно, и так, если мозгами шевелить лень...
Только коряво это. Оценку долго ждать придётся, если хоть мало-мальская точность нужна..

Несколько странно... Такой метод подойдет ведь для сигналов любой формы и их суммы. Если Вы точно знаете, что там одна чистая синусоида, то можно эффективнее... Хозяин - барин.