Файл: Григоркина Р.Г. Прикладные методы корреляционного и спектрального анализа крупномасштабных океанологических процессов.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