W (s) = |
|
KT1 |
|
|
1 |
|
+ |
|
KT2 |
|
|
1 |
, |
||
T |
|
|
T s +1 |
T |
|
|
T s +1 |
||||||||
|
|
−T |
|
|
|
−T |
|
|
|||||||
|
1 |
|
2 |
|
1 |
|
|
2 |
|
1 |
|
2 |
|
||
где коэффициенты перед дробями являются вычетами передаточной функции, вычисленными в соответствующих полюсах. Получаем представление передаточной функции в виде параллельного соединения двух апериодических звеньев, как показано на рис. 3.1.
K |
1 |
T1 −T2 |
s |
u |
1 |
|
y1 |
|||
|
y |
|||||
|
|
T1 |
||||
|
|
|
|
|
||
|
K |
1 |
|
y2 |
||
|
|
|
||||
|
T2 −T1 |
|
|
s |
|
|
|
|
1 |
|
|
||
|
|
|
T2 |
|
|
|
Рис. 3.1
Реакция на ступенчатое входное воздействие системы представляет собой сумму реакций апериодических звеньев, находящихся в параллельных ветвях соединения:
|
KT |
|
|
− |
t |
|
|
KT |
|
|
− |
t |
|
||
|
|
T |
|
|
T |
||||||||||
h(t) = |
|
1 |
|
1 − e |
1 |
|
+ |
|
2 |
|
1 − e |
2 |
|
||
|
|
|
|
|
|
|
|
|
. |
||||||
T |
−T |
|
|
T |
−T |
||||||||||
|
|
|
|
|
|
|
|
|
|||||||
|
1 |
2 |
|
|
|
|
|
|
2 |
1 |
|
|
|
|
|
Упрощая последнее выражение, получим
|
|
|
|
|
T |
|
−T |
|
|
T |
|
−T |
|
||
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
t |
|
h(t) = K |
|
− |
|
1 |
e |
1 |
+ |
|
2 |
e |
2 |
||||
1 |
|
|
|
|
|
|
|
|
|
. |
|||||
T |
−T |
|
|
T |
−T |
|
|
||||||||
|
|
|
|
1 |
2 |
|
|
|
2 |
1 |
|
|
|
|
|
На рис. 3.2 представлены графики переходных процессов в параллельных ветвях и их сумма для случая T1 > T2.
Из рисунка видно, что суммарный график h(t) имеет точку перегиба и представляет собой монотонно возрастающую функцию от нуля до установившегося значения, равного коэффициенту усиления K. Отметим, что в начальный момент времени не только выходной
41
сигнал равен нулю, но и его первая производная по времени также равна нулю, поскольку наклоны касательных в нуле к переходным характеристикам в ветвях параллельного соединения одинаковые по величине и противоположные по знаку. Это соответствует принятому соглашению о нулевых начальных условиях по умолчанию.
h(t) |
|
|
|
KT1 |
|
|
|
T1 −T2 |
|
|
h1 |
K |
|
|
h |
|
|
|
|
0 |
T2 |
T1 |
t |
|
|||
KT2 |
|
|
h2 |
|
|
|
|
T2 −T1 |
|
|
|
|
|
|
Рис. 3.2 |
Структурная модель в виде последовательного соединения
Для получения логарифмических частотных характеристик удобно использовать метод разбиения передаточной функции на простые множители, что означает представление ее в виде произведения типовых динамических звеньев. Именно в виде произведения двух апериодических звеньев и задана передаточная функция системы в рассматриваемом примере. Таким образом, структурная модель принимает вид последовательного соединения двух апериодических звеньев на рис. 3.3, причем не имеет значения порядок их соединения.
u K |
1 |
|
1 |
1 |
y |
|||
|
T1 |
|
s |
T2 |
|
s |
|
|
|
|
1 |
|
|
1 |
|
||
|
|
|
T1 |
|
|
|
T2 |
|
Рис. 3.3
42
Логарифмические амплитудную и фазовую частотные характеристики последовательного соединения получают путем суммирования соответственно амплитуд и фаз всех элементов соединения.
Результирующие ЛАФЧХ и годограф передаточной функции для рассматриваемого примера приведены на рис. 3.4, а и б.
|
L, дБ |
L(ω) |
|
|
θ, град |
|
20lgK |
-20 |
1 |
|
|
|
T2 |
|
|
||
|
0 |
|
0 |
||
|
1 |
|
|
||
|
|
|
|
|
|
а) |
θ(ω) |
T1 |
|
-40 |
-90 |
|
|
|
|
|
|
|
-40 |
|
|
|
-180 |
|
Im |
пл. W ( jω) |
|
|
|
|
|
|
ω →∞ |
|
|
б) |
|
ω = |
0 |
|
|
|
|
-1 |
1 |
K |
Re |
Рис. 3.4 |
|
|||
В области низких частот ω ≤ |
0.1 |
сдвиг по фазе между выход- |
||
T |
|
|||
|
|
|||
|
1 |
|
|
|
ной и входной синусоидами практически отсутствует, а отношение амплитуды выходной синусоиды к амплитуде входной равно коэф-
фициенту усиления K. В области высоких частот ω ≥ 10 синусоиды
T2
43
находятся почти в противофазе, а амплитуда на выходе во много раз меньше амплитуды входного сигнала.
Годограф передаточной функции W ( jω) на рис. 3.4, б также
показывает изменения амплитуды и фазы с ростом частоты от нуля до бесконечности.
Структурная модель системы с обратными связями по переменным состояния
Такой способ представления структурной математической модели известен как метод приведения к канонической форме.
Для системы второго порядка структурная схема показана на рис. 3.5, где координаты вектора состояния x1 и x2 выбраны в качестве выходных сигналов интеграторов.
u |
c |
+ |
– |
1 x2 |
1 x1 y |
|
s |
s |
|||
|
|
|
|
a
b
Рис. 3.5
Замыкая дважды контуры с обратными связями, находим выражение передаточной функции
W (s) = c . s2 + as +b
Разделив числитель и знаменатель этого выражения на коэффициент b, приведем передаточную функцию к виду колебательного устойчивого динамического звена
W (s) = K , T 2 s2 + 2ξTs +1
44
где коэффициент усиления K = bc , постоянная времени T = 1b , а
безразмерный коэффициент затухания ξ = ab . Корни характери-
стического уравнения колебательного устойчивого звена являются комплексно-сопряженными с отрицательной действительной частью, а это будет лишь в случае 0 <ξ <1.
Аналитическое выражение переходной характеристики колебательного устойчивого звена имеет вид
|
|
|
1 |
|
ξ |
|
1 − |
ξ |
2 |
|
1 - ξ |
2 |
|
|
|
|
|
|
|
−T t |
|
|
|
|
|
|
|
||||
h(t) = K |
1 |
− |
1 −ξ |
2 e |
sin |
T |
|
|
t + arctg ξ |
|
|
, |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где действительная часть корня − Tξ определяет затухание колеба-
ний, а мнимая часть |
1 −ξ2 |
частоту колебаний. |
|
T |
|||
|
|
На рис. 3.6 представлен график этого переходного процесса, на котором показаны установившееся значение hуст, два первых максимальных значения hm1 и hm2, разделенные по времени на период колебания Tк.
Указанные значения будут использованы в работе для идентификации параметров модели колебательного звена.
h
hm1
hm2
hуст
Tк |
t |
Рис. 3.6
45