Исследование оптимальных по степени устойчиоости решений при пропорционально-интегральном управлении Текст научной статьи по специальности «Математика»

УДК 62-83.681.325.2

ИССЛЕДОВАНИЕ ОПТИМАЛЬНЫХ ПН СТЕПЕНИ УСТОЙЧИВОСТИ РЕШЕНИЙ ПРИ ПРОПОРЦИОНАЛЬНО-ИНТЕГРАЛЬНОМ УПРАВЛЕНИИ

А.М. Шубладзе (1), В.Е. Попадько (2), С.И. Кузнецов (3), A.A. Якушева (2)

(1) Институт проблем управления им. В.А. Трапезникова РАН, г. Москва;

(2) Российский государственный университет нефти и газа им. И.М. Губкина, г. Москва;

(3) ОАО «НИИТеплоприбор», г. Москва

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

ВВЕДЕНИЕ

Работа продолжает цикл статей, связанных с исследованием систем управления максимальной степени устойчивости. Синтезированные по этому критерию системы управления отличаются от других систем сравнительно высоким быстродействием и грубостью, робастностью к изменениям параметров объекта управления, так как наиболее «медленные», крайние правые корни характеристических полиномов таких систем максимально удалены от мнимой оси.

Один из возможных способов синтеза систем максимальной степени устойчивости основан на применении метода ^-разбиения, предложенного Ю.И. Неймарком [1, 2] и развитого в работах Б.Т. Поляка, П.С. Щербакова, Е.Н Грязиной [3, 4] и др. В случае ПИ-закона управления объектом с полиномом нулевой степени в числителе его передаточной функции указанный метод позволяет с помощью итерационных процедур сравнительно просто определить значение максимальной степени устойчивости и соответствующие ей настроечные параметры регулятора.

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

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

твовать системе, находящейся на границе устойчивости, и никакими вариациями настроечных параметров регулятора систему нельзя сделать устойчивой. Как следует из работ [5, 6], при ПИ-управлении динамическим объектом в оптимальном по степени устойчивости случае крайними правыми корнями характеристического полинома системы могут быть не меньше трех и не больше шести корней.

В настоящей работе выполнено исследование каждого из возможных случаев распределения корней и представлены достаточные условия их оптимальности. Таких случаев семь. Отличаются они числом действительных и комплексно сопряженных пар крайних правых корней характеристических полиномов и видом частотного годографа Михайлова характеристического полинома в комплексной плоскости, мнимая ось которой проходит через эти крайние правые корни.

Полученные оптимальные решения будут сформулированы в виде условий на коэффициенты оптимального полинома, определители Гурви-ца этого полинома и поведение годографа в окрестности точек пересечения им начала координат.

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

П

x(n)(t) + X аix('' - 1}(t) = k0u(t), (1)

i = 1

где x(i)(t), i = 0, ..., n, n l 3 — фазовые координаты, a¡ и k0 — параметры, u(t) — ПИ-управление:

t

u(t) = unn(t) = kns(t) + kM j"е(т)йТ, (2)

0

s(t) = g(t) — x(t), g(t) — задающее воздействие.

Оптимизируемыми являются параметры кП и кИ. Указанные параметры выбираются так, чтобы имела место максимальная степень устойчивости [6]

1оп = —minmaxReAj.(кп, ки), (3)

где Aj — корни характеристического полинома

n

Dn + i(A) = An + 1 + X a¿A¿ + ko(knA + ки). (4)

1

После замены А на А1 — / при кП = кПоп и кИ = кИоп в выражении (4) характеристический полином принимает вид [6]:

А(п + 1)(А1, кПоп кИоп) =

п Нк1

= А1^ П (А1 + юг ) ^1[(п + 1)-п^-2пк] (А1) =

I = 1

п2

= А1 +1 + X а* а1- 1 + X а* (кпоп, киоп) А'Г 1, (5)

1 = 3 1 = 1

где

а* = 1 д Р(п + 1)(А> кПоп, кИоп)

* (/ - 1) дА1 _1

при А = - /оп, (6)

^1[п + 1 - п _ 2пк] (А1) — полином, не имеющий корней в правой полуплоскости, п§ и у. — кратности

§ 1

нулевого корня и комплексно сопряженных пар корней оптимального решения, расположенных на мнимой оси. Число пар комплексно сопряженных корней определяется равенством

пк1

пк = X

1 = 1

Рассматриваемую задачу исследования оптимальных решений будем решать на основе результатов, полученных в работах [6, 7], с помощью машинных методов оценивания устойчивости систем управления. При этом нам понадобится следующий модифицированный годограф Михайлова из работы [6] характеристического полинома (5), в

2

котором А1 = /ю, / = -1:

Зй(/ю) = Ке^1(п + 1)(/ю) + ^^(п + 1)(/ю) =

= Яе0(ю2) + / 1т0(ю2). (7)

Как и обычный модифицированный годограф (7) будет соответствовать устойчивому полиному (5) (полиному с корнями в левой полуплоскости), если при увеличении частоты ю от нуля до бесконечности годограф последовательно проходит против часовой стрелки в п + 1 квадрант.

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

1.СЛУЧАЙ 1:

ДЕЙСТВИТЕЛЬНЫЕ КОРНИ ОПТИМАЛЬНОГО РЕШЕНИЯ

В этом случае имеют место условия

п§ = 3, пк = 0 (8)

при которых крайним правым корнем полиномов

(4) и (5) является один действительный корень кратности 3. Согласно работе [4] условия (8) могут выполняться только при

а* (/оп, кПоп, кИоп) = 0,

40) (/.п ^поп ки„„> > ^ (9)

где 1 < у < п — 2, / = 1, 2, 3, а* определяются

выражением (6), а)0) (/оп, кПоп, кИоп) — определители Гурвица полинома (5). При выполнении условий (9) достаточность условия оптимальности следует из того, что значение крайнего правого корня полинома (5) равно значению крайнего правого корня его второй производной, т. е. тахЯеА1г.(кПоп, кИоп) = тахЯеАц 2, / = 1, ..., п, где А1;. — корни уравнения ^1п(А1, кПоп, кИоп) = 0, а А1/ 2 — корни второй производной полинома (5)

д — I ( п + 1 ) ( А 1 , кПоп, кИоп ) = 0

дА 2 .

Любой выбор значений кПоп и кИоп в управлении (2), не удовлетворяющий условиям (9), приводит к уменьшению степени устойчивости / (3),

так как при этом в полиноме (5) а3 (/оп) = 0, а * (/оп, кП, ки) ф 0 хотя бы для одного а *, / = 1, 2, что и нарушает необходимое условие устойчивости полинома (5).

Рис. 1. Годограф оптимальной системы для случая 1

Как следует из выражений (5) и (9), минимальная степень при ю в Яе0(ю ) — см. формулу (7) —

равна четырем, а в 1т0(ю ) — двум, поэтому годограф (7) в рассматриваемом случае имеет вид, изображенный на рис. 1.

Он получен для полинома (5), в котором

^ 1 [ п + 1 - п £ - 2 п к ] (А1) = ^1[6](А1) = (А1 + 1) ,

п = 8, п§ = 3, пк = 0.

Из рис. 1 следует, что любое изменение значений кПоп или кИоп приводит к появлению корней полинома (5) в правой полуплоскости, т. е. к уменьшению степени устойчивости / системы (1), (2). Таким образом, вычисленные по формулам (6) и (9) значения кПоп и кИоп гарантируют в рассматриваемом случае максимальную степень устойчивости системы (1), (2).

2. СЛУЧАЙ 2:

ДВА ДЕЙСТВИТЕЛЬНЫХ КОРНЯ И ОДНА КОМПЛЕКСНО СОПРЯЖЕННАЯ ПАРА КОРНЕЙ

В этом случае имеют место условия

п§ = 2, пк = 1, у1 = 1, (10)

при которых крайними правыми корнями полиномов (4) и (5) являются один действительный корень кратности 2 и одна комплексно сопряженная пара корней. Согласно работе [6] условия (10) могут выполняться только при

а* (/оп, кПоп, кИоп) = 0, Ап- 2 (/оп, кПоп, кИоп),

А)0 * (/.п кПоп. киоп) > I0, (11)

где 1 < у < п — 3, а * определяются выражением (6),

/ = 1, 2, А)0) (/оп, кПоп, кИоп) — определители Гурвица размерности 1 < у < п — 2 полинома (5). Три уравнения, из которых находятся значения /оп, кПоп и кИоп, и (п — 3) неравенства (11) являются необходимыми условиями максимальной степени устойчивости системы (1), (2). Подтверждением этому служит годограф (7), изображенный на рис. 2 и полученный решением уравнений (11). Такой годограф не соответствует системе максимальной степени устойчивости, так как перемещением путем вариации параметров кП и ки начала координат комплексной плоскости А1 в любую точку области А получаем годограф (7), соответствующий устойчивой системе. Поэтому степень устойчивости (3) системы (2) в рассматриваемом случае может быть увеличена.

Рис. 2. Годограф не оптимальной системы для случая 2

Рис. 3. Годограф оптимальной системы для случая 2

Годограф, представленный на рис. 2, получен для полинома (5), в котором

^1 [п + 1 - - 2пк] (А1) = ^1[6](А1) = (А1 + ,

п = 8, п§ = 2, пк = 1, у1 = 1, ю1 = 0,08.

Для получения достаточных условий оптимальности нужно к условиям (11) добавить условие на направления годографа в точках пересечения им начала координат:

<(^(0), ^(0)), (Rel(ю2), 1т1(ю?))) < 0, (12)

где Re1(ю2) = дКе °1 (п+1)(гю),

5(ю 2)

1т1(ю2) = д( ю- 1 1 т Д1 1 +1 )(^ю)), (^(0), -Ке1(0)) -

д(ю )

вектор, перпендикулярный направлению годографа Zn0(/ю) при ю = 0, ^е 1( ю 1), 1т1( ю 1)) — вектор, перпендикулярный направлению годографа Zn0(/ю) при ю = ю1.

Из рис. 2 следует, что для годографа условие (12) не выполняется.

Покажем, что условия (11) и (12) в рассматриваемом случае являются достаточными условиями оптимальности (3). Как следует из условий (11), на мнимой оси плоскости А1 расположены нулевой корень кратности 2 и одна чисто мнимая комплексно сопряженная пара корней, остальные корни лежат левее мнимой оси. Для того, чтобы два нулевых корня изменением значений кП и кИ сдвинуть влево и чтобы при этом не появилась новая чисто мнимая комплексно сопряженная пара корней, необходимо начало координат плоскости А1 поместить в область между частью годографа, проведенного до первого пересечения им мнимой оси, и самой мнимой осью. На рис. 2 это область А, где точкой пересечения является начало координат. Нетрудно заметить, что при выполнении условий (12) годограф Zn0(/ю) (7) в окрестности частоты ю = ю1 не попадает в оговоренную выше область. Поэтому любое перемещение корней влево приводит к сдвигу комплексно сопряженной пары чисто мнимых корней вправо. Следовательно, сдвинуть путем изменения значений кП и кИ все рассматриваемые корни влево при выполнении условия (12) невозможно. Из этого следует, что условия (11) и (12) являются достаточными условиями оптимальности в случае 2.

На рис. 3 изображен годограф (7), соответствующий оптимальному по степени устойчивости решению и удовлетворяющий условиям (11) и (12).

Он получен для полинома (5), в котором

^ 1 [п + 1 - п?- 2пк] (А1) = ^1[6](А1) = (А1 + 1) ,

п = 9, п8 = 2, пк = 1, У1 = 1, ю1 = 0,6.

Из рис. 3 следует, любое изменение значений кПоп или кИоп приводит к неустойчивости полинома (5). Таким образом, выполнение условий (11) и (12) обеспечивает в рассматриваемом случае максимальную степень устойчивости системы (1), (2).

3. СЛУЧАЙ 3:

ОДИН ДЕЙСТВИТЕЛЬНЫЙ КОРЕНЬ И ОДНА КОМПЛЕКСНО СОПРЯЖЕННАЯ ПАРА КОРНЕЙ

В этом случае имеют место условия

\ = Ь пк = 1 Уі = Ь (13)

при которых крайними правыми корнями полиномов (4) и (5) являются один действительный корень и одна комплексно сопряженная пара корней. Согласно работе [6] условия (13) могут выполняться только при

а* (1оп, кПоп, кИоп) = 0, 4п - 1 (1оп, кПоп, кИоп) = 0,

40) (^оп. *Поп. кИоп) > 0- (14)

где 1 < j < п — 2, а * — определяется выражением (6),

40) (1оп, кПоп, кИоп) — определители Гурвица размерности 1 < j < п — 1 полинома (5). Двух уравнений (14) недостаточно для нахождения трех неизвестных /оп, кПоп и кИоп. Третье недостающее уравнение может быть получено из соображений о невозможности изменения только значения кПоп в управлении (2) невозможно увеличить степень устойчивости системы (1), (2). В рассматриваемом случае такая ситуация будет, когда годограф (7) полинома (5) касается в начале координат мнимой оси. Такой факт имеет место при выполнении

условия Яе1( ю1) = 0 из которого следует, что полином

В1 п1 (ю2) = 1т0(ю2) + /юЯе1(ю2), (15)

где п1 = 2[(п + 1)/2] — 1, [•]— целая часть числа, имеет ту же комплексно сопряженную пару корней, что и исходный полином (5). Поэтому в оптимальном случае справедливо уравнение

<- 1 ('ол кПоп) = ^ (16)

где 4 п0)_ 1 (Т^, кПоп) — определитель Гурвица размерности п1 — 1 полинома (15).

Рис. 4. Годограф неоптимальной системы для случая 3

Рис. 5. Годограф оптимальной системы для случая 3

Из условий (14) и (16) можно определить три неизвестных параметра /оп, кПоп и кИоп. Но условия

(14) и (16), как и в предыдущем случае условий (11), являются необходимыми условиями оптимальности. Подтверждается это годографом (7) полинома (5), изображенном на рис. 4, для которого условия (14) и (16) выполнены, но очевидно, что увеличением параметров кП и кИ полином (6) может быть сделан устойчивым.

Годограф рис. 4 получен для полинома (5), в котором

В1 [п + 1 - п?-2пк] (А1) = В1[6](А1) = (А1 + 1) ,

п = 8, п% = 1, пк = 1, у1 = 1, ю1 = 0,32.

Достаточные условия оптимальности можно получить добавлением к условиям (14) и (16) условия на направления годографа в точке пересечения им на частоте ю1 начала координат. Это направление должно быть противоположно направлению мнимой оси, что имеет место при выполнении неравенства

1шД ю 1) < 0, (17)

где функция 1т1(ю ) определена выражением (12).

Условия (14), (16) и (17) являются достаточными условиями оптимальности в рассматриваемом случае. На рис. 5 изображен годограф (7), соответствующий оптимальному по степени устойчивости решению и удовлетворяющий условиям (14),

(16) и (17).

Он получен для полинома (5), в котором

В1 [п + 1 - п,-2пк] (А1) = В1[6](А1) = (А1 + ^

п = 8, п^ = 1, пк = 1, ю1 = 0,72.

4. СЛУЧАЙ 4:

ОДИН ДЕЙСТВИТЕЛЬНЫЙ КОРЕНЬ И ДВЕ ПАРЫ КОМПЛЕКСНО СОПРЯЖЕННЫХ КОРНЕЙ

В этом случае имеют место условия

п§ = 1, пк = 2, у1 = 1, у2 = 2, (18)

при которых крайними правыми корнями полиномов (4) и (5) являются один действительный корень и две комплексно сопряженные пары корней. Согласно работе [6] условия (18) могут выполняться только при

а* (1оп, кПоп, кИоп) = 0,

4п- 3 (1оп, кПоп, кИоп) = 4п~) 3 (1оп, кПоп, кИоп) = 0,

4)0) кИоп) > (I9)

где 1 < j < п — 4, а * определяется выражением (6),

4п- 3 (1оп, кПоп, кИоп), 4п- 3 (1оп, кПоп, кИоп) — основной и побочный определители Гурвица полинома (5), описанные в книге [7]. Три уравнения, из которых находятся /оп, кПоп и кИоп и (п — 4) неравенства в условиях (19) являются, как и в рассмотренных случаях 2 и 3, необходимыми условиями максимальной степени устойчивости системы (1), (2). На рис. 6 изображен годограф полинома (5), построенного из условий (19), но который не является годографом оптимальной системы. Очевидно, что увеличением кПоп и уменьшением кИоп полином (5) может быть сделан устойчивым.

Рис. 6. Годограф неоптимальной системы для случая 4

Рис. 7. Годограф оптимальной системы для случая 4

Годограф получен для полинома (5), в котором

В1 [п + 1 - п5- 2пк] (А1) = В1[6](А1) = (А1 + ^ п = 10,

= 1, пк = 2, У1 = у2 = 1, ю1 = 0,46, ю2 = 0,65.

Оптимальность такой системы имеет место при выполнении следующих неравенств:

Яе1( ю2) > 0, Яе1( ю 2) < 0,

Яе1( ю2) < 0, Яе1( ю 2) < 0,

где ю1 < ю2 — частоты, на которых годограф пересекает начало координат, функция Яе1(ю ) определена при рассмотрении случая 2. На рис. 7 изображен годограф полинома (5) оптимальной системы. Он получен для полинома (5), в котором

В1 [п + 1 - п^-2пк] (А1) = В1[6](А1^ = (А1 + ^) , п = 10,

п^ 1, пк 2, У1 У2 1 ю1 0,3, ю2 0,5.

5. СЛУЧАЙ 5:

ДВЕ ПАРЫ КОМПЛЕКСНО СОПРЯЖЕННЫХ КОРНЕЙ

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

пк = 2, у1 = 1, у2 = 2, (20)

при которых крайними правыми корнями полиномов (5) и (6) являются две комплексно сопряженные пары корней. Условия (20) могут выполняться только при

4п - 2 (1оп, кПоп, кИоп) = 4^7-) 2 (1оп, кПоп, кИоп) = 0,

4у°) кИоп) > 0. (21)

где 1 < j < п - 3. 4п-3 (^ ^оП ‘ноп^ 4п_ 3 (^.П кПоп, кИоп) — основной и побочный определители Гурвица. Из двух уравнений (21) невозможно найти три параметра /оп, кПоп и кИоп. Для получения необходимого третьего уравнения воспользуемся результатами работы [6], из которой следует, что при условиях (20) имеет место касание годографа в точках пересечения им начала координат на частотах ю1 и ю2. Поэтому существует полином

^(ю2) = vRe1(ю2) + 1т1(ю2),

где V — тангенс угла наклона касательной, корнями которого также являются ю = ю1 и ю = ю2, в силу чего полином от /ю

Вп0у(/ю) = 1т0(ю2) + /юS1>2), п0 = ^/2^ (22)

2

где 1т0(ю ) из выражения (7), имеет корни при тех же ю = ± ю1 и ю = ± ю2. В связи с этим оказываются справедливыми два следующих уравнения:

4!?.3.« ('оп кПоП V) - 4п0-3,, (^ V) = 0,(23)

где 4.3. ^ ('оп, *Поп, V), ' = 0, 1 — основной и побочный определители Гурвица полинома (22).

Таким образом, мы имеем четыре уравнения (21) и (23) с четырьмя неизвестными /оп, кПоп, кИоп и V. Как и в предыдущем случае, условия (21) и (23) являются необходимыми условиями максимальной степени устойчивости системы (1), (2). На

Рис. 8. Годограф неоптимальной системы для случая 5

Рис. 9. Годограф оптимальной системы для случая 5

рис. 8 изображен годограф полинома (6), построенного из условий (21) и (23), но который не является годографом оптимальной системы. Он получен для полинома (5), в котором

В [п +1 - пг- 2пк] (А1) = ВИ6](А1) = (А1 + 1)6, п = 9, п§ = 0, пк = 2, у 1 = у2 = 1, ю1 = 0,03, ю2 = 0,23.

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

<^е^ ю 1), 1т1( ю 1)), ^еД ю 2), 1т1( ю2))) < 0,

22

в котором функции Re1(ю ) и 1т1(ю ) определены

2 2

ранее (в условиях (12) и (17)), ^е^ю2), 1т1(ю)) и

22

(Re1( ю 2), 1т1( ю 2)) — векторы направления годографа ^^(/ю) при ю = ю1 и ю = ю2 соответственно, ю1 и ю2 — частоты, на которых годограф (7) проходит через начало координат. На рис. 9 изображен годограф полинома (5) оптимальной системы. Любое изменение положения этого годографа в направлении касательного вектора путем вариаций значений кПоп и кИоп приводит к переходу двух пар крайних правых корней в правую полуплоскость. Изменение в направлении, перпендикулярном этому вектору, приводит к переходу одной из пар в правую полуплоскость.

Годограф получен для полинома (5), в котором

В1 [п + 1 -пг-2пк] (А1) = В1[6](А1) = (А1 + ^

п = 9, п^ = 0, пк = 2, ю1 = 0,05, ю2 = 0,41.

1т х10 3

6. СЛУЧАЙ 6:

ОДНА КРАТНАЯ ПАРА КОМПЛЕКСНО СОПРЯЖЕННЫХ КОРНЕЙ

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

пк = 1, У1 = 2, (24)

при которых крайними правыми корнями полиномов (5) и (6) является кратная комплексно сопряженная пара корней. Условия (24) могут выполняться при

4п- 2 (1оп, кПоп, кИоп) = 4п - 2 (1оп, кПоп, кИоп) = 0,

4п0-2 (/„ кИоп) > 0. (25>

где 1 < j < п - 3. 4п0) 2 (/оп. ^По^ кИоп) и 4п1- 2

кПоп, кИоп) — основной и побочный определители Гурвица. Для нахождения трех параметров /оп, кПоп и кИоп, как и в работе [7], введем полином

вЦп -1) (/оп) = Rel(ю2) + /юГт^ю2), (26)

22

где функции Re1(ю ) и 1т1(ю ) определены ранее. Полином (26) имеет корень на той же частоте, что и полином (5). Поэтому справедливо условие

А 10п - 2) (/оп) = 0, (27)

где 41°п - 2) (/оп) — основной определитель Гурвица полинома (26). Трех уравнений (25) и (27) достаточно для определения параметров /оп, кПоп и кИоп,

Рис. 10. Годограф оптимальной системы для случая 6

причем условия (25) и (27), как и в случае 1, являются достаточными условиями оптимальности. Вначале из уравнения (27) находится оптимальное значение степени устойчивости /оп, затем из уравнений (25) — параметры кПоп и кИоп. На рис. 10 изображен годограф полинома (5) оптимальной системы. Он получен для полинома (5), в котором

В1 [п + 1 -п§-2пк] (А1) = В1[6](А1) = (А1 + 1) ,

п = 9, п, = 0, пк = 1, у1 = 2, ю1 = 0,2.

7. СЛУЧАЙ 7:

ТРИ ПАРЫ КОМПЛЕКСНО СОПРЯЖЕННЫХ КОРНЕЙ

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

пк = 3 У1 = ^2 = ^3 = 1 (28)

при которых крайними правыми корнями полиномов (4) и (5) являются три пары комплексно-сопряженных корней. Условия (28) могут выполняться при

4п-4 (/оп, кПоп, кИоп) = 4п-) 4 (/оп, кПоп, кИоп) =

= 4- 4 (/оп, кПоп, кИоп) = 0,

4)0) (/оп кПоп кИоп) > 0, (29)

где 1 < j < п - 5, 4п0-4 (/оп, кПоп кИоп) = 41-4 (/оп,

кПоп кИоп) = 42 4 (/оп, кПоп кИоп) — основной и ПО-

бочные определители Гурвица. Условия (29) являются необходимыми условиями максимальной степени устойчивости системы (1), (2). На рис. 11 изображен годограф полинома (5), построенного при выполнении условий (29), но который не является годографом оптимальной системы.

Он получен для полинома (5), в котором

В1 [п + 1 -п^-2пк] (А1) = В1[6](А1) = (А1 + 1) ,

п = п, п, = 0, пк = 3, У1 = у2 = у3 = 1, ю1 = 0,03, ю2 = 0,15, ю3 = 0,254.

Для оптимальности системы (1), (2) достаточно выполнение следующих условий на скалярные

произведения. Если ^е^ ю 1), 1т1( ю 1)) — направление годографа (7) в точке пересечения начала

22

координат на частоте ю1, ^е^ ю 2), 1т1( ю 2)) —

2

направление годографа на частоте ю2, (Re1( ю 3),

2

1т1( ю 3)) — направление годографа на частоте ю3,

Рис. 11. Годограф неоптимальной системы для случая 7

Рис. 12. Годограф оптимальной системы для случая 7

2 2

где Яе1(ш ) и 1ш1(ш ) определяются выражением

(15), то произведение скалярных произведений

<(-Іш1( ю2), Яе1( ю2)), (Яе1( ), 1ш1( ))> х

<(-ІШі( ю 1), Яеі( ю2)), (Яеі( ю 3), ІШі( ю 3 ))> < 0.(30)

Доказательство оптимальности условий (29) и (30) основано на том, что при сдвиге путем изменения значений кПоп и кИоп любых двух чисто мни-

мых комплексно сопряженных пар корней влево от мнимой оси плоскости А1 третья пара сдвигается вправо. При выполнении условия (38) отсутствует общая область устойчивости, аналогичная области В, которая была в случае, показанном на рис. 11.

На рис. 12 изображен годограф полинома (5) оптимальной системы. Он получен для полинома

(5), в котором

D1 [n + 1 - ng-2nk] (Ai) D1[6](Al) (A1 + !)

11, ng = 0, nk = 3, у = 1,

6

g

k

1 = 0,1, ш2 = 0,4, ш3 = 0,8.

ЗАКЛЮЧЕНИЕ

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

ЛИТЕРАТУРА

1. Неймарк Ю.И. Устойчивость линеаризованных систем — Л.: ЛКВВИА, 1949.

2. Неймарк Ю.И. Динамические системы и управляемые процессы. — М.: Наука, 1978.

3. Поляк Б.Т., Щербаков П.С. Робастная устойчивость и управление. — М.: Наука, 2002.

4. Грязина Е.Н. К теории D-разбиения // Автоматика и телемеханика. — 2004. — № 12. — С. 15—28.

5. Шубладзе А.М. Достаточные условия экстремума в системах максимальной степени устойчивости. I // Автоматика и телемеханика. — 1997. — № 3. — С. 93—105.

6. Шубладзе А.М. Достаточные условия экстремума в системах максимальной степени устойчивости. II // Автоматика и телемеханика. — 1997. — № 8. — С. 67—79.

7. Гантмахер Ф.Р. Теория матриц. — М.: Наука, 1967.

Ш (495) 334-88-81, e-mail: shub@ipu.ru

Статья представлена к публикации членом редколлегии А.Г. Бутковским. □

Читайте в следующем номере статью

Цыганова В.В., Кадымова Д.С. Механизмы управления политико-экономическими сис-