СИНТЕЗ ЦИФРОВОЇ СИСТЕМИ КЕРУВАННЯ ІЗ СТАНДАРТНИМ ПІ – АЛГОРИТМОМ КЕРУВАННЯ
4.1 Задача розрахунку
Управління технологічним об’єктом із використанням керуючої ЕОМ здійснюємо методом безпосереднього цифрового керування (БЦК).
Розрахунок безпосереднього цифрового керування (БЦК) зводиться до вибору таких параметрів налаштування системи, якій забезпечували наперед задану якість перехідних процесів. Регулятор в системі БЦК реалізується у вигляді програм на базі мікропроцесорної техніки.
При розробці системи БЦК необхідно розв’язати наступні задачі:
розробка алгоритму БЦК;
вибір оптимального кроку квантування (періодичність опитування давачів і видачі керуючих виливів);
вибір керуючої ЕОМ;
реалізація програмного забезпечення;
вибір вимог до пристроїв зв’язку з об’єктом при узгодженні ЕОМ із виконавчими механізмами.
В системах БЦК існують аналоги алгоритмів аналогових неперервних регуляторів і тому задача розрахунку регулятора БЦК заключається у визначенні коефіцієнта передачі, часу ізодрому, часу визначення регулятора.
Є і додатковий параметр налаштування – h – період дискретності, який впливає на якість перехідного процесу та інші показники системи з БЦК.
4.2 Вибір оптимального кроку квантування
В системах БЦК вивід сигналів із здавачів і видача керуючих впливів на об’єкт здійснюється в дискретні моменти часу.
Задача оптимального вибору періоду квантування відноситься до основної проблеми оптимального витрачання машинного часу, який витрачається на кожний контур БЦК.
Вибір оптимального періоду квантування зводиться до налагодження точності керування від величини періоду квантування керуючих впливів u [nh] і керуючої змінної x [nh]. При великому h погіршується якість керування і знижується стійкість системи. При малому h збільшується завантаження ЕОМ, не економічно витрачається машинний час. Тому приходять компромісного рішення.
Існують різні методи вибору періоду квантування. Період квантування обчислюємо трьома методами:
на основі метода Гурмана;
на основі В-характеристик;
на основі практичного досвіду.
Метод Гурмана
Період квантування визначається як
(4.1)
Метод В-характеристик
(4.2)
де (4.3)
де - абсолютна похибка вихідної величини, =
Отже ,
Практичний досвід
(4.4)
де -час, за який вихідна величина набула значення 95% від усталеного значення.
Вихідна величина змінюється за законом
тоді
Отже, нехай
Визначаємо які періоди квантування відносяться до великих, а які до малих. Малі кроки квантування, якщо
Великі:
Отже, крок і - великі і - малий.
4.3 Розрахунок параметрів налаштування цифрового регулятора
На базі ЕОМ можна реалізувати самі складні і ефективні закони регулювання. Із стандартних законів регулювання найбільш поширені – ПІ і – ПІД закони регулювання.
В пункті 4.2 було доведено доцільність використання ПІ-регулятора.
Методика визначення параметрів налаштування цифрового регулятора для великих і малих періодів квантування різна.
Для малих періодів квантування параметри цифрового регулятора розраховуються за значеннями Кр,Ті, неперервного ПІ-регулятора.
Для великих тактів квантування способи дискретизації неперервних диференційних рівнянь стають несправедливими.
Поведінка систем оцінюється на основі квадратичного критерію якості
(4.5)
де - ваговий коефіцієнт;
- відхилення керуючої змінною від її усталеного значення .
Параметри цифрового регулятора визначають в результаті мінімізації критерію ….
При малих кроках квантування справедлива нерівність а при великих кроках квантування –
4.3.1 Розрахунок параметрів налаштування ПІ-регулятора при великому кроці квантування.
Для класу об’єктів керування типу , існують декілька методів визначення параметрів налаштування дискретних регуляторів. Розглянемо дві методики для ПІ-регулятора:
, (4.6)
(4.7)
Тобто для об’єкта де і
Передавальна функція ПІ-регулятора:
. (4.8)
За методом Такахаші:
, (4.9)
(4.10)
Тобто для =20с і
,
,
Передавальна функція ПІ-регулятора:
,
Отже .
Визначаємо параметри регулятора для кроку квантування за методом Такахаші:
,
.
4.4 Дослідження стійкості системи
Досліджуємо одноконтурну цифрову систему керування на стійкість
Структурна схема системи приведена на рисунку (4.1)
y U(nT) U(t) y(t)
–
n
Рисунок 4.1 – Цифрова система керування
де Д(z) - передавальна функція цифрового регулятора;
- передавальна функція екстраполятора;
- передавальна функція об’єкта.
Визначаємо передавальну функцію неперервної частини
,
У z-формі:
,
де
Визначаємо передавальну функцію приведеної частини для і
Отримані результати означають, що початок перехідного процесу зміститься відповідно на 20с і 10с.
Знаходимо передавальні функції замкнутої системи для різних кроків квантування:
;
Знаходимо полюси функції:
Оскільки полюси знаходяться в середині кола радіуса R=1, то система стійка.
Д1(z) Wп1(z)=,
.
Полюси функції:
Z2 -1,5411z + 0,5675=0,
Z1,2 = 0,770550,16184,
Z 1= 0,9323; Z 2 =0,6081
Оскільки:
то система стійка.
Д3(z) Wп2(z)= ,
Ф3,
Полюси функції:
Z2 – 1б728z +0,7786=0,
Z1,2 =0,8640,179j,
<0
- система стійка.
,
,
Полюси функції:
,
Оскільки:
то система стійка.
4.5 Визначення якості перехідних процесів
Параметри якості перехідного процесу визначаємо із отриманих при моделюванні на ЕОМ кривих перехідних процесів (рисунок 1,2,3,4). Моделювання здійснюємо за допомогою…..
Результати зведемо в таблицю 4.1
Таблиця 4.1 - Параметри якості перехідних процесів
Метод визначення | Крок квантування,n | Передавальна функція регулювання, Д(z) | Показники ПП
перерегулювання , % | час регулю- вання, t, с
Продовження таблиці 4.1
Практичний |
20
10 |
~1
1,05 |
~ 200
250
Такахаші |
20
10 |
60 %
45 % |
~ 320
250
4.5 Дослідження системи із малим кроком квантування
Параметри цифрового регулятора для малого кроку квантування знаходяться за параметрами неперервних регуляторів. Отже визначаємо параметри неперервного Пі-регулятора. Використовуємо метод Циглера –Нікольса. Знаходимо критичну частоту.
Розв’язуючи рівняння:
Знаходимо пропорційне критичне налаштування, розв’язуючи рівняння:
Параметри налаштування регулятора:
;
тоді коефіцієнти пропорційності
,
стала ізодрому - .
Параметри цифрового регулятора для кроку квантування n=5c
Передавальна функція цифрового регулятора:
,
Параметри цифрового регулятора для кроку квантування n=2c
Для
Визначаємо передавальні функції розімкнутих і замкнутих систем для кроку n=2c
Для кроку n=5c
,
.
Досліджуємо ці системи на стійкість. Знаходимо полюси замкнутих функцій:
Z2 – 1,918z +0,9201=0,
Z1,2 =0,959,
,
Z2 – 1,7935z +0,80775=0,
Z1,2 =0,89675,
.
Отже системи стійкі.
Перехідні процесу змоделюємо на ЕОМ. Результати приведені на рисунку
3…, …., ….,