Вычисление жесткого масштабированного относительного графика линейных стационарных систем

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

• ID статьи: 2511.17297
• Название: Computing the Hard Scaled Relative Graph of LTI Systems
• Авторы: Julius P. J. Krebbekx, Eder Baron-Prada, Roland Tóth, Amritam Das
• Учреждения: Технологический университет Эйндховена, Австрийский институт технологий, ETH Zurich, Институт информатики и управления HUN-REN
• Классификация: eess.SY (Системы и управление), cs.SY, math.OC
• Дата подачи: 21 ноября 2025 г.
• Ссылка на статью: https://arxiv.org/abs/2511.17297

Аннотация

Масштабированные относительные графики (SRG) предоставляют новый графико-частотный метод анализа нелинейных систем, где линейные стационарные (LTI) системы являются основными строительными блоками. Для анализа замкнутых контуров обратной связи, содержащих нестабильные LTI компоненты, необходимо использовать жесткий SRG, поскольку он надлежащим образом отражает поведение входа/выхода в расширенном пространстве L₂. В данной статье разработан систематизированный метод вычисления для точного расчета жесткого SRG LTI систем, которые могут быть нестабильными и содержать интеграторы. Статья также исследует связь с критерием Найквиста, включая многомерный случай, и демонстрирует метод на нескольких примерах.

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

Исследовательская проблема

Статья направлена на решение центральной проблемы анализа нелинейных систем: как систематически вычислить масштабированный относительный график (SRG) замкнутых систем обратной связи, содержащих нестабильные LTI компоненты.

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

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

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

1. Мягкий SRG: Первоначально предложен в 3, применим только к стабильным SISO системам, имеет проблемы с нестабильными системами
2. Расширенный SRG: Разработан в 11 путем добавления информации критерия Найквиста для решения проблемы SISO нестабильных систем, но ограничен одномерным случаем
3. Жесткий SRG: Определен в 12, теоретически может обрабатывать нестабильные системы, но не имеет систематического метода вычисления

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

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

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

1. Систематизированный метод вычисления: Предложен неприближенный алгоритм (Теорема 2), позволяющий точно вычислить жесткий SRG LTI систем, включая нестабильные системы и системы с интеграторами
2. Доказательство теоретической эквивалентности: Доказано, что для SISO систем жесткий SRG полностью эквивалентен расширенному SRG (Теорема 3)
3. Расширение на MIMO: Метод распространен на многомерный (MIMO) случай с полной структурой вычисления для квадратных систем
4. Связь с критерием Найквиста: Показано, как жесткий SRG служит альтернативным методом для критерия MIMO Найквиста и обобщенного критерия Найквиста (GNC)
5. Практический алгоритм: На основе представления передаточной функции граница конструируется путем пересечения и удаления дисков, центрированных на вещественной оси в комплексной плоскости

Детальное описание метода

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

Дана LTI система G: L²ₑᵖ → L²ₑᵍ с передаточной функцией G(s) ∈ Rᵍˣᵖ(s), вычислить её жесткий SRG:

$\text{SRG}_e(G) = \bigcup_{T>0} \text{SRG}(G|_T)$

где G|_T — оператор G, усеченный на временном интервале 0,T.

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

1. Представление гиперболической выпуклости

Замыкание любого h-выпуклого множества C может быть представлено как:

$\text{cl}\,C = \bigcap_{\alpha \in \mathbb{R}} \left(D_{R_\alpha}(\alpha) \setminus D_{r_\alpha}(\alpha)\right) =: \mathcal{G}(\{r_\alpha\}, \{R_\alpha\})$

где $D_r(x)$ обозначает замкнутый диск с центром в x и радиусом r.

2. Определение сдвинутой системы

Для α ∈ ℝ определяется сдвинутая система:

$G_\alpha = \begin{pmatrix} G \\ 0_{(n-q)\times p} \end{pmatrix} - \begin{pmatrix} \alpha I \\ 0_{(n-p)\times p} \end{pmatrix}$

где n = max{p,q}.

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

Предложение 2: Вычисление радиуса максимального усиления Rα

$R_\alpha = \begin{cases} \bar{\sigma}(G_\alpha) & \text{если } G \in RH_\infty^{q\times p} \\ \infty & \text{в противном случае} \end{cases}$

Ключевая идея:

• Для стабильных систем используется норма H∞ (максимальное сингулярное значение)
• Для нестабильных или неправильных систем Rα = ∞

Предложение 3: Вычисление радиуса минимального усиления rα

$r_\alpha = \begin{cases} \underline{\sigma}(G_\alpha) & \text{если } G_\alpha \text{ имеет минимальную фазу} \\ 0 & \text{в противном случае} \end{cases}$

Ключевая идея:

• Через вычисление обратной системы $G_\alpha^{-1}$ и применение Предложения 2
• rα отлично от нуля только когда Gα имеет минимальную фазу (все передаточные нули в открытой левой полуплоскости)
• Использование формы Смита-Макмиллана: полюсы Gα являются нулями $G_\alpha^{-1}$ и наоборот

Основной алгоритм (Теорема 2)

Для квадратной LTI системы G: L²ₑᵖ → L²ₑᵖ с передаточной функцией G(s) ∈ Rᵖˣᵖ(s):

$\text{SRG}_e(G) = \mathcal{G}(\{r_\alpha\}, \{R_\alpha\})$

Этапы алгоритма:

1. Для каждого α ∈ ℝ построить Gα = G - αI
2. Использовать Предложение 2 для вычисления Rα
3. Использовать Предложение 3 для вычисления rα
4. Вычислить пересечение всех кольцевых областей

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

1. Обработка расширенного пространства L₂

• Традиционные методы работают в пространстве L₂, не могут обрабатывать расходящиеся сигналы
• Жесткий SRG определен в расширенном пространстве L²ₑ, включающем периодические и расходящиеся сигналы
• Обработка нестабильности через рассмотрение усеченных сигналов для всех T > 0

2. Обнаружение минимальной фазы

• Ключевая инновация: вычисление rα сводится к проверке, имеет ли Gα минимальную фазу
• Минимальная фаза ⟺ все передаточные нули в открытой левой полуплоскости
• Для неминимальнофазовых систем rα = 0, отражая потенциальные проблемы в обратной связи

3. Унификация с критерием Найквиста

Для SISO систем доказано: $\text{SRG}_e(G) = \text{SRG}'(G)$

где SRG'(G) — расширенный SRG, содержащий информацию мягкого SRG и критерия Найквиста.

4. Обработка специальных случаев

• Интеграторы: Автоматически приводят к Rα = ∞
• Неправильные системы: Обнаруживаются через форму Смита-Макмиллана при дефиците ранга
• Временные задержки: Приводят к rα = 0, Rα не затронут

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

Примеры систем

SISO системы

1. G₁(s) = 1/(s² + s + 1): Стабильная система второго порядка
2. G₂(s) = 1/(s(s+1)): Система с интегратором
3. G₃(s) = 3/((s-2)(s/10+1)): Нестабильная система (полюс в правой полуплоскости)

MIMO системы

1. G₄(s): Стабильная система 2×2 $G_4(s) = \begin{pmatrix} \frac{1}{s+1} & \frac{1}{s+2} \\ \frac{1}{s+4} & \frac{1}{s+3} \end{pmatrix}$
2. G₅(s): Нестабильная система 2×2 $G_5(s) = \begin{pmatrix} \frac{s+7}{s-1} & \frac{s-5}{(s+2)^2} \\ \frac{1}{(s+4)^3} & \frac{s}{(s+3)^2} \end{pmatrix}$
3. G₆(s): Система 3×3 с интегратором

Методы визуализации

• Построение границы жесткого SRG на комплексной плоскости (серая область)
• Наложение границы мягкого SRG (темно-серая область)
• Маркировка критических точек, таких как точка -1 (для анализа стабильности)

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

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

Результаты для SISO систем

1. G₁: Жесткий SRG совпадает с мягким SRG (система стабильна), образует ограниченную область
2. G₂: Из-за интегратора жесткий SRG распространяется на бесконечность, включает полуось [0,∞)
3. G₃: Жесткий SRG нестабильной системы включает дополнительные области, отражая информацию об охвате Найквиста

Проверка эквивалентности расширенного SRG:

• Для G₂, $\mathcal{N}_{G_2} \cap \mathbb{R} = [0,\infty)$
• Это полностью совпадает с жестким SRG на рисунке 2b
• Подтверждает корректность Теоремы 3

Результаты для MIMO систем

G₄ (стабильная система 2×2):

• Границы жесткого и мягкого SRG очень близки
• Обе отделены от точки -1, указывая на стабильность замкнутого контура

G₅ (нестабильная система 2×2):

• Жесткий SRG показывает два отверстия (рисунок 2e)
• Эти отверстия соответствуют областям стабильной конфигурации замкнутого контура
• Меньшее отверстие может использоваться для проектирования коэффициента усиления

G₆ (система 3×3 с интегратором):

• Жесткий SRG частично неограничен (из-за интегратора)
• Граница мягкого SRG стремится к ±j∞ в мнимом направлении

Применение анализа стабильности

Пример проектирования коэффициента усиления (на основе G₅)

Статья демонстрирует три стратегии выбора коэффициента усиления:

1. Чистое масштабирование: k₁ = 0, k₂ = 1/5
   • Интервал [-6, 0) ⊈ SRGₑ(G₅)
   • После масштабирования точка -1 отделена от SRG
2. Отражение и масштабирование: k₁ = 0, k₂ = -2
   • Отражает правое маленькое отверстие на точку -1
   • Использует геометрические свойства жесткого SRG
3. Трансляция: k₁ = -1.5, k₂ = 1
   • Сдвигает маленькое отверстие влево, чтобы оно включало точку -1
   • Гарантирует стабильность замкнутого контура

Каждый случай удовлетворяет условию стабильности: $\text{dist}(-1, \text{SRG}_e(k_1I + k_2G_5)) \geq r > 0$

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

1. Точность: Для SISO и нормальных матричных MIMO систем граница точна (неприближенна)
2. Численная стабильность: Конечное множество α приводит к численной границе, которая всегда является внешней аппроксимацией, не недооценивая SRG
3. Геометрическая интуитивность: Жесткий SRG обеспечивает интуитивное геометрическое объяснение, визуализирует запасы стабильности
4. Модульность: После вычисления жесткого SRG системы G его можно повторно использовать для различных проектов коэффициентов усиления без пересчета

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

Развитие структуры SRG

1. Исходный SRG 2,3: Предложен Ryu и др., Chaffey и др., основан на мягком SRG в пространстве L₂
2. Расширенный SRG 4,11: Расширен путем добавления информации критерия Найквиста на нестабильные SISO системы
3. Жесткий SRG 12: Предложен Chen, определен на L₂0,T, теоретически может обрабатывать нестабильные системы

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

1. Анализ нелинейности: Обобщенный круговой критерий 4, малая теорема усиления 3
2. Управление с переключением: Анализ 8 и проектирование 9
3. Распределенные системы: Условия стабильности 10
4. Частотные характеристики: Нелинейные диаграммы Боде и определение полосы пропускания 5
5. Фазовый анализ: Концепции опережения/отставания фазы 6

Связь с классическими методами

1. Критерий Найквиста: Жесткий SRG в SISO случае эквивалентен расширенному SRG, который содержит информацию Найквиста
2. Критерий MIMO Найквиста 13: Жесткий SRG предоставляет альтернативный метод без необходимости вычисления числа охватов
3. Обобщенный критерий Найквиста (GNC) 14: Жесткий SRG избегает проблемы ручного сшивания характеристических траекторий

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

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

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

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

1. Полнота алгоритма: Разработан первый систематизированный метод для точного вычисления жесткого SRG LTI систем (включая нестабильные системы)
2. Теоретическая унификация: Доказана эквивалентность жесткого SRG и расширенного SRG в SISO случае, унифицированы две структуры
3. Практическая ценность: Жесткий SRG предоставляет практическую альтернативу критерию MIMO Найквиста и GNC со следующими преимуществами:
   • Не требует вычисления числа охватов
   • SRG можно переиспользовать при изменении коэффициента усиления
   • Непосредственно обеспечивает границы L₂ усиления и запасы стабильности
4. Геометрическая интуитивность: Геометрическая конструкция на основе пересечения дисков обеспечивает четкую визуализацию и интуитивное понимание

Ограничения

1. Зависимость от модели

• Текущий метод основан на представлении передаточной функции, требует параметрической модели системы
• Преимущество критерия Найквиста в том, что он может использовать данные частотного отклика, полученные экспериментально
• Будущая потребность: Прямое вычисление rα и Rα из входных/выходных данных

2. Сложность обнаружения неминимальной фазы

• Предложение 2 требует только проверки наличия нестабильных полюсов (не требует подсчета)
• Предложение 3 требует обнаружения нулей неминимальной фазы
• Известная сложность: Обнаружение нулей неминимальной фазы из данных является сложной задачей 22

3. Неквадратные системы

• Текущая полная теория применима только к квадратным системам (p = q)
• Широкий случай (p > q): Gα всегда имеет нетривиальное алгебраическое ядро, приводящее к rα = 0
• Высокий случай (p < q): Можно использовать левый псевдообратный, но требует значительных дополнительных технических деталей

4. Предположение о корректности

• Анализ стабильности жесткого SRG требует предположения о корректности (well-posedness) системы обратной связи
• Мягкий и расширенный SRG помимо границ усиления гарантируют корректность
• Это является недостатком жесткого SRG относительно расширенного SRG

Будущие направления

1. Методы, управляемые данными

Прямое вычисление жесткого SRG из частотного отклика или входных/выходных данных, в частности:

• Обнаружение нестабильных полюсов из данных
• Идентификация нулей неминимальной фазы из данных

2. Расширение на неквадратные системы

• Завершение теории для высоких систем (p < q)
• Разработка альтернативных методов для широких систем

3. Применение к нелинейным системам

• Применение результатов к анализу нелинейных MIMO систем 7
• Интеграция с анализом нелинейной обратной связи жесткого SRG из 12

4. Оптимизация вычисления

• Разработка эффективных численных алгоритмов
• Обработка вычислительной сложности для высокомерных MIMO систем
• Адаптивный выбор множества α для баланса между точностью и вычислительной стоимостью

5. Анализ робастности

• Расширение на неопределенные системы
• Интеграция с теорией робастного управления

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

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

1. Теоретическая строгость

• Полные математические доказательства: Все основные результаты (Предложения 1-3, Теоремы 2-3) имеют строгие доказательства
• Геометрический фундамент: Основана на прочном теоретическом фундаменте гиперболической выпуклости
• Доказательство эквивалентности: Теорема 3 устанавливает четкую связь с классическими методами

2. Инновационность метода

• Унифицированная структура: Впервые предоставляет унифицированный метод вычисления для стабильных и нестабильных LTI систем
• Элегантный трюк с обратной системой: Вычисление rα через инверсию Gα является элегантным решением
• Геометрическая конструкция: Метод пересечения дисков одновременно интуитивен и легко реализуем

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

• Реализуемый алгоритм: Метод основан на стандартных инструментах линейных систем (сингулярные значения, обнаружение минимальной фазы)
• Модульное проектирование: Вычисленный SRG можно переиспользовать для различных проектов коэффициентов усиления
• Интуитивная визуализация: Геометрическое представление на комплексной плоскости облегчает понимание и применение

4. Достаточность примеров

• Охватывает SISO и MIMO, стабильные и нестабильные, системы с интеграторами и другие случаи
• Практическая демонстрация применения к проектированию коэффициентов усиления
• Четкая и эффективная визуализация

Недостатки

1. Теоретические ограничения

• Ограничение на квадратные системы: Основные результаты применимы только к квадратным системам, ограничивая область применения
• Корректность: Не может гарантировать корректность, как расширенный SRG
• Предположение о нормальных матрицах: Для общих MIMO систем граница может быть неточной

2. Практические вызовы

• Зависимость от модели: Требует точной передаточной функции
• Обнаружение неминимальной фазы: Обнаружение из данных является сложной задачей
• Вычислительная сложность: Для высокомерных MIMO систем требуется вычисление для непрерывного множества α

3. Недостаточность экспериментов

• Отсутствие анализа времени вычисления: Не сообщается об эффективности алгоритма
• Отсутствие обсуждения численной точности: Численная ошибка конечного множества α не квантифицирована
• Отсутствие практического сравнения с GNC: Только теоретическое сравнение, отсутствуют практические примеры сравнения

4. Проблемы представления

• Плотность символов: Большое количество математических символов может повлиять на читаемость
• Отсутствие псевдокода алгоритма: Хотя метод ясен, отсутствует описание алгоритма, готовое к прямой реализации
• Объяснение графики: Некоторые графики (например, численные артефакты на рисунке 2b) требуют дополнительного объяснения

Оценка влияния

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

1. Теоретическая полнота: Заполняет пробел в структуре SRG для обработки нестабильных систем
2. Методологический прогресс: Предоставляет новый инструмент для анализа нелинейных систем
3. Модернизация классической теории: Предоставляет современное геометрическое объяснение критерия Найквиста

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

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

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

• Теоретическая воспроизводимость: Все доказательства полностью приведены в приложении
• Реализуемость: Основана на стандартных инструментах, легко реализуется
• Ясность примеров: Предоставлены конкретные передаточные функции для проверки
• Отсутствие кода: Не предоставлена реализация с открытым исходным кодом (это сожаление)

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

Наиболее подходящие применения

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

Менее подходящие сценарии

1. Чистый анализ, управляемый данными: Когда невозможно получить точную модель
2. Неквадратные системы: Когда p ≠ q, метод неполный
3. Приложения реального времени: Если вычислительная сложность является критическим ограничением
4. Системы очень высокой размерности: Когда размерность велика, вычисление может стать сложным

Общая оценка

Это высококачественная теоретическая статья по теории управления, которая вносит значительный вклад в структуру масштабированного относительного графика. Основные преимущества:

• Предоставляет первый систематизированный метод для обработки жесткого SRG нестабильных LTI систем
• Теория строга, доказательства полны
• Устанавливает четкую связь с классическим критерием Найквиста
• Закладывает основу для анализа нелинейных систем

Основные области для улучшения:

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

Рекомендуемый рейтинг: ⭐⭐⭐⭐☆ (4/5)

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

Избранные ссылки

1. 2 Ryu et al., 2022: "Scaled relative graphs: Nonexpansive operators via 2D Euclidean geometry" — исходное математическое определение SRG
2. 3 Chaffey et al., 2023: "Graphical Nonlinear System Analysis" — введение SRG в теорию управления
3. 11 Krebbekx et al., 2025: "Scaled Relative Graph Analysis of General Interconnections" — определение расширенного SRG
4. 12 Chen, 2025: "Soft and Hard Scaled Relative Graphs" — первое определение жесткого SRG
5. 13 MacFarlane & Postlethwaite, 1977: Классическая литература по критерию MIMO Найквиста
6. 14 Desoer & Wang, 1980: Обобщенный критерий Найквиста (GNC)