Обработать результаты измерения Опции

[SergSit]

Честно говоря незнаю в какой раздел форума писать. Если не по адресу поправте. Задача стоит следующая: обработать результаты измерения(которые снимаются с частотой 100Гц) на предмет удаления из измерений влияние медленно меняющихся процесов (изменения температура и влажность среды скажем за время большее 3 сек). Размерность результатов измерения 8 или 16бит. Только просьба не отправляйте к чтению литературы по DSP. Меня интересует не теория , а практика построения такого медленного фильтра. Из того , что успел почитать понял, что мне вроде нужен БИХ фильтр ФНЧ 1-го порядка. Как это практичеки реализовать в программе не нашел. То, что попадалось касается построения более быстрых фильтров. Может кто поделиться ссылками, где обьяснено как это сделать. Заранее сильно не пинайте, делаю это первый раз и могу не понимать суть проблемы. Реализовать надо на AVR.

[xemul]

Занятный подход. Нужно удалить зависимость от температуры и влажности, но вместо этого опять хочется пофильтровать...
Не находите, что методически верным было бы определить зависимости измеряемой величины от влияющих факторов, а в процессе - измерять эти факторы, корректировать соответственно результат измерения измеряемой величины, а уже потом, при желании, фильтровать.
Т.е. набираете сначала достаточное для обработки количество экспериментальных данных (измеряемая величина (измеренная более другим средством измерения более высокого класса точности) + результат измерения Вашим измереметром + измерения влияющих факторов) в интересующих диапазонах изменения измеряемой и влияющих величин, строите модель, связывающую результаты Вашего измереметра с измеряемой величиной с учетом влияющих факторов, воплощаете эту модель в AVR.

[Палыч]

Занятный подход.

Всё зависит от того, что получаем при обработке данных... Например, при обработке ЭКГ для получения частоты сердечных сокращений (к примеру) не важно как изменяется (вследствие движения человека) средняя линия сигнала. Измерить напряжения мышц человека - затруднительно, да и не нужно - легче отфильтровать...

[xemul]

Всё зависит от того, что получаем при обработке данных... Например, при обработке ЭКГ для получения частоты сердечных сокращений (к примеру) не важно как изменяется (вследствие движения человека) средняя линия сигнала. Измерить напряжения мышц человека - затруднительно, да и не нужно - легче отфильтровать...

Я не специалист по обработке ЭКГ, но Вы описАли простейший случай, когда возможно разделить спектры измеряемой величины и влияющих факторов. И, насколько мне известно, обработка ЭКГ производится не по абсолютной амплитуде сигналов и в этом отношении носит скорее качественный, чем количественный характер.
Например, если нас интересует только изменение температуры, измеряемой термопарой, то мы можем просто стабилизировать температуру ее холодных концов. Если же нас интересует абсолютное значение температуры, то температуру этих концов придется измерять.
Автор пока умалчивает о интересующем диапазоне частот измеряемой величины, предоставляя нам возможность всласть потелепатировать.

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

Т.е. организовать два (кольцевых) буфера на 3-5-... значений Вы уже не сможете из-за хардварных ограничений?

Это все понимаю. У меня стоит задача другая: усреднить и отфильтровать.

Хотелось где нибудь именно о цифровой фильтрации медленно меняющихся процессах и использование их на практике.(именно на практика)

Т.е. Вы уверены, что фильтрации будет достаточно?
Порядок фильтра (и количество потребных для фильтрации отсчетов) будет определяться требуемым подавлением в полосе задержания, которое можно рассчитать, зная коэффициенты, связывающие результат измерения и влияющие факторы, и требуемую точность измерения.
Без теории будет тяжко.
А быстро или медленно... В этом мире, благодаря т. Эйнштейну, все относительно.

[SergSit]

Автор пока умалчивает о интересующем диапазоне частот измеряемой величины, предоставляя нам возможность всласть потелепатировать.

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

[Палыч]

В моем случае о диапазоне частот измеряемой величины говорить тяжело, т.к. измеряемая величина от отсчета к отсчету остается постоянной (за исключением шумов скажем в пределах 1%).

Как же Вы собираетесь фильтровать, если не известен диапазон частот полезного сигнала?

За какое время T величина должна измениться на 10%, чтобы был зарегистрирован факт изменения? Теперь найдите значение 1/T, и среднепотолочно можно будет считать это за требуемую полосу снизу.

Это - даже не среднепотолочно, а - "не знаю что". Имхо, так и порядок величины нижней частоты сигнала не корректно определять.

[xemul]

Это - даже не среднепотолочно, а - "не знаю что". Имхо, так и порядок величины нижней частоты сигнала не корректно определять.

Я на языке автора попытался объяснить (правда, предполагая ответ, - следующим предложением, которое Вы не процитировали), как все-таки можно _оценить_ потребную полосу снизу.
А на счет корректно или нет... По скорости нарастания сигнала Вы сможете определить частоту среза фильтра, чтобы сигнал искажался не более чем на ...? Именно это со среднепотолочными допусками я и предложил.

Согласен может я не понимаю сути проблемы, нехватает знаний. Рассуждаю так. Мне надо убрать из измеренного потока, медленную ошибку. Т.е. мне нужно построить ФВЧ с частотой среза соответствующей постоянной времени 3сек. В чем не прав поправте. Быстрые изменения фильтр буде пропускать, медленные( с какой-то точностью) не пропускать.

Я и спросил, насколько быстрые изменения сигнала Вы хотите регистрировать.

[Палыч]

Я и спросил, насколько быстрые изменения сигнала Вы хотите регистрировать.

С постановкой задачи у автора - проблема. Но, как я понял, это - что-то похожее на термостабилизацию, когда нужно предпринимать некие действия при выходе температуры за границы диапазона. Наверное, температура меняется не быстро при неизменности внешних факторов. Более важен вопрос: "Насколько медленные изменения сигнала Вы хотите регистрировать?" и "Корректно ли отбрасывать медленные изменения температуры при неизменности внешних факторов?"