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

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

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

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

Добавлен: 21.10.2024

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

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

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

+ kM (ty -

T) °x2 (t2 - Я)] +

kM [x2 Ui -

T) ф (t2 - Я)] -b

 

=

+

M [cp (Ц — t ) ф (t2— Я)] dxdX =

 

1 120 )! (t )

(Я) Ң

— — + t ) dxdX',

 

 

fl 1

 

 

 

 

 

 

 

о о

 

 

 

 

 

 

 

 

fl

f%

 

 

 

 

*22 Ui,

^2) — J I W2 (x) w2 (X) kn (-Ц — T,-

t2 — X) dxdX\

(9 0 )

 

 

о

0

12

 

 

 

 

 

 

 

 

 

 

 

*12 Ui,

^2) = J w 2 U) *11 (^i> ^2 — Я) йЯ;

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

*21 (<1. У =

J ш2 СО *11 (*1 — т-

^2) dx.

 

о

Математические ожидания, входящие в уравнения (89) и (90), при совмест­ ном нормальном распределении выходных координат для нелинейной функции (72) имеют вид:

 

 

Фо U) =

т1 (0 т 2 (0 +

d12 (0;

 

М

ixi Ui) ф

(^2)]~ т

2 (^2) *11 Uit

U)

mi Uz) *12

М [ф (Ц) Хі

(t2) ]= Ші (ti) k2i (t1,

t2j“b

tn2 (Ц) kn

M

[x2 (Ц) Ф

(t2) ]= tTlx Uz) *22 Ult

U)“b

mZ(^2) *21

M

о

о

 

m2 (Ц) k12 Ul, tè't

[ф (Ц) х2 (t2)] = mi Ul) *22 Ul, h) +

 

M

[ф (Ц) ф (t2) ] =

тг (Ц) mi Uè

*22 Ult U) +

 

 

+ m2 Ui) m2 (t2) kn Ui,

t2) +

+

m2 (Ц) mi Uz) *12 Ui,

t2) + mi Ui) m2 (t2) k21 (h* tz) +

 

~b *11 Ui, t2) k22 (/j,

t2) -f- ki2 (ti,

t2) k2i (ti, t2).

(91)

(Ц>и)< Ui<U)t Ultt2)'t

(92)

Уравнения (89), (90) с учетом выражений (91) и (92) были проинтегрированы

методом квадратурных формул на ЦВМ с шагом А =

0,1 для следующих пара­

метров системы:

Тг = 0,5, Т2 — 1,

k = 2, г =

1, Q — 0,58. Результаты

решения для

(^, t2), mXi и DXi

приводятся на рис. 8, 9, 10. На рис. 9 и 10 для сравнения показаны математическое ожидание и диспер­ сия выходной координаты Xi (t), определенные методом статистиче­ ских испытаний по тысяче реали­ заций (кривые 2).

thC

Рис. 8. Поверхность корреляцион­ ной функции процесса хг (t) нели­ нейной системы, содержащей дву­ мерный нелинейный элемент

30


Рис. 9. Математическое ожидание про-

Рис. 10. Дисперсия процесса хх (t)

цесса Xi (t) нелинейной системы:

нелинейной системы:

1 — метод квадратных

формул; 2 — метод

1 — метод квадратурных формул; 2

статистических

испытаний

метод статистических испытаний

4. Применение теории процессов Маркова к анализу, непрерывных систем управления

Широкий класс систем управления описывается дифферен­ циальным уравнением

x = f{t, х) + G (t, х) I (0,

(93)

где вектор-столбец фазовых координат системы х имеет размер­ ность п\ f (t, х) — вектор-функция п измерений; G (t, х) — ма­ трица размерности [п, т). Случайные возмущения, действующие на систему, представлены в формуле (94) в виде вектора-столбца I (t) т измерений. Случайный векторный процесс £ (t) является нормально распределенным белым шумом со следующими ста­ тистическими характеристиками:

М[1 (01 = 0;

 

M lU t, ) g*.(fa)] = Q (ti) б (Ц — t2).

(94)

Анализ системы, описываемый уравнением (93), предполагает вычисление некоторых показателей точности. Достаточно общей формой такого показателя является

т

1 — М j ф (t), t] dt,

о

где ф [х (t), t] — скалярная функция фазовых координат сис­ темы. Обозначим через р (t, х) плотность распределения вероят­ ностей случайного процесса х (t). Тогда

т

 

I = J dt j dx ф (х, t) р (t, х),

(95)

31


где J dx означает я-Кратный интеграл по всем компонентам век­ тора X в пределах от —оо до -foo, т. е.

СОJ dxt соJ dx2. .

В частном случае

Ф Ix (i), t] — X* (t) V (t) X (t) +

+ X* (0 Лх (t) 6 {T — t),

где V (t) я А — положительно

определенные матрицы размер­

ности [я,

я], а 8 (Т і) — дельта-функция Дирака.

 

Тогда

показатель

точности

системы

/ равен

 

 

 

/ =

М

JX *

(t)

V (t) X

(t)dt +

 

 

 

 

 

О

 

 

 

 

 

 

 

 

=

т

(t) V (01 dt

 

 

 

+ ** (Т) Ах (Т)

Jtr

+

tr (Т) Л].

(96)

 

 

 

о

 

 

 

 

 

 

В формуле (96) матрица Г (t) размерности

[я,

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

выражением

 

 

 

 

 

 

 

 

 

Г (0 = М [х (0 X* (01,

 

 

 

а

 

 

 

 

П

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

trC =

S сп

 

 

 

 

 

 

 

 

 

г=і

 

 

 

 

является следом матрицы С размерности

[я,

я]

с элементами

си ,

і, j = 1 , 2 , . . ., я.

 

 

 

 

 

 

 

 

Показатель точности, определяемый выражением (96), зависит только от второй начальной моментной функции Г (t) случайного процесса в то время, как в общем случае [см. формулу (95)] показатель / зависит от закона распределения случайного про­ цесса. Расчет плотности распределения вероятностей р (t, х) случайного процесса х (t) может быть проведен с применением

теории

процессов Маркова.

 

 

 

 

 

Случайный процесс х (t) называется процессом без последей­

ствия,

или процессом Маркова, если для любых

<* t2 <j • • •

• • • <j

tk

при произвольном

k

 

 

 

 

 

 

F {tin

xk \ t j,

Xj,

12, X2, • . M

Xft_i) —

где

 

 

F (tk, Xfe I 4-i,

 

 

(97)

 

F ('Jk,

Xk I ti,

 

12X2;

• • -j

 

 

=

 

 

Xi\

 

 

 

p lx (tk)

xk IX (t^j

Xt.,

/

1, 2,

* è • * & 1]«

32


Условная плотность распределения вероятностей марковского процесса х (t) обладает свойством, аналогичным выраженному формулой (97):

(98)

Приведенное определение марковского векторного процесса X (t) означает, что если производится предсказание поведения х (t)

в

момент 4

на

основании

известных

значений х (4), і

= 1,

2,

. . ., k

1, то

точность

предсказания

будет такой же,

как

в случае известного измерения только в последний момент вре­ мени 4-і- Иначе говоря, вся информация о марковском процессе сосредоточена в последнем его измерении. Можно сказать и так: будущее марковского процесса независит от его прошлого при

известном настоящем

значении.

 

Для марковских процессов

 

 

р (4 Xt IX

( т ) , х ^

s <

t) = р (i, xt \s, xs),

где под р (t, хДх(т),

т < s <

t)

понимается плотность распре­

деления вероятностей в момент t при известной реализации про­ цесса X (т), т < s < 4

Согласно известным свойствам многомерного дифференциаль­ ного закона распределения х (t) [82]

х 2; . . .; 4-2» Х/і- і )- ■■ Р (4> Х 2 І 4> -Н) Р (4> х і)-

(99)

Применяя к выражению (99) свойство (98) марковского про­

цесса, получим при 4

< 4 <( . . .

< 4

 

 

р

(4)

4. -^2»

• • •»

4.

x k)

р (4.

x k \ 4-1> x k-i) P

(4-1. Xk - A '

4-2.

xk-i)• •p (4. *2І 4 .

x d

p (4. x i)-

Таким образом, многомерный дифференциальный закон рас­ пределения вероятностей произвольного порядка k полностью

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

одномерными условными законами распределения

р (4 х) и р (4

х |т , у), которые содержат, следовательно, полное

статистическое описание марковского процесса.

имеет обоб­

Важное значение в теории марковских процессов

щенное уравнение Маркова

[8, 23]

 

 

p (4

X | t , y) =

 

 

 

T < s <

4

3 A. M. Батков

33


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

lim

Р [\х (t +

Д^) — x (t) I > е |х (t) =

х] = 0 ,

дг->о АГ

 

 

 

что означает

 

 

 

 

 

1іп1 Ж

I

dyp(t + A t , y \t , x ) =

0.

 

&t->o

I у_х

I ■> Е

 

Другое, более жесткое, определение непрерывности случайного процесса состоит в следующем: случайный процесс непрерывен, если для любого п-мерного вектора е, все компоненты которого положительны,

lim -jj

j

(у — x f p ( t + At, y\t, x)dy = 0.

At->°

I y—x |>

e

Непрерывные марковские случайные векторныелтроцессы х (t), или процессы диффузионного типа, характеризуются п-мерным вектором коэффициентов сноса с (t, х) и матрицей коэффициентов диффузии Ѳ (t, X) размерности [п, я], которые определяются выражениями:

с (t, х) ~ Ііш ~

M[x{t At) — х (t) \x(t) = x] =

 

Д ^О

 

 

= 1іт4т

f (у — x)p(t-\- At, y\t, x)dy;

(100)

Д£->0

J

 

Ѳ(t, x) — lim -rjM [(x(t -f At)—x (t)) (x) (t -j- At)

д/->о

— x{t))*\x{t) = x] = Y\m^j f (y — x) (y — x)*p(t + At, y\t, x) dy.

ДГ>0 Af J

( 101)

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

Процесс х (/), описываемый дифференциальным уравнением (93), обладает марковским свойством и является непрерывным [8,

34