Ненавязчивая автоматизация ректификационной установки

# Фрагмент кода регулятора текущей скорости отбора (curQ) в зависимости от 
# текущего значения уставки (Tset) и скорости изменения температуры в колонне
# (maTcDeriv). Базовые критерии выглядят так:
# 1. Если отбор включен (состояние "старт") и текущее значение температуры в 
#    колонне (cTc) превышает нижний порог температуры (Tset - qT, 
#    где qT - квант датчика) и производная температуры в колонне положительна
#    (maTcDeriv > 0), то выполняется критерий для выключения отбора ("стоп")
# 2. Если отбор выключен (состояние "стоп") и текущее значение температуры в
#    колонне ниже верхнего порога (Tset + qT) и производная температуры в 
#    колонне отрицательна (maTcDeriv < 0), то выполняется критерий для
#    включения отбора ("старт")
# Поскольку используются цифровые датчики, то в регулятор добавлена компенсация
# цифрового "дребезга" температуры на уровне кванта датчика. Для этого вводятся
# счетчики событий 1 и 2, перечисленных выше, ctrp и ctrm соответственно. Если
# критерии 1 или 2 выполняются 3 и более раз подряд, то отбор реально 
# выключается (для 1) или включается (для 2). В противном случае выполнение 
# критериев для "стоп" или "старт" считаются "цифровым дребезгом" датчиков и
# счетчики ctrp и ctrm сбрасываются.
#
# Во время паузы ("стоп") производится "взвешенная" коррекция базовой скорости
# отбора (bQ) в зависимости от текущего ее значения и температуры в кубе. 
# Оптимальное значение веса (wt = 0.8) определяется экспериментально и задается
# при инициализации системы.
...
    if not stopFlag: # В данный момент отбор производится
      if cTc >= Tset - qT and  maTcDeriv > 0 : ctrp += 1 # Выполнен критерий 1
      elif ctrp > 0: ctrp -= 1 # Увы, это - "дребезг" температуры - сброс счетчика
      if ctrp > 2: # Порог превышен более 2 раз подряд, значит реальный "стоп"
        curQ = Qmin # Выключаем отбор
        stopFlag = True # Флаг состояния "стоп"
        ctrp = 0 # Сбросим счетчик выполнения критерия
      else: # В противном случае сохраняем текущую скорость отбора
        curQ = bQ
    else: # Сейчас - пауза в отборе (т.е. "стоп")
      if cTc < Tset + qT and maTcDeriv <= 0 : ctrm += 1 # Выполнен критерий 2
      elif ctrm > 0: ctrm -= 1 # Увы, это dct-nfrb "дребезг" температуры
      if modes[mode].reg == 1: # В этом режиме работы базовая скорость отбора
                               # корректируется с учетом веса wt старого
                               # значения (bQ) и рассчитанного по температуре
                               # в кубе (cQ)
        cQ = cR*(100.0 - Ts[hid]) # Скорость отбора, рассчитанная по т-ре в кубе
        bQ = bQ*wt + cQ*(1.0 - wt)
      if ctrm > 2: # Критерий 2 выполнен более 2 раз подряд - значит "старт"
        curQ = bQ   # Старт
        stopFlag = False
        ctrm = 0 # Сброс счетчика
    sd.Q = curQ # Установка скорости отбора в контроллере клапана отбора
...