Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
_MU_LR_2НАДО.doc
Скачиваний:
1
Добавлен:
01.07.2025
Размер:
571 Кб
Скачать
☆

3. Объекты и средства исследования

В работе рассматривается система управления летательным аппаратом (ЛА) по курсу. Ее структурная схема показана на рисунке.

Линейная математическая модель, описывающая рыскание ЛА, имеет вид

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

.

Структурная схема системы стабилизации ЛА на курсе

Линейная модель привода (рулевой машины) представляет собой интегрирующее звено с передаточной функцией

,

охваченное единичной отрицательной обратной связью. На угол перекладки руля и скорость перекладки накладываются нелинейные ограничения

, .

Для измерения угла рыскания используется гирокомпас, математическая модель которого записывается в виде апериодического звена первого порядка с передаточной функцией1

,

В качестве управляющего устройства используется ПИД-регулятор с передаточной функцией

, где сек и сек.

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

Базовый регулятор, построенный по линейной модели, выделен серым фоном. Сигнал на его выходе представляет собой желаемый угол перекладки руля. Для блока нелинейной коррекции типа «насыщение» выбираются пределы, равные ограничениям на угол перекладки руля.

Если насыщения нет, разность сигналов равна нулю, и обратная связь не работает; используется закон управления, синтезированный для линейной системы. Если сигнал превышает допустимые пределы, разность подается (со знаком «минус») на вход интегратора через усилитель. Таким образом, при насыщении сигнал на входе интегратора ослабляется тем сильнее, чем больше разность между желаемым и допустимым углами перекладки. Такой метод коррекции получил в литературе название anti-windup (противодействие «наматыванию»).

В ходе выполнения работы для выбора оптимального значение коэффициента применяются процедуры численной оптимизации пакета NCD Blockset.

  1. Инструкция по выполнению работы

Этап выполнения задания

Команды Matlab

Сделайте свою папку рабочей папкой Matlab.

ЛКМ по кнопке справа от поля Current Directory

Создайте новый М-файл, в первой строчке наберите имя скрипта sysdata для загрузки исходных данных. Сохраните файл с именем lab6go.m и выполните его.

File – New – M-file

File – Save

клавиша F5

Откройте окно рабочей папки.

View – Current directory

Откройте модель, построенную в лабораторной работе № 5.

двойной щелчок на lab5.mdl

Сохраните модель в своей папке под именем lab6.mdl.

File – Save as ...

Откройте блок Нелинейная модель-Регулятор. Удалите линию, соединяющую вход блока с входом И-канала.

выделить ЛКМ, нажать Delete

Скопируйте сумматор, расположите его перед интегратором и измените знак у второго входа на «минус».

перетащить ПКМ

|+- в поле List of signs

Скопируйте сумматор еще раз, разверните его на 180 градусов и измените расположение входов.

перетащить ПКМ

нажать Ctrl+I

+-| в поле List of signs

Откройте блок Нелинейная модель-Привод и перетащите в окно Регулятор второй нелинейный блок (ограничение на угол перекладки руля). Разверните его на 180 градусов.

перетащить ЛКМ

нажать Ctrl+I

Добавьте блок Gain из группы Math Operations и установите коэффициент усиления, равный Kaw.

Kaw в поле Gain

Соберите схему, показанную на рисунке. Скопируйте ее в отчет.

Добавьте в скрипт lab6go команду, устанавливающую начальное значение для Kaw. Сохраните скрипт и выполните его.

Kaw = 1;

Далее мы сравним переходные процессы в этой системе с результатами моделирования линейной и исходной нелинейной систем. Закройте окна всех подсистем, кроме главного окна модели. Дайте блоку название Система с компенсацией и поменяйте цвет фона на светло-голубой (этим цветом будут изображаться на осциллографах соответствующие графики).

Скопируйте в новую модель подсистему Линейная модель из модели lab4.mdl и подсистему Нелинейная модель из модели lab5.mdl.

открыть нужную модель, перетаскивание ЛКМ

Добавьте в модель мультиплексор (блок Mux из группы Signal Routing) и установите для него количество входов 3.

двойной щелчок ЛКМ

3 в поле Number of inputs

Создайте копию мультиплексора и соберите схему по рисунку так, чтобы на каждом из осциллографов изображались по 3 сигнала с разных моделей. Входы f всех систем надо соединить с выходом блока Возмущение, а входы r – с выходом блока Заданный курс.

Перейдите в окно скрипта lab6go в редакторе. Добавьте строчки для запуска процедуры моделирования и для вывода графика, учитывая, что на каждой половине рисунка надо выводить по 3 графика, массивы phi и delta содержат по 4 столбца. Удобно скопировать и немного изменить команды для оформления графика из файла lab4graph.m.

sim ('lab6' )

figure(1)

...

Запустите скрипт. После отладки сохраните скрипт на диске и скопируйте его в отчет. Скопируйте в отчет построенный график. Опишите полученный эффект.

Добавьте в основное окно модели усилитель (блок Gain из группы Math Operations) и установите коэффициент усиления, равный 1/phiZad. Этот (нормирующий) блок будет служить для того, чтобы на выходе установившееся значение было равно 1 (нормировка).

View – Library Browser

перетаскивание ЛКМ

двойной щелчок по блоку

1/phiZad в поле Gain

Добавьте в модель блок NCD Outport из пакета NCD Blockset и соедините его через нормирующий усилитель (см. выше) с выходом phi системы с компенсацией насыщения. Сохраните схему в отчете.

Откройте блок NCD Outport. Цель оптимизации – выбрать некоторые коэффициенты системы так, чтобы переходный процесс оставался внутри черной зоны графика. Горизонтальные и вертикальные границы ограничений (красные полосы) можно сдвигать мышкой. Полоса выделяется нажатием ЛКМ, при этом она становится белой. Кнопка Split позволяет разбить выделенную полосу на две независимых части. Нажав ПКМ на красной полосе, можно редактировать ограничение в числовом виде.

Сдвиньте верхнее ограничение так, чтобы ограничить перерегулирование примерно значением 10%. Правую границу первого блока сдвиньте примерно на время 60-65 сек.

Выберите пункт верхнего меню Optimization – Parameters и введите Kaw в поле Tunable variables (настраиваемые параметры). В поле Lower bound (нижняя граница) введите 0 (нуль).

Запустите процедуру оптимизации. Она может занимать несколько минут. Признаком окончания будет появление текста Optimization Converged Successfully в командном окне Matlab.

ЛКМ по кнопке Start

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

Определите найденное оптимальное значение Kaw. Для этого надо ввести Kaw в командном окне Matlab. Запишите его в отчет и перенесите в скрипт lab6go.

Kaw

Запустите скрипт lab6go на выполнение и скопируйте полученный график в отчет.

клавиша F5

Определите перерегулирование и время переходного процесса, используя функцию overshoot.

[si,Tpp] = overshoot(phi(:,1),...

phi(:,4));

Измените в командном окне Matlab заданный угол курса на и повторите моделирование. Скопируйте полученный график в отчет. Определите перерегулирование и время переходного процесса.

в командном окне

phiZad = 30;

выделить нужные строки скрипта, нажать F9

Сделайте выводы о применении схемы с компенсацией насыщения.

Таблица коэффициентов

Вариант

, сек

, рад/сек

, сек

, сек

16.0

0.06

1

1

16.2

0.07

2

2

16.4

0.08

1

3

16.6

0.07

2

4

16.8

0.06

1

5

17.0

0.07

2

6

17.2

0.08

1

1

17.4

0.07

2

2

17.6

0.06

1

3

17.8

0.07

2

4

18.0

0.08

1

5

18.2

0.09

2

6

18.4

0.10

1

1

18.6

0.09

2

2

18.8

0.08

1

3

19.0

0.07

2

4

19.2

0.08

1

5

19.4

0.09

2

6

19.6

0.10

1

1

18.2

0.0694

2

6