СИНТЕЗ РОБАСТНОГО УПРАВЛЕНИЯ ИНТЕРВАЛЬНЫМИ СИСТЕМАМИ С ЗАПАЗДЫВАНИЕМ НА ОСНОВЕ ГЕТЕРОГЕННОЙ ДЕКОМПОЗИЦИИ ФУНКЦИОНАЛА ЛЯПУНОВА — КРАСОВСКОГО

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

Огромную роль в становлении теории систем с запаздыванием сыграли работы

Классические и структурные методы активно развиваются на базе частотных методов синтеза регуляторов, которые позволили адаптировать подходы на основе ПИД-регуляторов

Современный этап развития характеризуется переходом к робастным и адаптивным подходам

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

Работы Красовского и Колмановского заложили базовую форму функционала, предлагаемая статья является естественным развитием на случай априорно известной неравномерной чувствительности. Равномерное разбиение — частный случай гетерогенного, но в случае гетерогенного — «энергия» функционала Ляпунова — Красовского распределяется по оси запаздывания неравномерно. Приведен численный пример для системы, чувствительной к малому запаздыванию. Равномерное разбиение «размазывает» штраф, а гетерогенное – концентрирует его в нужной области, что снижает консервативность системы.

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

Рассматривается линейная динамическая система с интервальной параметрической неопределенностью и запаздыванием по состоянию

где x(t)∈Rn — вектор фазовых координат; u(t)∈Rm — управляющее воздействие; τ(t) —переменное запаздывание с ограниченной производной, удовлетворяющее условиям

Интервальные матрицы представимы в центрированной форме:

аналогичные представления справедливы и для матриц AdI и BdI.

Для более компактного описания коррелированных неопределенностей может использоваться следующее представление в виде политопа:

где Gi — постоянные матрицы, определяющие направления возмущений по параметрам.

Цель синтеза состоит в определении параметров обратной связи по полному вектору состояния (введение наблюдателя не требуется; в общем случае, коэффициент обратной связи может быть интервальной матрицей)

такой, чтобы замкнутая система

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

Проведем гетерогенную декомпозицию интервала запаздывания на отрезки

где точки разбиения располагаться неравномерно, адаптируясь к динамике системы (1). Такая параметризация позволяет снизить консервативность условий устойчивости в силу того, что функционал Ляпунова — Красовского точнее отражает расположение «энергии» системы по оси запаздывания.

Пусть функционал Ляпунова-Красовского имеет следующую структуру

где

Здесь P>0, Qj>0, Rj>0, Zj>0 — матрицы, подлежащие определению. Слагаемое (4) — стандартная квадратичная форма, как в системах без запаздывания. Слагаемое (5) — штраф за отклонение состояния на интервале запаздывания. Слагаемое (6) — штраф за отклонение управляющего воздействия на интервале запаздывания. Слагаемое (7) — учитывает «энергию» производной состояния на интервале запаздывания.

Интегральные неравенства, возникающие при вычислении полной производной функционала Ляпунова-Красовского, оцениваются неравенством Йенсена, которое дает грубую оценку. Для любого интервала [a,b] и любой положительно определенной матрицы справедливо неравенство

где Ωk — векторы, определяемые через ортогональные полиномы Лежандра Pk(s) на [a,b] :

При этом неравенство (4) соответствует неравенству Йенсена, при k=1 — неравенству Виртингера, при k≥2 — уточненные оценки высших порядков. В настоящей работе использовано m=3.

После подстановки всех оценок и выполнения преобразований производная функционала Ляпунова-Красовского будет ограничиваться сверху квадратичной формой

Здесь ξT(t)=(xT(t), xT(t-τ(t)) xT(t-τm) xT(t-τM) x(t)) - расширенный вектор размерности R5n; — размерность одного блока расширенного вектора; Φ∈R5n×5n — симметричная блочная матрица, имеющая 5 блоков Φij∈Rn×n, определенных ниже:

Здесь символ обозначает симметричные относительно диагонали блоки.

где: P>0, Q>0, R>0, Z>0 — положительно определенные матрицы из функционала Ляпунова-Красовского; X∈Rn×n — дополнительная матрица; τM — верхняя граница запаздывания; μ — ограничение на производную запаздывания.

Таким образом, матрица в формуле (9) является симметричной блочной матрицей размера 5nx5n с элементами, содержащими произведения неизвестных матриц, которые линеаризуются с помощью замены переменных X=P-1, Y=KX, с использованием преобразования конгруэнтности (умножение слева и справа на блочно-диагональную матрицу "blok" "diag" {X, X, X, X, X} и её транспонирование). Последнее позволяет привести матрицу привести к окончательным линейным матричным неравенствам.

где

При этом матрица Q в функционале Ляпунова-Красовского заменяется на матрицу

После решения линейных матричных уравнений искомые матрицы вычисляются как

Вычисление коэффициента обратной связи по формуле (10) обеспечивает асимптотическую устойчивость замкнутой системы (2) при всех допустимых параметрах из заданных интервалов.

Научная новизна предлагаемого подхода состоит в следующем — гетерогенная декомпозиция + интервальная робастная обратная связь + интегральные неравенства Лежандра m=3 в рамках одного LMI-подхода. В известной нам литературе (включая работы Fridman, Park, Lee) такая комбинация для интервальных систем с запаздыванием отсутствует.

Для верификации метода синтеза проведено тестирование предлагаемого подхода в среде matlab и сравнение его с двумя известными вариантами решения поставленной задачи.

Параметры проведенного эксперимента по предложенному алгоритму.

Порядок системы n = 2, Интервальные матрицы (численно):

AI = [[-2,1; -1,9], [0,9; 1,1]; [-1,1; -0,9], [-2,1; -1,9]]

AdI = 0,1·AI, BuI = [[0,8; 1,2]; [0,9; 1,1]], μ = 0,1, τm = 0.

Солвер: SeDuMi 1.3.4 + YALMIP R2023a в MATLAB.

Разбиение для варианта В: [0; 0,3; 0,7; 1,0; 1,23] (4 неравномерных отрезка).

При проверке предлагаемого подхода используется интервальная система, дополненная запаздыванием и управлением:

При проведении вычислительного эксперимента сравнивались два варианта:

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

- метод гетерогенного разбиения запаздывания на 4 неравномерных отрезка и использование усовершенствованной формулы оценки неравенств при m=3 (вариант В — предлагаемый подход).

В обоих случаях синтезировался стационарный регулятор (задача максимизации временного запаздывания). Рассмотрено три варианта алгоритмов:

Метод A — равномерное разбиение + Йенсен.

Метод B (предложенный метод) — гетерогенное разбиение + Лежандр m=3.

Метод C — метод Виртингера (Park, Lee, 2015) без гетерогенного разбиения.

Метод D — метод свободных матриц (free-weighting matrices) с равномерным разбиением.

В таблице приведены результаты расчетов τ_max:

В случае с шумом и дискретизацией инженерная значимость может быть ограничена. Полученное увеличение τmax на ~3–7% может быть существенно для систем на границе устойчивости, но требуются дополнительные исследования чувствительности к используемой реализации.

Проведены замеры времени выполнения алгоритмов для системы из примера (n=2, N=4 – количество отрезков разбиения запаздывания):

1. Вариант А (равномерное, Йенсен): время решения LMI ~ 0,23 с.

2. Вариант В (гетерогенное, m=3): время ~ 0,41 с.

3. Рост в ~1,8 раза.

4. Для n=4 рост составляет 2,5 раза.

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

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

К основным научным результатам относятся: разработка функционала Ляпунова-Красовского с гетерогенной декомпозицией, позволяющей учитывать неравномерную чувствительность системы к запаздыванию на разных временных отрезках, и применение усовершенствованных интегральных неравенств порядка m=3, что существенно снижает консервативность условий устойчивости по сравнению с классическими неравенствами Йенсена и Виртингера.

Проведённый вычислительный эксперимент показал, что гетерогенная декомпозиция интервала запаздывания позволяет увеличить максимально допустимое запаздывание по сравнению с классическим методом.

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