Метод оценки вероятностей больших выбросов случайных процессов Текст научной статьи по специальности «Математика»

Том XXXIV

УЧЕНЫЕ ЗАПИСКИ ЦАГИ 2 00 3

№ 3—4

УДК 519.2

МЕТОД ОЦЕНКИ ВЕРОЯТНОСТЕЙ БОЛЬШИХ ВЫБРОСОВ СЛУЧАЙНЫХ ПРОЦЕССОВ

В. П. КУЗЬМИН

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

пересечения

уровня

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

1. Постановка задачи и методика оценки вероятности. Рассматривается задача об оценке вероятности первого достижения случайным процессом х (ґ) заданного уровня х0 за время ґ. Будем считать, что уровень Хо достаточно высок, так что среднее время достижения такого уровня значительно превосходит время корреляции для данного процесса, которое может, например, определяться формулой [1]:

где Кх (т) — корреляционная функция случайного процесса.

Считается, что начальные условия находятся достаточно далеко от заданного уровня. В этом случае вероятность достижения заданного уровня за время ^ представляется в виде [1], [2]

и задача сводится к определению параметра д, который связан со средним временем достижения

В ряде случаев не трудно определить среднее число всех пересечений в единицу времени по формуле [1], [2]:

Р(ґ) * 1 - в~^

уровня М ] = — и равен среднему числу первых пересечений в единицу времени.

то

(1)

о

где у = ~^~ — производная случайной величины х ^), /(х, у) — совместная плотность

распределения величин х и у. Распределение производной в момент пересечения определяется соотношением:

А ( у ) = -

о

В случае стационарного гауссовского процесса из формулы (1) получаем [1], [2]:

\2 '

N =

С

у

2пс

ехр

( х0 - т х )

2с2

где тх = М [ х] — среднее значение переменной х (t), сх и су — средние квадратические отклонения переменной х ^) и ее производной.

Производная в момент пересечения будет иметь распределение Релея

Л (у ) = ^гехР

V 2су у

(2)

Пересечения, соответствующие одной реализации, будем называть серией [1], тогда среднее число всех пересечений N0 и среднее число первых пересечений (реализаций) д связаны соотношением

N

д = -

п*

(3)

где п* — среднее число пересечений в одной серии.

В работе [3] предложен способ определения среднего числа пересечений в серии, основанный на подсчете числа пересечений для нестационарного процесса, начинающегося от момента первого пересечения, и получено соотношение для определения величины п*

п* = 1 + п = 1 + |( N ^)- Ы0 ) Л,

(4)

где N ^) — среднее число пересечений в единицу времени для нестационарного случайного

процесса, начинающегося от момента первого пересечения; п — среднее число пересечений за первым пересечением.

Если нестационарный процесс является гауссовским, то среднее число пересечений с положительной производной определяется соотношением [2], [3]:

хо - тх ^)

С х(t)

А (М+),

(5)

где

f (M+) = exp

-1 m2

-2п М+Ф[М+ (t)], M+ =

xy

my(t) + x0 - mx(t).

C У (t) C У (t)

1 x

Ф(x) = ^= [ exp(-u2/2)du — функция Лапласа. J2n J v '

Для недифференцируемого процесса бесконечно большими оказываются как среднее число всех пересечений в единицу времени, так и среднее число пересечений в одной серии. В этом случае значение параметра д находится как предел при стремлении к нулю интервала корреляции и для гауссовского процесса определяется формулой [3]:

exp

( ( \2 >

(x0 - mx )

2с2

ц = -

где параметр щ определяется интегралом

Cx (t)

exp

(x0 - mx (t))

' 2с2 (t)

exp

(x0 - mx )2

2с2

(6)

Таким образом, решение задачи определения параметра д сводится к интегрированию, если известно изменение по времени условных моментов распределения для нестационарного процесса, следующего за первым пересечением.

Основная сложность в использовании предложенного способа заключается в необходимости выделения из всех пересечений первых в данной серии [3]. В данной работе предложен простой способ моделирования всех пересечений (случайных реализаций, пересекающих заданный уровень) и выделения среди них первых пересечений.

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

(л □ ткор ), связанного с начальными условиями, являются стационарными и гауссовскими.

Рассмотрим все реализации, пересекающие заданный уровень, от момента пересечения как функции интервала времени А = t - Л0), где Л0) — момент пересечения. Совокупность таких реализаций образует нестационарный случайный процесс, поскольку для него х(АЛ = 0) = х0 и, следовательно, В[х(0)] = 0, а при больших значениях | А | этот случайный процесс практически

не отличается от исходного. Определим вначале в общем случае соотношения, которым удовлетворяют моменты нестационарного процесса с заданными условиями пересечения уровня.

Пусть исходный стационарный случайный процесс со спектральной плотностью ^ (ю) задан рядом Фурье

x(t) = 14 (^г- cos Ю/ + Q sin ю/),

(7)

где А2 = ^ (ю г-)dю, £, г- и дг- — независимые стандартные (с нулевым средним и единичной

дисперсией) гауссовские случайные величины. Число членов ряда может быть любым, в том числе и бесконечным.

Отметим, что в общем случае стационарный случайный процесс задается стохастическим интегралом [6]. Запись в виде ряда (7) связана с тем, что полученные далее формулы используются для численных расчетов, где число членов ряда всегда конечно, а вывод формул

n

t

00

nt = И

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

Объединим все случайные величины в один вектор и запишем разложение (7) в виде скалярного произведения

x(t) = ( at 5),

(8)

где вектор \ состоит из стандартных гауссовских величин, а компоненты вектора at определяются формулами:

ati = Ai cos (ю it) для i=1, 2,..., n, ati = Ai sin (ю/) для i=n + 1, n + 2, ...2n.

Производная случайного процесса определяется соотношением

dx(t)

У(t) = -

dt

= ( b 5).

(9)

где bti =-Airni sin (rnit) для i = 1, 2,., n, bti = Airai cos (ю/) для i =n+1, n + 2, ...2n.

Пусть в некоторый момент времени to случайная реализация пересекает уровень xo с производной Уо . Эти условия запишутся в виде

(ao5) = ^ (bo5) = Уo,

(і0)

где а0 = ал (Л0), ь0 = Ъл (Л0)•

Векторы а0), Ъ), аЛ, Ъ( удовлетворяют очевидным соотношениям:

(а0Ъ0) = (аЛЪЛ) = 0 — условие некоррелированности переменной и ее производной для стационарного процесса;

(а0аЛ) = Кх (Л - ^ ) — корреляционная функция переменной х(Л);

(Ъ0ЪЛ) = Ку (л - Л0 ) — корреляционная функция производной у(Л);

/ \ \ dKx (л - Л0)

(а0ЪЛ) = - (Ъ0аЛ) =--—-----= Кх (Л - Л0) — производная от корреляционной функции.

Случайные реализации, удовлетворяющие условиям (10), определяются случайным вектором ^, в котором зафиксированы два ортогональных компонента, соответствующие

переменной х (Л) и ее производной в момент времени ^

~ ^0 uo

5 = 9 xo +-уУ0 +n-

ao (nao) bo (nbo)

(іі)

В соотношении (11) вектор п состоит из независимых стандартных гауссовских случайных величин.

Подставим случайный вектор ^ в соотношения (8), (9) и получим выражения для изменения переменных х (л) и у (л) на случайных реализациях, пересекающих заданный уровень с заданной производной:

:(t - to ) =

(aoat) , (boat)

(aobt) , (bobt)

У (t - t0 )^^“2 x0 +

У0

У0

f f n

v v f f n

v v

(aoat) (boat)

.2 bo U2

a

o

bt - ao

(aobt) , (bobt)

X- - bo

(і2)

(іЗ)

//

Используя соотношения (12) и (13), можно получить условные моменты случайной величины х (Л) и ее производной, как функции от интервала времени А = Л - ^:

»ГГ ,А м -х(Лґ) К'х(Лґ)

т(Лґ) = и[х(Лґ)]=—х0--К^у°,

/л ч ,гг /л м -У (Лґ) -х (Лґ)

тУ(Лґ)=и [ У (Лґ )]=—у0+—х„,

(14)

"2(Лґ) = В [ х(Лґ)] = Кх (0 -)-х(Лґ) - -х (Лґ), л 1 -1 х Кх (0) -у (0)

(15)

'(Лґ) = В[у(Лґ)] = Ку(0 -)

гху (Лґ )а х (Лґ )с у (Лґ) = -ху (Лґ) = К'х

КУ(Лґ) - -х 2(Лґ)

Ку (0) -х (0) ,

КУ (Лґ) -х (Лґ)

ху'

-у (0) -х(0)

где гху — коэффициент корреляции между переменными х и у,

-Л х ( ґ1, ґ2 ) = -х ( ґ2 ґ1)

-х ( ґ1 )-х (ґ2 ) -х(ґ1 )-х(ґ2 )

-х (0)

-У (0)

(16)

Если переменная х (ґ) недифференцируемая, то у нее нет производной, и выражения для

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

(а0^) = х0-

Случайный вектор, определяющий все пересечения, задается в виде

*0 (п«0 )

(17)

Случайный процесс определяется соотношением

(. , \ (а0аґ)

х (ґ - ґ0 ) = ^Т~ х0

( ґ

аґ а0 2 V V а02 )

(а0аґ)

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

и [ х(Лґ )] = ~~~~~ х0,

-х (0)

В [ х(Лґ )] = -х (0) -

-х2(Лґ) -х (0) ,

(18)

(19)

Для стационарного процесса значение момента времени Л0 несущественно, поэтому в дальнейшем будем считать Л0 = 0, т. е. будем отсчитывать время от момента пересечения.

В этом случае коэффициенты для косинусов и синусов в разложении случайного процесса (7) вычисляются независимо:

где П1 и П2 — векторы, состоящие из независимых стандартных гауссовских случайных величин.

3. Методика выделения первых пересечений. Приведенные выше соотношения позволяют моделировать все случайные пересечения, а выделение среди них первых сводится к проверке очевидного условия

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

Сравнивая соотношение (24) с формулой (3), легко установить, что /1 = 1/п*.

Для реализаций, соответствующих первым пересечениям, можно также определить оценки моментов случайного процесса при Л > 0 и с помощью формул (5) и (4) найти среднее число пересечений после первого п и среднее в серии п* = п + 1.

Проведение расчетов по определению среднего числа пересечений в серии может показаться лишним, поскольку известное значение /1 решает поставленную задачу. Однако при проведении численных расчетов часто оказывается, что величина п более устойчива к ошибкам вычислений, чем величина вероятности Р1. Кроме того, условие (23) проверяется на интервале времени, который не определен строго и практически определяется как интервал сходимости интегралов (4) или (6).

В случае недифференцируемого случайного процесса вероятность Р1 = 0, поскольку среднее число всех пересечений бесконечно велико. В этом случае предлагается проводить параметрические расчеты для нескольких, уменьшающихся по величине значений шага АЛ. Количество реализаций, для которых выполняется условие (23) и, следовательно, величина щ, определенная для этих реализаций, будет зависеть от величины шага и искомое значение является пределом

Очевидно также, что Нш Р^АЛ) = 0, поэтому успех данной процедуры определяется

АЛ ^0

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

с малым шагом по времени.

(21)

(22)

х(ґ) < х0 для -ткор < ґ < 0.

(23)

(24)

пґ = ііш пґ (Лґ).

Лґ——0

4. Методы формирования случайных реализаций. Приведенные выше формулы для случайных векторов (11), (17), определяющих все реализации, или формулы для корреляционных функций (16), (20) позволяют формировать случайные реализации различными способами. Рассмотрим два возможных способа.

4.1. Формирование случайной реализации в дискретные моменты времени. Задается дискретный набор моментов времени Л, для / = 1, 2, ..., т из заданного диапазона, например,

-ткор < Л <ткор. Случайная реализация задается вектором с компонентами X = Ах (). Для

данного набора моментов времени определяется корреляционная матрица К размера т х т по соотношению (16) для дифференцируемого процесса, либо по соотношению (20), если процесс недифференцируемый. Элементами этой матрицы являются значения корреляционной функции

в соответствующие моменты времени Ку = Ках (лг-, Л у ).

Матрица К является симметричной, и поэтому она может быть представлена в виде [5]:

Тогда случайный вектор Лх определяется через случайный вектор состоящий из стандартных гауссовских величин, по формуле

К полученным значениям необходимо прибавить математическое ожидание, определяемое соответственно по формулам (14) либо (18)

4.2. Формирование случайных реализаций с помощью ряда Фурье. Способ формирования случайных реализаций с помощью ряда Фурье представляется технически особенно простым.

Для заданного спектрального разложения (7) определяются два вектора а0 и Ь0 с компонентами ат = Аг-, Ът = Дюг-.

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

Для обеспечения условий в момент пересечения случайные векторы П1 и П2 преобразуются по формулам (21), (22). Если случайный процесс недифференцируемый, преобразуется только вектор гц. Случайная реализация в произвольный момент времени определяется с помощью ряда (7).

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

где с и С2 — стандартные гауссовские величины.

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

К = ЄЄ7.

Лх =

х (ґг ) = Шх (і, ) + Лхг-.

где Г (1/2 )Г(р -1/2)

задается таким образом, что | S (ю) й ю = 1

0

Случайный процесс с такой спектральной плотностью имеет конечную дисперсию при

1 3 3

р> 2, пр„че„ при р< 2 о„ -вляется н£диффбр6нциру6мым: а при р> 2 — дпффере„ц„руемым

1

с дисперсией производной, равной с у = -

2 (р-3/2)

Корреляционная функция данного процесса определяется соотношением

Кх (Л) =

где Ку (Л) — модифицированная функция Бесселя, Г(у) — гамма-функция, у = р-1/2. Интервал корреляции имеет значения

Ткор = ] Кх (Л)йЛ =

Г (1/2 )Г(у + У2 ) Г(у) .

Рассматривается задача определения вероятности первого пересечения уровня Х0 = Ясх. Из анализа корреляционных функций для нестационарного процесса следует, что лишь при целых значениях р = 1 и р = 2 нестационарный процесс является марковским. В этих случаях

Ках (^, Л2 ) = 0, если переменные ^ и Л2 имеют разные знаки. Этот факт можно легко проверить,

если учесть, что Кх = е 111 для р = 1 и Кх = е I(1 + | Л |) для р = 2.

Равенство нулю корреляционной функции означает, что задание одного или двух условий в момент времени Л = 0 делает процессы до и после пересечения некоррелированными. В частности, для р = 1 это означает, что выделение первых пересечений среди всех не изменит моментов распределения при Л > 0.

Для р = 2 некоррелированность процесса обеспечивается двумя условиями (10), и процессы при Л < 0 и Л > 0 связаны между собой одним параметром у — производной в момент пересечения. Следовательно, выделение первых пересечений среди всех приведет лишь к изменению распределения производной при первом пересечении, но не изменит зависимости моментов случайного процесса от времени для Л > 0 по сравнению со всеми пересечениями при заданном значении производной.

При дробных степенях р процесс не является марковским, и это означает, что значения

х (л) при Л > 0 зависят от всех значений случайной функции при Л < 0, и эта связь не может быть

сведена к конечному числу условий и, следовательно, конечному числу переменных.

В дальнейшем величины, определенные для всех пересечений, будем отмечать индексом «0», а вычисленные только для первых пересечений — индексом «1». Кроме того, для дифференцируемого процесса также будем использовать параметр п{, определяемый соотношением п = 2 пп*сх/су.

5.1. Случай недифференцируемого процесса. Рассмотрим вначале случай р<3/2, когда случайный процесс является недифференцируемым, и вычислим величину пЛ0 для всех пересечений по формуле (6), которая с учетом формул (18), (19) и условия Кх (0) = 1 имеет вид:

1

ехр

Я2 (1 - Кх (Л))2

2 1 - Кх2 (Л)

Я 2 ]

- ехр

2

(25)

Результаты численных расчетов по формуле (25) приведены на рис. 1.

Значение пЛ0 может быть вычислено по приближенной формуле, которая для высоких уровней получается из предположения, что основной вклад в интеграл (25) дает область малых значений Л.

При малых значениях Л из разложения функции Бесселя получим [4]

где

Кх (ґ) * 1 - Ьґ Г(1 -V)

2у

Ь =

2v '

Тогда

Г(1 + v)2

Я — тх (ґ) * ЯЬґ2v, с2 * 2Ьґ2v,

1 1

п

і

у * [-----й2— ехР (-Я^Ьґ2v/4)

ґ 0(2Ь)12 ґv V

1 1

2v^2A

' 4 " 2v 2 " 1 1'

_ Я 2Ь _ _ 2v 2 _

(26)

где Г( х) — гамма-функция.

В общем случае, если спектральная плотность нормирована (с2 = 1) и при больших

частотах ведет себя как ^(со) -

ю

2р

, для больших значений Я справедлива формула (26) с

параметром Ь, определяемым по формуле:

Ь = А

^Лг( 3/2 -р)

22р(р-1/2 )Г(р)'

Анализ формулы (26) показывает, что показатель степени при величине Я меняется от нуля

при р = 3/2 до бесконечности при р ^ 1/2 . При

( л Л

величина

р —— 1/2 более точно при —1— □ Я2

2р — 1

определяется формулой:

пу *(2р — 1)(Я2 +1)е_

Тогда среднее время достижения уровня Я составит

М [ґ1 ]*( 2р — 1)( Я 2 +1)

и будет стремиться к нулю при р ^ У2 независимо от величины уровня.

2

Для выделения первых пересечений используется метод формирования случайных значений х (7) в дискретных точках. Для этого задаются дискретные моменты времени, симметричные относительно значения ^ = 0. Заданные узлы сгущаются при приближении к нулю по геометрической прогрессии. Минимальное значение шага А^ш;п задается, а знаменатель прогрессии подбирается так, чтобы заданное число узлов примерно покрывало интервал корреляции. В расчетах было принято общее число узлов т, равное 160. Таким образом, условие (23) отсутствия пересечений при t < 0 проверяется в 80 дискретных точках.

\ V > 4х ** ■ — . V Ч. ’''■■‘V.— " о» _

ч ч Ч -V ч ч ч ч ч X N ч . > '"'1,48 'м4

ч ч ч ч ч ^ ■ '1,25 'х Р = 1 Ч

2 ' 4 6 8 ' " '11

•в"/!

-1нД/т

1,48

1,25

и,7

2 ' V '6 а 10

Рис. 3

Величина пл определяется численно по формуле (6), где в качестве моментов распределения используются оценки, полученные для первых пересечений, т. е. для случайных реализаций, в которых условие (23) выполняется во всех точках при ^ < 0.

При проведении статистических испытаний задается очень большое число реализаций,

порядка 106, но расчеты проводятся до появления 1000 первых пересечений.

Расчеты показали, что основное влияние на величину п^ оказывает смещение математического ожидания у первых пересечений по сравнению со всеми пересечениями. Примеры сравнений математических ожиданий для первых и для всех пересечений приведены на рис. 2. На этом рисунке приведено изменение по времени безразмерных математических ожиданий, отсчитываемых от заданного уровня: Атх ^) = Я -(х0 - тх ) Vе х (*). На рис. 3 приведены характерные процессы уменьшения вероятности /1 и процесс сходимости параметра п(1 при уменьшении величины минимального шага А^ш;п.

5.2. Случай дифференцируемого случайного процесса. В случае дифференцируемого

Рис. 2

случайного процесса ((р>3/2) для определения величины п* проводятся параметрические расчеты при различных значениях производной в момент пересечения, в результате чего определяются зависимости п (у0 ) и р (у0 ). Окончательные значения п и Р1 получаются путем осреднения по случайным значениям производной.

Распределение производной для первых пересечений будет иметь вид:

Р (Уо)

^ехр

/1(Уо ) = ;

Уо 2с 2

] р1 (У)

У

ехр

У

2с 2

(27)

а величины р и п определятся интегралами

Р1 = | Р1 (У ) /1 ( У) ^У’ п = |п (У ) /1 (У) ^у.

о о

При проведении численных расчетов задавалось 20 значений производной из диапазона (о - 4с У ).

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

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

В данном примере выделение первых пересечений приводит как к изменению распределения производной (рис. 4, уо = уо/су), так и изменяет моменты распределения при

^ > о, что приводит к существенному изменению среднего числа пересечений в серии, в первую очередь — при значениях р, близких

к величине 3/2.

Суммарные результаты расчетов для различных значений степени р приведены на рис. 5 для Я = 5

и Я = 1о. Здесь приведены значения параметра пп, определенные по описанной методике (символы *), результаты расчетов параметра п{о по всем пересечениям (сплошные линии) и приближенные оценки, которые для недифференцируемого процесса (р < 3/2) соответствуют формуле (26), а для дифференцируемого — простому условию Ц = ^о (штриховые линии).

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

Рис. 5

5.3. Случай дифференцируемого колебательного процесса. В качестве другого примера рассмотрим случайный процесс со спектральной плотностью вида

С

5 (ю) = -

(ю2 -1)2 + 4к2ю2 для к □ 1.

Значение константы С несущественно, поскольку рассматриваются безразмерные значения заданных уровней Я = ± —°. В данном примере контролируется выход за два уровня. При р = 1

с х

приведенная спектральная плотность соответствует дифференциальному уравнению второго порядка, рассмотренному в работах [1], [3].

Для формирования случайных реализаций задавалось 60 членов ряда Фурье. Интервалы по

частоте сгущаются к значению Юо = ^ 1 - 2к2 по геометрической прогрессии, максимальное

значение частоты задавалось равным 3. Минимальное значение интервала по частоте и интервал времени, в течение которого подсчитывались пересечения, определялись по формулам

Аюш;п = о,2 • к/ р, М* =р/к. Отметим, что в данном случае интервал времени для подсчета числа пересечений лишь примерно соответствует интервалу корреляции.

Расчеты проводились одной серией со случайным распределением производной в момент пересечения. Контроль пересечений производился путем последовательного поиска экстремумов

функции х (^) в окрестностях точек ^ = гП 1 - к2, вначале при г = —1, -2, -3,-А /п для определения типа траектории и, если пересечение оказывалось первым, для этой траектории подсчитывалось число пересечений по той же методике для г = 1, 2, 3,., А^ /п.

!ал

о

°0

18 Л

1Л

/ /£

л

.5 ' 4 ^ ", Ч N N ч . ч ,5 ’ 2

Ч Ч Ч N Ч V N X' Ч Ч X ' Ч \ 4 ч Ч ч Чч 'ч

ч V ■ч ч. . Ч* ч. Ч*

Я = 3

7,5

10

Рис. 6

В результате определялась как величина вероятности р так и среднее число пересечений в серии п*. Результаты расчетов приведены на рис. 6 для различных значений параметров R и р для

к = 10-2. Из приведенных результатов видно, что условие П* Р1 = 1 выполняется достаточно точно при всех значениях параметров, что объясняется в первую очередь протяженным интервалом для контроля пересечений.

Для оценки влияния схемы вычислений на точность определения величин п* и Р1 рассматривался случай р = 1, для которого спектральная плотность задавалась в виде

5 (ю) = - 4к

п

(ю2 - ,)’

-4к2ю2

так, что с х = с = 1.

Корреляционная

соотношением

Кх (*) = ке

-к

cos

(ю7) sin (ю7)

ю

функция определяется

(28)

где ю' = ^ 1 - к2 .

Расчеты проводились сериями для различных фиксированных значений производной у0, для которой задавалось 20 значений из диапазона (0 ^ 3с у ).

В этом примере значения х (^) при ^ < 0 и при ^ > 0 некоррелированы, и при заданном значении У0 процедуры выделения первых пересечений и определения среднего числа пересечений в серии проводятся независимо. Кроме того, процедура определения среднего числа пересечений может быть существенно упрощена. При малых значениях параметра к пересечения происходят

в окрестностях экстремумов математического ожидания и среднее число пересечений определяется как сумма вероятностей [3]:

-а®

где ^ — точки экстремумов функции

" Я " -Ф Я -\тх ()

_Сх _ _ сх (и) _

Я - К (^ )1

(29)

С х ( ^ )

соотношениями (14) и (15) по известной корреляционной функции (28).

В табл. 1 приведены результаты расчетов для различного количества точек т, в которых контролируется условие (23) для Я = 10.

Т аблица 1

к = 10-2

к = 10-3

к = 10-4

к = 10-5

т п* п*Р1 п* п*Р1 п* п *Р1 п* п*Р1

1 1,35 1,08 5,20 1,56 40,0 3,10 340 6,80

2 1,30 0,95 5,05 1,31 37,0 2,84 335 3,80

4 1,29 1,00 5,04 1,15 37,0 2,17 334 3,10

8 1,29 0,99 5,04 1,05 36,9 1,48 334 2,20

16 1,29 1,00 5,03 0,91 37,0 1,04 334 1,65

32 1,30 0,97 5,03 0,97 36,9 0,94 334 1,22

64 1,29 0,99 5,03 0,90 37,0 0,85 333 0,92

Из результатов табл. 1 видно, что при увеличении количества точек контроля уменьшается часть первых пересечений в серии, а произведение п*р стремится к 1, хотя при малых значениях

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

Это объясняется зависимостью вероятности р от числа точек контроля. При изменении числа точек контроля изменение величины вероятности может быть записано в виде р1 (Уо, т ) = р1 (Уо)с (т) и множитель С(т), исчезает при нормировке распределения для производной (27).

Данный результат позволяет просто определить распределение производной для первых пересечений из условия х(-п) < Хо:

/ (У0 ) =

Р(|х(-п)| < Я)Уо^■у0/2

|Р (х(-п) < Я) ус

Г уо/2

где Р (| х(-п) | < Я) = Ф

-Я - тх (-п) СТх (-п)

и для малых значений к —

Д-п) *\я2

сХ (-п) = ■! 1 - ехр (-2к п).

Как следует из результатов работы [3], при малых значениях параметра к среднее число

, - 1 1 ф(к, Я) п

пересечений в серии определяется формулой п* = 1 + п = 1 +----------—. Результаты расчетов данной

пкЯ 2

работы приведены на рис. 7 сплошными линиями для различных значений параметров Я и к. Из приведенных результатов видно, что незначительное отличие от универсальной зависимости наблюдается лишь при Я = 3. На этом же рисунке символами «*» отмечены результаты расчетов, в которых число пересечений определяется методом статистических испытаний (Монте-Карло). Видно, что суммирование вероятностей по формуле (29) дает меньшее (на ~10%) число пересечений. Для объяснения этого факта в табл. 2 приведено сравнение значений параметра п, вычисленных по формуле (29) и методом статистических испытаний по 1000 реализаций.

Данные табл. 2 показывают, что расхождение результатов не связано с определением распределения производной, а определяется предположением, на котором основана формула (29), о том, что у случайных реализаций точки экстремумов совпадают с точками экстремумов математического ожидания и, следовательно, подсчитывать пересечения можно только по точкам экстремумов. При конечных значениях Я это предположение выполняется неточно, что приводит к некоторому уменьшению числа пересечений в серии.

Работа выполнена при поддержке Российского фонда фундаментальных исследований (проект 01-01-00431).

ЛИТЕРАТУРА

1. Стратонович Р. Л. Избранные вопросы теории флуктуаций в радиотехнике.—

М.: Изд-во МГУ.— 1966.

2. Тихонов В. И. Выбросы случайных процессов.— М.: Наука.— 1970.

3. Кузьмин В. П., Ярошевский В. А. Оценка предельных отклонений фазовых координат динамической системы при случайных возмущениях. — М.: Наука, Физматлит.—

1995.

4. Я н к е Е., Э м д е Ф., Лёш Ф. Специальные функции.— М.: Наука.— 1977.

5. Корн Г., Корн Т. Справочник по математике.— М.: Наука.— 1974.

6. Королюк В. С., Портенко Н. И., Скороход А. В., Турбин А. В.

Справочник по теории вероятностей и математической статистике.— М.: Наука.— 1985.

Рукопись поступила 21/Х112001 г.