- •В.И. Гаркушенко теория автоматического управления
- •Введение
- •1. Задание на курсовую работу
- •1.1 Cиловые системы автоматического управления
- •1.1.1 Исходные данные
- •1.1.2 Техническое задание
- •1.2 Содержание пояснительной записки
- •2. Порядок выполнения работы
- •2.1 Расчет сау в линейном приближении
- •2.2 Расчет сау с учетом нелинейности
- •3. Практические рекомендации
- •3.1 Введение
- •3.2 Математическое описание функциональных элементов
- •3.2.1 Генератор постоянного тока
- •3.2.2 Тахогенератор
- •3.2.3 Электромашинный усилитель
- •3.2.4 Двигатель постоянного тока
- •3.3 Определение передаточных функций замкнутых систем
- •3.4 Определения коэффициента усиления электронного усилителя
- •3.5 Построение области устойчивости по коэффициенту усиления
- •3.5.1 Алгебраический метод построения области устойчивости
- •3.5.2 Частотный метод построения области устойчивости
- •3.6 Построение лах и лфх разомкнутой системы и определение запасов устойчивости замкнутой системы
- •3.7 Построение переходной характеристики замкнутой системы
- •3.8 Синтез последовательной коррекции с помощью лах и лфх
- •3.9 Электрическая схема корректирующего устройства
- •3.10 Синтез дискретного корректирующего устройства
- •3.10.1 Приближенный синтез дискретного корректирующего устройства
- •3.10.2 Синтез дискретного корректирующего устройства частотным методом
- •3.11 Определение параметров автоколебаний
- •3.12 Исследование системы на абсолютную устойчивость
- •Список литературы
- •Приложение а. Корректирующие цепи
- •Приложение б. Коэффициенты гармонической линеаризации нелинейных характеристик
- •Приложение в. Справочные данные
3.6 Построение лах и лфх разомкнутой системы и определение запасов устойчивости замкнутой системы
Построение ЛАХ и ЛФХ разомкнутой системы можно проводить ручным и машинным способом.
Для использования ручного способа построения асимптотической ЛАХ и ЛФХ необходимо, чтобы исходная передаточная функция была представлена в стандартном виде, т.е. в виде произведения типовых звеньев. При необходимости приведение к стандартному виду произвольной передаточной функции можно проводить с помощью функции zpk.Например, для передаточной функции
при выполнении командной строки
>> W=tf([3 6 2],[200 10 3 0]),W1=zpk(W)
Transfer function:
3 s^2 + 6 s + 2
----------------------
200 s^3 + 10 s^2 + 3 s
Zero/pole/gain:
0.015 (s+1.577) (s+0.4226)
--------------------------
s (s^2 + 0.05s + 0.015)
исходная передаточная функция представлена в виде Zero/pole/gain, т.е. в виде произведения типовых звеньев, у которых свободные члены отличны от единицы. Для приведения передаточной функции к стандартному виду необходимо свободные члены вынести за скобки.
Ручной способ построения асимптотической ЛАХ и ЛФХ по передаточной функции разомкнутой системы проводится в следующей последовательности:
исходная передаточная функция представляется в стандартном виде:
,
где – коэффициент усиления разомкнутой системы;,; порядок полинома числителяне превышает порядка знаменателя, что соответствует физически реализуемым системам;
определяется значение ;
определяются сопрягающие частоты , вычисляются;
строится асимптотическая ЛАХ:
через точку проводится слева направо прямая с наклономдБ/дек до первой слева сопрягающей частоты, где– число интегрирующих звеньев в передаточной функции (если, то прямая проводится параллельно оси частот; если вместо интегрирующих звеньев присутствуютдифференцирующих звеньев, то следует принятьи наклон асимптоты будет положительным);
в сопрягающей частоте ЛАХ терпит излом относительно предыдущего участка на дБ/дек илидБ/дек. Если сопрягающей частоте соответствует звено первого порядка, то излом составляетдБ/дек, для звена второго порядка (с комплексно-сопряженными корнями) соответственнодБ/дек. Знак "+" соответствует звеньям, расположенных в числителе, а знак "–" в знаменателе передаточной функции;
далее проводится прямая до следующей сопрягающей частоты, в которой ЛАХ терпит излом аналогично предыдущему;
в области высоких частот ЛАХ уходит в бесконечность с наклоном дБ/дек, где– порядок числителя,– порядок знаменателя передаточной функции.
строится ЛФХ в виде суммы ЛФХ типовых звеньев, входящих в передаточную функцию, при этом в области низких частот ЛФХ начинается со значения рад, а в области высоких частот ЛФХ стремится к значениюрад.
При ручном способе построения ЛАХ и ЛФХ следует использовать масштаб: по оси ординат 20 дБ - 40 мм и I град - I мм; по оси частот 1 дек - 100 мм. При построении ЛФХ воспользоваться шаблонами ЛФХ апериодического и колебательного звена с заданным коэффициентом затухания, построенных с помощью таблиц [1] и вырезанных из ватмана.
При машинном способе построения с использованием системы MATLAB строится точная ЛАХ, на которую следует нанести график асимптотической ЛАХ, которая в дальнейшем используется для синтеза последовательного корректирующего устройства.
Рассмотрим машинный способ построение ЛАХ и ЛФХ с помощью системы MATLAB для рассмотренной выше передаточной функции разомкнутой системы схемы А:
,
для принятых параметров ;;;;;.
Воспользуемся вспомогательной программой, составленной в Script-файле:
% передаточная функция разомкнутой системы
Wpas=75.5*tf([0.4 1],[0.04 1])*tf([1],[0.2 1])*tf([1],[0.96 1])*tf([1],[0.38 1]);
% сопрягающие частоты
omega=[0.1 1/0.96 1/0.4 1/0.38 1/0.2 1/0.04 1/0.001];
% построение ЛАХ и ЛФХ, запасов устойчивости по амплитуде и фазе
margin(Wpas); grid on; figure
% построение асимптотическая ЛАХ в заданной области частот по передаточной функции
L1=20*log10(75.5);L2=20*log10(75.5);L3=L2-20*log10(0.96/0.4);
L4=L3;L5=L4-20*log10(0.38/0.2);L6=L5-40*log10(0.2/0.04);L7=L6-60*log10(0.04/0.001);
L=[L1 L2 L3 L4 L5 L6 L7]; semilogx(omega,L); grid on; hold on
% точная ЛАХ для сравнения
[Lg,f,w]=bode(Wpas,{0.1,1000});Lg1=20*log10(squeeze(Lg)); semilogx(w,Lg1,'--')
Приведенный текст программы можно скопировать и вставить в рабочее поле m-file, которое открывается при нажатии левой верхней кнопки в среде MATLAB. Затем выполнить программу нажатием кнопки со стрелкой (Run) на панели инструментов m-file. При этом вычисленные переменные в Script-файле являются глобальными и доступны в любом другом Script-файле, а также в среде MATLAB.
Если программу оформить в виде m-файла с помощью оператора function, то переменные будут локальными в пределах данного m-файла.
В результате выполнения получим графики представленные на рис. 17, рис.18.
Рис. 17
Рис. 18
Согласно критерию Найквиста для устойчивости замкнутой отрицательной единичной обратной связью физически реализуемой разомкнутой системы с передаточной функцией ,необходимо и достаточно, чтобы АФЧХ охватывала точкус учетом знака в суммераз, гдечисло правых корней характеристического уравнения разомкнутой системы. При этом следует считать, что при наличиинулевых корней уравненияхарактеристикапридополняется дугой бесконечно большого радиуса с раствором угла, т.е. начинается на вещественной положительной полуоси с бесконечно большого значения. Положительный охват соответствует повороту радиус-вектора относительно точкина уголпротив часовой стрелки (по часовой стрелке – отрицательный охват). Принятые знаки поворота радиус-вектора объясняются тем, что углу, отложенному от вещественной положительной полуоси, в первом квадранте соответствуети, следовательно,, а в четвертом квадрантеи, следовательно,.
Подсчет охватов удобно проводить с помощью правила Цыпкина, согласно которому сумма охватов точки равна сумме переходов с учетом знакавещественной полуоси в соответствии с рис. 19, на котором показаны возможные случаи переходов вещественной полуоси. ЕслиАФЧХ касается вещественной полуоси , то это эквивалентно тому, что она совершает два полуперехода с обратными знаками, в сумме равных нулю.
При прохождении АФЧХ через точку замкнутая система имеет пару чисто мнимых корней. Это следует из условия, которому соответствует уравнение
,
где характеристическое уравнение замкнутой системы имеет пару мнимых корнейпри некотором значении.
Логарифмический аналог критерия Найквиста устанавливается с учетом связи АФЧХ и ЛАХ и ЛФХ передаточной функции .Критерий устойчивости формулируется следующим образом: для устойчивости замкнутой системы необходимо и достаточно, чтобы в области частот , соответствующих положительным значениям,сумма переходов ЛФХ через линии с учетом знаков (положительный – снизу вверх, отрицательными – сверху вниз и т.д.) равнялась , где– числоправых корней характеристического уравнения разомкнутой системы .
Пересечение ЛФХ линий соответствует точкам, принадлежащим вещественной полуосиплоскости АФЧХ.В большинстве случаев, в том числе для рассматриваемых схем САУ, ЛФХ пересекает только одну линию .
На рис. 17 представлены точные ЛАХ и ЛФХ с указанными значениями запасов устойчивости по амплитуде: (при частоте 12,7 рад/с) и фазеград (при частоте 18 рад/с).
На рис. 18 представлена асимптотическая ЛАХ (и точная ЛАХ для сравнения), которая может быть использована в дальнейшем для синтеза последовательной коррекции.
Графики рис. 17, рис. 18 следует представлять в наибольшем масштабе (растянуть и отредактировать график аналогично предыдущему) и вставлять в текст документа Word, выбирая параметры данной страницы "альбомная".
В соответствии с логарифмическим аналогом критерия Найквиста замкнутая система неустойчива. Для устойчивой замкнутой системы значения должны быть положительными.
Критический коэффициент усиления разомкнутой системы , при котором замкнутая система с отрицательной единичной обратной связью находится на границе устойчивости, определяется согласно критерию Найквиста из условия
,
при этом .
Значение критической частоты можно найти с помощью ЛФХ из условияград, при этом для критического коэффициент усилениявыполняется равенство, что может быть использовано при графическом определениипо ЛАХ и ЛФХ.
Для минимально-фазовых передаточных функций условие соответствует границе колебательной устойчивости.
Для рис. 14 при уточненном значении 12,6518 рад/с вычислениеосуществляется с помощью командной строки
>> Kpas=75.5;Kkp=Kpas/abs(freqresp(Wpas,[12.6518]))
В результате получим значение Kkp =35.3756 равное значению, найденномус помощью метода D – разбиения.
Для проверки правильности полученного результата можно также построить годографы движения корней характеристического уравнения замкнутой системы при изменении коэффициента усиления разомкнутой системы с помощью командной строки
>> WKpas=Wpas/Kpas; rlocus(WKpas);
Результат вычислений представлен на рис. 19. С помощью левой клавиши мыши отметить точку пересечения одной из ветвей годографа с мнимой осью комплексной плоскости корней и не отпуская клавишу подвигать точку в левую и правую сторону. При этом в открывшемся окне будут указаны значения коэффициента усиления и соответствующие ему корни (полюса). Так, например, на рис. 19 видно, что при коэффициентев замкнутой системе все корни устойчивые из них одна пара комплексно-сопряженных корней; при,пара комплексно-сопряженных корней становится неустойчивой:.
Таким образом, при замкнутая система находится на границе колебательной устойчивости ипри становится неустойчивой.
Вычислить корни характеристического уравнения замкнутой системы при выбранном коэффициенте усиления можно с помощью командной строки
>> Kkp=35.3756; p=esort(tzero(1+Kkp*Wpas/Kpas))
Рис. 19
в результате выполнения которой получим
p =
-0.0000 +12.6518i
-0.0000 -12.6518i
-2.4977
-31.1754