Интегрирование показаний датчика на С, как правильно? Опции

[smk]

Мне нужно проинтегрировать показания датчика. Сами показания измеряются АЦП контроллера с дискретностью в 2 миллисекунды. Как я понимаю это должно быть численное интегрирование по методу трапеций. Начальное значение (обычно "а" в формулах) равно 0, как я понимаю. Количество проделанных измерений я могу посчитать. Это, как я понимаю будет "b" и "N" одновременно. Соответственно "h" получится всегда равно 1. А значит само интегрирование сведется к подсчету среднеарифметических значений между соседними замерами и их последующее суммирование? Наверняка я что-то понимаю не так. Если кто уже делал подобное, подскажите как правильно. Код я наверняка сам напишу, мне бы быть уверенным что алгоритм верный. Спасибо.

[Tanya]

Мне нужно проинтегрировать показания датчика.
Если кто уже делал подобное, подскажите как правильно. Код я наверняка сам напишу, мне бы быть уверенным что алгоритм верный. Спасибо.

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

Лично мне за много лет работы с данными на практике хватает двух методов: усреднение и медианный фильтр. Но это может зависеть от данных, какова специфика датчика и какие помехи нужно побороть.

У Вас много лет не было практики со спектральными приборами или хроматографами там разными...
Это не помехи бороть, а весь полезный сигнал собрать.

[smk]

Это не помехи бороть, а весь полезный сигнал собрать.

Совершенно верно. Т.е. я так понимаю, что нужно суммировать предидущее и текущее значения, делить их на 2, а результат накапливать в переменной в виде суммы всех вышеупомянутых среднеарифметических?

Кстати, попутно вопрос. Как с помощью printf вывести unsigned long (напр. 12123123)? Или лучше объявить как float? Тогда выводится правильно.

[ReAl]

Как с помощью printf вывести unsigned long (напр. 12123123)?

Код

#include <inttypes.h> // само включит <stdint.h>
#include <stdio.h>

int main()
{
        uint64_t foo = 12123123ULL;
        printf( "foo = %" PRId64 "\n", foo);
        return 0;
}

Это даст правильный результат где угодно.

[smk]

Это даст правильный результат где угодно.

Спасибо.

По интегрированию еще бы так...

[Tanya]

По интегрированию еще бы так...

А блюдечко с голубой каемочкой приготовили?
Вы же обещали сами написать...
В чем суть вопроса Вашего?
Численное интегрирование не проблема. Проблема в том, чтобы правильно выделить пик, который нужно интегрировать. В общем случае...
Вот пока еще нет пика есть некий фон - какая-то неизвестная кривуля. Или прямуля.
Ее нужно правильно выделить. Второе - нужно отрезать хвост. Иначе наинтегрируем до бесконечности. И нос нужно правильно найти.
Есть еще небольшие тонкости связанные с конечной разрядностью и накоплением ошибки вычислений.