САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
ФАКУЛЬТЕТ ПРИКЛАДНОЙ МАТЕМАТИКИ – ПРОЦЕССОВ УПРАВЛЕНИЯ
Севостьянов Руслан Андреевич
Выпускная квалификационная работа аспиранта
Оптимизация алгоритмов цифрового управления
подвижными объектами
Направление 09.06.01 "Информатика и вычислительная техника"
Программа МК.3021.2014 "Системный анализ, информатика и управление"
Научный руководитель,
доктор физ.-мат. наук,
профессор
Веремей Е. И.
Рецензент,
кандидат физ.-мат. наук
Арефина А. И.
Санкт-Петербург
3
Содержание
Введение
3
Постановка задачи
7
Обзор литературы
10
Глава 1. Многоцелевые регуляторы 13
1.1. Общая структура закона управления 13
1.2. Поиск базового регулятора
15
1.3. Поиск коэффициентов асимптотического наблюдателя 21
1.4. Синтез корректирующего устройства
Глава 2. Учет запаздывания
23
28
2.1. Применения прогноза для компенсации запаздывания 28
2.2. Учет запаздывания для динамических регуляторов
Глава 3. Объект управления
31
37
3.1. Суда на воздушной подушке
37
3.2. Уравнения нелинейной динамики
39
3.3. Линейная модель объекта управления
42
Глава 4. Практическая реализация многоцелевого регулятора
4.1. Коэффициенты модели
43
43
4.2. Режим отработки командного сигнала
46
4.3. Компенсация ступенчатых возмущений
47
4.4. Коррекция морского волнения 48
Глава 5. Практическая реализация компенсирующего регулятора 53
5.1. Режим отработки командного сигнала
53
5.2. Компенсация ступенчатых возмущений
55
5.3. Коррекция морского волнения 57
Выводы
.. 60
Заключение
62
Список литературы
63
Перечень собственных публикаций по теме
65
4
Введение
В настоящее время с каждым днем возрастают потребности в
использовании
транспорта
различного
рода.
Растет
поток
клиентов
пассажирских перевозок, увеличивается объем транспортировки грузов.
Помимо требований к количеству единиц транспорта усиливаются требования
к качеству передвижения и транспортировки, будь то безопасность и комфорт
пассажиров, или аккуратность перевозки хрупких предметов.
Отдельный интерес представляют различные виды морских подвижных
объектов (МПО). В этой области существует очень большое разнообразие
моделей, принципиально различающихся по своим характеристикам, и
имеющих свои достоинства и недостатки. В качестве частного примера одного
из перспективных классов современных МПО, который будет использован
далее в работе для демонстрации применения используемых методов, является
класс судов на воздушной подушке (СВП). В принципе, одной из главнейших
задач при проектировании любого морского транспорта является уменьшение
сопротивления воды за счет уменьшения площади соприкосновения корпуса
судна с водой. СВП решают эту задачу своеобразным путем. Между корпусом
судна и поверхностью воды поддерживается область повышенного давления
воздуха (так называемая «воздушная подушка», отсюда и название класса
судов), так что прямой контакт корпуса с водой отсутствует. Следует отметить,
что существует два вида СВП: скеговые и амфибийные. В скеговых СВП
область воздушной подушки ограничивается жесткими боковыми стенками,
которые погружены в воду. Этот вид СВП по динамике больше похож на
классические морские суда, но может развивать большие скорости с меньшими
энергозатратами.
Область воздушной подушки амфибийных СВП (АСВП) ограничена
гибким ограждением, не имеющим прямого контакта с поверхностью воды.
Эта
особенность
обеспечивает
возможность
движения,
в
принципе,
практически по любой поверхности. За счет этого очень распространена
5
эксплуатация АСВП в труднодоступных местах, например, на болотах или
арктических льдах. Кроме этого АСВП часто используют в качестве транспорта
для переброски военного десанта с суши на морские суда.
Динамика АСВП в целом отличается от динамики классических морских
судов.
Отсутствие
прямого
контакта
корпуса
с
поверхностью
воды
минимизирует гидродинамическое сопротивление, за счет чего существенно
повышается инерционность бокового движения, оказывая влияние на динамику
дрейфа. Однако гидродинамическое сопротивление не исчезает полностью.
Под корпусом АСВП образуется провал водной поверхности. Рябь воды влияет
на распределение давления внутри воздушной подушки, за счет чего
образуется сила сопротивления. На малом ходу сопротивление возрастает с
увеличением скорости. Однако после преодоления определенного значения
скорости движения (т.н. «горбовая скорость») сопротивление резко падает. Это
определяет вторую важнейшую особенность динамики АСВП. Кроме этого,
АСВП в гораздо большей степени подвержены влиянию порывов ветра, чем
классические морские суда. Указанные особенности обосновывают наличие
специальных требований и ограничений к качеству управления АСВП в
различных режимах движения.
С
развитием
техники,
информационных
технологий
и
методов
прикладной математики появляется все больше возможностей облегчения
работы людей, пилотирующих различные транспортные средства, вплоть до
полной замены человека электронной системой управления. Помимо
улучшения качества функционирования, это позволяет увеличить гарантии
безопасности, поскольку довольно часто к транспортным авариям приводит т.н.
«человеческий
фактор».
Системы
автоматического
управления
могут
реагировать быстрее, просчитывать большое количество вариантов развития
событий, применять сложные алгоритмы и методы для достижения
поставленных целей. В целом автоматические устройства совершают гораздо
меньше ошибок, чем люди-операторы. В качестве примеров подобных систем,
существующих в настоящее время можно привести автопилоты самолетов и
6
кораблей, системы ADAS (advanced driver assistant systems) – помощь
водителям автомобилей.
В настоящее время подавляющее большинство систем автоматического
управления
реализуется
на
цифровых
устройствах:
компьютерах
и
микроконтроллерах. Использование цифровых устройств вносит целый ряд
особенностей в формирование управления. Во-первых, очевидно, управление
должно иметь дискретный характер, т.е. вычисляться и подаваться на объект
управления в отдельные точечные моменты времени, между которыми
управляющий сигнал остается постоянным. Во-вторых, работа с цифровыми
устройствами неразрывно связана с процессами хранения, обработки и
передачи информации. Это, в свою очередь, приводит к проблеме наличия
запаздываний в системе. Вообще, в любой реальной физической системе в той
или иной мере присутствуют запаздывания. Например, рулевой привод не
может мгновенно изменить положение руля. Измерения с датчиков тоже
происходят за некоторое ненулевое время. Эти запаздывания могут быть и
совсем незаметными, а могут быть достаточно большими, оказывая
существенное влияние на динамику объекта управления. Чаще всего наличие
запаздываний
приводит
к
ухудшению
динамики,
вплоть
до
потери
устойчивости, однако известны случаи, когда наличие запаздываний оказывает
положительное влияние на качество функционирования системы управления.
Цифровая природа бортовых вычислительных устройств добавляет еще
один тип запаздывания: вычислительное. Расчет управляющих воздействий
также занимает некоторое ненулевое время, которое может оказаться довольно
существенным.
Далее в данной работе будут рассмотрены две задачи. Первая относится к
автоматическому управлению произвольным подвижным объектом с учетом
ряда требований к качеству управления и ограничений за счет использования
специальной структуры закона управления. Вторая задача относится к учету
запаздывания при построении закона управления. При этом важным
требованием является сохранение конечного спектра собственных чисел и
7
передаточной матрицы исходной замкнутой системы без запаздывания, т.е.
сохранение динамических характеристик. Эта задача может быть решена на
основе т.н. компенсационного подхода, суть которого заключается в
определенной трансформации исходного регулятора для учета запаздывания в
контуре управления.
В
качестве
частного
иллюстративного
примера
использования
предложенных подходов к оптимизации законов управления за счет учета
запаздывания будет рассмотрена задача стабилизации курсового угла АСВП
при движении над морской поверхностью.
8
Постановка задачи
В рамках данной работы рассматриваются две задачи. Первая
относится к управлению движением произвольного подвижного объекта с
учетом ряда требований и характеристик при отсутствии запаздывания
командного сигнала. Вторая задача касается учета запаздывания при
проектировании системы управления, при этом важным моментом является
сохранение динамических свойств системы без запаздывания.
Рассмотрим нелинейные уравнения динамики подвижного объекта
управления в отклонениях от некоторого контролируемого движения:
G ( x, δ, d) ,
x
(1)
где x E – вектор состояния объекта, δ E m – вектор отклонений органов
n
управления (например, рулей) от нейтрального положения, d E nd – вектор
возмущающих воздействий.
Синтез управляющего сигнала будем проводить по выходу системы (1),
линеаризованной в окрестности положения равновесия при постоянной
продольной составляющей скорости хода и при отсутствии управляющего и
возмущающего воздействий. Результатом линеаризации является система
обыкновенных линейных дифференциальных уравнений с постоянными
коэффициентами (LTI-система) вида
(2)
характеризует отклонения вектора состояния объекта
x Ax Bδ Hd ,
где вектор
x
управления от положения равновесия, A , B , H – матрицы с постоянными
компонентами.
Будем считать, что полный вектор состояния недоступен для
измерения. Добавим к системе (2) уравнения выхода
(3)
– вектор измеряемых координат, С – матрица с постоянными
y Cx,
где y E
k
компонентами.
9
Кроме этого, в дополнение к приведенным выше уравнениям следует
рассматривать нелинейные уравнения приводов исполнительных органов,
используемых для управления движением
δ u σ , σ f u (u) , δ f (δ u ) ,
где u E
δu E m
m
(4)
– управляющий сигнал, σ E m – скорость отклонения рулей,
– вспомогательный вектор. Здесь векторные функции f u , f
представляют существенные нелинейности, например, типа «срезки»
f u (u) sat (u, u 0 ) и f (δ u ) sat (δ, δ 0 ) ,
где параметры u 0 и δ 0 задают величины соответствующих ограничений.
Для системы (2),(3),(4) требуется синтезировать закон управления,
удовлетворяющий следующим требованиям:
1) необходимо обеспечить асимптотическую устойчивость системы в
любом режиме движения;
2) качество отработки командного сигнала по курсовому углу должно
удовлетворять
определенным
требованиям
(чаще
всего
это
требования по времени перехода и величине перерегулирования),
формальное определение требований будет дано ниже;
3) при воздействии постоянного ступенчатого возмущения (например,
влияние ветра) отклонение курсового угла не должно превышать
заданное значение;
4) реакция системы на периодическое возмущение (например, морское
волнение) так же должна удовлетворять определенным требованиям,
при этом возможна реализация одного из двух режимов. В первом
режиме основной задачей является минимизация реакции выхода
системы (режим «точный»). Во втором режиме в центре внимания
находится минимизация интенсивности перекладки рулей (режим
«экономичный»). В рамках данной работы будет рассмотрен только
последний режим.
10
Далее добавим в систему (2),(3) постоянное запаздывание по входам:
(t ) Ax(t ) Bδ(t h) Hd (t h),
x
y Cx,
(5)
где в дополнение введенным выше обозначениям указано заданное
вещественное число h 0 , определяющее запаздывание по управлению и
возмущению. Здесь будем считать, что вектор измерений
состоянием
x , при этом векторы u ,
y совпадает с
d и y могут быть непосредственно
измерены.
Необходимо обеспечить асимптотическую устойчивость системы (5) с
учетом наличия запаздывания. При этом в центре внимания находится задача
сохранения динамических свойств системы без запаздывания (2),(3),
замкнутой регулятором, синтезированным в процессе решения предыдущей
задачи.
11
Обзор литературы
В современной литературе довольно глубоко описаны методы анализа и
синтеза стабилизирующих обратных связей с учетом ряда требований и
ограничений. Однако, как правило, эти методы предназначены для применения
к одному определенному режиму движения объекта управления и учитывают
только требования по отдельным динамическим характеристикам. В качестве
примера можно выделить работы [1-3].
В центре внимания данной работы находится задача стабилизации
движения подвижных объектов с учетом всего комплекса требований и
ограничений в целом, поэтому методы, описанные в указанных работах, не
могут быть применены с той же эффективностью. В качестве основной базы
для решения поставленной задачи можно выделить идею использования т.н.
многоцелевых структур, впервые описанную в работе [4]. Далее методология
синтеза была развита в работах [5;6], а в работах [7;8] показана
эффективность применения таких структур для реализации автопилота
морского судна. Одним из главных достоинств МС можно выделить тот факт,
что синтез отдельных частей регулятора может быть произведен в
определенном
смысле
последовательно
и независимо.
Это
позволяет
эффективно учесть все поставленные требования к динамике. Кроме этого,
независимость подразумевает возможность отключения различных компонент
регулятора в отдельных режимах движения, что может положительно
сказываться на вычислительной эффективности.
С другой стороны, при реализации систем управления на реальных
объектах с применением цифровых устройств неизбежно появляется проблема
наличия запаздывания. Вопросы анализа динамики систем с запаздыванием, а
также методы синтеза управления с учетом запаздывания также являются
одной из ключевых задач теории управления. Монография [9] является одной
из первых работ, посвященной развитию теории для решения указанной
12
проблемы. В книге [10] обсуждаются практические вопросы описанной теории.
Алгебраические методы исследования динамики систем с запаздыванием
приводятся в работе [11]. В монографии [12] рассмотрены вопросы
устойчивости по Ляпунову с точки зрения методов, развивающих идеология
функционалов
Ляпунова-Красовского.
теоретические
методы
анализа
В
систем
целом
с
можно
запаздыванием
сказать,
что
разработаны
достаточно глубоко.
В настоящее время активно развиваются вопросы, связанные с
практическим применением полученных теоретических результатов для
стабилизации систем с запаздыванием. Здесь можно выделить два подхода.
Первый подход систематизирован в книге [13]. Он основывается на
применении статических обратных связей для линейных систем. Вообще
говоря, у систем с запаздыванием бесконечно много собственных чисел.
Описанный подход преследует цель сдвига конечного числа крайних правых
собственных чисел системы в сторону левой открытой полуплоскости.
Движение осуществляется за счет итеративного процесса внесения небольших
изменений в вектор коэффициентов обратной связи. Явным недостатком этого
подхода является тот факт, что у предлагаемого процесса для каждой величины
запаздывания есть конечная граница, по достижении которой дальнейший
сдвиг крайних собственных чисел уже невозможен. При этом указанная
граница довольно часто находится в правой полуплоскости, так что
стабилизация возможна далеко не всегда. В работе также описан метод
преодоления указанной проблемы за счет использования асимптотических
наблюдателей особой структуры (по сути, за счет внесения дополнительной
динамики в регулятор). Однако зачастую на практике оказывается, что для
достижения
заданной
степени
устойчивости
системы
необходимо
синтезировать довольно значительной число таких наблюдателей. Иными
словами, динамика полученного регулятора может оказаться неоправданно
сложной и невыгодной с точки зрения вычислительной эффективности. При
этом спектр замкнутой системы остается бесконечномерным. Тем не менее,
13
применение указанного подхода может иметь ряд преимуществ с точки зрения
робастности, однако этот вопрос требует отдельного исследования.
Вторым подходом к синтезу стабилизирующих обратных связей для
линейных систем является компенсация запаздывания через использование
прогноза состояния системы. Впервые эта идея была описана в статье [14],
далее развита в монографии [15], а в работе [16] обобщена для нелинейных
систем и систем с распределенными параметрами. Одной из главных
особенностей компенсационного подхода является тот факт, что замкнутая
система имеет конечное число собственных чисел. В связи с этим
рассматриваемый подход в литературе именуется назначением конечного
спектра для систем с запаздыванием.
В настоящее время широкое развитие получили вопросы практической
реализации компенсирующих обратных связей. Кроме этого, важнейшей
задачей является обобщение компенсационного подхода для использования
исходных динамических регуляторов. В статье [17] описан подход к решению
этой задачи, при этом в центре внимания находится не только вопрос о
сохранении конечного спектра корней, но и проблема сохранения передаточной
матрицы эталонной системы. Указанная статья обосновывает подходы к
стабилизации подвижных объектов с учетом запаздывания и комплекса
требований к динамике, использованные в данной работе.
14
Глава 1. Многоцелевые регуляторы
1.1. Общая структура закона управления
В данной главе рассматриваются регуляторы специального класса,
получившие название многоцелевой структуры управления [4], позволяющие
эффективно решать задачи стабилизации линейных систем с учетом комплекса
требований и ограничений к динамике в различных режимах движения. Одной
из главных особенностей описываемой структуры является возможность
разбить исходную задачу на более простые подзадачи и решать их в
определенном смысле последовательно.
Рассмотрим линейные уравнения динамики подвижного объекта в общем
виде с учетом вектора измерений и добавлением линейного уравнения привода
рулей:
Ax Bδ Hd,
x
y Cx,
δ u.
(6)
Общая структура математических моделей законов управления для
реализации в составе системы управления движением подвижным объектом
включает следующие элементы:
z Az Bδ G ( y Cz )
ξ F ( p )(y Cz )
u μ z ν y ξ
– уравнение асимптотического наблюдателя;
– уравнение корректирующего устройства;
(7)
– уравнение скоростного управляющего сигнала.
(8)
Первый элемент решает две задачи. Во-первых, очевидно, наблюдатель
дает оценку состояния системы в тех случаях, когда полный вектор состояния
недоступен для непосредственного измерения, но система обладает свойством
наблюдаемости. Во-вторых, этот элемент влияет на свойство астатизма, т.е. на
качество подавления системой постоянных ветровых возмущений.
15
Корректирующее
устройство
(7)
отрабатывает
реакцию
на
периодические волновые возмущения. При этом, в зависимости от
поставленных
задач,
этот
элемент
может
выступать
либо
в
роли
оптимального динамического компенсатора, минимизирующего влияние
возмущения на выход системы, либо в роли оптимального динамического
фильтра, минимизирующего реакцию управления на возмущения, например,
для экономии топлива и сохранения ресурса износоустойчивости рулей.
Уравнение скоростного закона (8) является самым важным элементом,
определяющим значение управляющего сигнала для обеспечения общей
устойчивости системы, а также для удовлетворения поставленных требований
и ограничений по динамике переходных процессов. Вообще говоря, можно
перейти к использованию эквивалентного позиционного закона управления по
выходу наблюдателя в виде
u kz k 0δ ν 0 y ξ ,
где
(9)
k μ ( A GC) , k 0 μB , ν 0 μG ν .
(10)
Стоить отметить, что использования управляющего сигнала в виде (9)
является предпочтительным при переходе к учету запаздывания в следующей
главе. В связи с отмеченным обстоятельством далее будем использовать
управление именно в таком виде.
В представленных выше уравнениях
ν , а также передаточная матрица
ξ Em ,
F (s )
постоянные матрицы G , μ ,
корректирующего устройства
подлежат поиску в процессе оптимального синтеза, исходя из требований по
динамике объекта управления в условиях ветра и волнения.
Заметим, что уравнения (7) могут быть представлены в нормальной
форме вида
p αp β( y cz ),
ξ γ p m f ( y cz ),
n
где p E – вектор состояния корректирующего устройства,
1
(11)
ξ Em
– его
выходной вектор (обратим внимание на то, что его размерность должна
16
совпадать с размерностью вектора δ отклонений исполнительных органов), α
,
β
, γ,
mf
– постоянные матрицы, удовлетворяющие условию
γ (E n s α) 1 β m f F ( s)
1
С учетом (9) и (11), уравнения стабилизирующего закона управления с
многоцелевой структурой могут быть представлены в виде
z Az Bδ G (y Cz),
p αp β(y Cz ),
(12)
u kz k 0 δ ν 0 y ξ.
Отметим тот факт, что регулятор (12) подразумевает реализацию на
бортовом оборудовании, которое является цифровым устройством. В связи с
этим проведем дискретизацию уравнений (12) с некоторым периодом
дискретизации Td . В результате получим систему разностных уравнений
z[n 1] A d z[n] B d [n] G d (y[n] Cz[n]),
p[n 1] α d p[n] β d (y[ n] Cz[n]),
u[n] kz[n] k 0 δ[n] ν 0 y[n] ξ[n],
которая предназначена для реализации на борту подвижного объекта.
1.2. Поиск коэффициентов базового закона управления
В первую очередь при выборе настраиваемых параметров многоцелевой
структуры
необходимо
решить
задачу
поиска
коэффициентов
закона
управления в скоростном или позиционном варианте [18]. В свою очередь, эти
коэффициенты однозначно определяются матрицами базовой обратной связи
по состояния вида
u K x (x x * ) K δ ,
(13)
где x * E n – заданный фиксированный вектор командного сигнала, K x и K –
постоянные матрицы искомых коэффициентов. Необходимо отметить, что
обратная связь (13) реализуется непосредственно, однако при этом является
17
базовым законом управления для решения всех рассматриваемых ниже
функциональных задач.
Таким образом, первой задачей является поиск матриц K x и K ,
которые затем используются для однозначного нахождения матриц μ и ν в
скоростном законе (7) и далее для расчета матриц k , k 0 и ν 0 позиционного
закона управления по выходу асимптотического наблюдателя.
Поиск указанных матриц определяется требованиями к динамике
собственного движения объекта управления. Эти требования являются
наиболее приоритетными среди всех динамических требований, которые
предъявляются к качеству функционирования системы автоматического
управления, замкнутой обратной связью
с многоцелевой структурой.
Выполнение обозначенных требований является обязательным условием для
любого режима функционирования системы управления.
Будем
подвижного
понимать
объекта
под
такой
собственным
режим
его
движением
произвольного
функционирования,
осуществляется при нулевых начальных условиях
x(0) 0 , δ(0) 0
который
вектора
состояния и вектора органов управления соответственно, а также при
отсутствии внешних возмущений, т.е.
d(t ) 0 .
В этом режиме основной
задачей управления является переход из нулевого начального положения
e(0) 0
вектора управляемых переменных в некоторое заданное конечное
положение e z .
Для дальнейших выводов рассмотрим вспомогательную скалярную
функцию с нормой в смысле евклидова пространства
(t ) e(t )
заданную при
,
t [0, ) .
С целью формализации проблемы поиска матриц K x и K , на
движениях замкнутой нелинейной системы (1), (4), (13) введём функционал,
определяющий качество динамического процесса, в виде
J 1 J 1 (K )
m (K ) z
,
z
(14)
18
где
z ez
(t ) , K {K , K } , значениями которого являются
, m tsup
x
[ 0 , )
величины перерегулирования по вспомогательной переменной .
(t)
m
z
2
t
T
p
0
Рис. 1. Наглядное представление перерегулирования и быстродействия
Кроме этого зададим еще один функционал
T p T p (K ) inf t m : (t ) M ( z , ), t t m ,
(15)
где M ( z , ) : z / z , – фиксированное вещественное число
(обычно 0.03 ). Функционал (15) определяет быстродействие замкнутой
системы, по сути – время, за которое
(t )
окрестность
не
командного
сигнала
и
перейдет в некоторую покинет
ее.
На
рис.
1
продемонстрировано схематичное представление сущности приведенных
функционалов (14) и (15).
Выбор
настраиваемых
коэффициентов
должен
обеспечивать
одновременную минимизацию приведенных функционала на множестве k ,
на котором обеспечивается только устойчивость замкнутой системы, либо на
таком множестве, чтобы корни i характеристического полинома 0 ( s )
19
замкнутой линейной системы (2), (4), (13) находились в некоторой заданной
области C комплексной плоскости:
k K : i C , 0 ( i ) 0, i 1, n .
Существуют различные способы задания области C . В качестве
примера можно привести область [18], условно изображенную на рис. 2,
которую можно формально представить в виде
C C 2 s x yj C1 : x , 0 y ( x ) tg ,
где величина определяет заданную степень устойчивости, а величина
–
заданную степень колебательности замкнутой системы.
Im s
C
Re s
Рис. 2. Пример задания допустимой области расположения корней
Стоит отметить тот факт, что на практике чаще всего одновременное
достижение минимума функционалов (14) и (15) невозможно, поэтому в
действительности коэффициенты базового закона (13) выбираются с учетом
ограничения значений приведенных функционалов:
J 1 (K ) J 10 , T p (K ) T p 0 ,
где постоянные числа J 10 и
Tp 0
(16)
задаются из практических соображений.
Вопрос об обеспечении ограничений (16) за счёт выбора матрицы K
может быть, в принципе, сведён к задаче оптимального параметрического
синтеза [18]. В данном случае эта задача является одним из сложных вариантов
20
проблемы нелинейного программирования, отягощенной сложным заданием
целевой функции и допустимого множества искомых матриц.
В настоящее время для поиска указанных матриц с учетом заданных
требований чаще всего используются следующие методы:
– направленный перебор коэффициентов закона (13);
– поиск коэффициентов базового закона (13) как решения задачи
модального управления, например – с биномиальным распределением
собственных значений;
– поиск закона (13) как решения задачи LQR-оптимизации;
– поиск коэффициентов в процессе параметрической оптимизации с
применением специализированного пакета системы MATLAB.
После нахождения матриц K x и K осуществим переход к скоростному
закону управления. Этот переход обусловлен наличием требования астатизма
для динамики системы. В скоростном законе (8) присутствуют выходные
переменные асимптотического наблюдателя, но предварительным шагом здесь
является построение скоростного закона по фактическому состоянию системы:
(17)
u μ x ν y
Астатизм в классической теории управления достигался за счет введения
в регулятор дополнительной интегрирующей составляющей, как это делается,
например, в стандартном пропорционально-интегрально-дифференциальном
(ПИД) законе управления. Однако в таком случае могут существенно
ухудшиться динамические характеристики качества процессов управления.
Перебор коэффициентов интегральной составляющей, как правило, может
улучшить только отдельные характеристики, но в целом качество динамики
будет хуже.
Указанную проблему можно решить следующим образом. Приведем
общую процедуру вывода коэффициентов скоростного закона управления (17)
через матрицы базовой обратной связи (13) [16;18]. Пусть матрицы базового
стабилизирующего
управления
(13)
известны.
Рассмотрим
замкнутой им системы в линейном приближении в виде
уравнения
21
x Ax Bδ , e Lx , δ u ,
(18)
u K x x K δ ,
(19)
где e – вектор контролируемых переменных, при нулевом командном сигнале
и при отсутствии внешнего возмущения. Будем считать, что в данном случае
dim e dim u m , причем rank L m . Тогда выражение для вектора
e можно
переписать в виде
e L a x a L e x e , L L a
где dim x e m ,
L e , x x a
xe ' ,
(20)
dim x a n m , причем матрица L e является невырожденной,
благодаря чему можно выразить вектор x e в виде
x e Le1e Le1L a x a .
(21)
Теперь подставим (21) в уравнение (18) с учетом (20):
x ( A 2 A1Le1L a )x a Bδ A1Le1e .
(22)
Введем обозначения
γ x a
δ ,
T
S A 2 A1Le1L a B .
Тогда из (22) имеем
γ S 1 x S 1 A 1Le1e .
(23)
С учетом введенных обозначений уравнение (19) можно записать в виде
u Mγ me ,
(24)
1
где M K xa K xe L e L a K , m K xe Le1 . Далее подставив выражение (23) в
(24), получим скоростной базовый закон:
u MS 1x (MS 1A1Le1 m)e .
Если матрица
(25)
MS 1 A1Le1 m не вырождена, то (25) обеспечивает
астатизм для системы (18). Подставим (25) в уравнения привода:
δ MS 1x (MS 1A1Le1 m)e .
(26)
В статическом положении равновесия с учетом возмущения имеем x 0 ,
0 , тогда
δ
(MS 1 A1Le1 m)e 0 .
(27)
22
Рассматривая (27) как систему линейных уравнений и учитывая
невырожденность матрицы MS 1A1Le1 m , получим единственное решение
e 0 . Следовательно, закон управления (25) действительно обеспечивает
астатизм для линейной системы (18) по вектору e .
Отдельно отметим тот факт, что тождественность преобразований в силу
уравнений объекта при переходе к астатическому регулятору (25), обеспечивает
также и тождественность динамики собственного движения для закона
управления (13) и закона управления
u MS 1 x (MS 1 A1Le1 m)e K x x * .
Последнее свидетельствует о том, что ограничения (16) для базового
закона (13) сохраняются и для скоростного астатического закона (25).
В итоге можно заключить, что если контролируемый и измеряемый
вектор совпадают, т.е. если e y , то коэффициенты скоростного закона (8) или
(17) определяются по формулам
μ MS 1 , ν (MS 1 A1Le1 m) .
(28)
Стоит отметить, что общий подход к выводу формул (28) в частных
случаях может быть существенно упрощен. Далее в соответствующей главе
будет продемонстрирована процедура для поиска матриц μ и ν в задаче
стабилизации курсового угла СВП.
2.2. Поиск коэффициентов асимптотического наблюдателя
Следующим элементом, подлежащим поиску в процессе синтеза
многоцелевой
структуры,
является
матрица
G
коэффициентов
асимптотического наблюдателя. Здесь поиск определяется требованиями к
динамике объекта управления в режиме воздействия внешних ступенчатых
возмущений (как правило, за счет влияния ветра). При этом отметим, что
командный сигнал в этой ситуации не входит в закон управления, т.е. имеем
условие
x * 0 .
23
Требования к динамике в данном режиме имеют более низкий
приоритет по сравнению с требованиями к режиму собственного движения,
однако все еще остаются одними из наиболее важных условий.
Внешние возмущения, воздействующие на объект управления в
рассматриваемом режиме, определяются равенством
d(t ) d 0 1(t ) ,
где
1(t )
(29)
– ступенчатая единичная функция, d 0 E nd – постоянная векторная
величина. При этом, в данном режиме также выполняются равенства
δ(0) 0 ,
x ( 0) 0 ,
т.е. начальные условия по векторам состояния объекта управления и
исполнительных органов являются нулевыми.
Будем
считать,
астатическим
законом
что
обратная
управления
связь
с
реализуется
многоцелевой
позиционным
структурой
при
отключенном динамическом корректоре, при этом уравнения регулятора
принимают вид
z Az Bδ G ( y Cz ),
u kz k 0 δ ν 0 y ,
(30)
где матрицы k , k 0 и ν 0 вычисляются по формулам (10) через матрицы
скоростного управляющего сигнала μ и ν , найденные на предыдущем этапе.
Также будем считать, что величина d 0 задана, причем для нее
существует положение равновесия в системе (1),(4),(30).
Переходный процесс по вспомогательной переменной
(t )
в режиме
ступенчатых возмущений с регулятором (30) с учетом астатических свойств
последнего может быть схематично представлен на рис. 3.
0
t
(t , G )
p
24
Рис. 3. Наглядное представление переходного процесса при астатизме
Задача о поиске матрицы G также может быть сведена к задаче
оптимизации. Для этого зададим функционал на движениях системы (1),(30),
значения
которого
определяются
максимальным
отклонением
по
регулируемым координатам:
p p (G ) max (t , G ) .
(31)
t[ 0 , )
Фактически функционал (31) определяет максимальный «провал»
регулируемых координат при воздействии ступенчатого возмущения. В
соответствии с вышесказанным матрица G выбирается таким образом,
чтобы обеспечить минимум функционала (31) и устойчивость наблюдателя,
не ухудшив при этом качество управления в предшествующем режиме
движения. Решение указанной задачи может производиться методами,
описанными в предыдущем параграфе, а также путем синтеза наблюдателя
как фильтра Калмана.
2.3. Синтез корректирующего устройства
Последним элементом, подлежащим поиску при решении задачи
синтеза многоцелевой структуры, является передаточная матрица
F (s )
корректирующего устройства. Вообще говоря, поиск этой матрицы зависит от
поставленных требований к динамике: максимальная точность стабилизации
или минимальная интенсивность управляющего воздействия при движении
объекта управления в условиях воздействия периодических возмущений.
В данной работе рассматривается второе требование. Целью снижения
интенсивности функционирования рулей может являться экономия ресурса
привода и снижение потребления топлива. В качестве дополнительного
эффекта при этом может произойти выход из режима автоколебаний,
25
порожденных нелинейностями привода и, как следствие, – может быть
достигнуто дополнительное повышение точности стабилизации.
Для формализации указанных задач рассмотрим блок-схему [18]
замкнутой линейной системы управления (18),(12), изображенную на рис. 4.
Рис. 4. Блок-схема замкнутой линейной системы
Из
рис.
4
и
уравнений
(12)
очевидно,
для
системы
(18),
представляющей модель объекта управления, регулятором является канал
исполнительных органов, который в свою очередь является LTI-системой,
состоящей из уравнения приводов (4) и многоцелевой структуры (12). На
вход регулятора подаются текущие измерения y , а выходной переменной
является вектор
δ
отклонения рулей. При этом сама LTI-система,
являющаяся моделью канала, может рассматриваться, как система, замкнутая
регулятором в виде корректирующего устройства. В качестве объекта
управления
здесь
выступает
основная
часть
уравнениями
u,
δ
z Az Bδ G(y Cz),
u μ z ν y ξ ,
которые можно преобразовать к эквивалентному виду
канала,
описываемая
26
z Az Bδ G (y Cz ),
δ kz k δ ν y ξ,
0
где
(32)
0
k μ ( A GC) , k 0 μB , ν 0 μG ν .
Рис. 4 показывает, что входной переменной для основной части канала
является вектор
вектор
δ
,
ζ
где
y
.
ξ
Будем считать выходным вектором основной части
ζ y Cz .
Добавим к системе (34) уравнения выхода, тогда в
пространстве состояний модель основной части канала исполнительных
органов будет представлена уравнениями
z A GC B z G 0 nm y
,
k
k 0 δ ν 0 E m ξ
δ
δ 0 mn E m z 0 mk 0 mm y
.
C
0
E
0
ζ
δ
k m
k
k m ξ
(33)
Для вывода дальнейших результатов представим (33) в частотной
области, при этом получим систему
δ
y T ( s ) T12 ( s ) y
,
T( s ) 11
ζ
ξ T21 ( s ) T22 ( s ) ξ
где
T(s )
(34)
– передаточная матрица основной части канала исполнительных
органов. Эта матрица может быть в явном виде выражена через постоянные
матрицы системы (33) следующим образом:
E m
0
E nm s
T( s ) mn
C 0 km
0 mm
0
.
mk
0 km
Ek
В
соответствии
с
рис.
A GC B
k
k 0
4,
замкнем
1
G
ν0
0 nm
E m
систему
(33)
или
(34)
корректирующим устройством. При этом получим уравнения замкнутой
системы в виде
δ T11 ( s )
ζ T21 ( s )
ξ F ( s)ζ.
T12 ( s ) y
,
T22 ( s ) ξ
(35)
27
Fy (s )
Найдем передаточную матрицу
канала исполнительных органов
от входа y к выходу δ , используя уравнения (35). Для этого исключим
внутренние переменные
ζ
и
ξ
из этой системы и получим
Fy ( s ) T11 ( s ) T12 ( s )F( s ) E k T22 ( s )F( s ) T21 ( s ) .
1
При этом связь между вектором y измерений и вектором δ состояния
исполнительных органов определяется уравнением вход-выход
δ Fy (s ) y ,
которое является математической моделью в частотной области канала
исполнительных органов как LTI-системы.
Теперь рассмотрим функционал
J 2 J 2 (F ) Fy ( s, F )
заданный на множестве
,
F (s )
матриц
с правильными дробно-
рациональными компонентами, корни знаменателей которых лежат в левой
открытой
полуплоскости.
интенсивности
управления
Тогда
задача
может
о
быть
достижении
минимальной
сформулирована
в
виде
оптимизационной задачи
J 2 (F ) min , J 20 min J 2 (F ) ,
F
F
F( s) arg min J 2 (F) .
F
(36)
Стоит отметить, что решение указанных оптимизационных задач
можно быть найдено при помощи методов многоцелевого синтеза или Hоптимизации. Однако при практической реализации таких методов могут
возникать существенные трудности.
Как правило, в частных случаях можно существенно упростить задачу
оптимизации за счет введения ряда допущений. Например, в случае морских
подвижных
объектов
зачастую
в
качестве
приближенной
модели
периодических возмущений принимают гармонические функции. Решение
задачи синтеза корректора для этого случая будет показано далее в
28
соответствующей главе на примере задачи коррекции морского волнения для
стабилизации курсового угла СВП.
29
Глава 2. Учет запаздывания
Как было отмечено во введении, запаздывание неизбежно возникает при
реализации управления на реальных подвижных объектах. Особенно заметно
проявление запаздывания при работе с цифровыми системами управления. При
этом в настоящее время именно такие системы наиболее часто применяются на
практике. Как правило, запаздывание оказывает негативное влияние на
динамику системы, вплоть до потери устойчивости. Рассмотрим один из
подходов к учету наличия запаздывания в контуре управления. Основная идея
здесь состоит в компенсации запаздывания за счет прогноза состояния системы
и, соответственно, сохранения динамических характеристик некоторой
исходной системы управления без запаздывания.
2.1. Применение прогноза для компенсации запаздывания
Для начала рассмотрим предлагаемый подход в простейшем случае.
Введем в систему (2) постоянное запаздывание h 0 , для удобства будем
считать, что возмущение подается в систему также с запаздыванием. Тогда
имеем систему вида
(t ) Ax(t ) Bu (t h) Hd(t h),
x
y Cx.
(37)
Будем считать, что выходной вектор y совпадает с вектором состояния
системы
x , при этом векторы u ,
d и y подлежат непосредственному
измерению. Наряду с системой (39) будем рассматривать вспомогательный LTI
объект
Ax Bu Hd(t h),
x
y Cx.
(38)
Для системы (38) построим стабилизирующее управление вида
u Kx,
(39)
30
обеспечивая гурвицевость матрицы A BK . При этом уравнения замкнутой
системы (38), (39) приобретают вид
x A BK x Hd(t
h),
(40)
y Cx.
Далее обратимся к исходной системе с запаздыванием (37) поставим
задачу о нахождении управления
u ( y , d)
(41)
таким образом, чтобы замкнутая система (39), (41) приняла вид (40). В
научной литературе такая задача получила название задачи назначения
конечного спектра для систем с запаздыванием [13]. Очевидно, что решение
задачи можно получить [16], взяв управление (41) в виде
u(t h) Kx(t ) ,
или, что эквивалентно
u(t ) Kx(t h) .
(42)
В общем случае прямая реализация управления (42) физически
невозможна, поскольку для этого требуется знание будущего состояния
объекта управления в момент t h . Однако тот факт, что моделью исходного
объекта является LTI-система, позволяет использовать формулу Коши для
прогноза состояния системы через текущий вектор
d (t ) ,
заданные на отрезке
x(t )
и векторы
u (t )
и
t [t h, t ] :
Ah
x(t h) e x(t ) e
t h
Ah
e
A ( t )
Bu( h) Hd( h) d .
(43)
t
Теперь введем в рассмотрение вспомогательную динамическую
векторную переменную z E n :
t
z (t ) e At e A Bu() Hd() d .
0
Преобразуем формулу (43) с учетом тождества
t
t h
t
0
0
t h
A
A
A
e Bu() Hd() d e d e d ,
из которого следует, что
31
t h
x(t h) e Ah x(t ) e Ah e A ( t ) Bu( h) Hd( h) d
t
t
e Ah x(t ) e At e A Bu() Hd() d
t h
t
e Ah x(t ) e At e A d
0
t h
e d e
A
Ah
x(t )
0
t
t h
0
0
e At e A d e Ah e Ah e At e A d e Ah x(t )
t
e
At
e Bu() Hd() d
A
t h
Ah
e e
A (t h)
0
Ah
e Bu() Hd() d
A
0
Ah
e x(t ) z (t ) e z (t h).
Таким образом, формула (43) представлена в виде
x(t h) z (t ) e Ah x(t ) z (t h) ,
(44)
что позволяет переписать формулу обратной связи по прогнозу (42)
следующим образом:
u e Ah x Kz Ke Ah z (t h) .
(45)
Прямая реализация (45) сопряжена с рядом значительных проблем,
связанных
с
численным
вспомогательной переменной
взятием
интеграла,
представляемого
z , как отмечено, например, в работе [13]. В
соответствии с этим трансформируем полученное представление прогноза
(44). Найдем производную вспомогательной переменной
z:
t
z (t ) Ae At e A Bu() Hd() d e At e At Bu(t ) Hd(t )
0
(46)
Az (t ) Bu(t ) Hd(t ).
С учетом выражения (46), по сути, вместо обратной связи (45) получаем
регулятор с собственной динамикой
z Az Bu Hd,
u e Ah x Kz Ke Ah z (t h),
практическая реализация которого не представляет трудностей.
2.2. Учет запаздывания для динамических регуляторов
В предыдущем параграфе был рассмотрен подход к компенсации
запаздывания при стабилизации исходных объектов без запаздывания
32
статическими регуляторами. Однако в центре внимания данной работы
является применение многоцелевой структуры, которая по своей сути
является динамическим регулятором, как было показано выше. В связи с
этим рассмотрим обобщение метода синтеза компенсирующей обратной
связи. Ключевым моментом, обеспечивающим возможность применения
такого обобщения к задачам управления подвижными объектами, является
следующий факт. В большинстве случаев линейная система, приближенно
описывающая движение подвижного объекта (особенно в случае МПО)
является устойчивой, однако асимптотическая устойчивость отсутствует
ввиду наличия нулевых собственных чисел. Тем не менее, матрицу A можно
считать гурвицевой с малой степенью устойчивости, поскольку динамика
объекта компенсируется в результате использования прогноза. Далее так же
будем считать, что для измерения доступен полный вектор состояния
системы.
Рассмотрим вспомогательный объект без запаздывания в канале
управления с моделью
x Ax Bu Hd(t h),
η Cx.
(47)
Представим регулятор (12), полученный в результате синтеза элементов
многоцелевой структуры для объекта управления, в общем виде
u W ( p ) η ,
где
p d / dt ,
(48)
причем общие знаменатели дробно-рациональных элементов
передаточной матрицы W тождественно не равны единице.
В уравнениях пространства состояний регулятор (50) представляется
следующим образом:
ξ A k ξ B k η,
u C k ξ D k η,
где
ξ E –
(49)
вектор состояния регулятора. Матрицы представления (49)
удовлетворяют тождеству
W ( s ) C k E s A k B k D k ,
1
(50)
33
где
E –
единичная матрица размера .
Исключим из уравнений замкнутой системы (47), (49) переменную u :
x A BD k C x BC k ξ Hd(t h),
ξ B Cx A ξ.
k
(51)
k
Теперь вернемся к системе с запаздыванием (37) и построим для нее
обратную связь вида
u(t ) W ( p)Cx(t h) ,
(52)
которая приводит к устойчивой замкнутой системе (51). Как было показано
выше, прогноз можно выполнить по формуле (44):
x(t h) e Ah x(t ) z (t ) e Ah z (t h) ,
где векторная функция
z (t )
удовлетворяет уравнению
z Az Bu Ηd .
Подставляя прогноз (44) в уравнение (52), получаем
u(t ) W ( p )C e Ah x z e Ah z (t h) ,
т.е. приходим к компенсирующему регулятору по выходу.
Теперь представим полученную обратную связь в виде уравнений
пространства состояний с учетом тождества (50). Введем в рассмотрение
вспомогательную переменную
γ C e Ah x z e Ah z (t h) ,
будем считать ее входом регулятора (49). Дополним модель регулятора
z (t ) ,
в итоге получим обратную
γ Ce Ah x z e Ah z (t h) .
(53)
уравнением динамики векторной функции
связь вида
z Az Bu Hd,
ξ A ξ B γ ,
k
k
u C k ξ D k γ ,
Таким образом, получен обобщенный динамический регулятор,
позволяющий компенсировать запаздывание в системе управления и
получить динамические свойства исходной системы без запаздывания.
34
Заметим, что полученный регулятор требует измерения всех компонент
вектора состояния. Вообще говоря, можно показать, что существует метод,
позволяющий компенсировать запаздывание регулятором по выходу, однако
такой подход ориентирован на сохранение динамики вспомогательного
объекта в виде
x Ax Bu Hd (t h),
η C1x, C1 Ce Ah ,
замкнутого динамическим регулятором
~
u W( p) η .
Однако в центре внимания данной работы находится компенсация
запаздывания с целью сохранения динамики системы без запаздывания (38),
замкнутой
многоцелевым
управление
в
виде
(53).
регулятором,
Конкретная
поэтому
будем
структура
использовать
компенсирующего
многоцелевого регулятора подушке представляется в виде
z p Az p Bδ Hd,
z Az Bδ m G ( γ k Cz ),
δ m kz k 0δ m ν 0 γ k γp m f ( γ k cz),
p αp β( γ k Cz ),
γ k c e Ah x z p e Ah z p (t h) ,
δ z δ m ,
где δ m – модельное отклонение рулей, δ z – текущее вычисленное положение
рулей, необходимое для достижения поставленных целей управления, т.е. для
стабилизации с учетом запаздывания.
Обратим внимание на следующий факт. Уравнения (53) подразумевают,
что управляющий сигнал на выходе регулятора подается
напрямую на
объект управления. Однако в действительности при многоцелевом подходе
управлению напрямую подлежит скорость перекладки рулей, а не отклонение
непосредственно. В связи с отмеченным обстоятельством для обеспечения
асимптотического стремления фактического положения δ рулей к их
вычисленному положению δ z на привод рулей можно подать управляющий
сигнал в виде
u k u (δ z δ) ,
35
приводя его уравнение к виду
δ k u (δ z δ) ,
где k u 0 – диагональная матрица с достаточно большими вещественными
компонентами для быстрого завершения динамических процессов в
приводах.
Далее введем обозначения
A GC B
A c k m f C k 0
βC
0 21
G
B c ν 0 m f
β
032
γk ,
α
,
Cc 013 1 012 ,
D c 0.
С учетом введенных обозначений компенсирующий регулятор для
стабилизации произвольного подвижного объекта с учетом постоянного
запаздывания при многоцелевом подходе можно переписать в обобщенном
виде
z p Az p Bδ Hd,
ξ c A c ξ c B c γ k ,
γ k C e Ah x z p e Ah z p (t h) ,
δ z C c ξ c ,
u k u (δ z δ).
(54)
В соответствии с соотношениями (38), (54) на рис. 5 приведена блоксхема замкнутой системы управления с компенсирующим регулятором.
36
Рис. 5. Замкнутая система управления с компенсирующим
регулятором (58)
Далее, ориентируясь на тот факт, что компенсирующая обратная связь
должна быть реализована на бортовом оборудовании, которое является
цифровым устройством, проведем дискретизацию полученного регулятора с
некоторым периодом дискретизации Td . Здесь стоит отметить, что в
дискретном
варианте
экспоненты
при
вычислении
переменной
γk
необходимо заменить на степени матрицы A . В результате получим
дискретный регулятор
z p [ n 1] A d z p [ n] B d δ[ n] H d d[ n],
ξ c [ n 1] A cd ξ c [ n] B cd γ k [ n],
γ k [ n] C A n x[ n] z p [ n] A n z p [ n nh ],
h
h
δ z [ n] C c ξ c , [ n]
u[ n] k u (δ z [n] δ[ n]),
(55)
где nh [h Td ] .
В таком виде синтезированный регулятор подлежит непосредственной
реализации на бортовом оборудовании объекта управления. Отметим тот
факт,
что
дискретная
природа
полученных
уравнений
обеспечивает
достаточно простую и эффективную с вычислительной точки реализацию на
цифровых устройствах. Более того, можно говорить об использовании
полученной обратной связи в системах реального времени. Конечно,
необходимо принять к сведению доступные вычислительные ресурсы и
величину периода дискретизации, однако в целом для эффективного
37
функционирования дискретного регулятора (55) достаточно относительно
небольших вычислительных мощностей.
38
Глава 3. Объект управления
3.1. Суда на воздушной подушке
В
качестве
описанных
объекта
подходов
в
управления
частном
для иллюстрации
случае
рассмотрим
применения
морские
суда
специального класса, поднимающиеся над водной поверхностью при помощи
аэростатического давления между корпусом судна и поверхностью воды. В
англоязычной литературе такие суда относят к ACV (Air Cushion Vehicle) или
hovercraft типу, а в русскоязычных изданиях их именуют судами на
воздушной подушке (СВП). При этом в данной работе рассматриваются СВП
амфибийного типа, представляющие наибольший интерес.
Рис. 6. Образование избыточного давления под днищем СВП
Схема образования воздушной подушки представлена на рис. 6.
Подъемная сила, компенсирующая вес судна, создается вентиляторами,
нагнетающими воздух в зону повышенного давления 3 (воздушная подушка,
ВП). Для более равномерного распределения давления, воздушный поток
может предварительно проходить через зоны жесткого ресивера 1 и гибкого
ресивера
2.
Воздушная
подушка
замкнута
гибким
ограждением,
конструктивно объединенным с корпусом. Благодаря гибкому ограждению
39
жесткий корпус СВП поднят над водой на достаточную высоту, при которой
исключается непосредственный контакт даже при развитом волнении моря.
Тем самым обеспечиваются амфибийные свойства СВП, его способность
двигаться не только над водой, но и над твердой поверхностью, выходить на
лед, берег и т.п.
Поскольку СВП не имеют контакта с водной поверхностью, их
движители располагают на палубе. Они представляют собой воздушные
винты с изменяемым шагом. Изменение направления движения СВП может
осуществляться при помощи установки разного шага винтов (в случае, если
движителей несколько), а также вертикальных аэродинамических рулей. На
малых скоростях хода СВП более эффективно использование разности шага
винтов, но при наборе скорости эффективность снижается, выводя на первый
план использование аэродинамических рулей. В некоторых моделях СВП
реализована возможность поворачивать непосредственно сами движители,
однако такая практика распространена в меньшей степени. В рамках данной
работы будем считать, что имеется единственный движитель, а управление
поворотом осуществляется только при помощи аэродинамических рулей.
Так как судно движется в воздухе, присоединенные массы воды и
соответствующие моменты инерции себя не проявляют. Однако силы и
моменты
гидродинамической
природы
существуют
в
результате
взаимодействия между воздушной подушкой и поверхностью воды. При
малых скоростях под днищем СВП образуется углубление в водной
поверхности. После преодоления так называемой «горбовой» скорости
(обычно более 30 узлов) углубление в водной поверхности перестает
образовываться, вследствие чего уровень гидродинамических сил и моментов
заметно снижается [18].
Далее приведем математическую модель СВП, как в виде системы
обыкновенных дифференциальных нелинейных уравнений с учетом всех сил и
моментов, так и в линейном приближении, которое далее будет использовано
для синтеза управления.
40
3.2. Уравнения нелинейной динамики
Рис. 7. Система координат
Нелинейные дифференциальные уравнения динамики СВП в боковом
движении можно записать в виде [20; 21]:
m(V x Vz ω y ) Fxa Fxh T ,
m(V z Vx ω y ) Fza Fzh Fzc Fzr ,
x M xa M xh M θ M xr ,
I xω
y M ya M yh M yr ,
I yω
или в виде обыкновенных дифференциальных уравнений
V x ( Fxa Fxh T Vz ω y )/m,
V z ( Fza Fzh Fzc Fzr Vx ω y )/m,
x ( M xa M xh M θ M xr )/I x ,
ω
y ( M ya M yh M yr )/I y ,
ω
где
Iy
V x , V z , ω x
-
и
моменты
ωy –
(56)
линейные и угловые скорости, m – масса судна, I x и
инерции.
Индекс
a
обозначает
аэродинамическое
сопротивление, h – гидродинамическое сопротивление, r – силы и моменты,
действующие на рули, Fzc – сила реактивной струи потока воздуха,
выходящего
из-под
воздушной
подушки
в
случае
крена,
Mθ
–
восстанавливающий момент силы давления воздуха внутри воздушной
41
подушки, T – сила тяги движителя. Соответствующая система координат
представлена на рис. 7.
Выражения для сил и моментов аэродинамического сопротивления
представляются в виде
Fxa 0.5C xa (β a )ρ aV02 S ,
Fza 0.5C za (β a )ρ aV02 S ,
M xa 0.5Fza H ,
ω
M ya 0.5(C mya (β a ) C myya
ωy L
V0
)ρ aV02 SL,
где β a – угол аэродинамического дрейфа (относительно ветра), ρ a –
плотность воздуха,
V0
– скорость судна относительно ветра, S – площадь
боковой поверхности судна, L и H – длина и высота судна, соответственно;
ω
C xa (β a ), C za (β a ), C mya (β a ), C myya – аэродинамические коэффициенты.
Гидродинамическое сопротивление имеет форму
Fxh 0.5C xh ( Fr , β)ρ wV 2W 2/3 ,
Fzh 0.5C zh ( Fr , β)ρ wV 2W 2/3 ,
M xh 0.5Fzh H ,
ω
M yh 0.5(C myh ( Fr , β) C myyh ( Fr )
где
β
ω yW 1/3
V
)ρ wV 2W ,
– угол дрейфа судна, ρ w – плотность воды, V – скорость судна, W –
водоизмещение судна. C xh (Fr , β), C zh (Fr , β), C mxh (Fr , β),
ωy
C myh (Fr , β), C my
h
–
гидродинамические коэффициенты. Fr – число Фруда:
Fr
V
.
Lg
Число Фруда является безразмерной величиной, позволяющей в
судостроении сравнивать гидродинамику моделей разного размера.
Стоит отметить, что из-за сложного профиля корпуса для реальных
судов, как правило, аэродинамические и гидродинамические коэффициенты
не подлежат теоретическому расчету, а измеряются в процессе экспериментов
в аэродинамической трубе и бассейне.
Точный расчет сил и моментов, действующих на рули, так же затруднен
в связи со сложностью протекания аэродинамических процессов, поэтому
42
здесь используется некоторое приближение, полученное в результате
проведения экспериментов и описанное в книге [20]:
ρa 2
V0 S ,
2
Fzr x r ,
Fzr C r δ
M yr
M xr Fzr y r ,
C r C r1 [(1 k r )(1 s )
kr
S r'' S r'
],
Sr Sr
1.9π1r s
,
r (1.9π λ1r )( s 1)
1 2.53T
1,
(ρ a D 2V02 )
s
где C r – силовой коэффициент, δ – угол отклонения руля, x r и y r –
координаты центра руля, C r1 – коэффициент модели без движителей. S r –
площадь поверхности руля, S r' ,S r'' – области площади руля внутри и вне
воздушного потока, λr ,λr1 – соответствующие части длины рулей, D –
диаметр воздушного винта. k r и s – вспомогательные коэффициенты.
Для вычисления силы воздушного потока из-под гибкого ограждения
воздушной подушки, а также восстанавливающего момента воздушной
подушки в условиях ненулевого угла крена, вообще говоря, необходимо
строить и решать динамические уравнения давления в области воздушной
подушке. Однако по тем же причинам, что обосновывают использование
приближенных формул расчета сил, действующих на рули, обычно для целей
моделирования
формулы
с
достаточно
рассмотреть
коэффициентами,
упрощенные
получаемыми
в
приближенные
процессе
проведения
экспериментов с моделью:
Fzc mgkθ θ,
M θ mghθ N θ ω x ,
где θ – угол крена, h – метацентрическая высота, k θ – коэффициент силы
потока,
Nθ
– коэффициент восстанавливающего момента воздушной
подушки. Коэффициенты k θ и N θ для реальных судов, как правило,
измеряются
экспериментально,
так
гидродинамические коэффициенты.
же
как
и
аэродинамические
и
43
3.3. Линейная модель объекта управления
Линеаризуем систему нелинейных дифференциальных уравнений (6) в
окрестности положения равновесия V x const ,
Vz y r 0
полагая, что
внешние возмущения отсутствуют. В этом режиме очевидным образом можно
пренебречь динамикой СВП по продольной скорости. Также ввиду малости
можно пренебречь динамикой по углу крена. Кроме этого вместо боковой
скорости Vz в линейных моделях МПО принято рассматривать угол дрейфа,
V
z
который в данной ситуации представляется как β V .
x
С учетом вышесказанного линейная модель динамики СВП в боковом
движении, используемая для стабилизации курсового угла, представляется
следующими уравнениями:
β a11β a12 y b1 r h1d ,
y a21β a22 y b2 r h2 d ,
y ,
где d – скалярное внешнее возмущение (ветровое или волновое).
44
Глава 4. Практическая реализация многоцелевого
регулятора
4.1. Коэффициенты модели
В качестве конкретных значений параметров нелинейной модели СВП
использованы данные гипотетического судна, являющиеся осреднением
данных, принимаемых для различных судов в работах [19-22]. Основные
характеристики судна приведены в табл. 1.
Характеристика
Объёмное водоизмещение
Длина
Ширина
Высота
Масса
Момент инерции
относительно продольной
оси
Момент инерции
относительно вертикальной
оси
Момент инерции
относительно вертикальной
оси
Максимальный угол
перекладки рулей
Максимальная скорость
перекладки рулей
Максимальная скорость
Площадь руля
Продольное плечо руля
Вертикальное плечо руля
Плотность воды
Плотность воздуха
Размерность
м3
м
м
м
кг
кг м2
Величина
3 628.8
27
14
9.6
70 000
1 680 933.3
Обозначение
W
L
B
H
m
Jxx
кг м2
5 395 833.3
Jyy
кг м2
4 790 100
Jzz
град
40
max
град/с
15
d/dt
м/с
м2
м
м
кг/м3
кг/м3
28
5.4
13.5
2.4
1000
1.292
Vmax
Sr
xr
yr
w
a
Таблица 1. Характеристики рассматриваемой модели СВП.
45
Как было описано выше, линейная модель бокового движения судна в
отклонениях
от
положения
равновесия
Vx const ,
Vz ω y θ δ 0
и
отсутствии внешних возмущений и с учетом уравнения динамики привода
руля представляется в виде
β a11β a12 ω y b1δ h1 d ,
y a21β a22 ω y b2 δ h2 d ,
ω
(57)
ω y ,
δ u.
Для проверки описанных подходов была реализована компьютерная
модель нелинейной динамики СВП в среде Matlab-Simulink. В результате
идентификации линейной модели судна при прямолинейном движении с
продольной скоростью около 12 м/с были получены следующие значения
коэффициентов:
a11 0.3176, a12 0.852,
a12 0.0102, a22 0.1383,
b1 0.005, b2 0.0217,
h1 0.0006, h2 0.00025.
Сравним динамику линейной и нелинейной систем. На рис. 8 показана
реакция систем на единичное ступенчатое отклонение руля. В этой ситуации
модели совпадают практически идеально. Рис. 9 демонстрирует динамику
моделей при воздействии ступенчатого возмущения в виде бокового ветра со
скоростью 10 м/с. Здесь можно заметить, что на первых секундах
присутствует
небольшое
расхождение
по
боковой
скорости,
однако
рассмотрение динамики в целом позволяет сказать, что и в этой ситуации
имеем достаточное совпадение моделей.
Таким образом, можно сделать вывод о том, что полученная линейная
модель является довольно хорошим приближением для динамики исходной
нелинейной системы в принятом режиме движения и может быть
использована для дальнейшего решения задач синтеза многоцелевого и
компенсирующего регуляторов.
46
Рис. 8. Сравнение линейной и нелинейной моделей при единичном отклонении руля
Рис. 9. Сравнение линейной и нелинейной моделей при воздействии бокового ветра
47
4.2. Режим отработки командного сигнала
Синтезируем первый элемент многоцелевой структуры управления для
системы (57). Поскольку на этом этапе задача синтеза наблюдателя еще не
решена, будем использовать базовый закон управления в скоростном
варианте по состоянию системы.
Помимо очевидного требования по обеспечению асимптотической
устойчивости этот элемент должен обеспечивать желаемое качество
управления.
Прежде
всего,
величина
перерегулирования
не
должна
превышать 3% от заданного значения командного сигнала. Во-вторых, время
переходного процесса должно быть минимальным с учетом соблюдения
предыдущего
требования.
Указанным
требованиям
удовлетворяют
коэффициенты закона (13)
K x 0.3365 23.1742 7.0711 , k 1.006 ,
полученные в результате решения задачи LQR-оптимизации с матрицами
Q diag 0 1.5 0.5 0 , R 0.01.
Далее
необходимо
найти
коэффициенты
скоростного
закона
управления. Как было сказано ранее, в частных случаях процедура поиска
может
быть
существенно
упрощена.
Обратимся
к
системе
(57),
представляющей линейную модель бокового движения СВП. По сути,
необходимо решить систему (57) относительно переменных
y ,
и .
Считая, что K x k1 k 2 k 3 , результате получим
b
b2
]/ d,
[( a12 b2 a 22 )
1
,
y
]/ d,
a11
a 21
[( a11 a 22 a 21 a12 )
(58)
d a 21b1 a11 b2 .
Подставив
(58)
в
базовый
закон
управления
коэффициенты скоростного закона управления
k 3 , 1 ( a 21 k b2 k1 ) / d , 2 (b1 k1 a11 k ) / d ,
3 ((a12 b2 b1 a 22 )k1 (a11 a 22 a12 a 21 )k ) / d k 2 .
,
(13),
получим
48
Подставляя значения K x и K в полученные формулы, получим
μ 0.2757 46.2269 29.8032 , 7.0711 .
На рис. 10 представлена динамика судна в режиме отработки
командного сигнала с полученным скоростным законом управления при
повороте на 10 градусов по курсовому углу. Можно заметить, что требование
по величине перерегулирования соблюдается, при этом переход в зону 3%
осуществляется менее чем за 10 секунд, что является довольно эффективным
результатом.
Рис. 10. Динамика курсового угла при повороте на 10º
4.3. Компенсация ступенчатых возмущений
Перейдем к синтезу следующего элемента многоцелевой структуры –
асимптотического
требованием
наблюдателя.
является
Будем
обеспечение
считать,
отклонения
что
необходимым
курсового
угла,
не
превышающего 0.4 градуса при боковом ветре 10 м/с. С учетом этого
требования получим коэффициенты матрицы G в результате решения задачи
синтеза наблюдателя как фильтра Калмана:
49
G 1.4272 1.9543 1.977 T
.
Также теперь можем перейти к позиционному закону управления по
выходу наблюдателя:
k 0.3843 23.1742 148.8717 T , k0 1.0006, 0 155.9427 .
На рис. 11 продемонстрирована динамика системы при воздействии
бокового ветра. Как видно из графика, указанное требование соблюдается,
при этом курсовой угол сходится к нулевому значению.
Рис. 11. Динамика курсового угла при воздействии бокового ветра со скоростью 10 м/с
4.4. Коррекция морского волнения
Наконец, обратимся к синтезу динамического корректора. Напомним,
что в рамках данной работы корректор выступает в роли оптимального
динамического
фильтра,
интенсивности
перекладки
целью
рулей
которого
при
является
воздействии
минимизация
периодических
возмущений (морское волнение).
Для упрощения поиска решения задачи (36) возьмем в качестве
приближения для
модели
периодических
возмущений гармоническое
50
колебание с постоянной частотой. В таком случае задача сводится к тому,
чтобы обеспечить условие
Fy ( j0 ) 0
[18]. Необходимым и достаточным
условием, которому должна удовлетворять передаточная функция
F (s )
для
этого является равенство
~
~
F( j0 ) [T12 ( j0 ) TT22 ( j0 )] 1 T,
(59)
~
где T T11 ( j0 )T21 1 ( j0 ) .
f 0 F ( j0 )
Введем обозначение
для величины, вычисленной по
формуле (59). Зададим матрицу α для системы (11) в виде
0
α
0
Передаточная функция
1
.
1
F p (s )
системы (11) является рациональной
Fp ( s )
( s )
p ( s)
дробью вида
,
где
p ( s ) det(E 2 s ),
( s ) p ( s ) γ (E 2 s α ) 1 β 1 s 0 .
Условие
Fp (i0 ) f 0
можно обеспечить, если взять
0 Re[ f 0 p (i0 )],
1
Im[ f 0 p (i0 )]
0
s.
0
Возьмем γ 0 1 , тогда столбец β можно найти как решение
1
линейной системы
1 0 0
0
,
0 0 1 1
таким образом
0 0
.
1 1 0
Зададим частоту 0 0.65 , а также корни общего знаменателя
передаточной матрицы корректора 0 1 0.2 . Проводя вычисления в
51
соответствии с подходом, описанным в соответствующей главе, получим
значения коэффициентов корректора:
0
α
0.04
1
,
0.4
69.2093
,
β
83.98
γ 0 1
.
Наличие фильтрующих свойств замкнутой системы отражается
графиком частотной характеристики от возмущения к отклонению рулей с
фильтрацией и без, показанной на рис. 12.
Рис. 12. Частотная характеристика замкнутой системы (59), (12) от возмущения к
отклонению рулей
В качестве приближенной модели морского волнения, примерно
соответствующего 4 баллам по шкале Бофорта, возьмем полигармоническое
возмущение вида
d (t ) 0.2sin(0.65t ) 0.1sin(0.55t ) 0.1sin(0.75t ).
Динамика судна для этого режима представлена на рис. 13 и рис. 14.
Здесь первые 300 секунд корректор был отключен, при этом можно
наблюдать довольно значительные колебания при перекладке рулей.
Очевидно, что включение корректора после 300-й секунды значительно
52
снижается отклонения руля, при этом также снижается интенсивность
колебаний курсового угла. Таким образом, можно сделать вывод о том, что
синтезированный
динамический
корректор,
выступающий
в
роли
оптимального фильтра, довольно эффективно справляется с поставленной
задачей, снижая интенсивность перекладки руля на порядок. При этом можно
отметить, что снижение интенсивности достигается даже с учетом того
факта, что в модели волнения присутствуют побочные гармоники кроме
основной частоты. Это обосновывается тем фактом, что значения частотной
характеристики от возмущения к отклонениям рулей на принятых частотах
существенно ниже при включенном фильтре.
Рис. 13. Динамика курсового угла при воздействии морского волнения
53
Рис. 14. Динамика перекладки рулей при воздействии морского волнения
54
Глава 5. Практическая реализация компенсирующего
регулятора
5.1. Отработка командного сигнала
Рассмотрим теперь влияние запаздывания на динамику системы,
замкнутой многоцелевым регулятором, синтезированным в предыдущей
главе. Прежде всего, обратимся к режиму отработки командного сигнала при
условии
отсутствия
внешних
возмущений
и
с
выключенным
корректирующим устройством. Динамика судна при повороте на 10 градусов
по курсовому углу и запаздывании h=0.8с представлена на рис. 5. Очевидно,
что качество управления существенно ухудшилось, появилась довольно
значительная колебательность.
При запаздывании h=1.2с система практически теряет устойчивость,
переходя в режим постоянных колебаний (рис. 16).
Рис. 15. Динамика курсового угла при повороте на 10º и запаздывании h=0.8c
55
Рис. 16. Динамика курсового угла при повороте на 10º и запаздывании h=1.2c
Таким образом, очевидно, что наличие запаздывания необходимо
учитывать при построении систем управления для СВП в частности и для
подвижных объектов в целом. Для компенсации запаздывания далее
проведем трансформацию многоцелевого регулятора в соответствии с
подходом, описанным в главе 2. Подавать управление на привод руля будем в
виде
u ku ( z ), ku 100 .
Динамика
системы,
замкнутой
многоцелевым
компенсирующим
регулятором, приведена на рис. 17. Можно отметить, что качество управления
в этом случае практически идентично качеству функционирования исходной
системы без запаздывания, что соответствует ожидаемому поведению. Таким
образом, компенсирующий регулятор справляется со своей задачей, по сути
исключая влияние запаздывания на динамику системы.
56
Рис. 17. Динамика системы, замкнутой компенсирующей обратной связью при повороте
на 10º и запаздывании h=0.8c
5.2. Компенсация ступенчатых возмущений
Далее рассмотрим влияние запаздывания на стабилизацию системы при
наличии бокового ветра со скоростью 10 м/с. Рис. 18 демонстрирует, что уже
при запаздывании h=0.5с качество отработки ступенчатого возмущения
сильно ухудшается. При этом увеличение запаздывания до 0.8с, приводит
фактически к потере устойчивости (рис. 19). Очевидно, что в этом режиме
также
необходим
учет
запаздывания
при
проектировании
системы
управления.
Применим теперь компенсирующий многоцелевой регулятор. Анализ
рис. 20 позволяет сделать вывод, что и в этом режиме применение
компенсационного подхода позволяет избавиться от негативного влияния
запаздывания и получить качество процессов, практически совпадающее с
качеством исходной системы без запаздывания.
57
Рис. 18. Динамика системы при воздействии бокового ветра со скоростью 10м/с и
запаздывании h=0.5c
Рис. 19. Динамика системы при воздействии бокового ветра со скоростью 10м/с и
запаздывании h=0.8c
58
Рис. 20. Динамика системы, замкнутой компенсирующей обратной связью при
воздействии бокового ветра со скоростью 10м/с и запаздывании h=0.5c
5.3. Коррекция морского волнения
Перейдем к рассмотрению влияния запаздывания на функционирование
системы управления в режиме коррекции волновых возмущений. Из рис. 21 и
рис. 22 видно, что запаздывание h=0.8c ощутимо увеличивает амплитуду
колебаний системы при выключенном корректоре, однако его введение в
контур управления позволяет так же эффективно снизить интенсивность
управления, а также уменьшить амплитуду колебаний курсового угла, как и в
случае системы без запаздывания. Рис. 23 и рис. 24 демонстрируют тот факт,
что использование компенсирующего регулятора позволяет значительно
улучшить динамику системы при выключенном корректоре. При этом
эффективность включенного корректора остается практически неизменной.
Таким образом, можно отметить, что компенсирующий многоцелевой
регулятор так же может быть использован и в режиме морского волнения с
запаздыванием.
59
Рис. 21. Динамика курсового угла при воздействии волнения и запаздывании h=0.8c
Рис. 22. Динамика перекладки рулей при воздействии волнения и запаздывании h=0.8c
60
Рис. 23. Динамика курсового угла при использовании компенсирующего регулятора при
воздействии волнения и запаздывании h=0.8c
Рис. 24. Динамика перекладки рулей при использовании компенсирующего регулятора при
воздействии волнения и запаздывании h=0.8c
61
Выводы
Прежде всего, полученные результаты исследования позволяют сделать
вывод о том, что применение многоцелевых регуляторов для стабилизации
подвижных объектов по линейным моделям позволяет эффективно решать
поставленные задачи с учетом всего комплекса требований, предъявляемых к
динамике системы управления в различных режимах функционирования. При
этом относительная независимость отдельных элементов многоцелевой
структуры упрощает процесс синтеза, позволяя разбить исходную задачу на
отдельные более простые подзадачи. Также явным преимуществом является
возможность отключения отдельных элементов многоцелевой структуры в
зависимости от текущей ситуации.
Далее можно отметить существенное влияние запаздывания на динамику
системы, замкнутой многоцелевым регулятором, в сторону значительного
ухудшения качества управления. В частности, для движения СВП с
полученными коэффициентами регулятора запаздывания величиной порядка 1с
достаточно для полной потери устойчивости. Таким образом, очевиден тот
факт, что влияние запаздывания необходимо учитывать при проектировании
системы управления.
Наконец, можно отметить тот факт, что компенсировать запаздывание
позволяет особая трансформация многоцелевого регулятора. Важнейшим
моментом здесь является сохранения динамических свойств исходной системы
без
запаздывания.
Представленные
графики
наглядно
демонстрируют
практически полную идентичность динамических процессов. Кроме этого, в
качестве
существенного
преимущества
можно
выделить
отсутствие
необходимости решать отдельную задачу синтеза. Компенсация достигается за
счет
добавления
вспомогательной
векторной
функции
с
собственной
динамикой, а также трансформации входа регулятора с учетом прогноза.
62
Если говорить о недостатках использованного подхода, то здесь можно
отметить необходимость измерения возмущающего воздействия, что не всегда
возможно сделать точно. Также открытым остается вопрос робастности
полученного управления.
Однако в целом можно сделать итоговый вывод о том, что
представленные
подходы
действительно
обеспечивают
эффективную
стабилизацию подвижных объектов даже в условиях развитого морского
волнения, воздействия ветра, а также при наличии запаздывания в контуре
управления.
63
Заключение
В
работе
были
рассмотрены
задачи
управления
движением
произвольного объекта с учетом комплекса ограничений и требований к
реакции на внешние постоянные и периодические возмущения, а также учет
запаздывания в это задаче. Полученные результаты позволяют сделать вывод о
том, что многоцелевая структура является мощным инструментом для
достижения поставленных целей, при этом оставаясь удобной для реализации
на бортовых устройствах в режиме реального времени и гибкой в тех случаях,
когда учет всех требований не является обязательным условием работы.
Кроме
этого
было
продемонстрировано
негативное
влияние
запаздывания на стабилизацию подвижных объектов на примере задачи
стабилизации СВП при движении над морской поверхностью. Результаты,
полученные при учете запаздывания в задаче стабилизации, позволяют
говорить о явных преимуществах компенсационного подхода. При наличии
исходного базового регулятора для системы без запаздывания компенсация
позволяет не только минимизировать эффект запаздывания, но и получить те
же динамические характеристики, т.е. сохранить желаемое расположение
собственных чисел, а также передаточную матрицу системы. Кроме этого,
компенсирующий регулятор достаточно просто и эффективно реализуется на
бортовом оборудовании и пригоден для работы в режиме реального времени.
В качестве основного дальнейшего направления развития можно
выделить исследование робастных свойств использованных регуляторов, а
также сравнение компенсационного подхода с методами типа непрерывного
размещения собственных чисел (continuous pole placement [13]).
64
Список литературы
1.
Fossen, T.I., Guidance and Control of Ocean Vehicles, New York, John Wiley
& Sons, 1994.
2.
Perez, T., Ship Motion Control: Course Keeping and Roll Stabilization Using
Rudder and Fins, London: Springer–Verlag, 2005.
3.
Wang Xianzhou and Xu Hanzhen, Robust Autopilot with Wave Filter for Ship
Steering, Journal of Marine
4.
Веремей Е.И., Корчанов В.М. Многоцелевая стабилизация динамических
систем одного класса. Автоматика и телемеханика, 1988, №9. С. 126–137.
5.
Веремей Е.И., Еремеев В.В., Корчанов В.М. Синтез алгоритмов
робастного
управления
движением
подводных
лодок
вблизи
взволнованной поверхности моря. Гироскопия и навигация, 2000, № 2.
С. 34–43.
6.
Веремей
Е.И.,
Корчанов
В.М.,
Коровкин
М.В.,
Погожев
С.В.
Компьютерное моделирование систем управления движением морских
подвижных объектов. СПб.: НИИ Химии СПбГУ, 2002. 370 С.
7.
E.I. Veremey. Dynamical Correction of Control Laws for Marine Ships'
Accurate Steering // Journal of Marine Science and Application. – Volume 13,
Issue 2. – June 2014. – Pages 127-133. (DOI: 10.1007/s11804-014-1250-1).
8.
Veremey E.I. Optimization of filtering correctors for autopilot control laws with
special structures // Optimal Control Applications and Methods. Volume 37,
Issue 2, pages 323 – 339. March/April 2016. (DOI: 10.1002/oca.2170).
9.
Красовский Н.Н. Некоторые задачи теории устойчивости движения. М.:
Физматгиз, 1959. 211 с.
10. Янушевский Р.Т. Управление объектами с запаздыванием. М.: Наука, 1978.
416 с.
11. Жабко А. П., Харитонов В. Л. Методы линейной алгебры в задачах
управления. СПб.: Изд-во С.–Петерб. ун-та, 1993. 320 с.
65
12. Kharitonov V.L. Time-Delay Systems. Lyapunov Functionals and Martices.
Birkhäuser, 2013. 312 p.
13. W. Michiels and S.-I. Niculescu. Stability and Stabilization of Time-Delay
Systems. An Eigenvalue-Based Approach. SIAM, Philadelphia, 2007, 378 p.
14. Smith O.J.M. A Controller to overcome dead time // ISAJ. 1959. Vol. 6, No. 2.
P. 28–33.
15. Marshall J.E., Gorecki H., Korytowski A., Walton K. Time-Delay Systems.
Stability and Performance Criteria with Applications. Ellis Horwood, 1992. 244
p.
16. Krstic M. Delay Compensation for Nonlinear, Adaptive, and PDE Systems.
Birkhäuser, 2009. 466 p.
17. Веремей Е.И. Компенсирующие регуляторы по выходу для LTI систем с
запаздыванием по управлению. Современные методы прикладной
математики,
теории
управления
и
информационных
технологий
(ПМТУКТ-2017) Сборник трудов X международной конференции. 2017.
С. 106 – 110.
18. Веремей Е.И. Среднеквадратичная многоцелевая оптимизация: учеб.
пособие. СПб.: Изд-во С.-Петерб. ун-та, 2016. 408 с.
19. Лукомский Ю.А., Чугунов В.С. Системы управления морскими
подвижными объектами. Л.: Судостроение, 1988. 272 с.
20. Войткунский Я.И. Справочник по теории корабля: в трех томах. Том 3.
Управляемость
водоизмещающих
судов.
Гидродинамика
судов
с
динамическими принципами поддержания. Л.: Судостроение, 1985. 544 с.
21. Yun L., Bliault A. Theory and Design of Air Cushion Craft. London: Arnolds,
2000. 632 p.
22. Колызаев Б.А., Косоруков А.И., Литвиненко В.А. Справочник по
проектированию судов с динамическими принципами поддержания. Л.:
Судостроение, 1980. 472 с.
Перечень собственных публикаций по теме
66
1.
Sevostyanov R., Shayakhmetova L. Stabilization of the inverted pendulum
considering delay // Selected Papers of the II International Scientific
Conference "Convergent Cognitive Information Technologies" (Convergent
2017). Moscow, Russia, 2018. P. 318 – 324. (Scopus)
2.
Севостьянов Р.А. Цифровая стабилизация движения колесного робота //
Системы управления и информационные технологии. Т. 71, Н. 1.
Воронеж, 2018. С. 43 – 36. (ВАК)
3.
Севостьянов Р.А. Автоматизация исследований динамических систем //
Современные методы прикладной математики, теории управления и
информационных технологий (ПМТУКТ-2017) Сборник трудов X
международной конференции. Воронеж: Новая книга. 2017. С. 318 – 322.
(РИНЦ)
4.
Sevostyanov Ruslan A.
Air Cushion Vehicle Model for Educational
Purposes // Proceedings of the 2nd International Conference on Applications
in Information Technology. Aizu, Japan, 2016. P. 34 – 37.
5.
Sevostyanov R., Shayakhmetova L. Automation system for the research
process of the air cushion vehicle model // Selected Papers of the I
International Scientific Conference "Convergent Cognitive Information
Technologies" (Convergent 2017). Moscow, Russia, 2018. P. 127 – 134.
(Scopus)
6.
Sevostyanov Ruslan A., Veremey Evgeny I. Multipurpose stabilization of the
advanced marine surface crafts // 6th Seminar on Industrial Control Systems Analysis, Modeling and Computation, ITM Web of Conferences, Vol. 6, 2016.
(WoS CC)
7.
Веремей Е.И., Погожев С.В., Севостьянов Р.А. Стабилизация курса
морских судов в экономичном режиме движения // Сборник трудов VIII
международной
конференции
«Современные
методы
прикладной
математики, теории управления и компьютерных технологий», Воронеж:
Новая Книга, 2015. С. 92 – 95. (РИНЦ)
67
8.
Веремей Е.И., Погожев С.В., Севостьянов Р.А. Фильтрующая коррекция
многоцелевых законов управления движением морских судов // Системы
управления и информационные технологии, 2015, № 4 (62). С. 4 – 7.
(ВАК)
9.
Севостьянов Р.А.. Стабилизация движения робота при наличии
транспортного запаздывания // Навигация и управление движением:
материалы XVI конференции молодых ученых. СПб.: Электроприбор,
2014. С. 258 – 264. (РИНЦ)
10. Севостьянов Р.А. Управление движением мобильного робота с учетом
транспортного запаздывания // Процессы управления и устойчивость:
Труды 45-й международной научной конференции аспирантов и
студентов. СПб.: Издат. Дом С.-Петерб. гос. ун-та, 2014. С. 385 – 390.
(РИНЦ)
11. Севостьянов Р.А., Нефедов Д.Э. Система удаленного управления
мобильным роботом // Процессы управления и устойчивость: Труды 44-й
международной научной конференции аспирантов и студентов. СПб.:
Издат. Дом С.-Петерб. гос. ун-та, 2013. С. 479 – 484.
12. Севостьянов Р.А., Нефедов Д.Э., Скрябина Т.Д. Программная поддержка
процессов управления мобильным роботом // Навигация и управление
движением: материалы XV конференции молодых ученых. СПб.:
Электроприбор, 2013. С. 224 – 229. (РИНЦ)
13. Севостьянов Р.А. Программная поддержка процессов управления
гусеничным роботом в реальном времени // Процессы управления и
устойчивость: Труды 43-й международной научной конференции
аспирантов и студентов. СПб.: Издат. Дом С.-Петерб. гос. ун-та, 2012. C.
388 – 393.
Отзывы:
Авторизуйтесь, чтобы оставить отзыв