Файл: Пугачев, В. С. Основы статистической теории автоматических систем.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 208

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

mi = М [Л, {Y)],

mi (0) = mi0.

(8.39)

(t = 1,2,

. . . . п)

 

Здесь индекс /е заменен на i. Начальными условиями этого урав­ нения язляются математические ожидания переменных в начальный момент времени.

Для получения уравнений относительно вторых начальных мо-

ментов

переменных системы умножим обе части уравнения (8.33)

на yktji

и проинтегрируем по всем переменным в бесконечных пре­

делах. Проводя вычисления, аналогичные предыдущим, и изменяя индексы i = k, j = I, получим следующую систему обыкновенных дифференциальных уравнений относительно вторых начальных мо­ ментов переменных системы:

«,-/ =

М [УД- +

У Д + 25,;],

(8.40)

(i,

j = 1 , 2 , . .

., n)

 

где M — оператор математического ожидания, осредняющий пере­ менные системы. Начальными условиями системы (8.40) являются зна­ чения вторых начальных моментов переменных системы. Уравне­ ния (8.39), (8.40) являются обобщениями уравнений, полученных в [18].

Аналогично вычисляют моменты произвольного порядка. Вы­

числим смешанный начальный момент N-vo порядка:

 

ar1rt ...rn = M[Y'1'Y ? Y '> ...Y rn],

(8.41)

где гъ г.2, . . ., гп — положительные

целые числа, такие,

что их

сумма равна величине порядка смешанного момента:

 

£ rv = N >

0.

(8.42)

V=1

 

 

Умножая обе части уравнения (8.33) на уг'у£ . . . угп и интегри­

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

а,.,

 

 

 

 

 

 

 

 

1) М

B..Y-2

П У гр

+

Xl

rir

5..У г1Уг1П

Yrp (8.43)

[

“ 1 Р=1

р .

 

1 /=1

 

4 1

1 Р=1

Р

 

^2)

• •

•»

9, 1,

2, . . .,

Л-),

 

где штрих над суммой означает отсутствие в ней членов при i = j. Уравнения (8.43) являются достаточно общими, сараведливыми для всех автоматических систем, для которых можно вычислить

14*

211


коэффициенты сноса и диффузии. Из уравнений (8.43) следуют урав­ нения для первых и вторых начальных моментов. Полагая гк = О, при i =j=k и >'i = 1 из уравнения (8.43) получаем уравнения для пер­

вых начальных моментов: mi = aooo...r1.=i...oo = М [Л(.]. Анало­ гично можно получить уравнения для вторых начальных моментов: полагая гк = 0, при i =j= k и к ф j и = /у = 1 из уравнения (8.43) имеем

а ц = а оо...г£= 1.../-/= 1...оо = М [ ^ И / 4 “ ^ i^i ~Ь 2 5 £;-],

что соответствует уравнениям (8.40).

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

S гМ

[Л,УГ«+ BuY f

(rt - 1)] П У*

+

i=l

l

P — 1

 

+

2 '/ у . м a i

П Y rp = 0.

(8.44)

0=1

 

Установившиеся значения моментов первого и второго порядка

определяются уравнениями (8.39), (8.40) при т1 = а £/- = 0.

Если нелинейные функции ср (у) могут быть представлены в виде степенных рядов

фik( y ) = S

1А 'У 2 ---УЪ

(8-45)

•••.

чп = о

 

то коэффициенты сноса и диффузии, как это следует из формул (8.23), (8.31), будут также являться степенными рядами переменных. Учи­ тывая, что в уравнениях (8.43) коэффициенты сноса и диффузии ум­ ножаются на переменные и далее осредняются применением оператора математического ожидания, можно отметить, что моменты в правой части уравнения (8.43) будут более высокого порядка, чем в левой части. Для определения этих моментов необходимо записать новые дифференциальные уравнения соответствующего порядка. Однако они вновь будут зависеть от моментов еще более высокого порядка. Следовательно, для определения всех моментов необходимо соста­ вить бесконечную систему уравнений. Частный случай этой системы уравнений, соответствующий отсутствию мультипликативных воз­ мущений, получен В. С. Пугачевым [57].

Рассмотрим линейные системы. Подставив выражения для коэф­ фициентов сноса (8.25) и диффузии (8.26) в уравнение (8.43) и приме­ нив операцию математического ожидания, получим следующую си-

212


стему уравнений относительно смешанных начальных моментов N-to порядка [76]:

“ Г1■■■гп ~

.

S

 

 

+

+ ~2 Gfqqk) <Zrv ..ri-:l.. .rk+l.. ,rn +

+ ^~2 bkqG jkq

+

b t k t n f 'j a ri...

r f—l-

...rnJ +

- J .

t r i (r I

'— 1) [G fk iq 'X

ari...

r.-2...

rk+l...

rq+l

...rn + 2Gfkqb£qari...

rk+i...

ri-2

...rn +

 

 

 

 

г

 

 

 

1

" ,

 

 

 

 

H- b i j t b iq G kqCCr

.,

_ 2 .

. . r }

4 “ ~FT

i j

t i ? j

{ G i k j q X

 

X

. .r(—1.■■rj—1■■•r/.+I- ■•r9+[- • -Ot "Г" (Gikqbjq -f- Gjfcqbiq) X

 

Хау-.гь+1.-гг 1-..Г/-1".гл +

bikbiqGkqar^

...r-—I.. ,r;—I. ••/•„).

(8.46)

 

 

(rj, r 2, . .

 

rn =

0, 1, . .

N)

 

 

 

Знак штрих у суммы означает отсутствие слагаемого i =

/. Если

в системе отсутствуют мультипликативные возмущения, то из си­

стемы (8.46)

следуют уравнения

 

 

 

 

П

 

,k+i rn + biltm^arv..r.-i

rn) +

К

-

ri ( ^ r 1..,i- i ..

 

 

+ Т S й ( й - 1 ) ^ Ь ц Х

 

 

 

^ i ,

k , q—l

 

 

X G k q C l r

 

г - 2 . . . Г + 4

IS'

nrjbibbjfikjXr г 1...Г 1 .

. , . (8.47)

 

 

*

i , j , k, <7=1

'

 

Для моментов первого и второго порядка из системы (8.46) полу­ чаем

^

[ ( a‘ft

m‘k Т

mk 4~

ft, <7=1

 

 

 

+ - у b k q G f k q + biktrLk

(8.48)

acj= У] [(a ;A+ mfA+

у Gfw^

ос/г/+

+ ) m 7+

ft, <7=1

 

 

 

4 ” ^ a /ft 4"" m /ft 4 “ -g G j q q k ' j a fti 4~

b k q G j k q + b j ^ t n ’k ^

/I

+ y X

X ^

\ ( G i l i j q

G j k i q ) Clkq

( G i k q b j q

G j b q b i q ) П 1 к _ - \ - b i k b j q G k q ] - (8.49)

ft, <7=1

Анализ уравнения (8.48) показывает, что коррелированность ад­ дитивных и мультипликативных возмущений приводит к появлению

213


эффекта детектирования, т. е. при отсутствии математического ожи­ дания аддитизного возмущения = 0 математическое ожидание

выходного сигнала не равно нулю в установившемся режиме, а про­ порционально взаимной интенсивности аддитивных и мультиплика­ тивных возмущений. Таким образом, воздействие параметрических возмущений, коррелированных с аддитивными возмущениями, при­ водит к появлению систематических ошибок выходных переменных. Это явление в частности наблюдается в измерителях при недостаточ­ ной жесткости конструкции измерителя, в системах с переменной структурой и самонастраивающихся системах. Уравнения (8.48), (8.49) являются обобщением уравнений для моментов [79] на случай параметрических возмущений.

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

П

 

Ф/Л- (у) = Фо,А + Y kikp (Ур — тр),

(8.50)

р=1

завися­

где ср0гй, kikp— параметры статистической линеаризации,

щие от первых и вторых моментов переменных. Подставляя это вы­ ражение в формулы (8.24), (8.31) и далее подставляя коэффициенты сноса и диффузии в уравнения (8.39), (8.40) и применяя операцию

математического

ожидания,

получаем

 

 

 

 

 

 

 

 

П

 

 

 

 

Л

 

 

 

 

 

 

 

ПТ-i = Yli^ik

 

 

 

9

 

kikp

\^0PqGikPq

Ь pqGi)tq\ -{-

 

k = l

 

 

 

z

k , p, q—I

 

 

 

 

 

 

 

 

 

 

 

+

t b ikmN(ty,

 

 

 

 

(8.51)

 

 

п

 

 

Г

A=1

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

aij = Y i an +

niu)

фоп/п/+

Y

kUp (api — mpmj)

+

 

 

 

i=i

 

 

L

 

 

p=i

 

 

 

 

 

 

 

Yl

kilP

ф0Pqfftj

 

Y

k pqr (a

r j

n i r i r i j )

G u Pr +

~2

G fiNpm } j

+

l, p .q = l

 

 

■r —l

 

 

 

 

 

 

 

 

 

mkЩ +

Y

( au + n4i)

фоцтк +

£

kjip (aPi rriptrii) +

 

 

 

 

i=i

 

 

 

 

 

9 =

1

 

 

 

 

+

Y

kjtq

ф0pqfll-i

“i- 9

Yi

k pqr (pCr i

t1lrtH i)

G jlp r +

-o-Gy/p/Hi

f +

l. P. <7=

I

 

 

r= l

 

 

 

 

 

 

 

 

+

m fn ii+

S

фопфо/г +

Y

kiipkjrq {ctpq — mpmq)

G i ljr H~

 

 

l. r= 1

l L

 

 

p , q—l

 

 

 

 

 

 

 

 

 

 

 

 

Gifjfpou -j- Gju"фо//1+ Gij.

 

 

(8.52)

 

 

 

 

(i, / = 1 , 2 , . . . ,

n).

 

 

 

 

 

214


Нелинейные уравнения (8.51), (8.52) являются обобщениями урав­ нений для моментов [1]. Нелинейность уравнений (8.51), (8.52) обусловлена зависимостью параметров статистической линеаризации от первых и вторых моментов переменных.

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

Удобно построить систему уравнений относительно кумулянтов переменных. Число уравнений для кумулянтов в случае нелиней­ ных систем также будет бесконечным. Однако с увеличением номера кумулянта его значимость уменьшается, поскольку высшие куму­ лянты характеризуют отклонение закона распределения вероят­ ности от нормального. Для получения системы уравнений N-ro по­ рядка достаточно положить равными нулю кумулянты N + 1-го по­ рядка. Точность решения задачи анализа систем с параметрическими возмущениями можно оценить, если вычислить дополнительно N + 1-й кумулянт.

По определению, кумулянтом N-ro порядка случайного вектора

называется величина

 

 

(8.53)

где g (к; t) — одномерная характеристическая функция перемен­

ных системы, гъ г2, . . .,

гп — положительные целые числа такие,

что их сумма равна N >

0.

Как известно, кумулянт первого порядка совпадает с математи­

ческим ожиданием. Кумулянт второго порядка равен второму цен­ тральному моменту. Поэтому из уравнений (8.39), (8.40) получаем

«i =

М Д ];

(8.54)

( i = l , 2, . . ., п)

 

= М [У Д + У Д +

2В„] - xtM [А,] - щМ Д ] .

(8.55)

Уравнения (8.54), (8.55) образуют совместную систему. Аналогичным образом можно получить кумулянты высших поряд­

ков. Так, для кумулянтов третьего порядка получаем следующее уравнение:

В данном уравнении производные по времени от кумулянтов первого и второго порядков в правой части можно исключить, поль­ зуясь уравнениями (8.54), (8.55).

215