Метод повышения робастных свойств линейной замкнутой системы в условиях структурно-параметрической неопределенности описания передаточной функции объекта управления Текст научной статьи по специальности «Математика»

УДК 004

DOI 10.25205/1818-7900-2018-16-1-86-99

О. Ф. Данилов \ А. Н. Паршуков \ И. Л. Полянская 2

1 Тюменский индустриальный университет ул. Володарского, 38, Тюмень, 625000, Россия

2 Уральский государственный университет путей сообщения (филиал в г. Тюмени) ул. Калинина, 5, Тюмень, 625000, Россия

danilov_of@mail.ru, anparshukov@mail.ru,polyanskaya_il@inbox.ru

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

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

Ключевые слова: робастность, модальный регулятор, структурно-параметрическая неопределенность.

Введение

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

о (п, p) У ^) = Ъ (т, p) п ), an = 1, П > m, (1)

здесь u - входной (управляющий) сигнал, у - выходной (управляемый) сигнал, t - непрерывное время; о (п, p) и Ъ (т, p) - дифференциальные операторы 1.

Здесь и далее под записью f (/, p) понимается полиномиальный оператор в степени Ь

f 0,P) = • ^ ,

г = о

1 Далее предполагается, что операторы о (п, p) и Ъ (т, p) взаимно простые, т. е. не имеют общих нулей.

Данилов О. Ф., Паршуков А. Н, Полянская И. Л. Метод повышения робастных свойств линейной замкнутой системы в условиях структурно-параметрической неопределенности описания передаточной функции объекта управления // Вестн. НГУ. Серия: Информационные технологии. 2018. Т. 16, № 1. С. 86-99.

ISSN 1818-7900. Вестник НГУ. Серия: Информационные технологии. 2018. Том 16, № 1 © О. Ф. Данилов, А. Н. Паршуков, И. Л. Полянская, 2018

гдеЦ (г е 0,1) - постоянные коэффициенты, р1 = ё1 ¡й^ - оператор дифференцирования по времени 2. Полиномиальному оператору / (I, р) соответствует алгебраический полином

/ (I, г):

f (l, s) = £ f

i

s'

s i

i=1

где s - переменная преобразования Лапласа.

Оператор f (l, p) (и соответствующий полином f (l, s)) будем называть «точечным», подразумевая, что в пространстве коэффициентов данному оператору соответствует точка.

Качество управления назначается ограниченной односвязной областью S , определяющей допустимое расположение полюсов передаточной функции (ПФ) на комплексной плоскости C1 . Предполагается, что область S расположена слева от мнимой оси. В технологии синтеза модального регулятора (описанной в работе [1]) регулятор R ищется в виде дифференциального уравнения l-го порядка:

ß (i, p) u (t) = a (l, p) y (t) + X(i, p)g (t), ßi = 1, (2)

здесь g (t) - заданный эталонный сигнал. В результате ПФ замкнутой системы имеет вид

Wcl( s) = b (ms) -x(l, s) , (3)

a(n,s)-ß(l,s)-b(m,s)-a(l,s)'

Известно, что для объекта управления, заданного дифференциальным уравнением n-го порядка, любое заданное расположение полюсов ПФ (3) можно обеспечить динамическим регулятором (2) порядка l = n — 1 (и выше) [2]. Коэффициенты регулятора находятся из уравнения

aet (2n — 1,s) = a(n,s)-ß(n — 1,s) — b(m,s) •a(n — 1,s), ae2tn—1 = 1, (4)

где за aet. {2n — 1,s) принят характеристический полином эталонной системы. Полином aet. (2n — 1, s) выбирается произвольным образом при условии, что все множество его корней

Л (aet ) = { X: aet. (2n — 1, X) = 0} с C1, лежит внутри области S , т. е. выполняется

Л^ ) с int S . (5)

Приравнивая коэффициенты при одинаковых степенях переменной s в левой и правой частях уравнения (4), получаем систему линейных алгебраических уравнений: C • R = D , где C - матрица C = | cij j, i е 1,2n — 1, j e 1, 2 n — 1, составленная из коэффициентов полиномов a (n, p) и b (m, p) объекта (1) по правилу

0, i < j,

c>j =

ai—j, i e j, n + j, 0, i > n + j,

2 Подразумевается, что p0 = 1

при i е 1,2 n -1, j е 1, n -1, и

0, i < j - n +1,

с.. = s

ij

- b +. . ,, i е j - n +1, да + j - n +1,

n+l-J-1 ^ J ^ J '

0, i > m + j -n +1,

при i е 1, n1 +1, j е l +1,21 +1; R = col (P0, •••, Pn-2, a0, •••, an-1) - искомый вектор коэффициентов регулятора (2); далее иногда будет удобно представлять его состоящим из двух векторов: R1 = col(Р0, •••, Pn-2) и R2 = col(a0, an1), тогда вектор R символически

записывать в виде R = col^R1, R2); D = col(d2n-1) - вектор, составленный из коэффициентов полиномов aet' {2n -1,5) и a (n, 5),

d =1 aI<-l, i е1,n -1,

I 011 ~0,-п , * е п,2п -1.

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

1) принцип исключения нуля;

2) методы теории Н.

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

Методы теории Н, представленные в работах [10-13] и др., позволяют формулировать достаточное условие робастной устойчивости при наличии структурно-параметрической неопределенности описания, однако существенным недостатком таких критериев является использование Н метрики операторов, что сильно ограничивает классы операторов, для которых выполняются условия этих критериев.

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

3 Далее для краткости изложения будем использовать слова робастные свойства или робастность, подразумевая под ними робастную устойчивость и робастное качество управления.

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

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

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

Е"(д, р) • А(п, р) у (í) = Е"(V, р) • В(т, р) и (г), п > т, д > V, (6)

здесь А(п, р), В(т, р), Е' (д, р), Е"(V, р) - семейства (множества) дифференциальных операторов:

п—1 ^ •

А(п, р) = а0 (п,р) + АА(п,р), а0 (п,р) = ^аг • р' + рп, А) =ео1(а0, •••, Оп—1, 1),

1=0

ДА(п,р) = <5а(п,5А) = ^5а; • р1': 5А=со1 (5а0, •••, 5ап) е Л сЯп+\ — А0 е а| ,

т •

В(т,р) = 6) (т,р) + ДВ(т,р), 6)(т,р) = Х6 • р' , В) = со1 ((, •••, Ьт),

1'=0

Г т •

АВ(п2, р) = <56 (т,5В) = ^56'. • р': 5В=со1(5Ь0, •••, 56т )е В с Ят+\ — В0 е В

I '=0

Е"(д,р) = ^0(д,р) + АЕ'(д,р), е0 (д,р) = 1 + • р', Е0 = со!(, •••, еЦ),

'=1

АЕ '(д, р) = ^5 е' (д, 5Е ') = £ 5 е" • р' : 5Е '=со! (, •, 5е\) е Ех с Ед, — Е0 е Е1,

Е"(V,р) = е0'(V,р) + АЕ>,р), е0"(V,р) = 1 + £>'• р' , Е0" = со!((, •, еу"),

'=1

АЕ "(V,р) = <5е''(V,5Е") = ]Г 5е'г р' : 5Е "=со!(5е;, •, 5еу") е Е2 с Я*, — Е" е Е21.

Далее предполагаем, что области A , B, Е1 и Е2 в (6) связны.

В записи (6) подразумевается, что пара {E' {q, p), E '{v, p)) описывает ту часть объекта

управления, для которой уже выполняются свойства устойчивости и качества управления, заданного областью S; т. е. предполагается, что для всех операторов из множеств E' (q,p) и E "(v, p) выполняются

Л(E') с int S , Л (E '') с int S, (7)

где Л^' ) и Л(E') обозначены множества нулей всех дифференциальных операторов E' (q,p) и E"(v,p) соответственно. Очевидно, регулировать следует только оставшуюся часть объекта управления (т. е. операторы A(n,p), B(m,p)); в этом смысле пара {E ' (q, p), E ''{v, p)) далее будет именоваться «структурными возмущениями», а {A (n, p), B (m, p)) - «основной динамикой» (подлежащей регулированию).

Для объекта (6) регулятор (2) (порядка l = n -1) будем рассчитывать по классической схеме синтеза, изложенной в предыдущем разделе; в качестве «точечного» объекта управления (1) примем P0:

ао (n, p) y (t) = b0(m, p) u (t),

где операторы а0 (п,р) и Ь0 (т,р) принадлежат основной динамике объекта (6). Замыкая объект (6) рассчитанным регулятором, получим ПФ:

Е>, 5) • В(т,5) -х(п -1,5)

ЖС'( 5) =--.

Е'(ц, 5) • Л(п, 5) • Р(п -1,5) - Е"(V, 5) • В(т, 5) ^а(п -1,5)

Знаменатель ПФ Жс1. (5) является характеристическим полиномом замкнутой системы:

Лс'. (2п + ц -1,5) = Е' (ц, 5) • А(п, 5) • Р(п -1,5) - E''(v, 5) • В(т, 5) • а(п -1,5).

С учетом определений операторов А(п, р), В(т, р), Е'(ц, р) и Е'(V, р) в (6) получим окончательное выражение для Ас'. ( 2п + ц -1,5) :

Лс1. (2п + ц -1,5) = |ас'. (2п + ц -1,5) = е'(ц,5Е') • (ае" (2п -1,5) + 5а(п,5Л) • р(п -1,5) --5 Ь (т, 5В )^а(п -1,5))+(е'( ц,5Е')- е'( V, 5Е "))• Ь (т, 5 )а(п -1,5):

У5Л е А, 5В е В, 5Е'е £1,5Е''е Е2} . (8)

Следуя статье [3], для семейства полиномов (8) введем новое множество функций, определенных следующим образом: каждый элемент ас1. (2п + ц -1,5) множества (8) разделим

на полином е'(д,5Е'), получим новый элемент

а с1(2п -1,5) =

ас1(2п + ц-1,5) е'( ц, 5Е') ,

или, в развернутом виде,

ас1 (2п -1,5) = ае. (2п -1,5) + 5 а (п, 5Л) • р(п -1,5) -

(

-5 Ь (т, 5В) • а( п -1,5) +

»Л

1 -

е" (V, 5Е") е' (ц,5Е '' )

Ь(т,5) ^а(п -1,5).

Множество элементов ас1 (2п -1,5) обозначим за Лс1 (2п -1,5) , таким образом: Лс1. (2п -1,5) = | ас1 (2п -1,5)=ае. (2п -1,5) + 5а (п, 5Л) -р(п -1,5) - 5 Ь (т, 5В) -а(п -1,5) +

е ' (V, 5Е'')

- ■ (Е')

Ь(т,5) ^а(п -1,5) : У5Л е А,5В е В,5Е ' е £1,5Е '' е Е2

(9)

Как следует из определения (9), семейство Лс' (2п -1,5) представляет собой множество дробных функций ас1 (2п -1,5) , каждая из которых имеет (в общем случае) 2п + ц -1 нулей и ц полюсов. В этом смысле будем говорить, что и семейство Лс1' (2п -1,5) имеет 2п + ц -1 нулей и ц полюсов. В соответствии с условием (7), все полюса дробных функций ас. ( 2п -1,5 ) лежат внутри области Б .

Пусть из множества функций Acl {2n -1,5) найдется хотя бы одна такая, что все ее нули лежат внутри области S, тогда потеря количества нулей функции ad. {2n -1,5) , находящихся внутри области S, происходит только в следующих случаях:

1) при вариациях векторов 5A, 5B, 5E ', 5E '' в (6) часть нулей и полюсов функции acC. (2n -1,5) в (9) совпадает (поскольку совпадающие нули и полюса находятся внутри области S, то их сокращение не повлияет на вариацию аргумента функции acl (2n -1,5) );

2) при переходе (хотя бы) одного нуля функции acl {2n -1,5) через границу области S .

Утверждение 1. Пусть

1) все множество нулей семейства полиномов E' (q, p) лежит внутри области S;

2) семейство функций (9) содержит хотя бы одну такую, что все ее нули лежат внутри S. Тогда для того, чтобы (все) множество нулей функций (9) лежало внутри области S, достаточно, чтобы выполнялось условие:

Р- >0,

где введены обозначения:

р-= min р {5 ), (10)

V sedS V '

Р(5) = Р1 (5)-Р 2 (5) , (11)

р, {5)= min aet(2n-1,5) + 5a(n,5A)-ß(n-1,5)-5b(m,5B)-a(n-1,5) ,

1V ' V5AeA, 5BeB I I

Р2 {5) = Ц{5)• b+{5)-|a(n -1, 5) |, e ''(v, 5E '' )

ц (5 )= max

v ' V8E"e£,, SE"e£.

1 --

e ' (q,5E')

Ъ, (s) = max I b0(m,s) + 5b(m,5B) I.

+ V ' V5BeSl U I

Доказательство данного утверждения следует из результатов, изложенных в статье [3].

Утверждение 1 формулирует достаточное условие робастной устойчивости и робастного качества управления для (8).

Функция р_ может служить в качестве оценки свойств робастного качества управления.

Расчет функции р(s) для каждой точки s границы области S сводится к решению задач поиска экстремумов функций с ограничениями на векторы варьируемых параметров, заданных областями A, B, E1 и E2. Указанная задача поиска условного экстремума функций решается методом неопределенных множителей Лагранжа.

В настоящей работе мы воспользуемся свободой в выборе расположения нулей эталона aet' (2n -1, s) (ограниченной условием (5)) с целью максимизации р-. Ниже излагается алгоритм поиска вектора R коэффициентов регулятора (2), максимизирующего 4 значение р-.

Основной результат.

Алгоритм повышения робастных свойств

Рассмотрим процедуру повышения р _ как итерационный процесс; обозначим за k номер итерации алгоритма. Для удобства изложения номер итерации будем указывать в качестве первого из аргументов векторов и функций; например, вектор R , полином aet' (2n — 1, s),

4 Имеется в виду локальный максимум функции р_ в пространстве векторов R

функции р(я), р_ и т. д. далее будут записаны как R(k), aet (k,2n -1, s), р(, s) и p_ (k), соответственно.

Вектору R (k) соответствует характеристический полином эталона aet (k,2n - 1,s), который определяется формулой

aet(k,2n -1, s) = a0 (n, s) -ß (k,n -1, s) - b0 (m, s) -a (k,n -1,s), (12)

где a0 (n, s) и b0 (m, s) - полиномы, принадлежащие основной динамике объекта управления (6). На вариации вектора R (k) будем накладывать условие

Л(aet) с intS , (13)

здесь за Л (aet) обозначено множество нулей полинома aeL (k,2n - 1,s) . Требование выполнения условия (13) на каждой итерации алгоритма является существенным, поскольку в противном случае нельзя гарантировать выполнения всех условий утверждения 1.

На начальной итерации (k = 0) алгоритма вектор R (0) рассчитывается по формуле (12) для заданного полинома aet' (0,2n -1, s) (выбираемого из условия (13)). Для R (0) по формулам (10) и (11) определяются р (0, s) и р- (0) .

На очередной k-й итерации алгоритма будем искать вектор R (k) такой, что:

1) после подстановки его коэффициентов в (12) для aet' (k,2n - 1,s) выполняется (13);

2) выполняется

р. (k)> р. (k -1) . (14)

Обоснование алгоритма

Обозначим за Q все множество точек s границы б S, в которых функция р (k -1, s) принимает минимальное значение:

Q = | s (i): s (i) = arg min р (k -1, s) j с б S .

Предполагая, что число Щ(к) таких точек на границе д Б на к-й итерации алгоритма конечно, множество & может быть переписано в виде

Q = | s (i) = arg min р (k -1, s) : V i e I == 1,

Для каждой точки 5 (/) (5 (7 )е &) найдем нормированный обобщенный 5 градиент функции р(к -1,5) по Ё (обозначим его за 0(/)) 6:

в (i ) =

фСЪМ) , если Ф1 (i), Ф2(i), (i)-^2 (i )\' Л) Л)

0 T, если Ф1 (i ) = Ф2 (i);

5 Слово «обобщенный» означает, что в тех случаях, когда градиент не определен, он доопределяется как нулевой вектор [14].

6 Здесь и далее под записью 0 понимается вектор-столбец соответствующей размерности, состоящий из ну-

Ф (' ) =

Ш1П

V5AеЛ, 5ВеВ

если Ш1П

V5AеЛ, 5ВеВ

( VЯе • Я + хЯе) • VЯе + ((1Ш • Я + х1ш ) • Vь

)

)2

1( (

VЯe • Я + тКе ) +( • Я +Т1Ш

Я = Я(к—1) =5(')

VЯе • Я + хЯе) + ( • Я + х1ш

Я = Я(к—1)

5=5" ('')

* 0;

0 , если ш1п

V5AеA, 5ВеВ

((Яе • Я + хЯе)2 + ((1Ш • Я + хЬп)2

Я = Я(к—1)

5=5 ('')

= 0;

ф2 (' ) =

0Т, 6+(5 )•

(ЁЯе • Я )• ЁЯе + (Ё1Ш • Я )• £1

Яе • Я)2 +(Ё1Ш • Я)2

Я=Я(к —1)

=5(' )

я2

если а

(п — 1,5)Я2 =Я2(к—1)* о, и ц(5) 5('* 0, и Ь+(5) | * 0;

5=5 К

0 Т е Я2 п—1,

если а

(п —^5 )|| Я =Я (к—1)= а или Ц (5 ) Ц.}= 0,

5=5 ('')

или

6+(5 ) 1=5(' )= 0

у7 =

Vя = (Яе (ц), •, Яе (^2п—1)), Vьп = ( 1ш (ц), •, 1ш (^))

ЁЯе = (Яе (§1), •, Яе &)), Ё1ш = ( 1ш (§), •, 1ш (§я)),

( а0 (п, 5) + 5 а (п, 5А)) • —1,

—(б0(т, 5) + 5 6(т, 5В) )• 5

7 е 1, п — 1, 7п 7 е п,2п — 1,

§ 7 = —\

7 е 1, п ,

хЯе = Яе

((а0(п, 5) + 5 а (п, 5А)) • 5п 1),

х1ш = 1ш

(( а0(п, 5) + 5 а (п, 5А)) • и" 1) .

Множество обобщенных градиентов & (''), вычисленных для всех точек множества О обозначим как 0 : 0 = | 9(''), V'' е 11 с Я2п—1 . Индексное множество I разделим на два под-

множества

11 ={': & (' )* 0, ' е I, & (') е 0 ), 12 = {': & (' ) = 0, ' е I, & (' )е 0 ).

В случае если 11 не пусто, из множества 0 формируем два множества

0+={ & ('):' е 11, & (') е 0 ), 0 ={ — & ('):' е 11, & (' )е 0 ).

Множества 0 + и 0— определяют направления градиентов и антиградиентов для функции р (к -1,5) в точках 5 множества О .

Пусть существует вектор ¥ (¥ е Я2п—1) такой, что выполняются условия

¥Т • & + (')>!

&

¥Т • &— ('') < —£0, V' е 11, (в0 > 0)

(15)

ми-

тогда вектор ¥ определяет направление возрастания функции р (к -1,5) во всех точках

нимума О (разумеется, за исключением тех точек, где градиент &(1) обращается в ноль).

Проблема построения (поиска) вектора ¥, обеспечивающего выполнение условий (15), может быть решена одним из алгоритмов построения «фазового портрета» (например, алгоритмом Гаусса - Зайделя 7, изложенным в [14]).

Если ¥ существует, вектор Я (к) будем искать в виде

Я(к) = Я(к -1) + у • ¥4, (7> 0), (16)

где за ¥ обозначен модуль вектора ¥ , у - параметр, определяющий величину шага алгоритма вдоль выбранного направления ¥ . Шаг у будем выбирать из условий (13) и (14) (к вопросу выбора такого у мы еще вернемся).

Случай 11 = 0 означает, что во всех точках множества О градиент функции р (к -1,5)

по Я либо не определен, либо равен нулю 8. В этом случае в (16) в качестве вектора ¥ последовательно выбираются векторы ортонормированного базиса в Я2"-1 . Условия останова алгоритма (16):

1) не существует (не найден) вектор ¥, являющийся решением задачи (15) 9;

2) выполнение условий достижения локального максимума функции р (к, 5) по параметру у хотя бы в одной из точек 5(1) множества О 10:

ё р1 ( к, 5 ) ёр2(к,5)]

ё у у=0 ё у

= 0.

у=0

5=5 (1)

ё2 р1 (к, 5 ) ё2р2(к,5 )

ё у2 у=0 ё У2 ]

< 0, 1 е 12.

у=0

5 =5 (1)

(17)

Здесь использованы обозначения

ё Р1 (к, 5 ) ё у

шт УЗЛеА, ЗБеБ Ф )

= < ) 7=0

у=0 V 5=5(1) )

5=5(1) 0,

если

если

ф)|т=0 * 0,

5=5(1)

ф)) 7=0 = °

5=5(1)

7 Данный алгоритм позволяет построить вектор ¥ (конечно, если он вообще существует), за число шагов, не превышающее значения 4/в0 (здесь е0 предполагается заданным); если алгоритм за указанное число шагов не находит искомый вектор, следовательно, его не существует.

8 Иначе говоря, в соответствующих точках 5 (5 е О ) функция р (к -1,5) достигла локального экстремума

по Я .

9 Вектор ¥, являющийся решением задачи (15), определяет направление возрастания функции р(к -1,5) общее для всех точек из множества О ; если такого вектора не существует, то Я (к -1) соответствует локальному

максимуму функции р_ по Я.

10 Очевидно, нет необходимости проверять выполнение условий, входящих в выражения (17), для всех точек множества О - подозрительными на локальный максимум являются те точки из множества О , в которых вектор обобщенного градиента & (1) обращается в ноль, т. е. такие точки 5(1), где 1 е 12.

ёР2(к,5) ё у

у=0

5 =5 (/)

Ц (5 )• ь+(5 )•

Ку )

если

у=0

5 =5 (/)

ю(у) 7=0 * 0,

5=5(1)

если ю (у) у=0 = 0,

п(у)=(фЯе • Я (к — 1) + ФЯе • ¥ • У^ + хЯе )• ФЯе • ¥ • у^ + +(ф 1Ш • Я (к — 1) + Ф 1ш • ¥ • У^ + хЬп )• Ф 1ш • ¥ • , ф)=(фЯе • Я (к — 1) + ФЯе •¥ • +хЯе) + (фЬп • Я (к — 1) + ФЬп •¥ • ^ + х

а(у )=( ЁЯе • Я2 (к — 1) + ЁЯе •¥ 2 • ^) ЁЯе •¥2 • ^ + +(ё1ш • Я2(к — 1) + Ё 1ш •¥, • ^Ё^ ¥ • ^, ¥2 = со1 (Уп—1, ..., У2п—1), ю(у) =(ЁЯе • Я22(к — 1) + ЁЯе • ¥2 • У^)2 + (Ё- • Я22(к — 1) + Ё1Ш • ¥2 • ^)2,

ё 2Р1 (к, 5)

ё у2

у=0

5 =5 (/)

Ш1П

V5AеЛ,

Т^у)

п ; (у )—

п(у ) л/Щ

2

у=0

5 =5(/ )

ё2Р2(к,5) ё у2

у=0 5=5(г)

ц (5 )• ь+(5 )

а; (у )—

Ку )

2

7=0 5 =5 (/)

если

если

если

если

§(у) 7=0 * 0,

5=5(г)

§(у) 7=0 = 0,

5=5 (г)

ю(у) 7=0 * 0,

5=5(г)

ю (у) 7=0 = 0,

5=5(г)

здесь п \ (у) и а ^ (у) - производные функций п(у) и а(у) по переменной у:

Ч (У) =(ФЯе • ¥ • X, )2 +(Ф - • ¥ • )2, а; (у) =( ЁЯе • ¥2 • )2 +(Ё1Ш • ¥2 •)/ )2.

Если ни одно из перечисленных условий останова не выполняется, определяется шаг у алгоритма (16). По смыслу поиск такого у (удовлетворяющего условиям (13) и (14)) состоит из 2-х этапов.

Этап 1. Определяется открытый интервал 11 (обозначим его за Г):

Г =[0, 7 *) , 7* > 0,

где за 7 * обозначено такое значение параметра у в (16), при котором нули полинома ае. (к,2п —1,5) в (12) выходят на границу области S . Искомое значение 7* (если оно существует) определяется формулой

7* = пип 7+ (5 (г)),

(18)

V5BеB

11 Данный интервал не может быть пустым (в силу выполнения условия (13) на к - 1-й итерации алгоритма).

здесь за у+ (5) обозначено положительное решение системы уравнений

Re (aet (k -1,2« -1, s)) + у • Re (Aaet (k -1,2« - 2, s)) = 0, Im(aet (k -1,2« -1, s)) + у • Im ( Aaet' (k -1,2« - 2, s)) = 0,

(19)

в (19) введено обозначение

Aaet (k - 1,2« - 2, s) = a0 («, s) •

( 1 n-1 Щ ^ VJ

V Щ J=1

• s

J-1

- bo (m, s) •

( 1 2«-1 Щ ^ ^ VJ

V J=«

за (у = 1,2« -1) обозначены компоненты вектора ¥ .

Решение переопределенной системы уравнений (19) относительно параметра у существует только при условии

ra«k

Re (A aet (k -1,2« - 2, s)) - Re (aet (k -1,2« -1, s)) Im (A aet' (k -1,2« - 2, s)) - Im (aet' (k -1,2« -1, s))

=1 ,

(20)

при этом значение у определяется выражением

-Ие (аг' (к -1,2« -1,5))

Y(s ) =

Re (Aaet'(k -1,2« - 2, s)) - Im (aet'(k -1,2« -1, s)) Im (A aet'(k -1,2« - 2, s))

если

если

Re (A aet'(k -1,2« - 2, s)) 0, Im (A aet(k ,2« - 2, s)) 0,

(21)

Если значение (21) отрицательно или такого значения вообще не существует 12, значит, не существует такого положительного значения у в (16), при котором нули полинома ае'. (к,2« -1,5) выходят на границу области Б в точке 5 (5 е О ). Таким образом, в качестве

у+ (5) в выражении (18) можно принять любое сколь угодно большое значение.

Прежде чем перейти ко второму этапу процедуры нахождения у отметим, что вектор ¥ (являющийся решением задачи (15)) соответствует общему направлению возрастания функции р (к -1,5) во всех точках 5 множества О (и, возможно, убывания в других точках границы д Б). Следовательно, в силу непрерывной зависимости функции р- (к) от вектора Я (к) (а значит, и от у), существуют такие значения у (у> 0), при которых функция р- (к) возрастает.

Этап 2. Внутри интервала Г ищется такой вложенный интервал с закрытыми границами [0, у], что для у выполняется условие (14). Значение у может быть найдено одним из методов

прямого поиска экстремума функции на заданном интервале Г.

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

12 Это является нарушением условия (20).

Алгоритм итеративного повышения робастности

Для работы алгоритма необходимо задать параметры алгоритма и выполнить следующие предварительные расчеты:

1) задать число Ы+, определяющее максимально допустимое число точек минимума функции р- на границе области Б, и число в0, определяющее величину «зазора» при решении задачи (15);

2) в качестве начальной итерации (к = 0) задать характеристический полином эталона ае. (2« -1,5) , удовлетворяющий условию (13), и по классической схеме синтеза модального

регулятора рассчитать вектор Я (0) .

Шаг 1. Расчет функций р (к, 5) и р- (к) по формулам (10), (11). Проверка условия

N (к )< М+, (22)

если данное условие выполнено, переход на шаг 2, в противном случае переход на шаг 9 с формулировкой «превышение допустимого числа точек минимума функции р (к, 5) на дБ ».

Шаг 2. Формирование множеств О , 0 и индексных множеств 11, 12.

Шаг 3. Проверка условия 11 = 0. Если данное условие выполнено - переход на шаг 5, в противном случае переход на шаг 4.

Шаг 4. В качестве ¥ выбирается любой ранее не выбиравшийся вектор ортонормирован-ного базиса в Я2«-1; выбранный вектор (вместе с векторами, выбранными на предыдущих итерациях) запоминаются. В том случае если такой вектор есть, переход на шаг 6; если таких

векторов нет, переход на шаг 9 с формулировкой «достигнут локальный максимум функции р- ».

Шаг 5. Формируются множества 0 +, 0- и решается задача (15). В случае если вектор ¥, являющийся решением задачи (15), не найден, переход на шаг 6; в противном случае -на шаг 9 с формулировкой «достигнут локальный максимум функции р- ».

Шаг 6. Вектор Я (к) ищется в виде выражения (16). Проверяются условия (17). Если условия (17) не выполнены, переход на шаг 7, в противном случае - переход на шаг 9 с формулировкой «достигнут локальный максимум функции р- ».

Шаг 7. По формулам (18) и (19) определяется интервал Г, внутри которого затем ищется такое значение у, при котором выполняется условие (14) (методами прямого поиска экстремума функции на заданном интервале). Если искомое значение у не найдено, переход на шаг 4.

Шаг 8. Принять к ^ к +1 и перейти на шаг 1.

Шаг 9. Останов.

Данный алгоритм обеспечивает нахождение такого вектора Я настроек регулятора (2), при котором функция (10) достигает своего локального максимума (при условии, что ни на одной из итераций алгоритма не нарушается условие (22).

Заключение

В настоящей статье разработан алгоритм вариации настроек модального регулятора (2), последовательно повышающий меру робастности (10) замкнутой системы. Алгоритм не требует точного знания операторов структурных возмущений - достаточно знания функции ц (5), определенной выражением (11). Данный алгоритм может быть реализован на ЭВМ.

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

Список литературы

1. Кузовков Н. Т. Модальное управление и наблюдающие устройства. М.: Машиностроение, 1976.

2. Соловьев И. Г. Методы мажоризации в анализе и синтезе адаптивных систем. Новосибирск: Наука, 1992.

3. Паршуков А. Н. Критерий робастной устойчивости и качества управления линейной замкнутой системой в условиях структурно-параметрической неопределенности описания в передаточной функции объекта управления // Вестн. НГУ. Серия: Информационные технологии. 2017. Т. 15, № 1. С. 59-69.

4. Ackermann J. Robust control: systems with uncertain physical parameters. London: Springer, 1993.

5. Barmish B. R. New tools for robustness of linear systems. New York: MacMillan, 1994.

6. Поляк Б. Т., Цыпкин Я. З. Частотные критерии робастной устойчивости и апериодичности линейных систем // АиТ. 1990. № 9. С. 45-55.

7. Киселев О. Н., Поляк Б. Т. Критический коэффициент усиления для последовательности неопределенных звеньев // АиТ. 1995. № 9. С. 93-104.

8. Поляк Б. Т., Цыпкин Я. З. Устойчивость и робастная устойчивость однотипных систем // АиТ. 1996. № 11. С. 91-104.

9. Шарый С. П. Новые характеризации множества решений для интервальных систем линейных уравнений // Вычислительные технологии. 2016. Т. 21, № 5. С. 111-118.

10. Francis B. A. A course in H control theory. Lect. Notes Control Inf. Sci., V. 88, Berlin: Springer, 1987.

11. Doyle J. C., Francis B. A., Tannenbaum A. Feedback control theory. New York: MacMillan, 1992.

12. Поляк Б. Т., Цыпкин Я. З. Робастный критерий Найквиста // АиТ. 1992. № 3. С. 25-31.

13. Киселев О. Н., Поляк Б. Т. Синтез регуляторов низкого порядка по критерию H и по критерию максимальной робастности // АиТ. 1999. № 3. С. 119-130.

14. Вапник В. Н., Червоненкис А. Я. Теория распознавания образов (статистические проблемы обучения). М.: Наука, 1974. 416 с.

Материал поступил в редколлегию 26.01.2018 O. F. Danilov, A. N. Parshukov, I. L. Polyanskaya

1 Tyumen Industrial University 38 Volodarsky Str., Tyumen, 625000, Russian Federation

2 Ural State Transport University (a branch in Tyumen) 5 Kalinin Str., Tyumen, 625000, Russian Federation

danilov_of@mail.ru, anparshukov@mail.ru,polyanskaya_il@inbox.ru

THE METHOD OF INCREASING OF ROBUSTNESS FOR LINEAR CLOSED-LOOP SYSTEMS WITH STRUCTURAL-PARAMETRIC UNCERTAINTY INTO TRANSFER FUNCTION OF OBJECT CONTROL

The article is established measure of robustness for linear closed-loop systems with indeterminacy of structure and parametric uncertainty into transfer function of object control. Also, in the article presented the algorithm of variation of parameters of a modal controller that increases the robustness for closed-loop system. The results may be realized on computer.

Keywords: robustness, modal controller, indeterminacy of structure, parametric uncertainty.

References

1. Kuzovkov N. T. Modalnoye upravlenie i nablyudayushie ustroistva. Moscow, Mashino-stroenie, 1976, 184 p. (In Russ.)

2. Solovev I. G. Metody majorizatchii v analize i sinteze adaptivnyh sistem. Novosibirsk, Nauka, 1992, 188 p. (In Russ.)

3. Parshukov A. N. Kriteriy robastnoy ustoytchivosti i kachestva upravleniya lineinoy zamk-nutoy sistemoy v usloviyah strukturno-parametricheskoy neopredelennosti opisaniya v pereda-tochnoy funktchii ob'ekta upravleniya. Vestnik NGU. Seriya: Informatchionnye tehnologii, 2017, vol. 15, no. 1, p. 59-69. (In Russ.)

4. Ackermann J. et al. Robust control: systems with uncertain physical parameters. London, Springer, 1993, 413 p.

5. Barmish B. R. New tools for robustness of linear systems. New York, MacMillan, 1994, 394 p.

6. Polyak B. T., Tchipkin Ya. Z. Chastotnye criterii robastnoy ustoichivosti i aperiodichnosti lineinyh sistem. AiT, 1990, № 9, p. 45-55. (In Russ.)

7. Kiselev O. N., Polyak B. T. Kriticheskii koeffitchient usileniya dlya posledovatelnosti neopre-delennih zven'ev. AiT, 1995, № 9, p. 93-104. (In Russ.)

8. Polyak B. T., Tchipkin Ya. Z. Ustoichivost' i robastnaya ustoichivost odnotipnih sistem. AiT, 1996, № 11, p. 91-104. (In Russ.)

9. Sharii S. P. Novie harakterizatchii mnozhestva reshenii dlya intervalnih sistem lineinih urav-nenyi. Vichislitel'nye tehnologii, 2016, vol. 21, no. 5, p. 111-118. (In Russ.)

10. Francis B. A. A course in H control theory. Lect. Notes Control Inf. Sci., V. 88, Berlin, Springer, 1987.

11. Doyle J. C., Francis B. A., Tannenbaum A. Feedback control theory. New York, MacMillan, 1992.

12. Polyak B. T., Tchipkin Ya. Z. Robastnyi kriterii Nyquist'a. AiT, 1992, № 3, p. 25-31. (In Russ.)

13. Kiselev O. N., Polyak B. T. Sintez regulyatorov nizkogo poryadka po kriteriy H i po kriteriy maksimal'noy robastnosti. AiT, 1999, № 3, p. 119-130. (In Russ.)

14. Vapnik V. N., Chervonenkis A. Ya. Teoriya raspoznavaniya obrazov (statisticheskie proble-my obucheniya). Moscow, Nauka, 1974, 416 p. (In Russ.)

For citation :

Danilov O. F., Parshukov A. N., Polyanskaya I. L. The Method of Increasing of Robustness for Linear Closed-Loop Systems with Structural-Parametric Uncertainty into Transfer Function of Object Control. Vestnik NSU. Series: Information Technologies, 2018, vol. 16, no. 1, p. 86-99. (In Russ.)