- •Белорусский государственный университет
- •Исходные данные приведены в табл. 1
- •Выразим к через параметр т2.
- •Для получения уравнений состояний в нормальной форме используем дифференциальное уравнение замкнутой системы
- •Таким образом, исходная система не удовлетворяет заданным показателям качества, ее следует скорректировать.
- •Учет требований качества переходного процесса: tnn и σ, запасов устойчивости учитываются при формировании среднечастотной области Lж(ω). Здесь можно воспользоваться графиком (рис. 1.6).
- •Корректирующему звену соответствует функция :
- •Для нахождения переходных характеристик замкнутой системы с корректирующим звеном предварительно сформируем модель в пространстве состояний. Передаточная функция замкнутой системы имеет вид:
- •II исследование линейных импульсных сау
- •На рис.2.5 представлена диаграмма Боде исследуемой дискретно системы, с отмеченными на ней запасами устойчивости по амплитуде и фазе.
- •III исследование нелинейной непрерывной системы автоматического управления
- •Литература
Выразим к через параметр т2.
(1.6)
Зависимость К(Т2) приведена на рис. 1.3
х=[0:0.0001:3]
plot(x,(0.0492+1.096*x+1.047*x.*x)./((0.047+1.047*x).^2))
Рис. 1.3
При заданном параметре T2 находим граничное значение Кгр коэффициента передачи К.
Кгр=К(Т2=0.1)= 7.3554
3. Полагая К = 0.7КГР, записываем аналитическое выражение для ф() =axgW(j), L()= 20lg|W(j)| из W(s) при s = j.
К=0.7Кгр= 5.15
Передаточную функцию разомкнутой системы можно записать в виде
(1.7)
где
Тогда строим графики логарифмических характеристик разомкнутой системы, с помощью MATLAB (оператор bode или margin) Рис. 1.4 а. Предварительно с помощью функции paz=tf([K],[a0 al a2 аЗ 0]) найдем:
Transfer function:
5.15
---------------------------------------
0.0047 s^4 + 0.1517 s^3 + 1.147 s^2 + s
Рис. 1.4а
Строим график АФЧХ с помощью MATLAB (оператор nyquist) рис. 1.4 б для разомкнутой системы.
Рис. 1.4б
Запасы устойчивости по модулю и фазе определяются по логарифмическим характеристикам (см. рис. 1.4 а): на частоте среза ωс определяется запас по фазе — Δ φ, а запас по амплитуде ΔL - на частоте при которой φ (ω) = -р. Таким образом, ΔL ≈ 3.1дБ, Δ φ ≈ 7.41°, что является недостаточным.
4. Величина ошибки по скорости определяется как eск=V1 /K. Для ориентировочной оценки tnn и а следует построить переходной процесс h(t) (оператор step в MATLAB) при v(t) = l[t] и по нему определить tnn и а.
Для получения уравнений состояний в нормальной форме используем дифференциальное уравнение замкнутой системы
D(s)y(t) = Kv(t). Если D(s) = b0s4 + b1s3 + b2s2 + b3s + b4 = 0, то уравнение состояния имеет вид:
(1.8)
Для описания динамических систем в пространстве состояний в Matlab применяются модели подкласса ss, которые основаны на линейных дифференциальных или разностных уравнениях.
Модель непрерывной системы в подклассе ss имеет вид:
dx/dt = Ах + Bv;
y = Cx + Dv,
где: х - вектор состояния; u- вектор входа; у - вектор выхода.
Для формирования моделей в подклассе ss предназначена функция ss
sys =ss(A, В, С, D).
В результате под именем sys получаем ss-объект с числовыми характеристиками в виде четверки матриц {А, В, С, D}, которые должны иметь согласованные размеры. Матрицу D в данном случае полагаем равной 0.
Для построения переходного процесса h(t) воспользуемся оператором step в MATLAB.
еализация функций имеет вид:
sys=ss([0 1 0 0;0 0 1 0;0 0 0 1;-b4/b0 -b3/b0 -b2/b0 -b1/b0],[0 0 0 k/b0]',eye(4),zeros(4,1));
step(sys)
a =
x1 x2 x3 x4
x1 0 1 0 0
x2 0 0 1 0
x3 0 0 0 1
x4 -1096 -212.8 -244 -32.28
b =
u1
x1 0
x2 0
x3 0
x4 1096
c =
x1 x2 x3 x4
y1 1 0 0 0
y2 0 1 0 0
y3 0 0 1 0
y4 0 0 0 1
d =
u1
y1 0
y2 0
y3 0
y4 0
В результате получим графики представленные на рис. 1.5. Нас будет интересовать Out(l).
Величина ошибки по скорости определяется как eск=V1 /K = 2.4/5.15=0.466> eск зад = 0.07
Для ориентировочной оценки tnn и σ следует построить переходной процесс h(t) (оператор step в MATLAB) при v(t)=1[t] и по нему определить tnn и σ. Эти величины из графика Out(l) определяются следующим образом:
(4.2)
Время переходного процесса определяется с учетом следующих соотношений: εуст=v(t)/(l+K), где v(t) = l[t], а К=5.15 - общий коэффициент передачи разомкнутой системы. Тогда εуст=1/(1+5.15)= 0.1626 и следовательно tnn из графика Out(l) tn.n.≈ 20c > tnn зад= 1.5c.
Рис. 1.5