Многокритериальная параметрическая оптимизация трехканальной системы стабилизации летательного аппарата с перекрестными связями Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 681.513.1

МНОГОКРИТЕРИАЛЬНАЯ ПАРАМЕТРИЧЕСКАЯ ОПТИМИЗАЦИЯ ТРЕХКАНАЛЬНОЙ СИСТЕМЫ СТАБИЛИЗАЦИИ ЛЕТАТЕЛЬНОГО АППАРАТА С ПЕРЕКРЕСТНЫМИ СВЯЗЯМИ

А.С. Аксенов1, Е.М. Воронов2, К.К. Любавский2, С.И. Сычев1

1ОАО "Корпорация "Тактическое ракетное вооружение", Королев Московская обл., Российская Федерация e-mail: [email protected]

2МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: [email protected]

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

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

MULTI-CRITERIA PARAMETRIC OPTIMIZATION

OF THE TRIPLE-CHANNEL CROSS-COUPLING STABILIZING SYSTEM OF AN AIRCRAFT

АА Aksenov1, E.M. Voronov2, K.K. Lyubavskiy2, S.I. Sychev1

1OAO "Corporation "Tactical Missile Armament", Korolev, Moscow region,

Russian Federation

e-mail: [email protected]

2Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: [email protected]

Methods for optimizing control of multi-object and multi-criteria systems on the basis of stable-effective compromise concept are used in multi-criteria task to parametrically optimize the structurally-complex triple-channel stabilizing system of a statically stable aircraft with normal aerodynamic design under conditions of initial structural inconsistence of channels based on cross coupling between channels. The generalized stable-effective algorithm of multi-criteria optimization is formed on the basis of equilibrium-arbitral compromise. A complete set of cross couplings is revealed in the practically useful model of linearized stabilizing system of an anti-ship missile with "perturbed" reference trajectory of flight. Optimal modes of aircraft

stabilization are investigated taking into account the stability, performance, static accuracy and control response under exposure to perturbations and with changes in the cross-coupling conditions.

Keywords: multi-object and multi-criteria system, stable-effective compromise, reference trajectory, equilibrium-arbitral algorithm, Nash equilibrium, stabilizing system, criterion of efficiency, pattern search.

В настоящей работе развивается подход равновесно-арбитражной балансировки каналов в многосвязном регулировании и управлении [1] для системы стабилизации летательного аппарата (ЛА) на основе теории оптимизации управления многообъектными многокритериальными системами (ММС) [2], которая применяется для трех классов задач управления структурно-сложными системами в условиях исходной структурной несогласованности, конфликтной ситуации и неопределенности (среды, "активного" партнера, цели). Рассматриваемая задача относится к первому классу и обобщает результаты для системы стабилизации (ССт) ЛА, полученные в работе [1].

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

Математические модели каналов ССт и трехканальной ССт с перекрестными связями. Рассмотрено динамическое описание типовой противокорабельной ракеты. Структурные схемы контура стабилизации продольного и бокового каналов, а также канала крена [3-5] в виде соединения передаточных функций (ПФ) ЛА, рулевого привода, различных фильтров, регуляторов и цифровых делителей напряжения (ЦДН) приведены на рис. 1.

Управляющими параметрами продольного канала являются величины КШг, КПу, управляющими параметрами бокового канала — величины КШу, Knz, управляющими параметрами канала крена — КШх, KY.

Структурная схема трехканальной ССт с перекрестными связями представлена на рис. 2.

Формирование математического описания показателей и постановка задачи комбинированной оптимизации трехканальной ММС с перекрестными связями. Оценить качество регулирования можно непосредственно по расчетным или опытным переходным процессам или косвенно по другим динамическим свойствам [6, 7]. Таким образом, методы оценки качества систем регулирования подразделяют на прямые и косвенные. К первым относят определение характеристик переходного процесса, времени выхода на установившийся режим,

от от 2

63

о

3

Я

Ш СП

к

ё р

о

о

Я

ч к

СП

о ч

о

о £

о о я к

о

/ ±15

К,..

+ 18

к 1

и 0,05« + 1

0.015^+ 1

0,005^ + 1

+ 5

К,

\ Блок 2

к\ К2

0,332

Л/,

Фильтр 1 1

.2 2

0,005 " 5 " + 2 • 0,005 ■ 0,6 5 + 1

8,33

ЦДН«/,

т

РП

1

0,015+ 1

+ 30

А"

| Фильтр 2

1

0,005.? + 1

0,04

1

.2 2

0,005"«"+ 2 ■ 0,005- 0,5 л- + 1

со,,

57,3

Фильтр 3

1

0,03 2 .V2 + 2 ■ 0,4 • 0,03 5 + 1

V (Т35 + 1)(Т45 + 1)

А

Г, 5+1

Цз2 + 2Т2^+ 1

V________________у

<

^возм

Рис. 1 (начало). Структурные схемы контура стабилизации продольного (а) и бокового (б) каналов

О

ОТ ОТ

О

63

3

Я

и

СП

к

ё р

о >8

£

к

СП

о ч

о

43 о о X 8

О

±15

К8у г

£

0,332

М„

Блок 3

8,33

ЦДНу

РП

1

0,01 ^ + 1

+ 30

Ку

Туй + 1

к

57,3

А"

Фильтр 2

57,3

1 0,04 1 1 1 Шд

0,005 5 + 1 1 1 0,005252+2-0,005-0,5^+ 1

Рис. 1 (окончание). Структурная схема контура стабилизации канала крена (в):

М:возы, Мувози, Мхвози — возмущающие моменты, действующие на ЛА (относительно осей Z, У и X); 5В — отклонение руля высоты в продольном канале; И^ — ПФ ЛА изменения угловой скорости тангажа по углу 5Ъ; Т,1,Т,2,ц,/\ — коэффициенты этой ПФ, выражаемые через динамические коэффициенты а,: ,; — ПФ ЛА изменения перегрузки по угловой скорости тангажа; — коэффициен-

ты этой ПФ, выражаемые через динамические коэффициенты сц Кш:, Кп:, К&п:, К2 — коэффициенты регуляторов, зависящие от скоростного напора д; 5Н — отклонение руля направления в боковом канале; — ПФ ЛА изменения угловой скорости рыскания по возмущающему воздействию; Тх, Т2, ^, — коэффициенты этой ПФ, выражаемые через динамические коэффициенты 6г: , ; И''*'* — ПФ ЛА изменения перегрузки по угловой скорости рыскания; Тз, Т4 — коэффициенты этой ПФ, выражаемые через динамические коэффициенты Кшу,КПу,К&Пу — коэффициенты регуляторов, зависящие от скоростного напора д; дэ — отклонение элеронов в канале крена; ИЛ? — ПФ ЛА изменения угловой скорости крена по углу 6Э; КЭ,ТЭ — коэффициенты этой ПФ, выражаемые через динамические коэффициенты с^у, КШг, КШу, А 7, — коэффициенты регуляторов, зависящие от скоростного напора д; 1-2 — контуры по перегрузкам и углу крена

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

Далее использованы прямые методы оценки качества и установившейся точности по переходному процессу при реакции на типовое воздействие (оценку точности или статической ошибки) и косвенные частотные методы (определение границы и запасов устойчивости с помощью построения Б-разбиения в пространстве параметров системы для двухканальной системы [1] и экспериментальным образом для трехканальной ССт).

В работе рассмотрен вариант исследования ССт в режиме стабилизации: входом является возмущающий момент, а выходом — угловые скорости тангажа, рыскания и угол крена. Основная задача ССт в указанном режиме — парирование внешних возмущений, т.е. стабилизация углового положения ЛА (обеспечение заданного углового положения или угловой скорости) при воздействии возмущений. Поэтому в качестве одного из критериев выбрана статическая точность, которая характеризуется статической ошибкой.

Критерий статической точности. Статическая ошибка определяется по установившейся реакции на типовое ступенчатое воздействие. Требуется, чтобы статическая ошибка была не выше ошибки при отдельном расчете параметров каждого канала.

При ступенчатом воздействии на продольный и боковой каналы, а также на канал крена получаем следующие возмущающие моменты, Н м: МгВозм = 52620, Муъозм = 62212 и Ыхвозы = 10941. Требования к статической ошибке по перегрузке в продольном, боковом каналах и канале крена задаются тремя неравенствами вида

еПг < 0,003; еПу < 0,003; е7 < 0,05 °/с.

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

Jl^& (^уст пх ^треб пх) ^ ШШ;

ф (^уст пу ^треб пу) ^ Ш1П; (1)

^ (^уст 7 ^треб 7 ) ^ Ш1П .

Критерий колебательности ССт. Демпфирующие свойства автоматической системы оцениваются перерегулированием о, которое

Рис. 2. Структурная схема трехканальной ССт с перекрестными связями:

Wдyc — ПФ фильтра инерциального датчика угловой скорости (ДУС); блок 1,2,3 — блоки для обеспечения свойств устойчивости, качества, элементов адаптации во внешних контурах по перегрузкам (1,2) и углу крена (3), выделенные на рис. 1; РП — блоки соответствующих рулевых приводов, содержащих в общем случае нелинейные свойства вследствие люфта редуктора и сухого трения; , ^о, То — значения углов тангажа, рыскания, крена на опорной траектории, относительно которой проводилась линеаризация уравнений движения ЛА

представляет собой нормированное динамическое отклонение регулируемой величины [6, 7]. Перерегулирование определяется по реакции нормальной перегрузки на ступенчатое воздействие в продольном и боковом каналах и угла крена в канале крена:

а = Птах - Пуст • 100%;

пуст

а = ^тах - 7уст • юо%.

Туст

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

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

= (апг — атреб пг )2 + (Пп2 — Птребпг )2 ^ Ш1П;

¿2'ф = (аПу - атреб О + Ыпу - Птреб пу) ^ Ш1П; (2)

22

j2y (^Y ^треб y ) + (nY ^треб y ) ^

min,

где атребг = 0,2 (20%).

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

Критерий быстродействия системы стабилизации. Время регулирования ¿р — это минимальное время, по истечении которого (с момента подачи ступенчатого воздействия) отклонение выходной величины от установившегося значения не превышает некоторого заданного значения А. Математически такое определение можно записать следующим образом:

¿р = т1п{Тр : |к (¿) - к (то) | < А, ¿ > Тр}.

тр

Как правило, А = 0,05.

Минимизируемые показатели сформированы в виде квадратичных невязок времени регулирования и требуемого значения:

Jы = (¿п2 — ¿треб пг )2 ^ т1п;

^ф = (¿пу — ¿треб пу) ^ т1п; (3)

J3Y = (¿7 — ¿треб 7) ^ т1П,

где ¿треб — время, сформированное в соответствии с результатами моделирования на множестве управляемых параметров.

Следует отметить, что оценку времени переходного процесса для двухканальной ССт можно выполнить аналитически [1, 6], но в силу сложности математического описания трехканальной системы требуемое время переходного процесса предложено оценивать по результатам многофакторного моделирования с учетом предварительной оценки значения времени в усеченной двухканальной ССт с перекрестными связями.

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

J4& = {КШг - KTреб)2 + {Knz - Кребnz )2 ^ min;

,1А>ф = (КШу - КтребШу) + (Kny - Ктребnj ^ min; (4)

J47 = (КШх - Ктреб Wx )2 + (К7 - Ктреб y)2 ^ min .

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

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

Метод Хука-Дживса предназначен для поиска безусловного локального экстремума функции и относится к прямым методам, т.е. опирается непосредственно на значения функции. Алгоритм делится на две фазы: 1) исследующий поиск; 2) поиск по образцу [10].

Рис. 3. Иллюстрация первой (а) и второй (б) фаз алгоритма для двух координат

1. Исследующий поиск. На начальном этапе задается стартовая точка (точка 1) и шаги к по координатам (рис. 3, а). Затем замораживаются значения всех координат кроме первой, вычисляются значения функции в точках хх + кх и хх — кх (хх — первая координата точки; кх — значение шага по этой координате) и выбирается точка с наименьшим значением функции. В этой точке замораживаются значения всех координат кроме второй, рассчитываются значения функции в точках

+ к2 и — к2, выбирается точка с наименьшим значением функции и т.д. Когда шаги к по всем координатам станут меньше соответствующих значений точности ^, алгоритм завершается и точка признается точкой минимума. Таким образом, проведением исследующего поиска по всем координатам была получена новая точка с наименьшим значением функции в окрестности (точка 2). Теперь можно осуществить переход ко второй фазе алгоритма.

2. Поиск по образцу. На начальном этапе поиска по образцу откладывается точка 3 в направлении от точки 1 к точке 2 на том же расстоянии. Ее координаты получают по формуле

Хз = Х1 + Л(х2 - Х1), (5)

где х — точка с номером г; Л — параметр алгоритма, принимающийся равным 2. Затем в точке 3 проводится исследующий поиск, как на первой фазе алгоритма, за исключением того, что шаг на этой фазе не уменьшается.

Если на указанной фазе в результате исследующего поиска удалось получить точку 4 (рис.3,б), отличную от точки 3, то точка 2 становится точкой 1, а точка 4 — точкой 2 и повторяется поиск по образцу. Если не удается найти точку 4, отличную от точки 3, то точка 2 переобозначается как точка 1 и повторяется первая фаза алгоритма — исследующий поиск.

С учетом специфики поставленной задачи блок-схема алгоритма поиска точки для формирования критерия устойчивости будет иметь вид, представленный на рис. 4 (показано для четырех параметров).

Критерий оптимизации в каждом канале рассматриваемой ССт является векторным и предъявляет требования к точности, демпфи-

Рис. 4. Блок-схема первой фазы алгоритма — исследующий поиск (справа показана процедура получения интервала устойчивости для г-го параметра)

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

■ Ф^ ■ Ф1

Ф = Ф^ = Ф2

. Ф^ . Ф3

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

4 4

Ф1 = Ф^ = J J-r& = J Jji ^ min;

j=i j=i 44

Ф2 = Ф^ = j J = ^2 aj2 Jj2 ^ min;

j=1 j=1 44

Ф3 = Ф7 = aJY JJY = Y1 aJ3 JJ3 ^ min-

j=1 j=1

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

В перечисленных принципах конфликтного взаимодействия заложены три фундаментальных понятия теории оптимального управления ММС (ТОУ ММС) [9]: стабильность; эффективность; стабильно-эффективный компромисс (СТЭК).

Стабильность ММС — обеспечение межобъектно устойчивых, сбалансированных по эффективности процессов функционирования и проектирования многообъектных структур на основе уравновешивания в условиях исходной структурной несогласованности, конфликта и неопределенности.

Эффективность ММС — достижение максимального целевого качества объектов, коалиции и ММС в целом на основе устойчивого и рационального коалицинирования.

Стабильно-эффективный компромисс — объединение стабильности и эффективности на множестве решений - от полного совпаде-

ния данных свойств в одной точке множества значении вектора показателей Ф (или множества управляющих параметров ф) до обеспечения возможной степени сближения в условиях информационно-тактических расширений компромиссов.

Прикладная задача рассматривается в условиях исходной структурной несогласованности, когда каналы стабилизации при наличии перекрестных связей необходимо бескоалиционно сбалансировать по эффективности (минимизации потерь Ф^, г = 1, 2, 3, в свертке по устойчивости, качеству, точности и быстродействию в каждом канале). Затем результат следует спроецировать в точку парето-границы множества значений отображения Ф(ф) (предельной эффективности ММС — трехканальной ССт), наиболее близкую к сбалансированной точке уравновешивания Фг = (Ф1, Ф2>, Ф3), где г — индекс равновесия; ф — общее множество управляющих параметров.

При бескоалиционном взаимодействии основной принцип оптимальности (стабильности) — равновесие по Нэшу.

Набор решений дг = (дг1,..., q) является равновесным по Нэшу относительно скалярных показателей ФС = ^^ , каждый из

которых представляет собой функцию потерь коалиции К^, если

€ г € Мк = (1, 2,...,Шк), ФС [с[ > ФС (qr),

где | qг) = {qr'1,..., qr'г-1, q\ qr'г+1,..., qr'mfc}. Здесь тк € Мк — множество коалиций. В прикладной задаче каждая коалиция состоит из одного канала ССт и тк = 3.

Для трехканальной ССт определение равновесия превращается в следующую систему неравенств при qг =

Ф2(qr1, q2, qr3) > Ф2(qr1, qr2, qr3); (6)

Ф3^г1^г2^3) > Ф3(qr1,qr2,qr3).

Другими словами, тот канал, в котором параметры отклонились от равновесия, теряет эффективность (потери Ф^ увеличиваются).

Простейшим необходимым условием (6) является система равенств при ^ е [2]:

дФ1 (q1, qr2, qr3)

%

дФ2 (qr1,q2,qr3)

дФз (qr1,qr2,q3)

= 0

= 0;

= 0,

где qr = (qr1, qr2, qr3) — равновесие по Нэшу.

Пусть множество индексов Mk = {1}, K = KMMC, Ф = = (Ф1,..., Фт) — вектор эффективности ММС, Ф» — скалярные показатели ММС. Вектор q0 £ Q оптимален по Парето, если из условия q £ Q, Ф(д) > Ф(д0), следует либо равенство Ф(д) = Ф(д0), либо система неравенств несовместима и хотя бы одно из неравенств противоположно по смыслу.

Простейшее необходимое условие оптимальности по Парето — следующая система равенств при q £ Q:

дФ1(д) =0; дФМ = 0; дФМ = 0

п 0; о 0; п 0.

dq dq dq

Исходная точка арбитражной схемы Нэша (АСН) — равновесное решение [2], а результат дает одну из точек области оптимальности по Парето, наиболее близкую к равновесной точке по значениям показателей. Таким образом, удается достичь объединения свойств стабильности (обеспечиваемой равновесным решением по Нэшу) и эффективности (обеспечиваемой оптимальностью решения по Парето), а также получить СТЭК.

Арбитражная схема Нэша при q £ Q [2] и нечетном числе показателей, существующая в трехканальной ССт, принимает вид

Фа = [Ф1 (q) - Ф1 (qr)] [Ф2 (q) - Ф2 (qr)] [Фз (q) - Фз (qr)] ^ min

q

при Ф^) < Ф»(qr), i = 1, 2, 3, Фа — арбитражная функция Нэша.

Запишем этапы метода равновесно-арбитражной балансировки.

Этап 1. Получение параметров системы стабилизации в каждом канале ССт (начальные приближения для равновесно-арбитражной оптимизации каналов ССт).

Этап 2. Нахождение балансировочного параметрического решения ССт на основе равновесия по Нэшу.

Этап 3. Парето-оптимизация параметрического решения ССт на основе арбитражной схемы Нэша.

На этапе 1 применяется методика расчета ССт в одном канале, сформулированная в работах [1, 8], которая включает в себя:

• метод стандартных коэффициентов, позволяет получить переходный процесс заданного вида (с требуемым перерегулированием) и значения параметром канала [8];

• оптимизацию внутреннего контура демпфирования ССт с ДУС (выбор параметра KWz, (KWy ,K>7)) по критерию минимизации резонансного пика частотной характеристики замкнутого контура демпфирования, что снижает колебательность канала [8];

• параметрический синтез во внешнем контуре ССт (выбор параметра Knz, (Kny, K7)) с использованием частотной характеристики для достижения требуемых запасов устойчивости и точности [8].

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

Разработка сетевого алгоритма поиска равновесия по Нэшу на этапе 2 равновесно-арбитражного метода. Для получения равновесия по Нэшу в трехканальной ССт по шести параметрам (два управляющих параметра в каждом канале) разработан алгоритм поиска на основе определения равновесия по Нэшу.

Для компактной иллюстрации разработанного алгоритма без ограничения общности рассмотрим работу алгоритма на примере двухко-алиционной игры К1,К2, где первая коалиция соответствует каналу тангажа, а вторая — каналу рыскания, с вектором параметров

Q =

K < K < K

-Ivwz min _ А vwz _ А vwz max

K < K < K

nz min _ A vnz _ A vnz max

K < K < K

min — 1 vWy — 1 vWy max

K < K < K

J1ns min — 1 *-ny — 1 *-ny max

В многофакторном анализе применяется алгоритм для трехканальной системы.

Предполагается, что на описанных выше интервалах сформирована сеть параметров с определенными шагами:

Q1 = \_Kwz min — K^z —Kuz max] = [Kwz min, Kwz min + shagi , . . . , KWz maj ; Q2 [Knz min —Knz —Knz max] [Knz min, Knz min + shag2 , ... , Knz max] j Q3= [KUy min —KWy —KLVy max] = \_K^y min, KWy min + shag3 , . . . , Kwy max] ; Q4 \_Kny min — Kny —Kny max] \_Kny mln, Kny min + shag4 , ... , Kny max] .

Пусть векторный целевой показатель сформирован в соответствии с формулами (1)-(4) в виде Ф = (Ф1, Ф2),

( (1,1) ••• (1, йш^з ^ш^)) \

Ф1=

у (dim(Qi) dim(Q2), 1) • • • (dim(Qi) dim(Q2), dim(Q3) dim(Q4))/

ф,=

/ (1,1) ••• (1, dim(Q3)dim(Q4)) \

v(dim(Qi) dim(Q2), 1) • • • (dim(Qi) dim(Q2), dim(Q3) dim^))/

где diш(Q2) — число различных комбинаций параметров кана-

ла тангажа; diш(Q3) diш(Q4) — число различных комбинаций канала рыскания.

Блок-схема для алгоритма поиска представлена на рис. 5.

Многофакторный анализ эффективности оптимизированной линеаризованной трехканальной ССт. На основе структурных схем контуров стабилизации бокового и продольного каналов, канала крена была построена программная модель в имитационной системе МЛТЬЛБ 81шиКпк, в которой все каналы связаны полным набором перекрестных связей. По экспериментальным данным было выделено три режима полета ЛА со следующими параметрами:

Рис. 5. Блок-схема алгоритма поиска равновесия по Нэшу (а) и простой отладочный пример (б)

Режим:

1.......................... М = 1, Н = 550 м, д = 6,6423е+04 Н/м2

2.......................... М = 1,5, Н = 8000 м, д = 5,6070е+04 Н/м2

3.......................... М = 1,06, Н = 15000 м, д = 9,4733е+03 Н/м2

Примечание. М — число Маха; Н — высота; д — скоростной напор.

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

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

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

Результаты оптимизации

Оптимизируемый параметр Критерий

Этап Кш, Knz Шу Kn ny Knx KY Ф^ Ф^ Ф7

Для режима 1

Равновесие по Нэшу 0,2940 1,6338 0,3427 1,9959 0,0528 0,3558 0,0665 0,0919 0,0308

СТЭК 0,2940 1,6338 0,3954 1,9959 0,0447 0,3558 0,0619 0,0747 0,0126

Для режима 2

Равновесие по Нэшу 0,4872 1,9922 0,4393 1,4867 0,0457 0,3849 0,0291 0,0393 0,0139

СТЭК 1,9922 1,4867 0,3849 1,4867 0,0457 0,3849 0,0213 0,0393 0,0139

Для режима 3

Равновесие по Нэшу 0,7410 1,33669 0,7345 1,2149 0,2970 2,5483 0,2100 0,0738 0,1651

СТЭК 0,7410 1,8576 0,7345 0,7072 0,2970 2,5483 0,2096 0,0731 0,1651

Пример области показателей для режима 2 приведен на рис. 6. На рисунке точка 2, полученная в результате программного применения

0,7

Рис. 6. Область показателей трехканальной ССт между каналами тангажа, рыскания и крена для режима 2:

1 — точка равновесия по Нэшу; 2 — точка СТЭК; 3 — точки области значений Ф

алгоритма арбитражной схемы Нэша [1], соответствует решению с предельной эффективностью и находится в области Парето.

Заключение. В результате проделанной работы получены оптимальные значения искомых параметров трехканальной ССт с перекрестными связями при трех режимах полета. С помощью этих значений можно аппроксимировать коэффициенты в функциональных зависимостях регуляторов для наилучшего функционирования ЛА в нерассмотренных режимах полета.

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

По сравнению с решением, равновесным по Нэшу (см. таблицу), удалось добиться улучшения показателей в точке СТЭК в форме равновесно-арбитражного решения, т.е. методы теории оптимизации

Рис. 7. Переходные процессы в исходной (а) и оптимизированной (б) системах

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

ЛИТЕРАТУРА

1. Воронов Е.М., Мелехина Ю.В., Веселовская О.А., Мусин Е.Р. Равновесно-арбитражная многокритериальная балансировка каналов в многосвязном регулировании и управлении // Вестник МГТУ им. Н.Э. Баумана. Сер. Приборостроение. 2007. № 4. С. 99-119.

2. Воронов Е.М. Методы оптимизации управления многообъектными многокритериальными системами на основе стабильно-эффективных компромиссов. М.: Изд-во МГТУ им. Н.Э. Баумана, 2001. 576 с.

3. Кузовков Н.Т. Системы стабилизации летательных аппаратов. М.: Высш. шк., 1976. 304 с.

4. Лебедев А.А., Карабанов В.А. Динамика систем управления беспилотными летательными аппаратами. М.: Машиностроение, 1965. 528 с.

5. Боднер В.А., Козлов И.М.Стабилизация летательных аппаратов и автопилоты. М: Оборонгиз, 1961. 508 с.

6. Бесекерский В.А., Попов Е.П. Теория систем автоматического регулирования. М.: Наука, 1975. 768 с.

7. КимД.П. Теория автоматического управления. Т. 1. М.: Физматлит, 2003. 288 с.

8. Веселов А.П., Евстифеев В.В., Карабанов В.А. Расчет и проектирование систем стабилизации летательных аппаратов. М.: МВТУ, 1979. Вып. 5. 22 с.

9. Коньков В.Г. Устойчивость управляемых технических систем / под ред. В.И. Сивцова. М.: Изд-во МГТУ им. Н.Э. Баумана, 1990. 78 с.

10. Аттетков А.В., Галкин С.В., Зарубин B.C. Методы оптимизации / под ред. B.C. Зарубина, А.П. Крищенко. М.: Изд-во МГТУ им. Н.Э. Баумана, 2003. 440 с.

REFERENCES

[1] Voronov E.M., Melekhina Yu.V., Veselovskaya O.A., Musin E.R. Equilibrium-arbitral multicriterial channel balancing in multivariate regulation and control. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Priborostr. [Herald of the Bauman Moscow State Tech. Univ., Instrum. Eng.], 2007, no. 4, pp. 99-119 (in Russ.).

[2] Voronov E.M. Metody optimizatsii upravleniya mnogoob'ektnymi mnogokriterial'nymi sistemami na osnove stabil'no-effektivnykh kompromissov [Control optimization methods of multicriterial multiobjective systems based on stably efficient compromises], Moscow, MGTU im. N.E. Baumana Publ., 2001. 576 p.

[3] Kuzovkov N.T. Sistemy stabilizatsii letatel'nykh apparatov [Stabilizing systems of airborne vehicles]. Мoscow, Vysshaya Shkola Publ., 1976. 304 p.

[4] Lebedev A.A., Karabanov V.A. Dinamika sistem upravleniya bespilotnymi letatel'nymi apparatami [Dynamics of control systems of airborne vehicles]. Moscow, Mashinostroenie Publ., 1965. 528 p.

[5] Bodner V.A., Kozlov M.S. Stabilizatsiya letatel'nykh apparatov i avtopiloty [Airborne vehicles stabilization and autopilot systems], Moscow, Oboroniz Publ., 1961. 508 p.

[6] Besekerskiy V.A., Popov E.P. Teoriya sistem avtomaticheskogo regulirovaniya [Theory of automatic regulation systems]. Moscow, Nauka Publ., 1975. 768 p.

[7] Kim D.P. Teoriya avtomaticheskogo upravleniya. T.1. Izd. 2 [Automatic control theory. In 2 volumes. 2nd ed.]. Moscow, Fizmatlit Publ., 2003. 288 p. (vol. 1).

[8] Veselov A.P., Evstifeev V.V., Karabanov V.A. Raschet i proektirovanie sistem stabilizatsii letatel'nykh apparatov [Calculation and design of stabilizing system of airborne vehicle]. Moscow, MVTU Publ., 1979. 22 p. (iss. 5).

[9] Kon'kov V.G., Sivtsov V.I. Ustoychivost' upravlyaemykh tekhnicheskikh system [Stability of controlled engineering system]. Moscow, MGTU im. N.E. Baumana Publ., 1990. 78 p.

[10] Attetkov A.V., Galkin S.V., Zarubin B.C., Krishchenko A.P. Metody optimizatsii [Optimization methods]. Moscow, MGTU im. N.E. Baumana, 2003. 440 p.

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

Алексей Сергеевич Аксенов — инженер-конструктор ОАО "Корпорация "Тактическое ракетное вооружение". Автор ряда работ в области теории управления и ее приложений в системах управления летательных аппаратов.

ОАО "Корпорация "Тактическое ракетное вооружение", Российская Федерация, 141080, Московская область, Королев, ул. Ильича, д. 7.

A.S. Aksenov — design engineer of OAO "Corporation "Tactical Missile Armament". Author of a number or publications in the field of control theory and its applications in the aircraft control systems.

OAO "Corporation "Tactical Missile Armament", ul. Il'icha 7, Korolev, Moscow region, 141080 Russian Federation.

Евгений Михайлович Воронов — д-р техн. наук, профессор кафедры "Системы автоматического управления" МГТУ им. Н.Э. Баумана. Автор более 200 научных работ в области теории управления, систем управления и ее приложений. МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.

E.M. Voronov — Dr. Sci. (Eng.), professor of "Automatic Control Systems" department of the Bauman Moscow State Technical University. Author of more than 200 publications in the field of control theory and its applications, control systems. Bauman Moscow State Technical University, Vtoraya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.

Кирилл Константинович Любавский — студент кафедры "Системы автоматического управления" МГТУ им. Н.Э. Баумана. Соавтор докладов международных конференций по теории управления и ее приложений.

МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.

K.K. Lyubavskiy — student of "Automatic Control Systems" department of the Bauman Moscow State Technical University. Co-author of reports at international conferences on control theory and its applications.

Bauman Moscow State Technical University, Vtoraya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.

Станислав Игоревич Сычев — начальник отдела опытно-конструкторского бюро ОАО "Корпорация "Тактическое ракетное вооружение". Автор ряда работ в области теории управления и ее приложений в системах управления летательных аппаратов. ОАО "Корпорация "Тактическое ракетное вооружение", Российская Федерация, 141080, Московская область, Королев, ул. Ильича, д. 7.

S.I. Sychev — head of department the Research and Development Bureau of OAO "Corporation "Tactical Missile Armament". Author of the number of publications in the field of control theory and its applications in aircraft control systems. OAO "Corporation "Tactical Missile Armament", ul. Il'icha 7, Korolev, Moscow region, 141080 Russian Federation.