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

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

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

фициенты разложения. Далее осуществляется оценка

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

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

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

А . О ц е н к а с в е р х у и з м е н е н и я с р е д н е г о к в а д р а т н ч е с к о го о т к л о н е н и я с л у ч а й н о го п р о ц е с с а в л и н е й н о й с т а ц и о н а р н о й си с те м е за сч е т с л у ч а й н о го р а з б р о с а е е п а р а м е т р о в

Среднее квадратическое отклонение о(а) установив­ шегося случайного процесса в линейной нестационарной системе при стационарном случайном воздействии и кон­ кретном значении вектора а параметров системы опреде­

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

+00

а2 (а )= J /с (/'to; а)/с (— /со; a)G(w)du>.

(5.9)

—00

 

Здесь К (jeо; а)-— частотная характеристика системы при конкретном значении вектора а параметров системы,

aG(co)— спектральная плотность воздействия.'

Ві[29] показано, что для некоторого класса нестацио­ нарных воздействий, выраженных через неслучайные и стационарные случайные функции времени, определение средних квадратических отклонений установившихся

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

Среднее квадратическое отклонение а с учетом слу­

чайного разброса параметров будет находиться из соот­ ношения

з2 = | 32(a)f(a)da,

(5.10)

21 0

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

Математическое ожидание вектора и обозначим через

<хо. Тогда среднее квадратическое отклонение (То процес­ са в системе без учета случайного разброса параметров будет

ао=<т(ао). (5.11) Произведем оценку сверху относительного изменения б среднего квадратического отклонения процесса за счет

разброса параметров.

Значение б определим по формуле

 

8 = 1 = - ^

(5.12)

°0

 

и найдем оценку свеху абсолютной величины б. Предполагая компоненты Aaq случайного отклонения

Да вектора параметров малыми настолько, что величину о2(а) можно разложить в ряд по степеням этих откло­

нений и ограничиться членами второго порядка малости, из (5.10) и (5.11) получим

2

О

 

дг°2К)

^

 

а

 

++SS daqod<tro

4 j “r ’

(5.13)

 

 

0=1 г- 1

 

 

где п — число

компонент

вектора а

(число

случайных

параметров); щ 0 и а>0— математические ожидания q-то

и г-го параметров;

К а2аг— корреляционный момент слу-

чайных отклонений ^-го и г-го параметров.

 

Введем корреляционный момент относительных от­

клонений <7-го и г-го параметров:

 

 

 

 

Да5

Лаг

 

к

(5.14)

 

а 9о

 

а го

 

а 2оа ,*о

 

 

 

 

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

 

д2а2(а„)

 

 

 

dqr

 

 

 

(5.15)

2

да2оЙа

 

а 7яа г о-

 

О

 

 

 

 

Тогда из (5.12) — (5.15)

найдем,

что

 

 

 

П

П

 

 

 

 

/

1 +

S

E

а яqri\г К

qr 1.

(5.16)

 

 

 

 

 

 

q~rl Га!

 

 

 

 

14:

211

Отсюда с учетом малого влияния разброса параме­ тров нетрудно получить, что

ПП

(5 Л 7 >

<7 = 1 г= 1

Абсолютная величина б удовлетворяет неравенству

П П

(5-18)

q=\ r=l

Обозначим через A qr оценку сверху абсолютной вели­ чины коэффициента aqr:

|а<7г|^АдГ. (5.19)

Сучетом (5.18) получим оценку сверху величины б:

ПП

 

 

і8!<4 SS^r|K?r|-

(5-20)

 

 

 

 

q=\ г—1

 

 

 

 

Для практического использования этой оценки необ­

ходимо знать величины A qr.

 

 

 

 

Перейдем к определению величин A qr, удовлетворяю­

щих неравенству

(5.19).

 

 

 

 

 

В соответствии с (5.9), (5.11) и (5.15) коэффициенты

ciqr

можно

представить

 

в следующем

виде:

 

 

 

 

+ 0

О

 

 

 

 

 

 

 

J

фдг (/<о;

«о) G (“ ) da»

 

 

 

 

 

---------------------------------- (5.21)

 

 

 

^ Ко (/со;

а0) Ко (— /со; а 0) G (со) rfco

 

где

 

а,0аг

 

 

 

 

 

 

*^4 (/«>; аі):

 

 

 

 

 

 

 

(/“*: *о)А'о ( — /®;

а о) +

+

А '> ( 4

а 0) Д 0аго( - /с о ;

а 0) + ;А 'оГО( 4

а оЖоЯо( —/®; ®о)+

 

 

+ 4 0 4

*o)^oeo“re(-/®;*o)]-

(5.22)

 

 

 

 

В формуле (5.22)/С”30 г° ( 4

а„), K0q°(i<°l а0)

и т- Д-

обозначают

производные

от

К0(/«*;

а0) по

величинам

212

Для получения оценки A qr удобно формулу (5.21)

записать следующим образом:

 

+ 00

 

 

 

 

 

“0)/С0(— /и; a0)G(<o)'du

 

 

 

J

ФдД/ю;

«о) /С„ (/«;

 

(5.23)

<Ѵ=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

/Со (/со: а0) /С0 ( — /со; «0)]G (и )

Ао

 

 

 

 

Здесь

 

 

 

 

 

 

 

Щг{Іи; ао)_____

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ф?г (/®;

во)== /<0 (/со;

«о) /Го (— /«о;*а0)

 

 

 

 

 

а2оаго

/С0Ѵ ™(/«а;

«о)

І- '

Ко {—іи;

« о ) ,

 

 

 

 

2

 

 

К0(/со;

а0)

 

а„)

1 -

 

,

/ Г р 10 ( / с о ;

а 0)

/ С

о

(— /to;

а

0)

/

С ^

Д — / с о ;

“о)

К І г 0

Ц и ;

а 0 )

1

/С0 (/öj;

 

а0)

/Со

(— /со;

а0)

Ко {—Іи;

а„)

/Со (/со;

а„)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.24)

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

Go(co) = K O(/Ö ; ао)Ка(— /со; ao)G(co).

(5.25)

Оценка второго типа учитывает, что спектральная плотность Go(со) начиная с некоторой частоты £2 прак­ тически равна нулю, т. е.

Go((o)~OnpM | c ö | ^ Q .

(5.26)

Для обеих оценок из (5.23) и (5.24) можно получить

I ^-Aqr — Bqr^rBnBr-

(5.27)

Здесь для оценки первого типа

В,

■■max

 

К

0 Я °

г0 { І и ; а . )

С С:

 

 

 

Jq r -

 

ілл

Го

/Со (/со;

а.)

 

ПО со

 

 

 

 

 

 

 

(5.28)

ВдР* max

/Cp80 (/со;

«о)

 

1,0 /Со

(/<о;

«о)

 

ПО w

213

для оценки второго

типа

 

 

 

Вq r

max

К

(/со;

«о)

а<7оаГо Ко (/to;

а0)

 

I

П О

( О

 

т

I s a

 

 

 

(5.29)

 

 

 

К*8» (/« ;

а0)

 

ß97 max

 

 

Ко (/со;

а0)

 

I

П О

to

 

(О l ^ ß

 

 

 

 

Величины ß g,-, ß g

и ß,- достаточно

просто находятся

[37], если передаточную функцию системы Л 'о (р ; ао) вы­

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

Покажем, как это делается для случаев., когда

(р; <*) = 1

К (р; а)

а)

(5.30)

+

к

( р \

и

 

 

 

1

 

 

(р ; а )

 

+

 

а) *

(5.31)

1

к

( р ;

 

 

где К(р\ а) — передаточная функция разомкнутой систе­

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

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

чала оценку первого типа.

Итак,

пусть Ко(р\ а) имеет вид (5.30). Тогда

 

K qo (/<о; а»)

к“80(/со;

а0)

(5.32)

К0 (/со;

а0)«а?о= [1 — Ко(Іт’> а о)1

К (/со;

а„) а 9о-

Так как К{р\ ао) представляет собой произведение

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

/((/со; a0)=Pq(/co; aq0)Qq(jor, aw, ■■ ag-)-io, . • ., ctno) j

где Pq(jсо; aqo) — передаточная функция звена, зависящая от q-ro параметра;

aq-w,

(5.33)

простейшего

Qq(P> ■■■) = Pq ^p’-

— передаточная функция, не за­

висящая от 9-го параметра,

214

Из (5.32) и (5.33) полупим

 

 

К 80 (/«о;

«„)

P4qo(/со; а2о)

(5.34)

Л' (/со; а.)

— Яв (/со;

а3о) аЧо-

 

Пусть

 

 

 

 

Со

:гаах

я“80 (/со; аЧо)

(5.35)

Яд (/со; <*зо)

*9о

 

П Э О )

 

Значения величин Cq приведены в табл. 5.1.

Имея частотную характеристику системы без учета

разброса параметров, легко найти величину

 

 

 

 

М, =

max 11

~ К 0 Ощ а0) |.

 

 

 

(5.36)

В силу (5.32), (5.35) и (5.36) получим

 

 

 

 

 

До«» (/со;

а0)

80

 

\мхсч.

 

 

(5.37)

 

 

П О ш

Ко (/со; ««)

 

 

 

 

max

 

 

 

 

 

 

 

 

 

 

Сопоставляя (5.28)

с

(5.37),

видим,

что

величину B q

можно определить следующим образом;

 

 

 

 

 

 

 

 

Bq — MyCq.

 

 

 

 

 

 

(5.38)

Найдем B qr. Из

(5.30) получаем

 

 

 

 

 

 

V

0 г 0 ( / “ ; “ о)

 

___д “ « » “ '-»

(/со;

а0)

 

 

 

 

Ко (/со;

а0)

а Ч о а

г 0

 

к ( ; Ч

о .

а о

)

а

? о а

г

„ л .

 

1_____

, У

80 (/со;

“ о )

_

У

г0

(/со;

а

0)

X 1 +

К (/со;

0С0)

 

Я (/со; а0)

а0„ --,TT—■—

- аГо Д

 

"’о

Д (/со;

«„)

^ ■

 

 

 

Хі

К

(/со;

а0)

 

 

 

 

 

 

(5.39)

 

 

 

I -Ь Д (/со; а0)]*‘

 

 

 

 

Учитывая вид передаточной функции К(р\ а0), можем

написать,

что

 

а0)=Pqr(j(o\ aq0, aro) X

 

 

 

 

K(jw,

 

 

 

 

XQqr(/со; ( Х м ,

■ .

CCg-10,

Q

. q +

10,

■ ■

 

 

 

 

 

. .,

(Xr—10, Ctr+10,

• • •,

Ctno),

 

 

 

(5.40)

где Pgr (p;

aqo,

ar0) — передаточная

функция,

зависящая

от q-то и

r-го

параметров, Qqr(p\

. • . ) — передаточная

функция,

не зависящая

от q-то и r-го

параметров.

215

Отметим, что Pqr(p; ago, aro) может быть либо про­

изведением двух передаточных функций простейших звеньев, зависящих от q-то и r-го параметров соответст­

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

В силу (5.40)

К 90

г0 (/со; о)

Pqr° Г0 (/ы> “ЗО’

“ го)

-г . . .

Л’ (/<о; «„)

а?»аг» —

Р іг (/to; а9о, аг0)

<5-4 >

Если

/V (р;

а?0, аГо) — произведение

передаточных

функций простейших звеньев,

то

 

 

Ф“" (/«; ад„, ar„)

Ядг (/со; ago, a r„) a?o^ro —

 

я ,90 (/<■>;

“«о)

Nx

prr0 (/“; aro)

(5.42)

 

~~ Яд (/со;

ago) a’o

X

P r

(/со; a r0) a»v

 

 

Из (5.41) и (5.42)

в силу

(5.35)

получим

 

шах

^ “Зо“го (/»; а0)

 

 

 

я990(/to; aSo)

 

Я (/со; а0)

 

 

шах

X

ПО СО

 

 

П О

С О

Яд (/со; ago) a ’ °

 

я“г0(/со;

аг„)

a Н---п п

(5.43)

 

X max

яг(/со;

аг0)

 

 

о — ^9^г*

 

П О с о

 

 

 

Для рассматриваемого случая произведение CgCr

удобно

обозначить через

 

 

 

 

 

 

CqT= C qCr.

 

 

 

(5.44)

Если Pqr{p\ ctgo, aro)— передаточная функция одного

простейшего звена, то можно ввести величины Cqr, вы­ числив их по формуле

я“9»“г0(/со; ago, ar0)

 

С?г = шах

ЯдД/со;

ago, a r0)

a Voa r<>

. (5.45)

 

 

' д г О ш .

“ go.

 

Значения

Cqr для этого случая приведены в табл. 5.1

и 5.2.

 

 

 

(5.45)

получим

Из соотношений (5.39), (5.43) или

шах

К0> “г0 (/«;

*») - -

^ С чгМ ^ 2 С чСгМ ,М г,

П О с о

 

 

 

 

 

где

 

 

 

 

(5.46)

М2=

шах IК0(/да; а0)

 

(5.47)

 

 

П О О )

216

Т а б л и ц а 5.1

К

о,

+

ся СЯЛ\

СЯ СЯ

V/ А

ся ся

Я 5

О. g*

с с

к

:1л

о,

+г

+к,

ил Ч ОА ил

СЯ СЯ

СЯ СЯ

V/

А\

V/

л\

ся

ся

ся

ся

к

s

к

а

Си Си

a,

а.

в

в

с

а

 

 

■а

w

Е>,

MJ>

 

 

к

 

+

fc

+ P

L

s-ч

\

к

V/

і Х Г і

'

S?

JJ

I

fljji

+

 

 

 

a ся

V/ А\

СЯ СЯ

5 к Э* о.

с с

'2, « тз и <м

а СЯ

ѵ /А\

ся а s м

с

■«

^ ІСЧ

СЧ

+

2)7

 

 

 

Т а б л и ц а 5.2

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

Параметры

Сіг

С

“«о и “го

 

С Ч г

1 + 2\Тр + Г Ѵ 2

2IT и Т2

0

0

Ѵ~2

К2

1

2 IT и Т%

 

2 bd при

£2 < 2 8

1 + 2 6 Tp+T'p '

/

 

f при

S ГЗ: 2 3

Ѵ~2.

 

 

 

 

 

 

К2

Пр и м е ч а н и е . В табл. 5.1 и 5.2.:

ЙГ

 

 

4(2£3— 1)з + 9 (2 ^а— 1) + [4(2 £ а— 1 ) а + 3] +

V ( 2 ^ — 1 ) 4 - 3 .

V1+ WT*

 

 

(2£3—1)» + 1+ (2^2—І)К(2^_ 1)2 + 3;

Ь= ,

 

:

й, = -^г- ;

 

 

 

К(1—8>Л)а+(2уя)а

 

'

 

 

 

с= --- 1

:

 

йа= --- :

¥

 

 

24, V I — £3

 

 

тѴ\ —2

 

 

d _______№ ______

0 _ V (2^»—1) + Ѵ(2^—Гр+ 3

 

Ѵ(і—я^т-игУв)»

 

а~

т

'

 

Из сравнения

 

(5.28) с (5.46) видно, что

в качестве

Bqr можно принять

 

 

 

 

Bqr— CqyMi + 2CqCrMiMz-

 

(5.48)

Используя (5.27), (5.38) и (5.48), из (5.20) получим оценку первого типа для рассматриваемого случая, когда передаточная функция системы определяется соотноше­ нием (5.30):

ПП

і8і<4 S S {Cqr+{ Ш *+Мі)CqCr]1Kqr 1 (5-49)

<7= I ur = l

Аналогично получим оценку первого типа для случая, когда передаточная функция системы определяется по формуле (5.31):

Пп

I3I<"ГS S Мг [c,rt Ш f£r\ I КчтI- (5.50)

ч—\ г=]

218

Из (5.4Ö) и (5.50) видно, что найденные сверху оцен­ ки увеличения среднего квадратического отклонения про­ цесса достаточно просты. Эти оценки определяются через корреляционные моменты Kqr относительных отклонений

параметров, максимумы амплитудных частотных харак­ теристик системы при номинальных значениях параме­ тров и коэффициенты Cq и Cqr, приведенные в табл. 5.1

и 5.2.

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

реписать выражение (5.23)

для коэффициента aqr с уче­

том (5.26):

 

 

 

+ S2

 

 

 

J Ф<гг (/<*>;

®о) К> (/и;

ао) К0(— /“>;

aoVß (<*>)dw

а?г = —— ур----------------------------------------------. (5.51)

I Д

0 ( / с о ; а 0) Д " 0 (— / с о ; а 0) G

( со) d a

Я

 

 

 

Далее, поступая таким же образом, как и при полу­ чении оценки первого типа, легко показать, что оценка второго типа будет иметь аналогичный вид, с той лишь разницей, что величины М и М2, Cq и Cqr необходимо за­ менить на М М'г, С'д и C'qr. Причем

М \ =

max I 1 — Д0(/со; а0)|,

 

 

П О

С О

 

 

 

I 00

|S=8

 

(5.52)

/И'2 =

max |/С0 (/<*>; а0)|,

 

 

П О

с о

 

 

 

I с о

 

 

 

а С'д и C'qr имеют значения,

приведенные в табл. 5.1

и 5.2.

 

 

 

 

Б. О ц е н к а с в е р х у в л и я н и я с л у ч а й н о й в з а и м о с в я з и ка н а л о в

 

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

у п р а в л е н и я

на е е то ч н о с ть

 

В некоторых случаях исследуемая система может быть отнесена к классу двухканальных систем автомати­ ческого управления [4, 27].

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

219

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