Рубцов В. И.
ВЫПОЛНЕНИЕ ДОМАШНЕГО ЗАДАНИЯ №2 ПО КУРСУ ТЕОРИИ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ ЛИНЕЙНЫЕ СИСТЕМЫ С ПРИМЕНЕНИЕМ ПРОГРАММЫ MATLAB
ВВЕДЕНИЕ
При выполнении сложных расчетных заданий в курсе теории автоматического управления не обойтись без применения специализированных компьютерных программ. Одной из таких программ является Matlab и его приложения Control System Toolbox, Simulink и Symbolic Toolbox. В данном пособии, на примере выполнения типового домашнего задания, приведены основные методы работы с указанными приложениями. Также приведены основные теоретические положения для определения передаточные функций, для расчета переходной, импульсно переходной функций, определения установившееся и начальное значения переходной функции.
СОДЕРЖАНИЕ ЗАДАНИЯ
Для системы автоматического управления, структурная схема которой представлена следующее:
1. Найти передаточные функции разомкнутой системы, замкнутой системы и передаточную функцию замкнутой системы по ошибке.
2. Написать дифференциальное уравнение замкнутой системы для переменных x(t) и ε(t).
3. Построить ЛАЧХ и ФЧХ разомкнутой системы. Определить запасы устойчивости системы по модулю и по фазе.
а) ассимптотические амплитудные и приближенные фазочастотные характеритики с помощью средств ручного счета.
б) точные характеристики рассчитать с помощью ЭВМ.
4. Определить установившееся значение ошибки при входном воздействии
а) g(t)=g0=const
б) g(t)=g0+g1×t
5. Определить устойчивость замкнутой системы используя критерии Гурвица и Найквиста (в последнем случае по ЛАЧХ и ЛФЧХ).
6. С помощью ЭВМ найти переходную функцию системы. Построить график этой функции
7. С помощью ЭВМ найти импульсную переходную функцию системы. Построить график этой функции.
ПРИМЕНЕНИЯ ПРОГРАММЫ MATLAB ДЛЯ ВЫПОЛНЕНИЯ ДОМАШНЕГО ЗАДАНИЯ
При анализе и проектировании современных систем автоматического управления не обойтись без использования специализированных компьютерных программ. Одной из таких программ является Matlab, содержащей ряд специализированных приложений, облегчающие анализ и синтез систем автоматического регулирования. Control System Toolbox является одним из таких приложений. Это приложение предназначено для описания непрерывных и дискретных систем управления с постоянными параметрами. При реализации этого пакета использованы принципы объектно-ориентированного программирования. Введен новый класс объектов: линейные системы с постоянными параметрами. Этот класс включает линейные непрерывные и и дискретные, одномерные с одним входом и одним выходом и многомерные со многими входами и (или) многими выходами. Дополнительные параметры позволяют указать, является ли данная система непрерывной или дискретной, имеет ли запаздывание на входе системы. В пакете прикладных программ Control System Toolbox введен ряд функций, позволяющих определять динамические параметры системы: нули и полюса, динамические параметры системы, собственные частоты, коэффициенты демпфирования, и.т.д. В пакете введен ряд функций позволяющих осуществлять синтез регулятора на основе задания желаемого расположения полюсов или минимального значения квадратического критерия качества.
Другим, не менее полезным приложением Matlab, является Symbolic Toolbox. Данное приложение позволяет решать обыкновенные и дифференциальные уравнения, системы уравнений. Применять прямое и обратное преобразования Фурье и Лапласа, что позволяет получать переходные и импульсно переходные функции систем автоматического управления.
ВЫПОЛНЕНИЕ ЗАДАНИЯ
При описании выполнения домашнего задания №2 для вариантов 1-16 и 17-32, будет различаться только выполнение пункта 1. При выполнении пунктов со 2 по 7 мы будем пользоваться полностью идентичными командами.
ВЫПОЛНЕНИЕ ЗАДАНИЯ: ВАРИАНТЫ ЗАДАНИЯ 1-16
Пункт 1
Структурная схема изображена на Рис. 1.
Для выполнения первого пункта задания воспользуемся приложением Control System Toolbox. Зададим передаточные функции отдельных звеньев. Для этого воспользуемся командой tf (transfer function)
W1=tf([τ 1],[T 1]);
W2=tf([k1],[T1 1 0]).
W3=tf([k2],[1 0])
Для упрощения схемы воспользуемся правилами структурного преобразования.
Для этого применим правило структурного преобразования для систем с местной обратной связью и команду feedback.
W4=feedback(W2,1).
Полученная упрощенная схема приведена на рис. 2. Далее воспользуемся правилом упрощения для последовательно соединенных звеньев:
Рис. 1
Рис. 2
Рис. 3
W=series(W1,W4,W3)И получим следующую эквивалентную схему, изображенную на рис. 3. Как видно из рис. 3, мы получили передаточную функцию разомкнутой системы. Для вычисления передаточной функции замкнутой системы по управляющему воздействию используем команду feedback.
Fx=feedback(W,1)
Для вычисления передаточной функции замкнутой системы по ошибке воспользуемся соотношением: Fε=1-Fx
Fe=1-Fx.
ВЫПОЛНЕНИЕ ЗАДАНИЯ: ВАРИАНТЫ ЗАДАНИЯ 17-32
Пункт 1
Структурная схема изображена на Рис. 4.
Для выполнения первого пункта задания воспользуемся приложением Control System Toolbox. Зададим передаточные функции отдельных звеньев. Для этого воспользуемся командой tf (transfer function)
W1=tf([τ 1],[T 1]);
W2=tf([k1],[T1 1 0]).
W3=tf([k2],[1 0])
W4=tf(1,[1 0])
Для упрощения схемы воспользуемся правилами структурного преобразования.
Для этого применим правило структурного преобразования для систем с местной обратной связью и команду feedback.
W24=feedback(W2,W4).
Полученная упрощенная схема приведена на рис. 5. Далее воспользуемся правилом упрощения для последовательно соединенных звеньев:
W=series(W1,W24,W3)
И получим следующую эквивалентную схему, изображенную на рис. 3. Как видно из рис. 6, мы получили передаточную функцию разомкнутой системы. Для вычисления передаточной функции замкнутой системы по управляющему воздействию используем команду feedback.
Fx=feedback(W,1)
Для вычисления передаточной функции замкнутой системы по ошибке воспользуемся соотношением: Fε=1-Fx
Fe=1-Fx.
Пункт 2
Пункт 2 сложности не представляет и для его выполнения компьютер не требуется
Пункт 3 б)
Для построения ЛАЧХ и ЛФЧХ системы с разомкнутой главной обратной связи воспользуемся следующими командами:
bode(W) – данная команда позволяет построить ЛАЧХ и ЛФЧХ произвольной передаточной функции
margin(W)- данная команда позволяет построить ЛАЧХ и ЛФЧХ произвольной передаточной функции и оценить запасы устойчивости по амплитуде и по фазе.
Пункт 4 а)
Для выполнения данного пункта воспользуемся приложением Symbolyc Toolbox.
Сначала определим символьные переменные.
Syms s
Зададим передаточную функцию системы по ошибке в символьном виде:
Fe=(a0*s^4+a1*s^3+a2*s^2+a3*s+a4)/(b0*s^4+b1*s^3+b2*s^2+b3*s+b4)
Рис. 4
Рис. 5