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

книги из ГПНТБ / Абгарян К.А. Матричные и асимптотические методы в теории линейных систем

.pdf
Скачиваний:
97
Добавлен:
25.10.2023
Размер:
17.28 Mб
Скачать

2 8 0

Д И Н А М И Ч Е С К И Е

Х А Р А К Т Е Р И С Т И К И

СИ СТЕ М

 

[ГЛ. X

Значит,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

М?1=

1,

 

M!] =

ö2 ^

,

 

 

 

 

 

 

 

 

= ß4

dA^

/rt

dA^

 

d\„

 

 

 

 

 

 

 

 

dt

 

dt

 

dt

+ ^ ( 1 , - 4 )

 

 

 

 

 

*!? -

1.

 

A Ü W ^ b - ,

 

 

 

t& =

ß‘

dA.,

/ о dA,

 

dAt \

d2A,

 

. .

 

 

 

 

 

Hf

16

Hi

 

rtt

) I

dt2

 

 

 

 

 

 

 

_ .

m[I] _

„3

 

 

 

 

 

 

 

m 12 a ,

m \2 — Ü —J f

 

 

 

 

 

 

 

 

m $ — — ß5

 

 

 

 

 

 

 

 

 

 

 

[0]

= a,

 

[li

 

,

dA,

 

 

 

 

 

Щ2

m22 = — ß3

—jf-

 

 

 

 

 

ПІѴ2

= Q5

 

 

 

 

 

 

 

 

 

 

л[ 0 ] _ л

т [ 1 ] _

 

d t

,

1 І 2 ] _

 

„з

dht

dA,

'

Л ,

A x ,

A l

ß

A l

 

f l

A )

dt

,

j[0]2 _

12,

 

2

— а

d(

,

 

32[2]___— а3

df

di

 

 

 

di

,

 

A

 

A

A

 

 

 

 

 

а

 

 

d X \

dh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где а =

1/(AX— A2).

 

 

 

 

 

 

 

 

 

 

Подставляя найденные величины в выражение

(г =

= 0,1, 2), получаем соответственно

 

 

 

 

 

 

 

g f ’ (Л £) =

а (I) І^ехр

1 'kxdt — exp I X2dlj

,

 

gl1>(^£) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= a ® '

! +

a2(öJ M

| L + fl2 (0 .A W

X

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

X

exp

( (A X — f l - ^ - j dt

 

 

 

1 + f l 2 l

d t

e x p ] ^A2 - V a - dll \ u i } ,

 

§ 7] П О С Т Р О Е Н И Е И М П У Л Ь С Н О Й П Е Р Е Х О Д Н О Й Ф У Н К Ц И И 281

s f

( ',а =

«(?)

k<?

 

 

і

 

 

 

 

 

й _

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

 

 

 

 

 

 

 

 

 

 

Af (f, 0 exp

j [%2 + a

^ + a3- ^ -

dt

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А'2»/, 0 = 1 +

«2 (0 - ^ L

 

+ а2 (0

 

 

-

 

 

 

•ß4(0

(*, (I) -

h, (0)

 

 

-

3

dh Ш j

+

 

 

 

 

+

о2( 0 - ^ - а 2(Ю^к Й ) ^

 

 

 

 

 

 

 

 

 

 

 

 

 

dl

 

 

 

+

</ (0

dM 0

 

л

»1')

 

dXs (t) \ .

а' %

( ()

(^s (0

(0)

dt

 

 

 

dt

 

dt

I

Ч-----ПЛ

 

 

 

 

 

 

1

dr

 

 

 

 

 

 

 

 

(a, s =

1,2;

s+= er).

 

 

 

 

n

 

 

d:o

,

1

dg

 

1

 

 

 

 

 

 

П р и м е р . _ +

T - - 2 L _ _ (7 = w.

 

 

 

 

В данном случае

 

 

 

 

 

 

 

 

 

 

A, = .

5 — 1

 

X0 —

/ 5 + 1

 

а

=

 

 

 

 

2^

 

 

 

 

 

2t

 

 

 

 

 

Учитывая

это,

получаем

 

 

 

 

 

 

 

 

/5-1 3- /5

g!01(Л0 = — =-(*

2 g

2

-t

 

• h

 

 

 

 

; о

 

 

 

/ 0 =

y j - (^

l_0- r

^

1+ß)

g ii2)(/,0=

І 36( / E 1- 3 r

 

P g , + P )

( )

 

 

 

 

/5+1

3+ /5

2

É 2 ),

(ß =

- ^ » 0 , 8 9 5

 

/ 5

/ ;

( ß =

11

: 0,985).

5 / 5

 

 

Для сравнения приведем точное выражение импульсной переходной функции, известное для данного примера:

Как видим, отличие приближенного выражения g(2) (/, 0 от точного Si (/, 0 незначительно.

2 8 2

Д И Н А М И Ч Е С К И Е Х А Р А К Т Е Р И С Т И К И С И С Т Е М

[ГЛ.

X

 

Расчеты по формуле (7.20) приводят практически

к

тем же результатам, некоторое различие имеется только в значениях коэффициентов перед скобкой. Так, например,

согласно

(7.20)

 

 

 

 

~(2) н

м

25 г/Ре1—5

/—

Іо

\

§ 8. Реакция системы на показательное возмущение. Передаточная функция

Реакция системы (по всем выходам) на сигнал в виде по­ казательной функции exp (kt), действующий (на промежут­ ке (—оо, t)) на систему по /-му входу, согласно (3.4) пред­ ставляется в виде

t

* / ( М ) = J

— оо

При замене переменных t t' = s имеем

о о

Xj (к, і) = £ gj (s, t s) ex {t—s)ds.

о

Реакция системы по выходу і на входной сигнал в виде показательной функции, поданный на /-й вход,

xiI(X,t) = wil(k,t)e>-‘,

(8.1)

где

 

со

 

Щі {К 0 = 1gii (s, t — S)e-^sds.

(8.2)

6

 

Функция i&ij (к, t), определенная соотношением (8.2), называется передаточной функцией системы (соответствую­ щей /-му входу и і-му выходу)*). Столбцовая матрица

СО

Wj (k, t) = j gj (s, t s) e-Xsds

о

*) Передаточная функция w u (k,t) определена только в области схо­ димости интеграла в соотношении (8.2); во многих практических случа­ ях (но не всегда) возможно путем аналитического продолжения область определения передаточной функции распространить на всю А-плоскость, исключая некоторые особые точки.

§ 8]

Р Е А К Ц И Я С И С Т Е М Ы НА В О З М У Щ Е Н И Е

2 8 3

представляет передаточные функции системы, отвечающие /-му входу и всем ее выходам. Полный набор передаточных функций системы, имеющей I входов и п выходов, дается п X /-матрицей

№ (М ) = (<М М ) <М М ) ••• М М )),

которая связана с матрицей импульсныхпередаточных функций системыследующими эквивалентными соотноше­ ниями:

оо

Wfi,t) = §G(s,t — s)e-tods,

(8.3а)

О

 

 

t

 

 

W {к, 0 = J

G(/ — г, Г) er-ь «-Wdf.

(8.36)

— о о

 

Согласно (8.1) /-й столбец матрицы W (к, t)

 

W/(k,t) = Xj(k,t)e-Xt,

(8.4)

где Xj (к, t) — решение

дифференциального уравнения

A (t)^ j -

= B(t)x, + hl (t)e>-‘,

(8.5)

отвечающее нулевому состоянию системы (т. е. имеется в виду то решение уравнения (8.5), которое отвечает тривиаль­ ному (нулевому) решению соответствующего однородного уравнения).

Столбцовая матрица w,- (к, t) является решением диффе­ ренциального уравнения

■$Г = іи (0 - tân 1O' + Л“ 1(t) h, (0,

(8.6)

отвечающим нулевому состоянию системы, в чем можно убе­ диться путем подстановки в это уравнение выражения (8.4). Значит, матрица передаточных функций W (к, t) является решением дифференциального уравнения

= [U (t) hE„\ W + A~l (t) H (/),

соответствующим нулевому состоянию системы.

284

Д И Н А М И Ч Е С К И Е Х А Р А К Т Е Р И С Т И К И С И СТ Е М

[ГЛ. X

§ 9. Связь между входными и выходными сигналами системы посредством передаточной функции

Между матрицами входных сигналов и (f) и выходных сигналов X (t) имеет место соотношение (см. (3.4))

і

 

х ( 0 = J G(t — t',t')u(t')di'.

(9.1)

00

Предположим, что к входным сигналам можно применить преобразование Лапласа и L (и) — преобразование Лапла­ са матрицы входных сигналов. Тогда

 

 

 

С+іоо

 

 

 

и ® = - Щ- і L ^ )e u dk

( с > с а)

 

 

 

С — /оо

 

 

(са — абсцисса абсолютной сходимости

 

преобразования

Лапласа).

 

 

 

 

 

Подставим выражение и (і) в (9.1) и поменяем порядок

интегрирования. Получим

 

 

 

 

 

of-too

*

G(t t', t') eu'dt'

 

*(9 = 2пі

с — і о о

j

L (и) dk,

или

 

 

 

 

 

'+ic

 

 

 

 

 

Х{0 =

\

G(t — t', t') е-ь "- r)dt’ L (и) eud\.

2ni

 

c—loo

 

 

 

 

Отсюда

 

 

 

 

 

 

(0 =

- 2 s r

-f

 

 

Преобразование Лапласа обеих частей последнего ра­

венства

приводит к соотношению

 

 

 

 

L(x) = W (X, t) L (и).

 

(9.2)

Из (9.2) следует, что передаточная функция линейной системы есть отношение преобразования Лапласа выходно­ го сигнала к преобразованию Лапласа входного сигнала.

§ ІО] П О С Т Р О Е Н И Е П Е Р Е Д А Т О Ч Н О Й Ф У Н К Ц И И 2 8 5

§ 10. Построение передаточной функции

10.1. Передаточная функция стационарной системы.

Учитывая (2.1) и

(8.36),

в случае

стационарной

системы

(А — const,

Н =

const)

имеем

 

 

 

 

№(Х) =

j G(t— Г)е~хи~п сй' =

J eUline~Ku- n dfA~lH- =

 

— o o

 

 

 

 

 

— oo

 

 

 

 

 

 

 

oo

 

 

 

 

 

 

 

 

 

 

= J

eUse-%sdsA~lH = L(eUs) A - ]H,

или, поскольку

 

 

6

 

 

 

 

 

 

L (eUs) = (XE — U)-'

 

 

 

 

 

 

(C M . § 4 г л .

VII), то

 

 

 

 

 

 

 

 

 

 

W (X)^(XE — U)-xA - xH.

(10.1)

Пусть J

= diag (/, (X,), J2

(X2),

..., Jp (Xp)) — жорданова

форма матрицы U, а К

=

(Кѵ К2>КР) — соответствую­

щая преобразующая матрица,

так

что

 

 

 

и = KJM

/ м

= к ~ 1=/ Мі

 

 

 

 

 

 

V

 

 

 

'■мр

 

Тогда (10.1) можно записать в виде

 

 

 

IV (X) =

21 Ко (XEka -

J ar ' M 0A - lH,

 

причем

 

 

 

а=1

 

 

 

 

 

 

 

 

 

 

 

н.

 

 

ңка~1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nk

(XEka

Ja)

Х

— Х д

( Х -

Х а ) 2

+ ■

ка

( x - x j а

10.2. Передаточная функция нестационарной системы. Рассмотрим некоторые из возможных путей построения пе­ редаточной функции нестационарной системы.

10.2.1. И с п о л ь з о в а н и е

в ы р а ж е н и я и м ­

п у л ь с н о й п е р е х о д н о й

ф у н к ц и и. Учиты­

вая (2.3), из (8.3)

имеем

 

t

Gu (/ — 1\ t’) A

(/') H (V) e Kll~n dt’

\Ѵ(Х, 0 = j

2 8 6

Д И Н А М И Ч Е С К И Е Х А Р А К Т Е Р И С Т И К И

С И С Т Е М

[ГЛ. X

И

 

СО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W (X, t) = J

G0(s, t -

s) A~x (t s) H (t s) e~Xsds,

где

 

0

 

 

 

 

 

 

 

G0 = X ( t ) X - l ( n = X ( t ) X ~ l ( t - s ) ,

 

 

 

 

а X

(t) — фундаментальная

матрица

однородного уравне­

ния

(2.1).

 

 

 

 

 

 

 

Имея

точное

или приближенное

выражение

матрицы

G0, можно, пользуясь этими формулами, довести построение

передаточной функции до конца.

 

 

 

10.2.2.

П р е д с т а в л е н и е

 

п е р е д а т о ч н о й

ф у н к ц и и н е с т а ц и о н а р н о й

с и с т е м ы ч е ­

р е з

п е р е д а т о ч н у ю

ф у н к ц и ю с и с т е м ы

п р и

з а м о р о ж е н н ы х

п а р а м е т р а х .

Импульс­

ную

переходную

функцию

G (s, t s)

разложим в ряд

Тейлора в окрестности точки (s, t):

 

 

 

Заменяя в (8.3а) G (s,

t — s) ее разложением в ряд Тей'

лора,

имеем

 

 

 

 

 

 

«7 (X, о =

СО

 

с о

 

 

 

 

VG(S, о e~Ksds -

J s

dG(a;‘]

e-^ds +

 

 

 

0

 

о

 

 

 

 

о о

о

Здесь

J G (s, t) é~Xsds - R (X, i) = Roo(Я,, t),

о

со

0

и, вообще,

со

[sl

* ° M 4 - e 'u ds =

J

dt1

D

 

§ 101

П О С Т Р О Е Н И Е

П Е Р Е Д А Т О Ч Н О Й

Ф У Н К Ц И И

2 8 7

Учитывая это, получаем

 

 

 

W (X, t) = R (К, t) +

Rn (X, /)+ -!■

R22 (X, t) +

. - • (10.2)

Матрицу R (X, t) молено трактовать как матрицу переда­ точных функций системы в условиях, когда ее параметры в момент времени t заморожены (т. е. имеют постоянные значения, соответствующие моменту времени t).

Если импульсная переходная функция в качестве вто­ рого аргумента имеет не t, а медленное время т = st, то тогда

 

 

со

 

 

W (X, т) =

j o (s, х — es) e~Xsds

 

 

о

 

и разложение (10.2) принимает вид

 

W (X, х) =

R (X, т) +

sRn (X, X) +

S2R22 (X, т) +

10.2.3.

П е р е д а т о ч н а я

ф у н к ц и я к а к р е ­

ш е н и е д и ф ф е р е н ц и а л ь н о г о у р а в н е н и я . Передаточная функция может быть построена и как решение дифференциального уравнения (8.6). Для решения уравне­ ния (8.6) или по крайней мере упрощения этой задачи мож­ но воспользоваться методом расщепления дифференциаль­ ной системы на подсистемы меньшего порядка. С этой целью наряду с (8.6) введем в рассмотрение уравнение

= \Ѵ № - *Еп] W + Л-1 (Т) h, (t)

(X=

st), (10.3)

которое при е = 1 совпадаете (8.6).

U (т)

уравнения

Пусть собственные значения матрицы

(10.3) разбиваются на непересекающиеся группы Х\а), Х ^ ,...

..., X{kJ (а =

1, 2,

...,

р; Уka = п), так что она может быть

представлена

в

форме

 

 

 

 

 

 

 

 

U =

2 КаКМа.

 

 

 

 

 

 

 

 

а=\

 

 

В этих условиях собственные значения матрицы U ХЕ

такл(е

разбиваются

на

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

группы

вида

Х\а)~

X, w

-

X, ....

Х{°1

- Х ( а = 1 ,2 .......р;

Уka =

п), а

288 Д И Н А М И Ч Е С К И Е Х А Р А К Т Е Р И С Т И К И С И С Т Е М [ГЛ . X

сама матрица U %Еп представима в форме

 

U hEn — 2

Ка (Ла — Я,Ека) Ма.

 

а=

1

Формальное решение уравнения (10.3) определяется

соотношениями

 

 

2

K a (i, e)yaj,

 

< 7 = 1

dUai

= [A<J (т, е) — XEka\ ijaj + М„ (т, е) А ' (т) lij (t)

di

 

(0 = 1

, 2 , . .. ,/?).

Удерживая в разложениях матриц Ка, Аа, Ма некото­ рое число первых членов, получим соотношения, определяю­ щие приближенное решение уравнения (10.3):

w f =

2

KP (т>і е) yPh

dtM)

а=в1

т

а] уР, + М р (т, е) А~1(т) h, (t)

= (Аа0 (Т, е) -

(а = 1 , 2 , .. ., р).

Полагая е = 1, отсюда находим соотношения, опреде­ ляющие приближенное решение уравнения (8 .6 ):

w(p

= Ѣ к Р (t) yPh

 

duV)

0 = I

 

Ь £ *а| yp, +

(t) A~' (/) hj (t)

—З Г = [ A (a° (О -

( o = l , 2 , . . ., p).

Допустим, что YP (К t) — невырожденная квадратная матрица порядка k0, удовлетворяющая матричному урав­ нению

т г = |Л<«’ « -

Тогда

t

ypj = J 1/!,, ) ( ^ 0 1 /аГ І ( ^ / ' ) ^ а , ( П ^ ' ( П / і / ( П ^

—ОС

iS 10]

ПОСТРОЕНИЕ ПЕРЕДАТОЧНОЙ ФУНКЦИИ

2 8 9

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

w)r) (*-. 0 =

рL

= 2 к Р (t) \

 

(х, t) Y P

(х, t') м Р

Л-

(t')1

h, (t') dt'.

0 =1

 

-CO

 

 

 

 

 

 

 

 

В более компактной форме

 

 

 

 

 

 

w p {К 0

=

<

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= /Си

 

( 0j’

У(ГЧ М )У (0 ~ \М ') М (Ѵ

М

_ У ) М

Г) ^ '.

Здесь

 

—'СО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/Ог) =

(/СГ

/СГ . .. О

,

м ю =

|

••• |,

 

 

 

 

 

 

 

 

\MlpJ

 

 

 

 

У( г 1 =

diag (У<,'\ Y p .........Yp).

 

 

 

В частном случае, когда собственные значения матрицы U простые, разбивая эти собственные значения на п «групп» (по одному собственному значению в каждой «группе»), будем иметь

W(P = 2

Ka' (t) уа],

0

= 1

<4и{г)

 

 

dt

=

( Х Р

(0 -

Х ) у р

- м Р (О л - 1( 0 h ,

( I )

 

а '

 

 

 

 

 

 

 

(о = 1,2,

. . . , п).

 

Тогда

 

 

 

 

 

 

w p

(X, і) =

 

 

 

 

 

 

 

<

 

 

 

=

2 KP

J

exp J (Хр (Г) -

X) dt"Mp (t') A ~ \f) h, (t') dt',

О—1

—oo

t'

 

 

 

 

или, более компактно,

 

 

 

wp(X, t) =

 

 

 

 

 

 

 

t

 

t

 

 

 

 

= Kin (t) J

exp J (A{r\t") XEn) dt"Mir) (t') Л- 1 (t') h, (t') dt',

 

—CO

t'

 

 

 

 

где

 

A<'> =

diag(tf\A ir>..........XP).

 

 

 

 

10 К. А. Абгарян

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