Материал: бархоткин системы автоматического управления

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

где 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.

1.6. Частотные характеристики линейных динамических систем

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 + ejωt

 

 

 

 

 

 

 

 

2

и входное воздействие можно записать в следующем виде:

 

 

 

z(t) = Z cosωt =

Z

e+ jωt +

Z

ejω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