Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Рогов Е.Ф. Основы теории автоматического регулирования учебное пособие

.pdf
Скачиваний:
24
Добавлен:
30.10.2023
Размер:
12.3 Mб
Скачать

Пусть

M /i( 0 1 - f - /i ( 0 и i: f/s (0 1 -f*/a (O'-

Т о г д а

l I / , ( 0 1 L \f, ( 0 1

j / x (t - x ) / a ( 0 r f t =

f f x( x ) / 2 ( * - x) rfx.

 

о

0

Значит, здесь произведению двух изображений будет соответство­ вать не произведение оригиналов, а их свертка.

6 . Предельные свойства преобразований Лаплрса. Если

то

1 / ( 0 1 « = о =

! р / 7 0 > ) ! р = - ,

( * )

 

1 / ( 0 ! * = « =

| p / =' ( p ) | p = 0 J

( * * )

т. е. предельные значения функции

/ ( 0

можно найти, подста­

вив в ее изображение, умноженное на р,

значение р = 0 или р = со.

Формулой (**) можно воспользоваться при определении величины статической ошибки или установившегося значения регулируемой величины системы.

Изображения

основных

функций,

встречающихся в задачах

регулирования,

приведены

в

таблице:

 

 

 

 

 

О р и ги н а л

 

И зо б р а ж е н и е

 

 

II

 

 

 

A

 

 

 

 

 

 

 

P

 

а

w

_ r f : 1

W

 

 

1

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

sin

0) t

 

 

 

0)

 

 

 

 

-|_

ш 2

 

 

 

 

 

 

 

COS <ii t

 

 

 

P

 

 

 

 

 

 

 

p 2 - j -

(O2

 

t"

 

 

 

n\

 

 

 

 

 

 

pn 11

 

 

 

 

 

 

 

e

—at

a >

0

1

 

 

 

 

п р и

P +

 

 

 

 

 

 

 

 

at

 

 

 

 

CO

 

e

 

sin ы t

 

 

( p -f - a ) * -f~

e

—at

 

 

 

p +

«

 

cos a t

 

 

(p - j -

a ) 2

- J - со2

362

§ 2. Решение линейного дифференциального уравнения операционным методом

Рассмотрим дифференциальное уравнение с правой частью:

 

d n x ( t )

я,

+

••• +

«„ А' ( 0

— у (t) .

 

го " dp

' +

Умножим левую и правую части уравнения на

и проинтег­

рируем их

в

пределах от 0

д о о о ;

 

 

 

 

 

dnx(t)

 

rf" - 1 x(t)

• +

anx{t)

e pt dt = j у (t) e

pt d t

U"

di"

A' lh

dt"

 

 

 

 

о

 

п

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

Г dn x (t)

 

 

 

 

 

 

 

 

 

ao I

dt" ~ e

pt d 14- a

 

P‘^

+

--

A(Oe

ptdt =

о

=>'(0 e~pt^ .

0

На основании прямого преобразования Лапласа предыдущее уравнение можно записать в виде

a0L

dn х (t)

+ a\L

da- l x(t)

+ ... -\ - a„L\x(t)]=L\y{t)\ .

(2)

 

~dF~

 

dr- - '

 

 

(Обычно сразу к каждому члену уравнения (1) применяют прямое преобразование Лапласа, опуская промежуточные преобразова­

ния, и получают уравнение 2 ).

 

 

изображе­

Теперь заменим по приведенным выше формулам

ния производных

через изображения

первообразной

функции и

соответствующие

начальные условия:

 

 

 

a0pnL [а ( 0 1 а0 pn-i х (0 )

dx(0)

 

рп~ 2

 

 

+ агр"-1 L [а (01 ~ «г рп~2 х (0 )

j - рп~ 8

^ а (О)

 

dt

 

 

 

 

 

”*■••• + ал ^ I*(01 —L\y (01•

363

Если все выражения, содержащие начальные условия,

перенести

в правую часть, а в левой части вынести затем L [лг(0 |

за скобки,

то получим

 

 

ia0 р« + а, р" - 1 4 -.,. -f а„)Цх (01 = L [у (О] + г (р),

(3)

где г (р) =c1pn~1-j-.. ,-j- c„-ip-j-c„~—сумма всех членов, содержащих начальные условия.

Коэффициенты многочлена г(р) зависят от начальных условий и равны

с, = а„ х (П).

dx (0 )

 

 

 

 

(4)

"

d»-'x_{0 )

+ Й1

+ • • • + a„-i х (0 ).

 

0 dt*-1

 

 

 

 

Если все

начальные

условия

нулевые:

 

 

_<1х\0 )

dl х {0 )

_

__

 

 

~

dt

dt2

dt"—' ~

то г ( р ) ~ 0 ,

Следовательно, уравнение

(3)

примет

вид

где

 

 

Я (р )/.[* (0 ] =

М.У(0 ],

 

 

(5)

 

 

 

 

 

 

 

 

 

 

 

 

Rip) — 4 P a -i « 1 р" - 1 +

. .. + а„.

 

 

Решим уравнение (3)

относительно

L [-т(0].

тогда

 

 

 

L\x(t) }:

( 0 1

, г { р )

 

 

(6 )

 

 

 

 

RiP)

^ Rip) '

 

 

 

Теперь

необходимо

перейти

от

полученного

 

изображения

L [х(0 ]

к оригиналу

.с(0- Положим

y{t) = 0,

тогда из

уравнения

 

(6 )

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7)

364

Найдем функцию x(t),

для которой

 

 

L 1 ( 0 1

г (р) _ И Р " '1 4- са Рп~ 2 4 • ■■4- С„

(8)

' я (Р)

а0Р" 4 аг

4 • • • 4 ап

 

 

С этой целью разложим правую часть равенства на простые дроби

(рассматривается случай,

когда нет кратных корней). Тогда

 

 

L и о |

 

 

 

г(р)

 

 

 

 

 

'(P - P i)(P “

 

Ра) •••(*> — Рп)'

 

 

 

 

 

 

 

 

 

. _____ I

_ 4 2 __(_

 

I

4 я

(9)

 

 

Р

 

Р\ ‘

Р

Рг

"

 

Р — .

 

 

 

 

 

где ри р2,->;Рп— корни

характеристического уравнения/?(р)=0 ;

 

А1г Л2 ,

А„ — коэффициенты

(числители) простых дробей.

Коэффициенты

Лк

вычисляются

 

по

формуле

 

 

 

 

 

 

Лк =

г (рО

 

 

 

( 10)

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

Я'(Рк) ’

 

 

 

 

 

 

 

 

d R {p)

 

 

 

 

 

 

 

(Рк)

 

 

 

 

 

 

R '

 

d р

р

Рк

 

 

 

 

 

 

 

 

 

 

Подставляя

в (9)

эти

значения

Лк, находим

выражения для

L\x

(0 1 :

 

 

 

 

 

 

 

 

 

 

 

 

 

_

г(р) ___

г(рг)

__ i____.

r ip*)

 

 

L [л- (П

 

 

 

 

 

p -

pi

R’(p2)

P — P2

 

 

~

R

( p )

R ’ i P

i )

 

 

 

 

 

 

Г_[РЯ)_ .

 

 

 

 

 

 

 

 

 

R '

iPn)

P - Pn ’

 

И Л

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L [x(0|

 

V

 

? (Pk)

 

( H )

 

 

 

 

Й /?'(Pk)(p-pk)

 

 

 

 

 

 

 

 

Если учесть, ч т о -----------»-ет

 

,

То для изображения L \x{t)j

из

предыдущей

 

Р -

Рк'

'

 

 

 

 

 

x(t)\

формулы

получим ее оригинал

 

 

 

 

 

 

П

f(Pk)

Рк*

 

(12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/?'.(Рк)

365

 

Следовательно,

для нахождения оригинала необходимо:

 

1)

определить

R(p);

 

 

R(p) — 0;

 

2)

найти

корни уравнения

по

3)

определить коэффициенты многочлена ?(р)> вычислив их

формулам

(4);

 

 

R' (р);

 

 

4)

найти

производную

 

 

5)

подставить последовательно ръ р2, ... , рп в R' (р) и в г(р)

и

наити значения

г(р к)

 

(где Л = 1 , 2 , Я

-р-——;

 

 

6 )

 

 

R Ы

( 1 2 ).

 

 

 

составить сумму

 

1

 

Если среди корней имеются комплексные сопряженные:

то

 

 

 

Pk=

*-f/'S,

pk + 1 = a —/3,

 

 

 

 

 

 

 

 

 

 

 

 

r (Pk)

_'>~f~ J Ti__'

i ,•

 

 

 

 

я - м

 

~

, + у з

 

 

 

 

'(Pk+i)_

■Jri

- j 3

 

 

 

 

(Pk+1)

 

 

 

 

 

 

 

.Запишем комплексные числа в векторной форме:

'(Р О

С

/?'(Рк)_

r (Pfc+i)

- ,Л с Ну,

/ ?

' ( Р к + 1

)

 

где

 

 

<р= arctg 4-.

Л =

]/53 -f о2 ;

Среди слагаемых в

У r(pk) Pkt

ЙR'lPk)

будут содержаться слагаемые

ЛеУге(*-МЗН + Л е-f? ек1—

Воспользовавшись

тождеством

Эйлера

ejz = cos z -j- j sin z ,

получим формулу для определения оригинала

функции по се изо­

366

бражению (для случая, когда характеристическое уравнение имеет комплексные корни):

 

п —2S

s

 

x(t) =

Z c j ' k -f- ^ 2 Лк c"ktcos (8к ^ -f- 9к) ,

(13)

где (п 2 ,s') —

число

действительных корней характеристиче­

 

ского

уравнения;

 

s— число нар сопряженных комплексных корней ха­ рактеристического уравнения;

 

Ск

 

г ( Р к )

Лк — 1 о2 +

a2;

-зk = ак tg ~ ;

 

 

 

R' ы

 

 

 

 

 

 

о

 

 

рк — действительный корень характеристического урав­

^■к и

нения;

 

 

действительная и

мнимая

части

соответственно

 

 

 

комплексных корней характеристического уравне­

о и

 

 

ния;

 

 

действительная и

мнимая

части

з — соответственно

 

 

 

выражения

'(Pk)

в случае,

когда

Рк яв-

 

 

R' М

 

 

ляется

 

 

 

 

 

 

 

 

комплексным корнем.

 

 

 

Пусть теперь в уравнении (6 ) у (t) щ 0,

но

все начальные

условия равны нулю и в силу

этого г (р) =

0 .

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L 1x{t)}

L[y{t)]

 

 

 

 

(14)

 

 

 

Rip)

'

 

 

 

 

 

 

 

 

 

 

 

 

В теории регулирования в качестве функций

у (t)

применяются

сигналы,

изображение

которых

является

дрооно-рациональными

функциями вида

М{р)

где М ip) и Т(р)

полиномы от р.

 

 

 

 

Т( Р)’

 

 

 

 

 

 

 

Таким образом, вместо уравнения (14) получим

 

 

 

L\xit)\

Mjp)

. .

Mjp)

 

 

(14а)

 

 

Т ip) Rip)

S{ P)

 

 

 

 

 

 

 

где

5 ip) = Tip) Rip).

367

Рассуждая совершенно так же, как н для случая y ( t ) ~ 0 и

г(р) £ 0 , получим искомый

оригинал

x(t):

 

х ( о = 2

М(рь)

е

Pkt

(15)

т

(Pk)

 

 

S '

 

 

 

где pk — корни уравнения S(p)-=0.

Если

в уравнении

S ( p ) ~ 0

есть один нулевой корень,

т. с

5 (р) =

р S,(p) =

'\

а значит

рг —0,

то оригинал x(t)

для

уравнения

(14а)

можно определить как по формуле (15), так

и по

формуле

 

 

 

 

 

 

 

 

 

 

 

 

М (0)

V

М (Pk)

Pkt

(КЗ)

 

 

■ * ( 0 = ^ ( ° )

k ^ P -V (p k )e

 

 

 

где рк

корни

уравнения

S (р) = 0.

 

 

г

 

М (0 )

определяет установившееся значение пере-

Слагаемое

о~7

т

менной

 

с>1

(.и)

 

 

 

 

 

 

х (/) после окончания переходного процесса.

(16)

Когда

Т(р) = р,

S 1(p) = R(p).

Поэтому формула

примет

вид

 

 

 

 

 

 

 

 

 

x(t)

Л ! ( 0 )

у

М ( pk)

-Pkt

(16а)

 

 

лч (»

г ,

/’к А” о>«)

0

 

 

 

 

 

где рк — корни характеристического уравнения замкнутой систе­ мы R(p) 0 .

В случае, если часть корней уравнения S (р) () будут чис­ то мнимыми или комплексными не кратными, а остальные корни— не кратными вещественными, можно применять дли определения оригинала формулу (15), а при одном нулевом корне — формулу (16). В результате вычислении среди слагаемых полученного ори­ гинала x(t) будут слагаемые, Содержащие комплексные числа. Для получения окончательного решения необходимо от алгебраи­ ческой формы комплексного числа в слагаемом оригинала, полу­ ченном по формуле (15) или (16) при подстановке комплексного

корня, перейти к

векторной (показательной) форме

комплексного

числа, как это показано в книге

|23| при решении задачи № 142.

При наличии комплексных корней в уравнении 5 (/?) = ( 1 для

определения оригинала

x{t)

можно пользоваться

и специаль­

ными формулами разложения, например формулой (13).

В случае кратных корней в уравнении S(p) — 0

для опреде­

ления

оригинала

х (t)

используются другие формулы разложе­

ния,

приведенные

в книге. [6 ],

ч. I ,

 

Пр и ме р . Найти решение уравнения

 

 

I ю d*x (0

,

о,

dx(t)

3 0 -v(O =

dJdT -- + y{t),

(17)

dt*

dr-

 

'

1,1

dt

 

если

начальные условия

при

^ = 0

равны

 

 

 

 

 

 

 

 

х (0) — Л'0 =

3 и

rfJ?(2 L _ n

1

 

 

а

входной

сигнал

_у(7)=1|/‘|.

 

30,

/;0 =

=

1).

примере а 0=

1 /

 

 

10,

п2 = 31, а3 =

 

Применим к каждому члену уравнения (17)

прямое преобра­

зование Лапласа,

у

результате получим

 

 

 

 

L

сР х (t)

+ 10 L

 

dl x

 

dx (t)

- i- 30 I.

'

 

 

x (7)

 

 

 

" d r

 

 

 

dP

 

dt

 

 

 

 

 

 

 

 

=

 

L

dy{t)

+ l

У it)

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

Выполняя дальнейшие преобразования в этом уравнении и учиты­ вая, что вместо /-[-t(0 ] можно записывать X(р), найдем:

f X { p )

- Р*х(0) + р

d.r ( 0 )

tf+ ( 0 )

} + 10 [ X(p)

 

 

 

dt

+

dP

 

Px{0 ) +

dx{0 )

( 4 31 j p X ( p ) - x (0)

+ 30 Ar(p) = p K(p) + Y{p),

 

dt

 

 

 

 

 

или с учетом начальных условий и входного сигнала:

+ Х(р) 3

4 10 j р*Х(р) - 3 р j + 31 j р Х ( р ) - 3 j -f 30 А (р) =

Решая это уравнение относительно Х(р),

получим

Х (Р) =

у ( Л + 1 ) + + Р )

( 18)

у у io+ у - з у у у з о ’

 

где г (р) = Ър'1+ 30р +

93.

 

369

С учетом уравнения (6 ) запишем уравнение (18) в следующем виде:

 

У(пЛ = _______ р_±1_______ ,

Зр» + 30р + 93_

 

КР} р (f + ю р а■- f з 1р + зо) 'г р» + ю р а + 3i р + зо

или

* (р )

М( р ) , г(р)

 

S i p y R i p ) ’

 

 

 

 

где

Л*(Р) — Р + 1 ;

S(p) = T{p)R(p) = pi f Юр’ -f 31ра + 30р;

Т{р) = р\

R(P) = P' + Юр2 + 31 р Ч- 30;

г(р) = 3 р 2+30р+93.

Для определения оригинала x{t)

по изображению (19) не­

обходимо найти

отдельно

оригинал за

счет действия на систему

входного сигнала

y{t),

т. е.

 

х х( 0

М{р) ■S{p) :

и оригинал за счет учета начальных условии, т. е.

х ч( 0

г(р)

 

■RiP)

 

 

 

Решение уравнения (19), следовательно, будет равно сумме

х х (t),

Хъ(0. Т. е.

 

 

 

x(t) = x 1( 0

+ *s(0 -

(20)

Если в системе начальные условия будут равны нулю, то вто­ рое слагаемое в уравнении (19) исчезает и решение уравнения будет

л: (t) — х х(()•

370

Найдем оригинал х г (t) для изображения

(р) = Щр)

(21)

s (p)

 

по формуле (16), так как корни уравнения S(p) = 0

следующие:

р3 = 0, р2 = — 2, р3 = — 3, р4 = — 5.и среди них ость один нулевой корень.

Определим полиномы, которые необходимо подставлять в фор­

мулу

(16).

Так, М(р) = р +

1; З Д = р 5 1(р)=рЧ-Ш р8+31р2+30р.

Значит, S^p) — ^ у 1 0 /Я +

3! р + 30.

В нашем примере Sx(p)

совпадает с полиномом R(p)-

Определим производную для поли­

нома S 1(p),

т. е.

 

 

 

 

 

d_Sy(p)

=

3ра + 20р + 31 .

 

 

d р

 

 

 

Теперь

подставим значения

корней в формулу (16) и образуем

сумму

для

определения х г (t) :

 

 

M{Pl)

* 1 ( 0

Si(Pi)

,e

P iSi'iPt)

 

М(р*)

pP ^,

M(p9)

pP»«

 

>2 ^ 1' (p2)

^ Р з^ /С Р з )

 

_

J_+ l c~ 2,_ i e “Л . A < T 51

( 22)

3 0 ^ 6

3

' 1 5

 

 

Если в системе начальные условия не нулевые, то необходимо

определить оригинал

x^{t)

для

второго слагаемого

уравнения

(19).

Оригинал

х%(t)

определяется по формуле

(12), так как

корни

уравнения

7?(р) = 0

все вещественные и

не

кратные:

P i = - 2, р2 = — 3, р3 = — 5 .

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

3

г(Ръ)

pkt

 

 

 

 

Хо( 0 = £

 

(23)

 

 

е

 

 

 

 

к~\

RrЫ

 

 

 

Для определения оригинала х 2 (t) имеем:

г (р) = 3 р- 4- 30 р г 93;

R(p)

10pa-f 31р + 30;

dR(p) = 3pa-f-20p-f 31, d р

371

Соседние файлы в папке книги из ГПНТБ