Моделирование нелинейной системы управления мобильными объектами в среде Matlab
Вид материала | Документы |
СодержаниеNCD Blockse |
- П. П. Порешин московский инженерно-физический институт (государственный университет), 23.75kb.
- Моделирование системы управления стеклоочистителем на базе нечеткого контроллера, 16.79kb.
- Математическое моделирование системы управления для управляемого выпрямителя на igbt-транзисторах, 26.76kb.
- Моделирование автоматизированной системы контроля и управления технологическими объектами, 35.31kb.
- Программа спецкурса "Компьютерное моделирование нелинейных волновых процессов" Специальность, 27.11kb.
- Математическая модель тиристорного асинхронного электропривода в среде matlab, 183.43kb.
- Типография Издательства Мордовского университета 430000, Саранск, ул. Советская,, 365.09kb.
- Тема диссертационной работы, 45.86kb.
- Лабораторная работа на тему "Решение оду и систем оду в среде Simulink", 27.38kb.
- При выполнении сложных расчетных заданий в курсе теории автоматического управления, 83.98kb.
Моделирование нелинейной системы управления
мобильными объектами в среде Matlab
Шарипбаев А.А. д.т.н., профессор, заведующий кафедрой "Вычислительная техника" ЕНУ им Л.Н.Гумилева,
Атанов С.К., к.т.н, доцент кафедры «Компьютерные системы» КазАТУ им. С.Сейфуллина
Управление мобильными объектами с помощью любой инерциальной системы может рассматриваться как взаимодействие двух процессов: решения навигационной задачи и решения задачи стабилизации. Рассмотрим эту проблему на примере летательного аппарата, где данные задачи наиболее актуальны в связи с возможностью их автономной работы и высокими скоростями движения. Первая задача заключается в определении требуемой траектории летательного аппарата и в вычислении фактической траектории, вторая — в управлении аппаратом для поддержания требуемого курса с заданной точностью [1]. Так на рис 1. приведена типичная схема бесплатформенной системы управления космическим аппаратом для выхода на геостационарную орбиту. Подобный механизм используется для управления крылатыми ракетами. С бортовой цифровой вычислительной машиной 1 соединены три группы датчиков, условно обозначенных через Д1 , Д2 и Д3.
![](images/239557-nomer-6eaab877.png)
Рис. 1 - Схема бесплатформенной системы управления ориентацией:
1 – бортовая цифровая вычислительная машина;
2 – блок согласования;
3-исполнительные органы
Вырабатываемые в машине сигналы управления преобразуются должным образом в блоке согласования 2, после чего поступают на исполнительные органы системы ориентации 3. Воздействуя на динамику ракеты (в зависимости от работы исполнительных органов, изменяется его угловое движение, и на входе вычислительной машины появляются измененные сигналы датчиков системы ориентации). На приведенной схеме все датчики условно разбиты на три группы в зависимости от основной задачи, выполняемой ими в полете.
При этом задача системы стабилизации — обеспечить управление рулями и тягой таким образом, чтобы выполнить задаваемую программу полета с требуемой точностью. Ее структурная схема показана на рис. 2.
![](images/239557-nomer-m79e9a2d2.gif)
Рис 2.Структурная схема системы стабилизации ракеты на курсе
Линейная математическая модель, описывающая движение аппарата, имеет вид[2]:
![](images/239557-nomer-91b2ccb.gif)
где
![](images/239557-nomer-3861e54c.gif)
![](images/239557-nomer-m6023856d.gif)
![](images/239557-nomer-25883904.gif)
![](images/239557-nomer-m37b6fad4.gif)
![](images/239557-nomer-da2014c.gif)
![](images/239557-nomer-m31fc7245.gif)
Линейная модель привода (рулевой машины) представляет собой интегрирующее звено с передаточной функцией
![](images/239557-nomer-6ddf5d7.gif)
охваченное единичной отрицательной обратной связью. На угол перекладки руля и скорость перекладки накладываются нелинейные ограничения
![](images/239557-nomer-5364ab08.gif)
![](images/239557-nomer-1c6b9310.gif)
Для измерения угла рыскания используется GPRS-навигатор, математическая модель которого записывается в виде апериодического звена первого порядка с передаточной функцией [2]:
![](images/239557-nomer-m3d501120.gif)
В качестве управляющего устройства рассмотрим ПИД-регулятор с передаточной функцией для конкретного примера из[3]:
![](images/239557-nomer-5159a9f9.gif)
![](images/239557-nomer-34761563.gif)
![](images/239557-nomer-m7d9249cf.gif)
Для компенсации эффекта насыщения, вызванного ограниченным углом перекладки руля, используем схему с внутренней нелинейной обратной связью, охватывающей интегратор в составе регулятора.
![](images/239557-nomer-m3248e840.gif)
Базовый регулятор, построенный по линейной модели, выделен серым фоном. Сигнал
![](images/239557-nomer-mcbbf10b.gif)
Если насыщения нет, разность сигналов
![](images/239557-nomer-m8540632.gif)
![](images/239557-nomer-mcbbf10b.gif)
![](images/239557-nomer-m8540632.gif)
Для выбора оптимального значение коэффициента
![](images/239557-nomer-3ce8c238.gif)
Движение судна описывается линейной математической моделью в виде передаточной функции и для конкретного частного случая примет вид:
![](images/239557-nomer-5522b930.gif)
![](images/239557-nomer-m3bc5d69e.gif)
![](images/239557-nomer-67fb8f5d.gif)
Линейная модель привода представляет собой интегрирующее звено с передаточной функцией
![](images/239557-nomer-6ddf5d7.gif)
![](images/239557-nomer-m33641ee1.gif)
охваченное единичной отрицательной обратной связью. На угол перекладки руля и скорость перекладки наложим нелинейные ограничения
![](images/239557-nomer-5364ab08.gif)
![](images/239557-nomer-1c6b9310.gif)
Измерительное устройство (гирокомпас, GPRS-навигатор или иное) моделируется как апериодическое звено с передаточной функцией
![](images/239557-nomer-m3d501120.gif)
![](images/239557-nomer-m11bf1bdd.gif)
В качестве управляющего устройства используется ПИД-регулятор с передаточной функцией
![](images/239557-nomer-5159a9f9.gif)
где
![](images/239557-nomer-11f3821e.gif)
![](images/239557-nomer-m2538a921.gif)
![](images/239557-nomer-34761563.gif)
![](images/239557-nomer-m437bbb56.gif)
Для компенсации эффекта насыщения, вызванного ограниченным углом перекладки руля, используем схему с нелинейной обратной связью, охватывающей интегратор в составе регулятора. В ходе моделирования рассчитаем оптимальное значение коэффициента усиления в обратной связи с помощью пакета NCD Blockset [3].
В этом случае блок компенсации насыщения будет иметь вид:
- подсистема «Регулятор»
![](images/239557-nomer-52f4ff7.png)
- пределы для блока Saturation
- номинальное значения
- модель для сравнения трех типов систем
![](images/239557-nomer-2a0b9d1f.png)
- переходные процессы при
![](images/239557-nomer-m7d9d2736.gif)
При введении компенсации насыщения интегрирующего звена видим устойчивое и оптимальное управление при заданном изменении курса.
Проведем анализ оптимальной системы управления
![](images/239557-nomer-3ce8c238.gif)
![](images/239557-nomer-m4c3397ea.png)
- переходные процессы при
(поворот на 90 градусов)
![](images/239557-nomer-29040292.gif)
- перерегулирование
, время переходного процесса
сек.
- переходные процессы при повороте на 30 градусов
![](images/239557-nomer-m7651ab2c.gif)
- перерегулирование
, время переходного процесса
сек.
Таким образом, применение оптимальной компенсации позволило выйти на заданный курс полета на 20-40 % быстрее линейной системы и нелинейной системой управления даже в экстремальных условиях полета, а также избежать перерегулирования, что могло на высоких скоростях полета привести к турбулентности и потери управляемости ракеты.
Литература
1. Киреев Н.Г. Системы управления беспилотных летательных аппаратов. – К.: УМК ВО, 1993. – 160 с.
2. Пельпор Д.С. Гироскопические системы ориентации и стабилизации. Справочное пособие. - М.: Машиностроение, 1982, 165 с.
3. Потемкин В. Г. MATLAB 6: среда проектирования инженерных приложений.—
М.:Диалог–МИФИ, 2003.— 448с.