Основы теории автоматического регулирования. В 2 частях. Ч.1. Учебное пособие
.pdf
6.2. Функции и характеристики
Рассмотрим функции и характеристики теоретического анализа систем автоматического управления для трех пространств:
–физического;
–операторного пространства Лапласа;
–частотного.
Под физическим пространством будем понимать временное пространство. Основные положения операторного пространства Лапласа рассмотрены в разделе 6.1. Частотным пространством будем считать пространство анализа установившихся режимов при гармоническом воздействии заданной частоты.
6.2.1. Физическое пространство
Переходная характеристика. Это реакция звена, вызванная подачей на его вход единичного ступенчатого воздействия 1(t) (рис. 6.2). Переходную характеристику обозначают h(t).
Весовая функция, или импульсная переходная характеристика.
Весовая функция представляет собой реакцию звена на единичный импульс δ(t )(рис. 6.3). Весовую функцию обозначают w(t).
Поскольку δ(t )=1′(t ), то w(t)= h′(t). Отсюда следует, что
h(t)= ∫w(t)dt. |
|
6.2.2. Операторное пространство |
|
Характеристикой операторного про- 1(t) |
h(t) |
странства является передаточная функция. |
|
Для определения этой характеристики ис- |
|
ходноедифференциальноеуравнениеопи- |
|
сания звена |
|
|
a |
dnx |
+a |
dn x |
|
+ |
|||
|
|
n |
|
n-1 |
|||||
|
0 |
dt |
1 |
dt |
|
||||
|
|
|
|
|
|
|
|
||
+a x =b |
dm f |
+b |
dm-1 f |
+ +b f |
|||||
|
m |
|
|
m-1 |
|||||
n |
0 |
dt |
1 |
dt |
m |
||||
|
|
|
|
|
|
|
|||
следуетперевестивоператорноепространство. Будем исходить из того, что все
δ(t) w(t)
Рис. 6.3.
К определению весовой функции
41
начальныеусловиянулевые.Этоозначает,чтоеслиполученыизображениявыходногоивходногосигналов x(t) X (p)и f, (t) F (p)
то для производных входного и выходного сигналов будем иметь
dxdt pX (p), dfdt pF (p),
(6.1)
dnx |
pnX (p), |
dm f |
pmF (p). |
dtn |
|
dtm |
|
После подстановки выражений (6.1) в исходное дифференциальное уравнение описания звена получим описание звена в операторном пространстве:
a pnX (p)+a pn−1X (p)+ +a X (p)= |
|
||
0 |
1 |
n |
|
=b0 pmF (p)+b1pm−1F (p)+ +bmF (p). |
(6.2) |
||
Из слагаемых левой части уравнения (6.2) вынесем изображение выходного сигнала, а из слагаемых правой части – изображение входного сигнала. Содержимое круглых скобок обозначим полиномами R(p) и V(p):
(a pn +a1pn−1 + +an )X(p)=(b0 pm +b1pm−1 + +bm )F (p),
0
R(p) |
V (p) |
гдеR(p)–полиномприизображенииреакции;V(p)–полиномпри изображении воздействия.
В итоге
R(p)X(p)=V(p)F(p). (6.3)
Передаточной функцией звена называется отношение изображения выходного сигнала к изображению входного сигнала:
W (p)= |
X (p) |
|
(6.4) |
|
F (p). |
||||
|
||||
Введенная для операторного пространства характеристика звена названа так потому, что по известному изображению входного сигнала позволяет получить изображение выходного сигнала:
X (p)=W (p)F (p). (6.5)
42
Из уравнений (6.4) и (6.3) следует, что
W(p)= |
X (p) V(p) |
(6.6) |
|||
|
|
|
|||
F (p) = |
R(p). |
||||
|
|||||
То, что передаточная функция, наряду с определением, может быть представлена и в виде отношения полиномов V(p) и R(p), является обоснованием совпадения внешнего вида передаточной функции и сокращенной формы записи дифференциального звена, поскольку в соответствии с (5.5) сокращенная форма записи дифференциального полинома определяется как отношение совпадающих по внешнему виду полиномов V(p) и R(p).
Установим связь между передаточной функцией и характеристиками физического пространства.
Будем считать, что на вход звена в физическом пространстве подается δ-функция. С одной стороны, как отмечалось в подразделе 6.2.1, выходным сигналом звена является весовая функция w(t).Сдругойстороны,δ(t ) 1,ивсоответствиис(6.5)X (p)=W (p).
Это означает, что весовая функция – оригинал передаточной функции или передаточная функция есть изображение весовой функции.
6.2.3.Частотное пространство
Ив теории цепей, и в теории автоматического регулирования широко используется пространство установившихся режимов при гармоническом воздействии заданной частоты. Это пространство логично было бы назвать гармоническим, однако его характеристики называют частотными. Поэтому пространство установившихся режимов при гармоническом воздействии заданной частоты будем называть частотным.
Общейчастотнойхарактеристикойявляетсяамплитудно-фазо- вая частотная характеристика звена. Известно, что при наличии
всистеме автоматического управления гармонического воздействия известной частоты в установившемся режиме все переменные имеют вынужденные колебания этой же частоты со своими амплитудой и фазовым сдвигом по отношению к источнику
43
гармонического сигнала. Таким образом, вводя определение ам- плитудно-фазовой частотной характеристики звена, надо обеспечить сохранение такой информации.
Дляопределенияамплитудно-фазовойчастотнойхарактеристи- ки звена исходное дифференциальное уравнение описания звена
a |
dnx |
+a |
dn−1x |
+ +a x =b |
dm f |
+b |
dm−1 f |
+ +b f |
|||||
|
n |
|
n−1 |
|
m |
|
m−1 |
||||||
0 |
dt |
1 |
dt |
n |
0 |
dt |
1 |
dt |
m |
||||
|
|
|
|
|
|
|
|
|
|
||||
следует перевести в частотное пространство.
Будемисходитьизтого,чтонавходзвенапоступаетгармоническоевоздействие f (t )= fm sin (ωt +ψf )(ψf =0)иизображениегар-
монического воздействия на входе звена имеет вид f (t) fmejωt. В установившемся режиме на выходе звена наблюдается гармонический сигнал x (t )= xm sin (ωt +ϕ), а его изображение имеет вид
x(t) xmej(ωt+ϕ).
Поскольку исходное дифференциальное уравнение описания звена содержит слагаемые с производными входного и выходного сигналов, то соответствующие им изображения будут такими:
|
df |
jωfmejωt, ..., |
dm f |
( jω)m fmejωt; |
|
dt |
dtm |
||
|
|
|
||
dx |
jωxmej(ωt+ϕ), ..., |
dnx |
( jω)n xmej(ωt+ϕ). |
|
dt |
|
|
dtn |
|
Подставляя в исходное дифференциальное уравнение описаниязвенаизображениявходногоивыходногосигналов,атакжеих
производных, получим |
|
|
|
|
|
|
|||
a |
( jω)n x |
m |
ej(ωt+ϕ) +a |
( jω)n−1 x |
m |
ej(ωt+ϕ) + +a x |
m |
ej(ωt+ϕ) = |
|
0 |
|
1 |
|
n |
|
|
|||
|
=b0 ( jω)m fmejωt +b1 ( jω)m−1 fmejωt + +bm fmejωt. |
(6.7) |
|||||||
Изкаждогослагаемоголевойчастиуравнения(6.7)вынесемсо-
множитель xmej(ωt+ϕ), а из каждого слагаемого правой части – со- |
||||||||||||||
множитель fmejωt: |
|
|
|
|
|
|
|
ej(ωt+ϕ) |
|
|||||
a |
|
( jω)n +a |
( jω)n−1 |
+ +a |
x |
m |
= |
|||||||
0 |
|
|
1 |
|
n |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
( jω) |
m |
+b1 ( jω) |
m−1 |
+ |
|
|
fme |
jωt |
. |
(6.8) |
||
b0 |
|
|
+bm |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
44
Содержимое квадратных скобок левой части уравнения (6.8) обозначим R(jω), а правой части – V(jω): R( jω)xmej(ωt+ϕ) =
=V ( jω) fmejωt.
Амплитудно-фазовой частотной характеристикой (АФЧХ) называют отношение
W ( jω)= |
V ( jω) |
. |
(6.9) |
|
|||
|
R( jω) |
|
|
Покажем, что данное определение АФЧХ позволяет получить необходимуюинформациюобамплитудеифазовомсдвигевыходного сигнала. Действительно, из (6.9) следует
W ( jω)= |
V ( jω) |
= |
xm |
ej(ωt+ϕ) |
= |
xm |
ejϕ. |
R( jω) |
|
|
|||||
|
|
fm ejωt |
|
fm |
|||
Видно, что амплитудно-фазовая частотная характеристика содержит информацию об изменении амплитуды и фазовом сдвиге в результате прохождения сигнала через звено.
При определенном заданном значении частоты амплитуднофазовая частотная характеристика, по сути, является комплексным числом:
W ( jω)= Amejϕ = Am cosϕ+ jAm sinϕ=W ′(ω)+ jW ′′(ω). (6.10)
Поскольку в выражении для передаточной функции (6.6) и в выражении для амплитудно-фазовой частотной характеристики (6.10) участвуют одни и те же полиномы, то это означает, что способом получения аналитического выражения для АФЧХ будет использование аналитического выражения для передаточной функции: необходимо провести замену оператора p на jω и привести комплексное число к каноническому виду.
Пример 6.2. Пустьзаданоаналитическоевыражениепередаточной функции звена
W (p)=TpK+1.
45
+j |
|
|
|
ω = ∞ |
ϕ(ωk) |
ω = 0 |
|
|
|
|
+1 |
|
|
ωk |
Годограф |
|
Am(ωk) |
|
АФЧХ |
|
|
|
Рис. 6.4. Годограф амплитуднофазовой частотной характеристики
Требуется найти аналитическое выражение для АФЧХ звена:
W ( jω)= |
K |
|
|
= |
|
|
|
K |
|
−T(jω)+1 |
= |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||||||
T(jω)+1 |
T(jω)+1 −T(jω)+1 |
||||||||||||||||
|
|
|
|
|
|
||||||||||||
|
K 1− jωT |
) |
|
|
|
|
|
|
K |
|
KωT |
|
|
||||
( |
|
|
|
|
|
|
|
|
|
|
|||||||
= |
|
|
= |
|
|
− j |
|
. |
|||||||||
(1+ jωT )(1− jωT ) |
|
1+ω2T 2 |
1+ω2T 2 |
||||||||||||||
Поскольку при определенном значении частоты амплитуднофазовая частотная характеристика, по сути, является комплексным числом, то можно предложить графическую интерпретацию АФЧХ. Если частоту монотонно изменять от значения ω = 0 до значения ω = ∞, каждую полученную точку наносить на комплексную плоскость и затем их последовательно соединить, то получим годограф АФЧХ (рис. 6.4).
Помимо амплитудно-фазовой частотной характеристики выделяют амплитудную частотную характеристику А(ω) и фазовую частотную характеристику ϕ(ω):
( ) |
( |
′ |
( )) |
2 |
|
( |
|
′′ |
( |
|
)) |
2 |
|
|
A ω = |
W |
ω |
|
+ W |
|
ω |
|
|
; |
(6.11) |
||||
|
|
|
|
|
|
|||||||||
|
|
|
|
W ′′ |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||
ϕ(ω)=arctg |
|
. |
|
|
|
|
|
|
||||||
|
|
|
|
W ′ |
|
|
|
|
|
|
||||
Первая характеристика (6.11) показывает, как изменяется ам- плитудасигнала,проходящегочереззвено,вторая–какизменяет- сяфазовыйсдвиг.Обехарактеристикиимеютграфическуюинтерпретацию в функции от частоты.
46
6.3.Элементарные типовые звенья
6.3.1.Пропорциональное звено
Рассмотрим функции и характеристики пропорционального звена.
1.Уравнение звена: x(t)= kf(t).
2.Передаточная функция звена: W(p)= k.
3.Весовая функция: w(t). Поскольку весовая функция является оригиналом передаточной функции, то, представляя передаточную функцию в виде W(p)= k = k 1 и используя свойство линейности преобразования Лапласа, для весовой функции получим выражение w(t)= kδ(t).
|
t |
t |
4. |
Переходная характеристика звена: h(t)= ∫w(τ)dτ= ∫kδ(τ)dτ= |
|
t |
0 |
0 |
= k∫δ(τ)dτ= k 1= k (рис. 6.5). |
|
|
0 |
Амплитудно-фазовая частотная характеристика звена: |
|
5. |
||
W(jω)= k ( W ′(ω)= k, W ′′(ω)=0). Графическая |
интерпретация |
|
АФЧХпропорциональногозвена–этоточканавещественнойпо- ложительной полуоси комплексной плоскости.
6. Амплитудная частотная характеристика (АЧХ):
A(ω)= 
W ′2 (ω)+W ′′2 (ω) = k. В процессе прохождения через звено гармонический сигнал любой частоты изменяет амплитуду в k раз.
7. Фазовая частотная |
характеристика |
|
||||
(ФЧХ): |
|
|
|
h(t) |
||
|
′′ |
|
0 |
|
||
|
|
|
|
|||
ϕ(ω)=arctg |
W (ω) |
=arctg |
=0. |
k |
||
′ |
|
|
||||
|
W (ω) |
|
k |
|
||
В процессе прохождения через звено |
|
|||||
гармонический сигнал любой частоты не |
|
|||||
получает фазового сдвига. |
|
|
|
t |
||
6.3.2. Инерционное звено первого порядка |
Рис. 6.5. Переходная |
|||||
Рассмотрим функции и характеристики |
характеристика |
|||||
пропорционального |
||||||
инерционного звена первого порядка. |
звена |
|||||
47
1.Уравнение звена:T dxdt + x = kf(t).
2.Передаточная функция звена:W(p)=Tpk+1.
3.Весовая функция звена:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
||
|
k |
|
|
|
|
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
||||||||
−1 |
−1 |
k |
|
−1 |
k |
|
|
|
|
|
k |
− |
|
t |
|
||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||
w(t)= L |
|
|
|
= L |
|
|
|
|
|
|
|
= L |
|
|
|
|
|
|
|
|
= |
|
|
e |
T |
|
. |
|
|
|
|
1 |
|
|
|
|
|
1 |
|
T |
|
||||||||||||||
|
Tp+1 |
|
|
T |
|
p+ |
|
|
|
T |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
p− |
− |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Получим весовую функцию, используя теорему разложения. Характеристическое уравнение звена, которое получают в результате приравнивания к нулю знаменателя передаточной функции D(λ)=Tλ+1=0, имеет единственный вещественный корень
λ = −T1 .
Определим постоянный коэффициент
C1 = HG(′(λλ)) =Tk .
Имеем
k
W (p)= pT+λ w(t)=Tk eλt =Tk e−T1t .
Видно, что единственный вещественный корень характеристического уравнения определяет экспоненциальную составляющую
|
движения весовой функции (рис. 6.6). |
w(t) |
4. Переходная характеристика звена: |
t
Рис. 6.6. Весовая функция инерционного звена первого порядка
|
|
t |
|
t |
k |
|
− |
1 |
τ |
|
k |
t |
− |
1 |
τ |
|
|
|
|
|
||||
h(t)= ∫w(τ)dτ= ∫ |
e |
T |
dτ= |
∫e |
T |
d |
τ= |
|||||||||||||||||
T |
|
|
|
T |
|
|
||||||||||||||||||
|
|
0 |
|
|
0 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
k |
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
||
= |
(−T)e−T |
τ |
|t0 = −ke−T t |
+k = k 1 |
−e |
−T t |
. |
|||||||||||||||||
|
|
|||||||||||||||||||||||
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Переходная характеристика инерционного звена первого порядка представлена на рис. 6.7.
48
Переходная характеристика может быть получена в процессе проведения эксперимента. Такая характеристика может использоваться при определении типа звена и его параметров.
Параметрами инерционного звена пер- t вогопорядкаявляютсякоэффициентпере-
дачи k и постоянная времени Т.
Если экспериментальная переходная характеристика соответствует переходной характеристикеинерционногозвена,тоопре-
деление параметров звена осуществляется следующим образом:
1)определяют установившееся значение переходной характеристики. Это значение – коэффициент передачи k;
2)на переменной части переходной характеристики регистрируют ее значение h* и момент времени t*. Далее с использованием аналитического описания переходной характеристики инерционного звена определяют его постоянную времени:
|
|
|
|
|
− |
1 |
t* |
|
* |
|
||||
|
|
−e |
|
T |
|
|
|
|
||||||
k 1 |
|
|
|
|
= h |
; |
|
|||||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
t* =1− h* ; |
|
|
|||||||||
|
e− |
|
|
|
||||||||||
T |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
k |
* |
|
||
− |
1 |
t |
* = ln |
1− h |
; |
|||||||||
|
|
|
||||||||||||
|
|
|||||||||||||
T |
|
|
|
|
|
|
|
k |
|
|
||||
|
|
|
|
|
|
|
|
|
||||||
T = − |
|
|
|
t |
|
|
. |
|
||||||
ln |
|
1 |
|
* |
|
|
||||||||
|
|
|
|
|
|
− h |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
5. Амплитудно-фазовая частотная характеристика: |
|
|||||||||||||
W(jω)= |
k |
|
|
= |
k |
|
1− jωT |
= |
k(1− jωT) |
= |
||||
|
|
|
|
|
|
|
|
|||||||
T(jω)+1 |
T(jω)+11− jωT |
1+ω2T 2 |
||||||||||||
|
|
|
|
|
|
|||||||||
|
|
k |
|
|
|
|
kωT |
|
|
′ |
|
′′ |
|
|
= |
|
2 2 |
− j |
2 2 |
|
|
|
|
||||||
|
=W (ω)+ jW (ω). |
|
||||||||||||
1 |
+ω T |
|
|
1+ω T |
|
|
|
|
|
|
||||
49
W′′(ω) |
|
A(ω) |
|
ω = ∞ |
ω = 0 |
|
|
|
W′(ω) |
0 |
ω |
Рис. 6.8. Амплитудно- |
Рис. 6.9. Амплитудная |
||
фазовая частотная |
частотная |
|
|
характеристика |
характеристика |
|
|
инерционного звена |
инерционного звена |
||
первого порядка |
|
|
|
ϕ(ω)
0
− π2
ω
Рис. 6.10. Фазовая частотная характеристика инерционного звена первого порядка
Приизменениичастотыотω = 0 до ω = ∞годографАФЧХбудет иметь такой вид, как представлено на рис. 6.8.
6. Амплитудная частотная характеристика звена (рис. 6.9):
|
|
|
|
|
|
|
|
|
k |
2 |
|
|
|
kωT |
2 |
|
|
|
W ′2 |
(ω)+W ′′2 (ω) = |
|
||||||||||||||
A(ω)= |
|
|
|
|
|
|
+ |
|
|
= |
|||||||
1 |
2 2 |
2 2 |
|||||||||||||||
|
|
|
|
|
|
+ω T |
|
|
|
|
1+ω T |
|
|
||||
|
|
= |
k2 (1+ω2T2 ) |
|
= |
|
|
k |
|
|
|
. |
|
|
|
||
|
|
(1+ω2T 2 )2 |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
1+ω2T 2 |
|
|
|
||||||||||
С увеличением частоты входного гармонического сигнала имеет место уменьшение амплитуды выходного сигнала.
7.Фазовая частотная характеристика звена (рис. 6.10):
−kωT
ϕ(ω)=arctg |
W ′′(ω) |
=arctg |
1+ω2T 2 |
= −arctg(ωT ). |
′ |
k |
|||
|
W (ω) |
|
|
|
1+ω2T 2
С увеличением частоты входного гармонического сигнала имеет место появление фазового сдвига, приближающегося кϕ= − π2 .
6.3.3. Колебательное звено
Рассмотрим функции и характеристики колебательного звена.
1. Уравнение звена: a0 d22x +a1 dx + x = kf(t). dt dt
50
