
книги из ГПНТБ / Пугачев В.Н. Комбинированные методы определения вероятностных характеристик
.pdfфициенты разложения. Далее осуществляется оценка
сверху коэффициентов разложения. Коэффициенты раз ложения сложным образом выражаются через характе ристики упрощенной системы в отсутствие разброса па раметров. Однако иногда абсолютную величину их мож но оценить сверху и эту оценку в достаточно простом виде выразить через обобщенные характеристики упро щенной системы в отсутствие разброса параметров. Оценка сверху коэффициентов разложения позволяет по лучить искомую оценку изменения вероятностной харак теристики за счет разброса параметров.
Поскольку при получении оценок сверху допустимо существенное упрощение системы, то ниже рассматрива ются случаи, когда упрощенная система является линей ной и стационарной.
В настоящем параграфе дано применение изложен ного метода оценки для двух довольно общих случаев.
А . О ц е н к а с в е р х у и з м е н е н и я с р е д н е г о к в а д р а т н ч е с к о го о т к л о н е н и я с л у ч а й н о го п р о ц е с с а в л и н е й н о й с т а ц и о н а р н о й си с те м е за сч е т с л у ч а й н о го р а з б р о с а е е п а р а м е т р о в
Среднее квадратическое отклонение о(а) установив шегося случайного процесса в линейной нестационарной системе при стационарном случайном воздействии и кон кретном значении вектора а параметров системы опреде
ляется соотношением
+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яа г о- |
||||
|
2зО |
|
|
|
|
||
Тогда из (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