Управление объектом с запаздывающим управлением по выходу в условиях параметрической неопределенности

МАТЕМАТИКА

УДК 62-50

А. М. Цыкунов Астраханский государственный технический университет

УПРАВЛЕНИЕ ОБЪЕКТОМ С ЗАПАЗДЫВАЮЩИМ УПРАВЛЕНИЕМ ПО ВЫХОДУ В УСЛОВИЯХ ПАРАМЕТРИЧЕСКОЙ НЕОПРЕДЕЛЕННОСТИ

Введение

Задача управления объектами с запаздыванием в управлении в условиях априорной неопределённое™ их параметров занимает особое место в теории управления. Это обусловлено тем, что устойчивость и качество переходных процессов в замкнутой системе зависят от величины запаздывания. Вследствие этого для улучшения качественных показателей замкнутой системы часто используют различные прогнозирующие устройства [1, 2]. Однако это не снимает всех проблем, а именно таких, как требования к устойчивости разомкнутой системы. Хорошо известно [3], что для неустойчивого объекта управления можно построить систему управления только в том случае, когда запаздывание по сравнению с доминирующей постоянной времени объекта управления является малой величиной. Эта проблема усугубляется в случае априорной неопределённости параметров.

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

В данной работе предлагается одна из возможных схем построения адаптивной системы управления без прогнозирующего устройства в сочетании с хорошо известной схемой расширенной ошибки.

Постановка задачи

Рассмотрим объект управления, динамические процессы в котором описываются дифференциальным уравнением с запаздывающим аргументом

где Q(р), Я(р) — нормированные линейные дифференциальные операторы с постоянными коэффициентами; у(V) — скалярная регулируемая переменная; и(^) — скалярное управляющее воздействие; р = & / & — оператор дифференцирования; к > 0 ; к — постоянное время запаздывания. Эталонная модель задана уравнением

где ут (V) — выход модели; г(V) — задающее воздействие; кт > 0, Qm (р), Ят (р) — нормированные линейные дифференциальные операторы с постоянными коэффициентами.

Требуется построить систему управления, обеспечивающую выполнение целевого условия

Q(р)у(0 = кЯ(р)ы(* — к), и (у) = 0,5 е [—к, 0),

Qm (р)Ут ({) = ктКт (р)г(* — к\

Нш(у(^) — Ут (V)) = 0

при следующих ограничениях:

Предположения: 1. Полиномы О(1), От(1), Л(1), Ят(1) - гурвицевы, где 1- комплексная переменная в преобразовании Лапласа. 2. Коэффициенты полиномов О(1), ^(1) и величина к зависят от вектора неизвестных параметров X. Известны порядки полиномов О(1), Qm (1), Л(1), Ят (1) - п, п, т, т соответственно, п -1 > т .

Метод решения

Преобразуем уравнение (1) к виду

у

+ 8

Здесь 5^) — экспоненциально убывающая функция, зависящая от начальных условий. Выделим в выражении Qm (р) / Q(р) целую часть, т. к. deg Qm = п, deg Q = п . Тогда получим

у

От

У

где АQ(р) = Q(р) — Qm (р), deg АQ = п — 1.

Если все параметры объекта управления известны, то закон управления

н^) = ^

Q(P) ^т {рШр) к

обеспечивает следующий вид уравнения замкнутой системы:

У

Qm

т. е. переходные процессы в замкнутой системе после затухания составляющей 5(^) будут такие же, как в эталонной модели.

Так как параметры объекта управления неизвестны, то возьмем настраиваемый фильтр

состояния (&и (О = ^тСн (О + Ьи (V) + *рг (ОС и (V), где £ и е Яп ; ЬТ = [0, ..., 0, 1]; Ат — гурвицева матрица в форме Фробениуса с характеристическим многочленом Qm (1); Р(^) — настраиваемый вектор.

Так как матрица Ат задана в форме Фробениуса, то существует вектор Р0 такой, что матрица Ат + будет иметь характеристический полином Q(1). Поэтому уравнение фильтра преобразуется к виду

£ и

Ц(() = дТ £и (V) + ^0С и» (V),

где £и (V) е Яп; д, д0 — коэффициенты полинома Qm (1), £ит (V) — производная п -й компоненты вектора £и (V). Если (Р(^) — Р0)Т £и (V) будет стремиться к нулю при V , то в пределе

си (() =

1 Рп-1

-и

, ^

_ Q

ления (6) может быть реализована в виде 0Т(0Си(V), где 0^) — настраиваемый вектор параметров. Для формирования второй составляющей возьмём фильтр

£

где £(0 е Ят ; ^ — матрица в форме Фробениуса, с характеристическим полиномом Ят(1). Тогда имеем

Г (t) =

u

Rm

что позволяет реализовать вторую составляющую в (6) в виде (t)v(t), где m(t) - настраивае-

мый вектор. Возьмём следующий закон управления:

u

где a(t) — настраиваемый параметр. Тогда из уравнения (5) получим

y(t) = kRm-(P) ((0(t—h)—0o)T Z u (t—h)+(m(t—h)—mo )T V(t—h)+

Qm

+

где вектор 0o составлен из коэффициентов полинома DQ(l), а вектор m0 - из коэффициентов

полинома DR(l) с противоположными знаками, записанными в обратном порядке; a0 = km /k .

В идеальном варианте вектор b(t) должен стремиться к b0, а вектор 0(t) - к 00, но b0 = —00, поэтому будем настраивать вектор 0(t), а вектор b(t) в (8) будем брать b(t) = —0(t). Тогда уравнение фильтра (8) примет вид

z u

h(t) = qT Zu (t) + q0Z um (t).

В результате все фильтры оказались независящими от неизвестных параметров объекта управления, а следовательно, они реализуемы.

Введем векторы CT (t) = [0T (t), mT (t), a(t)] , c0 = [0T , m0 , a0] , WT (t) = [ZT (t), XT (t), r(t)] и из уравнений (3) и (11) получим уравнение ошибки e(t) = y(t) — ym (t) без учета экспоненциально убывающей функции d(t) .

e

Qm

Получили обобщенный настраиваемый объект, для которого применима любая известная схема расширения [6], например следующая:

Є1

Qm (P) ,

g

где t(t) - настраиваемый параметр, c > 0. Уравнение расширенной ошибки

eP (t) — e(t) + ej (t) будет иметь вид

eP (t) = W (P)(k (c(t — h) — C0 )T J(t — h) + (t(t) — k) g (t — h) —

— Cep

T

Выберем полиномы T(p), Rm (p), Qm (p) так, чтобы передаточная функция m

^

имела вид —1—, где а > 0 . Тогда уравнение для расширенной ошибки запишется следующ

1 + а

образом:

ep (t) = ~aep (t) + (k(c(t - h) - c0 )T J(t - h) + (r(t) - k)g(t - h) -

-Xep

Утверждение. Пусть выполнены условия предположений и алгоритмы настройки имеют вид

ep

i+e2

¿

где р > 0, Г = ГТ > 0 . Тогда выполнены целевые условия (3) и все сигналы в замкнутой системе ограничены.

Доказательство утверждения. Проинтегрируем первое уравнение в

[t - h; t] c(t) _ c(t - h) - I —----^------ds . Тогда уравнение (15) можно записать в виде

Л 1 + e:

T г ep (s)PJ(s - h) T

e,

Л 1 + e2

t-h

где c (t) _ c(t) - c0, t(t) _ t(t) - k .

Возьмём функционал Ляпунова - Красовского

V(t) _ ln(1 + ep(t)) + kcT(t)r-1c(t) +112(t) + p P

0 t

+ Jdv J s

-h t+v

и вычислим полную производную на траекториях системы (16), (17):

V(t) _ -2ac(t)e2 (t) + 2c(t)ep (t) JT (t - h) J s(s)ep (s)J(s - h)ds

t-h

t

- hs(t)ep (t) JT (t - h)rj(t - h) - J a(s)e2 (s)JT (s - h)J(s -

t-h

Ce p

p

Выделив полный квадрат и учитывая неравенство o(t) £ 1, получим

t

V(t) - 2c(t)e2 (t)(a + JT (t - h)(ßI - rh)J(t - h)) - J gT (s, t)Г-1 g(s, t)ds,

t-h

где g (s, t) _ (ep (t)JT (t - h)r- e (s)JT (s - h)r). Тогда, выбрав c и Г из условия %I - hr> 0,

JT (t — h)G-ep (- *T‘

получим V(t) < — ac(t)ep , откуда следует lim ep (t) = 0 . Далее доказательство того, что выполнено

целевое условие (3), полностью аналогично доказательству в [6] и поэтому здесь не приводится.

Пример. Уравнения объекта управления и эталонной модели имеют вид: (p3 + 2p2 + 2p + l)y(t) = 4u(t-5), (p3 + 3p2 + 3p +1)ym(t) = 2r(t-5). В фильтре (12) Zu G R3, Zu(t) = A,Zu (t) + (t)h(t) + ¿a(t)r(t), h(t) = qTZu (t) + qoZuffl (t), матрица 4m имеет характери-

стический полином l3 + 3l2 + 3l +1, wT (t) = [Z U (t), h(t), r(t)]. Схема расширения имеет вид

e (t) = Tp\ (g (t - h) x(t) - Cep (t )JT (t - h) J(t - h)), g (t) = (cT (t )T-1 (p) - T-1 (p)cT (t ))Jt),

Qm

T(p)J (t) = w; (t) , где T(p) = p2 + 2p +1. На рисунке приведены графики переходных процессов в системе при нулевых начальных условиях, когда параметры системы управления имели значения: Г = diag{0,3 0,3 0,3 0,3}, р = 0,5, % = 20 .

0,6

0,4

0,2

-0,2

11(f) n

/Т

J

; t, С

100 200 300

Переходные процессы в системе СПИСОК ЛИТЕРА ТУРЫ

1. Паршева Е. А., Цыкунов А. М. Адаптивное управление объектом с запаздывающим управлением со скалярными входом-выходом // Автоматика и телемеханика. - 2001. - № 1. - С. 142-149.

2. Цыкунов А. М. Адаптивное управление с компенсацией влияния запаздывания в управляющем воздействии // Изв. РАН. Теория и системы управления. - 2000. - № 4. - С. 78-81.

3. Колмановский В. Б., Носов В. Р. Устойчивость и периодические режимы регулируемых систем

с последействием. - М.: Наука, 1981. - 448 с.

4. Фуртат И. Б., Цыкунов А. М. Адаптивное управление объектами с запаздыванием по выходу // Изв.

вузов. Приборостроение. - 2005. - Т. 48, № 7. - С. 15-19.

5. Niculescu S. I., Annaswamy A. M. An adaptive Smith - controller for time delay systems with relative degree

n £ 2 // Systems and Control Letters. - 2003. - Vol. 49, N 5. - P. 347-358.

6. Narendra K. S., Annaswamy A. M., Singh R. P. A general approach to the stability analysis of adaptive sys-

tems // Int. J. Control. - 1985. - Vol. 41, N 1. - P. 193-215.

Статья поступила в редакцию 1.10.2006

MANAGEMENT OF AN OBJECT WITH TIME-DELAY MANAGEMENT ON OUTPUT IN CONDITIONS OF PARAMETRICAL UNCERTAINTY

A. M. Tsykunov

The problem of construction of an adaptive control system for an object with time-delay management and unknown parameters is considered in the paper. The initial model is transformed, and algorithms of management are synthesized by means of augmented error method.