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

книги из ГПНТБ / Цейтлин Я.М. Проектирование оптимальных линейных систем

.pdf
Скачиваний:
2
Добавлен:
24.10.2023
Размер:
9.35 Mб
Скачать

сходящегося, по крайней мере, при малых значениях s (что со­ ответствует большим значениям t), то из сравнения двух выраже­ ний для Y (s) следует, что

Си = С* •

Коэффициенты С* находятся путем деления полинома М (s), являющегося числителем передаточной функции, на полином N (s), являющийся ее знаменателем. Деление осуществляется по обыч­ ному правилу деления полинома на полином.

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

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

деляется / +

1 моментами функции

веса

( f i 0 , ц^, . . .,

так

как члены ряда (1.69), начиная с / +

1-го,

становятся (точно или

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

равными нулю.

 

 

 

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

зована

в функцию

 

 

уй (t) = Н(р)х (t),

 

где Н

(р) — заданный линейный стационарный оператор,

назы­

ваемый

оператором идеального преобразования.

 

Обозначим передаточную функцию синтезируемой системы

через F (s) и соответствующую ей функцию веса через Р (і).

Функ­

цию веса, соответствующую оператору идеального преобразова­

ния, обозначим через Р0

(t).

 

Функция на выходе синтезируемой системы в установившемся

режиме согласно (1.69) может быть представлена в виде

</(0 = £с,*<*>(0,

 

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

 

M 0 = 2 № ( f e ) ( 0 .

 

 

k

 

Условием равенства

выражений у (t) и у0 (t)

является равен­

ство коэффициентов Ck

и СІ, т. е. равенство

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

моментов синтезируемой

функции веса и функции веса, соответ-

ствующей оператору идеального преобразования. Такие условия имеют вид

со

lxkP(x)dx=[i0k,

О

со

 

 

 

 

 

 

 

 

 

 

 

 

 

где fife = J x(k)Po

(т) dx — заданные

константы,

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

10. ЧАСТОТНЫЕ

 

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

 

 

 

Предположим, что на вход нашей стационарной системы посту­

пает гармоническое

колебание

 

 

 

 

 

 

 

 

 

 

 

х (/) = A cos со/.

 

 

 

 

Уравнение (1.62)

прлмет

вид

 

 

 

 

 

 

 

 

N

(р) у (t)

=

М

(р) A cos

со/.

 

 

(1.70)

Так как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cosarf

= -i-(e'a <

 

+е~ш),

 

 

 

уравнение (1.70) эквивалентно

следующим двум

уравнениям:

 

N(p)y2

 

 

 

У)=*М(р)±<Гш,

 

 

 

где

 

 

Ух (0 + у* (0 = у

(t).

 

 

 

 

 

 

 

 

 

Будем искать

частные

решения

ух

(t)

и у2

(t)

этих уравнений

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

ух

it) =

Д е ' <

M

i

> =

Л і Є ' * ' еш;

 

1

(1.71)

у2

(/) = Д,е~'(

ш Н

* 2

)

=

Л 2 е - ' * ' е - ' ш ' . J

 

Подставив (1.71)

в (1.70)

и учтя,

что

 

 

 

 

 

р е

=

 

; =.(гсо) е

,

 

 

получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v * w

( и е ' " ' = 4 м

( и е ' ю / ;

 

 

 

 

( - /со)е - і а /

 

 

 

/со) e - ' w ' (

 

откуда

Л а Є

_

2 i V ( - t o ) ~ 2 Ґ

<

т > -

Функции

 

 

 

 

 

 

F (/со) =

F (s) | S = 1 M

и

F (—/со) =

F (— s) | s = i u )

можно представить

в виде

 

 

 

 

F (/со) = Re [F (/со) ] +

Im

[F (/со) ] = \F (/со) | е1'ф

где

 

 

 

 

 

 

| F (/со) | =

V Re2 [F(гсо] + I m 2

[F(гсо)];

Очевидно, что

F (— /со) = Re [F (— /со) ] + Im [F (— /со) ] = = Re [F (ico) ] — Im [F (/со) ] = \F (/со) | е~('ф

Подставив (1.72) и (1.73) в (1.71), получим:

(1.72)

(1.73)

 

Лі Є '*' =

^-|F(/co)|e'>('w );

 

Ля е-'*« =

- 4 - 1 F (ico) 1 е~'ф ( *№ > ,

откуда

 

 

 

 

Фі = Ч>2 = ф (/©)•

(1.74)

 

)

С учетом

(1.74) выражения (1.71) примут вид:

 

У і ( 0 = 4 і ^ ( И І е ' [ в ' + , , ( ' в ) ] ;

 

у. (О =

! z7 (to) | е-^

 

Поэтому

 

 

/

У (0 =

У! (О + Уг (0

= Л | F (/со) | cos

[со/ + Ф -(/со)].

Таким образом, амплитуда гармонического колебания на вы­ ходе линейной стационарной системы равна произведению ампли­ туды обусловившего его гармонического колебания на входе на модуль того комплексного выражения, в которое обращается пере-

52

даточная функция системы F (s) в результате подстановки в нее

s = /СО.

Фаза гармонического колебания на выходе отличается от фазы соответствующего гармонического колебания на входе на величину аргумента того комплексного выражения, в которое обращается передаточная функция системы F(s) в результате подстановки

внее s =^т

Всвязи с изложенным функция

F (/со) = F (s) | s = 1 ( 0 = Re [F (/со)] + Im [F (ico)]

называется амплитудно-фазовой частотной характеристикой си­ стемы.

Модуль этой

функции

 

 

 

 

| F (їй) | =

F (s) F ( -

s) \ 1

= ш = 1/Re2 [F (/со)] + Im 2

[F (/со)]

называется амплитудно-частотной характеристикой системы.

Аргумент

этой функции

 

 

 

 

 

 

cp(/co) = arctg ^ g g g j

 

 

называется фазовой частотной характеристикой системы.

Функции

Re

[F

и Im [F(/co)], равные

вещественной и

мнимой

частям

амплитудно-фазовой частотной характеристики,

называются вещественной

и

мнимой амплитудно-частотными ха­

рактеристиками

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

 

 

Из

того,

что знаменатель

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

совпадает

с оператором

левой части дифференциального

уравнения (если

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

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

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

Иными словами, передаточная функция устойчивой системы не должна иметь ни одного полюса в правой полуплоскости, т. е. должна быть аналитической функцией в правой полуплоскости (включая мнимую ось). Существуют алгебраические (Гурвица, Раусса) и частотные (Найквиста, Михайлова и др.) критерии

устойчивости. Сущность упомянутых критериев достаточно хо­ рошо изложена в обширной литературе'тю теории автоматического регулирования и управления.

Не будем касаться здесь этих вопросов, так как при синтезе оптимальных систем требование устойчивости учитывается авто­ матически: оптимальная система (т. е. система, удовлетворяющая заданному критерию оптимальности) ищется в классе устойчивых систем — таких систем, передаточные функции которых не имеют полюсов в правой полуплоскости.

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

В соответствии с выражением (1.68) функция на выходе выра­ жается через производные функции х (t) на входе и моменты функ­ ции веса. Для того чтобы функция у (t) была ограниченной, тре­ буется ограниченность моментов функции веса.

Таким образом, требование

устойчивости

системы совпадает

с требованием

 

 

 

со

 

 

 

_ [т*Я(т)гіт <со,

k =

0, 1, 2,

. . .,

о

 

 

 

если

 

 

 

dv

0,

v>N.

 

x(t) =

 

dtv

Глава II

СЛУЧАЙНЫЕ ФУНКЦИИ В ЛИНЕЙНЫХ СИСТЕМАХ

И . С Л У Ч А Й Н Ы Е Ф У Н К Ц И И

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

Вреальных системах автоматического управления сигналы на их входе и выходе весьма часто являются не детерминированными,

аслучайными функциями, т. е. функциями, конкретный вид ко­ торых меняется от опыта к опыту.

Иногда случайные изменения характера протекания процессов

вотдельных конкретных опытах невелики и ими можно прене­ бречь. В этих случаях допустима идеализация, сводящаяся к тому, что соответствующие процессы считаются детерминированными функциями. Однако во многих очень важных ситуациях указанная идеализация оказывается недопустимой и соответствующие про­ цессы следует считать случайными функциями. Система автомати­ ческого управления осуществляет, таким образом, преобразование случайных функций, в связи с чем возникает проблема установле­ ния закономерностей такого преобразования.

Прежде чем перейти к установлению этих закономерностей, необходимо дать определение случайной функции и ее основных характеристик [10].

Под случайной функцией х (t) понимается функция, значение которой при каждом фиксированном значении аргумента t яв­ ляется случайной величиной, а при каждом фиксированном опы­ те — некоторой неслучайной функцией (данной реализацией слу­ чайной функции).

Это определение иллюстрируется рис. 9, где изображено не­ сколько реализаций случайной функции х (t). Реализации обо­

значены

) , 2, 3, 4, . . . Действительно,

при фиксированном

t =

tx

случайная функция х (t) является случайной

величиной

х

(^),

принимающей значения

xtl, х12,

. .

.,

xlk, . .

. При фиксирован­

ном t =

t2

случайная функция

х (t)

является

случайной величи­

ной х (t2),

принимающей значения х21,

х22, . .

.,

х2ь . . .

х

При

фиксированном

опыте

1

случайная

функция

является детерминированной функцией— кривой 1, при фикси­ рованном опыте № 2 — детерминированной функцией — кривой 2 и т. д.

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

Под одномерной функцией распределения случайной функции х (t) понимается функция

Ft (х, t) = Р (0 < * ] ,

где Р [. . . ] — вероятность.

Рис. 9

Одномерная плотность

вероятности

 

 

 

 

 

 

 

 

 

 

 

 

fl(*. ')=•

dFx(x,

t) .

 

 

 

 

 

 

 

 

 

 

 

дх

 

 

 

 

 

 

двумерная

функция

распределения

 

 

 

 

 

 

 

F2

и

хг,

tlt

t3)

=

Р

(tj

< л ; ь

х

{t2)

 

<.х%];

 

двумерная

плотность

вероятности

 

 

 

 

 

 

 

 

 

 

 

t

(v

v

і

і

\

 

2 (Х\,

х2 ,

і\,

t2)

,

 

 

л-мерная функция распределения

 

 

 

 

 

 

 

 

F/i (xl> •

•' xn>

11>

•>

tn)

=

P

[x ( ^ l )

<3

#1»

і і

'і X (І„)

< < Л Я ] ' ,

n-мерная

плотность

вероятности

 

 

 

 

 

 

 

 

 

f i

Y

 

Y

t

 

'

і

\

dnFn(xi

 

 

Xn,

^i

in)

 

UK*!,

• • .»

xn,

н> • • •»

Ы) —

 

 

dxxdxr

 

. .

dxn

 

и т.

д.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

55

Если

рассматриЁать

не одну случайную

функцию х

(t), а

две — х

(t) и у (t), то

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

совместные

законы

распределения ординат этих двух случайных функций в различные

моменты

времени,

например:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F*(x,

У,

Н

t')

=

P[x(t)<x,

 

 

y(t')<y];

 

 

 

 

 

F4(xu

 

x2,

yu

y2,

 

U,

t2j> U,

U) =

P[x{ti)<ixh

 

 

 

 

x{t2)<x2,

 

 

 

 

 

 

 

 

 

 

УІи)<Уі,

 

 

 

y{h)<y^;

 

 

 

 

 

 

 

f

(x

у

 

-tj

1,

f

 

t

і

 

i'\

=

 

d V 4 ( * i ,

x2,

yb

y2,

 

tb

t2,

t\, t'2) ш

ГЛ*и x2,

yu

У г ,

tu

 

h,

h,

t2)

 

 

 

 

дХідХгдУіду2

 

 

 

'

 

 

F2n(xi,

 

 

. . .,

xn,

 

 

yi,

. . .,

yn,

tu

. . .,

tn,

 

t\, . .

., tn)

=

 

=

P{x{U)<xu

 

. . ., x{tn)<xn,

 

 

y(t[)<yu

 

 

. . . ,

 

 

y(Q<yn];

 

 

f2n(XU

• . .,

Xn,

 

 

У\,

.

. .,

Уп,

tl,

...

.,

tn,

h, . .

.,

tn)

=

 

 

 

_ d 2 " ( * i

 

 

 

xn,

yb .

. .. yn,

tu

.

.,

tn't'x

 

 

 

Q

 

а также

 

 

 

 

 

 

 

dxt.

 

. . дхпдУі.

. . дуп

 

 

 

 

 

'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fkt,kAxU

 

• • •>

 

 

 

Уи

•»

Ука,

h,

. . .,

tkl,

tx,

. . .,

tk2)—

=

P[x(tx)<xx,

 

. . .,

x(tkJ<xkl,

 

 

y(t'x)<yx,

 

...

.,

 

 

y(t'kt)<yk,];

 

 

/fe,ft

(

x,

. .

.,

Xk

 

 

Ух, . . .,

Укг,

U, . . .,

tkit

ti,

. . .,

thl)

=

 

 

2

x

 

 

 

 

lt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;_

 

d k l + k 2 F k l k 2 ( x i

 

 

 

 

 

Уг

 

 

 

 

'і- • •• V

 

h • •

*kt)

 

 

 

 

 

 

 

 

 

 

dxv

 

. . дХкхдУі

 

. . .

ду^

 

 

 

 

 

 

для любых значений

 

kx

и k2

и т. д.

 

 

 

 

 

 

 

 

 

 

 

В литературе встречаются понятия «случайная функция»,

«случайный

процесс»

 

и

«случайная последовательность».

Первые

два понятия тождественны, хотя иногда случайными процессами называют случайные функ-ции времени t.

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

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

оо

М [х (t)\ = тх (t) = \ xh (х, t) dx.

(II. 1)

—со

 

Математическое, ожидание-случайной функции представляет собой неслучайную функцию, которая при каждом значении аргу-

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

На практике, когда плотность вероятности f1 (х, t) неизвестна, использование формулы ( I I . 1) не представляется возможным. В этом случае оценка математического ожидания может быть опре­ делена в результате статистической обработки полученных реали­

заций X[ (t) случайной функции х (t).

В качестве такой оценки ис­

пользуется

 

 

 

 

mx(t)

=

-jrtxl(t).

 

(П.2)

Центрированная случайная

функция

х° (t) определяется как

х° (t)

= х

(0 -

тх

(/).

Начальные и центральные моменты порядка к, характери­ зующие степень разброса случайной функции х (t) в любой мо­ мент времени t, являются неслучайными функциями и опреде­ ляются по формулам:

со

ak(t)=

J х * М * .

t)dx =

M[^(t)];

со

—со

 

 

 

 

 

М 0 = I

U-mx(t)]kfi(x,

t)dx =

M[[x\t)]k}:

—со

 

 

 

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

со

a i

(t) = М[х

(t)] =

тх (t) =

J хЬ (х,

t) dx

 

 

 

 

. — с о

 

 

 

и ее дисперсия

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

^2 (*) = М {[хй

(О]2 \=Dx(t)

=

D [х (01 =

J

[х -

тх

(*)] h (х, t) dx.

 

 

 

 

 

—со

 

 

 

Очевидно,

что

 

 

 

 

 

 

 

Действительно,

 

V-i (t) = о.

 

 

 

 

 

 

 

 

 

 

 

 

ii1(t)

= M[x°(t)]=M[x(t)~mx(t)}

 

=

 

=

М[х

(О] -

М [тх

(t)] = тх

(0 -

тх

(t) =

0.

Оценкой

для дисперсии

случайной

функции

х (t) является

 

"

1 ( = i

 

 

 

где тх (t) определяется по

формуле (II.2).

 

случайной

Степень статистической связи между ординатами

функции х

(t) в моменты времени tu t2,

. . ., tk,

• • .

характери­

зуется неслучайными функциями — начальными и центральными моментами данной случайной функции. Начальный момент по­

рядка

k

определяется

как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

«* (tu

h,

 

tk)=M

х (tt)

. . . х {tk)\

=

 

со

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= J

• • •

І

ЧЧ

• • • 4h

 

(Ч,

Ч,

•• • Ч, к,

к,

• • ., ik)d4

d 4

• • • d 4 .

— 0 0

 

—-со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Центральный

момент

порядка

k

определяется

как

 

 

 

 

Н

(h, U, ...,h)

=

M [х° (h) х° (/,) ...хй

(tk)]

=

 

=

М\[х(^)-тх

 

(h)] [х (t2) -

тх

(Q1 ...[х

(tk) -

тх {tk)\}

=

 

 

со

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

і

• • •

і

І*і ~

т

* ('i)l [ Х 2 — тх (t2)] ...\xk

— mx

(tk)]

х

 

 

— со

 

— со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X /* ( х ъ

. . . , xk, tu

 

. . . tk) dxl...

 

dxk.

 

 

При

k =

1 получим

рассмотренные

ранее

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

«!

(t) = М[х

(t)] =

j х\х

(х,

t) dx = mx

(t);

 

 

 

 

 

 

 

 

 

 

 

— CO

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Их (t) =

M

[x° (01 =

0.

 

 

 

 

 

При

k =

2

будем

иметь:

со

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

«2

 

 

*я) =

Л1 fo) х (4)] =

j

|

Вдг/а

( х ь

х2 ,

tu

t2) dxx

dx2;

 

 

 

 

 

 

 

 

 

— с о

C O

 

 

 

 

 

 

 

 

 

\it

(tlt

t2)

= M\

\x ft) — mx

(У ] [x (ta) mx {t2)\}

=

 

 

 

CO

CO

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

J

J

l4

— mx(t1)][x2

 

 

 

 

mx(t2)]f2{x1,x2,t1,t2)dx1dx2.

 

— C O

— C O

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Второй центральный момент ц2 (tlt t2) носит название корре­ ляционной (автокорреляционной) функции и для каждых двух фиксированных моментов времени tx и t2 равен корреляционному моменту (моменту связи) случайных величин х (tx) и x{t2). Итак, корреляционная функция

 

Кхх

*,) = Ц2

(*ь t2) = М [х°

х0 (4)]

=

 

=

М {[х (tt) -

тх ft)] ft) -

тх ft)] \

=

со

со

 

 

 

 

= J

J [x1 — tnx(ti)][x2

mx(t2)]f2(x1,x2,t1,t2)dx1dx2.

—W— со

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