Файл: Методы оптимизации в статистических задачах управления..pdf

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

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

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

Добавлен: 21.10.2024

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

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

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

Рис. 24. Приближенная ма­ тематическая модель случай­ ного процесса:

а — схема формирования муль­

типликативной помехи;

б — ма­

тематическая

модель

процесса

а (0;

в — схема

реализации

мультипликативной

помехи

 

а (i)

x t

(t)

 

■alt)

6лох

умножения a(t)xL(t)

а)

+N

t p ---- тъ2 a(t)

' * г

-/V u

S)

+ H

■a(t)xi(t)

- N

&)

Процесс а (f) является случайным, стационарным с равным нулю математическим ожиданием М (01 и корреляционной функцией

R ia(0 =

' I.

(197)

Таким образом, система (196) содержит мультипликативную помеху а (t) xt (t) (рис. 24, а) по і-й координате выходного про­ цесса. Для случайного процесса а (t) может быть построена сле­ дующая модель (рис. 24, б). Реле в случайные моменты времени переключается из положения 1 в положение 2 и обратно. Пред­ положим, что положения 1 и 2 равновероятны, а распределение моментов переключений подчиняется закону Пуассона со средней частотой р:

P(n, 0 =

где р (n, t) — вероятность того, что в течение интервала времени t произойдет точно п переключений. Построенный таким способом процесс имеет корреляционную функцию (197) и математическое ожидание, равное нулю, что доказывает возможность представле­ ния процесса а (і) моделью рис. 24, б. Это позволяет заменить схему рис. 24, а схемой рис. 24, в, Содержащей блоки постоянных коэффициентов + N и —N , а также контакты случайным образом переключающегося реле. Таким образом, система управления с мультипликативной помехой заменена системой со случайным скачкообразным изменением состояний. В первом состоянии она описывается дифференциальным уравнением

X == Ä qX —j- Â ^ x —I- H f ,

75


где

0 0

.

0

 

 

 

 

 

'

0

 

 

 

6

 

 

N

 

>

 

—— — Nx,

 

0

 

 

 

0

 

 

0

’ .

 

 

 

 

 

 

'

0

 

6

 

во втором — уравнением

 

 

 

 

 

X — А ах + А гх +

Bf,

 

 

где

 

 

 

 

 

 

0

 

 

 

 

0

 

.

'

 

0

 

 

 

 

' 0

 

 

 

ö

 

4 =

+ N

 

f

 

Nxt

 

 

 

0

 

 

0

 

 

0

' .

 

 

 

 

 

 

 

0

 

Ö

 

Поскольку интенсивности

переходов

из

первого

состояния

во второе и обратно одинаковы и равны

р,

матрица

А (t) имеет

следующее выражение:

 

 

 

 

 

 

 

Л(і) =

— I*

М- I

 

 

 

р

— р II

 

 

 

 

 

 

Математическбе ожидание и дисперсионная матрица выходных координат рассматриваемой системы могут быть определены изло­ женным выше методом. На рис. 25 представлена схема модели­ рования М [х (t)}.

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

76



Рис. 25. Схема моделирования М [х (t) ] в системе с муль­ типликативной помехой

мости статистических характеристик выходных координат си­ стемы от старших моментов мультипликативной помехи. Для учета старших моментов мультипликативной помехи следует применить более сложные модели.

3. Линейные системы со случайными коэффициентами типа «белый» шум

В предыдущих параграфах было проведено исследование ли­ нейных динамических систем со случайными возмущениями коэф­ фициентов, когда эти возмущения представляют собой случайные величины или марковские чисто разрывные процессы (например, телеграфный сигнал). В первом случае вся мощность процесса отклонения параметра сосредоточена на нулевой частоте его спектральной характеристики, во втором случае случайный про­ цесс имеет более богатый спектральный состав. В этом параграфе будут рассмотрены возмущения параметров, имеющие равную интенсивность всех гармоник, или процессы типа «белый» шум. Статистическому анализу линейных динамических систем со слу­ чайными коэффициентами типа «белый» шум посвящен ряд работ [118, 121, 123, 131 ], в которых для решения задачи анализа приме­ нен аппарат теории процессов Маркова и уравнения Колмогорова.

Рассмотрим линейную динамическую систему, вектор выход­ ных координат которой х (/) удовлетворяет системе уравнений в векторно-матричной форме

X = А (t) X + f (t),

 

где X— вектор-столбец п измерений; матрица

А (t) размерности

[п, п] с элементами als (t), i, j =

1,

2 , . . ., n,

характеризующи­

мися статистическими свойствами

 

 

 

М Іаи (/)]

=

аи (t)\

( 198)

М [аи (t) akl (т) 1=

 

 

 

ki (і) б (t — г);

77


/ (/) — и-мерное входное возмущение со статистическими харак­ теристиками

M[ f ( t ) ] = 7 ( 0 ;

I

(199)

М [/ (/) /* (т)] = N (/) б (/ — т). (

Случайные процессы Л (/) и f (t) коррелированы:

М lau (i) °fk (т) ] = G(i, k (t) б (t — т).

(200)

При принятых предположениях относительно свойств случай­ ных процессов А (/) и f (/) [см. формулы (198)—(200)] случайный вектор X (/) является марковским и его одномерная плотность распределения вероятностей р (/, х) удовлетворяет уравнению Колмогорова (108):

Зр

( ж ) '

*Ъ + - Т ІГ

 

(201)

dt

дхдх*Р® V' ХУ

 

Если известно начальное распределение вероятностей вектора выходных координат х (/):

Р (0, х) = ро (X),

(202)

то уравнение (201) также имеет начальное условие (202). Определим вектор коэффициентов сноса с (/, х) и матрицу

коэффициентов диффузии Ѳ (/, х) рассматриваемого процесса х (/):

 

c(t, х) — lim -r-M [Дх(/)|х(/) = х] —

 

 

 

 

 

Д->0 л

 

 

 

 

 

 

 

f

 

 

 

 

)

 

= й 4 -'м J [Л (т) * w + / м

( о = * =

 

lim_J_

М

і

} dx

А (т)

х(/) + ]^Я[Л(Я)х(Я)+/(Я)]

+

д^о д

 

<

 

L

Lі

 

 

\

/( т) X (/) =

xj

= Ига -д- I |^Л (/) X +

/ + - у Д м (/) X -ф-

л/(0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(203)

В формуле (203) для вектора коэффициентов сноса с (/, х)

матрица Л (/) имеет элементы

atj (і), / , / = 1 , 2 , . . ., п;

матрица

(/) и вектор d^ (/)

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

 

 

 

 

 

 

jД_

г+Д

т

 

 

 

Daa (t) =

lim

М J dx

J dЯЛ (x) А (Я);

 

 

 

 

 

Д->0

I

t

t

 

 

 

dAf (0 =

Hm

tf+A

T

 

 

 

T

M J dx

j dXA (x) f (Я),

 

78