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

книги из ГПНТБ / Основы авиационной автоматики учеб. пособие

.pdf
Скачиваний:
50
Добавлен:
22.10.2023
Размер:
13.63 Mб
Скачать

и понятие моментов J,, * этой функции, определяемых выраже­ нием

“ (/) да

(5.61)

Ji,k = \ z ( t ) z [ t ) d t .

о

Можно показать, что изображение по Лапласу этой функции

имеет вид

 

______________1

 

 

 

 

 

 

1

 

L[z{t)\

 

 

 

 

 

 

 

 

 

( Р) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

КР" + К-\Р"~І + ■• •+ ^ l P + K

 

t

\ p

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=0

 

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

применяя

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

Лапласа

к

(5.62)

уравне­

нию

(5.59) и учитывая

его

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

(5.60),

находим

 

 

L [Х„ 2

(t)

 

Хп_, 2 (t) + . . . +

Хі 2 (f) +

Х0

Z(t)] =

 

 

= K \ P n z { p )

- 2 ( 0 ) ]

+

К

-

І р л ~ 1 г { р )

+

. . . + \

p z { p )

+

Х0г(/>)=

 

 

Рп г(Р) — \

1

+

Хл_, рп~' z{ p) + . .. + \ p z (р) +

Х0z(p)=

 

 

= 2 (р) \ъврп +

Хл_1/7"-1+

. • . +

Xj р +

X0J — 1

=

о,

 

откуда получаем для

z(p)

 

выражение

вида (5.62).

Отметим

также, не приводя доказательств, ряд свойств моментов

У,-,*

функции z(t):

моментов

 

с симметричными

индексами

 

а)

равенство

 

 

 

 

 

 

 

 

 

Ji,k — Jk,h

 

 

 

 

 

 

(5.63)

б)

положительность моментов

с одинаковыми индексами

 

 

 

 

 

 

 

 

 

0;

 

 

 

 

 

 

 

(5.64)

в)

свойство моментов с четной суммой индексов

 

 

 

 

 

 

Уі, к == У(й+2і)), к == (

1 )1' У(*+т|), (Я-т)) >

 

 

(5.65)

г)

 

свойство

моментов

с

нечетной

суммой индексов,

если

і + k < 2га— 1:

 

 

 

 

 

Уг, А — о,

 

 

 

 

 

 

(5.66)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

і + k < 2п — 1 — нечетное число;

д)

 

свойство момента У,-,*

= У(П-і),л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У(л -1),„=

 

 

 

 

 

 

 

(5.67)

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

Л / =

Ап, (/+п

(5.68)

 

2ХЛД

 

268

где Xn— коэффициент при р п знаменателя г(р), или коѳффи-

(«)

циент при z(t) уравнения (5.59);

Д— детерминант п-го порядка, составленный из коэффи-

циентО'В h

0

0 - Х 0

Д =

6 0

0 0

К

- 0

 

- х 3 . .

. 0

 

х ,

 

 

:

. \

 

 

\

 

 

 

\

о

 

1

1 С

 

 

 

 

 

го

0

0

0

(5.69).

- X

К

Дп, 1+1 — алгебраическое

дополнение Д

для

элемента

п-ой

строни и +

1)-го столбца.

 

следующие

преобра­

Для доказательства

(5.68)

 

проделаем

зования:

уравнение

(5.59)

последовательно

на z(t),

1) умножим

(л-1)

результате

чего

получаем

систему

уравнений:

z (/),..., z(l),в

(л)

(л-1)

 

 

 

 

 

 

 

 

 

 

 

К z (t)z ( 0 +

K - i z ( 0

z ( 0

+

• ■ +

^ i z ( 0

z ( 0

+

x 0 z 2

( 0

=

(л) .

(л-1) .

 

 

 

 

+ \

 

 

 

 

 

К z (t)z (0 ~h K-i z (0z(t)~h ■••+ ^i z < 2 (0

z

{

t ) z { t )

(л) (л—1)

(/.-1)

(i)

+ •

■-f- \

. (л-1)

 

 

(л-1)

 

\ nz (t) z ( t )

Xn_j z2

z(t) z(é)-\-~k0z (t)z(t) — 0;

 

 

 

 

 

 

 

 

 

 

 

 

(5.70)'

2) проинтегрируем по времени от 0 до сю систему (5.70), в ре­ зультате чего, с учетом (5.61), получаем систему уравнений от­ носительно У;, к-

Х„ Ул,о + ^„_і У(л—1), о + . • • + /ц Уі.о + Х3 А о = 0;

Х„Ул,і + ^,і-і-Ал- d. 1+ • •

Уі, 1+ Х0Л>, 1“ 0; (5.71)

Х„ /л.(л—1) + \ - 1 У(л-1), (л-1)

+• ■• +^іУі,(л-1)-|-0,(л-1)= 0;

3)преобразуем систему (5.71), используя свойства момен­

тов У,-, * (5.63) -+- (5.67), в результате чего получаем систему

269

уравнений относительно У/,

 

 

^0 У), О— X, -А, 1Н~ ^4 J%2 —Xß Уз, 3

+ Xg -У. 4 —

.......................... =

ОУо,о + Xt Уl, 1 — Х3 Уг, 2 + ^5 Уз,з

^7 Уі, 4 + ........................

= 0;

О Уо, 0 \ К х + X, У*. 2 -- ^4 Уз, 3

+ У),4 — .............................

= 0 ;

I. ..........................................................................

^Я-1 У(Л-1),(П-1)--- —— •

 

^Ал

 

 

(5.72)

Система уравнений (5.72) представляет линейную систему из п

уравнений с и неизвестными

У/, решение

которой определяет­

ся формулами Крамера:

 

 

 

 

 

 

 

У .,=

\

 

 

 

(5.73)

 

 

 

Д

 

 

 

 

где Д — детерминант системы (5.72):

 

 

 

 

^0

^2

 

- > - С

Х 8

. .

0

 

 

 

 

 

0

 

^3

\

- х 7 .

. .

0

 

 

 

 

 

 

0

~~ ^0

^'2

- к

V

 

.

0

А

 

 

 

 

 

 

(5.74)

\

ОО

Л“ 1

Д; — детерминант, полученный из А путем замены + 1) -го столбца А коэффициентами правой части системы

О

(5.72), т. е. столбцом

О

■ Поэтому

 

 

 

 

 

 

1/2ХП

 

 

 

А

2Х„

 

 

(5.75)

 

 

 

 

■где Д„,(,+і) — алгебраическое дополнение Д для

элемента

п-ой

строки и + 1)-го столбца.

 

 

У, ,• выраже­

С учетом (5.74), (5.75) из (5.73) находим для

ние в виде (5.68).

 

квадратичной

оценки

J —

Для вычисления интегральной

со

= J i 2(É)dtустановим ее связь с моментами У/,/ функции z(t).

о

270

Для этого 'Прежде заметим, что если изображение функции ^(і)

имеет вид

 

 

£

\ р 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЧІР)

і-0____

 

x,„ Р т +

- 1 Рт 1 +

 

■ +

*1 Р +

хо, т<.п,

 

 

£

\ р ‘

 

*„/>" +

\j- l і0,1-1

+

 

• +

Р +

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

і- о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то,

учитывая вид изображения z(£)

(5.62),

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

 

 

 

 

 

 

т

 

 

 

 

 

т

 

 

 

 

 

 

 

 

т (у) =

S

^

"н--------= 2

 

 

 

Z7' * (/>)>

 

 

 

 

 

 

 

і—О

£

^ У*

/=0

 

 

 

 

 

 

 

 

 

 

 

 

 

і=0

 

 

 

 

 

 

 

 

 

 

откуда находим связь -[(7)

с z(t):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l ( t ) = L

 

£

-'чР12{Р)

 

 

 

 

 

 

 

 

 

 

 

 

 

і=О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

(m—1)

 

(т)

 

=

£

*і 2 (0 = *о Z (0

+ У-! 2 ( 0 +

• • • +

хт-1

2 (0 + xmz (0- (5-76)

 

/-0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выражение (5.76) позволяет с учетом (5.61),

(5.63) — (5.67)

за­

писать:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ =

[ т 2{ t ) d t =

\[y.Qz[t) + *,z(0 + .-H- xm- U ( b + xmz(OJ2^

=

 

 

6

 

 

О

 

 

 

 

 

 

 

 

 

 

 

 

 

=

~

[V z2(O + 2^

 

.

. . .

 

 

 

 

 

 

 

im)

 

j

1z(0 z ( / ) + 2x0x2z(i)

2(O+... + 2*0*eiz(/) z(*) +

 

o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

.

 

 

. ..

 

 

 

 

 

 

 

. ( m )

 

 

 

 

'^і22Г" (t) +

2Xj x2z(/) Z (£)+ - . ■+

2xtv.mz (t) Z (0 +

 

 

 

 

 

 

 

+

x22 2:2 (0 +

■-

 

 

 

(m)

 

 

 

 

 

 

 

 

• + 2*2 y-mz {t)z {£)+••■

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[m)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ *т*- (*)] dt

 

 

 

V

-7o,о +

2-/-0

/о, 1 +

2x0 x2 У0, 2 -f-. •

+

2'/.0 xm

y0, /л +

 

 

 

 

 

y-\ J\, 1

2Xj x-2 J1, 2 +

 

• +

2 x, y.m У!, m -f-

 

 

 

 

 

 

 

 

+

x22

Jt, 2

+

 

• . . +

2x2 xm У2,m + . •

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ X 2

jm, m

 

271

 

— х02 ■Л), 0

2Xq*2 J\, 1-{- 2*0 X^ •/2, 2

............................

 

 

4“

 

 

+

xi2

J\, 1 — 2xjXj / j,! +

.......................................

 

 

+

 

 

 

 

 

 

 

'*■2 'J2, 2 — ............................

 

 

't-

 

 

 

 

 

 

 

 

 

 

 

' Ь Xm J 71*ui

 

где

 

BqJo, 0 +

J\, l

+^2^2, 2 +

• •

+

Ba Jm, m,

 

(5.77)

 

 

 

ßn

*

2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*0 >

 

 

 

 

 

 

 

 

 

 

5, = X,2 2x0 x2;

 

 

 

 

 

 

 

 

 

B 2= X22 - 2xj X3 +

2x0 X4;

 

 

 

 

 

 

 

Д =

*2.

 

 

 

 

 

 

 

 

 

 

 

/72

 

 

 

 

 

 

 

Учитывая полученное ранее выражение (5.68)

для

момен­

тов J;,i,

из (5.77)

находим искомое выражение для J:

 

 

 

Алл

+В,

 

Дл,3

,

,

о

Ал, /л+1

 

 

J = B ,о —

 

+ ß 2 u"-3

 

 

 

 

 

 

 

2Х..Д

2Х„Д

5 Г І + - " + в ” ! Г Х

 

[ ß 0 V l +

Bt^n.2 +

б

2 й я ,3 +

• •

• +

BmД п , m + 1 ] = -

° .

2Х„ Д

 

 

 

 

 

 

 

 

 

2Х„ Д

где

Aß — детерминант, полученный

из

Д путем

замены его

 

 

последней строки коэффициентами

Bt:

 

 

*0

0

0

ßo

1 to

 

О

0

X, •

-

• •

0

- \

\

.

\

\*

0 . . . х„_2 - К

ßl . . . Вп-2 Ва- 1

3.

Примеры использования интегральных

 

квадратичных оценок

П р и м е р

1.

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

с передаточной функцией Ф(р)= -— ——, если ? ( / ) = / / ( / ) —

( Р+ ' Г

- Я ( о о ) .І ^ ) .

2 7 2

Ре ш е н и е .

1.Определение ■\(р) = L [H(t) — Н ( оо)<1 (а?)]:

а) определение L[H{t)]:

 

L[H{t))=H{p):

Ф (

р

) _

1

.

 

р

 

 

 

 

 

 

 

 

 

 

р (

р + 1

) 2

б)

определение L [ Н(оо)-1(£)]:

 

 

 

 

 

L[H(cx>).\{t)} = L {< P ](p )p^

-1(0]

 

Р

 

 

 

 

 

 

 

 

 

 

 

в)

определение

f

(р):

 

 

 

 

 

 

 

 

тг(/0 =

 

1

 

 

_L= 1 - р 2- 2р - I

 

 

 

 

 

р

 

р ( р + I)2

 

р

( р

+

! )

2

 

 

= — р — 2 __

X,р + х0

 

 

где

 

(Р ■+•

I)2

 

^22 Р2 + \ Р + ^0

 

 

 

 

 

 

 

 

 

 

 

 

тп=

1,

 

V.Q—

2,

Xq=

1,

 

п = 2;

 

Xj =

— 1;

X, =

2,

 

 

 

 

 

 

 

 

 

Х2 — Хп — 1.

2. Определение J =

1 -[2{ t ) d t ~ — —

 

 

 

 

 

 

J

 

 

 

 

2Х2Л

 

 

 

 

 

 

о

 

 

 

 

 

 

 

а)

определение

Д:

 

 

 

 

 

 

 

 

 

Д —

Xq

— Х2 =

Х0 X] =

2;.

 

 

 

 

 

 

0

 

X,

 

 

 

 

б)

определение

Дя :

 

 

 

 

 

 

 

 

 

 

 

 

 

Bq=

у-о2 =

4)

 

 

 

 

 

 

 

ß , =

v =

i,

 

 

 

А* =

—х2

=

X0ßj +

X2ß 0 =

1' +

4 = 5;

 

А,

 

А

 

 

 

 

 

 

 

в)

определение J:

 

 

 

 

 

 

 

 

 

 

/ =

 

2Х,Л

 

 

 

= 1,25-

 

 

 

 

 

 

2-1-2

 

 

18. И зд. № 5312

 

 

 

 

 

 

 

 

 

27$

П р и м е р 2.

Определить интегральную квадратичную оценку / для си­ стемы, структурная схема которой показана на рис. 5.14, если

4(0 = g(0-

Р и с. 5.14. Структурная схема линейной стационар­ ной системы

Ре ш е н и е .

1.Определение у (/?) = /.[^(01 = Ф(/7):

5(/ >+ 1)

/7(2/7+ 1)(0,5/7+1)

ТГ(р) = Ф(р) =

j,______5 + 1)

Р(2/?—|—1) (0,5р -1-1)

5 1/7 +

1)

 

 

 

 

х ,/7 +

х„

Р (2Р + 1) (0,5/7 +

1)

+

Ъ(р +

1 )

1'з/73+ ^2/?2 +

КР + К

где

 

 

х0 = 5,

 

 

 

т = 1,

 

Л0 = О,

 

л =

3;

 

хі =

5;

=

6,

 

 

 

 

 

 

Х2 =

2,5,

 

 

 

 

 

 

 

X =

1.

2.

Определение

а)

определение Д:

д =

 

0

0

К - К

 

0

 

б)

определение Дв:

/ = I 72 (t) dt = —

J

2Xg Д

о

 

— XQ(Xj X2

XqXj) — 5 (6 • 2,5 — 5 • 1) = 50;

 

 

ß 0 =

x02= 25,

 

X2

5, =

x,* = 25,

 

0

~ 5' 1•25+2,5* 1*25 = 187.5

0

 

-X , =

Bo

Bl

0

 

274

в) определение J:

J =

s

1 8

7 ,5

 

 

 

= 1 ,8 7 5 .

 

П р и м е р 3.

2Х3 Д

2 - 1

- 5 0

 

 

 

 

 

 

 

1

Для системы с передаточной функцией

Ф{р)

/?2 + 2 . ^ + 1

выбрать значение

5 =

50

 

 

 

из условия минимума интегральной

квадратичной оценки /,

если к (£) = Н (t)

7/(оэ)- 1(f),

Ре ш е н и е .

1.Определение -\(р) = L \Н {£) Н (со) -1 (/)]:

а)

определение L[H{t)\.

 

 

 

 

 

 

 

L [H(t)\ —Н (р) =

Ф(Р)

 

1

 

 

 

 

 

 

р

 

/?(/72+ 2|/7+1)

 

б)

определение L [Н(со)-1(£)]:

 

 

 

1

 

L [Я(оо) - 1 (0] -

L [Ф (/7)р.о • 1 (0] = L

[1 • 1 (0]

 

= — ;

в)

определение т(р):

 

 

 

 

 

 

р

 

 

 

 

 

 

 

 

1

 

 

1

 

- 0 - 2 1

*1Р + у'0

7 (Р) =

 

1)

/7

 

/72+

2 5 /? + 1

Х*/72+ Х ,Р + Х с

где

/?(/72+ 2 $ /? +

 

т = 1,

 

*о = - 2 3 ;

 

 

 

 

 

Xq— 1>

 

 

Л = 2;

 

*і = — 1;

 

Xj =

25;

 

 

 

 

 

 

 

 

 

Х9 — Х„ = 1.

 

2.

Определение

J — \~^{t)dt==

2Х , Д

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

определение

Д:

 

 

 

 

 

 

 

 

 

 

Д =

 

Ч)

г'2

=

Х0 Х1

= 2 |;

 

 

 

 

 

0

X,

 

 

 

 

 

 

 

 

 

 

 

б)

определение

 

 

Д0 =

ѵ

=

^ ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ß i - V =

Ь

 

 

 

 

Xq

Х2 Х0Д, +

Х2Д0 — 1+ 452;

 

 

д 0

в х

 

 

 

 

 

 

в)

определение J:

 

1 +

4?2

1 +452

 

 

J =

 

 

 

 

 

 

45

 

 

 

 

2 Х ,Д

2 - 1 - 2 5

 

 

 

18*

275

3. Определение Е— S0 из условия У(Е0)—шіпУ:

dJ(i)

(1 -f 4Е2)-4 — 4?-85

dl «-«.

Ш 2

1 +

4So2- 8 S 02 = O - So = 0,5.

§ 5.4. ТОЧНОСТЬ ЛИНЕЙНЫХ СТАЦИОНАРНЫХ CAP, ПОДВЕРЖЕННЫХ ДЕЙСТВИЮ СЛУЧАЙНЫХ УПРАВЛЯЮЩИХ

ИВОЗМУЩАЮЩИХ СИГНАЛОВ

1.Постановка задачи. Статистические характеристики

случайных сигналов

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

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

В данном изложении рассмотрение поставленной задачи бу­ дет проведено при следующих допущениях:

1) все действующие на систему случайные сигналы являют­ ся стационарными, т. е. такими случайными сигналами, стати­ стические характеристики которых не зависят от сдвига во времени;

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

3)сигналы ошибок рассматриваются в установившихся ре­ жимах.

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

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

276

Напомним кратко физический и математический смысл этих ха­ рактеристик.

М а т е м а т и ч е с к о е о ж и д а н и е mx [t) случайного ста­ ционарного сигнала x(t) есть постоянная величина, равная среднему значению этого сигнала, и определяется .выражением:

т х (і) = тх = М [х Н)].

(5.78')

Разность случайного сигнала и его математического ожидания является случайным сигналом с нулевым математическим ожи­

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

о

сигналом и обозначается x(t):

 

 

 

 

x{t)

= x { t ) -

mx {t),

 

 

 

 

(5.79)

 

 

 

 

M[x{t)]=-0.

 

 

 

 

 

(5.80)

 

К о р р е л я ц и о н н а я

фу н к ц и я

Кх (х)

случайного

ста­

ционарного

сигнала

x(t)

есть

неслучайная

функция,

характе­

ризующая

взаимосвязь последующего

и

предшествующего се­

чений сигнала х ( і ),

разнесенных

на

время

х,

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

выражением:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K A ') = M[x(t)x{t + x)].

 

 

(5.81)

Д и с п е р с и я

D x(t)

случайного

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

сигнала

x(t)

есть постоянная

величина, характеризующая степень

раз­

броса случайного сигнала

относительно его среднего значения,

и определяется выражением:

 

 

 

 

 

 

 

 

 

 

 

Dx {t) = Dx =

М [°2М].

 

 

(5.82)

Из

(5.81) и

(5.82)

можно установить важное свойство:

 

 

 

 

 

DX = K J 0),

 

 

 

 

(5.83)

т. е. дисперсия равна корреляционной функции при х = 0.

 

С п е к т р а л ь н а я

п л о т н о с т ь

 

(ш)

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

случайного

сигнала

x(t)

есть неслучайная

функция,

характе­

ризующая мощность

составляющих

случайного

сигнала

x(t),

спектральный состав

‘которых

заключен

в

диапазоне

частот

<в -ь- to -f- d и,

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

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

Фурье

от

кор­

реляционной функции:

оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

» =

j/C x ( x ) e - ^ d x .

 

 

(5.84)

— со

Учитывая свойство четности корреляционной функции, мож­ но показать, что Sx (io) есть также четная и вещественная функция, определяемая выражением:

оо

 

S x (a>) = 2 J Кх (^) cos сох dx.

(5.85)

Cl

 

277

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