где Mдо - крутящий момент двигателя на частоте wо при управляющем напряжении lо .
Значения всех производных вычисляются в рабочей точке ( wо , lо ) и, следовательно, являются константами. В соответствии с первым методом малых
отклонений Ляпунова при малых значениях приращений |
|
ω, |
|
l можно пренебречь |
|||||||||||||||||
членами ряда Тейлора с порядком выше первого. |
|
|
|
|
|
|
|
|
|||||||||||||
После подстановки линеаризованных зависимостей (1.2) и (1.3) в исходное |
|||||||||||||||||||||
уравнение (1.1) получим: |
|
|
¶M |
|
|
|
|
¶M |
|
|
|
|
|
|
|
|
|
|
|||
|
dw |
|
|
æ |
ö |
|
æ |
д |
|
ö |
|
|
|
æ |
¶M |
ö |
|
||||
j |
|
|
= M |
до |
+ ç |
|
д |
÷ |
|
Dw - ç |
|
|
÷ |
Dl - M |
cо |
- |
ç |
|
c ÷ |
Dw |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
dt |
|
ç |
|
÷ |
|
ç |
¶l |
|
|
÷ |
|
|
è |
¶w ø |
|
|||||
|
|
|
è |
¶w ø |
0 |
è |
|
|
ø |
0 |
|
|
0 . |
||||||||
Уравнение статики M дo = Mсo |
|
|
|
|
|
|
|
|
|
|
|||||||||||
описывает поведение двигателя в |
|||||||||||||||||||||
установившемся режиме. Оно следует из уравнения динамики (1.1), если приравнять dw
нулю производную dt .
С учетом того, что M дo = Mсo , проведем сокращения в уравнении (1.4) и
сгруппируем в его левой части члены с |
|
|
ω : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
dw |
é |
¶ |
|
|
ö |
|
æ ¶M |
|
|
ö ù |
æ ¶M ö |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
j |
æ |
M |
c |
|
- ç |
|
|
|
д |
÷ |
|
úDw = ç |
|
|
д |
÷ |
|
|
Dl |
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
+ êç |
|
|
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
dt |
|
|
|
|
¶w |
|
|
¶l |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
è |
¶w |
ø |
0 |
ç |
|
÷ |
|
ú |
ç |
|
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
ê |
|
|
|
|
|
è |
|
|
|
|
|
ø |
0 |
è |
|
|
ø |
0 |
|
|
|
|
|
|
|
|
|
|
(1.5) |
||||||||
|
|
|
|
ë |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
û |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
В результате имеем линеаризованное уравнение двигателя в абсолютных |
|||||||||||||||||||||||||||||||||||||||||
приращениях (или отклонениях). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Упростим уравнение (1.5). После деления его правой и левой частей на |
|||||||||||||||||||||||||||||||||||||||||
выражение, стоящее в квадратных скобках, получим: |
|
|
|
|
¶M |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
æ |
|
д |
ö |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ç |
|
|
|
÷ |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
|
dw |
|
|
|
|
|
|
|
ç |
|
|
¶l |
|
÷ |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ Dw = |
|
|
|
|
|
è |
|
|
|
ø0 |
|
|
|
|
Dl |
||||||||
|
|
éæ |
¶M |
|
ö |
|
|
æ ¶M |
|
ö ù |
|
|
dt |
|
|
éæ |
¶M |
|
|
ö |
|
æ |
|
¶M |
д |
ö ù |
|||||||||||||||
|
êç |
|
c |
÷ |
|
- ç |
|
д |
÷ |
|
ú |
|
|
|
|
|
|
|
|
êç |
|
|
|
c |
÷ |
- ç |
|
|
÷ |
|
ú |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
¶w |
|
|
ç |
¶w |
÷ |
|
|
|
|
|
|
|
|
|
¶w |
|
|
ç |
|
¶w |
÷ |
|
|
|||||||||||||||||
|
è |
ø |
0 |
|
è |
ø |
|
ú |
|
|
|
|
|
|
|
|
è |
|
ø |
0 |
è |
|
ø |
|
ú |
|
|||||||||||||||
|
ê |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
ê |
|
|
|
|
|
|
|
|
|
|
0 |
|
|||||||||||
|
ë |
|
|
|
|
|
|
|
|
|
|
|
|
û |
|
|
|
|
|
|
|
|
ë |
|
|
|
|
|
|
|
|
|
|
|
|
|
û . |
||||
Каждый член уравнения имеет определенную размерность. Однако для исследования системы желательно получить уравнение:
∙в относительных (нормированных) единицах;
∙с безразмерными коэффициентами или с коэффициентами, имеющими размерность времени в степени, равной порядку производной, при которой стоит данный коэффициент.
Введем относительные изменения для входной (управляющей) и выходной (управляемой) величин:
Dl = z(t) |
Dw = x(t) |
, причем ω = |
ω(t) . |
||||||
lo |
; |
wo |
|
||||||
Для ω при wо = сonst имеем: |
|
|
|
|
|
|
|
|
|
|
dw |
= |
d |
(w + Dw) = |
d(Dw) |
|
|
||
|
|
|
|
|
|||||
|
dt |
|
dt |
o |
|
|
dt . |
|
|
|
|
|
|
|
|
||||
Для перехода в уравнении (1.6) к величинам z(t) и x(t) |
выполним следующие |
||||||||
действия: |
|
|
|
|
|
|
ω величиной wо , для чего |
||
1) пронормируем приращение угловой скорости |
|||||||||
разделим правую и левую части уравнения на wо ; |
|
|
|
||||||
(1.4)
(1.6)
21
2) пронормируем l величиной lо , для чего умножим и разделим правую часть
уравнения на lо . Путем соответствующего группирования членов перейдем к относительным управляющим воздействиям:
j
éæ |
¶M |
c |
ö |
|
æ |
¶M |
д |
ö |
|
ù |
êç |
|
÷ |
|
- ç |
|
÷ |
|
ú |
||
|
|
|
|
|
||||||
¶w |
|
ç |
¶w |
÷ |
|
|||||
è |
ø |
0 |
è |
ø |
|
ú |
||||
ê |
|
|
|
|
|
0 |
||||
ë |
|
|
|
|
|
|
|
|
û |
|
Введем обозначения:
x(t) = Dw ; wо
æ |
Dwö |
|
|
|
|
|
|
|
|
|
æ ¶M |
|
ö |
|
|
l |
о |
|
|
|
|
|
|
||||||||
ç |
|
÷ |
|
|
|
|
|
|
|
|
|
|
ç |
|
|
|
д |
÷ |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
dç |
|
÷ |
|
|
|
Dw |
|
|
|
|
|
ç |
|
¶l |
|
÷ |
|
|
wо |
|
|
|
|
Dl |
|
||||||
è wо |
ø |
|
+ |
= |
|
|
|
è |
|
|
ø0 |
|
|
|
|
|
|
||||||||||||||
|
dt |
|
|
|
wо |
|
éæ |
¶M |
c |
ö |
|
|
|
|
æ |
|
¶M |
д |
ö ù lо |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
êç |
|
|
÷ |
|
- ç |
|
|
|
|
÷ |
|
ú |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
¶w |
|
|
|
|
ç |
|
|
¶w |
÷ |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
êè |
ø0 |
|
|
|
è |
|
|
ø |
0 |
ú |
|
(1.7) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
ë |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
û . |
|||
T = |
|
|
|
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
; |
|
|
|
|
|
|
|
||||
éæ |
¶M |
c |
ö |
|
|
|
æ ¶M |
д |
|
ö ù |
|
|
|
|
|
|
|
||||||||||||||
|
|
êç |
|
|
÷ |
|
|
- ç |
|
|
|
÷ |
|
|
ú |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
¶w |
|
|
|
ç |
¶w |
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
è |
|
ø |
0 |
|
è |
ø |
|
|
ú |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
ê |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
ë |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
û |
|
|
|
|
|
|
|
|
|
|
||||
|
|
æ ¶M |
ö |
|
l |
о |
|
|
|
|
|
|
|||
|
|
ç |
|
|
д |
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
¶l |
|
|
w |
|
|
|
|
|
|
|||
|
|
ç |
|
÷ |
|
|
|
|
|
Dl |
|
||||
k = |
|
è |
|
|
ø0 |
|
|
о |
|
|
|
|
; |
||
éæ |
¶M |
c |
ö |
æ |
¶M |
д |
ö |
|
ù |
lо |
|||||
êç |
|
÷ |
- ç |
|
|
|
÷ |
|
ú |
|
|
||||
|
|
|
|
|
|
|
|
||||||||
¶w |
ç |
|
¶w |
÷ |
|
|
|
||||||||
êè |
ø0 |
è |
|
ø |
0 |
ú |
|
|
|||||||
ë |
|
|
|
|
|
|
|
|
|
|
û |
|
|
||
Вычислим производную x'(t) :
z(t) = Dl . lо
|
|
|
|
æ |
Dw |
ö |
|
|
|
|
|
|
ç |
÷ |
|
|
|
¢ |
d[x(t)] |
|
dç |
|
÷ |
|
|
|
|
è wо |
ø |
|
|
||||
x (t) = |
|
|
= |
|
|
|
|
|
dt |
|
dt |
. |
(1.9) |
||||
|
|
|
||||||
Подставив (1.8) и (1.9) в (1.7), получим дифференциальное линеаризованное уравнение в относительных величинах, связывающее относительные угловую скорость
двигателя x(t) и изменение управляющего воздействия z(t) : |
|
′ |
(1.10) |
Tx (t) + x(t) = kz(t) . |
Коэффициент T называется постоянной времени двигателя (объекта регулирования) и имеет размерность времени (с). Он определяет быстродействие двигателя, т.е. как быстро устанавливается заданная угловая скорость при изменении управляющего воздействия.
Коэффициент k носит название передаточного коэффициента и является безразмерной величиной. Он характеризует зависимость между изменением выходной
величины x(t) и управляющим воздействием z(t) в установившемся режиме. Следует подчеркнуть, что дифференциальное уравнение (1.10) описывает
динамику линеаризованной системы управления двигателем в окрестности рабочей
точки ( wо , lо ). Если в исходной нелинейной системе (1.1) изменится рабочая точка, то
это вызовет изменение коэффициентов T и k .
Рассмотрим геометрический смысл линеаризации (рис.1.16). Линеаризация нелинейной зависимости двух переменных означает:
∙замену исходной кривой AB отрезком ее касательной A'B' в точке O′ , соответствующей установившемуся режиму;
∙параллельный перенос начала координат в точку O′ .
22
В дальнейшем будем рассматривать стационарные линейные динамические системы (ЛДС), т.е. системы, которые представляются следующим дифференциальным уравнением:
|
|
an x |
(n) |
(t) + an−1x |
(n−1) |
|
&& |
& |
|
|||||
|
|
|
|
|
|
(t) + ... + a2 x(t) |
+ a1x(t) + ao x(t) = |
|
||||||
|
|
= bm z |
(m) |
(t) + bm−1z |
(m−1) |
t |
& |
(1.11 |
||||||
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
(t) + ... + b2 z(t) + b1z(t) + bo z(t) , |
||||||||
где x(t) и z(t) - соответственно выходной и входной сигналы, причем |
|
|||||||||||||
выполняется условие физической реализуемости n ³ m . |
|
|||||||||||||
В операторной форме дифференциальное уравнение (1.11) имеет вид: |
|
|||||||||||||
A( p)x(t) = B( p)z(t) , |
|
|
|
|
|
|
|
|
|
|
||||
где A( p) , B( p) - дифференциальные операторы левой и правой частей |
|
|||||||||||||
уравнения; |
|
pn + ...a p + a |
|
B( p) = b pm + ...b p + b |
|
|||||||||
A( p) = a |
n |
о ; |
|
|||||||||||
|
1 |
|
|
|
|
|
|
n |
1 |
о ; |
|
|||
p = dpd - оператор дифференцирования.
Mc Mc
B 
B
O
Δω
A
A
0 |
|
|
ω |
||
Рис.1.16.Графическая интерпретация процесса
линеаризации динамической системы
Для математического описания ЛДС широко используются частотные характеристики, позволяющие судить о реакции системы на входные гармонические воздействия различных частот.
Пусть на вход ЛДС поступает гармоническое воздействие в виде бесконечной косинусоиды: z(t) = Z cos(ωt) . Определим x(t) - реакцию системы на этот сигнал.
Известно, что согласно формуле Эйлера |
e+ jωt = cosωt + j sin ωt , поэтому |
|||||||||
cosωt = |
e+ jωt + e− jωt |
|
|
|
|
|
|
|
|
|
2 |
и входное воздействие можно записать в следующем виде: |
|||||||||
|
||||||||||
|
|
z(t) = Z cosωt = |
Z |
e+ jωt + |
Z |
e− jωt = z (t) + z |
2 |
(t) |
||
|
|
|
|
|||||||
|
|
2 |
|
2 |
1 |
. |
||||
|
|
|
|
|
||||||
Для всех ЛДС справедлив принцип суперпозиции: реакция системы на
несколько одновременных входных воздействий равна сумме реакций на каждое воздействие (рис.1.17).
23
|
|
z1(t) |
|
|
|
|
|
x1(t) |
|
|
|
|
|
|
|
|
|
Линейная |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
z(t) |
|
z2(t) |
|
|
|
|
|
x2(t) |
|
|
|
|
|
|
x(t) |
|
|
|
|
динамическая |
|
|
|
|
|
|
|
||||||
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
... |
|
|
|
|
система |
... |
|
|
|
|
|
|
|
||
|
|
zn(t) |
|
|
|
|
|
xn(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
||||||||
Рис.1.17. К пояснению принципа суперпозиции
Отсюда следует практический вывод: если на вход ЛДС поступает сложный сигнал, который можно представить в виде совокупности нескольких простых сигналов
zi (t) , то результирующим сигналом на выходе системы будет сумма реакций этой системы на каждую составляющую xi (t) .
|
|
|
N |
|
Другими словами, если на вход ЛДС подается сигнал |
z(t) = å zi (t) |
, то на ее |
||
i=1 |
||||
|
N |
|
|
|
|
x(t) = å x (t) |
|
|
|
выходе будет сигнал |
i=1 i |
, где N - количество компонентов, на которые |
||
раскладывается входной сигнал z(t) . Для частного случая, когда входной сигнал имеет
вид z(t) = z1(t) + z2 (t) , выходной сигнал вычисляется по формуле: xi (t) = x1(t) + x2 (t) . Определим эффекты, создаваемые каждым из двух экспоненциальных воздействий, на которые разложен исходный входной сигнал.
Для ЛДС справедливо следующее правило: реакцию системы на входное
воздействие можно вычислить путем умножения этого сигнала на функцию параметров системы:
x1(t) = z1(t)W ,
где W - некоторая функция параметров системы. Например, для радиоэлектронной системы функция W зависит от номиналов резисторов, конденсаторов, индуктивностей, составляющих данную систему. При этом для стационарной системы W не является функцией времени (W ¹ W(t)).
Продифференцируем m раз входной сигнал
x (t) = |
Z |
e+ jwtW |
|
|
|
|
|
|
|
2 |
, в результате получим: |
||||||||
сигнал 1 |
|
||||||||
|
|
|
& |
(t) = |
Z |
( jω)e |
+ jwt |
; |
|
|
|
|
|
|
|||||
|
|
|
z1 |
2 |
|
||||
|
|
|
|
|
|
|
|
||
&z&1(t) = Z2 ( jω)2 e+ jwt ;
….
z (t) = Z e+ jwt |
|
|
1 |
2 |
и n раз выходной |
|
||
x&1(t) = Z2 ( jω)e+ jwtW ; &x&1(t) = Z2 ( jω)2 e+ jwtW;
….
z (m) |
(t) = |
Z |
( jω)m e+ jwt ; x(n) (t) = |
Z |
( jω)n e+ jwtW. |
|
|
||||
1 |
2 |
1 |
2 |
|
|
|
|
|
|||
После подстановки выражений для входного и выходного сигналов и их производных в дифференциальное уравнение (1.11) получим:
[a |
n |
( jω)n + a |
n-1 |
( jω)n-1 +... + a |
2 |
( jω)2 + a ( jω) + a |
0 |
] |
Z |
e+ jwt |
× W = |
|||||||||||||||||
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|
|
|
|||||||
= |
[b |
( jω)m + b |
|
|
|
( jω)m-1 |
+ ...+ b ( jω)2 |
+ b ( jω) + b |
|
] |
Z |
e+ jwt . |
|
|||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||
Откуда |
|
m |
|
|
m-1 |
|
|
|
|
|
|
2 |
|
1 |
|
0 |
2 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W =W ( jω) = |
b |
|
|
( jω)m + b |
|
|
( jω)m-1 +...+ b ( jω)2 + b ( jω) + b |
, |
||||||||||||||||||||
m |
|
|
m-1 |
|
|
|
|
2 |
|
|
|
1 |
|
|
0 |
|||||||||||||
|
( jω)n + a |
|
|
|
|
|
|
+ a ( jω) + a |
|
|||||||||||||||||||
|
|
|
|
|
a |
n |
n-1 |
( jω)n-1 +...+ a |
2 |
( jω)2 |
0 |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|||||||||
24
причем для реальных динамических систем m ≤ n .
Если в числителе и знаменателе раскрыть скобки, то получим:
W ( jω) = A(ω) + jB(ω) C(ω) + jD(ω) ,
где A(ω) = bо − b2 ω2 + b4 ω4 − b6 ω6 + ...; B(ω) = b1ω − b3ω3 + b5ω5...; C(ω) = aо − a2ω2 + a4ω4 − a6ω6 + ...; D(ω) = a1ω − a3ω3 + a5ω5...,
причем A(ω) и C(ω) - четные функции; B(ω) и D(ω) - нечетные.
W ( jω) носит название амплитудно-фазовой частотной характеристики или просто амплитудно-фазовой характеристики (АФХ), иногда ее называют комплексной
частотной характеристикой (КЧХ). Неправильно называть W ( jω) передаточной функцией, так как данный термин закреплен в ТАУ за другим важным понятием.
С целью упрощения записи формул введем обозначения: A = A(ω) ; B = В(ω) ; C = C(ω) ; D = D(ω). Тогда имеем:
W ( jω) = |
A + jB |
= |
(A + jB)(C − jD) |
= |
||
C + jD |
(C + jD)(C − jD) |
|||||
|
|
|
||||
= AC + BD |
+ j BC − AD = U (ω) + jV (ω) = H (ω)e jϕ(ω) , |
|||||
C2 + D2 |
C2 + D2 |
|
||||
где U (ω) = Re[W ( jω)] - вещественная частотная характеристика ЛДС; V (ω) = Im[W ( jω)] - мнимая частотная характеристика.
Функция U (ω) является четной в силу того, что
U (−ω) = Re[W (− jω)]= Re[W ( jω)]= U (ω) .
Здесь используются следующие соотношения:
∙произведение двух четных функций является четной функцией;
∙произведение двух нечетных функций является четной функцией;
∙произведение нечетной и четной функции дает нечетную функцию.
Функция V (ω) - нечетная, так как
V (−ω) = Im[W (− jω)]= − Im[W ( jω)]= −V (ω) .
В полярных координатах на комплексной плоскости АФХ имеет:
W ( jω) = H (ω)e jϕ(ω).
Величина H (ω) называется амплитудно-частотной характеристикой (АЧХ) и
|
|
|
|
определяется по формуле: H (ω) = |
|
U 2 (ω) + V 2 (ω) . Так как H (ω) = H (−ω) , то АЧХ |
|
является четной. |
|
|
|
Функция ϕ(ω) называется фазово-частотной характеристикой (ФЧХ) |
|||
|
ϕ(ω) = arctg V (ω) |
||
динамической системы. При этом |
|
U (ω) , если U (ω) > 0 (первый и четвертый |
|
квадранты). Формулы для вычисления ϕ(ω) для различных квадрантов представлены на рис.1.18.
25