Контрольная работа по ТАУ в MatLab
ДОМАШНЕЕ ЗАДАНИЕ К ТЕМЕ 10
ВАРИАНТ 1
Выбрать регулятор (закон регулирования) для статического объекта, передаточная функция которого приведена ниже.
Дано:
;
· Максимально-допустимое значение динамического отклонения регулируемой величины εmax = 9 %;
· Принятый типовой переходный процесс – апериодический переходный процесс;
· Максимально-возможное значение возмущения на объект, выраженное в % хода регулирующего органа Δхвх = 4%.
· Допустимое остаточное отклонение, т.е. статическая ошибка, εст = 5%
· Допустимое время регулирования Тр = 20 сек.
Для упрощения расчетов и исследования характеристик САУ используем рекомендованный пакет прикладных программ Control System Toolbox, который является приложением системы MatLab и ориентирован на решение задач по теории автоматического управления.
Для использования функций Control System Toolbox в первую очередь следует ввести в среду MatLab все полученные ПФ. Наиболее удобно это делать в так называемой tf –форме.
Для создания tf-модели данной ПФ используется функция с одноименным идентификатором
H=tf ( [bm, bm-1, …, b0],[an, an-1,…,a0] ),
где H – имя созданной ПФ, […] – массивы коэффициентов полиномов числителя и знаменателя ПФ.
%Создадим tf-модель объекта управления:
Wou=tf([2.5], [8 1],'inputdelay',4);
% Определим переходный процесс при максимально-возможном значение возмущения на объект, выраженное в % хода регулирующего органа Δхвх = 4%.
step (0.04*Wou) ; grid on
Рис. 1. Переходный процесс исходного объекта.
% Определим качество переходного процесса
stepinfo(0.04*Wou)
ans =
RiseTime: 17.5778
SettlingTime: 35.2975
SettlingMin: 0.0901
SettlingMax: 0.1000
Overshoot: 0
Undershoot: 0
Peak: 0.1000
PeakTime: 80.6826
Переходный процесс удовлетворяет заданным критериям:
· по максимально-допустимому значение динамического отклонения регулируемой величины εmax = 9 %;
· по типу переходного процесса – апериодический переходный процесс;
· по допустимому остаточному отклонению, т.е. статическая ошибка, εст = 0%.
Но не удовлетворяет по допустимому времени регулирования Тр = 20 сек., которое равно 35,2975 с.
Введем корректирующее форсирующее звено и подберем его параметры.
Wkor= tf([7 1], [0 1]);
W=Wou*Wkor;
%W=feedback(Wou, Wpi)
step (0.04*W) ; grid on
Рис. 2. Переходный процесс скорректированного объекта.
% Определим качество переходного процесса
stepinfo(0.04*W)
ans =
RiseTime: 2.0607
SettlingTime: 18.6626
SettlingMin: 0.0901
SettlingMax: 0.1000
Overshoot: 0
Undershoot: 0
Peak: 0.1000
Таким образом, время переходного процесса уменьшилось до 18.6626 с, что соответствует заданным требованиям.