
- •2.4 Определение устойчивости линейной сау по критерию
- •1.2 Функциональная схема сау
- •1.3 Структурная схема сау
- •2 Исследование линейной сау
- •2.1 Расчет передаточной функции
- •2.2 Определение устойчивости линейной сау по критерию Гурвица
- •2.5 Переходный процесс линейной системы и определение показателей качества
- •2.6 Построение ачх линейной системы и определение показателей качества
- •2.7 Построение лачх и определение запасов устойчивости
- •2.8 Аппроксимация лачх и определение передаточной функции
- •2.9 Вывод по исследованию линейной системы
- •3 Исследование нелинейной системы
- •3.1 Преобразование системы
- •3.2 Построение фазового портрета нелинейной системы
- •3.3 Вывод по исследованию нелинейной системы
- •4 Исследование дискретной системы
- •4.2 Определение устойчивости дискретной системы по критерию Шур-Кона
- •4.3 Переходный процесс дискретной системы
- •4.4 Билинейное преобразование дискретной системы
- •4.5 Построение лачх дискретной системы и определение запасов устойчивости
- •4.6 Вывод по исследованию дискретной системы
2.5 Переходный процесс линейной системы и определение показателей качества
Используя математический редактор Matlab, построим переходный процесс системы по передаточной функции (2).
A(w)
t, с
Рисунок 6 - Переходный процесс системы
2.6 Построение ачх линейной системы и определение показателей качества
Используя математический редактор MathCAD, построим амплитудно-частотную характеристику системы.
Заменим в уравнении (3)
:
(30)
(31)
ω2
Amax
A0
A(w)
w, Гц
ωp
Amax
ω1
Рисунок 7 - Амплитудно–частотная характеристика линейной системы
Определение показателей качества:
А0 = 0;
Аmax = 285;
резонансная частота
=
0,2, Гц;
показатель колебательности
;
полоса пропускания
=
201,525
2.7 Построение лачх и определение запасов устойчивости
Используя математический редактор MATLAB, построим ЛАЧХ системы.
>> w=tf([48746880],[1,6 4,9 18 35 40 34 22 7,5 1])
Transfer function:
48746880
-----------------------------------------------------------------------------------------
1,6 s^8 + 4,9 s^7 + 18 s^6 + 35 s^5 + 40 s^4 + 34 s^3 + 22 s^2 + 7,5 s +1
>> bode(w)
∆L

∆φ
L(w), Дб
φ(w), deg
w, рад/с
Рисунок 8 - ЛАЧХ Линейной системы
∆L=150о
∆φ=-520о
2.8 Аппроксимация лачх и определение передаточной функции
w, рад/с
Рисунок 9 - Аппроксимация ЛАЧХ линейной системы системы
Передаточная функция аппроксимированной ЛАЧХ имеет вид
(32)
С учетом коэффициентов выражение (32) примет вид:
(33)
2.9 Вывод по исследованию линейной системы
В ходе исследования линейной системы была получена передаточная функция системы. Определили, что система является устойчивой по критериям устойчивости Гурвица, Михайлова и Евсюкова. Построили переходный процесс и АЧХ системы, определили показатели качества. При помощи математического редактора MATLAB построили ЛАЧХ системы и определили запасы устойчивости.
3 Исследование нелинейной системы
3.1 Преобразование системы
U1(t)
U1(t)
U1(t)
U1(t)
U1(t)
U(t)
Y(t)

Рисунок 10 - Функциональная схема САУ с нелинейным элементом
Функциональная схема состоит из следующих элементов:
ЗТ – заквасочный танк для производства сыра;
ВС – ванна сыродельная;
ФА – формовочный аппарат для производства сыра;
П – пресс для сыра;
УП – Установка для посола.
Н.Э. – нелинейный элемент.
График, описывающий нелинейный элемент приведен на рисунке 11.
Рисунок 11 - Cтатическая характеристика нелинейного элемента
Применяя правила преобразования структурных схем, упростим схему, изображенную на рисунке 10.
Wобщ (p)
Н.Э.
Рисунок 12 – Структурная схема нелинейной системы
Введем вынужденную обратную связь.
Рисунок 12 - Итоговое преобразование САУ с нелинейным элементом
.
(34)
С учетом всех преобразований Wобщ (p) примет вид:
(35)
3.2 Построение фазового портрета нелинейной системы
Передаточная функция есть
или
,
(27)
где
-передаточная
функция линейной системы;
Подставляя в формулу (7) значение передаточной функции получим:
Приведенную формулу можно записать в виде:
(29)
Введем замену
и
и исключим из правой части уравнения
производную:
Для того чтобы построить фазовый портрет, необходимо, чтобы степень числителя и знаменателя не превышала вторую степень, поэтому элементы выше второй степени исключаем. Тогда получим:
(31)
(32)
Так как в качестве нелинейного элемента
используется реле со статической
характеристикой, то подставляя значение
для трех участков, получим систему
уравнений:
(33)
Создадим матрицу для решения дифференциального уравнения.
Зададим матрицы для трех начальных условий:
(34)
Возьмем количество точек равным 1000 и конечное время интегрирования 1, то матрица решений запишется как:
(37)
Построим фазовый портрет нелинейной системы :
ν
Рисунок 13 - Фазовый портрет нелинейной системы
h(t)
t
Рисунок 14 - Переходный процесс нелинейной системы