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

книги из ГПНТБ / Болдырев, В. С. Методы математической статистики в гидрографии и кораблевождении учеб. пособие

.pdf
Скачиваний:
7
Добавлен:
19.10.2023
Размер:
6.33 Mб
Скачать

■ é r ) K ( ^ ) - K C O )

А .

Сравнивая полученное выражение с (4 .5 ),

видим, что

эта формула дает приближенное значение интеграла,полу­ ченного из (4 .5) путем замены 7 на 7 и вычис­ ленного методом трапеций.Это обстоятельство позволяет применить формулы для остаточного члена при численном интегрировании.Обозначим через diJfXj изменение диспер­ сии оценки X вследствие замены непрерывной реализации дискретными отсчетами. Как известно, ошибка, возникаю­ щая при численном интегрировании методом трапеций с магом А , моает быть представлена так:

* " T r Z f V ‘ - y ) >

где

h °

. N

« в а т я ютервале

.0,1). Применяя эту формулу

к нашему случаю, получим

.2 т [2

-у)к(г)<і-Z . (4.10)

Предположим, что интервал выбран достаточно малым. Тогда сушу в (4.10) можно рассматривать как приближен­ ное значение интеграла от

Ы Ь -

2 0 7

взятого методом прямоугольников. То есть приближенно можно написать: Т+&

о

Т+й Г

* j t f c  h ^ - т Ц 1- т )к №

.

о

о

 

Разлагая правую часть полученного равенства по сте­

пеням Л и сохраняя

только линейные и квадратичные

члены, получим

^

 

â I ) [ x ] ^ ( A ä - 2 ß ü ) j T 1

(ч.П)

О

K ( * ) - - g r K ( T ) ■

Так как выражение (4 .II) является полиномом второй степени относительно Л , то при

 

 

Л

В

 

 

<Ул[Х]

^

^ m i n

А

 

 

ю»е«т минимальное

значение, равное

 

 

 

 

 

с * . і я

При достаточно

больном

Т

знак отноиения ~

опре

деляѳтся

знаком числителя.

Поэтому приЬ > 0 можно вы­

брать шаг дискретности Д - й опт , при котором

S D

2 0 8

будет

отрицательный и, следовательно, дисперсия оценки

X

,

полученной на

(4 .2 ),

будет меньше дисперсии,по­

лученной по (4 .4 ).

Результаты

расчетов,

выполненных

по

изложенной выше методике,

показывает,

что выигрыш

в точности определения^

уменьшается с

ростом Г , а

общий характер зависимости нормированной средней квад­ ратической ошибки оценки математического ожидания вид, представленный на рис.24.

Перейдем к нахождению оценки корреляционной функции X ( t ) • Случайную ф ункциях ^) пока не будем считать стационарной.

По определению,

n(t,t,)= /* {[ x ( t ,) - W ] [ х и ,) - X ( t , )]j >

поэтому для нахождения оценки K(jt~t2) мохно исходить

из тех же соображений, которые были положены в основу

14

209

наховдения X ( t ) • При наличии п

реализаций длиной Г

оценка К (t, t 2 ) может быть найдена по формуле

/='

ЕслиX ( t i ) неизвестно, то в предыдущей формуле вместо

него можно поставить его оценку и тогда

) - * ( & ] > -is )

1

где, как и обычно в подобных задачах, множитель — за-

менен на

У

п

 

для того, чтобы оценка была несмещенной.

п-1

Выполняя перемножение, получим

Г 1

Рассмотрим теперь методы нахожденія оценки корреля­ ционной функции стационарной случайной функции. Пусть в нашем распоряжении имеется одна реализацияx ( t ) дли­ тельностью (О 7 Т ) , а математическое ожидание X ( t ) из­

вестно. Теперь

т-т

K { Z ) = j ^ z [ [ * ( 0 - X ] [ x ( t + T ) - X ] c t t ,

 

а если X ( t ) - неизвестно, то

 

г-г

 

I [ c ( t ) - X ] [ x ( t + Z ) - X ] c / t .

(,4. Ь )

2ГО

Можно вести обработку опытного материала и иначе: сначала получить оценку математического ожидания произ­

ведения x ( t ) x ( t + T ) I а затем перейти к оценке корреля­

ционной функции, воспользовавшись известным соотношени­ ем между центральным и начальным моментами второго по­ рядка:

K( T ) = M[ x ( t ) x ( t + T ) ] -( X ) 2

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

 

T-Z

_

 

 

1

 

 

x(t)x(ti- Z)dt - (X ) 2

(4.16)

 

О

 

 

 

о

 

 

Формула (16)

не совпадает

с (15), однако при больших

Т результаты,

полученные

по этим формулам,

будут

практически одинаковыми. Проверим несмещенность оценки /С (Т ) . Вычислим математическое ожидание обеих частей

(4.16)«

м [к і т )]= к (т ) - я [х ] ■

Для эргодических случайных процессов (€ігг>К('і)=-О)

г ~

1

 

г -— о « .

 

дисперсия D[X

\

»

 

 

 

по доказанному ранее, стремится к

нулю при росте

Т .

Следовательно,

 

€ i m M [ K ( z j \ = K ( Z ) »

 

Отметим, что

 

в отличие от X математическое

ожида­

ние сценки К (Z) ра0НО К ( Т ) только в пределе.

Чтобы

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

нии X

используется

не

вся реализация, а только за­

пись за

время Т - Т

,

то есть

 

 

 

Г»ІІ

г-г

X- т-г

x ( t ) d t .

 

 

 

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

 

 

T-Z T-Z

 

 

 

К(Th(г- t r J

 

.

(4.17)

О О

 

 

 

Определяя дисперсно К ( 7 )

по общим правилам,

будем

иметь

,2

ф(т)]-/>і{[#(т)]2}-{м [Х(т )]]

Подставляя

сюда H ( Z ) из (4 .17),

получим

 

 

т-г

 

 

 

 

О

 

 

 

[*(?3 * 7)-*(U

 

dtlclt2alt3alt<f-

 

- { K ( z ) - M

\ - £

 

r ) K < a J t f

 

о

 

 

 

 

Полученная

формула показывает, что для нахождения

 

)]J

 

 

 

1)[К(Т)] надо знать

моменты до четвертого порядка

включительно. Если же процесс нормальный, то эти мо­ менты можно выразить через математическое ожидание и корреляционную функцию. Так как для любых четырех нор­ мальных величин X, і х 2 і х 3 ; х ^

212

M \ x t х г х 3 x ^ J

К)2 K3L/ + К/3

+Kfq Кгз +Xt Х 2

Кзі +

+Хг Х3 КЧі +

W

2^ X

z X 4 K/3+ \ X , K 12+ X ,X ,K 23+

+ X , X g X s X t .

 

 

 

 

Понимая под

х ,

, х 2

, х 3 ,

соответственно x ( t , ) ,

x ( t 3) y \ x ( t ! + T ) - x ( t 2) ]

И [ х ( ^ г ) - х ( ^ ) ]

. п о ­

лучим

 

 

 

 

 

D [K (r)]

(тг р[г/((т,)-к(г-г,)-/<(('Т,

 

 

V

 

 

 

 

 

О

 

 

 

 

г_

t

^ ( Т г Ъ ) [ к * ( Ъ ) +к ( т + Ъ ) к (? - т , ) \ ^

V

в

 

K ( t f t l)K(É3~t/ - T ) o t t /c/t2 clt3 ,

(4 .18)

J J

 

 

 

О

 

 

Т,= т-

т .

 

Если X известно

точно, то аналогичные рассуадения

дают

т-і

я[Х(г)]-^^{{T-z-r^^bKCT-zYiztz^dz,.

(4 .19)

213

Из формул (4.18) и (4.19) следует, что при наличии эргодичности нормального процесса x ( t )

-é’im D [Я ('Z ) \ = 0 у

Т ——с*0 1

т .е . оценка К ( Т ) является состоятельной и несмещенной. Вид формулы (4.19) показывает, что точность ординат

К ( Т )

с увеличением Z

падает. Происходит эго потому,

что интервал осреднения

функции x ( t ) x ( t

+ Т )

умень­

шается.

 

 

 

 

 

 

На практике обычно^

 

неизвестно наместо него,как

отмечалось

ранее,

используется оценка X

Кроме того,

перед

определением К(Т)

реализация центрируется,т.е.

 

 

 

 

 

г

 

 

 

 

 

 

 

О

 

 

 

 

 

/- С

 

 

 

 

 

 

о

 

 

 

 

Таким образом, можно записать, что

 

 

 

т-т

 

 

г

 

 

K ( T ) - y j ^ c ( t [ x ( t ) x ( t + T ) - j - x ( t + T > ( t ' ) c l t ' -

 

Т

О

 

Г

0 т

 

 

-yjx(t)x(t')cltl- уъ dt'

 

 

 

О

 

 

о

О

 

 

Проведем осреднение по многим реализациям!

_

_

2

СТ

 

і

J

X

Kx (Z)=K(Z) + j 2

J ( r - T ) K ( t ) c l t -

214

(

(J-Z-t)K(t-T)olt-

 

0fz T-T-t)K(t)dt+ \0 r

 

 

+ ^(T-?)K(t+T)dit+^ - і) K ( t ) d t ~

 

( r - t ) K ( t ) c l t +

\ (T~ t ) K ( t + T ) c / t 1

(*.20)

Из этой формулы следует, что даже

при большом числе

реализаций KX (Z) не

сходится кА'СГ)

и остается

зав и ся*.

щей от длины записи Т . При<£-~0 *з (*.20) получается выражение для дисперсии измеряемой величины в зависи­ мости от времени осреднения;

Kx ( 0 ) = K ( 0 ) - j t ( T - t ) K ( t ) c t t .

Рассмотрим особенности оценки корреляционной функп»и по дискретному числу ординат случайного процесса.Пусть

Х = 0 , а X {t) - нормальный процесс. Тогда

К(Z U — V —

т- £ +і

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

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

215

Кроме того, для различных Z величина Л

должна

быть различной, а расчеты для ее нахождения очень

громоздки. По указанным соображениям, обычно

стараются

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

T-Z

K(Z) =^ - ^ [ x ( t ) - x ] [ x ( t +Z )-X ]c lt -

О

Некоторые корреляторы вычисляют значение интеграла

Принципиальная схема таких приборов состоит в сле­ дующем. С подвижной ленты, на которой записана реализа­ ция, снимают значения ординат случайной функции, соот­ ветствующих двум моментам времени, сдвинутым относи­ тельно друг друга на интервал • Величину 2Г можно менять в соответствующих .пределах. Затем значения ор­ динат центрируются и перемножаются. Полученное произ­ ведение запоминается в сумматоре и складывается с ос­ тальными парными произведениями. Далее находится сред­ нее значение суммы, которое и будет соответствовать заданной величиной 'ZT , ординате корреляционной функ­ ции. Перечисленные действия могут выполняться с по-

216

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