VII

Стационарные процессы. Дискретное и непрерывное время

[32/0%]
LaTeX
Задача VII.1

Покажите, что пуассоновская случайная мера (см. §9 главы I), заданная на полукольце ограниченных борелевских множеств в Rn\mathbb {R}^{n}, является ортогональной (нецентрированной) случайной мерой.

?
Задача VII.2

Пусть {Zλ,λ∈R}\left\{ Z_{\lambda }, \lambda \in \mathbb {R}\right\} — комплекснозначный L2L^{2}-процесс такой, что

  1. E[∣Zλ−Zν∣2]→0\mathbb {E}\left[\left|Z_{\lambda }-Z_{\nu }\right|^{2}\right] \rightarrow 0 для любого ν∈R\nu \in \mathbb {R} при λ↓ν\lambda \downarrow \nu,

  2. приращения ортогональны, т.е. для всех λ1<λ2<λ3\lambda_{1} < \lambda_{2} < \lambda_{3}

    E[(Zλ2−Zλ1)(Zλ3−Zλ2)‾]=0 \mathbb {E}\left[\left(Z_{\lambda _{2}}-Z_{\lambda _{1}}\right) \overline{\left(Z_{\lambda _{3}}-Z_{\lambda _{2}}\right)}\right] = 0

На полукольце K={(a,b],−∞<a⩽b<∞}\mathscr {K} = \left\{ (a, b],-\infty < a \leqslant b < \infty \right\}, где (a,a]=∅(a, a] = \varnothing, введем семейство случайных величин Z((a,b]):=Z(b)−Z(a)Z((a, b]): = Z(b)-Z(a). Докажите, что ZZ ортогональная случайная мера на K\mathscr {K}.

Если мера порождается указанным выше процессом Zλ,λ∈RZ_{\lambda }, \lambda \in \mathbb {R}, то вместо интеграла ∫Rf(λ)Z(dλ)\int_{\mathbb {R}} f(\lambda ) Z(d \lambda ) используют также запись ∫Rf(λ)dZλ\int_{\mathbb {R}} f(\lambda ) d Z_{\lambda }.

?
Задача VII.3

Пусть {Zλ,λ∈R}\left\{ Z_{\lambda }, \lambda \in \mathbb {R}\right\} — центрированный гауссовский процесс с ортогональными приращениями, непрерывный справа в среднем квадратическом (выясните, обязан ли тогда процесс быть непрерывным слева в среднем квадратическом). Пусть при каждом t∈Tt \in T функция g(t,⋅):R→Cg(t, \cdot ): \mathbb {R} \rightarrow \mathbb {C} удовлетворяет условию

∫R∣g(t,λ)∣2μ(dλ)<∞ \int _{\mathbb {R}}\left|g(t, \lambda )\right|^{2} \mu (d \lambda ) < \infty

где μ\mu — структурная мера, отвечающая ортогональной случайной мере ZZ, порожденной процессом {Zλ,λ∈R}\left\{ Z_{\lambda }, \lambda \in \mathbb {R}\right\} (см. упражнение 2). Докажите, что тогда процесс

Y={Y(t)=∫Rg(t,λ)dZλ,t∈T} Y = \left\{ Y(t) = \int _{\mathbb {R}} g(t, \lambda ) d Z_{\lambda }, t \in T\right\}

является центрированным гауссовским процессом.

?
Задача VII.4

(замена меры в интеграле) Пусть ZZ — ортогональная случайная мера, заданная на Ω×G\Omega \times \mathscr {G} и имеющая σ\sigma-конечную структурную функцию (меру) μ\mu на измеримом пространстве (Λ,A)(\Lambda , \mathscr {A}) (система множеств G={A∈A:μ(A)<∞})\mathscr {G} = \left\{ A \in \mathscr {A}: \mu (A) < \infty \right\} ). Пусть h:Λ→C,h∈L2(Λ,A,μ)h: \Lambda \rightarrow \mathbb {C}, h \in L^{2}(\Lambda , \mathscr {A}, \mu ). Введем

H={B∈A:∫B∣h(λ)∣2μ(dλ)<∞} \mathscr {H} = \left\{ B \in \mathscr {A}: \int _{B}\left|h(\lambda )\right|^{2} \mu (d \lambda ) < \infty \right\}

и положим

V(B)=∫Λ1B(λ)h(λ)Z(dλ),B∈H V(B) = \int _{\Lambda } \mathbf{1}_{B}(\lambda ) h(\lambda ) Z(d \lambda ), \quad B \in \mathscr {H}

Объясните, почему VV является ортогональной случайной мерой (центрированной, если центрирована мера ZZ ) с σ\sigma-конечной структурной мерой

ν(B)=∫B∣h(λ)∣2μ(dλ)<∞,B∈H \nu (B) = \int _{B}\left|h(\lambda )\right|^{2} \mu (d \lambda ) < \infty , \quad B \in \mathscr {H}

Докажите, что если g:Λ→C,g∈L2(Λ,σ{H},ν)g: \Lambda \rightarrow \mathbb {C}, g \in L^{2}(\Lambda , \sigma \left\{ \mathscr {H}\right\} , \nu ), то

∫Λg(λ)V(dλ)=∫Λg(λ)h(λ)Z(dλ) \int _{\Lambda } g(\lambda ) V(d \lambda ) = \int _{\Lambda } g(\lambda ) h(\lambda ) Z(d \lambda )

В связи с теоремой Карунена (см. также [38,58,134][38,58,134] ) вновь обратимся к вопросу о канонических представлениях случайных процессов.

Рассмотрим пространство L2[a,b]L^{2}[a, b], состоящее из комплекснозначных функций, квадрат модуля которых интегрируем по Лебегу на [a,b][a, b]. Скалярное произведение в этом гильбертовом пространстве определим, как обычно, формулой

(f,g)=∫abf(t)g(t)‾dt,f,g∈L2[a,b] (f, g) = \int _{a}^{b} f(t) \overline{g(t)} d t, \quad f, g \in L^{2}[a, b]

Следующий пример является хорошей иллюстрацией теоремы 3.

Пример 2. Для винеровского процесса WW на отрезке [0,2π][0,2 \pi ] справедливо представление (Карунена)

W(t)=12π∑k=−∞∞1−eiktikzk,t∈[0,2π] W(t) = \frac{1}{\sqrt{2 \pi }} \sum _{k = -\infty }^{\infty } \frac{1-e^{i k t}}{i k} z_{k}, \quad t \in [0,2 \pi ]

где zkz_{k} — ортонормированные центрированные случайные величины, ряд сходится при каждом t∈[0,2π]t \in [0,2 \pi ] в среднем квадратическом (при k=0k = 0 считаем (1−eikt)/ik=−t)\left(1-e^{i k t}\right) / i k = -t).

Действительно, для s,u∈[0,2π]s, u \in [0,2 \pi ] имеем

1[0,s](u)=12π∑k=−∞∞ck(s)eiku \mathbf{1}_{[0, s]}(u) = \frac{1}{\sqrt{2 \pi }} \sum _{k = -\infty }^{\infty } c_{k}(s) e^{i k u}

Функции ( 1/2π)eiku,k∈Z1 / \sqrt{2 \pi }) e^{i k u}, k \in \mathbb {Z}, образуют полную ортонормированную систему в пространстве L2[0,2π]L^{2}[0,2 \pi ]; при каждом s∈[0,2π]s \in [0,2 \pi ] ряд (69) сходится в этом пространстве и коэффициенты фф урье

ck(s)=12π∫02π1[0,s](u)e−ikudu=12π(1−e−iks)/ik,k∈Z c_{k}(s) = \frac{1}{\sqrt{2 \pi }} \int _{0}^{2 \pi } \mathbf{1}_{[0, s]}(u) e^{-i k u} d u = \frac{1}{\sqrt{2 \pi }}\left(1-e^{-i k s}\right) / i k, \quad k \in \mathbb {Z}

(c0(s)=s/2π)\left(c_{0}(s) = s / \sqrt{2 \pi }\right). Равенство Парсеваля показывает, что при s,t∈[0,2π]s, t \in [0,2 \pi ]

min⁡{s,t}=⟨1[0,s],1[0,t]⟩=12π∑k=−∞∞(1−e−iks)(1−eikt)k2 \min \left\{ s, t\right\} = \left\langle \mathbf{1}_{[0, s]}, \mathbf{1}_{[0, t]} \right\rangle = \frac{1}{2 \pi } \sum _{k = -\infty }^{\infty } \frac{\left(1-e^{-i k s}\right)\left(1-e^{i k t}\right)}{k^{2}}

Представление (68) следует из теоремы 3 , в которой надо положить Λ=Z\Lambda = \mathbb {Z},

f(t,k)=12π⋅1−e−iktik,t∈[0,2π],k∈Λ f(t, k) = \frac{1}{\sqrt{2 \pi }} \cdot \frac{1-e^{-i k t}}{i k}, \quad t \in [0,2 \pi ], \quad k \in \Lambda

и выбрать меру μ\mu считающей мерой на Λ\Lambda, т.е. такой, что μ({k})=1,k∈Z\mu (\left\{ k\right\} ) = 1, k \in \mathbb {Z}. □\square

?
Задача VII.5

Докажите, что представление (68) справедливо (п.н.) с независимыми случайными величинами zk∼N(0,1),k∈Zz_{k} \sim \mathrm{N}(0,1), k \in \mathbb {Z}. Точнее говоря, докажите, что при таком выборе величин zk(k∈Zz_{k}(k \in \mathbb {Z} ) правая часть (68) имеет непрерывную модификацию, являющуюся центрированным гауссовским процессом с ковариационной функцией min⁡{s,t},s,t∈[0,2π]\min \left\{ s, t\right\} , s, t \in [0,2 \pi ].

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

C\mathrm{C} непрерывной на [a,b]×[a,b][a, b] \times [a, b] ковариационной функцией r=r(s,t)r = r(s, t) свяжем ограниченный линейный оператор Фредгольма AA, который действует в L2[a,b]L^{2}[a, b] и определяется следующим образом:

Af(s)=∫abr(s,t)f(t)dt,f∈L2[a,b] A f(s) = \int _{a}^{b} r(s, t) f(t) d t, \quad f \in L^{2}[a, b]

Функция r=r(s,t)r = r(s, t) называется ядром интегрального оператора AA.

?
Задача VII.6

Пусть X={X(t),t∈[a,b]}X = \left\{ X(t), t \in [a, b]\right\} — комплекснозначный L2L^{2}-процесс, непрерывньй в среднем квадратическом на отрезке [a,b][a, b]. Докажите, что это свойство равносильно тому, что на [a,b][a, b] непрерывна функция m(t)=E[X(t)]m(t) = \mathbb {E}\left[X(t)\right] и на [a,b]×[a,b][a, b] \times [a, b] непрерывна ковариационная функция r(s,t)=cov⁡(X(s),X(t))r(s, t) = \operatorname {cov}(X(s), X(t)).

?
Задача VII.7

Докажите, что оператор AA, заданный формулой (70), компактен, т. е. всякое ограниченное по норме множество оператор AA переводит в предкомпактное (т.е. множество, замыкание которого есть компакт).

Введенньй оператор AA — оператор Гильберта-Шмидта (см. [35; гл. IX, § 2, п. 1]), поскольку

∫ab∫ab∣r(s,t)∣2dsdt<∞ \int _{a}^{b} \int _{a}^{b}\left|r(s, t)\right|^{2} d s d t < \infty

Заметим также, что оператор AA является самосопряженным, поскольку (см. [35; гл. IX, § 2, п. 2]) его сопряженный оператор A∗A^{*} задается эрмитово сопряженным ядром r(t,s)‾\overline{r(t, s)}, а для ковариационной функции r(s,t)r(s, t) справедливо равенство r(s,t)=r(t,s)‾r(s, t) = \overline{r(t, s)} при всех s,ts, t.

Напомним (см. [35[35; гл. 4,§6]4, \S 6] ) классический результат.

Теорема 17 (Гильберт-Шмидт). Пусть AA — компактный самосопря женный оператор в (комплексном) гильбертовом пространстве НН. Тогда существует ортонормированная система {ϕn,n∈J}\left\{ \phi_{n}, n \in J\right\} собственных векторов оператора AA, отвечающих ненулевым собственным значениям, такая, что любой элемент h∈Hh \in H записывается единственным образом в виде

h=∑n∈Jcnϕn+u h = \sum _{n \in J} c_{n} \phi _{n}+u

где cn∈C,n∈J,u∈Ker⁡A(m.e.Au=0)c_{n} \in \mathbb {C}, n \in J, u \in \operatorname {Ker} A(m . e . A u = 0). Множество JJ всегда конечно или счетно (в последнем случае ряд (72) сходится по норме ∥⋅∥=(⋅,⋅)1/2\left\| \cdot \right\| = (\cdot , \cdot )^{1 / 2} ). Более того, в силу компактности оператора AA число собственных значений λn\lambda_{n} во внешности любого круга с центром в 0 конечно и поэтому их можно перенумеровать в порядке невозрастания модулей ∣λ1∣⩾∣λ2∣⩾⋯\left|\lambda_{1}\right| \geqslant \left|\lambda_{2}\right| \geqslant \cdots. Если J счетно, то lim⁡n→∞λn=0\lim_{n \rightarrow \infty } \lambda_{n} = 0.

Отметим, что собственные значения самосопряженного оператора AA в HH действительны, а собственные векторы, отвечающие разным собственным значениям, ортогональны, в частности, (ϕn,u)=0,n∈J\left(\phi_{n}, u\right) = 0, n \in J.

Для оператора AA из (70) в силу сепарабельности H=L2[a,b]H = L^{2}[a, b] в подпространстве Ker AA можно выбрать ортонормированный базис {ϕk}k∈M\left\{ \phi_{k}\right\}_{k \in M}, где MM — конечно или счетно (считаем множества MM и JJ непересекающимися; если Ker⁡A=0\operatorname {Ker} A = 0, то M=∅M = \varnothing ). Дополнив систему {ϕn,n∈J}\left\{ \phi_{n}, n \in J\right\} до ортонормированного базиса {ϕn,n∈J∪M}\left\{ \phi_{n}, n \in J \cup M\right\} в HH, перепишем (72) в виде разложения по собственному базису оператора AA :

h=∑k∈J∪Mckϕk,ck=(h,ϕk),k∈J∪M h = \sum _{k \in J \cup M} c_{k} \phi _{k}, \quad c_{k} = \left(h, \phi _{k}\right), \quad k \in J \cup M

Если семейство индексов J∪MJ \cup M счетно, то ряд Фурье в (73) сходится по норме L2[a,b]L^{2}[a, b].

В силу упражнений 6 и 7 для ковариационной функции rr непрерывного в среднем квадратическом процесса X={X(t),t∈[a,b]}X = \left\{ X(t), t \in [a, b]\right\} имеем согласно (73), что при каждом t∈[a,b]t \in [a, b]

r(t,⋅)‾=∑k∈J∪Mck(t)ϕk(⋅) \overline{r(t, \cdot )} = \sum _{k \in J \cup M} c_{k}(t) \phi _{k}(\cdot )

где

ck(t)=∫abr(t,s)‾ϕk(s)‾ds=λkϕk(t)‾=λkϕk(t)‾ c_{k}(t) = \int _{a}^{b} \overline{r(t, s)} \overline{\phi _{k}(s)} d s = \overline{\lambda _{k} \phi _{k}(t)} = \lambda _{k} \overline{\phi _{k}(t)}

Тем самым,

r(t,⋅)=∑k∈Jλkϕk(t)ϕk(⋅)‾ r(t, \cdot ) = \sum _{k \in J} \lambda _{k} \phi _{k}(t) \overline{\phi _{k}(\cdot )}

Если JJ — счетно, то ряды (73), (74) сходятся в L2[a,b]L^{2}[a, b] при любом t∈[a,b]t \in [a, b]. Отметим, что здесь все λk\lambda_{k} действительны и λk=0,k∈M\lambda_{k} = 0, k \in M.

?
Задача VII.8

Докажите, что (в силу неотрицательной определенности функции rr ) λk>0\lambda_{k} > 0 при k∈Jk \in J.

Следующая теорема показывает, что соотношение (74) допускает усиление (см., например, [60; с. 263]).

Теорема 18 (Мерсер). Пусть r=r(s,t)r = r(s, t) — непрерывная на [a,b]×[a,b]κова[a, b] \times [a, b] \kappa о в а риационная функция. Тогда для всех s,t∈[a,b]s, t \in [a, b]

r(s,t)=∑k∈Jλkϕk(s)ϕk(t)‾ r(s, t) = \sum _{k \in J} \lambda _{k} \phi _{k}(s) \overline{\phi _{k}(t)}

причем, если JJ — счетно, то на [a,b]×[a,b][a, b] \times [a, b] ряд (75) сходится абсолютно и равномерно. (По-прежнему, λkuϕk\lambda_{k} u \phi_{k} обозначают положительные собственные значения и соответствующие собственные функции оператора AA из (70).)

В качестве следствия отсюда вытекает

Теорема 19 (Карунен-Лоэв). Пусть X={X(t),t∈[a,b]}X = \left\{ X(t), t \in [a, b]\right\} — непрерывный в среднем квадратическом на [a,b][a, b] центрированный L2L^{2}-процесс, имеющий ковариационную функцию r=r(s,t)r = r(s, t). Тогда при каждом t∈[a,b]t \in [a, b]

X(t)=∑k∈Jλkϕk(t)zk n. н.,  X(t) = \sum _{k \in J} \sqrt{\lambda _{k}} \phi _{k}(t) z_{k} \quad \text{ n. н., }

где λk,ϕk\lambda_{k}, \phi_{k} — те же, что и в (75), а {zk,k∈J}\left\{ z_{k}, k \in J\right\} — центрированные ортонормированные случайные величины. Для счетного J ряд (76) сходится в среднем квадратическом.

Доказательство. Положим T=[a,b],Λ=JT = [a, b], \Lambda = J и f(t,k)=λkϕk(t)f(t, k) = \sqrt{\lambda_{k}} \phi_{k}(t) для t∈Tt \in T, k∈Λk \in \Lambda. Тогда свойство (75) превращается в соотношение (18), где μ\mu — считающая мера на Λ\Lambda, т.е. μ({k})=1\mu (\left\{ k\right\} ) = 1 для k∈Λk \in \Lambda и требуемый результат о справедливости (76) непосредственно вытекает из теоремы Карунена (теорема 3 ).

?
Задача VII.9

Оценив величину

E[∣X(t)−∑k⩽nλkϕk(t)zk∣2] \mathbb {E}\left[\left|X(t)-\sum _{k \leqslant n} \sqrt{\lambda _{k}} \phi _{k}(t) z_{k}\right|^{2}\right]

докажите, что если выполнены условия теоремы 19 , то в разложении (76) в качестве величин zkz_{k} можно выбрать центрированные ортонормированные величины

zk=1λk∫abX(t)ϕk(t)dt z_{k} = \frac{1}{\sqrt{\lambda _{k}}} \int _{a}^{b} X(t) \phi _{k}(t) d t

Интеграл здесь определяется (и существует) как предел в среднем квадратическом соответствующих интегральных сумм Римана. Объясните, почему для непрерывной на [a,b]×[a,b][a, b] \times [a, b] ковариационной функции r=r(s,t)r = r(s, t) собственные функции ϕk∈C[a,b],k∈J\phi_{k} \in C[a, b], k \in J. Кроме того, докажите, что среднеквадратическая сходимость в (76) с данными zk,k∈Jz_{k}, k \in J, будет равномерна по t∈[a,b]t \in [a, b].

Пример 3. Найдем, пользуясь теоремой 19, разложение Карунена-Лоэва для винеровского процесса W(t)W(t) на отрезке [0,1][0,1].

Для этого определим собственные значения и собственные функции интегрального оператора, отвечающего непрерывной на [0,1]×[0,1][0,1] \times [0,1] ковариационной функции r(s,t)=min⁡{s,t},s,t∈[0,1]r(s, t) = \min \left\{ s, t\right\} , s, t \in [0,1]. Согласно упражнению 8 величины λk>0\lambda_{k} > 0 для k∈Jk \in J и, как показывает упражнение 9,ϕk∈C[0,1]9, \phi_{k} \in C[0,1].

Рассмотрим относительно λ>0\lambda > 0 и ϕ∈C[0,1]\phi \in C[0,1] уравнение

∫01r(s,t)ϕ(t)dt=λϕ(s),s∈[0,1] \int _{0}^{1} r(s, t) \phi (t) d t = \lambda \phi (s), \quad s \in [0,1]

т. е. уравнение

∫0stϕ(t)dt+s∫s1ϕ(t)dt=λϕ(s) \int _{0}^{s} t \phi (t) d t+s \int _{s}^{1} \phi (t) d t = \lambda \phi (s)

Поскольку левая часть (78) дифференцируема в каждой точке s∈[0,1]s \in [0,1] (в концах отрезка надо рассматривать соответствующие односторонние производные), то функция ϕ\phi дифференцируема и

∫s1ϕ(t)dt=λϕ′(s),s∈[0,1] \int _{s}^{1} \phi (t) d t = \lambda \phi ^{\prime }(s), \quad s \in [0,1]

Снова дифференцируя по ss обе части этого равенства, приходим к тому, что требуемая функция ϕ(s)\phi (s) должна удовлетворять следующему уравнению:

ϕ′′(s)=−(1/λ)ϕ(s) \phi ^{\prime \prime }(s) = -(1 / \lambda ) \phi (s)

Обшее решение этого уравнения задается формулой

ϕ(s)=Acos⁡(s/λ)+Bsin⁡(s/λ) \phi (s) = A \cos (s / \sqrt{\lambda })+B \sin (s / \sqrt{\lambda })

где A,BA, B — константы. Учитывая, что ϕ(0)=0,ϕ′(1)=0\phi (0) = 0, \phi^{\prime }(1) = 0, находим, что

A=0,(B/λ)cos⁡(1/λ)=0 A = 0, \quad (B / \sqrt{\lambda }) \cos (1 / \sqrt{\lambda }) = 0

откуда

λk=((k+1/2)π)−2,k∈Z+={0,1,…} \lambda _{k} = ((k+1 / 2) \pi )^{-2}, \quad k \in \mathbb {Z}_{+} = \left\{ 0,1, \ldots \right\}

Условия нормировки ∫01ϕk2(s)ds=1\int_{0}^{1} \phi_{k}^{2}(s) d s = 1 и равенства (79) дают следующие представления для функций ϕk(s)\phi_{k}(s) :

ϕk(s)=2sin⁡((k+1/2)πs),s∈[0,1],k∈Z+ \phi _{k}(s) = \sqrt{2} \sin ((k+1 / 2) \pi s), \quad s \in [0,1], \quad k \in \mathbb {Z}_{+}

Легко проверить, что все найденные величины λk\lambda_{k} и ϕk,k∈Z+\phi_{k}, k \in \mathbb {Z}_{+}, удовлетворяют исходному уравнению (78). Таким образом, согласно представлению (76) находим, что

W(t)=2∑k=0∞sin⁡((k+1/2)πt)(k+1/2)πtzk W(t) = \sqrt{2} \sum _{k = 0}^{\infty } \frac{\sin ((k+1 / 2) \pi t)}{(k+1 / 2) \pi t} z_{k}

где {zk,k∈Z+}\left\{ z_{k}, k \in \mathbb {Z}_{+}\right\}- последовательность ортонормированных и центрированных (в силу упражнения 9) случайных величин; ряд в (81) сходится в среднем квадратическом равномерно на [0,1][0,1].

Замечание 6. Упражнение 4 главы II показывает, что L2[W(t),t∈[0,1]]L^{2}[W(t), t \in [0,1]] — гауссовская система. Согласно упражнению 9 величины {zk,k∈Z+}\left\{ z_{k}, k \in \mathbb {Z}_{+}\right\}, фигурирующие в (77), где X(t)=W(t),a=0X(t) = W(t), a = 0 и b=1b = 1, образуют гауссовский процесс. Поскольку для центрированных действительных гауссовских величин ортогональность равносильна независимости, то, применив упражнение 11 главы IV, получаем, что ряд (81) при каждом t∈[0,1]t \in [0,1] сходится с вероятностью единица, где независимые величины zk∼N(0,1),k∈Z+z_{k} \sim \mathrm{N}(0,1), k \in \mathbb {Z}_{+}, определены на том же самом вероятностном пространстве (Ω,F,P)(\Omega , \mathscr {F}, \mathbb {P}), что и данный винеровский процесс.

?
Задача VII.10

Найдите на отрезке [0,c][0, c] разложение Карунена-Лоэва для процесса Орнштейна-Уленбека, имеющего ковариационную функцию r(s,t)=e−α∣s−t∣r(s, t) = e^{-\alpha \left|s-t\right|}, где α>0,s,t∈R\alpha > 0, s, t \in \mathbb {R}.

?
Задача VII.11

Пусть X={X(t),t∈R}X = \left\{ X(t), t \in \mathbb {R}\right\} — непрерывный в среднем квадратическом процесс с ковариационной функцией r(s,t)=R(s−t)r(s, t) = R(s-t). Пусть XX — периодический процесс с периодом aa, т.е. X(t+a)=X(t)X(t+a) = X(t) п. н. для каждого t∈Rt \in \mathbb {R}. Найдите для XX разложение Карунена-Лоэва.

Применения различных канонических разложений случайных функций даны, например, в монографии [58]. Важное значение канонические представления имеют для приложений теории случайных процессов к решению задач обнаружения сигналов в присутствии помех. В частности, большой интерес представляет вопрос, в каком случае в той или иной задаче можно ограничиться конечным числом членов, участвующих в каноническом разложении случайной функции. Иллюстрирующие примеры такого рода можно найти в [58;§174][58 ; \S 174].

Между теоремой 3 главы VII и теоремой 6 главы II имеется тесная связь.

?
Задача VII.12

Докажите, что если ковариационная функция r=r(s,t)r = r(s, t) допускает представление (18), то гильбертово пространство HH с воспроизводящим ядром r=r(s,t)r = r(s, t) состоит из функций

h(t)=∫Λg(λ)f(t,λ)‾μ(dλ),g∈L2[f],t∈T h(t) = \int _{\Lambda } g(\lambda ) \overline{f(t, \lambda )} \mu (d \lambda ), \quad g \in L^{2}[f], \quad t \in T

где L2[f]L^{2}[f] — замыкание в L2(Λ,A,μ)L^{2}(\Lambda , \mathscr {A}, \mu ) линейной оболочки функций f(t,⋅),t∈Tf(t, \cdot ), t \in T. Скалярное произведение в HH задается формулой

⟨h1,h2⟩=∫Λg1(λ)g2(λ)‾μ(dλ) \left\langle h_{1}, h_{2} \right\rangle = \int _{\Lambda } g_{1}(\lambda ) \overline{g_{2}(\lambda )} \mu (d \lambda )

где hkh_{k} и gk,k=1,2g_{k}, k = 1,2, связаны соотношением (82). Пусть ZZ — ортогональная случайная мера на G×Ω\mathscr {G} \times \Omega, имеющая структурной функцией меру μ\mu (см. теорему 2), и пусть L2[f]=L2(Λ,A,μ)L^{2}[f] = L^{2}(\Lambda , \mathscr {A}, \mu ). Объясните, почему для процесса X={X(t),t∈T}X = \left\{ X(t), t \in T\right\} вида (19), процесс YY, фигурирующий в теореме 6 главы II, задается формулой

Y(h)=∫Λg(λ)Z(dλ) Y(h) = \int _{\Lambda } g(\lambda ) Z(d \lambda )

где hh определяется по gg согласно (82).

?
Задача VII.13

Докажите, что (действительное) гильбертово пространство HH с воспроизводящим ядром r(s,t)=min⁡{s,t},s,t∈R+r(s, t) = \min \left\{ s, t\right\} , s, t \in \mathbb {R}_{+}, являющимся ковариационной функцией винеровского процесса, состоит из функций

h(t)=∫0tg(λ)dλ,t∈R+, для которых ∫0∞g2(λ)dλ<∞ h(t) = \int _{0}^{t} g(\lambda ) d \lambda , \quad t \in \mathbb {R}_{+}, \quad \text{ для которых } \int _{0}^{\infty } g^{2}(\lambda ) d \lambda < \infty

(интегрируем по мере Лебега). Объясните, почему L2[W]L^{2}[W], т.е. замкнутая в L2(Ω,F,P)L^{2}(\Omega , \mathscr {F}, \mathbb {P}) линейная оболочка W(t),t∈R+W(t), t \in \mathbb {R}_{+}, состоит из величин вида (83), где ZZ — ортогональная случайная мера на G×Ω\mathscr {G} \times \Omega, порожденная самим винеровским процессом (см. пример 1 ), а функция gg удовлетворяет условию (84).

Рассмотрим ряд конкретных примеров стационарных и нестационарных процессов (подразумевая под стационарностью стационарность в широком смысле).

?
Задача VII.14

Пусть X={X(t)=ξei(ηt+ν),t∈R}X = \left\{ X(t) = \xi e^{i(\eta t+\nu )}, t \in \mathbb {R}\right\}, где действительные случайные величины ξ,η,ν\xi , \eta , \nu таковы, что ν\nu равномерно распределена на [0,2π][0,2 \pi ] и не зависит от (ξ,η)(\xi , \eta ), а E[ξ2]<∞\mathbb {E}\left[\xi^{2}\right] < \infty. Убедитесь, что XX дает пример центрированного стационарного процесса с ковариационной функцией вида (36), где мера G=E[ξ2]FG = \mathbb {E}\left[\xi^{2}\right] F, а FF — распределение вероятностей величины η\eta. Положим теперь Y={Y(t)=ξcos⁡(ηt+ν),t∈R}Y = \left\{ Y(t) = \xi \cos (\eta t+\nu ), t \in \mathbb {R}\right\}. Проверьте, что YY — стационарный процесс. Если X={X(t),t∈R}X = \left\{ X(t), t \in \mathbb {R}\right\} — стационарный комплекснозначный процесс, то можно ли утверждать, что {Re⁡X(t),t∈R}\left\{ \operatorname {Re} X(t), t \in \mathbb {R}\right\} и {Im⁡X(t),t∈R}\left\{ \operatorname {Im} X(t), t \in \mathbb {R}\right\} стационарные процессы?

?
Задача VII.15

(телеграфный сигнал). Пусть N={Nt,t⩾0}N = \left\{ N_{t}, t \geqslant 0\right\} — пуассоновский процесс, имеющий постоянную интенсивность λ>0\lambda > 0. Пусть случайная величина ζ\zeta не зависит от процесса NN и P(ζ=−1)=P(ζ=1)=1/2\mathbb {P}\left(\zeta = -1\right) = \mathbb {P}\left(\zeta = 1\right) = 1 / 2. Введем процесс, называемый телеграфным сигналом или телеграфной волной, положив

X(t,ω)=ζ(ω)(−1)Nt(ω),t⩾0,ω∈Ω X(t, \omega ) = \zeta (\omega )(-1)^{N_{t}(\omega )}, \quad t \geqslant 0, \quad \omega \in \Omega

Найдите cov⁡(X(s),X(t)),s,t⩾0\operatorname {cov}(X(s), X(t)), s, t \geqslant 0 (сопоставьте ответ с упражнением 27 главы III).

?
Задача VII.16

(обобщение упражнения 15). Пусть {ζn,n⩾1}\left\{ \zeta_{n}, n \geqslant 1\right\} — последовательность независимых одинаково распределенных действительных случайных величин, независящая от пуассоновского процесса N={Nt,t⩾0}N = \left\{ N_{t}, t \geqslant 0\right\}, имеющего с вероятностью 1 единичные скачки в моменты τ1<τ2<⋯\tau_{1} < \tau_{2} < \cdots (см. теорему 2 главы II). Положим X(t,ω)=ζk(ω)X(t, \omega ) = \zeta_{k}(\omega ) при τk(ω)⩽t<τk+1(ω)\tau_{k}(\omega ) \leqslant t < \tau_{k+1}(\omega ), где k=0,1,…(τ0=0)k = 0,1, \ldots \left(\tau_{0} = 0\right). Если в последовательности {τk(ω),k⩾0}\left\{ \tau_{k}(\omega ), k \geqslant 0\right\} найдутся совпадающие точки (что возможно лишь с вероятностью 0 ), то пусть X(t,ω)=0X(t, \omega ) = 0 для t⩾0t \geqslant 0. Считая E[ζ02]<∞\mathbb {E}\left[\zeta_{0}^{2}\right] < \infty, найдите ковариационную функцию процесса X={X(t),t⩾0}X = \left\{ X(t), t \geqslant 0\right\}.

?
Задача VII.17

Для процесса дробового шума (см. определение 14 главы VI) с ковариационной функцией r=r(s,t)r = r(s, t) докажите, что возникает асимптотическая стационарность в том смысле, что существует lim⁡s→∞r(s,s+t)\lim_{s \rightarrow \infty } r(s, s+t).

Разнообразные обобщения процесса дробового шума (в том числе поля дробово20 шума) играют важную роль в построении различных моделей более сложных процессов; см. по этому поводу, например, [13,111][13,111].

?
Задача VII.18

Найдите спектральную плотность процесса Орнштейна-Уленбека.

?
Задача VII.19

Приведите пример (в связи с теоремой 8) неотрицательно определенной функции R(t),t∈RR(t), t \in \mathbb {R}, непрерывной всюду, кроме точки t=0t = 0.

?
Задача VII.20

Пусть h1(t),…,hn(t)h_{1}(t), \ldots , h_{n}(t) — комплекснозначные функции, заданные на множестве TT. Докажите, что r(s,t)=∑j=1nhj(s)hj(t)‾,s,t∈Tr(s, t) = \sum_{j = 1}^{n} h_{j}(s) \overline{h_{j}(t)}, s, t \in T, является неотрицательно определенной функцией. Дайте пример негауссовского процесса, имеющего эту ковариационную функцию.

?
Задача VII.21

Пусть {Xn,n∈Z}\left\{ X_{n}, n \in \mathbb {Z}\right\} — стационарный процесс со средним значением aa и ковариационной функцией R=R(n),n∈ZR = R(n), n \in \mathbb {Z}. Докажите, что при N→∞N \rightarrow \infty

1N∑k=0N−1Xk⟶L2(Ω)a⟺1N∑k=0N−1R(k)→0 \frac{1}{N} \sum _{k = 0}^{N-1} X_{k} \stackrel{L^{2}(\Omega )}{\longrightarrow } a \Longleftrightarrow \frac{1}{N} \sum _{k = 0}^{N-1} R(k) \rightarrow 0
?
Задача VII.22

Пусть в условиях предыдущего упражнения

1N2∑k=1N∑m=1NR(k−m)≡1N∑∣k∣⩽N−1R(k)(1−∣k∣/N)⩽cN−γ \frac{1}{N^{2}} \sum _{k = 1}^{N} \sum _{m = 1}^{N} R(k-m) \equiv \frac{1}{N} \sum _{\left|k\right| \leqslant N-1} R(k)(1-\left|k\right| / N) \leqslant c N^{-\gamma }

для некоторых c,γ>0c, \gamma > 0 и всех N∈NN \in \mathbb {N}. Докажите, что тогда

1N∑k=0N−1Xk→0 п.н. при N→∞ \frac{1}{N} \sum _{k = 0}^{N-1} X_{k} \rightarrow 0 \text{ п.н. при } N \rightarrow \infty

Покажите, что условие (85) выполнено, если R(n)=O(n−γ)R(n) = O\left(n^{-\gamma }\right) при n→∞n \rightarrow \infty.

Теорема 20 (Гапошкин; [14]). Пусть X={X(n),n∈Z}X = \left\{ X(n), n \in \mathbb {Z}\right\} — стачионарный процесс со средним a=0a = 0 и ковариационной функцией R=R(n),n∈ZR = R(n), n \in \mathbb {Z} (пусть R(0)=1)R(0) = 1). Тогда (86) имеет место в том и только том случае, когда G({0})=0uG(\left\{ 0\right\} ) = 0 u

lim⁡n→∞∫0<∣λ∣⩽2−nZ(dλ)=0n.μ. \lim _{n \rightarrow \infty } \int _{0 < \left|\lambda \right| \leqslant 2^{-n}} Z(d \lambda ) = 0 \quad n . \mu .

где ZZ — ортогональная случайная мера, фигурирующая в спектральном представлении (35), а GG — ее структурная функция (мера).

Отметим, что доказательство этого результата основано на представлении сумм k−1∑j=1kXjk^{-1} \sum_{j = 1}^{k} X_{j} для 2n⩽k<2n+1(n=0,1,…)2^{n} \leqslant k < 2^{n+1}(n = 0,1, \ldots ) в виде ∫∣λ∣⩽2−nZ(dλ)+ψk\int_{\left|\lambda \right| \leqslant 2^{-n}} Z(d \lambda )+\psi_{k}, где lim⁡k→∞ψk=0\lim_{k \rightarrow \infty } \psi_{k} = 0 п.н. и в L2(Ω)L^{2}(\Omega ).

?
Задача VII.23

Пусть X={Xn,n∈N}X = \left\{ X_{n}, n \in \mathbb {N}\right\} — центрированная действительная гауссовская стационарная последовательность с ковариационной функцией R=R(n),n∈∈ZR = R(n), n \in \in \mathbb {Z}. Докажите, что условие

1N∑k=0N−1R2(k)→0,N→∞ \frac{1}{N} \sum _{k = 0}^{N-1} R^{2}(k) \rightarrow 0, \quad N \rightarrow \infty

необходимо и достаточно для того, чтобы оценка R^N(m)\widehat{R}_{N}(m), введенная в (48), была состоятельной в среднем квадратическом, т.е. для m∈Zm \in \mathbb {Z}

E[∣R^N(m)−R(m)∣2]→0 при N→∞.  \mathbb {E}\left[\left|\widehat{R}_{N}(m)-R(m)\right|^{2}\right] \rightarrow 0 \text{ при } \quad N \rightarrow \infty \text{. }
?
Задача VII.24

Пусть X={Xn,n∈Z}X = \left\{ X_{n}, n \in \mathbb {Z}\right\} — процесс скользящего среднего вида (39). Потребуем, чтобы ∑k=−∞∞∣ck∣<∞\sum_{k = -\infty }^{\infty }\left|c_{k}\right| < \infty и εk,k∈Z\varepsilon_{k}, k \in \mathbb {Z}, были независимыми одинаково распределенными величинами, имеющими E[εk4]<∞\mathbb {E}\left[\varepsilon_{k}^{4}\right] < \infty. Пусть f^N(λ)\widehat{f}_{N}(\lambda ) — периодограмма (49). Найдите lim⁡N→∞cov⁡(f^N(λ),f^N(ν))\lim_{N \rightarrow \infty } \operatorname {cov}\left(\widehat{f}_{N}(\lambda ), \widehat{f}_{N}(\nu )\right) для частот λ,ν∈[−π,π]\lambda , \nu \in [-\pi , \pi ].

?
Задача VII.25

Пусть ZZ — ортогональная случайная мера, отвечающая процессу Орнштейна-Уленбека с параметрами α,β>0\alpha , \beta > 0 (см. упражнение 27 главы III). Докажите, что процесс

W(t)=∫−∞∞eitλ−1iλ⋅iλ+βαZ(dλ),t⩾0 W(t) = \int _{-\infty }^{\infty } \frac{e^{i t \lambda }-1}{i \lambda } \cdot \frac{i \lambda +\beta }{\alpha } Z(d \lambda ), \quad t \geqslant 0

(по непрерывности полагаем (eitλ−1)/iλ=t\left(e^{i t \lambda }-1\right) / i \lambda = t при λ=0\lambda = 0 ) является винеровским процессом. (Существование непрерывой модификации этого процесса легко получается с помощью теоремы Колмогорова (теорема 18 главы II).)

Интересно отметить, что представление (87) дает возможность построения винеровского процесса, отправляясь от процесса Орнштейна-Уленбека, а формула (III.88) показывает, как процесс Орнштейна-Уленбека выражается через винеровский процесс. В главе VIII мы рассмотрим также стохастическое дифференциальное уравнение Ланжевена, связывающее эти два замечательных процесса.

Пользуясь упражнением 4, представление (87) можно переписать в следующем виде (сравните с примером 2 )

W(t)=∫−∞∞eitλ−1iλV(dλ),t⩾0 W(t) = \int _{-\infty }^{\infty } \frac{e^{i t \lambda }-1}{i \lambda } V(d \lambda ), \quad t \geqslant 0

где ортогональная случайная мера VV на полукольце K\mathscr {K}, фигурирующем в упражнении 4 , задается формулой (66) с функцией h(λ)=(iλ+β)/α,λ∈Rh(\lambda ) = (i \lambda +\beta ) / \alpha , \lambda \in \mathbb {R}.

Замечание 7. Согласно теореме 1 главы III траектории винеровского процесса п. н. недифференцируемы ни в одной точке t∈R+t \in \mathbb {R}_{+}. Докажите, что в каждой точке t∈R+t \in \mathbb {R}_{+}не существует производной W(t)W(t) даже в смысле сходимости по вероятности. Тем не менее, если формально продифференцировать по tt равенство (88), получим так называемый "бельй шум"

W˙(t)=∫−∞∞eitλV(dλ) \dot{W}(t) = \int _{-\infty }^{\infty } e^{i t \lambda } V(d \lambda )

Это название связано с тем, что формально получена смесь гармонических колебаний eitλe^{i t \lambda }, каждое из которых входит в (89) с одинаковой интенсивностью (поскольку структурная для VV мера есть мера Лебега). Точный смысл понятию "белый шум" придается в теории обобщенньх случайных процессов (сравните с определением 10); введение в эту область содержится в [77].

Результаты, изложенные в настоящей главе, распространяются также и на процессы с векторными значениями, а также на случайные поля. Остановимся на некоторых из этих результатов.

Матричная функция C(s,t)=(Cjk(s,t))j,k=1m,s,t∈TC(s, t) = \left(C_{j k}(s, t)\right)_{j, k = 1}^{m}, s, t \in T, с комплекснозначными элементами называется неотрицательно (или положительно) определенной, если для всех n⩾1n \geqslant 1, любых t1,…,tn∈Tt_{1}, \ldots , t_{n} \in T и вектор-столбцов z1,…,zn∈Cmz_{1}, \ldots , z_{n} \in \mathbb {C}^{m}

∑j,k=1nzj∗C(tj,tk)zk⩾0 \sum _{j, k = 1}^{n} z_{j}^{*} C\left(t_{j}, t_{k}\right) z_{k} \geqslant 0

где символ * обозначает транспонирование и комплексное сопряжение всех элементов.

Векторный процесс X={X(t),t∈T}X = \left\{ X(t), t \in T\right\} со значениями в Cm\mathbb {C}^{m} называется L2L^{2}-nроцессом, если E[∥X(t)∥2]<∞\mathbb {E}\left[\left\| X(t)\right\|^{2}\right] < \infty для всех t∈Tt \in T. Здесь и далее норма ∥⋅∥\left\| \cdot \right\| в комплексном (в частности, действительном) евклидовом пространстве определяется обычным образом:

∥z∥=(z,z)1/2, где (z,w)=∑k=1mzkwˉk для z,w∈Cm \left\| z\right\| = (z, z)^{1 / 2}, \quad \text{ где } \quad (z, w) = \sum _{k = 1}^{m} z_{k} \bar{w}_{k} \quad \text{ для } \quad z, w \in \mathbb {C}^{m}

Определим векторную функцию среднего и матричную ковариационную (корреляционную) функцию L2L^{2}-процесса XX, положив

a(t)=E[X(t)]∈Cm,r(s,t)=E[(X(s)−a(s))(X(t)−a(t))∗] a(t) = \mathbb {E}\left[X(t)\right] \in \mathbb {C}^{m}, \quad r(s, t) = \mathbb {E}\left[(X(s)-a(s))(X(t)-a(t))^{*}\right]
?
Задача VII.26

Докажите, что матричная функция R(s,t),s,t∈TR(s, t), s, t \in T (с комплексными элементами) является ковариационной функцией некоторого векторного L2L^{2}-процесса X={X(t),t∈T}X = \left\{ X(t), t \in T\right\} тогда и только тогда, когда она неотрицательно определена.

Пусть TT — линейное пространство. L2L^{2}-процесс X={X(t),t∈T}X = \left\{ X(t), t \in T\right\} со значениями в Cm\mathbb {C}^{m} называется стационарным в широком смысле, если

a(t)=a∈Cm,r(s,t)=R(s−t),s,t∈T a(t) = a \in \mathbb {C}^{m}, \quad r(s, t) = R(s-t), \quad s, t \in T

Матричная функция C(t),t∈TC(t), t \in T, называется неотрицательно определенной, если неотрицательно определена функция r(s,t):=C(s−t),s,t,s−t∈Tr(s, t): = C(s-t), s, t, s-t \in T.

Как и прежде, рассматриваемые процессы (и поля) считаем центрированными.

Говорят, что последовательность векторов εn,n∈Z\varepsilon_{n}, n \in \mathbb {Z}, со значениями в Cm\mathbb {C}^{m} является бельм шумом, если в (90)

a=0,R(0)=I,R(n)=0 при n≠0, a = 0, \quad R(0) = I, \quad R(n) = 0 \text{ при } n \neq 0,

здесь II — единичная матрица mm-го порядка. Заметим, что для процессов с непрерывным временем этот термин, как отмечалось выше, имеет более сложный смысл.

Векторный процесс скользящего среднего вводится формулой

X(n)=∑k=−∞∞Akεn−k,n∈Z X(n) = \sum _{k = -\infty }^{\infty } A_{k} \varepsilon _{n-k}, \quad n \in \mathbb {Z}

где ε={εn,n∈Z}\varepsilon = \left\{ \varepsilon_{n}, n \in \mathbb {Z}\right\} — белый шум, Ak,k∈ZA_{k}, k \in \mathbb {Z}, — матрицы (операторы), отображающие Cm\mathbb {C}^{m} в Cm\mathbb {C}^{m}. Ряд в (91) понимается как сходящийся в среднем квадратическом.

?
Задача VII.27

Докажите, что для сходимости в L2(Ω,F,P)L^{2}(\Omega , \mathscr {F}, \mathbb {P}) ряда (91) достаточно, чтобы ∑k=−∞∞∣Ak∣2<∞\sum_{k = -\infty }^{\infty }\left|A_{k}\right|^{2} < \infty, где ∣A∣:=(Tr⁡AA∗)1/2=(∑j,q=1m∣ajq∣2)1/2\left|A\right|: = \left(\operatorname {Tr} A A^{*}\right)^{1 / 2} = \left(\sum_{j, q = 1}^{m}\left|a_{j q}\right|^{2}\right)^{1 / 2} для A=(ajq)j,q=1mA = \left(a_{j q}\right)_{j, q = 1}^{m}. Найдите ковариационную функцию векторного процесса скользящего среднего.

Если TT — метрическое пространство с метрикой ρ\rho, то L2L^{2}-процесс X={X(t),t∈T}X = \left\{ X(t), t \in T\right\} со значениями в Cm\mathbb {C}^{m} называется непрерывным в среднем квадратическом в точке t∈Tt \in T, когда

E[∥X(s)−X(t)∥2]→0 при ρ(s,t)→0. \mathbb {E}\left[\left\| X(s)-X(t)\right\| ^{2}\right] \rightarrow 0 \text{ при } \rho (s, t) \rightarrow 0.

Непрерывность в среднем квадратическом на TT означает выполнение свойства (92) в каждой точке t∈Tt \in T.

Ознакомившись с приложением 4 , докажите, что справедлива следующая

Теорема 21. Функция R=R(t),t∈RdR = R(t), t \in \mathbb {R}^{d}, является ковариационной функцией стационарного в широком смысле непрерывного в среднем квадратическом комплекснозначного случайного поля X={X(t),t∈Rd}X = \left\{ X(t), t \in \mathbb {R}^{d}\right\} тогда и только тогда, когда

R(t)=∫Rdei⟨t,λ⟩G(dλ),t∈Rd R(t) = \int _{\mathbb {R}^{d}} e^{i \left\langle t, \lambda \right\rangle } G(d \lambda ), \quad t \in \mathbb {R}^{d}

где GG — некоторая конечная неотрицательная мера на B(Rd)\mathscr {B}\left(\mathbb {R}^{d}\right). При этом мера GG определяется однозначно.

Аналогичным образом может быть переформулирована и теорема 4.

Комплекснозначное поле X={X(t),t∈Rd}X = \left\{ X(t), t \in \mathbb {R}^{d}\right\} называется изотропным, если его ковариационная функция r(s,t)r(s, t) зависит только от tt и ∥s−t∥\left\| s-t\right\|. Если поле еще и однородно (иначе говоря, является стационарным в широком смысле), то r(s,t)=R(∥s−t∥)r(s, t) = R(\left\| s-t\right\| ) для s,t∈Rds, t \in \mathbb {R}^{d}. В этом случае для R(t)R(t) можно дать иное представление.

Теорема 22 (см. [16; т. 1, с. 262]). Д.ля того чтобы функция R=R(u)R = R(u), где u∈R+u \in \mathbb {R}_{+}, была ковариационной функцией однородного, изотропного и непрерывного в среднем квадратическом комплекснозначного поля X={X(t),t∈Rd}X = \left\{ X(t), t \in \mathbb {R}^{d}\right\}, необходимо и достаточно, чтобы

R(u)=2(m−2)/2Γ(m2)∫0∞I(m−2)/2(λu)(λu)(m−2)/2Q(dλ),u∈R+ R(u) = 2^{(m-2) / 2} \Gamma \left(\frac{m}{2}\right) \int _{0}^{\infty } \frac{I_{(m-2) / 2}(\lambda u)}{(\lambda u)^{(m-2) / 2}} Q(d \lambda ), \quad u \in \mathbb {R}_{+}

здесь Iν(x)I_{\nu }(x) — бесселева функция первого рода, Γ\Gamma — гамма-функция, QQ — неотрицательная конечная мера на B(R+)\mathscr {B}\left(\mathbb {R}_{+}\right), причем Q(R+)=G(Rm)=R(0)Q\left(\mathbb {R}_{+}\right) = G\left(\mathbb {R}^{m}\right) = R(0).

Упомянем также понятие однородного изотропного векторного поля XX, определенного на Rd×Ω\mathbb {R}^{d} \times \Omega и принимающего значения в Rd\mathbb {R}^{d}. Помимо условия (90) здесь требуется инвариантность распределений относительно группы вращений, т.е. требуется, чтобы (S−1X(St1),…,S−1X(Stn))=P(X(t1),…,X(tn))\left(S^{-1} X\left(S t_{1}\right), \ldots , S^{-1} X\left(S t_{n}\right)\right) \stackrel{\mathscr {P}}{ = }\left(X\left(t_{1}\right), \ldots , X\left(t_{n}\right)\right) для любого вращения SS в Rd\mathbb {R}^{d} и любых точек t1,…,tn∈Rd(n∈N)t_{1}, \ldots , t_{n} \in \mathbb {R}^{d}(n \in \mathbb {N}).

Матричную функцию множества G(B)=(Gjk(B))j,k=1mG(B) = \left(G_{j k}(B)\right)_{j, k = 1}^{m}, где B∈KB \in \mathscr {K} и K\mathscr {K} — полукольцо множеств Λ\Lambda, называют неотрицательно определенной, если матрица G(B)G(B) неотрицательно определена при каждом B∈KB \in \mathscr {K}.

Для непрерывного в среднем квадратическом однородного случайного поля X={X(t),t∈Rd}X = \left\{ X(t), t \in \mathbb {R}^{d}\right\} со значениями в Cm\mathbb {C}^{m} и вектора τ∈Cm\tau \in \mathbb {C}^{m} введем скалярное поле Y={Y(t)=⟨X(t),τ⟩,t∈Rd}Y = \left\{ Y(t) = \left\langle X(t), \tau \right\rangle , t \in \mathbb {R}^{d}\right\}. Можно показать, что это поле также является однородным и непрерывным в среднем квадратическом. Пользуясь этим свойством, докажите, что верна

Теорема 23. Для того чтобы функция R=R(t),t∈RdR = R(t), t \in \mathbb {R}^{d}, была матричной ковариационной функцией непрерывного в среднем квадратическом однородного векторного поля X={X(t),t∈Rd}X = \left\{ X(t), t \in \mathbb {R}^{d}\right\} со значениями в Cm\mathbb {C}^{m}, необходимо и достаточно выполнение представления (93), где RR понимается как m×mm \times m матричная функция, а GG — как m×mm \times m матричная счетно-аддитивная функция на B(Rd)\mathscr {B}\left(\mathbb {R}^{d}\right).

Векторной (со значениями в Cm\mathbb {C}^{m} ) ортогональной случайной мерой ZZ на некотором полукольце K\mathscr {K} подмножеств множества Λ\Lambda называется L2L^{2}-процесс {Z(B),B∈K}\left\{ Z(B), B \in \mathscr {K}\right\} такой, что

E[Z(B)Z(C)∗]=G(B∩C) \mathbb {E}\left[Z(B) Z(C)^{*}\right] = G(B \cap C)

где GG — матричнозначная счетно аддитивная мера на K\mathscr {K}, называемая структурной (матричной) мерой. Конструкция интеграла по скалярной ортогональной случайной мере легко обобщается на векторный случай.

Из теорем 3 и 23 выводится (см., например, [16; т. 1, с. 297])

Теорема 24. Поле XX, фигурирующее в теореме 23 , допускает представление вида

X(t)=∫Rdei⟨t,λ⟩Z(dλ),t∈Rd X(t) = \int _{\mathbb {R}^{d}} e^{i \left\langle t, \lambda \right\rangle } Z(d \lambda ), \quad t \in \mathbb {R}^{d}

где ZZ — векторная ортогональная случайная мера на B(Rd)\mathscr {B}\left(\mathbb {R}^{d}\right) со структурной функцией (мерой) GG, отвечающей спектральному представлению ковариационной функции поля. При этом между пространством L2[X]L^{2}[X], т. е. замкнутой в L2(Ω,F,P)L^{2}(\Omega , \mathscr {F}, \mathbb {P}) линейной оболочкой величин X(t),t∈RdX(t), t \in \mathbb {R}^{d}, и пространством L2(ν)=L2(Rd,B(R),ν)L^{2}(\nu ) = L^{2}\left(\mathbb {R}^{d}, \mathscr {B}(\mathbb {R}), \nu \right), где ν=Tr⁡G\nu = \operatorname {Tr} G, устанавливается изометрическое соответствие, при котором

  1. X(t)↔ei⟨t,λ⟩X(t) \leftrightarrow e^{i \left\langle t, \lambda \right\rangle },

  2. если Yk↔gk(λ)Y_{k} \leftrightarrow g_{k}(\lambda ), где Yk∈L2[X],gk∈L2(ν),k=1,2Y_{k} \in L^{2}[X], g_{k} \in L^{2}(\nu ), k = 1,2, то

    Yk=∫Rdgk(λ)Z(dλ),E[Y1Y2∗]=∫Rdg1(λ)g2(λ)∗G(dλ) Y_{k} = \int _{\mathbb {R}^{d}} g_{k}(\lambda ) Z(d \lambda ), \quad \mathbb {E}\left[Y_{1} Y_{2}^{*}\right] = \int _{\mathbb {R}^{d}} g_{1}(\lambda ) g_{2}(\lambda )^{*} G(d \lambda )

Получите отсюда, что справедлива

Теорема 25 (Котельников-Шеннон). Пусть однородное, непрерывное в среднем квадратическом поле X={X(t),t∈Rd}X = \left\{ X(t), t \in \mathbb {R}^{d}\right\} имеет ограниченный спектр, т.е. структурная мера сосредоточена на некотором параллелепипеде (−u1,u1)×⋯×(−ud,ud)\left(-u_{1}, u_{1}\right) \times \cdots \times \left(-u_{d}, u_{d}\right) с ненулевыми длинами ребер. Тогда

X(t)=∑n=(n1,…,nd)∈Zd∏k=1dsin⁡(uktk−πnk)uktk−πnkX(πn1u1,…,πndud) X(t) = \sum _{n = \left(n_{1}, \ldots , n_{d}\right) \in \mathbb {Z}^{d}} \prod _{k = 1}^{d} \frac{\sin \left(u_{k} t_{k}-\pi n_{k}\right)}{u_{k} t_{k}-\pi n_{k}} X\left(\frac{\pi n_{1}}{u_{1}}, \ldots , \frac{\pi n_{d}}{u_{d}}\right)

где при каждом t=(t1,…,td)∈Rdt = \left(t_{1}, \ldots , t_{d}\right) \in \mathbb {R}^{d} ряд сходится в среднем квадратическом.

Смысл представления (94) состоит в том, что оно дает возможность восстанавливать значения поля в произвольной точке t∈Rdt \in \mathbb {R}^{d} по его значениям в узлах решетки вида (πn1/u1,…,πnd/ud)\left(\pi n_{1} / u_{1}, \ldots , \pi n_{d} / u_{d}\right).

В связи с каноническими представлениями случайных процессов упомянем еще одно направление исследований, относящееся к обновляющим процессам. Рассмотрим L2L^{2}-процесс X={X(t),t∈T},T⊂RX = \left\{ X(t), t \in T\right\} , T \subset \mathbb {R}, и подпространства Ht(X),t∈TH_{t}(X), t \in T, состоящие из замыкания в среднем квадратическом линейной оболочки величин X(s)X(s), s⩽t(s,t∈T)s \leqslant t(s, t \in T). Спрашивается, можно ли найти обновляющий процесс Y={Y(t),t∈T}Y = \left\{ Y(t), t \in T\right\}, т. е. процесс с ортогональными приращениями, такой, что

Ht(X)=Ht(Y) при всех t∈T.  H_{t}(X) = H_{t}(Y) \text{ при всех } t \in T \text{. }

Вообще говоря, ответ отрицательный. Однако как показано Г. Крамером, можно видоизменить постановку задачи и получать (при каждом t∈Tt \in T ) представление Ht(X)H_{t}(X) в виде суммы ортогональных друг другу пространств Ht(Yn)H_{t}\left(Y_{n}\right), порожденных должными процессами Yn={Yn(t),t∈T},n=1,…,N(N⩽∞)Y_{n} = \left\{ Y_{n}(t), t \in T\right\} , n = 1, \ldots , N(N \leqslant \infty ). Такого рода представления играют важную роль при исследовании вопросов об абсолютной непрерывности и сингулярности распределений случайных процессов.

Стохастическое спектральное представление (38) можно получить, опираясь на теорию операторов в гильбертовом пространстве HH. Для этого напомним, что изометрический оператор UU (т.е. такой, что (Ux,Uy)=(x,y)(U x, U y) = (x, y) при всех x,y∈Hx, y \in H ), отображающий HH на себя, называется унитарным. Справедлив следующий классический результат.

Теорема 26 (Стоун). Пусть {Ut,t∈R}\left\{ U_{t}, t \in \mathbb {R}\right\} — группа унитарных операторов в гильбертовом пространстве H(UtUs=Us+t,s,t∈R)H\left(U_{t} U_{s} = U_{s+t}, s, t \in \mathbb {R}\right) и функции t↦(ξ,Utη)t \mapsto \left(\xi , U_{t} \eta \right) непрерывны при всех ξ,η∈H\xi , \eta \in H. Тогда для каждого t∈Rt \in \mathbb {R}

Ut=∫−∞∞eitλE(dλ) U_{t} = \int _{-\infty }^{\infty } e^{i t \lambda } E(d \lambda )

где E[⋅]\mathbb {E}\left[\cdot \right] — операторно-значная мера со значениями в ортогональных проекmopax.

О том, как определяется этот интеграл можно прочитать, например, в [64].

Для стационарного в широком смысле процесса X={X(t),t∈R}X = \left\{ X(t), t \in \mathbb {R}\right\}, применив теорему Стоуна к группе унитарных операторов сдвига (St)t∈R\left(S_{t}\right)_{t \in \mathbb {R}} на гильбертовом пространстве H=H(X)H = H(X) (см. доказательство теоремы 12 ), получаем, что

X(t)=StX(0)=(∫ReitλE(dλ))X(0)=∫ReitλZ(dλ) для t∈R, X(t) = S_{t} X(0) = \left(\int _{\mathbb {R}} e^{i t \lambda } E(d \lambda )\right) X(0) = \int _{\mathbb {R}} e^{i t \lambda } Z(d \lambda ) \quad \text{ для } t \in \mathbb {R},

где Z(B)=E[B]X(0)Z(B) = \mathbb {E}\left[B\right] X(0) — ортогональная случайная мера, заданная на борелевских множествах. Тем самым, получено иное (нежели в теореме Крамера (теорема 9 )) доказательство спектрального представления стационарных процессов (сравните с (38)). Приведенный способ получения спектрального представления, опирающийся на теорему Стоуна, был предложен А.Н. Колмогоровым (см. [64]). Аналогичным образом изучаются стационарные процессы с дискретным временем, а также случайные поля на Zd\mathbb {Z}^{d} и Rd\mathbb {R}^{d}.

Отметим, что использованный в главе V IIподход к каноническому представлению случайных функций, основанньй на теореме Карунена, действует для процессов X={X(t),t∈T}X = \left\{ X(t), t \in T\right\}, заданных на произвольном параметрическом множестве TT.

Различным аспектам исследования векторных процессов и полей посвящены монографии [6,44,64,80,87,134,186][6,44,64,80,87,134,186].

Рассмотрим кратко некоторые вопросы, связанные с преобразованиями (главным образом, стационарных) процессов.

Пусть имеется устройство ("система") AA, на вход которой подается процесс X={X(t),t∈T}X = \left\{ X(t), t \in T\right\} и на выходе получается процесс Y={Y(t),t∈S}Y = \left\{ Y(t), t \in S\right\}, который образно можно было бы представить в виде Y=AXY = A X. Если считать, что AA — отображение одного функционального пространства в другое, то надо потребовать, чтобы (почти все) траектории процесса XX лежали бы в области определения DAD_{A} отображения AA и чтобы действие этого отображения на эти траектории приводило также к случайному процессу (в смысле выполнения должных требований измеримости). Будем предполагать, что действие преобразования (отображения) AA на процесс XX приводит в каждый момент t∈St \in S к случайной величине Y(t)Y(t).

Введем операцию дифференцирования в среднем квадратическом. Производную X′(t)X^{\prime }(t) (комплекснозначного) L2L^{2}-процесса XX, заданного в окрестности точки t∈Rt \in \mathbb {R}, определим как производную функции X(t)X(t) со значениями в банаховом (гильбертовом) пространстве L2(Ω)=L2(Ω,F,P)L^{2}(\Omega ) = L^{2}(\Omega , \mathscr {F}, \mathbb {P}), т.е. как такой случайный элемент (обозначаемый X′(t)X^{\prime }(t) ), для которого

E[∣X′(t)−(X(t+h)−X(t))/h∣2]→0 при h→0 \mathbb {E}\left[\left|X^{\prime }(t)-(X(t+h)-X(t)) / h\right|^{2}\right] \rightarrow 0 \text{ при } h \rightarrow 0

Данный подход позволяет создать содержательное исчисление (с аналогом классической формулы Ньютона-Лейбница). Отметим, что если вместо сходимости в среднем квадратическом (95) требовать лишь сходимость по вероятности, то соответствующее исчисление неприводит бы к плодотворной теории как показывает утверждение следующего упражнения.

?
Задача VII.28

Докажите, что для пуассоновского процесса N={Nt,t⩾0}N = \left\{ N_{t}, t \geqslant 0\right\} постоянной интенсивности λ>0\lambda > 0 производная по вероятности в каждой точке t∈R+t \in \mathbb {R}_{+} равна нулю п. н.

В дальнейшем мы рассматриваем лишь производные в среднем квадратическом, не оговаривая этого специально. При этом, как обычно, производная kk-го порядка в точке tt определяется (если она существует) по значениям процесса X(k−1)(t)X^{(k-1)}(t) в окрестности этой точки.

?
Задача VII.29

Пусть X={X(t),t∈[a,b]}X = \left\{ X(t), t \in [a, b]\right\} — это L2L^{2}-процесс с ковариационной функцией r(s,t)r(s, t). Докажите, что X′(t)X^{\prime }(t) существует в точке t∈(a,b)t \in (a, b) тогда и только тогда, когда существует обобщенная вторая производная ∂~2r(s,t)∂s∂t\frac{\widetilde{\partial }^{2} r(s, t)}{\partial s \partial t} в точке (t,t)(t, t) (с теорией обобщенных функций это никак не связано), где по определению

∂~2r(s,t)∂s∂t:=lim⁡h,u→0(r(s+h,t+u)−r(s+h,t)−r(s,t+u)+r(s,t))/hu. \frac{\widetilde{\partial }^{2} r(s, t)}{\partial s \partial t}: = \lim _{h, u \rightarrow 0}(r(s+h, t+u)-r(s+h, t)-r(s, t+u)+r(s, t)) / h u .

Аналогичное утверждение справедливо и для производной в концевых точках, если иметь в виду односторонние производные.

Для доказательства этого и подобных утверждений ключевую роль играет следующая элементарная

Лемма 5. Функция f(t),t∈[a,b]f(t), t \in [a, b], принимающая значения в гильбертовом пространстве НН со скалярным произведением ⟨⋅,⋅⟩\langle \cdot , \cdot \rangle, имеет предел в точке t0∈[a,b]t_{0} \in [a, b] тогда и только тогда, когда существует предел lim⁡s,t→t0⟨f(s),f(t)⟩\lim_{s, t \rightarrow t_{0}}\langle f(s), f(t)\rangle.

?
Задача VII.30

Докажите, что у стационарного процесса X={X(t),t∈R}X = \left\{ X(t), t \in \mathbb {R}\right\} производная X(k)(t)X^{(k)}(t) существует в каждой точке tt тогда и только тогда, когда

∫−∞∞λ2kG(dλ)<∞ \int _{-\infty }^{\infty } \lambda ^{2 k} G(d \lambda ) < \infty

где GG — спектральная мера XX. Будет ли процесс X(k)X^{(k)} стационарен? Если будет, то каково его спектральное представление?

Наряду с дифференцированием важный пример линейной операции AA, т.е. такой, что

A(αX+βV)=αAX+βAV A(\alpha X+\beta V) = \alpha A X+\beta A V

для α,β∈C\alpha , \beta \in \mathbb {C} и X,V∈DAX, V \in \mathrm{D}_{A}, дает интегрирование в среднем квадратическом случайных процессов.

Будем для L2L^{2}-процесса X={X(t),t∈[a,b]}X = \left\{ X(t), t \in [a, b]\right\} под интегралом ∫[a,b]X(t)dt\int_{[a, b]} X(t) d t (или ∫abX(t)dt\int_{a}^{b} X(t) d t ) понимать интеграл Римана (если он определен) по отрезку [a,b][a, b] от случайной функции XtX_{t}, т.е. берется предел в пространстве L2(Ω)L^{2}(\Omega ) интегральных сумм

Римана. Подробнее об интегрировании и дифференцировании в среднем квадратическом см., например, [12; гл. 2].

Несобственные интегралы по всей прямой R\mathbb {R} вводятся стандартным образом как пределы (если таковые существуют) в соответствующем смысле интегралов по отрезкам [a,b][a, b], когда a→−∞a \rightarrow -\infty и b→∞b \rightarrow \infty.

Для L2L^{2}-процесса X={X(t),t∈R}X = \left\{ X(t), t \in \mathbb {R}\right\} с ковариационной функцией r=r(s,t)r = r(s, t) и комплекснозначной функции hh, заданной на R×R\mathbb {R} \times \mathbb {R}, положим

Y(t)=∫−∞∞h(t,s)X(s)ds,t∈R Y(t) = \int _{-\infty }^{\infty } h(t, s) X(s) d s, \quad t \in \mathbb {R}

Нетрудно убедиться, что процесс Y(t)Y(t) в (96) определен и конечен (п. н.), если конечен несобственный интеграл Римана

∫−∞∞∫−∞∞h(t,s)r(s,u)h(t,u)‾dsdu \int _{-\infty }^{\infty } \int _{-\infty }^{\infty } h(t, s) r(s, u) \overline{h(t, u)} d s d u

Функция h=h(s,t)h = h(s, t), фигурирующая в (96), называется в приложениях импульсной переходной функцией системы. Если формально подставить в интеграл (96) вместо X(s)X(s) функцию Дирака δs\delta_{s}, то получим Y(t)=h(t,s)Y(t) = h(t, s). Другими словами, h(t,s)h(t, s) есть отклик системы в момент tt на δ\delta-импульс, подаваемый на ее вход в момент ss.

Частным случаем процессов вида (96) являются однородные по времени процессы вида

Y(t)=∫−∞∞h(t−s)X(s)ds,t∈R Y(t) = \int _{-\infty }^{\infty } h(t-s) X(s) d s, \quad t \in \mathbb {R}

Соответствующие системы ("фильтры"), на выходе которых возникает процесс Y={Y(t),t∈R}Y = \left\{ Y(t), t \in \mathbb {R}\right\}, определяются импульсной переходной функцией h=h(s)h = h(s), s∈Rs \in \mathbb {R}, или функцией

H(iλ)=∫−∞∞h(s)e−iλsds,λ∈R H(i \lambda ) = \int _{-\infty }^{\infty } h(s) e^{-i \lambda s} d s, \quad \lambda \in \mathbb {R}

называемой частотной характеристикой или коэффициентом передачи. Если функция h=h(s)h = h(s) интегрируема по Лебегу, то функция H=H(iλ)H = H(i \lambda ), очевидно, определена. Частотная характеристика имеет следующий наглядный смысл. Если h∈L1(R)h \in L^{1}(\mathbb {R}), то для каждого λ∈R\lambda \in \mathbb {R} неслучайная функция X(s)=eiλs,s∈RX(s) = e^{i \lambda s}, s \in \mathbb {R}, будет собственной функцией преобразования (97), отвечающей собственному значению H(iλ)H(i \lambda ). Это делает понятным смысл следующего упражнения.

?
Задача VII.31

Пусть на вход фильтра (97) подается стационарный в широком смысле процесс X={X(t),t∈R}X = \left\{ X(t), t \in \mathbb {R}\right\}, имеющий спектральное представление (38) с ортогональной случайной мерой ZZ, структурная мера которой есть GG. Докажите, что если частотная характеристика фильтра H(iλ)∈L2(R,B(R),G)H(i \lambda ) \in L^{2}(\mathbb {R}, \mathscr {B}(\mathbb {R}), G), то ковариационная функция процесса Y={Y(t),t∈R}Y = \left\{ Y(t), t \in \mathbb {R}\right\} имеет вид

R(t)=∫−∞∞eiλt∣H(iλ)∣2G(dλ),t∈R R(t) = \int _{-\infty }^{\infty } e^{i \lambda t}\left|H(i \lambda )\right|^{2} G(d \lambda ), \quad t \in \mathbb {R}

а сам процесс YY допускает следующее спектральное представление:

Y(t)=∫−∞∞eiλtH(iλ)Z(dλ),t∈R Y(t) = \int _{-\infty }^{\infty } e^{i \lambda t} H(i \lambda ) Z(d \lambda ), \quad t \in \mathbb {R}

Вспоминая энергетическую интерпретацию спектральной меры, видим, что функция ∣H(iλ)∣2\left|H(i \lambda )\right|^{2} показывает, во сколько раз изменяется энергия простых гармонических составляющих процесса XX при прохождении через фильтр с частотной характеристикой H=H(iλ)H = H(i \lambda ).

?
Задача VII.32

Какая импульсная переходная функция отвечает полосово́му фильтру, т.е. фильтру, пропускающему без изменения гармонические составляющие процессов с частотами в полосе [a,b][a, b] (это значит, что H(iλ)=1[a,b](λ),λ∈RH(i \lambda ) = \mathbf{1}_{[a, b]}(\lambda ), \lambda \in \mathbb {R} )?

Интересно применение спектральных представлений процессов к решению дифференциальных уравнений вида

Pn(ddt)Y=Qm(ddt)X P_{n}\left(\frac{d}{d t}\right) Y = Q_{m}\left(\frac{d}{d t}\right) X

где Pn(z)=a0zn+a1zn−1+⋯+an,Qm(z)=b0zm+b1zm−1+⋯+bmP_{n}(z) = a_{0} z^{n}+a_{1} z^{n-1}+\cdots +a_{n}, Q_{m}(z) = b_{0} z^{m}+b_{1} z^{m-1}+\cdots +b_{m} — многочлены с постоянными коэффициентами, XX — заданный стационарный процесс, а неизвестный процесс Y={Y(t),t∈R}Y = \left\{ Y(t), t \in \mathbb {R}\right\} предполагается дифференцируемым в среднем квадратическом nn раз в каждой точке t∈Rt \in \mathbb {R}. Равенство (98) считается выполненным п. н. при каждом t∈Rt \in \mathbb {R}. О том, как решение YY ищется с помощью некоторого фильтра см., например, [17; с. 282-284].

Аналог уравнения (98), рассмотренного вьше для процессов с непрерывным временем, приводит в случае процессов с дискретным временем (t∈T⊂Z)(t \in T \subset \mathbb {Z}) к изучению так называемых процессов "авторегрессии — скользящего среднего", см. подробнее [85[85; т. 2 , с. 586]].

В заключение приведем еще один классический результат, относящийся к теории процессов, стационарных в узком смысле.

Введем на пространстве R∞\mathbb {R}^{\infty } оператор сдвига TT, положив Tx=yT x = y, где

x=(…,x−1,x0,x1,…),y=(…,y−1,y0,y1,…) x = \left(\ldots , x_{-1}, x_{0}, x_{1}, \ldots \right), \quad y = \left(\ldots , y_{-1}, y_{0}, y_{1}, \ldots \right)

причем yn=xn+1,n∈Zy_{n} = x_{n+1}, n \in \mathbb {Z}. Назовем множество B∈B(R∞)B \in \mathscr {B}\left(\mathbb {R}^{\infty }\right) инвариантным, если T−1(B)=BT^{-1}(B) = B. Легко видеть, что совокупность инвариантных множеств образует в B(R∞)σ\mathscr {B}\left(\mathbb {R}^{\infty }\right) \sigma-алгебру, которую обозначим I\mathscr {I}.

Пусть X={Xt,t∈Z}X = \left\{ X_{t}, t \in \mathbb {Z}\right\} — стационарный в узком смысле процесс, заданный на вероятностном пространстве (Ω,F,P)(\Omega , \mathscr {F}, \mathbb {P}). Положим IX=X−1(I)\mathscr {I}_{X} = X^{-1}(\mathscr {I}), где XX рассматривается как измеримое отображение (Ω,F)(\Omega , \mathscr {F}) в (R∞,B(R∞))\left(\mathbb {R}^{\infty }, \mathscr {B}\left(\mathbb {R}^{\infty }\right)\right). Тогда IX\mathscr {I}_{X} будет σ\sigma-алгеброй в F\mathscr {F}, называемой σ\sigma-алгеброй инвариантных событий процесса XX.

Теорема 27 (Биркгоф-Хинчин). Пусть X={X(t),t∈Z}X = \left\{ X(t), t \in \mathbb {Z}\right\} — стачионарный в узком смысле действительный процесс и E[∣X0∣]<∞\mathbb {E}\left[\left|X_{0}\right|\right] < \infty. Тогда

lim⁡N−M→∞1N−M∑t=M+1NX(t)=E[X0∣IX] \lim _{N-M \rightarrow \infty } \frac{1}{N-M} \sum _{t = M+1}^{N} X(t) = \mathbb {E}\left[X_{0} \mid \mathscr {I}_{X}\right]

где предел существует как почти наверное, так и в L1(Ω)L^{1}(\Omega ).

Говорят, что процесс XX является әргодическим, если P(A)\mathbb {P}\left(A\right) равно 0 или 1 для любого A∈IXA \in \mathscr {I}_{X}. В этом случае (при выполнении условий теоремы 27 ) в правой части (99) условное математическое ожидание E[X0∣IX]=E[X0]\mathbb {E}\left[X_{0} \mid \mathscr {I}_{X}\right] = \mathbb {E}\left[X_{0}\right] п.н.

Теория стационарных в узком смысле процессов теснейшим образом связана с теорией динамических систем (с изучением преобразований, сохраняющих меру). В качестве введения в эту область исследований можно прочитать гл. V в [85], где, в частности, дано простое доказательство соотношения (99) (при M=0M = 0 ), предложенное А. Гарсиа. Обобщение теоремы 27 на случайные поля содержится в работе [15]. Эргодической теории и динамическим системам посвящены, например, книги [36,70][36,70].

?