Файл: Методы оптимизации в статистических задачах управления..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