Рыжиков И.С. (ryzhikov-88@yandex.ru) - Сибирский государственный аэрокосмический университет им. академика М.Ф. Решетнева (аспирант), Красноярск, Россия, Семенкин Е.С. (styugin@rambler.ru) - Сибирский государственный аэрокосмический университет им. академика М.Ф. Решетнева, г. Красноярск, Россия | |
Ключевые слова: реле., терминальное управление, управление, динамическая система, эволюционные стратегии |
|
Keywords: relay, terminal control, control management, dynamic system modeling, evolution strategies |
|
|
Задача нахождения программного управления для систем, исполнительным механизмом в которых является идеальное или многоуровневое реле, изучена недостаточно. На сегодняшний день не существует метода нахождения управления в виде кусочно-постоянной функции, что обусловлено не только ее недифференцируемостью, но и тем, что нет единого метода нахождения управления для нелинейных систем. В [1] излагается метод коррекции точек переключения, при этом рассматривается линейная динамическая система, известно число точек переключений и начальное приближение предполагается достаточно близким к решению задачи. В [2] исследуется синтез реле для линейных динамических систем. Задача терминального управления является актуальной и описывается в классических учебниках по оптимальному управлению и в современных работах, например в [3, 4]. В данной работе рассматриваются метод решения экстремальной задачи, к которой была сведена терминальная задача управления, и система нахождения программного управления. Постановка задачи Для системы, описываемой дифференциальным уравнением в общем виде
где Также, если гамильтониан [5] линеен по функции управления, то есть H(x, p, u)=j1(x, p)+j2(x, p)u, а функция управления имеет ограничения по модулю, êu(t)ê Если реле является многоуровневым, функция управления представима в аналогичном виде:
где Пусть Вопрос о выборе числа точек переключения остается открытым, данная постановка задачи и предложенный далее метод решения предполагают неявную настройку количества точек переключения. Прежде всего, если некоторая точка переключения больше конечного времени правления, r>T, то она не влияет на процесс управления. Неявная настройка числа точек переключения обусловлена применением численной схемы интегрирования системы (1). Пусть h – шаг интегрирования, тогда при выполнении условия êri+1–ri–1ê В общем случае решение терминальной задачи для фиксированного времени управления совпадает с решением экстремальной задачи:
где
Существуют также ограничения в виде неравенств
которые обеспечивают выполнение условия нахождения каждой точки переключения внутри интервала [0, T]. В случае использования множества Для решения задачи условной оптимизации введем функцию штрафа
Такие постановки задачи рассматривались в [6]. Добавим еще несколько новых вариантов. Для ограничений ri для любого из критериев (4), (5), (7), (8). В этом случае получаем новые целевые функции, например для (5):
Таким образом, задача терминального управления сводится к задаче оптимизации в пространстве вещественных чисел или вещественных и натуральных чисел с целевыми функциями (3), (4), (6)–(8). Указанные функции многоэкстремальны и не являются аналитически заданными. Учитывая сложность решаемых задач, размерность и неоднородность входных переменных, в качестве алгоритма решения задачи целесообразно выбрать метод эволюционных стратегий [7]. Рассмотрим некоторые особенности работы алгоритма и внесенные модификации. Алгоритм решения экстремальной задачи Пусть каждый индивид задан кортежем
где Особенности задач (3), (4), (6)–(8) задают общий вид набора объективных параметров Для решаемых задач эффективность алгоритма возрастает при использовании измененной операции мутации. Пусть
Пусть Следующая модификация алгоритма обеспечивает выполнение неравенств (5). Каждый объективный параметр задает точку переключения следующим образом: Алгоритм случайного поиска дополняется случайным покоординатным спуском. Из N1 случайно выбранных индивидов случайно выберем N2 объективных хромосом и осуществим N3 шагов локального спуска величиной hl для определения лучшего решения. Гибридизация алгоритма позволяет значительно повысить точность найденных решений и эффективность в целом. В рассмотренных ниже экспериментах эмпирически выбранные настройки алгоритма оптимизации были следующими: размер популяции – 50 индивидов, количество популяций – 50, турнирная селекция с объемом 20 %, дискретное скрещивание, модифицированная мутация с вероятностью для каждого гена 1/úspê. Настройки локального спуска: N1=2×úopê, N2=úopê и N3=5 с шагом hl=0,05. Предложенный алгоритм был проверен на десяти специальных тестовых задачах, что позволило определить наиболее эффективные настройки. Применяя предложенный алгоритм решения оптимизационной задачи для выбранного представления структуры функции управления, может быть решена терминальная задача программного управления. Данный подход реализован в среде OCTAVE/MATLAB и может использоваться для решения любой из представленных задач управления. Так же реализована система решения задачи управления с выбранной пользователем постановкой, дружественным интерфейсом. Система разработана на языке программирования С++, результат решения задачи представляется в виде графиков выхода системы и функции управления. Аналогичные результаты сохраняются в отдельные файлы решения задачи терминального управления, где сохраняются численные решения. Параметры численной схемы интегрирования выбираются пользователем, при этом в программе по умолчанию выбран шаг интегрирования. Примеры решения задач управления Рассмотрим задачу перевода спутника с одной геостационарной орбиты на другую. Координаты спутника описываются нелинейным дифференциальным уравнением вида
Необходимо перевести объект из начального положения x(0)=(1, 0, 1, 0,785) в конечное x(T)= =(1, 0, 1, T) за конечное время T.
Рассмотрим задачу нахождения программного управления для случая, когда амплитуда реле является искомой, то есть задачу синтеза регулирующего устройства. Решение таких задач может быть полезно на этапе проектирования систем управления. Выберем следующие условия для терминальной задачи идеального релейного управления: –x¢¢+x¢×x2+x×u(t)=0, x(0)= По условиям задачи необходимо определить и точки переключений, и амплитуду реле. Воспользуемся критерием (8) при k=10, значение весового коэффициента было принято равным a=10. Усредненное решение задачи управления по 20 запускам алгоритма
Подводя итоги, отметим, что в работе был рассмотрен специальный подход к решению задачи перевода некоторой динамической системы из заданного состояния в желаемое притом, что исполнительный механизм может работать только в некоторых заранее известных или проектируемых режимах. Подход состоит в переводе задачи управления в задачу поиска оптимума, условного или нет, и в решении полученной оптимизационной задачи. Универсальность предложенного метода позволяет решать задачи в различных постановках, то есть для каждого отдельного семейства задач не требуется новый подход. Предложенный алгоритм, основанный на методе эволюционных стратегий, при нахождении экстремума тестировался на определенном наборе задач управления для различных систем, после чего были определены его наиболее эффективные настройки. Алгоритм разработан таким образом, что поиск происходит одновременно на пространствах с различным типом данных: вещественными и целыми числами, без перевода всех типов данных в бинарную строчку. Таким образом, выполняются все необходимые условия для решаемой задачи. На основании проведенных исследований был реализован и тестирован программный продукт. Представленные в работе примеры решения задач позволяют сделать вывод о том, что с помощью предложенного метода можно не только находить искомое решение задачи терминального управления для управления в виде реле, но и корректировать точки переключения. Дальнейшие исследования направлены на изучение влияния количества точек переключений на эффективность решаемых задач. Литература 1. Aida-zade K.R., Rahimov A.B., Journ. of Autom. and Inform. Sc., 2010, Vol. 42, Iss. 7, pp. 54–64. 2. Кучеров Д.П., Василенко А.В., Иванов Б.П. Алгоритм адаптивного терминального управления динамической системой с элементами дифференцирования // Автоматика. Автоматизация. Электротехнические комплексы и системы. 2009. № 23. C. 166–171. 3. Трушкова Е.А. Синтез управления в окрестности приближенного решения задачи с частично закрепленным правым концом // Программные системы: теория и приложения. 2011, № 2 (6). C. 31–35. 4. Зубер И.Б. Терминальное управление по выходу для нелинейных нестационарных систем // Дифференциальные уравнения и процессы управления. 2004. № 2. C. 36–42. 5. Алексеев В.М., Тихомиров В.М., Фомин С.В. Оптимальное управление. М.: Физматлит, 2005. 384 с. 6. Ryzhikov I.S., Semenkin E.S., Proc. Intern. Conf. «ICINCO», 2012, Vol. 1, pp. 333–337. 7. Schwefel H.-P. Evolution and Optimum Seeking: New York, Wiley & Sons., 1995, 456 p. |
http://swsys.ru/index.php?id=3408&lang=%E2%8C%A9%3Den&like=1&page=article |
|