Быстрое управление, прототипирование, моделирование оптимизационного управления отступным трекингом роя частиц для вращающегося инвертированного маятника

Краткое содержание

Этот протокол устанавливает стандартизированную процедуру быстрого прототипирования управления управлением для оценки оптимизации с оптимизацией Particle Swarm управление отслеживанием назад в условиях симуляции с фиксированным шагом в реальном времени.

Аннотация

Основная цель этого протокола — предоставить воспроизводимую модель симуляции с фиксированным шагом для оценки настройки усиления на основе оптимизации роя частиц (PSO) в нелинейных системах управления. Реализация начинается с формулировки модели маятника типа Фуруты, за которой следует интеграция контроллера обратного шага в среду выполнения с фиксированным шагом длительностью 2 мс. Методология включает систематический четырёхэтапный процесс: характеристику неидеальных ограничений реализации, определение многоцелевого пространства поиска для оптимизации роя частиц, выполнение автоматизированной настройки офлайн и оценку полученных параметров через стандартизированный набор сценариев отслеживания траекторий и отклонения возмущений. Эта система, использующая высокопроизводительные промышленные рабочие станции и стандартизированные интерфейсы сигнала, поддерживает последовательные повторные сравнения в рамках одной и той же архитектуры управления. Дизайн сравнивает базовый вручную настроенный контроллер обратного шага с оптимизированным для PSO вариантом, который использует точно такую же структуру управления, тем самым изолируя влияние выбора усиления. Производительность контроллера оценивается по трём различным операционным сценариям: шаговое следование по траектории, смешанное синусоидальное отслеживание и отклонение возмущений. Статистический анализ 10 повторных испытаний показал, что оптимизация на основе PSO снизила RMSE с 0,065 до 0,050 рад и ослабила пиковые отклонения маятника на 33,1%. Эти улучшения были достигнуты наряду с сокращением усилий по контролю RMS на 22,1%, что свидетельствует о более эффективном распределении энергии в рамках системы Ляпунова. В конечном итоге эта методология предоставляет структурированную модель моделирования для оценки нелинейных стратегий управления перед любым последующим физическим аппаратным внедрением.

Введение

Вращающийся инвертированный маятник типа Фурута представляет собой фундаментальный ориентир для валидации нелинейных алгоритмов управления благодаря своей нестабильности в открытом контуре и сложной недостаточной^{динамике 1,2,3}. Прежде чем рассматривать оценку, ориентированную на реализацию, эти теоретические конструкции требуют строгого промежуточного тестирования на основе моделирования. Таким образом, этот протокол создаёт стандартизированную, быстро управляющую модельную симуляцию прототипов управления для систематической оценки изменений производительности, вызванных оптимизацией роя частиц (PSO) на контроллерах с отслеживанием обратного шага в симулированной среде с фиксированным шагом.

В более широкой литературе многочисленные исследования изучали структурные модификации и настройку параметров для улучшения контроля вращающихся инвертированных маятников. Дизайны на основе оптимизации широко распространены; например, PSO использовался для выбора параметров контроллера и интегрирован в архитектуры управления fuzzy-hybrid ^4,5. По сравнению с другими био-вдохновлёнными метаэвристиками, PSO специально выбран для этой структуры из-за быстрой сходимости в низкомерных непрерывных поисковых пространствах и минимальных требований к гиперпараметрической настройке. Современная литература всё чаще подчёркивает необходимость использования интеллектуальных алгоритмов оптимизации в разнообразных и сложных сценариях управления. Например, продвинутые методы оптимизации эффективно сочетаются с адаптивным управлением на основе моделей (MRAC) и фреймворками дробного порядка для повышения точности отслеживания нелинейных^{сервозаводов 6,7}. Кроме того, оптимизационная настройка доказала свою высокую выгоду для управления связанной динамикой и внутренними ограничениями сложных электромеханических^{систем 8,9}.

Недавние исследования подтверждают, что PSO и его гибридные варианты значительно повышают эффективность отслеживания максимальной мощности в фотоэлектрических решётках, демонстрируя надёжную идентификацию параметров при условиях частичного^{затенения 10}. В робототехнике PSO успешно используется для оптимизации дополненных линейных и нелинейных проектов управления пропорциональными производными для параллельных манипуляторов, минимизируя ошибки отслеживания^{траектории 11}. Кроме того, интеграция PSO с адаптивным управлением скользящим режимом обратного шага оказалась критически важной для подавления вибраций в подвесных массах, приводимых в действие пневматическими искусственными^{мышцами 12}. Помимо фундаментального выбора параметров, интеграция современных стратегий обработки сигналов и надёжной оптимизации критически важна для поддержания стабильности замкнутого цикла в реалистичных, шумных физических^{условиях 13,14}.

Недавние исследования также были сосредоточены на одновременном отслеживании угла суставов и стабилизации маятника в неопределённых условиях с использованием устойчивой обобщённой динамической^{инверсии 15}, а также адаптивной нейронной оценки¹⁶. Кроме того, архитектуры управления шагом назад были значительно эволюционированы для обработки сложных возмущений в различных механических системах, таких как внедрение скользящих режимов для подавления вибраций зданий, использование квази-скользящих наблюдателей для электронных дроссельных клапанов и интеграция нелинейных наблюдателей возмущений для высокоточного регулирования скорости постоянного тока (17,18,19). Эти разнообразные приложения подчёркивают универсальность дизайнов обратного степа в сочетании с надёжными стратегиями оценки или оптимизации.

Критическим недостатком современной литературы является смешение модификаций структурных контроллеров с преимуществами настройки параметров. Многие сравнительные исследования противопоставляют совершенно разные архитектуры управления, что делает невозможным определить, связано ли прирост производительности с фундаментальным алгоритмом или просто с выбором превосходящего^{усиления 20,21}. Кроме того, существующие симуляционные исследования часто предполагают идеальные условия работы и сосредоточены исключительно на базовой стабилизации. Они часто не справляются с ухудшением производительности, вызванным реалистичными ограничениями реализации. Эта методология напрямую устраняет эти пробелы. Вводя имитируемые задержки дискретизации, шум датчиков и несоответствие демпфирования в рамках динамической парадигмы отслеживания траектории, предлагаемый протокол оценивает эффект оптимизации PSO на фиксированную управляющую структуру.

В отличие от традиционного численного интегрирования, этот протокол отличается тем, что эффективность настройки параметров отличается от вариаций структурных контроллеров, одновременно применяя ограничения по времени с фиксированным шагом. Вместо внедрения новой архитектуры управления этот метод оценивает один контроллер отслеживания назад в условиях симуляции с фиксированным шагом в реальном времени. Сравнивая вручную настроенную базовую линию, с оптимизированным для PSO вариантом того же контроллера, протокол разрабатывается так, чтобы приписывать наблюдаемые различия в трекинге шагов, синусоидальном отслеживании и отстранении возмущений процессу усиления-оптимизации. Основной вклад этого исследования — разработка ориентированного на бенчмарк протокола моделирования реального времени, который изолирует влияние выбора усиления на основе PSO от вариаций структурных контроллеров.

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

Протокол

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

1. Строительство станций и создание сигнальной конвенции

1. Постройте высокоточную модель моделирования вращающегося инвертированного маятника типа Фурута. Определим две обобщённые координаты системы: угол вращающегося рычага, θ, и угол отклонения маятника от вертикальной вертикали α.
2. Определим вектор состояния системы как x = [θ, α, θ̇, α̇]T, обозначая α = 0 как нестабильное вертикальное равновесие. Сформулировать математическую модель растения с использованием уравнений Эйлера-Лагранжа для явного определения нелинейной инерционной матрицы, сил Кориолиса и гравитационных векторов, связывающих активный вращательный рычаг и пассивный^{маятник 22}.
3. Определите, что положительные θ и положительные α последовательно используют одинаковую полярность энкодера между установкой, контроллером и скриптами постобработки.
4. Установите решатель моделирования в реальном времени с фиксированным шагом, строго используя номинальные механические и электромеханические параметры, описанные в таблице 1. Настройте систему замкнутого цикла в назначенной среде выполнения в реальном времени, как указано в Таблице материалов. Установите период обновления контроллера на 2 мс и убедитесь, что интервал логирования данных фиксирован на уровне 20 мс.
5. Включить в модель растения три конкретных неидеальных эффекта для имитации реалистичных ограничений реализации. Ввести шум измерения белого при нулевом среднем уровне с отклонениями 0,003 рад и 0,004 рад в каналы обратной связи рычага и маятника соответственно.
6. Введите задержку транспорта 4 мс в измеренный маятниковый сигнал. Одновременно увеличить имитируемые коэффициенты демпфирования установки на 8% по сравнению с номинальными значениями, используемыми в проектировании контроллера, чтобы имитировать несоответствие моделирования.
7. Ограничите выход команд мотора строго ±10 В.
   ВНИМАНИЕ: Настройте программные остановки безопасности для немедленного завершения испытания, если |θ| превышает 0,70 рад или если |α| Превышает 0,35 рад. Кроме того, реализуйте таймер watchdog для запуска сброса системы, если команда управления остаётся на пределе насыщения ±10 В более 100 мс.

2. Начальное условие состояния и предлог

1. Инициализируйте систему в нисходящем состоянии покоя перед каждым тестом.
2. Выполните задачу предпозиционирования, чтобы поднять маятник в перевёрнутое положение и включить стабилизирующий контроллер. Исключить эту фазу предпозиционирования из формальной оценки производительности.
3. Формальный сбор данных начинается только после подтверждения того, что отклонение маятника постоянно оставалось в пределах ±0,05 рад от вертикального равновесия не менее 1,0 с, сохраняя таким образом динамические начальные условия практически сопоставимыми во всех оценённых исследованиях.
4. Сбросьте пробное время до 0.0 с после выполнения этого условия и сразу начинайте выполнение назначенного профиля заданной точки. Убедитесь, что все контроллеры используют одинаковую процедуру предварительного размещения и соответствуют одинаковым стандартам приёмки.

3. Реализация базового контроллера обратного шага

1. Определите ошибку отслеживания eθ(t) согласно командной ссылке на угол руки θr(t) с помощью уравнения
   eθ(t) = θr(t) − θ(t) (1)
2. Постройте составную скользящую поверхность s(t) согласно уравнению 2. В рамках этой формулировки определим λ как строго положительный коэффициент наклона поверхности с отрицательной ошибкой.
   s(t) = ėθ(t) + λeθ(t) (2)
3. Сформулируйте референсный регулятор обратного степа путём интеграции номинального компенсирующего члена на основе модели с пропорциональным членом, стабилизирующим ошибки. Определите переменные ошибки отслеживания математически как z1 = eθ(t) и z2 = s(t).
4. Проведите анализ устойчивости Ляпунова для оценки сходимости ошибок отслеживания. Дифференцировать выбранную функцию Ляпунова V = 1/2z1² + 1/2z2² по времени, чтобы получить V̇ = z1ż1 + z2ż2. Подставим сформулированный управляющий закон в эту производную, чтобы проверить, что V̇ ≤ −k1z1² − k2z2² ≤ 0. Это отрицательно-полуопределённое условие теоретически гарантирует, что ошибки замкнутых систем сходятся в начало координат.
5. Ограничить последующий алгоритм PSO на исследование многомерного пространства параметров в пределах границ, выбранных для сохранения этого условия устойчивости на основе Ляпунова.
6. Замените прерывистую функцию знаков непрерывной функцией насыщения sat(x/φ), чтобы смягчить высокочастотные щебетания в равновесии, следуя стандартной практике обратного шага по^{пограничному слою 23}. Определите φ как ширину пограничного слоя, используемую для сглаживания переходов управления в уравнении 3.
   sat(x/φ) = 1 для x/φ > 1; sat(x/φ) = x/φ для |x/φ| ≤ 1; sat(x/φ) = −1 для x/φ < −1 (3)
7. Оцените угловые скорости по измерениям положения, применяя алгоритм фильтрации обратной разности первого порядка с частотой среза 25 Гц.
8. Применим фиксированные базовые коэффициенты усиления, указанные в таблице 2 (k1 = 3,60, k2 = 1,95, λ = 2,10 и φ = 0,12). Не меняйте эти значения после начала официального тестирования.

4. Оптимизация усиления с помощью возвратного шага с помощью оптимизации роя частиц (PSO)

1. Выполните офлайн-тюнинг перед оценкой в реальном времени. Оптимизировать четыре параметра контроллера (k1, k2, λ и φ) с использованием PSO, настроенного на 20 частиц, и строго ограничить выполнение до 35^{итераций 24}. Завершить оптимизацию на этом пороге, чтобы избежать ненужных вычислительных накладных расходов, поскольку эмпирическая офлайн-настройка показала, что рой стабильно сходится к стабильному глобальному наилучшему значению приспособленности в первые 20 итераций.
2. Инициализировать частицы равномерно в пределах определённых границ поиска: k1 ∈ [2.5, 5.5], k2 ∈ [1.2, 3.0], λ ∈ [1.2, 2.8] и φ ∈ [0.05, 1.20].
3. Оценивайте каждую частицу с помощью 20-секундного офлайн-тренировочного сценария. Убедитесь, что эта среда оценки строго воспроизводит пределы насыщения, шум датчиков, задержки и несоответствия демпфирования, установленные для формальных испытаний в реальном времени.
4. Сформулируйте целевой функцию минимизации как J = 0.55RMSE + 0.25max|α| + 0,20 RMSU.
5. Нормализуйте эти три целевых компонента с помощью фиксированных базовых эталонных значений, вычисленных как среднее пяти начальных базовых запусков, выполненных перед началом поиска роя. Поддерживайте эти нормируемые константы жёстко на протяжении всего процесса оптимизации.
6. Линейно уменьшайте вес инерции с 0,90 до 0,40 в ходе оптимизации. Зафиксируйте оба коэффициента ускорения на уровне 1,50.
7. Назначать заранее определённую штрафную стоимость любой конфигурации частиц, которая запускает остановку безопасности во время её оценки.
8. Извлеките лучшее глобальное решение сразу после завершения финальной итерации и назначите его как оптимизированный контроллер.
9. Реализуйте оптимизированные усиления, описанные в Таблице 2 , для всех формальных PSO-оптимизированных испытаний в реальном времени (k1 = 4,19, k2 = 2,44, λ = 1,83 и φ = 0,92).

5. Конфигурация платформы исполнения в реальном времени

1. Выполнять логику управления и математические модели завода одновременно в среде быстрого прототипирования управления (RCP) с фиксированным шагом в реальном времени. Используйте высокопроизводительную промышленную рабочую станцию с достаточной вычислительной мощностью, чтобы служить целевой машиной в реальном времени. Взаимодействовать алгоритмы управления с имитируемой динамикой завода через виртуальные аналоговые выходные каналы и входные каналы энкодера для имитации ограничений сигнала на уровне реализации.
2. Назначить строгий интервал в 2 мс для размера шага в реальном времени. См. таблицу 3 для подробного обзора аппаратного обеспечения и настроек выполнения.
3. Сбросьте все состояния контроллера перед запуском каждого испытания.
4. Установите эталонное значение на ноль и соберите исходный сегмент данных 2.0 с в вертикальном положении. Проверьте это окно, чтобы убедиться в отсутствии аномальных событий насыщения или дрейфа датчика перед запуском назначенного профиля отслеживания.
5. Регенерируйте реализацию шума измерения с помощью отдельного семя, выбранного из заранее определённого массива семян для каждого нового испытания.

6. Выполнение тестов с отслеживанием шагов

1. Примените последовательность шагов по кусочкам в 20 с, определённую как θr = 0 рад от 0,0 с до 2,0 с, θr = 0,50 рад от 2,0 с до 10,0 с и θr = 0,20 рад от 10,0 с до 20,0 с.
2. Проведите 10 валидных испытаний для базового контроллера и 10 действительных испытаний для контроллера, оптимизированного для PSO.
3. Записывайте прошедшее время, опорный угол, измеренный угол рычага, ошибку отслеживания, отклонение маятника и управляющее напряжение непрерывно в каждом испытании.
4. Оценивайте переходные характеристики непосредственно после входа 0,50 радиационного шага, возникающего при t = 2,0 с.
5. Рассчитайте время роста, определяемое как длительность, необходимая для перехода измеренного отклика угла руки с 10% до 90% от целевого значения.
6. Вычислите перестрелку как процент, на который пиковое угловое отклонение превышает цель 0,50 радиации.
7. Рассчитайте время осаждения, определяемое как первый случай, после которого ответ остаётся ограниченным в пределах ±2% от цели 0,50 рад минимум 1,0 с.
8. Отмечайте любое судебное разбирательство, которое не соответствует этому условию до отметки 10,0 секунд, как нерешённое. Исключите его строго из расчётов среднего времени осядения, сохраняя при этом для всех остальных метрик.
9. Вычислите ошибку в стационарном состоянии, усредняя ошибку отслеживания за временное окно от 9,0 до 10,0 с.

7. Выполнение синусоидального теста отслеживания

1. Введём 20-секундную смешанную синусоидальную эталонную траекторию, регулируемую уравнением 4.
   θr(t) = 0.26sin(2π·0.10t) + 0.12sin(2π·0.30t + 0.40) rad (4)
2. Выполните 10 действительных пробных запусков для каждой конфигурации контроллера.
3. Извлеките RMSE по полному испытанию, максимальную абсолютную ошибку отслеживания, управляющее RMS-напряжение, максимальное абсолютное отклонение маятника и фазовую задержку для каждого отдельного испытания.
4. Оцените фазовую задержку, выполняя перекрестную корреляцию между измеренной траекторией угла руки и командной ссылкой в заданном окне поиска ±0,50 с. Зафиксируйте задержку как положительное значение, когда измеренный выход временно следует за эталонным сигналом.

8. Выполнение теста отклонения возмущения

1. Держите эталонный угол руки строго на уровне 0 радиационного значения на протяжении всех 20 секунд.
2. Введите аддитивный импульс крутящего момента со стороны мотора ровно на 6,0 с и подайте второй идентичный импульс на 12,2 с.
3. Настройте амплитуду возмущения на 0,030 Н·м и установите ширину импульса на 0,12 с для обоих событий инъекции. Выполните 10 действительных пробных процессов на каждый контроллер.
4. Вычислите максимальную абсолютную ошибку отслеживания исключительно в окне 0,50 с сразу после каждого начала импульса. Определите и удержите более крупное из двух значений как репрезентативную пиковая ошибка после возмущения для этого испытания.
5. Вычислите время восстановления, определяемое как прошедшая длительность от начала импульса до первого записанного момента дискретизации, где |eθ| ≤ 0,010 радиация. Не интерполируйте время восстановления субвыборок, так как логирование данных происходит строго с интервалом в 20 мс.
6. Интегрируйте абсолютную ошибку отслеживания и рассчитайте управляющее RMS-напряжение кумулятивно за весь пробный запуск в течение 20 секунд.

9. Экспорт данных и статистическая сумма

1. Экспортируйте сырые наборы данных временной области из тестов шагового, синусоидального и дисмутенционного тестов в отдельные файлы данных.
2. Организуйте метрики повторных испытаний в структурированный матричный формат, присваивая одну строку на каждую пробу и одну выделенную колонку на каждую оцененную метрику.
3. Строго исключайте любые испытания, отмеченные как незавершённые или прерванные застраховочными остановками. Заменить каждый исключённый запуск проведением дополнительного испытания при идентичных условиях, чтобы гарантировать итоговый набор данных, состоящий ровно из 10 действительных испытаний на каждый контроллер для каждого тестового сценария.
4. Сообщите агрегированные результаты повторных испытаний как среднее ± стандартного отклонения.
5. Оцените нормальность распределения извлечённых метрик с помощью теста^{Шапиро-Уилка 25}.
6. Проанализируйте статистические различия между базовыми и оптимизированными контроллерами с помощью двухсторонних независимых t-тестов для нормально распределённых данных. Применим U-тест Манна–Уитни для любых параметров с ненормальными распределениями.
7. Установите порог статистической значимости на p < 0,05.

Результаты

Алгоритм оптимизации роя частиц (PSO) продемонстрировал быстрое начальное сокращение композитной целевой функции, за которым последовал период постепенной сходимости. В частности, оптимальное значение приспособленности снизилось до 1,8757 в первом вычислительном цикле, тогда как среднее значение приспособленности роя снизилось с 2,0573 до 1,2518 к 35-й итерации. Большая часть этой конвергенции происходила в первые 15–20 итераций. После этой стадии траектория лучшего решения по всему миру сходилась, что указывало на то, что рой достиг стабильной области в заранее определённом пространстве поиска, а не блуждал случайно. Окончательный набор усиления, полученный из PSO, сохранял структурную целостность базового закона обратного шага, но перераспределял относительную весовую нагрузку между стабилизирующими членами (Таблица 2). Как показано на рисунке 1, алгоритм оптимизации способствовал быстрому затуханию функции приспособленности, переходя от начальной фазы исследования к стабильному глобальному лучшему решению за 35 итераций. По сравнению с ручной настройкой базовой линии, оптимизированный контроллер принял более крупные значения для k1 и k2, уменьшенный наклон поверхности ошибки λ и более широкий φ пограничного слоя.

Производительность пошагового трекинга

Сравнительный анализ выявляет явные различия в производительности между двумя контроллерами во время эксперимента с трекингом шагов, иллюстрируя ожидаемые положительные и неоптимальные результаты этого протокола. Базовый контроль обратного шага представляет собой типичный субоптимальный результат; Хотя она сохраняла общую стабильность после команды 0,50 радика, введённой на 2,0 с, она демонстрировала быстрый начальный рост, преждевременно отклоняющийся от целевой траектории, что приводило к стойким остаточным колебаниям и требовало длительного времени для сходов. В то же время контроллер, оптимизированный для PSO, демонстрирует репрезентативный положительный результат, характеризующийся быстрым, затухающим подходом к заданной точке без срабатывания насыщения команд. В течение десяти последовательных испытаний оптимизированный контроллер снизил RMSE с 0,065 ± 0,002 рад до 0,050 ± 0,002 рад (тест Welch t, p < 0,001, Таблица 4). Временные метрики также улучшились: перераспределение снизилось с 4,644 ± 0,604% до 3,627 ± 0,620% (Mann–Whitney U, p = 0,009), а время роста было сокращено с 0,788 ± 0,025 с до 0,548 ± 0,025 с (Mann–Whitney U, p < 0,001). Эти различия в производительности показаны на рисунках 2 и 3, где оптимизированный отклик демонстрирует улучшенное демпфирование и более быструю сходимость по заданным точкам. Эта повышенная точность была достигнута параллельно снижением управляющего напряжения по среднеквадратичному корню на 22,1% (с 0,669 ± 0,004 В до 0,521 ± 0,003 В). Одновременно максимальное колебание маятника было снижено на 33,1%, снизившись с 5,812° до 3,889° (Таблица 4). Кроме того, хотя ни одно базовое испытание не удовлетворяло критерию осаждения 1,0 с в пределах окна 10,0 с, оптимизированный контроллер добился успешного осаждения в 70% испытаний, зафиксировав среднее время осаждения 1,057 ± 0,348 с.

Производительность синусоидального отслеживания

Обе стратегии контроля поддерживали ограниченное отслеживание в межсмешанной частоте синусоидального референса на протяжении всего периода 20-х тестов. Хотя визуально разрыв в производительности был менее заметен, чем в тестах с пошаговым ответом, статистические результаты повторных исследований показали улучшение результатов при оптимизированном подходе. Базовый управляющий сигнал имел более широкую огибающую ошибок отслеживания и заметно отставал от целевого сигнала. В отличие от этого, оптимизированный контроллер отслеживал заданную точку более точно как при низкочастотных, так и при высокочастотных вариациях, обеспечивая более точное распределение пиковых ошибок вдоль эталонной траектории. В результате RMSE в полном испытании был снижен с 0,061 ± 0,003 рад до 0,040 ± 0,003 рад после оптимизации. Как указано в таблице 5, оптимизированный контроллер также способствовал снижению управляющего напряжения RMS (с 0,557 ± 0,002 В до 0,402 ± 0,002 В) и ограничил максимальное перемещение маятника с 2,156 ± 0,079° до 1,464 ± 0,061°, что указывает на улучшение фазовой синхронизации и энергоэффективности (p < 0,001, Таблица 5; Рисунок 4).

Отторжение нарушений и устойчивость между сценариями

Введение аддитивных крутящих импульсов на 6,0 и 12,2 с вызывало большие отклонения угла рычага и более широкие отклонения маятника в отклике базового контроллера. Хотя оптимизированный контроллер также сразу фиксировал возмущения, его последующие отклонения были ниже. Процесс восстановления к равновесию с нулевой опорой показал качественно более быстрый тренд (от 0,014 с до 0,000 с); Однако этот конкретный показатель не достиг статистической значимости (p = 0,078), главным образом потому, что разрешение логирования данных в 20 мс вызывало измеримый эффект пола. В ходе многократных попыток оптимизированная параметризация снижала пиковую ошибку после возмущения. Интеграл абсолютной ошибки (IAE), вычисленный за 20 с, был уменьшен с 0,377 ± 0,012 до 0,248 ± 0,015, а также соответствующее снижение управляющего RMS-напряжения с 0,656 ± 0,003 В до 0,508 ± 0,005 В. Максимальная амплитуда колебаний маятника была аналогично ослаблена с 11,456 ± 0,117° до 7,561 ± 0,063° (Таблица 6). Рейтинг контролёров оставался стабильным во всех трёх тестовых профилях. Стратегия обратного движения, оптимизированная для PSO, приводила к снижению временных ошибок, снижению неточностей в стационарном состоянии и минимизации пассивных отклонений маятника (стр. < 0.001, Таблица 6; Рисунок 5) без необходимости увеличения общей мощности привода. Распределение основных показателей эффективности по межопытным испытаниям изложено на рисунке 6.

ДОСТУПНОСТЬ ДАННЫХ:

Поддерживающий набор данных, используемый для генерации фигур и таблиц, включая репрезентативные данные, метрики на уровне исследования, заранее заданные случайные семена и записи настройки PSO, был депонирован в Zenodo и доступен публично под именем DOI 10.5281/zenodo.20395966.

Рисунок 1: Оптимизационная сходимость роя частиц при настройке контроллеров. (A) Лучшие значения приспособленности роя за 35 итераций; (B) Средние значения приспособленности для каждой временной точки. Композитный объект быстро распадался на начальной фазе разведки, но постепенно замедлялся, когда сходился к стабильной области оптимизированного для PSO закона управления обратным шагом. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 2: Репрезентативные ответы пошагового отслеживания базового и оптимизированного для PSO контроллера обратного степпинга. (A) Командный отсчёт по углу руки и измеренный отклик угла руки для полной последовательности шагов длиной 20 секунд; (B) Расширенный вид реакции сразу после шага 0,50 рад, приложенного на 2,0 с; (C) Расширенный вид отклика после начала второго сегмента команды 0,20 радиации на 10,0 с. После оптимизации контроллер, оптимизированный для PSO, показал более быстрое время роста, меньшую пиковую ошибку и лучшее демпфирование по сравнению с базовым контроллером. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 3: Ошибка отслеживания по углу рычага и поведение управления под двумя контроллерами. (A) Ошибка отслеживания угла руки во время полного теста пошагового отслеживания; (B) Сигнал управляющего напряжения для того же испытания; (C) Угол отклонения маятника во время шагового отклика. Оптимизированный контроллер показывал меньшие ошибки, более узкий диапазон колебаний и относительно более плавное движение при более низком уровне команды RMS. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 4: Репрезентативная синусоидальная система отслеживания при смешанном частотном ориентире. (A) Командовал синусоидальными эталонными сигналами и измерял реакции угла руки в рамках эксперимента с отслеживанием 20 секунд; (B) ошибка отслеживания в том же временном диапазоне; (C) Угол отклонения маятника в процессе синусоидального отслеживания. Оптимизированный PSO-контроллер отслеживал композитную команду более точно, чем базовый контроллер, демонстрируя снижение фазовой задержки и уменьшенный диапазон ошибок в последовательных циклах. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 5: Репрезентативные реакции на нарушение возмущений при нулевой регуляции. (A) ошибка отслеживания по углу рычага во время теста на отстранение возмущения, включая нарушения крутящего момента и импульса на 6,0 и 12,2 с; (B) Угол отклонения маятника во время того же испытания; (C) Управляющий сигнал напряжения при отстранении возмущения. По сравнению с базовым контроллером, оптимизированное для PSO управление показало меньшее максимальное отклонение после возмущения и более быстрое восстановление к заданной точке. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 6: Межпробное распределение основных показателей эффективности во всех трёх ситуациях тестирования. (A) RMSE с отслеживанием шагов; (B) Время роста с отслеживанием шагов; (C) Синусоидальное отслеживание RMSE; (D) Синусоидальная фазовая задержка; (E) Пиковая ошибка отслеживания после возмущения; (F) Максимальное смещение маятника при условии возмущения-отклонения. Каждая панель показывает распределение 10 действительных испытаний при каждом состоянии контроллера. Центральная тенденция и дисперсия базовых и оптимизированных для PSO систем управления шагом назад сравниваются напрямую. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Параметр	Символ	Ценность	Подразделение
Длина вращающегося рычага	Lr	0.215	m
Длина центра масс маятника	Лп	0.168	m
Масса вращающегося рычага	МР	0.254	кг
Маятниковая масса	mp	0.097	кг
Вязкое демпфирование рычага	br	0.0031	N·m·s/rad
Вязкое демпфирование маятника	bp	0.0018	N·m·s/rad
Постоянный крутящий момент двигателя	Кт	0.053	N·m/A
Разрешение энкодера	Nenc	4096	Counts/rev
Период обновления контроллера	Тс	0.002	s

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

Контроллер	Усиление	Ценность	Описание
Базовый шаг назад	k1	3.6	Ручная настройка начального стабилизационного усиления
Базовый шаг назад	k2	1.95	Ручная настройка виртуального усиления управления
Базовый шаг назад	λ	2.1	Наклон поверхности с ошибкой
Базовый шаг назад	φ	0.12	Ширина пограничного слоя
Оптимизированный для PSO шаг назад	k1	4.19	Лучшая частица в мире после 35 итераций
Оптимизированный для PSO шаг назад	k2	2.44	Лучшая частица в мире после 35 итераций
Оптимизированный для PSO шаг назад	λ	1.83	Лучшая частица в мире после 35 итераций
Оптимизированный для PSO шаг назад	φ	0.92	Лучшая частица в мире после 35 итераций

Таблица 2: Базовые усиления контроллера с обратного шага и оптимизированные для PSO усиления контроллера обратного шага .В таблице представлены комбинации усиления для рабочих состояний обоих контроллеров. Базовые значения фиксировались при ручной настройке; Оптимизированные по PSO приросты относятся к итоговой глобальной лучшей частице, сохраняемой после 35 итераций оптимизации.

Предмет	Технические характеристики
Компьютер целей в реальном времени	Промышленный ПК, процессор класса Intel i7, 16 ГБ оперативной памяти
Режим выполнения	Симуляция с фиксированным шагом в реальном времени
Период обновления контроллера	2 мс
Интервал экспорта сигнала	20 мс
Интерфейс сигнала	Вход и аналоговый выход энкодера ±10 В
Конфигурация завода	Нелинейная модель вращающегося маятника типа Фурута
Введённые неидеальные эффекты	Шум измерения, задержка маятникового канала и несоответствие демпфирования
Шум измерения	0,003 радиационного угла рычага; 0,004 радиация на углу маятника
Задержка с маятниковым каналом	4 мс
Несоответствие демпфирования заводов	+8% относительно номинальной модели контроллера
Насыщение команд	±10 В
Критерии безопасности остановки	|θ| > 0,70 радиация, |α| > 0,35 радиация, или насыщение > 100 мс
Предварительное принятие позиционирования	|α| ≤ 0,05 радиации минимум 1,0 с до начала испытаний
Профиль опоры шагов	0 рад (0,0-2,0 с), 0,50 рад (2,0-10,0 с), 0,20 рад (10,0-20,0 с)
Синусоидальный эталонный профиль	0,26 sin(2π·0,20t) + 0,12 sin(2π·0,30t + 0,40) rad
Протокол возмущения	Два аддитивных импульса крутящего момента на 6,0 с и 12,2 с
Амплитуда и ширина возмущения	0,030 Н·м, 0,12 с
Действительные испытания по контроллеру и сценарию	10

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

Метрика	Базовый шаг назад n=10	Оптимизированный PSO-шаг назад n=10	Значение p	Статистический тест
Отслеживание RMSE (радиостанция)	0,065 ± 0,002	0,050 ± 0,002	<0.001	Тест Уэлча t-тест
Ошибка устойчивого состояния, 9,0–10,0 с (радиация)	-0,007 ± 0,001	-0,004 ± 0,001	<0.001	Университет Манн–Уитни
Пересечение (%)	4,644 ± 0,604	3.627 ± 0.620	0.009	Университет Манн–Уитни
Время подъёма	0,788 ± 0,025	0,548 ± 0,025	<0.001	Университет Манн–Уитни
Время урегулирования	Нет действительных урегулированных судебных разбирательств (0/10)	1.057 ± 0.348 (7 из 10 действительных урегулированных судебных процессов)	NA	Не сравнивать
Управляющий RMS (V)	0,669 ± 0,004	0,521 ± 0,003	<0.001	Тест Уэлча t-тест
Максимальное отклонение маятника (°)	5.812 ± 0.156	3,889 ± 0,109	<0.001	Тест Уэлча t-тест

Таблица 4: Метрики производительности с отслеживанием шагов по 10 действительным испытаниям на один контроллер. В этой таблице представлена статистика повторных испытаний эксперимента с трекингом шагов, включающая отслеживание RMSE, ошибку в стационарном состоянии, переброс, время роста, скорость сходимости заданной точки, управляемое RMS-напряжение и максимальное отклонение маятника. Значения указываются как средние ± стандартных отклонениях, если не указано иное.

Оптимизированный PSO-шаг назад n=10	Значение p	Статистический тест
0,040 ± 0,003	<0.001	Тест Уэлча t-тест
0,088 ± 0,006	<0.001	Тест Уэлча t-тест
0,232 ± 0,019	<0.001	Тест Уэлча t-тест
0,402 ± 0,002	<0.001	Тест Уэлча t-тест
1.464 ± 0.061	<0.001	Тест Уэлча t-тест

Таблица 5: Показатели производительности синусоидального отслеживания по 10 действительным испытаниям на один контроллер. В этой таблице представлены резюме повторных испытаний теста синусоидального отслеживания, которые включают полную пробную ошибку среднего квадрата корня (RMSE), максимальную абсолютную ошибку отслеживания, фазовую задержку и управляющее RMS-напряжение. Также приведено наибольшее отклонение маятника. Положительная фазовая задержка указывает на то, что измеренный выход отстаёт от заданной точки.

Значение p	Статистический тест
<0.001	Тест Уэлча t-тест
0.078	Университет Манн–Уитни
<0.001	Тест Уэлча t-тест
<0.001	Тест Уэлча t-тест
<0.001	Тест Уэлча t-тест

Таблица 6: Метрики производительности отклонения помех по 10 действительным испытаниям на один контроллер. Данные повторных испытаний экспериментов по отторжению возмущений включают пиковую ошибку отслеживания после возмущения, время восстановления, интеграл абсолютной ошибки (IAE), управляющее напряжение RMS и максимальное отклонение маятника. Согласно настройкам экспорта сигнала, измеримое время восстановления ограничено периодом логирования в 20 мс.

Обсуждение

Основная цель данного исследования — не просто объявить один контроллер превосходящим другой в изолированных условиях, но и показать, что оцененная здесь структура оптимизации усиления на основе PSO может повысить производительность контроллера обратного шага в определённых симулированных неидеальных сценариях. Вращающийся инвертированный маятник служит отличным ориентиром для этой оценки благодаря своему высоко нелинейному, неминимальному фазовому и недостаточно активируемым^{характеристикам 26}. Помимо метрик производительности, техническая целостность этого протокола зависит от целенаправленной характеристики среды симуляции в реальном времени. Критический этап реализации включает контролируемое введение 8% несоответствия демпфирования и задержки датчика в 4 мс, поскольку эти специфические неидеальные эффекты позволяют алгоритму PSO оптимизировать усиления при тех же возмущениях, что и при формальной оценке, а не только при идеализированных условиях.

Тестируя стратегии управления с помощью протоколов пошагового трекинга, синусоидального отслеживания и отклонения возмущений в единой установке фиксированного шага выполнения, результаты симуляции показали стабильную динамику производительности. В частности, улучшения, наблюдаемые как в ступенчатой, так и в синусоидальной реакции — в частности, снижение ошибок отслеживания, минимизация пересечения и снижение фазовой задержки — были достигнуты вместе с заметным снижением управляющего напряжения RMS. Это указывает на то, что процесс оптимизации не просто повысил точность отслеживания за счёт насыщения приводного блока; вместо этого он перераспределял усилия по контролю в существующей системе обратного степа на основе Ляпунова. Согласованность этого перераспределения параметров существенно зависит от критериев принятия предпозиционирования, определённых в протоколе Step 2.4. Строго следить за тем, чтобы маятник оставался в пределах определённого окна равновесия в течение длительного времени; В противном случае непредсказуемые начальные переходы могут привести к преждевременной застойной функции затрат или неконвергентному поведению роя в последующем офлайн-режиме. Для недостаточно активируемых механических систем это перераспределение важно, поскольку оно определяет, как распределяются усилия управления для управления связью между активной координатой (вращающийся рычаг) и пассивной координатой (маятник)²⁷.

Кроме того, одновременное снижение ошибки отслеживания по углу рычага и внутреннего смещения маятника подчёркивает преимущество предлагаемого протокола настройки. В управлении недостаточно активируемыми системами достижение точного отслеживания координат за счёт волатильных внутренних состояний обычно считается неоптимальной или даже неудачной конструкцией. Оптимизированный закон управления снизил амплитуду реакции неактивированного маятника на 33,1%, одновременно улучшая точность отслеживания. Это совокупное улучшение свидетельствует о том, что алгоритм PSO определил комбинацию усиления, которая лучше сбалансировала нелинейную связь между активными и пассивными координатами, а не полагалась на обратную связь с высоким усилением, способную усиливать паразитные колебания. Такая сбалансированная эффективность является более информативным показателем эффективности управления, чем чистое отслеживание траектории, поскольку указывает на то, что динамика внутреннего состояния была адекватно ограничена в процессах преобразования энергии²⁸. Аналогично, результаты отклонения возмущений подтверждают практическую устойчивость оптимизированных усилений в условиях смоделированного теста. Хотя интервал логирования в 20 мс вызывал эффект пола, который заглушал точное время восстановления субпроб, одновременные снижения пиковой ошибки после возмущения, накопленной абсолютной ошибки и максимального смещения маятника указывают на то, что оптимизированные параметры ограничивали степень распространения внешних импульсных инъекций по системной^{динамике 29}.

Успешное выполнение этого протокола требует устранения распространённых узких мест реализации. Например, если траектория роя не стабилизируется в течение 35 итераций, исследователям следует сначала проверить, достаточно ли границы поиска, указанные в протоколе Step 4.2, охватывают стабильную область Ляпунова. Кроме того, если целевая машина в реальном времени часто срабатывает на этапе отклонения возмущений, сокращение периода обновления контроллера, если это поддерживается платформой, или уточнение частоты среза фильтра обратной разницы могут помочь сохранить целостность замкнутого контура без ущерба логике оптимизации.

Несмотря на эти многообещающие результаты, методология остаётся ограниченной абстракциями, присущими средам симуляции. Реализованный протокол учитывал задержки дискретизации, насыщение команд, несоответствие демпфирования и шум измерения — факторы, которые часто игнорируются при чисто офлайн-численной интеграции. Однако настоящее исследование не включало физические эксперименты с вращающимся инвертированным маятником, физическую валидацию аппаратного обеспечения в цикле или тестирование физического развертывания. Поэтому результаты следует интерпретировать как доказательства на основе моделирования эффекта производительности настройки усиления на основе PSO в условиях фиксированного шага и неидеального моделирования. Будущие работы могут расширить смоделированный набор неопределённости, включив дополнительные эффекты трения, актуатора, соответствия и квантования в симуляции до того, как будет рассматриваться отдельное аппаратно-ориентированное^{исследование} 30,31,32.

В данном исследовании устанавливается протокол моделирования прототипирования быстрого управления для оценки оптимизации усиления на основе PSO для управления обратной скоростью. Количественные результаты 10 повторных испытаний показывают, что оптимизированный контроллер снижает RMSE с пошагового отслеживания с 0,065 до 0,050 rad и достигает снижения усилий по контролю RMS на 22,1% по сравнению с ручной настройкой базового уровня. Затухание в 33,1% при смещениях маятника указывает на улучшенную внутреннюю стабильность при смоделированных неидеальных условиях. Хотя эта среда моделирования включает шум, задержки, несоответствие демпфирования и насыщенность, она не полностью воспроизводит сложные кулоновские трения, мертвые полосы двигателя, структурную совместимость или квантование датчиков, характерные для физических систем. Будущие работы будут сосредоточены на расширении модели моделирования с включением дополнительных неидеальных эффектов и на проверке, остаётся ли тот же ранжировочный паттерн стабильным при более широких моделируемых условиях неопределённости.

Материалы

Список материалов, использованных в этой статье

Имя	Компания	Каталожный номер	Комментарии
Высокопроизводительная рабочая станция (Windows 11 Pro)	Различные / Настраиваемый вариант	N/A	Хост-машина для симуляции в реальном времени с фиксированным шагом и последующей обработки.
MATLAB (Версия R2024a)	MathWorks	https://www.mathworks.com/products/matlab.html	Моделирование нелинейной установки, кодирование контроллера, экспорт данных и управление параметрами.
Simulink (Версия R2024a)	MathWorks	https://www.mathworks.com/products/simulink.html	Конструирование блочно-схемной модели для маятника типа Фурута и выполнение контроллера.
Simulink Desktop Real-Time (Версия R2024a)	MathWorks	https://www.mathworks.com/products/simulink-desktop-real-time.html	Ядро реального времени с фиксированным шагом для выполнения модели управления на настольном компьютере.
Python (Версия 3.11)	Python Software Foundation	https://www.python.org/	Вторичная обработка данных, статистическая обработка и подготовка рисунков.
NumPy (Версия 1.26)	Разработчики NumPy	https://numpy.org/	Операции с числовыми массивами для экспортированных данных испытаний.
pandas (Версия 2.2)	Разработчики pandas	https://pandas.pydata.org/	Организация данных повторных испытаний и создание сводных таблиц.
SciPy (Версия 1.13)	Разработчики SciPy	https://scipy.org/	Статистическое тестирование и утилиты анализа сигналов.
Matplotlib (Версия 3.8)	Разработчики Matplotlib	https://matplotlib.org/	Создание графиков для сходимости, отслеживания и распределения.