Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Розанов Ю.А. Теория обновляющих процессов

.pdf
Скачиваний:
2
Добавлен:
23.10.2023
Размер:
4.01 Mб
Скачать

§ 3]

НЕКОТОРЫЕ ПРИМЕРЫ

ОБНОВЛЯЮЩИХ

ПРОЦЕССОВ

29

с тем

же структурным

типом dF (t)

и той же

крат­

ности р, что и максимальная система %j, / = 1, .

. р,

может быть дополнена выбранным ранее элементом х,

и, таким образом,

в Я (|) найдена система уи . . . ,

ур, х

кратности

р + 1 ,

большей,

чем

максимальная

крат­

ность р,

что невозможно

(см.

неравенство

(2.6)).

Следовательно, предположение о том, что Я (г|)=И=Я (£),

является неверным, и тем

самым теорема

доказана.

§ 3. Некоторые

примеры

обновляющих

процессов-

1. Марковские процессы. Мы рассмотрим ниже

два важных класса процессов — марковских

и ста­

ционарных в широком смысле.

 

 

Многомерный

случайный процесс КО =

(!г (0}"'

называют марковским в широком смысле, если для любого t проекции Pfo (t + h) величин £г (t + /г), /г^ О, на все «прошлое» Ht{|) принадлежат замкнутой

линейной оболочке значений

|;(0> г = 1 ,

>п.

 

Ограничимся

расмотрением

лишь конечномерных

процессов. В

этом случае для любого s < K

 

 

 

 

т

 

 

 

/ =

1, • • •,

 

(3.1)

Psh ( t ) =

h

Cij(t,

s)l,(s),

m,

 

 

/=i

 

 

 

 

 

 

 

где коэффициенты Cij(t, s)

легко определяются

по

корреляционным

функциям

 

 

 

 

 

s) =

E g ;(0 M s)’

0

/ =

1 , • ••,

т .

 

Например, в

 

случае

л и н е й н о

н е з а в и с и м ы х

значений | ; (s),

/ = 1, . . . , т ,

когда матрица В (s,

s) =

= [Bij{s, s)} является невырожденной и существует

обратная матрица В (s,

s)-1 =

{bkj (s,

s)}; коэффициенты

в выражении (3.1) можно найти по формуле

Cij(t, s) =

2 Blk (t,

s) bk, (s,

s),

i, j — l, . . . , in.

 

k

1

 

 

Согласно

общему предположению (1.1)

limPsg,(0 =

E<(0.

 

 

S~>t

 

 

 

30

 

ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ

[ГЛ. I

и

в

случае линейно независимых значений

£; (0,

г =

1,

т , имеется некоторый интервал

 

такой, что величины

 

 

 

 

Л^ (*^) == ^ sh (0* ^

1) • • • >

 

образуют

базис в линейной оболочке значений h (s),

i = 1,

in.

Поэтому

ортогональное

дополнение

Ht&) Q HtAl)

к подпространству Я ;,(|)

в

простран­

стве Я,(£), порождаемое всеми величинами

 

li(s) — Ptlli(s),

г = 1,

. . . .

т ,

совпадает с замкнутой линейной оболочкой величин Л; (s) — Л/ (0 = Лг ( s ) ~ Pt^t (s).

Случайный процесс ц (s) = {ц,- (s)}"\ ^ s ^ t, является

процессом с некоррелированными приращениями и обладает тем свойством, что

0

Я „ (I) = Я 5з (п) ©

я 5, (т,)

(3.2)

при любых

и в этом

смысле является

о б н о в л я ю щ и м

процессом для

| (s) =

{|,- (s)}"1,

(Можно затем перейти к обновляющему процессу с некоррелированными компонентами и упорядоченными структурными типами, как это сде­ лано в примере на стр. 22.)

Простейшим примером марковского в широком смысле процесса может служить многомерный про­

цесс I (t) = {£( (/))!”, t > tQ, удовлетворяющий системе

стохастических дифференциальных, уравнений вида

т

d £ i{t)= 2 Лг/(0 1 /(0 ^ + *1г(0>

i =

in, (3.3)

/ = i

или, в матричной форме,

£$(/)= л(ОКО Л + *1(0.

где A{i) = {Atj(t)} — непрерывная матричная функция

(неслучайная), а л (0 = (Лг (0)Г>

— процесс с не­

коррелированными приращениями.

 

5 3] НЕКОТОРЫЕ ПРИМЕРЫ ОБНОВЛЯЮЩИХ ПРОЦЕССОВ

31

 

Решение системы (3.3) при заданных начальных

значениях |

(/0) =

{£г (to)}",

точнее,

решение интеграль­

ного уравнения

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

I ( 0

=

I (/о) +

J А (s) I (s) ds +

[Г)(/) -

Т) (/о)], t >

t0. (3.30

 

 

 

и

 

 

 

 

 

 

 

 

 

может

быть представлено

в виде

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

= R (t, to)

l (to) +

J R (t,

S) dr] (s),

(3.4)

 

 

 

 

 

 

 

 

t„

 

 

 

 

где

m X m-матричная

функция

R(t, s) — {RiJ (t,

s)},

t ^ s , — так

называемая

р е з о л ь в е н т а — есть

ре­

шение обыкновенного

 

дифференциального уравнения

 

 

 

dR (t,

s)

A(t)R(t, s),

 

t > s,

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с начальным условием матрица). Действительно,

имеем

t

R(t, s) — I-f [ s

R(s,

s) — I (/ — единичная

по

самому определению

A(u)R(u, s)du

И

 

 

 

 

 

 

 

 

 

J R (t,

 

t Г

 

t

 

 

 

 

s) rfi] (s) =

J

/

+ J A (u) R (u,

s) du

dr] (s) =

tti

 

fa

 

S

 

 

 

 

 

= [л (0 Л (to)] +

 

t Г t

 

 

 

 

|

J A (u) R (u,

s) du

dr] (s) =

 

 

 

^0

S

j

t

Г и

-

 

 

 

 

 

 

 

 

 

= h (t) — Л Vo)] +

A (u)

J R(u,

s) dr\ (s) du.

Видно,

что

 

 

 

 

 

^0

 

 

 

 

 

 

 

 

 

 

 

!(/)

= {

R(t,

s) dr](S),

t ^ t 0,

 

32 ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ [ГЛ. I

есть решение интегрального уравнения (3.3') с на­

чальным

значением

|(/0) = О;

разность

же

x(t) =

= 1 ( 0 l{t),

t > t 0,

как решение однородного

диф­

ференциального

уравнения

dx (t)/dt = А (/) х {t)

с на­

чальным условием х (/0) = 1

(0 ) имеет вид

 

 

 

 

 

 

x(t) =

R(t, tQ)l(t0).

 

 

 

 

Предположим, что приращения г|г (t)—г)г (s),

s,

/ ^ / 0;

г == 1 , .. ., in,

н е к о р р е л и р о в а н ы с начальными

значениями 1У(/0), / =

 

1, . ..,

т .

Тогда, если переопреде­

лить процесс 1] (/) = {т|,- (О}',", t <

t0, точнее, рассматривая

вместо первоначального

процесса г)(0 ,

I <

to,

процесс

с некоррелированными

приращениями

г|(0 — il(0 ) +

-фКО), можно (не меняя обозначений) считать

 

 

Л; (0) =

£; (0).

 

 

 

 

 

 

 

(3.5)

и переписать

выражение (3.4) в

 

виде

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

1

(0

=

J

R(t, s)dp(s).

 

 

 

(3.6)

 

 

 

 

 

*0

 

 

 

 

 

 

 

 

 

Очевидно,

подпространства

Н ,{г|)

порождаются

значениями

1 ,-(0 )

и

приращениями

rp (s)— т)г (tQ),

/ 0 <

Используя тот факт,

что резольвента R{t, s)

удовлетворяет условию

 

 

h, t) ■ R (t, s)

 

 

 

 

R(t +

h, s) =

R{t +

 

 

 

при всех

t ^ s и /г^О ,

из

формулы (3.6)

получаем,

что

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R[t + h, t) l it) =

J Rit +

h, t) R it,

s) dri is) =

 

 

 

 

 

 

 

 

 

 

 

=

J R{t +

h, s) dц is)

H t-\r h) —

R it +

h,

t) l it) =

J

R it +

h,

s) dr\ (s).

Видно, что разности

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h it +

 

m

R ti it +

 

t) 1, it),

 

1,

. .. ,

 

h) -

2

h,

i =

m,

 

 

/=1

 

 

 

 

 

 

 

 

 

 

 

 

§ 3] НЕКОТОРЫЕ "ПРИМЕРЫ ОБНОВЛЯЮЩИХ ПРОЦЕССОВ

33

ортогональны

подпространству

//,(£)>

т. е. величины

т

 

 

 

t

т

 

 

 

У) R i i

(/'+

h, i) ls (t) =

J

У) Rit (t +

h, s) dx\j (s)

M

 

 

 

<=■ /=■

 

h), i =

1 ,. . m,

являются проекциями значении h (l +

на подпространства

Ht(£) s H, (ц).

 

 

Таким

образом,

процесс

| (/) — {h (/)}"'» t^ to ,

является м а р к о в с к и м

в широком

смысле, а фор­

мула (3.6)

дает нам

его

к а н о н и ч е с к о е

предста­

вление с помощью описанного выше (см. (3.5)) про­

цесса Л (0 =

{11г (0 )| . t > t 0,

с

некоррелированными

приращениями.

 

 

 

 

 

 

Очевидно,

Я ,(т])е

Я ,(|)

(см. (3.3')),

т. е. этот про.

цесс является

о б н о в л я ю щ и м

для

%(/.) — {£,• (t)}™,

t ^ t 0, а именно,

 

 

 

 

 

 

при всех

 

 

я ,Ш = я,(л)

 

 

 

t ^ t 0.

 

 

 

 

 

 

2. Стационарные процессы. Рассмотрим стацио­

нарный

в

широком

смысле

процесс

| (() =

(t)}™,

— оо < / <

оо,

т. е. такой,

что корреляционные функ­

ции Bkj (t, s) =

Egfc (t)

(s)

зависят

лишь от разности

t — s:

 

 

 

 

 

k ,j = 1.........m.

 

Bk,(t, s) = Bkl(t — s),

 

Стационарный процесс|(/) =

{£<(*)}",

— °° < t <

oo,

называют регулярным, если

 

 

 

 

 

 

 

П Я Д |) =

0.

 

 

(3.7)

 

 

 

 

t

 

 

 

 

 

Как известно :i!), регулярный процесс имеет спек­ тральную плотность /(Я) = (/Й/(Я)); это положительно определенная, матричная функция от А,, — оо < Я < оо, компоненты /^(Я) которой связаны .с корреляцион­ ными функциями соотношением

 

 

оо

 

 

 

Bkj{t) =

J eaifkj(X)dlj, k, j =

Г, . . . . m.

*)

См'., например,

Ю. А. Р о з а н о в ,

Стационарные слу­

чайные

процессы,

АЛ.,

Физматгиз, 1963. -

 

2 Ю. А. Розанов

34 ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ [ГЛ. 1

Более того, условие (3.7) равносильно тому, что спектральная плотность /(Я) почти всюду имеет один

и тот же ранг п ^ т

и представима

в виде

 

f М =

Ф М • Ф М * п.

в.,

(3.8)

где т X «-матричная функция ф(Я) = {ф*/ (Я)} удовле­ творяет условиям

оо

 

 

I II ф (Я) |\2dX <

оо,

 

“Г

 

(3.9)

| еш ф (Я) dX =

0

при t < О,

ОО

аф (Я)‘ обозначает матрицу, сопряженную к ф (Я).

Обозначим Н2 класс всех матричных функций Ф(Я), удовлетворяющих условиям (3.9). Преобразо­ вание Фурье

со

C{t) = ^ \

elMy (l)d l

(3.10)

—00

 

 

функций из Н2 обращается

в 0 при / < 0

и

оо

Ф (Я) = [ е~ш с (t) dt o'

аналитически продолжается в нижнюю полуплоскость 1 т Я < 0 комплексного переменного Я.

Отправляясь от спектрального представления

оо

 

 

 

U (t )= | е ш ЙФИЯ),

/2 =

1 , . . . , т ,

(3.11)

— ОО

 

 

 

где Ф(Я) = {ФА(Я)}™, — о о < Я < о о ,

— процесс

с не­

коррелированными приращениями, такой, что

 

р

 

 

 

Е Ф *(ц )Ф /(ц )= |/*/(Я )й Я ,

k , j = l , . . . , m ,

§ 3] НЕКОТОРЫЕ ПРИМЕРЫ ОБНОВЛЯЮЩИХ ПРОЦЕССОВ

35

можно получить следующее выражение для рассма­ триваемого стационарного процесса:

 

 

t

п

 

Sfc(0

=

J

/=1

(3.12)

 

 

—00

 

где

т

X л-матричная функция с (f) == {cft} (t)},

О,

есть преобразование Фурье матричной функции q?(A) = = {фАу(А)} класса Я 2, удовлетворяющей условию (3.8), а сам стохастический интеграл в правой части вы­ ражения (3.12) отвечает н е к о р р е л и р о в а н н ы м приращениям

Ч у (*)-Л /(5 )= { ^ = ^ - ф ( Я Г '^ Ф ( Я ) , (3.13)

оо

/= 1 , . . . , п.

Уточним здесь, что,

во-первых,

<р(А,)- 1 есть

п У т -

матричная функция,

о б р а т н а я

к ф(А.) (т.

е. про­

изведение ф(А,)- 1 cp (X)

есть п X «-единичная матрица),

и, во-вторых, для непересекающихся интервалов (s,, г1,)

и (s2,to) величины Л /(0)— Л /^ )

и Л/ (0) — Л/ (5г) не"

коррелированы, причем

 

Е| Л/(О Л/(5) \2== t s,

/ = 1 , . . . , « , (3.14)

а величины лй (А) — Л* (si) и Л/(0) — Л /Ы некоррелированы, для л ю б ы х интегралов, если k ф /. Можно рассматривать выражение (3.12) как неупреждающее линейное преобразование /г-мерного обобщенного про­

цесса л(/) = (л/ (0 }" типа ^белого шума», переписав

(3.12) в виде

-

1 ( 0 = J с (t s) л (s) ds.

(3.15)

 

— оо

 

Обозначим Ht{л) подпространство, порожденное всеми величинами Л /(0— Л/ С5)» / = 1 , . .. , / г, и Q, — оператор проектирования на Я,(т)). Поскольку

Я Д л ) э Я <(|), — o o < t < o o ,

2*

36- . ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ [ГЛ. I.

имеют место следующиенеравенства:

для--любого

!г > О

 

 

 

 

h

т

п

т

 

 

J

- S

S

1СЫ (s)-р d s = = S

Е 111* + /г) -

(* + Л) I2 <

о

f c = l = 1у

= 1

й

 

т

< 2=1 е | ы * + а) - Л £ * (* + а)Р,

а

где Р,, как и раньше, обозначает оператор проекти­ рования на подпространство #,(£)• Если допустить, что среди матричных функций ср (А.) класса Н2, удо­ влетворяющих условию (3.8), найдется функция ср°(А), для которой соответствующие приращения r)°(f) — r|°(s),

/ — 1 , . . . , п, дают равенство

 

 

 

 

 

=

 

 

 

 

(3.16)

 

 

 

 

 

 

 

оо

 

 

 

то

матричная

функция с°(*) =

- ^

|

еш ц>°(А) rfA

по от-

ношению

ко

всем

матричным

— СО

 

 

c(t) —

функциям

 

 

СО

 

 

 

 

 

 

 

 

=

2^- J еш ф(А)с!А

будет м а к с и м а л ь н о й

в том

 

 

— СО

что

 

 

 

 

 

 

 

смысле,

 

 

 

 

 

 

 

h

т

п

 

 

т

+/i)-w*+A)p>

J S 2K/wPds=2 E|^

о

k=\

1=1

 

 

k=\

 

h

rn

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

\

S

J l \ckl(s)?d s

 

 

 

 

 

 

 

0

k=] /=1

 

при всех

/г^О . Такая функция ф °(А )еЯ 2 действи­

тельно существует и определяется Teivx условием, что она м а к с и м а л ь н а в прямом смысле этого слова:

п

т

п

т

11ф°(А)||2= 2

2

|ф^ ( а) | * > 2

2 1 ф й/ М1 2 =||ф(А)|р

k—\ j= l

k=\ f= l

(3.17)

при всех комплексных А, 1 т А < 0 .

Фактически мы уже описали обновляющий про­ цесс, задав соответствующие приращения n]°(i) — ri°(s),

'§ 31 НЕКОТОРЫЕ ПРИМЕРЫ ОБНОВЛЯЮЩИХ ПРОЦЕССОВ

37

j — 1, т . Как показывает общая формула (3.14), кратность М обновляющего процесса равна п — рангу спектральной плотности /(Я):

М— п,

ав упорядоченной системе структурных типов

dF[(t)i> . . . > dFM(t) (см.

(2.4))

все dFt (t)

эквива­

лентны лебеговской мере

 

 

 

d F){t)~dt,

/ = 1.........

М.

(3.18)

Несколько иначе обстоит дело, если стационарный процесс I (/) = {£( {t)}"\ — оо < t < оо,' рассматривается лишь на к о н ё ч н о м интервале t0< t < Т (или при

t > t0).

Поясним это на примере одномерного стационар­ ного процесса £(/) со спектральной плотностью вида

f W = 2n\P(i7.)\2 ’

(З Л 9 )

П

 

где

Р (z) =

2 Pkzk — полином степени

п, все корни

которого

к

в

левой полуплоскости

R e z < 0

лежат

(в частности, р0 Ф 0).

В

этом случае

максимальной

функцией, удовлетворяющей условию (3.8),

будет

 

 

 

ф М — рцх)

 

 

а формула

(3.13)

дает

 

 

 

 

 

’ еш —eiXs

Р (г'Я) dO (Я) =

 

 

й (t)

1l(s) =

i%

 

 

„ш

a

p0 dO (Я) +

 

n—I

+| ( е ш ^ е ^ ) 5 ] р А+,(/Я)^Ф (Я) =

—oo

&=0

1

n—1

:P o| Uu)du +

% PH.l W *>(t)-Zw (s)l- (3-20)

k = 0

38

ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ

[ГЛ. I

Если это соотношение переписать в дифференциаль­ ной форме

dr\ (t) = S Pkl{k) (0 dt + pn

(0,

A=0

 

то легко заметить, что многомерный процесс ( it (t))nir t^ s t0, с компонентами

=k = l .........n,

удовлетворяет системе стохастических дифференциаль­ ных уравнений типа (3.3), а именно,

dU {t) =

lk+i(t)dt,

k = \ .........п 1 ,

 

 

 

 

(3.21)

»

= ("ЙГ £

w

dt +

— dr\ (/).

 

Pn

'/;=1

Кроме того, как и в случае общего стационарного процесса, определяемые выражением (3.20) величины

т](/) — г)(^о). t ^ t 0,

ортогональны к замкнутой

линей­

ной оболочке всех значений £(s),

и,

в част­

ности, к величинам

 

 

 

 

 

lk(t0) =

 

lik- l)(t0),

k = \ .........п.

(3.22)

Таким образом, многомерный процесс (tfc(0}">

принадлежит к тому типу,

что был уже рас­

смотрен нами, и

 

из

полученных

выше результатов

(см. стр. 33) следует,

что в рассматриваемом случае

о б н о в л я ю щ и м

для

[£й(/)}р

будет процесс

(Л* (*)}", t > t 0, с компонентами

 

 

T)k(t) =

U ( t o ) ,

k =

l .........п — 1,

 

Л» (0 =

Л (0 — Л (*о) +

In (to)-

 

Очевидно, этот же процесс будет

о б н о в л я ю щ и м

и для исходного

 

стационарного процесса |(/),

t^to-

Соседние файлы в папке книги из ГПНТБ