Файл: Григоркина Р.Г. Прикладные методы корреляционного и спектрального анализа крупномасштабных океанологических процессов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.06.2024

Просмотров: 234

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

•При этом следует иметь в виду, что условия стационарности в реальных процессах могут выполняться только с тем или иным приближением, так ка к физические процессы, в том числе и про­ цессы в океане, всегда имеют некоторую «динамику» развития. Использование теории нестационарных случайных функций в прикладных методах анализа сталкивается со значительными трудностями, обусловленными сложностью и громоздкостью ее приложения . Возникает проблема использования аппарата тео­ рии стационарных случайных процессов дл я описания заведомо нестационарных процессов. Это возможно для определенных

классов

нестационарных процессов,

которые

могут быть

выде­

л е н ы по

энергии

и мощности

этих

процессов

(Малахов,

1968).

П е р в а я группа

включает

процессы, полная

энергия которых

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

 

 

 

со

 

 

 

 

 

 

Ех= j x 2 ( 0 d i < o o .

 

(1.1)

 

 

—оо

 

 

 

iKo второй группе относятся процессы, у которых энергия бес­ конечна, а мощность — конечная величина

S x = l i m - £ = - J x a ( * ) < f t .

(1.2)

Это процессы, ординаты которых не о б р а щ а ю т с я в тождествен­

ный

нуль при

безграничном возрастании

продолжительности

процесса. К этой

группе относятся

п р е ж д е всего все

стационар­

ные

процессы, периодические процессы, в том числе

гармониче­

ские колебания, постоянная величина и т. п.

 

 

К

третьей

группе

п р и н а д л е ж а т

процессы,

не удовлетворяю ­

щие

условиям

(1.1) и

(1J2) (импульсные процессы, быстро зату­

х а ю щ и е процессы и другие нестационарные

процессы) .

Н е с м о т р я

на то, что процессы первой группы нестационарны,

а вторая группа

включает ка к стационарные, та к и нестационар­

ные процессы, корреляционные функции процессов обеих групп о б л а д а ю т свойствами корреляционных функций стационарных случайных процессов. В частности, косинус-преобразование Фу­ рье от корреляционных функций дл я процессов первой группы

дает

спектральную плотность

энергии,

а д л я процессов

второй

группы — спектральную

плотность мощности. Это

справедливо

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

процессов.

Таким образом, дл я процессов первой

и второй групп существу­

ет независимая

от времени

функция

спектральной

плотности.

Д л я

процессов

третьей

группы м о ж н о

поставить т а к ж е

вопрос

о существовании спектральной плотности дл я тех частот, на ко­

торых мощность процесса отлична от

бесконечности (Малахов,

1968).

 

 

 

Таким образом, р а с с м а т р и в а е м ы е

в

книге

прикладные мето­

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

анализа

стационарных про-

11


цессов могут быть т а к ж е использованы при исследовании неко­ торых классов нестационарных процессов, в том числе детерми­ нированных, гармонических, полигармонических и квазиперио­ дических процессов (см. гл. I I I ) .

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

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

ральных вопросов вероятностного анализа данных

наблюдений .

Этой проблеме посвящена

обширная литература .

 

 

 

Оценки д о л ж н ы удовлетворять

определенным

требованиям .

1. Оценка а г

д о л ж н а

быть состоятельной, т. е. д о л ж н а

схо­

диться по вероятности Р к оцениваемой величине

а при

неогра­

ниченном увеличении длины анализируемой реализации

Ти

 

 

l i m P ( | a r - f l | < e ) = l ,

 

(1-3)

где е —• наперед з а д а н н а я м а л а я

величина.

 

 

 

2. О ц е н к а д о л ж н а быть несмещенной. Несмещенной

оценкой

называют такую оценку, математическое ожидание

которой

р а в ­

но оцениваемой

величине

 

 

 

 

 

 

 

М(ат)=а.

 

 

(1.4>

Если это условие достигается только при Тп—>-оо,

то оценка

н а ­

зывается асимптотически

несмещенной

 

 

 

 

l i m M ( a r ) = a .

 

(1-5)>

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

12


тивность

оценки.

Если

при увеличении длины реализации

эффек ­

тивность

оценки

стремится к единице, то

т а к а я оценка

называ ­

ется асимптотически

эффективной. Д л я

случайных

процессов

качество оценки может

быть определено не только по

отношению

к оптимальной дисперсии, но и по отношению к дисперсии само­ го процесса. Отношение дисперсии оценки к дисперсии процесса называют показателем эффективности. Чем меньше этот показа­ тель, тем меньше дисперсия и тем эффективнее оценка.

Когда выборочное распределение оценки известно, м о ж н о найти доверительный интервал, в котором с заданной степенью достоверности '(доверительной вероятности) будет находиться оценка малого параметра . М е т о д ы определения доверительных интервалов и проверки статистических гипотез в большинстве случаев основаны на предположении о том, что реализации, из которых извлекаются выборки, подчиняются нормальному зако ­ ну. Это ограничение часто оказывается весьма существенным. З а к о н ы распределения океанологических характеристик обычно отличаются от нормального. Выборочное распределение самих оценок статистических характеристик в океане изучено еще очень

мало . Исключение, вероятно, составляют исследования ветрово­

го волнения. В связи с этим перспективными представляются

ме­

тоды оценивания, не з а в и с я щ и е от формы распределения,

т а к

называемые непараметрические методы.

 

§ 2. Автокорреляционная функция как характеристика внутренней структуры процесса

во временной области

К а к

отмечалось

в § 1, из-за

трудностей определения много­

мерных

законов

распределения

океанологических процессов

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

центральной

смешанной моментной

функции

второго

порядка .

Н а ч а л ь н а я моментная

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

мате­

матическим

ожиданием

случайного

процесса

 

 

 

 

 

 

mi(t)=M[X(t)],

 

 

 

(2.1)

где М — символ

математического ожидания,

X(t)—случайный

процесс.

 

 

 

 

 

 

 

Ц е н т р а л ь н а я

моментная функция

второго

порядка,

или

дис­

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

вокруг его математического о ж и д а н и я

и определяется следую­

щим

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

 

 

 

D[X(t)]=M[X{t)-mi(t)]\

(2.2)

где

D[X(t)]

— дисперсия случайного

процесса.

13


Оценки математического о ж и д а н и я и дисперсии широко при­ меняются при описании свойств океанологических процессов. В частности, оценки математического ожидания океанологиче­ ских характеристик используются при составлении атласов, на­

вигационных

пособий

и т.

д. Однако математическое

ожидание

и дисперсия

не содержат

информации о внутренней

структуре

процесса. Процессы

могут

иметь одинаковые математические

о ж и д а н и я и дисперсии, но различную внутреннюю структуру и различный характер протекания во времени, обусловленные кон­ кретными особенностями р е ж и м а отдельных районов Мирового океана. Например, в динамически активных районах темпера­ турные условия быстро изменяются и могут быть не связаны с предысторией.

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

ti)=M{[Xit1)-ml{tl)][X{t2)-ml{t2)]}. (2.3)

Автокорреляционная функция и математическое ожидание полностью определяют статистические закономерности случай­

ного стационарного процесса

в том случае, если распределение

ординат этого процесса подчиняется нормальному закону.

З н а я

оценку автокорреляционной

функции,

можно

решить широкий

круг задач, связанных с исследованием

и прогнозированием

из­

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

условий. >К числу

этих" з а д а ч отно­

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

Решение этих з а д а ч играет существенную роль в разработке физически обоснованных и статистически оправданных методов прогноза. Удачные опыты прогноза океанологических условий с помощью методов линейной экстраполяции случайных функций (Беляев, Болдырев, 1963; Овсянникова, 1967; Дмитриева, 1968) позволяют надеяться, что эти методы получат широкое распро­ странение. Это приводит к необходимости подробного исследова­ ния автокорреляционных функций различных океанологических процессов.

Д л я стационарного эргодического случайного процесса оцен­ ка автокорреляционной функции находится осреднением парных произведений ординат процесса

14


 

 

 

/?*(т) =

- r ^ r J

x(t)x(t+r)dt,

 

 

(2.4)

где

(т) —

оценка

истинного

значения

автокорреляционной

функции;

x(t)

—центрированная

 

реализация

процесса;

[О, Тп\

интервал

времени, на котором

з а д а н а реализация процесса (про­

должительность наблюдении);

 

т — временной сдвиг,

изменяю ­

щийся

в

пределах от

0

до

т т ;

т,„ — м а к с и м а л ь н а я

величина

ин­

тервала

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

автокор ­

реляции

'(максимальный

сдвиг

автокорреляционной

функции) .

В

случае

дискретизации исходной реализации

(квантовании

по времени) вычисление автокорреляционной функции осущест­

вляется по

формуле

 

 

 

 

 

 

 

 

 

 

 

 

R i < ш )

= - ^ ^ а м )

Х [ а + 1 ) м ] ,

 

(2.5)

где x(jAt).

х[(/+/)At]

— члены

последовательности

наблюде ­

ний;

jAt—t;

(j-\-l)At—t-\-x\

At

интервал

дискретности;

N

 

 

 

 

 

 

At=

Т

/ = 0 , 1 , • • •, m,

т

число

членов последовательности;

—гг-;

 

 

 

 

 

 

 

Nj

 

At(N—t)=Tn—x.

число

определяемых

ординат

оценки

Ш = т ;

Функции

(2.4) и

(2.5)

называют т а к ж е

автоковариационны ­

ми. Как следует из

(2.4), н а ч а л ь н а я

ордината

оценки

автокор­

реляционной

функции при

т = 0

соответствует

оценке

дисперсии

процесса

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г"

 

 

 

 

 

 

 

D*=d%

= R*(0)

=

-^-)^(t)dt.

 

 

 

(2.6)

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

У к а ж е м на основные свойства автокорреляционной функции.

15