Особенности структурно-параметрической адаптации в цифровых системах мониторинга и управления

УДК 519.711.3

ОСОБЕННОСТИ СТРУКТУРНО-ПАРАМЕТРИЧЕСКОЙ АДАПТАЦИИ В ЦИФРОВЫХ СИСТЕМАХ МОНИТОРИНГА И УПРАВЛЕНИЯ

В. Я. Карташов, С. С. Самойленко

PECULIARITIES OF STRUCTURAL AND PARAMETRIC ADAPTATION IN DIGITAL SYSTEMS OF MONITORING AND MANAGEMENT V. Ya. Kartashov, S. S. Samoilenko

В статье рассматриваются вопросы подстройки адаптивных систем управления в случаях изменения свойств объекта с течением времени. Приводится принципиальная блок-схема для таких систем управления. Выбирается метод порожденный алгоритмом непрерывных дробей на основе многочастотности, который позволяет системе адаптироваться к изменениям в режиме реального времени. Полученные результаты имеют важное практическое значение при разработке и реализации цифровых систем управления.

Questions of tuning adaptive control systems in case of changes in the properties of an object over time are discussed in the paper. A schematic block diagram for such control systems is provided.The method of continued fractions generated by the algorithm based on multifrequency is chosen, it allows the system to adapt to changes in real time. The received results have important practical value for developing and realization of digital control systems.

Самой сложной из типовых задач технической кибернетики согласно классификации данной в [1] является разработка и эксплуатация адаптивных систем управления. На начальном этапе создания таких систем, как правило, предполагается неизменность свойств объекта управления и отсутствие значительных изменений внешних условий, в которых он функционирует. Однако фактические свойства объекта и внешние условия его функционирования с течением времени изменяются, причем эти изменения могут привести к непредвиденным последствиям. В таких ситуациях необходимо разрабатывать такие системы мониторинга и управления, которые могли бы и внешние, и внутренние изменения идентифицировать и преждевременно вносить определенную коррекцию в систему управления. Для решения подобных задач и создаются адаптивные системы управления, реализованные на уникальных возможностях и достижениях цифровой техники.

Способность системы управления определять изменение и изменять структуру системы управления и составляет содержание понятия ее адаптивности. Исследованию и построению адаптивных систем управления объектами различной сложности посвящено значительное число публикаций [1; 2; 3; 5; 11; 13; 14; 15], в которых достаточно подробно рассматриваются классификация, особенности и возможный математико-алгоритмический аппарат. В то же самое время остается ряд основных нерешенных и неисследованных проблем в этой области знаний. Сначала заметим, что в данной работе не рассматривается случай изменений внешних условий - рассматриваются изменения свойств самого объекта. В этом случае будем рассматривать аналитические активные замкнутые динамические системы [3]. Как известно [15], изменения внешних возмущений могут быть исключены системой разомкнутого управления. Для устранения влияния непредсказуемых возмущений (особенно внутреннего характера для объекта) требуется управление с обратной связью. Для рассмотрения более подробно особенностей цифровых систем управления с элементами адаптации приведем высказывания исследователей об изменчивости свойств объектов управления.

Так в [1] авторы отмечают, что «в подавляющем большинстве технологических процессов непрерывная адаптация объекта в течении всего времени работы не нужна, поскольку непредсказуемые (непредвиденные) изменения объекта происходят относительно медленно». Однако эта оценка дается на качественном уровне в силу того, что для запуска алгоритма адаптации необходимо знать об изменении свойств или характеристик объекта управления. Будем считать, что управляющее устройство, реализованное на элементах цифровой техники, фактически не изменяется.

С другой стороны часто встречаются высказывания о довольно частом изменении характеристик объекта [1], причем частота этих изменений связана со сложностью объекта: чем сложнее объект, тем чаще происходят изменения его свойств.

Приведем наиболее часто встречающуюся [4; 6; 12; 15] блок-схему самонастраивающегося регулятора или схему адаптацию по замкнутому циклу (рис. 1).

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

практически случайно при значительных интервалах альном их функционировании, либо текущем его анавремени (при периодическом осмотре объекта) в силу лизе в режиме on-line. Хочется сразу заметить, что

дорогостоящих, сложных и т. д. процедур при специ- такие способы слабо теоретически разработаны.

Рис. 1. Блок-схема адаптивной системы с настройкой регуляторов по замкнутому циклу: Я - саморегулирующийся регулятор, ЯЯ - расчет регулятора; ИО - идентификация объекта; О - объект; п> - задающее воздействие; и - управляющее воздействие; е- ошибка регулирования; у - управляемая переменная

Таким образом, блок-схема адаптивной системы (рис. 1) имеет два контура, в которых кроме общей прямой части «Я^-О» имеется два обратных контура: один по управлению от у к элементу сравнения и второй по определению и идентификации модели объекта и настройки регулятора ЯЯ: (и,у)^-ИО^ЯЯ^-Я. Так как рассматривается цифровой регулятор, то на выходе его формируется напряжение определенного

уровня, которое сохраняется до следующего отсчета. Обозначим этот интервал - интервалом управления , тогда в зависимости от изменения ошибки управления е(^) на выходе цифрового регулятора формируется напряжение (рис. 2):

иОО = и(^ 8(д), при \\Ди < 1 < 0 + 1)Д^

Рис. 2. Графическое представление функционирования цифровой замкнутой системы

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

{у(0), у(Д:и),у(2Д:и),...,у(пД:и),...} .

С учетом суждений о частоте изменений свойств объекта необходимо изменить блок-схему адаптивной системы (рис. 1), поэтому в некоторых исследованиях вводится соответствующий элемент, учитывающий обсуждаемую неопределенность временного обнаружения изменения.

В настоящее время достаточно интенсивно развивается теория и реализация систем адаптивного управления на основе цифровой техники, в которых большое значение уделяется текущей идентификации

в реальном времени [6]. Согласно блок-схеме (рис. 1) синтез цифрового регулятора в большинстве случаев осуществляется с использованием упомянутых выше параметрических моделей объектов и возмущений. Для параметрической идентификации в режиме реального времени разработаны и для ряда типов объектов управления хорошо зарекомендовали себя рекуррентные методы оценивания параметров моделей объектов и сигналов [4; 6; 15]. Приведем наиболее распространенную общую блок схему такой адаптивной системы (рис. 3).

Если проанализировать блок схему (рис. 3), то можно сформулировать следующие замечания:

- такой принцип адекватности используется при параметрической настройке параметров модели объекта управления, структура которой заранее известна

и остается постоянной на все время функционирования системы;

- в каждый момент времени 1анализируется ошибка е^) = у^) - у(0; исследуется последовательность значений параметров с помощью различных методов: МНК и его модификаций, стохастической аппроксимации, фильтрация на основе ортогональных систем функций, рекуррентный метод максимального правдоподобия, рекуррентный метод вспомогательных переменных и т. д.

В блоке настройки модели определяются оптимальные значения параметров , которые затем

преобразуются в соответствующий набор параметров управляющего алгоритма. Значительное число исследований [4; 11] посвящено параметрическим адаптивным системам.

Рис. 3. Блок схема адаптивной системы с параметрической идентификацией объекта управления (ОУ): ММО - математическая модель объекта, БНМ - блок настройки модели; ВУ - вычислительное устройство, ИО - исполнительный орган. Я - регулятор, {С}, {С^ } - параметры ММО и их оценки,

{ ^ } - параметры регулятора

Однако актуальной и наиболее перспективной с практической точки зрения является задача построения адаптивных систем управления, в которых во время функционирования объекта появляются структурно-параметрические изменения, которые обусловливают развитие способов решения задач структурнопараметрической идентификации. Одним из таких методов является структурно-параметрическая идентификация (SP-идентификация) на основе дробнорациональной аппроксимации алгоритмами непрерывных дробей [7; 8]. Основа данного подхода определяется фундаментальным понятием теории автоматического управления [5] - передаточной функции динамического объекта, допускающего линеаризацию его математической модели. При моделировании используется согласованное z-преобразование, позволяющее с помощью правильных C-дробей из теории непрерывных дробей получать дискретные передаточные функции по измерениям вход-выходных переменных с периодом дискретизации At = const.

Кратко способ 8Р-идентификации представляет последовательность этапов, приводящих к таким основополагающим принципам осуществимости свойства адаптации математических моделей динамических систем. Прежде всего предполагается, что объект исследования имеет конечную память в динамике, т. е. если этот объект допускает линеаризацию для описания его поведения, то в качестве математической модели с достаточной точностью аппроксимации может быть использовано обыкновенное дифференциальное уравнение п-го порядка, где п определяет память объекта. Изменение п приводит к изменению структуры объекта. Коэффициенты дифференциального уравнения определяют динамические свойства объекта, такие как, устойчивость, апериодичность, колебательность выходной переменной и ряд других. Они получили название параметров объекта. 8Р-идентифи-кация предполагает определение на основе измеренных значений вход-выходных переменных в равноотстоящие моменты времени всех структурных и параметрических значений объекта. Погрешность измерений уменьшается за счет фильтрации, описанной в работе [7], поэтому полагаем, что погрешность мала и не влияет на практическую реализуемость результатов работы.

Первый этап разработки способа SP-идентификации заключается в определении взаимно-однозначного отображения объекта в непрерывном и дискретном представлении объекта исследования при фиксированном периоде дискретизации At = const. Достаточно тщательное исследование этой задачи, включающее теоретические, многочисленные модельные и реальные изыскания, показали, что таким отображением является согласованное z-преобразование [9]:

Z = epAt (1)

Х0 Х1 X2 X3 . .. X n

Уо У1 У 2 Уз . .. Уп

а10 a11 ai2 ai3 . .. ain

a20 a21 a22 a23 . .. a2n

am 0 am1 am 2 a m3 . .. amn

в которой первые две строки определяются исходными данными, а последующие элементы определяются рекуррентно по соотношению [10]:

где 7,р- комплексные переменные в дискретном представление объекта, в непрерывной форме интегрального преобразования Лапласа, соответственно, А1 - период дискретизации.

Сложность этого отображения обусловлена его многозначностью, которая привела к выделению основной полосы на комплексной плоскости р и к необходимости введения многолистной римановой поверхности. Свойства согласованного 7-преобразования (1) определяют значимость вариации периода дискретизации, в частности для каждого динамического объекта существует интервал значений Дt е (Д^т; Д^ах) сохранения взаимно-однозначного соответствия между основной полосой р-плос-кости и нулевым листом 7-плоскости. Более подробно это свойство приведено в [7]. Второй задачей, которая относится к базовой, является задача получения дискретной передаточной функции по совокупности измерений вход-выходных переменных в равноотстоящие моменты времени. По определению дискретная передаточная функция является отношением 7-преобразования выходной переменной к 7-преобразованию входной переменной при нулевых начальных значениях у(1):

am-2,n+1 am-1,n+1

am-2,0 am-1,0

где а-х п = хп; а0 п = уп. Расчет идентифицирующей

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

G(z, At) =

G (Г, At) =

= У (0) + y(At) z-1 +... + y(nAt) zn +... x(0) + x(At) z-1 +... + x(nAt) z ~n +...

= У( z, At) x( z, At)

Так как в работе рассматриваются динамические системы с конечной памятью, то для нахождения Сд (и, ДЦ используется итерационный алгоритм теории непрерывных дробей [9] - метод В. Висковатова, который определяет идентифицирующую матрицу и определяет дробно-рациональную аппроксимацию соотношения (2):

am0Z~ 1 + ...

Свойства согласованного 7-преобразования (1) и алгоритм дробно-рациональной аппроксимации (3) порождают два важных принципа моделирования дискретных систем:

^ е (Д^п; Д^ах) ;

Таким образом, 8Р-идентификация динамических объектов и выше сформулированные принципы позволяют привести обобщенную функциональную блок-схему системы адаптивного управления и мониторинга с возможностью определения структурных изменений динамического объекта (рис. 4).

Для реализации адаптивной цифровой системы управления при возникновении структурных изменений динамического объекта одной из основопола1

гающих задач является определение этих изменений. Для этого используем многочастотное квантование в следующей схеме: для управления выбираем большой интервал дискретизации Дt е (Д^,п; Д^ах), при котором существует дискретная передаточная функция С0(и, Д^) [8], для идентификации структуры и параметров модели выбираем период дискретизации Дt << Дtu, Дt е (Дtmin; Дtmax). Тогда на интервале ^ + Дtu) осуществляем опрос в М моменты времени [^;^ + Дt;...; tk + тД^...;^ + MДt] вход-выходных переменных на указанном временном множестве {х(^ + тД^,у(^ + тД^>, где т = 0;М. Это

множество измерений, является исходными данными для проведения 8Р-идентификации, решая которую определяем, что изменений не произошло, появились только параметрические или структурнопараметрические изменения, по сравнению с математической моделью объекта управления на предыдущем интервале управления ^к-1; ^-1 + Д^). Решение этой задачи является ключевым моментом в процессе адаптации. Для исследования этого вопроса необходимо иметь моделирующую систему, включающую имитацию входного ступенчатого сигнала и реакции объекта на этот сигнал (рис. 2).

Рис. 4. Обобщенная блок схема адаптивной цифровой СУ при структурных изменениях модели ОУ

Назовем моделирующую систему генератором вход-выходных сигналов при условии, что известна передаточная функция объекта управления ОУ (рис. 4), причем входной сигнал представляет ступенчатый сигнал, характеризующий функционирование цифрового регулятора Я, а выходной сигнал ОУ порождается входным сигналом и(1) и имитируемой передаточной функцией СОУ (р). Для анализа такого фрагмента системы управления (рис. 4) воспользуемся следующим представлением этой части системы (рис. 5).

Таким образом, входной сигнал имитируется сигналом:

иоо = £ дда - а -1)ду - 1а - iду]. (4)

Соотношение (4) можно представить в следующей форме, которая изображена на рис 5 пунктирной линией (при условии, что Д0 = 0 )

иОО = (Д; - Дi-l)l(t - с -1)ду =ддоо +

п-1 ;=1 (5)

+£ (Д+1 - дi)1(t - ;ди - дn1(t - nДtu).

Таким образом, нетрудно получить входной и выходной сигналы для любого линеаризуемого объекта конечного порядка. Предположим, что на входе объекта управления с передаточной функцией в(р) генерируется сигнал х(1) имеющий вид (5) (таблица).

Рис. 5. Принцип функционирования имитационной системы

Таблица

Время 1 < t [0, лу [ЛЪи,2ЛЪи) [(п _ 1)луплу t > пЛ^

х(1) 0 Д2 Д 0

Тогда преобразование Лапласа представляется в виде соотношения:

х(р) = 1 [Д1 + £ (д ;+1 - А; )е- ; рД1и - Дпв-прД1и ]. (6) р ;=1

Заметим, что временной отрезок [0, п Д^] отражает функционирование цифрового регулятора и его особенности выбираются в зависимости от целей имитации процесса управления путем изменения

Д,-.1^]

Ді1-1[Срр)] + (Д2 _ Д1)1-1[Срр)(

-1г6(Р)П

_і[6(Р),

Д^1^] + (Д2 _ Ді)!- [:р

+ (Д; _ Д: 1)-_1[^(Р)Є_(і_1)РЛСи] Р

Ді-_1 [Срр)] + ^ (Д: _ Д _1 )^1 [ е

р і=1 р

Дtu, п, {Д ! >п и самое главное - изменения передаточной функции в(р) практически на каждом или некоторой совокупности интервалов управления Д^ .

Учитывая свойства преобразования Лапласа, нетрудно получить сигнал, описывающий реакцию объекта на сигнал (таблица):

, если t < 0 , если 0 < t < Д^

, если Дtu < t < 2Д^

.. +

, если (i-1)Дtu < t < ;Д^

] , если (п-1)Дtu < t < п Д^

п+1

Д^1^] + £ (Ді _ Д _і)-_1[.

іг6(р)

-Дп-_1[^б(р),

, если t > пЛtu

На рис. 5 приведены схемы структурных особенностей вход-выходных сигналов, полученных путем имитации функционирования автоматической цифровой системы во внутреннем контуре.

При такой генерации разнообразных сигналов возникает задача идентификации передаточной функции на каждом интервале управления Д^ . То, что это

очень сложная проблема довольно часто обсуждается в публикациях. Так, например, В. Я. Ротач [1] считает такое представление ошибочным. Трудно не согласиться с В. Я. Ротачем, но всегда ли это так, скорее всего это зависит от динамических свойств объекта и от возможностей цифровой техники.

В конце данной работы приведем пример работы такого генератора, моделирующего некоторые изменения в объекте.

Пример.

Возьмем апериодический объект первого порядка, имеющий характеристики:

Тб +1

где к- коэффициент усиления,

Т - постоянная времени,характеризующая инерционность объекта.

Входной сигнал задается выражением:

хОО = £ Д;[1^ - 0 - 1)ДЦ - 1^ - ^)],

где п - количество ступенчатых воздействий,

Дt - интервал управления.

Для определенности возьмем количество ступенчатых воздействий равное 5, коэффициент усиления равный 1, постоянную времени равную 1, ступенчатые воздействия равные: Д1 = 3, Д2 = 1, Д3 = 2,

Д4 = 0, Д5 = 0 и интервал управления равный

Дt = 5.

Рис. 6. Сигнал на входе

тНтНтНтНтНтНтНтНтНтНС\\1С\\1С\\1С\\1С\\1С\\1

время

Рис. 7. Сигнал на выходе

Выводы

Данную работу следует рассматривать как обзорную, в которой обсуждены принципы адаптации в автоматических системах. Наиболее разработанной и предлагаемой на рынке цифровых систем является параметрическая адаптация, в то же время в реальных условиях объект управления по каким-либо причинам претерпевает структурно-параметрические изменения. В статье анализируются особенности таких заЛитература

- Кемерово, 2011. - 214 с.

Информация об авторах:

Карташов Владимир Яковлевич - доктор технических наук, профессор, заведующий кафедрой автоматизации и технической кибернетики КемГУ, 8-905-949-2136, v.kartashov.aitk@gmail.com.

Vladimir Ya. Kartashov - Doctor of Technical Sciences, Professor, Head of the Department of Research Automation and Technical Cybernetics, Kemerovo State University.

Самойленко Сергей Сергеевич - аспирант кафедры автоматизации и технической кибернетики КемГУ, 8-923-522-6442, my niserg@mail.ru.

Sergey S. Samoilenko - post-graduate student at the Department of Research Automation and Technical Cybernetics, Kemerovo State University.

Статья поступила в редколлегию 09.04.2014 г.

дач, предлагается подход реализации адаптационных процессов, включающий способ SP-идентификацию динамического объекта, принцип эквивалентности дискретных моделей непрерывного объекта и многочастотное квантование по времени. Для исследования потенциальных возможностей такого адаптационного механизма в статье предложен генератор вход-выходных сигналов объекта управления, работа которого проиллюстрирована на примере.