2025-11-10T02:47:53.458764

KLAP: KYP lemma based low-rank approximation for $\mathcal{H}_2$-optimal passivation

Nicodemus, Voigt, Gugercin et al.
We present a novel passivity enforcement (passivation) method, called KLAP, for linear time-invariant systems based on the Kalman-Yakubovich-Popov (KYP) lemma and the closely related Lur'e equations. The passivation problem in our framework corresponds to finding a perturbation to a given non-passive system that renders the system passive while minimizing the $\mathcal{H}_2$ or frequency-weighted $\mathcal{H}_2$ distance between the original non-passive and the resulting passive system. We show that this problem can be formulated as an unconstrained optimization problem whose objective function can be differentiated efficiently even in large-scale settings. We show that any minimizer of the unconstrained problem yields the same passive system. Furthermore, we prove that, in the absence of a feedthrough term, every local minimizer is also a global minimizer. For cases involving a non-trivial feedthrough term, we analyze global minimizers in relation to the extremal solutions of the Lur'e equations, which can serve as tools for identifying local minima. To solve the resulting numerical optimization problem efficiently, we propose an initialization strategy based on modifying the feedthrough term and a restart strategy when it is likely that the optimization has converged to a non-global local minimum. Numerical examples illustrate the effectiveness of the proposed method.
academic

KLAP: Низкоранговая аппроксимация на основе леммы КЯП для H2\mathcal{H}_2-оптимальной пассификации

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

  • ID статьи: 2501.05178
  • Название: KLAP: KYP lemma based low-rank approximation for H2\mathcal{H}_2-optimal passivation
  • Авторы: Jonas Nicodemus, Matthias Voigt, Serkan Gugercin, Benjamin Unger
  • Классификация: math.OC (математическая оптимизация и управление)
  • Дата публикации: 14 октября 2025 г. (препринт arXiv)
  • Ссылка на статью: https://arxiv.org/abs/2501.05178

Аннотация

В данной работе предложен новый метод пассификации под названием KLAP для линейных стационарных систем на основе леммы Калмана-Якубовича-Попова (КЯП) и связанных уравнений Люре. В рамках предложенного подхода задача пассификации соответствует поиску возмущения заданной непассивной системы, которое делает систему пассивной, одновременно минимизируя расстояние в норме H2\mathcal{H}_2 или частотно-взвешенной норме H2\mathcal{H}_2 между исходной непассивной системой и результирующей пассивной системой. Показано, что данная задача может быть сформулирована как задача безусловной оптимизации, целевая функция которой допускает эффективное дифференцирование даже в крупномасштабных постановках. Доказано, что любой минимизатор безусловной задачи порождает одну и ту же пассивную систему, и при отсутствии прямого канала каждый локальный минимизатор также является глобальным минимизатором.

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

Важность проблемы

  1. Требования к моделированию физических систем: Системы в электротехнике, энергетике, механике и пороупругости требуют гарантий пассивности для получения физически значимых результатов моделирования
  2. Требования к сетевым взаимосвязям: Пассивные системы служат строительными блоками крупномасштабных сетевых моделей, и сохраняющие мощность взаимосвязи пассивных систем порождают глобально пассивную систему
  3. Практические вызовы моделирования: Несмотря на то, что физические процессы по своей природе пассивны, модели, полученные посредством неструктурированных методов редукции порядка или методов идентификации систем на основе данных, часто оказываются непассивными

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

Существующие методы пассификации подразделяются на три основные категории:

  1. Методы на основе линейных матричных неравенств (ЛМН) с использованием леммы КЯП: Вычислительные затраты быстро растут с увеличением размерности системы из-за необходимости существования матрицы Ляпунова
  2. Методы на основе спектральных характеристик матрицы Гамильтона: Отсутствуют гарантии сходимости, может потребоваться множество итераций
  3. Методы на основе дискретных частот: Гарантируют пассивность только в определённых диапазонах частот

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

Данная работа направлена на разработку эффективного метода пассификации, способного:

  • Обрабатывать крупномасштабные системы
  • Обеспечивать гарантии сходимости
  • Находить оптимальные решения в смысле нормы H2\mathcal{H}_2

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

  1. Явная параметризация: Использование существования решения с минимальным рангом неравенства КЯП для получения явной параметризации любой пассивной системы с nmnm переменными решения
  2. Переформулировка в безусловную оптимизацию: Преобразование задачи выпуклой оптимизации с ограничениями в невыпуклую задачу безусловной оптимизации с установлением разрешимости, единственности и методов вычисления градиента
  3. Теория глобальной оптимальности: Доказательство того, что при кососимметричном прямом канале (D+DT=0D + D^T = 0) любой локальный минимизатор также является глобальным минимизатором
  4. Обнаружение локальной оптимальности: Предоставление нового критерия использования экстремальных решений неравенства КЯП для проверки того, является ли локальный минимизатор глобальным минимизатором
  5. Практические стратегии алгоритма: Предложение стратегий инициализации и перезапуска на основе модификации прямого канала

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

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

Дана линейная стационарная динамическая система: Σ:{x˙(t)=Ax(t)+Bu(t)y(t)=Cx(t)+Du(t)\Sigma : \begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = Cx(t) + Du(t) \end{cases}

Целью является нахождение модифицированной системы: Σ^(C^):{x˙(t)=Ax(t)+Bu(t)y(t)=C^x(t)+Du(t)\hat{\Sigma}(\hat{C}) : \begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = \hat{C}x(t) + Du(t) \end{cases}

такой, что Σ^(C^)\hat{\Sigma}(\hat{C}) является пассивной и минимизирует расстояние в норме H2\mathcal{H}_2 от исходной системы.

Основы теории

Лемма КЯП и параметризация пассивности

На основе леммы КЯП система является пассивной тогда и только тогда, когда существуют матрицы L∈Rn×mL \in \mathbb{R}^{n \times m} и M∈Rm×mM \in \mathbb{R}^{m \times m} такие, что: C=BTL−1(−LLT)+MLTC = B^T\mathcal{L}^{-1}(-LL^T) + ML^TD+DT=MMTD + D^T = MM^T

где L\mathcal{L} — оператор Ляпунова: L(X)=ATX+XA\mathcal{L}(X) = A^TX + XA.

Целевая функция и градиент

Целевая функция может быть представлена как: J(L)=tr((C−C^(L))P(CT−C^(L)T))J(L) = \text{tr}((C - \hat{C}(L))P(C^T - \hat{C}(L)^T))

где PP — грамиан управляемости. Градиент имеет вид: ∇J(L)=2XL−2P(CT−C^(L)T)M\nabla J(L) = 2XL - 2P(C^T - \hat{C}(L)^T)M

Архитектура алгоритма

Процедура алгоритма KLAP

  1. Инициализация: Получение начального L0L_0 с использованием алгоритма 1
  2. Оптимизация: Решение безусловной задачи методом L-BFGS
  3. Обнаружение глобальности: Проверка собственных значений Y∗=A−B(D+DT)−1M(L∗)TY^* = A - B(D+D^T)^{-1}M(L^*)^T
  4. Стратегия перезапуска: При обнаружении локальной оптимальности выполнение шага градиента и перезапуск

Стратегия инициализации

Возмущение прямого канала DD для обеспечения пассивности системы:

  • Вычисление λmin⁡=min⁡ωλmin⁡(Φ(iω))\lambda_{\min} = \min_\omega \lambda_{\min}(\Phi(i\omega))
  • Установка Dpert=D−(λmin⁡/2−ϵ)ImD_{\text{pert}} = D - (\lambda_{\min}/2 - \epsilon)I_m
  • Решение соответствующего алгебраического уравнения Риккати для инициализации

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

Тестовые системы

  1. Задача ACC: Система малого размера (n=4,m=1n=4, m=1)
  2. Рычаг CD-проигрывателя: Система среднего размера (n=120,m=2n=120, m=2)
  3. Высокоскоростная линия межсоединений смартфона: Крупномасштабная система (n=800,m=4n=800, m=4)

Методы сравнения

  • ЛМН: Стандартный метод на основе леммы КЯП
  • ЛМН-ТП: Метод ЛМН с параметризацией следа
  • Метод Гамильтона: Метод на основе возмущения собственных значений матрицы Гамильтона

Показатели оценки

  • Ошибка в норме H2\mathcal{H}_2: ∥G−G^(⋅;C^)∥H2\|G - \hat{G}(\cdot; \hat{C})\|_{\mathcal{H}_2}
  • Время вычисления и количество итераций
  • Частота сходимости к глобальному оптимуму

Результаты экспериментов

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

Сравнение вычислительной эффективности

МодельМетодИтерацииОбщее время (с)Время на итерацию (с)Ошибка H2\mathcal{H}_2
ACCKLAP122.29×10⁻⁴1.91×10⁻⁵8.71×10⁻¹
ACCЛМН134.61×10⁻³3.54×10⁻⁴8.71×10⁻¹
ACCЛМН-ТП113.59×10⁻²3.26×10⁻³8.71×10⁻¹
CD-проигрывательKLAP305.44×10⁻¹1.81×10⁻²1.06×10⁶
CD-проигрывательЛМН-ТП1166.04×10²5.21×10⁰1.00×10⁶
СмартфонKLAP22081.46×10²6.63×10⁻²8.32×10⁵

Ключевые выводы

  1. Вычислительная эффективность: Метод KLAP на 1-2 порядка быстрее традиционных методов ЛМН
  2. Глобальная сходимость: При отсутствии прямого канала все локальные оптимумы являются глобальными
  3. Эффективность стратегии перезапуска: Стратегия перезапуска успешно восстанавливает сходимость из негло­бальных локальных оптимумов
  4. Применимость к крупномасштабным системам: Метод эффективно работает на системах размерности 800

Анализ конкретных случаев

Задача ACC

  • Без прямого канала: все инициализации сходятся к глобальному оптимуму
  • С прямым каналом: 40% случайных инициализаций сходятся к негло­бальному локальному оптимуму
  • После применения стратегии перезапуска: все инициализации сходятся к глобальному оптимуму

Высокоскоростная линия межсоединений смартфона

  • Улучшение ошибки H2\mathcal{H}_2 на ~31% по сравнению с методом сравнения
  • Благодаря диагонализирующему преобразованию время решения одного уравнения Ляпунова сокращено с 550 мс до 4 мс

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

Классификация методов пассификации

  1. Методы на основе леммы КЯП: Порождают задачи выпуклой оптимизации, но имеют высокие вычислительные затраты
  2. Методы на основе спектра матрицы Гамильтона: Отсутствуют гарантии сходимости
  3. Методы на основе частотной дискретизации: Эффективны только в определённых диапазонах частот

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

  • Избегаются крупномасштабные вычисления ЛМН
  • Обеспечиваются теоретические гарантии сходимости
  • Применимо к крупномасштабным системам
  • Имеются явные критерии глобальной оптимальности

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

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

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

Ограничения

  1. Для нетривиальных прямых каналов могут существовать множественные локальные оптимумы
  2. Требуется предположение об асимптотической устойчивости системы
  3. В настоящее время ориентировано главным образом на оптимизацию в норме H2\mathcal{H}_2

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

  1. Расширение на ограниченную вещественную лемму для нахождения ближайшей сжимающей системы
  2. Применение к параметризованным системам и дифференциально-алгебраическим уравнениям
  3. Исследование задачи H∞\mathcal{H}_\infty-оптимальной пассификации

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

Достоинства

  1. Надёжный теоретический вклад: Обеспечивается полный теоретический анализ, включая существование, единственность и глобальную оптимальность
  2. Высокая инновационность метода: Умелое использование низкоранговой факторизации неравенства КЯП позволяет избежать вычислительных узких мест традиционных методов
  3. Выдающаяся практическая применимость: Алгоритм легко реализуется и применим к крупномасштабным системам
  4. Достаточная экспериментальная валидация: Метод проверен на эталонных системах различных размеров

Недостатки

  1. Проблема локальных оптимумов: Для общих прямых каналов возможно застревание в локальных оптимумах
  2. Зависимость от инициализации: Производительность метода в определённой степени зависит от качества инициализации
  3. Неполнота теоретического анализа: Анализ случая D+DT⊁0D + D^T \not\succ 0 недостаточно полный

Влияние

  1. Академическая ценность: Предоставляет новую теоретическую перспективу и методы решения задачи пассификации
  2. Практическая ценность: Особенно применимо к потребностям пассификации крупномасштабных инженерных систем
  3. Воспроизводимость: Код и данные открыты, что облегчает проверку и применение

Области применения

  1. Пассификация крупномасштабных линейных систем
  2. Восстановление пассивности после редукции порядка модели
  3. Постобработка при идентификации систем на основе данных
  4. Проектирование сетевых взаимосвязанных систем

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

Статья цитирует 58 связанных работ, охватывающих главным образом:

  • Основы теории диссипативных систем Willems, 1972
  • Лемма КЯП и теория положительной вещественности Anderson & Vongpanitlerd, 1973
  • Обзор методов пассификации Grivet-Talocia & Gustavsen, 2016
  • Методы численной оптимизации Boyd et al., 1989