Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лаб_ТАУ_ЛС 9.doc
Скачиваний:
0
Добавлен:
01.05.2025
Размер:
339.46 Кб
Скачать

Лабораторна робота № 9. Проектування регулятора для лінійної системи

Мета роботи. Освоєння методів проектування регулятора для одновимірної лінійної безперервної системи за допомогою середовища Matlab

Програма роботи

  • навчитися будувати моделі з'єднань лінійних ланок

  • навчитися використовувати модуль SISOTool для проектування простих регуляторів

Оформлення звіту

Звіт по лабораторній роботі виконується у вигляді зв'язного (читаного) тексту у файлі формату Microsoft Word (шрифт основного тексту Times New Roman, 12 пунктів, через 1,5 інтервалу, вирівнювання по ширині). Він повинен включати

  • назву предмету, номер і назва лабораторної роботи

  • прізвище і ініціали авторів, номер групи

  • прізвище і ініціали викладача

  • номер варіанту

  • короткий опис досліджуваної системи

  • результати виконання всіх пунктів інструкції, які виділені сірим фоном (див. нижчий): результати обчислень, графіки, відповіді на питання.

При складанні звіту рекомендується копіювати необхідну інформацію через буфер обміну з робочого вікна середовища Matlab. Для цих даних використовуйте шрифт Courier New, в якому ширина всіх символів однакова.

Опис системи

У роботі розглядається система стабілізації судна на курсі. Її структурна схема показана на малюнку.

Рис. 9.1. Структурна схема системи автоматичного управління

Лінійна математична модель, що описує рискання судна, має вигляд

де – кут відхилення, – кутова швидкість обертання навколо вертикальної, осі – кут повороту щодо положення рівноваги, – постійна часу, – постійний коефіцієнт, що має розмірність рад/сек. Передавальна функція відношення кута повороту до кута відхилення запишеться у вигляді

.

Привід (рульова машина) приблизно моделюється як інтегруюча ланка

,

охоплена одиничним негативним зворотним зв'язком.

Для вимірювання кута рискання використовується гірокомпас, математична модель якого записується у вигляді аперіодичної ланки першого порядку з передавальною функцією

,

Інструкція по виконанню роботи.

Основна частина команд вводиться в командному вікні середовища Matlab. Команди, які треба застосовувати в інших вікнах, позначені іконами відповідних програм.

Етап виконання завдання

Команди Matlab

  1. Введіть передавальну функцію моделі судна як об'єкт tf.

P = tf ( K, [Ts 1 0] )

  1. Введіть передавальну функцію інтегруючої ланки .

R0 = tf ( 1 [TR 0] )

  1. Побудуйте передавальну функцію рульового пристрою, замкнувши інтегратор одиничним негативним зворотним зв'язком.

R = feedback ( R0, 1 )

  1. Побудуйте передавальну функцію послідовного з'єднання об'єкту з приводом.

G = P * R

  1. Побудуйте перехідну характеристику для отриманої моделі і скопіюйте її в звіт через буфер обміну. Поясніть, чому функція нескінченно зростає і прямує до прямої. Який коефіцієнт нахилу цієї прямої? Закрийте вікно з графіком.

step ( G )

  1. Побудуйте передавальну функцію вимірювального пристрою .

H = tf ( 1 [Toc 1] )

  1. Побудуйте передавальну функцію розімкненого контуру.

L = G * H

  1. Побудуйте ЛАФЧХ розімкненої системи1.

bode ( L )

  1. Відзначте крапки, що визначають перетин ЛАЧХ з прямою 0 дБ і перетин ЛФЧХ з прямою -1800.

ПКМ – Characteristics – Stability (Minimum Crossing)

  1. Визначите, чи є замкнута система стійкою? Які запаси стійкості за амплітудою (Gain margin) і фазою (Phase margin)? Який регулятор неявно використовується в цьому випадку? Скопіюйте графік ЛАФЧХ в звіт.

ЛКМ на мітках-кружках

  1. Знайдіть максимальний коефіцієнт посилення розімкненої системи. Поясните цей результат.

ПКМ – Characteristics –

Peak Response

  1. Закрийте вікно з ЛАФЧХ і запустите модуль SISOTool.

sisotool

  1. Імпортуйте передавальну функцію G як модель об'єкту (Plant) і H як модель датчика (Sensor). Блоки F (передфільтр) і C (регулятор) залиште без змін (рівними 1).

File - Import

  1. Відключіть зображення кореневого годографа так, щоб у вікні залишилася тільки ЛАФЧХ.

View – Root Locus (відключити)

  1. Для того, щоб відразу бачити зміни перехідних процесів, запустіть LTIViewer 2з верхнього меню вікна SISOTool. Розташуєте два вікна поряд, щоб вони не перекривали одне одного.

Analysis –

Response to Step Command

  1. Залиште тільки графік перехідного процесу на виході, відключивши виведення сигналу управління.

ПКМ – Systems –

Closed loop r to u

  1. Визначите перерегулювання і час перехідного процесу 3. Скопіюйте графік в звіт.

ПКМ – Characteristics –

  • Peak Response

  • Settling Time

  1. Перейдіть у вікно SISOTool. Визначите коефіцієнт посилення, при якому перерегулювання приблизно дорівнює 10%. Як змінився час перехідного процесу? Які запаси стійкості в цьому випадку? Скопіюйте графік в звіт.

перетягання мишею ЛАЧХ

редагування в полі Current Compensator

  1. Перейдіть у вікно середовища Matlab і введіть передавальну функцію пропорційно-диференціального (ПД) регулятора

де сек

а – постійна часу судна.

Cpd = 1 + tf ( [Ts 0], [Tv 1] )

  1. Перейдіть у вікно SISOTool. Імпортуйте регулятор Cpd як базову модель для блоку C.

File – Import, Cpd -> C

  1. Визначите додатковий коефіцієнт посилення, при якому перерегулювання приблизно дорівнює 10%. Знайдіть час перехідного процесу і запаси стійкості. Порівняєте пропорційний і ПД-регулятори. Скопіюйте в звіт графік перехідного процесу.

перетягання мишею ЛАЧХ

редагування в полі Current Compensator

  1. Визначите додатковий коефіцієнт посилення, при якому час перехідного процесу мінімальний. Скопіюйте в звіт графік перехідного процесу.

перетягання мишею ЛАЧХ

редагування в полі Current Compensator

  1. Експортуйте отриманий регулятор в робочу область Matlab.

File – Export

у стовпці Export as змінити ім'я Cpd на C

кнопка Export to workspace

  1. Побудуйте передавальну функцію отриманої замкнутої системи. Подумайте, чому вийшов такий громіздкий вираз. Який повинен бути порядок передавальної функції?

W = C*G / (1 + C*G*H)

  1. Побудуйте мінімальну реалізацію передавальної функції W.

W = minreal(W)

  1. Визначите полюси передавальної функції замкнутої системи. Що означає близькість деяких полюсів до уявній осі? Чи вірно, що в цьому випадку буде малий запас стійкості?

pole ( W )

  1. Знайдіть коефіцієнт підсилення системи в усталеному режимі. Поясніть результат. Чи є у такої системи статична помилка при відстежуванні постійного сигналу? Чому? А для лінійно зростаючого сигналу?

dcgain ( W )

  1. Як зміниться статичний коефіцієнт посилення, якщо модель датчика прийме вигляд

?

  1. Побудуйте мінімальну реалізацію передавальної функції замкнутої системи від входу до сигналу управління (виходу регулятора).

Wu = minreal(C/ (1 + C*G*H))

  1. Побудуйте зміну сигналу управління при одиничному ступінчастому вхідному сигналі і скопіюйте графік в звіт. Поясните, чому сигнал управління прямує до нуля.

step ( Wu )

Таблиця коефіцієнтів

Варіант

сек

рад/сек

сек

сек

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]