Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,118

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

Корнюшин Ю.П. 1
1 Калужский филиал ФГБОУ ВПО «Московский государственный технический университет имени Н.Э. Баумана (национальный исследовательский университет)»
В статье рассматриваются вопросы синтеза оптимальных по быстродействию программных управлений. Рассматриваются объекты, относящиеся к классу нелинейных аффинных систем со скалярным управлением. В качестве математического аппарата используется линеаризация Ньютона – Канторовича, аппарат матричные операторы и абстрактная конечномерная L-проблема моментов М.Г. Крейна и А.А. Нудельмана. Алгоритм состоит из следующих основных этапов. 1. Линеаризация нелинейной математической модели объекта по схеме Ньютона – Канторовича. 2. Преобразование дифференциальных уравнений к форме с матричными операторами. 3. Синтез управления, согласно заданному критерию с использованием L-проблемы моментов. 4. Организация итерационного процесса, согласно схеме Ньютона – Канторовича. На 3 этапе используются методы функционального анализа, которые сводят задачу к классической вариационной задаче – максимизации системы моментных соотношений при ограничении на управление, которое рассматривается как равенство. Приводится пример иллюстрирующий метод синтеза.
синтез
оптимальный
быстродействие
нелинейный
матричный
оператор
проблема моментов
1. Александров В.М. Оптимальное по быстродействию управление в реальном времени линейными системами с возмущениями // Вестник НГУ. Серия: Математика, механика, информатика. – 2008. – Т. 8, вып. 3. – C. 3–25; URL: http://www.nsu.ru/xmlui/bitstream/handle/nsu/2702/01.pdf (дата обращения: 26.11.2016).
2. Ибрагимов Д. Н. Оптимальное по быстродействию управление движением аэростата // Труды МАИ. – 2015. – Выпуск № 83. – C. 1–31; URL: http://www.mai.ru/upload/iblock/bf3/ibragimov_rus.pdf (дата обращения: 26.11.2016).
3. Карагодин В.В., Горин В.А., Вишняков Е.П. Математические основы построения оптимального по быстродействию управления возбуждением синхронного генератора // Вопросы электромеханики. – 2013. – Т. 134. – C. 3–8.
4. Корнюшин Ю.П., Мельников Д.В. Концепция синтеза регуляторов следящих систем для класса нелинейных объектов управления на основе метода матричных операторов. // Автоматика. – 2008. – Т. 15. – С. 268.
5. Корнюшин Ю.П., Егупов Н.Д., Корнюшин П.Ю. Идентификация нелинейных объектов и систем управления с использованием аппарата матричных операторов // Электронное научно-техническое издание. Инженерный журнал: наука и инновации. – 2014. – № 6 (30). – С. 1–14; URL: http://engjournal.ru/articles/1315/1315.pdf (дата обращения: 26.11.2016).
6. Корнюшин Ю.П., Мельников Д.В., Егупов Н.Д., Корнюшин П.Ю. Синтез дополнительного регулятора для стабилизации угловой скорости ротора паровой турбины. // Вестник Московского государственного технического университета им. Н.Э. Баумана. Серия: Естественные науки. – 2015. – № 5 (62). – С. 100–112.
7. Корнюшин Ю.П., Егупов Н.Д., Акименко Д.А., Корнюшин П.Ю. Синтез нелинейных робастных регуляторов методом матричных операторов // Известия Тульского государственного университета. Технические науки. – 2011. – № 5–1. – С. 99–108.
8. Крейн М.Г., Нудельман А.А. Проблема моментов Маркова и экстремальные задачи. – М.: Наука, 1973. – 551 с.
9. Крендлер Е. Развитие теории оптимальных по быстродействию процессов управления / Кибернетический сборник. – М.: Издательство «Мир», 1975. – № 3. – С. 111–148.
10. Ловчаков В.И. Функции переключения оптимального по быстродействию регулятора для четырехкратного интегратора «Мехатроника, автоматизация, управление». – 2014. – № 9(162). – С. 3–6; URL: http://novtex.ru/mech/mech2014/annot09.html (дата обращения: 26.11.2016).
11. Мельников Д.В., Чжо Т.А., Окар М.М., Мин Ч.Т. Синтез систем регулирования первичных двигателей синхронных генераторов. // Фундаментальные исследования. – 2016. – № 10–3. – С. 509–515.
12. Моисеев Н.Н. Математические задачи системного анализа. – М.: Наука, 1981. – 488 с.
13. Понтрягин Л.С., Болтянский В.Г., Гамкрелидзе Р.В., Мищенко Е.Ф. Математическая теория оптимальных процессов – М.: Физматгиз, 1961. – 384 c.
14. Сухинин Б.В., Сурков В.В., Соловьев А.Э. Синтез оптимальных по быстродействию систем на основании использования теоремы об интервалах управления // Вести высших учебных заведений черноземья. – 2010. – № 2(20). – C. 57–63; URL: http://www.stu.lipetsk.ru/files/materials/3945/2010_2_012.pdf (дата обращения: 26.11.2016).
15. Ding M.-F., Liu Y., Gear J.A. A Modified Centered Climbing Algorithm for Linear Programming // Appl. Math. – 2012. – V. 3. – Р. 1423–1429; URL: http://file.scirp.org/pdf/AM20123000030_42233952.pdf (дата обращения: 26.11.2016).

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

Основной проблемой является необходимость решения краевых задач и наличие ограничений на управление. В случае нелинейных систем – разнообразие видов нелинейностей. Универсального подхода нет. Эти особенности приводят к необходимости в каждом конкретном случае разрабатывать свой метод решения задачи. Разработано множество методов, имеющих свои достоинства и недостатки [1–3, 10, 14, 15]. Список крайне неполон. Предлагается один из возможных подходов к решению задачи. Он основан на аппарате матричных операторов.

Аппарат матричных операторов нашел свое широкое применение для решения самых различных задач при исследовании и проектировании систем управления [5–7, 11, 12]. Его применение совместно с другими методами позволяет существенно расширить круг решаемых задач. Одной из таких является задача синтеза оптимальных по быстродействию программных управлений для нелинейных объектов. Желательно иметь решение задачи в замкнутом виде, однако это не всегда возможно. Если найдено оптимальное программное управление, то задача синтеза регулятора для нелинейного объекта упрощается. Она сводится, например, к задаче стабилизации движения объекта относительно программной траектории и программного управления, так называемая задача двухэтапной оптимизации [12].

В работе предлагается использовать аппарат матричных операторов совместно с абстрактной конечномерной L-проблемой моментов М.Г. Крейна и А.А. Нудельмана [8]. В работе [9] Е. Крейндлер показал, как можно использовать положения L – проблемы моментов для синтеза оптимального по быстродействию управления для класса линейных объектов.

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

Постановка задачи. Задан объект, описываемый уравнением вида

kor01.wmf (1)

где kor02.wmf, kor03.wmf.

Полагаем, что объект является вполне управляемым. Необходимо синтезировать скалярное управление u(t), переводящее объект из заданного начального состояния X(t0) в конечное состояние X(T) и доставляющее минимум следующему функционалу

kor04.wmf (2)

Алгоритм синтеза

Алгоритм состоит из следующих основных этапов.

1. Линеаризация нелинейной математической модели объекта по схеме Ньютона – Канторовича.

2. Преобразование дифференциальных уравнений к форме с матричными операторами.

3. Синтез управления, согласно заданному критерию с использованием L-проблемы моментов.

4. Организация итерационного процесса, согласно схеме Ньютона – Канторовича.

Линеаризация нелинейной математической модели

Данная процедура достаточно подробно рассмотрена [6, 7, 11] и в других работах автора, посвященных задачам синтеза. Отметим лишь, что в процессе линеаризации осуществляется переход от нелинейных дифференциальных уравнений (1) к последовательности линеаризованных уравнений вида

kor05.wmf (3)

kor06.wmf,

где

kor07.wmf

kor08.wmf (4)

Матрица kor09.wmf определена следующим образом:

kor10.wmf (5)

Линеаризованные дифференциальные уравнения (3) имеют переменные коэффициенты. Их структура соответствует той, с которой обычно выполняется синтез управления для неавтономных объектов.

Преобразование дифференциальных уравнений к необходимой форме

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

Для этого перейдем от дифференциального уравнения (3) к интегральному уравнению Фредгольма 2 – рода (индекс k для краткости записи опустим):

kor11.wmf

kor12.wmf, (6)

где матричные ядра kor13.wmf, kor14.wmf, kor15.wmf имеют вид

kor16.wmf,

kor17.wmf,

kor18.wmf.

Представим сигналы X(t), X(t0), u(t) и ядра уравнения (6) kor19.wmf, kor20.wmf, kor21.wmf в виде конечномерного разложения по заданному ортогональному базису kor22.wmf:

kor23.wmf, kor24.wmf,

kor25.wmf,

kor26.wmf,

kor27.wmf,

kor28.wmf, (7)

где kor29.wmf, kor30.wmf, kor31.wmf, kor32.wmf, kor33.wmf, kor34.wmf – спектральные характеристики соответствующих функций, kor35.wmf, kor36.wmf, kor37.wmf, kor38.wmf – клеточные матрицы, построенные на элементах ортогонального базиса.

Теперь можно перейти от интегральной формы описания объекта к спектральной форме

kor39.wmf

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

kor40.wmf, (8)

где

kor41.wmf,

kor42.wmf, kor43.wmf,

kor44.wmf, kor45.wmf.

Выражение (8) можно записать также в виде

kor46.wmf (9)

Здесь

kor47.wmf.

Полагая t = T, получаем систему моментных соотношений

kor48.wmf (10)

В зависимости (10) W(t, τ) – системой моментных функций, kor49.wmf – система моментов.

Синтез управления

Управление оптимальное u(t) по критерию (2) всегда ограничено. Это может быть ограничение по величине kor50.wmf, по энергии kor51.wmf или по площади kor52.wmf. Все они объединяются одним выражением

kor53.wmf. (11)

При p = ∞ и q = 1 имеет место ограничение управления по величине.

Левая часть неравенства (11) есть норма u(t) в пространстве Lp. Поэтому при синтезе управления можно воспользоваться методами функционального анализа. В [9] отмечается, что «управление будет оптимальным только в случае, если оно удовлетворяет неравенству (11) со знаком равенства». Кроме того для оптимальности управления необходимо и достаточно, чтобы оно максимизировало следующее выражение

kor54a.wmf, (12)

при некотором ненулевом векторе ∧.

Как видим, имеет место классическая вариационная задача – максимизация выражения (13) при ограничении (12), которое рассматривается, как равенство. Решение задачи известно – [11, 18]. Оптимальное управление определяется следующим выражением:

kor55f.wmf (13)

При условии

kor56.wmf (14)

Как видим из (13), оптимальное управление является функцией времени и зависит от значений элементов ненулевого вектора – ∧, минимального времени – T и величины ограничения – Cp: kor57.wmf.

Для значений p = ∞ и q = 1 имеет место известный результат [8, 13]

kor58.wmf (15)

С учетом вида управления (15), условие (14) становится следующим:

kor59.wmf (16)

Если длину вектора ∧ выбрать произвольно, например, так, чтобы

kor61.wmf, (17)

то условие (17) упрощается:

kor62.wmf (18)

Выражение (18) есть решение поставленной задачи на максимальное быстродействие.

Наименьшее T, которое удовлетворяет данному уравнению для заданного значения C при условии (17) есть искомое минимальное время для текущего шага итерационного процесса.

Организация итерационного процесса

Линеаризация по схеме Ньютона – Канторовича предполагает организацию итерационного процесса по уточнению решения уравнения (3) до тех пор, пока решение линеаризованного уравнения не будет близко к решению исходного нелинейного уравнения. Как правило, начальные приближения выбираются исходя из некоторых практических соображений, например, при решении данной задачи в качестве них можно использовать начальные состояния, из которых переводится объект. На текущем шаге итерационного процесса находится система моментных функций – W(t, τ) и система моментов – kor63.wmf (смотри систему моментных соотношений (10)), находится оптимальное время управления T, удовлетворяющее уравнению (18) и искомое управление u(t). Найденному управлению соответствует своя фазовая траектория, являющаяся основой для следующего шага итерационного процесса. Итерационный процесс продолжается до тех пор, пока с заданной точностью не будет выполнено ограничение на управление и не будет меняться время управления T.

Пример. Объект управления описывается уравнением

kor64.wmf,

где kor65.wmf.

Требуется синтезировать управление, переводящее объект из начального состояния kor66.wmf в конечное состояние kor67.wmf за минимальное время при ограничении на управление kor68.wmf.

kornyh1.wmf

Рис. 1. Управление u(t)

kornyh2.wmf

Рис. 2. Изменение координаты x1(t)

kornyh3.wmf

Рис. 3. Изменение координаты x2(t)

Начальные приближения kor69.wmf, kor70.wmf. Найденное время управления T = 0,7216256 с. Полученное конечное состояние: kor71.wmf. Ограничение на управление: kor72.wmf. Данная точность получена за 4 итерации.

На рис. 1 показан график изменения управления; на рис. 2–3 графики изменения фазовых координат. На рис. 2–3 точками показаны точные конечные значения фазовых координат.

Заключение

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

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований и Правительства Калужской области (грант № 16-41-400701).


Библиографическая ссылка

Корнюшин Ю.П. ПРИМЕНЕНИЕ АППАРАТА МАТРИЧНЫХ ОПЕРАТОРОВ ДЛЯ СИНТЕЗА ОПТИМАЛЬНЫХ ПО БЫСТРОДЕЙСТВИЮ СИСТЕМ УПРАВЛЕНИЯ // Фундаментальные исследования. – 2016. – № 11-5. – С. 925-930;
URL: http://fundamental-research.ru/ru/article/view?id=41278 (дата обращения: 21.11.2018).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.252