Full Text
О семинаре по проблемам нелинейной динамики и управления при Московском государственном университете имени М.В. Ломоносова1
А. В. Ильин, А. С. Фурсов, Ю. М. Мосолова (МГУ ВМК, Москва, Россия) О задаче стабилизации переключаемой интервальной линейной системы с соизмеримыми запаздываниями 18.12.2023
В настоящей работе исследуется задача цифровой стабилизации переключаемой интервальной линейной системы в случае, когда её режимы имеют различные запаздывания в управлении, а именно, рассматривается непрерывная скалярная переключаемая интервальная линейная система с соизмеримыми запаздываниями в управлении
(1)
где полунепрерывная справа кусочно-постоянная функция (ненаблюдаемый переключающий сигнал); множество переключающих сигналов , точки разрыва которых принадлежат множеству , некоторое положительное число, а ; вектор состояния; измеряемый скалярный выход; управляющий вход; композиция отображения и переключающего сигнала ; , и аналогичные композиции для отображений и , ( , , ( ) интервальные матрицы соответствующих размеров). Здесь величины постоянных запаздываний, причём рациональное число для любой пары .
Значение функции в каждый момент времени определяет активный режим переключаемой системы (1), описываемый интервальной линейной стационарной системой с запаздыванием в управлении
Решением уравнения состояния системы (1) при фиксированных тройках ( , , , ), заданном управлении (считаем, что при ), переключающем сигнале и начальном условии является решение соответствующей линейной нестационарной системы с запаздыванием
Требуется для переключаемой линейной системы (1) с заданным числом и ненаблюдаемыми переключающими сигналами построить цифровой регулятор вида
(2)
(3)
обеспечивающий глобальную равномерную асимптотическую устойчивость замкнутой непрерывно-дискретной системы
(4)
где
Система (4) записана при условии, что моменты времени и согласованы, , т.е.
Здесь период квантования по времени (считаем, что , существует такое, что , и для любого найдётся такое , что ), целая часть действительного числа, , , , ( порядок регулятора), , , дискретные функции, определённые на последовательности , формирующий элемент представлен фиксатором нулевого порядка [1, c. 25].
Замкнутую непрерывно-дискретную систему (4) называем глобально равномерно асимптотически устойчивой, а регулятор (2), (3) стабилизирующим, если при любых фиксированных тройках ( , , , ) для любых , и для соответствующего решения выполнено:
Для решения поставленной задачи предлагается применять подходы, разработанные в статьях [24], где исследовались проблемы цифровой стабилизации переключаемых интервальных линейных систем и переключаемых линейных систем с запаздыванием в управлении, одинаковом для всех режимов рассматриваемой переключаемой системы. Для решения первой проблемы в [3, 4] используется метод дискретизации переключаемой интервальной системы с дальнейшим поиском динамического дискретного регулятора на основе линейных матричных неравенств. В работе [2] предлагается подход для решения второй проблемы, включающий два основных шага переход от исходной непрерывной системы к её точной дискретной модели (вообще говоря, более высокого динамического порядка) и поиск дискретного динамического регулятора для полученной переключаемой дискретной системы. Модифицируя указанные подходы с учётом наличия соизмеримых запаздываний в режимах исходной переключаемой системы (1), разработана общая схема построения цифрового регулятора (2), (3), включающая следующие основные шаги:
1) переход от непрерывной системы (1) к её точной дискретной модели с учётом, что на её входе используется фиксатор нулевого порядка (точная дискретная модель фактически представляет собой семейство переключаемых дискретных систем с режимами, описываемыми системами разностных уравнений, вообще говоря, различных динамических порядков, не содержащих запаздываний);
2) построение интервального расширения для полученной дискретной модели;
3) приведение режимов полученной переключаемой линейной интервальной дискретной системы к единому порядку на основе метода расширения динамического порядка [5, c. 205];
4) построение процедуры численного поиска стабилизирующего дискретного регулятора вида (3) с использованием достаточного условия устойчивости переключаемых интервальных дискретных систем на основе метода функций Ляпунова [3].
Работа выполнена при финансовой поддержке Министерства науки и высшего образования Российской Федерации в рамках реализации программы Московского центра фундаментальной и прикладной математики по соглашению 075-15-2019-1621.
Литература. 1. Поляков, К.Ю. Основы теории цифровых систем управления: учеб. пособие / К.Ю. Поляков. СПб. : СПбГМТУ, 2002. 154 с. 2. Фурсов, А.С. Построение цифрового стабилизатора для переключаемой линейной системы с запаздыванием в управлении / А.С. Фурсов, С.И. Миняев, В.С. Гусева // Дифференц. уравнения. 2018. Т. 54, 8. С. 11321141. 3. Фурсов, А.С. Синтез цифрового стабилизатора по выходу для переключаемой интервальной линейной системы / А.С. Фурсов, С.И. Миняев, Ю.М. Мосолова // Дифференц. уравнения. 2019. Т. 55, 11. С. 15451559. 4. Фурсов, А.С. Цифровая сверхстабилизация переключаемой интервальной линейной системы / А.С. Фурсов, Ю.М. Мосолова, С.И. Миняев // Дифференц. уравнения. 2020. Т. 56, 11. С. 15161527. 5. Фурсов, А.С. Одновременная стабилизация: теория построения универсального регулятора для семейства динамических объектов / А.С. Фурсов. М. : Аргамак-Медиа, 2016. 238 с.
В. Е. Хартовский (ФИТМ ГрГУ имени Я. Купалы, Гродно, Беларусь) К вопросу назначения конечного спектра линейной системе нейтрального типа 16.10.2023
Объект исследования линейная автономная система нейтрального типа
(5)
(6)
Здесь решение уравнения (5); управление, выбираемое из класса кусочно-непрерывных функций; наблюдаемый выход; оператор сдвига, определяемый для заданного правилом , ;
где , , , .
Пусть единичная матрица, характеристическая матрица уравнения (5). Спектр (множество собственных значений) разомкнутой ( ) системы (5) совпадает с множеством корней уравнения , которое в подробной записи имеет вид , где некоторые полиномы, . Из этого уравнения видно, что спектр системы (5) состоит из бесконечного числа чисел.
Задача назначения конечного спектра заключается в построении регулятора с обратной связью такого, чтобы замкнутая система имела наперёд заданный конечный набор собственных значений, состоящий, как правило, из произвольных чисел с отрицательными действительными частями. Такая задача исследовалась многими авторами, наиболее важные результаты в этом направлении, а также решение более общей задачи модального управления (управления при помощи обратной связи коэффициентами характеристического уравнения) обсуждаются в работах [1, 2] (см. Введение). В статье [3] для систем запаздывающего типа с одним входом предложен алгебраический подход, позволяющий с единой позиции подойти к решению задач назначения конечного спектра и полной управляемости. При исследовании проблемы управления коэффициентами характеристического уравнения системы нейтрального типа (5) установлено [4], что в случае выполнения условия модальной управляемости исходную задачу можно свести к аналогичной задаче для системы запаздывающего типа с одним входом и дальнейшее исследование провести методами, описанными в [3]. Там же показано, что задача назначения конечного спектра разрешима тогда и только тогда, когда разрешима задача модальной управляемости. Систематизация подхода [4], а также исследование ослабленного варианта задачи модальной управляемости выполнены в статье [5, c. 321]. В настоящем сообщении методы [4, 5] обобщены на задачу назначения конечного спектра. Однако, в отличие от работ [4, 5], регулятор построен по неполным измерениям, определяемым наблюдаемым выходом (6). Основная идея заключена в модификации структуры регуляторов [4; 5, c. 321] посредством встраивания в их внутренний контур системы точной оценки решения, основанной на использовании конструкций финитных наблюдателей, разработанных в [5, 6]. Перейдём к постановке задачи.
Определим регулятор с обратной связью по выходу следующего вида:
(7)
Здесь вспомогательная переменная, операторы определяются на множестве непрерывных функций , имеющих кусочно-непрерывную производную, по правилу
(8)
где , заданные полиномиальные матрицы, заданные кусочно-непрерывные матричные функции, . Каждому оператору вида (8) поставим в соответствие матрицу
Выпишем характеристическую матрицу замкнутой системы (5)(7):
(9)
Задача 1. Пусть задан произвольный полином с вещественными коэффициентами. Требуется построить регулятор вида (7) такой, что определитель характеристической матрицы (9) замкнутой системы (5)(7) совпадает с полиномом , т.е. .
Получен критерий существования регулятора (7), обеспечивающего решение задачи 1.
Теорема 1. Задача 1 разрешима для любого полинома , , имеющего по крайней мере два вещественных корня (возможно, равных между собой), тогда и только тогда, когда выполняются условия:
1) , ;
2) , ;
3) , ;
4) , .
Условия 1), 2) представляют собой критерий модальной управляемости системы (5). Условия 3), 4) необходимы и достаточны для существования финитного наблюдателя для системы (5), (6). Эти же условия определяют критерий финальной наблюдаемости системы (5), (6) условия существования непрерывной операции восстановления “отрезка решения” , , при некотором достаточно большом , по результатам наблюдения , (при известном управлении , ).
Литература. 1. Хартовский, В.Е. Управление линейными системами нейтрального типа: качественный анализ и реализация обратных связей / В.Е. Хартовский. Гродно : ГрГУ, 2022. 500 с. 2. Зайцев, В.А. Модальное управление и стабилизация линейных систем статической обратной связью по выходу / В.А. Зайцев, И.Г. Ким. Ижевск : Изд. центр «Удмуртский университет», 2022. 184 с. 3. Метельский, А.В. Алгебраический подход к стабилизации дифференциальной системы запаздывающего типа / А.В. Метельский // Дифференц. уравнения. 2018. Т. 54, 8. С. 11191131. 4. Метельский, А.В. Критерии модальной управляемости линейных систем нейтрального типа / А.В. Метельский, В.Е. Хартовский // Дифференц. уравнения. 2016. Т. 52, 11. С. 15061521. 5. Метельский, А.В. Синтез финитного наблюдателя для линейных систем нейтрального типа / А.В. Метельский, В.Е. Хартовский // Автоматика и телемеханика. 2019, 12. С. 80102. 6. Метельский, А.В. О точном восстановлении решения линейных систем нейтрального типа / А.В. Метельский, В.Е. Хартовский // Дифференц. уравнения. 2021. Т. 57, 2. С. 265285.
Е. И. Атамась (МГУ ВМК, Москва, Россия) Один метод реализации интервальных систем 20.11.2023
Задача реализации управляемой линейной системы является классической и для стационарных систем давно имеет полное решение. Построение реализации тривиально в случае скалярных систем и немногим более сложно для векторных систем. В данном сообщении мы зададимся вопросом о методах построения реализации для интервальных передаточных матриц, т.е. матриц, элементами которых являются рациональные выражения с интервальными коэффициентами [1], например, в случае матриц размера будем иметь
где , вещественные интервалы.
Из-за привлечения интервальных вычислений непосредственный перенос многих методов построения реализации, известных для систем с точечными параметрами, оказывается невозможен. В работе [2] был предложен подход для скалярных систем, основанный на использовании канонической формы управляемости/наблюдаемости. К сожалению, обобщить его на векторные системы крайне затруднительно: построение известных аналогов таких канонических форм сопряжено с заменами координат, что в интервальном случае практически невозможно.
Попробуем применить для решения поставленной задачи алгоритм реализации Гилберта. Суть его состоит в разложении передаточной матрицы в сумму элементарных дробей, после чего каждая из них реализуется отдельно тривиальным образом. Однако в интервальной арифметике классический способ разложения в элементарные дроби на основе метода неопределённых коэффициентов также оказывается малоприменим, поскольку приводит к решению интервальной системы линейных алгебраических уравнений. Для получения более простого подхода сделаем следующие предположения о передаточной функции . Будем считать, что физически реализуема, т.е. степень числителя каждого её элемента строго меньше степени знаменателя, а все полюсы системы простые и вещественные. В этом случае можно записать равенство
где интервальные корни многочлена Сделанное нами предположение о простоте корней позволяет использовать относительно простые алгоритмы для их поиска (см., например, [3]), а также вычислять коэффициенты , для этого применим хорошо известную формулу где . Тогда нам потребуется лишь вычислить значения интервальных полиномов в интервальных точках, что является простой операцией. Так, каждой дроби вида будет соответствовать часть реализации
где единичная матрица соответствующего размера.
Таким образом, предложенный алгоритм позволяет строить реализацию векторной интервальной системы при указанных ограничениях.
Литература. 1. Прикладной интервальный анализ / Л. Жолен, М. Кифер, О. Дидри, Э. Вальтэр; пер. с англ. С.И. Куликова. М.; Ижевск: Институт компьютерных исследований, 2015. 467 с. 2. Атамась, Е.И. О переходе между различными представлениями интервальных управляемых систем / Е.И. Атамась // Вестн. Московского ун-та. Серия 15: Вычислит. математика и кибернетика. 2023. 4. С. 38. 3. Hansen, E.R. Sharp bounds on interval polynomial roots / E.R. Hansen // Reliable Computing. 2002. 8. P. 115122.
А. И. Астровский (БГЭУ, Минск, Беларусь) Наблюдаемость линейных нестационарных систем с ограничением на выходную функцию 11.12.2023
В классической постановке задачи наблюдаемости [1] для линейных систем предполагалось, что выходная функция системы в момент линейно связана с состоянием системы в этот же момент и её измерения всегда доступны наблюдателю. Вместе с тем в приложениях (например, в теории электрических цепей, в навигационной теории, в медицинских исследованиях и т.д.) есть задачи наблюдения для линейных систем, у которых выходная функция может быть измерена тогда и только тогда, когда фазовый вектор принадлежит некоторому заданному множеству из пространства Задачи наблюдения такого типа относят к специальному нелинейному классу задач наблюдения и называют задачами наблюдения с ограничениями на выходную функцию. Задачи наблюдения с геометрически ограниченными наблюдениями исследовались, например, в [24] для стационарных систем при ограничениях на фазовый вектор в виде конусов в конечномерном пространстве.
Рассмотрим для линейных нестационарных систем обыкновенных дифференциальных уравнений задачи наблюдаемости при условии, что выходная функция системы может быть измерена тогда и только тогда, когда она положительна:
(10)
Здесь -вектор-столбец состояния системы (10) в момент ; -матричная функция и -вектор-строка непрерывны на ; для каждого .
Говорят, что начальные состояния и из системы (10) положительно различимы, если соответствующие им выходные функции и не совпадают тождественно на отрезке , т.е. , Система (10) положительно наблюдаема, если любые два начальные состояния из положительно различимы.
Свойство положительной наблюдаемости накладывает довольно жёсткие требования на систему наблюдения. Например, при любая стационарная система вида (10) не является положительно наблюдаемой (хотя она наблюдаема при ). Скалярная нестационарная система наблюдения со знакопеременной функцией , , служит примером положительно наблюдаемой системы.
Отождествим систему , , , с парой . Пусть , , являются компонентами -вектор-строки , где фундаментальная матрица системы (10): , ,
Система (10) задаёт отображение “начальное состояние” “выходная функция”, которое действует по правилу , . Понятно, что система (10) положительно наблюдаема тогда и только тогда, когда отображение инъективно. Очевидно, что наблюдаемость пары является необходимым и достаточным условием положительной наблюдаемости системы (10).
Хорошо известно [1], что необходимым и достаточным условием наблюдаемости пары является линейная независимость системы функций на
Говорят, что система функций вполне линейно независима на отрезке , если любая нетривиальная линейная комбинация этих функций может быть равна нулю только на множестве меры ноль. Точку называют корнем-узлом для непрерывной функции , если и при переходе через эту точку значения функции меняют знак.
Система непрерывных функций на отрезке называется [5, 6] системой функций Чебышёва порядка , если любая нетривиальная линейная комбинация этих функций имеет не более корней на
Непосредственно из определения положительной наблюдаемости следует
Лемма 1. Для положительной наблюдаемости системы (10) необходимо, чтобы каждая нетривиальная линейная комбинация по системе функций имела хотя бы один корень-узел на промежутке .
Из определения положительной наблюдаемости и леммы 1 вытекает
Теорема 1. Система (10) положительно наблюдаема на отрезке , если для пары функции , ..., вполне линейно независимы на и у каждого нетривиального многочлена по этим функциям существует хотя бы один корень-узел на интервале .
Рассмотрим на отрезке систему наблюдения второго порядка: , , . Для этой системы функции , вполне линейно независимы на отрезке . Любой нетривиальный многочлен по этим функциям имеет точку в качестве корня-узла. Следовательно, данная система в силу теоремы 0 положительно наблюдаема на
Покажем, что ни одно из условий теоремы 0 нельзя опустить. Как следует из леммы 0, наличие хотя бы одного корня-узла у каждого многочлена , является необходимым условием положительной наблюдаемости системы (10). Ниже приведём пример системы (10), для которой функции , ..., не вполне линейно независимы на , у каждого нетривиального многочлена существует по крайней мере один корень-узел на и система (10) не является положительно наблюдаемой.
Пусть и , , , где непрерывные функции , имеют вид
Несложно заметить, что в приведённом примере и система функций линейно независима, но не вполне линейно независима на Функции можно представить в виде
Легко проверяется, что каждый нетривиальный многочлен имеет корень-узел. Для начального момента времени начальные состояния и порождают выходные функции
Понятно, что , , т.е. начальные состояния порождают одну и ту же выходную функцию. Поэтому данная система положительно ненаблюдаема.
Простые примеры показывают, что только одно условие вполне линейной независимости системы функций на отрезке не гарантирует положительную наблюдаемость системы (10). Следовательно, положительная наблюдаемость и дифференциальная наблюдаемость [57] существенно различные свойства систем наблюдения и одно из другого не следуют.
Лемма 2. Для того чтобы у каждого нетривиального многочлена по системе функций существовал хотя бы один корень-узел на промежутке , необходимо и достаточно, чтобы одновременно выполнялись два неравенства:
Доказательство леммы 2 следует из свойств непрерывных функций.
Теорема 2. Если функции , ..., образуют систему функций Чебышёва порядка на отрезке , то система (10) не является положительно наблюдаемой.
Доказательство теоремы 0 следует из того факта, что в линейной оболочке всякой системы функций Чебышёва существуют [8, 9] как строго положительные, так и строго отрицательные многочлены на .
Исходя из теоремы 2, несложно привести примеры систем наблюдения, которые не являются положительно наблюдаемыми. Например, система (10) вида , , , наблюдаема на любом отрезке , , но она не является положительно наблюдаемой по выходу на отрезке , так как соответствующие этой системе наблюдения функции образуют систему функций Чебышёва порядка на .
Пусть множество всех невырожденных при каждом квадратных -матриц с непрерывно дифференцируемыми элементами. Действие группы на паре определим стандартным образом: ,
Непосредственно из леммы 2.3 монографии [7] следует, что свойство положительной наблюдаемости инвариантно относительно действия группы на множестве систем наблюдения с непрерывными элементами.
Рассмотрим задачу положительной наблюдаемости для равномерно наблюдаемых [7] систем . Всякая равномерно наблюдаемая пара является [7] дифференциально наблюдаемой, и поэтому система функций вполне линейно независима на . Далее нам понадобится отображение, определённое по правилу
которое, как показано в [7], является полным инвариантом относительно действия группы на множестве равномерно наблюдаемых систем. Обозначим через , , компоненты -вектор-функции
Теорема 3. Пусть пара равномерно наблюдаема на отрезке . Система (1) положительно наблюдаема на тогда и только тогда, когда каждое нетривиальное решение скалярного дифференциального уравнения имеет хотя бы один корень-узел на интервале
Доказательство теоремы 3 основано на возможности преобразования равномерно наблюдаемой пары с помощью преобразования , , к системе вида (10) с матрицей в форме Фробениуса [7] и вектором .
Теорема 4. Пусть пара равномерно наблюдаема на отрезке и её полный инвариант , . Система (10) положительно наблюдаема тогда и только тогда, когда уравнение не имеет действительных корней.
Доказательство теоремы 4 следует из соотношения .
Оказывается, что задача положительной наблюдаемости тесно связана с проблемой неосцилляции [1].
Теорема 5. Если пара равномерно наблюдаема на отрезке и дифференциальное уравнение неосцилляционно, то система (1) не является положительно наблюдаемой системой.
Доказательство теоремы 5 следует из теорем 0, 0 и свойств неосцилляции [10].
Литература. 1. Красовский, Н.H. Теория управления движением / Н.H. Красовский. М. : Наука, 1968. 476 с. 2. Brammer, R.F. Geometrically constrained observability / R.F. Brammer // SIAM J. Control. 1974. V. 12, 3. P. 449459. 3. Sontag, E.D. Mathematical Control Theory / E.D. Sontag. Berlin : Springer-Verlag, 1990. 4. Астровский, А.И. Положительная наблюдаемость линейных нестационарных систем / А.И. Астровский // Изв. НАН Беларуси. Сер. физ.-мат. наук. 1999. 2. С. 3339. 5. Weiss, L. The concepts of differential controllability and differential observability / L. Weiss // J. Math. Anal. and Appl. 1965. 10. P. 442449. 6. Silverman, L.M. Controllability and observability in time-variable linear systems / L.M. Silverman , H.E. Meadows // SIAM J. Control. 1967. V. 5, 1. P. 6473. 7. Астровский, А.И. Линейные системы с квазидифференцируемыми коэффициентами: управляемость и наблюдаемость движений / А.И. Астровский, И.В. Гайшун. Минск : Беларус. навука, 2013. 214 с. 8. Карлин, С. Чебышевские системы и их применение в анализе и статистике / С. Карлин, В. Стадден; пер. с англ. под ред. С.М. Ермолова. М. : Наука, 1976. 568 с. 9. Крейн, М.Г. Проблема моментов Маркова и экстремальные задачи / М.Г. Крейн, А.А. Нудельман. М. : Наука, 1973. 551 с. 10. Левин, А.Ю. Неосцилляция решений уравнения / А.Ю. Левин // Успехи мат. наук. 1969. Т. 24. Вып. 2 (146). С. 4396.
1 Семинар основан академиками РАН С.В. Емельяновым и С.К. Коровиным.