Одной из важнейших задач функционирования системы автоматического контроля и управления технологическими процессами является оперативное обнаружение отклонений параметров технологического процесса от заданного поля допусков. Информация от датчиков поступает на управляющую вычислительную машину, где обрабатывается, после чего выдаются управляющие сигналы на исполнительные механизмы в случае отклонения от нормального хода технологического процесса или аварийных ситуаций. Одна из главных задач состоит в своевременном обнаружении отклонений технологического процесса от заданного поля допусков.
Следует отметить, что параметры контролируемых технологических процессов носят случайный характер, а при возникновении тревожных и аварийных ситуаций, носят ярко выраженный нестационарный характер, как по математическому ожиданию, так и по дисперсии.
Существующий интегральный метод обнаружения нестационарности в случайном контролируемом процессе не позволяет своевременно обнаруживать отклонения от заданного поля допусков, так как требует накопления показаний датчиков за определенный промежуток времени, что приводит к задержке обнаружения отклонения технологических параметров. Сущность интегрального метода заключается в том, что в качестве критерия обнаружения нестационарности в случайном контролируемом процессе X(t) принимается площадь S(t), ограниченная X(t) и сглаженной реализацией 
. Если величина площади отклонения S(t) больше некоторого критического значения ξ, то принимается решение о наличии отклонения (нестационарности) от заданного поля допусков.
В работе предлагается использование оценки дисперсии случайного процесса в качестве критерия обнаружения нестационарности, который значительно быстрее реагирует на изменения в контролируемом процессе, чем интегральный метод.
В реальных условиях контролируемые процессы являются случайными, содержащими полезный сигнал и помехи. для решения задачи оперативного обнаружения нестационарности необходимо вычислять сглаженные значения и дисперсии случайных процессов изменения контролируемых параметров с целью выделения полезного сигнала от помех.
Один из известных способов, используемых для получения оценок среднего значения и дисперсии на основе n независимых измеренных значений контролируемого случайного процесса, заключается в вычислении оценок по формулам:
(1)
Недостатком вычислений параметров по формулам (1) является необходимость накапливания значений данных за определенный промежуток времени, что приводит к задержке обнаружения существенных нестационарностей в контролируемом случайном процессе.
Оценку сглаженного значения можно вычислять с помощью известного оператора экспоненциального сглаживания [1]
, (2)
где α = i1/n.
Выражение по формуле (2) позволяет вычислять сглаженное значение 
в режиме реального времени без накопления промежуточных данных, т.е. без задержки. для этого требуется только знание текущего значения контролируемого параметра и сглаженное значение вычисленное на предыдущем такте опроса датчиков.
Для определения оценки дисперсии без накопления сумм можно использовать выражение полученное в [2]:
. (3)
При моделировании вычислялись оценки сглаженного значения и дисперсии по формулам (2), (3). В случае, когда замеры Xi снимаются через равные промежутки времени (∆t = const), имеем 
= const.
На рисунке 1 показаны результаты моделирования процесса обнаружения нестационарности в случайном контролируемом процессе. Сглаживание необходимо для отделения полезного сигнала от помех.
Вычисление дисперсии Di (см. рис. 1) позволяет определить момент наступления существенной нестационарности в условиях аддитивных помех гораздо быстрее, чем интегральным методом. Решение о начале нестационарности принимается при условии Di > ξ. Величина ξ уточняется в процессе эксплуатации системы для каждого контролируемого объекта в реальных условиях.
Например, после очередного вычисления оценки дисперсии по формуле (3) при n=57 (n – количество независимых измеренных значений контролируемого случайного процесса) получено значение Di > ξ (см. рис. 1). В результате принимается решение, что обнаружено отклонение технологического процесса от заданного поля допусков. для повышения надежности обнаружения отклонения следует формировать управляющие сигналы после трех и более подряд идущих превышений дисперсии допустимой нормы.
Кроме того, для повышения быстродействия обнаружения отклонений от заданного поля допусков следует использовать адаптивную дискретизацию контролируемых процессов (∆t ≠ const), при которой частота опроса контролируемых датчиков меняется в зависимости от изменения их показаний.
При изменении интервала опроса датчиков меняется коэффициент сглаживания 
. В этой связи необходимо определить зависимость коэффициента 
для выражений (2), (3) от параметров ∆t, T. Оператор (2) является дискретным аналогом апериодического звена первого порядка с постоянной времени T и коэффициентом усиления равным единице. Решение дифференциального уравнения, описывающего такое звено, имеет вид:
, (4)
где x(τ), 
– процессы, соответственно на входе и выходе звена.
При подаче на вход звена единичной ступенчатой функции x(τ)=1 из выражения (4) получим
= 1 – e-t/T
Переходя к дискретному представлению при tn = n*∆t и раскрывая рекурентную формулу (2) получим
(5)
Важность полученной формулы (5) заключается в том, что она характеризует зависимость коэффициента сглаживания α от заданной постоянной времени сглаживания T, при изменяющемся интервале времени между измерениями показаний датчика ∆t.
Следует отметить, что использование формулы (5) сопряжено с определенными трудностями, так как она требует подключения подпрограммы вычисления экспоненты. для устранения этой проблемы разложим функцию 
в степенной ряд вида:
Моделирование процесса обнаружения нестационарности
После разложения функции в степенной ряд и несложных преобразований получим более простую формулу:
, (6)
где L – нормирующий коэффициент (0 ≤ L ≤ 1).
В процессе моделирования на ЭВМ было получено значение L= 0,582, при котором погрешность вычисления 
по формуле (6) по сравнению с (5) не превышает 1 %. Таким образом, мы получили зависимость
. (7)
Полученная зависимость (7) позволяет использовать оператор экспоненциального сглаживания для оперативного вычисления сглаженных значений и дисперсий контролируемых процессов при изменяющемся интервале опроса датчиков. Один из вариантов алгоритма адаптивного опроса датчиков предложен в работе [2].
Оперативное обнаружение отклонений параметров контролируемого процесса от заданного поля допусков дает возможность для своевременной выработки управляющих воздействий на исполнительные механизмы, что повышает эффективность и надежность функционирования автоматической системы контроля и управления технологическими процессами.