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

книги / Матричные и асимптотические методы в теории линейных систем

..pdf
Скачиваний:
1
Добавлен:
19.11.2023
Размер:
26.59 Mб
Скачать

t

g l* (< . 0 =

a ©

k?

 

(t, S exp

j‘

 

 

 

 

d l -

 

 

_

k f (t, ?) exp j

(x. + a

^

+ a

^ - ^

где

 

 

 

 

 

 

 

 

 

 

 

 

 

- > 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k f

(t, i) =

1 +

а2 (i)

 

 

+

a2 (t)

 

 

_

 

 

- a‘(I) [(J, ©

-

к

© )

 

-----3 ( - ^

- ) 2] +

 

 

 

+

а

2 0

 

Л<Д0

а(Ю

ж

+

 

 

 

 

 

(

- ^ - а

 

 

+ a1 (0

ЛоЮ

и «МО

 

л , и \

+

<|*х.

 

 

 

Л

■('

 

 

Л

 

 

dt

)

^~

 

dt2-(М О -М О )]

 

 

 

 

 

(a, s = l , 2 ;

s ^ o ) .

 

 

 

n

 

£f"(7 .

 

1

do

 

I

=

 

 

 

 

П р и м е р .

 

+

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

я ,,=

I 5 — 1

 

 

Яд --

 

/ 5 + 1

 

a =

 

 

2/

 

 

 

 

 

 

 

2/

 

 

 

У5 *

Учитывая это,

получаем

 

 

 

 

 

 

 

 

 

/ s - i

з - Уб

- (

Уб+1

з+ УЪ

 

 

(1, м

- V << 2 £ 2

 

2 I 2 ),

 

 

I о

 

 

 

 

 

 

 

 

 

 

 

 

J1)

 

 

 

 

 

 

 

<“е1'+р)

(Р = у ъ = * о,895),

гг да 1) = -p i Л 1_э -

 

«р

(0 1) =■ - p i Л

' -3 -

 

r V

+|i)

if» =

y p g - * 0 ,9 8 5 ).

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

г ,д а 9 =

<P= 1 ).

 

Как видим, отличие приближенного выражения

(t, g)

от точного gi (t, |)

незначительно.

 

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

25

ivPfc1—Р

/HM+Pi

 

W

b [ n

1 1

т У

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

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

t

* / ( М =

I g ,{ t - t',t') ^ d t'

 

 

— оо

 

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

 

хI(К t) =

оо

 

$ gf(s, t s)eь «~s>ds.

 

 

о

 

Реакция системы по выходу i на входной сигнал

в виде

показательной функции, поданный на /-й вход,

 

xu(X,t) = wt!(X,t)eut

(8.1)

где

00

 

 

 

Щ1 0 = j gij (s, t s) e-^ds.

(8.2)

 

6

 

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

оо

Wf (X, t) = j* gf (S, t — s) e~Ksds

о

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

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

 

W(Я, t) =

(НУ]. (Я, t)

оуа (Я, t)

wt (Я, /)),

 

которая

связана с

матрицей

импульсныхпередаточных

функций

системы

следующими

эквивалентнымисоотноше

ниями:

 

со

 

 

 

 

 

 

 

 

 

W(X,t)= [о (s, t — s) e-^ds,

(8.3а)

 

 

О

 

 

 

 

 

t

 

 

 

 

W (X, t) = j G(t —t', f) e~>-«-‘‘Uf.

(8.36)

 

 

—oo

 

 

 

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

 

 

 

Wj(Kt) =

 

 

(8.4)

где Xj (Я, t) — решение дифференциального уравнения

 

A { t) - ^ ^ B ( t) x , + h,(t)^,

(8 .5 )

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

Столбцовая матрица Wj (Я, /) является решением диффе­ ренциального уравнения

VJ(0 - ХЕп) W + А (/) h, (t),

(8 .6)

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

- ^ = (У (/) - ХЕ„\ W + А~‘(<) Я (t),

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

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

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

t

 

x(t)= J Q(t — t',t')u(t')dt'.

(9.1)

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

“ w

= “2^ г 5 £(«)«“ *

( t > « j

 

(—ioо

 

(са — абсцисса

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

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

Лапласа).

Подставим выражение и (/) в (9.1) и поменяем порядок интегрирования. Получим

С—ioo

или

c+ioo г t

L (и) dk,

I — 00

x(f) = 2~ г I

I

L

i

L (u) eudX.

Cloo

 

— 00

 

Отсюда

 

 

c+loo

 

 

 

 

 

x ® = ~ S r I W(X,t)L(u)e>-‘dl.

 

 

 

c—ioo

 

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

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

 

 

 

L(x) = W(k,t)L(u).

(9 .2 )

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

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

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

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

(8.36),

в случае

стационарной системы

(А = const, Н =

const)

имеем

 

t

 

t

 

W(X) = J

 

= \

eua- i)e-u '- ,">dt'A-'H =

— GO

 

■— CO

 

 

 

OQ

 

= ]' eVse~uds A~'H = L(eUs) A~'H,

0

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

L (eUi) = Q,E— U)-'

(C M. § 4 гл. VII), TO

W (X) = {%E-U)-yA~xH.

(10.1)

Пусть J = diag (У2 (X^), Jz (Я2), ..., Jp (Xp)) — жорданова форма матрицы U, a К = (Kv Kt, •••, KP) — соответствую­ щая преобразующая матрица, так что

 

/

 

 

(М х

 

и = K J M

= / г 1 =

\М Р

 

 

\

 

 

 

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

 

W(Л) = 2

Ка

 

 

 

 

а=-1

 

 

 

 

 

причем

 

 

 

 

ka1

Ekn

 

Hk

 

 

*

 

 

Hk

кО

*п

+

+

Ъ

(XEka — Ja)_1 =

+

 

 

( X - X j

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

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

нестационарной системы.

 

10 .2 .1 . И с п о л ь з о в а н и е

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

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

ф у н к ц и и .

Учиты­

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

 

 

t

 

 

 

U?(X,t) = ]

— V) A~l(i')Hе K{t

ndt'

И

оо

Г (%, t) = [ G„fs, t — s) A~l (t — s)H ft— s) e~uds,

О

где

е0=хft) х г ' (п =х(ох-1ft— в),

а X (t) — фундаментальная матрица однородного уравне­ ния (2 .1 ).

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

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

 

 

10.2.2.

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

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

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

с и с т е м ы ч е ­

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

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

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

п а р а м е т р а х . Импульс­

ную переходную функцию

G (s,

t s)

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

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

t):

 

G (s,t — s) = G (s, t) — s - ^ - +

дЮ

dP

 

s3 dsG ,

31 dP +

Заменяя в (8.3a) G (s, t s) ее разложением в ряд Тей­ лора, имеем

сю оо

W(К, t) = j G(s, t) e~ksds -

j s —

e~Ksds +

 

о

0

 

 

 

 

d2G (s, t)

e k<>ds —

 

 

dP

 

Здесь

 

 

 

j G (s, t) e~k‘ds = R (Я, 0 = %o (Я, 0,

0

 

 

 

ac (s, /)

e k5ds =

d*R (l,t)

tfu (M )

dt

dh dt

и, вообще,

dlG (s, t)

e ksds d2lR (X, Q

Rli (K *)’

dt1

dtidt

 

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

 

W (X, <) = я (X, t) + Ru (*,,0 + 4 - Яи (Я, о +

(10.2)

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

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

рого аргумента имеет не t, а медленное время т =

et, то

тогда

 

 

 

 

 

00

 

 

W(А, т) =

| О(s, т — es) e~ksds

 

 

 

о

 

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

 

W(А, х) =

R (А, х) +

&RU (А, х) -f- -gp- «*#22 (k,т) +

 

10.2.3.

П е р е д а т о ч н а я ф у н к ц и я

к а к р е ­

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

= [UС*) — №п] w + А~1(х) hj (t)

(т = ei), (10.3)

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

U (х) уравнения

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

(10.3) разбиваются на непересекающиеся группы A(i0), А(20>, ...

=2. .... р; £ *< ,= я), так что она может быть

представлена в форме

и = S К„А„Мо.

0—1

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

также

разбиваются на

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

группы

вида

Х[а>-

X . W — X.......*,£> -

X (а = 1 ,2 ........ р;

=

в), а

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

U---кЕ„ = 2 Ко (Л(Г ---' кЕfe0) Мд. а—\

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

р~

Wj= HiKo(t,e)yoi,

а = ]

= [Ао (т, е) — %Ek<J\уа, + Ма(т, е) А~х(т) hj (t)

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

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

w{p = 2

0=1

к!тг) (Т, е)у%

лЛ')

-jjf- = [Л? (Т, 8) — XEk(J] уд} + м р (т, е) А~1(т) Л/ (IО

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

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

w{p = i ; к'р (t) у$,

du(r)-

= 1Л?’(t) -

I уЧ] + MX (t) A-' (t) hj (t)

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

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

т г - [ л р с о - ^ г

Тогда

t

уЩ= J Уа' (К t) Y^~' (X, t’) А# 1 (Г) А~' (Г) hj (Г) <#'

и соответственно w f (Я, t) =

= iltf'W . f Y V (К 0 Vf> V , t')M V [t')A - '(t')h ,(t')d t'.

-CO

В более компактной форме w f (Я, t) =

1

i

 

= K U) ( t )

j Y : n

( % , Г ) М ( Г ) ( П А - ' ( O h / V ' ) # ' .

— CO

Здесь

 

 

'MP'

Klr) = {K\r

K t

Kf), Mf,) =

 

 

Ж'<

Y,r) =

diag (Y\'\ V f .........Yf).

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

 

 

 

2

*!,'

(О Л

 

 

 

а=1

 

 

 

 

=

< * ? ( 0 -

ч

й ?

-

/И ? ’ (0 л - 1 (/) А , ( 0

 

 

(ст =

1 , 2 ......... п).

Тогда

 

 

 

 

 

 

ч>У (К I) =

 

 

 

 

 

 

;

л'

 

 

 

 

= 2 К » 1 (0

j

« Р J ( t f

(Г ) -

А) Л ’ М ? 1(/') <4_ l (<') A, (<') d i \

a=\

—oo

r

 

 

 

 

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

(Я, 0 =

//

= K{r)(t) j exp j (A(r)(i" )~ l £ n)d tM m

где

A ^ d i a g ^ i 'U ^ , . . . , Я<г)).

Г л а в а XI

ПРИБЛИЖЕННОЕ ИНТЕГРИРОВАНИЕ УРАВНЕНИЙ УПРАВЛЯЕМОГО ПРОЦЕССА

§ 1. Интегро-дифференциальная система уравнений управляемого процесса

Будем рассматривать управляемый процесс, течение ко­ торого представляется некоторыми параметрами (коорди­ натами) xv х2, ..., хп, удовлетворяющими системе уравне­ ний

П П I

2

ац (*) -^r = £

 

Ьи ® х>+ £ huui

0 • О

/=1

/=1

 

/—I

 

 

(i = 1 , 2 , . . . ,

п\

det (ац) ф 0).

 

Управляющие воздействия uj, рассматриваемые как вы­ ходные сигналы регулятора, предполагаются линейными

функциями входных сигналов регулятора

о/, которые фор­

мируются как линейные комбинации координат xlt x2t

хп:

v}-= l 1

tuxk

(j = 1 ,2, . . . ,

m).

(1.2)

A=1

 

 

 

 

Допустим, что

связь

между входными сигналами

vv v2...... vm и выходными сигналами wlt и2,

ы* регулято­

ра представлена посредством импульсных переходных функ­

ций ga (t — Г, t') (i = 1, 2,

l\ j = 1, 2 ,

m), так что

« i

 

(i.3)

\

 

/—1 — C O

 

 

Итак, рассматриваемый здесь процесс полностью опи­ сывается системой уравнений (1.1), (1.2), (1.3). Запишем эту

Соседние файлы в папке книги