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

Основная информация

• ID статьи: 2511.12634
• Название: Approximate Tracking Controllability of Systems with Quadratic Nonlinearities
• Авторы: Manuel Rissel (ShanghaiTech University), Marius Tucsnak (Université de Bordeaux)
• Классификация: math.OC (Оптимизация и управление)
• Дата публикации: 16 ноября 2025 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2511.12634

Аннотация

В данной работе исследуется проблема приближённой управляемости отслеживанием для конечномерных систем управления с непрерывным временем. Для заданного ε>0 изучается вопрос о существовании управления, при котором траектория состояния остаётся в ε-окрестности произвольно заданного пути в пространстве состояний. Авторы проводят систематическое исследование линейных и квадратичных стационарных систем. Для линейных систем получены отрицательные результаты: за исключением тривиального случая, когда пространство управления изоморфно пространству состояний, приближённая управляемость отслеживанием полного состояния невозможна даже в слабом смысле. Опираясь на эти отрицательные результаты для линейных систем, авторы переходят к нелинейной динамике и доказывают слабую приближённую управляемость отслеживанием в произвольных временных интервалах для общего класса систем с произвольной линейной частью и квадратичными нелинейными членами. Рассматриваемое понятие слабой приближённой управляемости отслеживанием включает метрику релаксации. Путём разработки приложений к связанным системам (включая задачи планирования движения) авторы подчёркивают релевантность этого слабого подхода и указывают на препятствия, возникающие в более сильных нормах. Данный подход даёт глобальные результаты даже в случаях, когда неуправляемая динамика может иметь особенности за конечное время.

Исследовательский контекст и мотивация

1. Исследуемая проблема

Традиционная теория управляемости сосредоточена на задачах точка-в-точку, то есть на исследовании того, может ли входной сигнал перевести траекторию состояния из произвольного начального состояния в предопределённое конечное состояние. В отличие от этого, управляемость отслеживанием (tracking controllability) сосредоточена на управлении полным состоянием или выходной функцией системы на всём временном интервале 0,τ таким образом, чтобы следовать предопределённой траектории. Конкретно, в данной работе исследуется следующая система управления:

$\dot{x}(t) + Ax(t) + f(x(t)) = Bu(t), \quad t \in [0,\tau]$

где x — состояние, u — управление, A — линейный оператор, B — оператор управления, f — (возможно) нелинейный член.

2. Значимость проблемы

Проблема управляемости отслеживанием имеет важное значение в нескольких областях приложений:

• Управление сложными сетями: необходимость поддерживать состояние сети вблизи желаемой траектории
• Робототехника: теоретическая основа для отслеживания пути и планирования движения
• Машинное обучение: анализ управляемости непрерывных моделей глубокого обучения, таких как нейронные ОДУ

3. Ограничения существующих методов

• Исследование управляемости выхода относительно хорошо развито (García-Planas и Domínguez-García, Zamorano и Zuazua и др.), но управляемость отслеживанием полного состояния недостаточно изучена
• Существующие исследования управляемости отслеживанием нелинейных систем в основном ограничены аффинными системами управления без дрейфа
• Для конечномерных систем, особенно с квадратичными нелинейностями, отсутствует систематическая теоретическая база

4. Исследовательская мотивация

• Управляемость отслеживанием линейных систем очень ограничена (что доказано в данной работе), что мотивирует исследование нелинейных систем
• Квадратичные нелинейности естественным образом возникают во многих математических моделях (например, система Лоренца, конечномерные аппроксимации уравнений гидродинамики)
• Разложение второго порядка может обеспечить более точное приближение нелинейных систем, чем простая линеаризация

Основные вклады

Основные вклады данной работы включают:

1. Отрицательные результаты для линейных систем: систематически доказано, что слабая приближённая управляемость отслеживанием линейных систем имеет место только в тривиальном случае, когда оператор управления B является сюръекцией (то есть пространство управления изоморфно пространству состояний) (Proposition 1.5)
2. Положительные результаты для квадратичных систем: для квадратичных нелинейных систем, удовлетворяющих определённым предположениям насыщения (Assumption 3.3), доказана слабая приближённая управляемость отслеживанием в произвольных временных интервалах (Theorem 3.5)
3. Ключевая роль релаксационной нормы: введена релаксационная норма |||·|||τ как инструмент измерения, позволяющий получить положительные результаты в этой слабой топологии, при этом контрпримеры показывают, что аналогичные результаты невозможны в более сильной норме L²
4. Приложения к связанным системам: основные результаты применены к задачам управления связанными системами ОДУ, включая задачи динамического планирования движения (Theorem 7.1)
5. Глобальные результаты: данный подход гарантирует глобальную управляемость даже в случаях, когда неуправляемая динамика может взорваться за конечное время

Подробное описание методов

Определение задачи

Определение 1.1 (Приближённая управляемость отслеживанием): Система (1.1) является приближённо управляемой отслеживанием на 0,τ относительно нормы Nτ на L²(0,τ;X), если для каждого ψ∈W^{1,2}((0,τ);X) и ε>0 существует функция управления u∈L²(0,τ;U) такая, что решение x с начальным условием x(0)=ψ(0) удовлетворяет:

$N_\tau(x-\psi) \leq \varepsilon$

Релаксационная норма (ключевой инструмент): для v∈L¹(0,τ;E) определяется как:

$|||v|||_\tau := \sup_{t\in[0,\tau]} \left\|\int_0^t v(s)ds\right\|_E$

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

Анализ линейных систем (Section 2)

Основная идея: посредством анализа двойственности задача управляемости отслеживанием преобразуется в характеризацию нулевого пространства оператора.

Ключевые этапы:

1. Введение оператора Fτ: L²(0,τ;U)→L²(0,τ;Y): $(F_\tau u)(t) = C\int_0^t \exp((t-\sigma)A)Bu(\sigma)d\sigma$
2. Вычисление его сопряжённого оператора Ψτ: $(\Psi_\tau g)(\sigma) = B^*\left[\int_\tau^\sigma \exp((t-\sigma)A^*)C^*g(t)dt\right]$
3. Доказательство Proposition 2.5:
   • Если B не является сюръекцией, существует η∈X{0} такой, что B*η=0
   • Конструируется специальная функция g(t)=-α(t)A*η-α̇(t)η, где α(t)=(τ-t)η
   • Доказывается, что Ψτg=0, но g≠0, следовательно, образ Fτ не плотен
4. Для слабой приближённой управляемости отслеживанием рассматривается расширенная система и получаются аналогичные выводы (Proposition 1.5)

Основные результаты для квадратичных систем (Sections 3-6)

Предположение насыщения (Assumption 3.3):

• Определяется рекурсивная последовательность: E₀=Range(B), Eᵢ=ℱ(Eᵢ₋₁)
• где ℱ(E) — максимальное подпространство такое, что любой γ∈ℱ(E) может быть представлен как γ=ξ₀-∑f(ξᵢ), ξᵢ∈E
• Предполагается, что существует n_X такой, что E_=X

Упрощённое предположение (Assumption 1.7): для каждого γ∈X существуют u,ξ∈U такие, что γ=Bu-f(Bξ)

Стратегия доказательства (основана на методе расширенных систем)

Шаг 1: Расширенная система (Section 4) Вводится система с двумя входами: $\dot{w}(t) + A(w(t)+\zeta(t)) + f(w(t)+\zeta(t)) = \gamma(t)$

где ζ — мультипликативное управление, γ — аддитивное управление. При ζ=0 система вырождается в исходную.

Ключевая лемма (Proposition 4.2): устанавливается непрерывная зависимость решения от данных, обеспечивая устойчивость для последующего приближения.

Шаг 2: Приближение траектории (Section 5, центральная техника)

Идея доказательства Proposition 5.1:

1. Для γ∈E_ℓ (кусочно-постоянные функции, принимающие значения в E_ℓ), предполагается γ=η-f(ξ), ξ,η∈E_{ℓ-1}
2. Используется свойство квадратичности (1.4): $f(w)-\gamma = \frac{1}{2}(f(w+\xi)+f(w-\xi))-η$
3. Конструируется высокочастотное колебательное управление: $\zeta_n(t) = \mathcal{E}_n\xi, \quad \text{где} \quad (\mathcal{E}_n\xi)(t) = \chi(nt/\tau)\xi$
   χ — 1-периодическая функция: χ(t)=ξ на [0,1/2), χ(t)=-ξ на [1/2,1)
4. Лемма 5.2 (ключевая техническая лемма): доказывается, что для непрерывного φ: X→X, $h_n := \phi(w_1+\zeta_n) - \frac{1}{2}(\phi(w_1+\xi)+\phi(w_1-\xi)) + A\zeta_n$
   удовлетворяет: $\lim_{n\to\infty}\|\mathcal{I}h_n\|_{C([0,\tau];X)}=0$
   где $(\mathcal{I}z)(t)=\int_0^t z(s)ds$
5. Доказывается сходимость в релаксационной норме: $|||\zeta_n|||_\tau \to 0 \quad \text{при} \quad n\to\infty$
   Это достигается анализом интеграла $\int_0^t\zeta_n(s)ds$ и использованием периодичности для сокращения членов

Шаг 3: Итеративная конструкция (Section 6)

1. Для целевой траектории ψ определяется γ=ψ̇+Aψ+f(ψ)
2. Применяется Proposition 5.1 для приближения от E_ℓ к E_{ℓ-1}
3. Итерация повторяется n_X раз, в конечном итоге получая управление в E₁=Range(B)
4. Поскольку B:U→Range(B) является сюръекцией, можно найти u∈U, реализующее это управление

Технические инновации

1. Искусное использование релаксационной нормы:
   • Слабее нормы L², позволяет высокочастотные колебания
   • Согласуется с физическим понятием "эффекта усреднения"
   • Контролирует ошибку в интегральном смысле
2. Конструкция высокочастотного колебательного управления:
   • Заимствует метод Agrachev-Sarychev (первоначально применённый к PDE жидкостей)
   • Использует специальную структуру квадратичной нелинейности
   • Реализует "принцип выпуклости" посредством быстрых колебаний
3. Стратегия многоуровневого приближения:
   • Использует рекурсивно определённую последовательность подпространств E_i
   • Каждая итерация снижает уровень на один, в конечном итоге достигая управляемого подпространства
   • Избегает прямых предположений о глобальной управляемости
4. Анализ устойчивости:
   • Proposition 4.2 гарантирует непрерывную зависимость при возмущениях
   • Позволяет приближать негладкие управления гладкими
   • Обеспечивает сходимость итеративного процесса

Экспериментальная установка

Проверка теоретических свойств на примерах

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

Пример 3.1 (Система Лоренца): $\dot{x} = \begin{bmatrix}\sigma & -\sigma & 0\\-\rho & 1 & 0\\0 & 0 & \beta\end{bmatrix}x + \begin{bmatrix}0\\x_2x_3\\-x_1x_2\end{bmatrix} + \begin{bmatrix}1 & 0\\0 & 1\\0 & 1\end{bmatrix}u$

Проверяется, что Assumption 1.7 выполнено, следовательно, система является слабо приближённо управляемой отслеживанием.

Пример 1.9 (Отрицательный результат): Рассматривается система: $\dot{x}_1 + x_1x_2 = u_1, \quad \dot{x}_2 + x_3^2 - x_1^2 = 0, \quad \dot{x}_3 - x_3x_2 = u_2$

• Удовлетворяет предположениям Proposition 1.8, следовательно, является слабо приближённо управляемой отслеживанием
• Однако не является (сильно) приближённо управляемой отслеживанием

Идея доказательства: выбирается целевая траектория ψ=ψ₁,0,0^T, где ψ₁ тождественно равна 1000 на 1/3,2/3. Если предположить существование последовательности приближений в L², то x₁≥999 на этом интервале, откуда ẋ₂≥998000, что противоречит ограниченности нормы L².

Пример 3.4 (Нетривиальность Assumption 3.3):

• X=ℝ⁶, U=ℝ³
• f(x)=0,0,0,x₁x₂,x₁x₃,x₂x₃^T
• B=I₃,O₃^T

Можно проверить, что Assumption 3.3 выполнено, но Assumption 1.7 не выполнено, что показывает, что первое действительно более общее.

Результаты исследования

Основные теоретические результаты

Теорема 3.5 (Основная теорема): Пусть выполнено Assumption 3.3, τ>0. Для любых ψ∈W^{1,2}((0,τ);X) и ε>0 существует управление u∈C^∞(0,τ;U) такое, что решение x с начальным условием x₀=ψ(0) удовлетворяет: $|x(\tau)-\psi(\tau)| + |||x-\psi|||_\tau < \varepsilon$

Следствия:

1. Глобальная приближённая управляемость (частный случай)
2. Даже если неуправляемая система может взорваться за конечное время, управляемая система остаётся глобально определённой

Приложения к связанным системам (Section 7)

Теорема 7.1: Рассматривается связанная система $\dot{x} + Ax + f(x) = Bu, \quad \dot{z} + \tilde{\Gamma}(z,x) + \mathcal{B}(z) = 0$

Если (z̄,x̄) удовлетворяет второму уравнению, и f,B удовлетворяют Assumption 3.3, то для любого ε>0 существует u такое, что: $\sup_{t\in[0,\tau]}|z(t)-\bar{z}(t)| + |x(\tau)-\bar{x}(\tau)| + |||x-\bar{x}|||_\tau < \varepsilon$

Сценарии приложений:

1. Динамическое управление: первое уравнение служит динамическим контроллером для второго
2. Планирование движения (Пример 7.3): управление x для достижения отслеживания траектории z

Конкретная установка Примера 7.3:

• X=Z=ℝ^d, Γ̃(a,b)=a₁b₁,...,a_db_d^T
• Цель: отслеживание эталонной кривой z_ref, остающейся в определённом ортогональном квадранте
• Конструкция: x_ℓ(t)=(𝒷_ℓ(z̄(t))-ż̄_ℓ(t))/z̄_ℓ(t)
• Результат: sup_|z(t)-z_ref(t)|<ε

Теоретические открытия

1. Фундаментальное различие между линейным и нелинейным:
   • Линейные системы: управляемость отслеживанием крайне ограничена
   • Квадратичные системы: богатая управляемость в слабой топологии
2. Критичность выбора нормы:
   • Норма L²: может не работать даже для квадратичных систем (Пример 1.9)
   • Релаксационная норма: позволяет использовать высокочастотные колебания
3. Достаточность предположения насыщения:
   • Assumption 3.3 обеспечивает проверяемое алгебраическое условие
   • Можно определить за конечное число итераций
4. Глобальность vs локальность:
   • Традиционная локальная управляемость за малое время требует сложных условий на скобки Ли
   • Данный метод напрямую даёт глобальные результаты

Связанные работы

Основные направления исследований

1. Классическая теория управляемости:
   • Монография Coron: систематическое изложение управляемости нелинейных систем
   • Метод расширения Jurdjevic-Kupka: методологическая основа данной работы
2. Управляемость выхода отслеживанием:
   • García-Planas & Domínguez-García 12: отслеживание выхода линейных систем
   • Zamorano & Zuazua 26: управляемость отслеживанием конечномерных линейных систем и двойственная характеризация
   • Данная работа: расширение от отслеживания выхода к отслеживанию полного состояния
3. Отслеживание в аффинных системах:
   • Liu & Sussmann 25: приближённое отслеживание в аффинных системах без дрейфа
   • Haynes & Hermes 13: нелинейная управляемость через теорию Ли
   • Данная работа: допускает дрейф (линейный + квадратичный)
4. Локальная управляемость за малое время:
   • Beauchard & Marbach 5: квадратичные препятствия для систем со скалярным входом
   • Данная работа: глобальная управляемость на всём временном интервале
5. Бесконечномерные системы:
   • Agrachev & Sarychev 2,3: управляемость уравнений Эйлера и Навье-Стокса
   • Shirikyan 24, Nersesyan 20,22: приближённая управляемость уравнений жидкостей
   • Koike et al. 15: релаксационное усиление в несжимаемых жидкостях
   • Данная работа: адаптация этих техник к конечномерным квадратичным системам

Преимущества данной работы

1. Систематичность: первое систематическое исследование управляемости отслеживанием полного состояния конечномерных систем
2. Полнота: как отрицательные (линейные), так и положительные (квадратичные) результаты
3. Проверяемость: обеспечены алгебраически проверяемые достаточные условия
4. Применимость: развиты приложения к связанным системам и планированию движения

Заключение и обсуждение

Основные выводы

1. Ограничения линейных систем: слабая приближённая управляемость отслеживанием ⟺ оператор управления является сюръекцией (Proposition 1.5)
2. Возможности квадратичных систем: при предположении насыщения слабая приближённая управляемость отслеживанием имеет место в произвольных временных интервалах (Theorem 3.5)
3. Значимость нормы:
   • Релаксационная норма: достижима
   • Норма L²: в общем случае недостижима (Пример 1.9)
   • Нормы L^p (p≥1): аналогичные препятствия (Замечание 1.10)
4. Прикладная ценность: применимо к динамическому управлению и планированию движения (Theorem 7.1)

Ограничения

1. Существенность квадратичного свойства:
   • Метод сильно зависит от квадратичной структуры
   • Неквадратичные возмущения обрабатываются только при "малых" предположениях (конец Section 1)
   • Расширение на полиномы высших степеней нетривиально
2. Необходимость слабой топологии:
   • В более сильных нормах (например, L²) результаты не имеют места
   • Физические приложения могут требовать более сильного управления
3. Проверка предположения насыщения:
   • Assumption 3.3 требует вычисления рекурсивной последовательности
   • Для высокомерных систем может быть вычислительно сложным
4. Требования гладкости:
   • Целевая траектория должна быть ψ∈W^{1,2}
   • Обработка негладких целей не рассмотрена
5. Стоимость управления:
   • Не обсуждается норма L² управления или энергия
   • Высокочастотные колебания могут привести к высокому энергопотреблению

Направления будущих исследований

Направления, намеченные в статье:

1. Полиномы высших степеней: расширение на кубические или полиномы более высокого порядка
2. Бесконечномерный случай: приложение к Galerkin-аппроксимациям PDE
3. Ограничения на управление: рассмотрение ограниченного управления или других ограничений
4. Оптимальность: исследование задач минимальной энергии или оптимального управления
5. Робастность: отслеживание, робастное к неопределённостям модели

Глубокая оценка

Преимущества

1. Теоретическая глубина:
   • Полные результаты как в положительном, так и отрицательном направлениях
   • Строгие математические доказательства
   • Глубокое раскрытие фундаментального различия между линейным и нелинейным
2. Методологические инновации:
   • Искусная адаптация техник управления PDE (метод Agrachev-Sarychev) к ODE
   • Введение релаксационной нормы имеет как математическое, так и физическое значение
   • Стратегия многоуровневого приближения элегантна и операциональна
3. Ясная структура:
   • От простого к сложному: линейное→квадратичное
   • От частного к общему: Assumption 1.7→Assumption 3.3
   • От теории к приложениям: основная теорема→связанные системы
4. Богатство примеров:
   • Система Лоренца (Пример 3.1)
   • Контрпримеры (Пример 1.9)
   • Нетривиальные предположения (Пример 3.4)
   • Планирование движения (Пример 7.3)
5. Технические детали:
   • Доказательство Lemma 5.2 изящно (использует относительную компактность + поточечную сходимость)
   • Proposition 4.2 гарантирует устойчивость процесса приближения
   • Анализ двойственности (Section 2) ясен и прозрачен

Недостатки

1. Практические соображения:
   • Отсутствуют численные эксперименты для проверки теоретических предсказаний
   • Не обсуждается практическая реализуемость управления (энергия, частота и т.д.)
   • Осуществимость высокочастотного колебательного управления в реальных системах вызывает сомнения
2. Ограничивающие предположения:
   • Хотя Assumption 3.3 проверяема, она может не выполняться для некоторых систем
   • Ограничение квадратичного свойства довольно сильно
   • Не обсуждается необходимость предположений
3. Сила результатов:
   • Получены только результаты в слабой топологии
   • Пример 1.9 показывает, что нельзя обобщить на норму L²
   • Управление конечным состоянием только приближённое
4. Полнота анализа:
   • Не обсуждается связь n_X (числа итераций) с размерностью системы
   • Отсутствуют оценки нормы L² управления
   • Зависимость от ε не уточнена
5. Глубина приложений:
   • Приложения в Section 7 относительно просты
   • Отсутствует сравнение с практическими алгоритмами планирования движения
   • Не проведён анализ вычислительной сложности

Влияние на область

Вклад в область:

1. Первое систематическое исследование управляемости отслеживанием полного состояния конечномерных систем
2. Установлена полная теория невозможности для линейных систем
3. Обеспечены операциональные достаточные условия для квадратичных систем
4. Связь между конечномерной теорией ODE и бесконечномерной теорией PDE

Практическая ценность:

1. Теоретическая основа для планирования движения
2. Новые идеи для проектирования динамических контроллеров
3. Возможные приложения в робототехнике и автоматическом управлении

Воспроизводимость:

• Теоретические результаты полностью воспроизводимы (чистые математические доказательства)
• Примеры проверяемы напрямую (алгебраические вычисления)
• Численная реализация требует дополнительной работы (численное приближение высокочастотного управления)

Ограничения:

• Практическое значение слабой топологии требует дальнейшего уточнения
• Квадратичное свойство ограничивает прямую применимость
• Требуются численные эксперименты для поддержки теоретических результатов

Применимые сценарии

Теоретически подходящие сценарии:

1. Механические системы с квадратичной нелинейностью: например, с демпфированием, содержащим квадратичные члены
2. Сети химических реакций: кинетика реакций второго порядка
3. Биологические системы: модели типа Lotka-Volterra
4. Низкомерные модели жидкостей: Galerkin-усечение

Сценарии, требующие осторожности:

1. Строгие требования к управлению: когда требуется отслеживание в L² или L^∞, метод неприменим
2. Системы с быстрым откликом: высокочастотные колебания могут быть нереализуемы
3. Системы с ограниченной энергией: стоимость управления может быть чрезмерной
4. Системы с доминирующей нелинейностью высокого порядка: квадратичное приближение недостаточно

Потенциальные направления расширения:

1. Интеграция с численной оптимизацией для разработки практических алгоритмов
2. Рассмотрение ограниченного управления
3. Исследование робастности и чувствительности
4. Расширение на системы с изменяющимися параметрами или стохастические системы

Библиография

Ключевые цитируемые работы:

1. 2,3 Agrachev & Sarychev: управляемость уравнений 2D Euler и Navier-Stokes, основной источник вдохновения для методов данной работы
2. 26 Zamorano & Zuazua: управляемость отслеживанием конечномерных линейных систем, обеспечивает основу для анализа двойственности
3. 15 Koike et al.: релаксационное усиление в несжимаемых жидкостях, демонстрирует применение релаксационной нормы в PDE
4. 24 Shirikyan: приближённая управляемость уравнения 3D Navier-Stokes, систематическое изложение принципа выпуклости
5. 14 Jurdjevic & Kupka: системы управления под действием групп, теоретическая основа метода расширения
6. 8 Coron: классический учебник по нелинейной управляемости

Общая оценка: Это высокого качества теоретическая работа по теории управления, систематически устанавливающая теоретическую базу управляемости отслеживанием полного состояния конечномерных систем. Путём введения релаксационной нормы и использования специальной структуры квадратичной нелинейности авторы искусно адаптировали техники управления бесконечномерными PDE к конечномерным ODE, получив глубокие положительные и отрицательные результаты. Основная ценность работы заключается в теоретическом вкладе и методологических инновациях, однако практическая применимость и численная проверка требуют дополнительной работы. Для исследователей, работающих в области нелинейной теории управления, это обязательная к прочтению статья; для прикладных исследователей необходимо оценить применимость к конкретным задачам.