Симплектическое снижение порядка модели портогамильтоновых систем

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

• ID статьи: 2203.07751
• Название: Symplectic model order reduction of port-Hamiltonian systems
• Авторы: Silke Glas, Mir Mamunuzzaman, Hongliang Mu, Hans Zwart
• Классификация: math.OC (Оптимизация и управление), math.CA (Классический анализ)
• Дата публикации: Март 2022 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2203.07751

Аннотация

В данной работе предложен новый метод снижения порядка модели (MOR), сохраняющий структуру, для линейных стационарных портогамильтоновых (pH) систем. Целью исследования является построение редуцированной pH-системы, которая остаётся интерпретируемой в физической области исходной полнопорядковой модели. В частности, если исходная высокомерная pH-система представляет собой электрическую цепь, желательно, чтобы редуцированная модель также интерпретировалась как электронная цепь. Для систем масса-пружина-демпфер (MSD) существуют методы MOR, гарантирующие сохранение этой специфической структуры. В работе показано, что новый метод MOR, основанный на симплектическом подходе, восстанавливает известный метод Арнольди второго порядка (SO-Arnoldi) для систем MSD. Однако для моделей pH цепей и других моделей с аналогичной блочной структурой метод даёт новые инновационные методы снижения порядка модели.

Научный контекст и мотивация

Описание проблемы

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

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

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

Научная мотивация

При переформулировании физических моделей в форме pH матрицы J, R и H обычно имеют специальную блочную структуру. Целью данной работы является сохранение этой блочной структуры в процессе редукции, обеспечивая, чтобы редуцированная модель оставалась интерпретируемой как физическая система того же типа, что и полнопорядковая модель.

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

1. Введение симплектико-pH метода: Предложен новый метод MOR для pH-систем на основе симплектического подхода, применимый к блочной форме (2), обеспечивающий, что редуцированная модель сохраняет форму pH-системы
2. Теоретическая связь: Доказано, что для частного случая систем MSD предложенный метод восстанавливает метод SO-Arnoldi
3. Численная верификация: Представлены численные результаты на электронных цепях (LRCR-цепь) с сравнением другими известными методами из литературы
4. Преимущества стабильности: Выявлено, что метод SO-Arnoldi в некоторых случаях производит нестабильные редуцированные модели, тогда как предложенный метод сохраняет стабильность

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

Постановка задачи

Дана линейная стационарная pH-система:

ẋ(t) = (J - R)Hx(t) + Bu(t)
y(t) = B*Hx(t)

где x(t) ∈ C²ⁿ — состояние, u(t) ∈ Cᵖ — входной сигнал, y(t) ∈ Cᵖ — выходной сигнал. Целью является нахождение системы той же формы, но с меньшей размерностью пространства состояний, такой что передаточные функции двух систем совпадают в предопределённых точках частоты.

Основная теоретическая база

Основы симплектической геометрии

Определение 1 (Симплектическая форма): Пусть V — чётномерное векторное пространство над C, симплектическая форма Ω: V × V → C удовлетворяет:

1. Антисимметричность: Ω(v₁, v₂) = -Ω(v₂, v₁)
2. Невырожденность: для всех v₁ ∈ V имеем Ω(v₁, v₂) = 0 тогда и только тогда, когда v₂ = 0

Определение 2 (Симплектическое отображение): Матрица Q называется (J₂ₙ, J₂ₖ)-симплектической, если:

Q*J₂ₙQ = J₂ₖ

Симплектическое MOR для бездиссипативного случая

Теорема 4: Для бездиссипативной pH-системы (R = 0), если Q ∈ C²ⁿˣ²ᵏ является (J₂ₙ⁻¹, J₂ₖ⁻¹)-симплектической матрицей и удовлетворяет:

• ran(B) ⊆ ran(Q)
• (sₘI₂ₙ - J₂ₙH)⁻¹Bμₘ ∈ ran(Q)

то редуцированная pH-система совпадает с исходной системой в заданных точках частоты.

Симплектическое MOR для диссипативного случая

Теорема 5: Для диссипативной pH-системы строится расширенная матрица:

QE := [Q₁  0  ]
      [0   Q₂]

где Q₂ является (J₂ₙ, J₂ₖ)-симплектической, Q₁* является левым обратным к Q₂. При надлежащих условиях редуцированная система удовлетворяет условиям интерполяции.

Алгоритм сохранения блочной структуры

Алгоритм 1 (Симплектико-pH MOR):

1. Построить V₁ на основе данных интерполяции
2. Вычислить V₂ := HV₁
3. Разделить V₁ и V₂ согласно блочной структуре
4. Построить Ψ₁, Ψ₂, Φ₁, Φ₂ по формулам (27)
5. Построить Q₁ и Q₂
6. Применить проекцию Петрова-Галёркина

Формулы построения:

Ψ₂ = orth(V₂₂)
Ψ₁ = V₁₂(Ψ₂*V₁₂)⁻¹
Φ₂ = V₂₁(Ψ₂*Ĵₙ*V₂₁)⁻¹Ĵₖ*
Φ₁ = V₁₁(Φ₂*V₁₁)⁻¹

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

1. Сохранение симплектической структуры: Через симплектические отображения сохраняется фундаментальная геометрическая структура гамильтоновой системы
2. Поддержание блочно-диагональной структуры: Специальное построение Q₁ и Q₂ обеспечивает сохранение блочной структуры физического смысла
3. Условия интерполяции: Гарантирует, что редуцированная система точно совпадает с исходной в предопределённых точках частоты
4. Единая теория: Включает метод SO-Arnoldi для систем MSD как частный случай

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

Численный пример: LRCR-цепь

Рассмотрена LRCR-цепь, показанная на рис. 1, содержащая n=50 элементов с параметрами:

• Ёмкости: Cᵢ = 1 пФ
• Индуктивности: Lᵢ = 1 нГн
• Сопротивления: Rᵢ,₁ = 1 мОм, Rᵢ,₂ = 1 кОм

Критерии оценки

• H₂-ошибка:

  H₂error = √(1/2π ∫ Trace((G(iω) - Gred(iω))*(G(iω) - Gred(iω)))dω)
  

• H∞-ошибка:

  H∞error = max‖G(iω) - Gred(iω)‖₂
  

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

1. pH-Arnoldi: Метод сохранения структуры на основе теоремы 2
2. SO-Arnoldi: Метод Арнольди второго порядка
3. Симплектико-pH: Предложенный в работе метод

Выбор точек интерполяции

Точки интерполяции выбираются равномерно на логарифмической шкале в интервале 10⁻², 10², например, 6 точек: 10⁻²i, -10⁻²i, 10⁰i, -10⁰i, 10²i, -10²i

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

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

Из рис. 3 можно наблюдать:

1. Симплектико-pH vs pH-Arnoldi: Симплектико-pH обеспечивает более низкие H₂ и H∞-ошибки в большинстве случаев
2. Симплектико-pH vs SO-Arnoldi:
   • При редуцированной размерности <24 производительность сопоставима
   • При редуцированной размерности =32 производительность SO-Arnoldi резко падает
   • При редуцированной размерности 28 и 36 SO-Arnoldi производит нестабильные системы

Анализ стабильности

На рис. 4 показаны нестабильные полюса, производимые методом SO-Arnoldi при r=28 и r=36, тогда как метод симплектико-pH сохраняет стабильность при всех тестируемых размерностях.

Теоретическая верификация: эквивалентность для систем MSD

Теорема 7 доказывает, что для систем MSD метод симплектико-pH эквивалентен методу SO-Arnoldi через унитарное преобразование: существует унитарная матрица Θ такая, что zred(t) = Θẑred(t)

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

Развитие MOR для pH-систем

1. Ранние работы: Методы пассивных систем на основе положительно вещественных передаточных функций
2. Согласование моментов: Работы Polyuga и van der Schaft
3. Аппроксимация структур Дирака: Альтернативный подход к сохранению структуры
4. Методы спектрального разложения: Методы на основе функции Попова
5. Оптимизационные методы: Методы на основе параметрической оптимизации и римановой оптимизации

Развитие симплектического MOR

1. Автономные гамильтоновы системы: Пионерские работы Peng и Mohseni
2. Нелинейные методы: Симплектическая редукция на основе многообразий
3. Методы, управляемые данными: Вывод операторов и идентификация систем

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

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

1. Успешно разработан метод MOR, сохраняющий pH-структуру и блочную структуру
2. Теоретически объединены различные методы MOR для систем MSD
3. Численные эксперименты подтверждают эффективность метода и преимущества стабильности

Ограничения

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

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

1. Согласование производных: Расширение на согласование производных передаточной функции для повышения качества аппроксимации
2. Адаптивный выбор: Разработка адаптивных стратегий выбора точек интерполяции
3. Нелинейное расширение: Расширение на нелинейные pH-системы
4. Крупномасштабные приложения: Оптимизация алгоритма для обработки систем большего масштаба

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

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

1. Теоретическая строгость: Прочная математическая база на основе симплектической геометрии
2. Сохранение структуры: Одновременное сохранение pH-структуры и физически значимой блочной структуры
3. Единство: Объединение существующих методов как частных случаев в единую теорию
4. Практичность: Метод имеет чёткую алгоритмическую реализацию и численную верификацию

Недостатки

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

Влияние

1. Теоретический вклад: Предоставляет новую теоретическую перспективу для MOR pH-систем
2. Практическая ценность: Имеет прямую ценность в инженерных приложениях, таких как анализ цепей
3. Методологическое значение: Демонстрирует важность геометрических методов в MOR

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

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

Список литературы

Статья содержит 31 ссылку, охватывающую важные работы в области теории MOR, моделирования pH-систем, симплектической геометрии и теории управления, обеспечивая прочную теоретическую базу для исследования.