1.10.3 Интегрирование данных

Интегрирование сигналов реализуется рекурсивными цифровыми фильтрами. Рассмотрим примеры анализа интегрирующих операторов.

Алгоритм интегрирования по формуле трапеций при нулевых начальных условиях:

y_k+1 = y_k+(s_k+1+s_k)/2. (10.3.1)

Принимая s_k = exp(jt) и y_k = H()exp(jt), подставляем сигналы в (10.3.1) при t_k = kt, t = 1 и решаем относительно H(). Получаем:

H() = (exp(j)+1)/[2(exp(j)-1)].

H() = cos(/2)/[2j sin(/2)].

Истинное значение интеграла равно (1/j)exp(jt). Отсюда:

K() = H()exp(jt)/[(1/j)exp(jt)].

K() = cos(/2)[(/2)/sin(/2)]. (10.3.2)

Интегрирование по формуле прямоугольников (интерполяционное среднеточечное). Оператор:

y_k+1 = y_k+s_k+1/2. (10.3.3)

После аналогичных подстановок сигнала и преобразований получаем:

K() = (/2)/sin(/2).

При численном интегрировании по формуле Симпсона уравнение фильтра имеет вид:

y_k+1 = y_k-1+(s_k+1+4s_k+s_k-1)/3. (10.3.4)

Частотный анализ фильтра проведите самостоятельно. Контроль:

Рис. 10.3.1. Коэффициенты соответствия.

Графики функций К(ω) приведены на рисунке 10.3.1. При интегрировании происходит накопление результатов по всему предыдущему циклу суммирования и в этих условиях значение коэффициента K() является более представительным и информационным, чем передаточная функция оператора для одной текущей точки.

Наиболее простые формулы цифро- вого интегрирования, трапеций и прямоугольников, ведут себя различным образом в главном частотном диапазоне. Формула прямоугольников завышает результаты на высоких частотах, а формула трапеций - занижает. Эти особенности легко объяснимы. Для одиночной гармоники площадь трапеции по двум последовательным отсчетам всегда меньше, чем площадь с выпуклой дугой гармоники между этими отсчетами, и разница тем больше, чем больше частота. В пределе, для гармоники с частотой Найквиста, отсчеты соответствуют знакочередующемуся ряду (типа 1, -1, 1, -1, ... или любые другие значения в зависимости от амплитуды и начального фазового угла) и при нулевых начальных условиях суммирование двух последовательных отсчетов в формуле (10.3.1) будет давать 0 и накопления результатов не происходит. Интегрирование по площади прямоугольников с отчетом высоты по центральной точке между двумя отсчетами всегда ведет к завышению площади прямоугольника относительно площади, ограниченной выпуклой дугой гармоники.

Формула Симпсона отличается от формул трапеций и прямоугольников более высокой степенью касания единичного значения, что обеспечивает более высокую точность интегрирования в первой половине главного диапазона. Однако на высоких частотах погрешность начинает резко нарастать вплоть до выхода на бесконечность на конце диапазона (полюс в знаменателе передаточной функции рекурсивного фильтра на частоте Найквиста).

Эти особенности интегрирования следует учитывать при обработке данных сложного спектрального состава.