Анализ периодических компонент интенсивности точечных процессов. Любушин А.А. - 13 стр.

UptoLike

ВУЗ: 

РГГРУ | Москва

Составители: 

Рубрика: 

13
max min
min
lg( ) lg( )
10 , lg( ) ( 1) , 1,...,
(1)
j
jj p
p
TT
TTj j
N
ϕ
ϕ
N
−
==+−⋅ =
−
(41)
4.
Задать режим вычислений – по всей выборке или в скользящем временном окне.
5. Если задан режим оценки в скользящем окне, то необходимо ответить на
дополнительные вопросы:
5.1.
Задать длину скользящего временного окна в размерном времени.
L
T
5.2.
Задать смещение скользящих временных окон в размерном времени.
L
TΔ
5.3.
Задать смещение
s
hift
T временных меток в выводном grd-файле.
Если задан режим оценки по всей выборке, то в текущей директории создается выводной
файл с именем "
PPP_out.dat", который представляет собой таблицу из 3-х колонок длиной
строк. В первой колонке последовательно идут тестированные значения периодов , во
второй – разность между максимумами логарифмических функций правдоподобия, а
в 3-й – значение безразмерного параметра
p
N
j
T
ln( )LΔ
,0 1aa
<
≤ амплитуды гармонического колебания
интенсивности процесса в модели (31).
Если задан режим оценки в скользящем временном окне то, в текущей директории
создается выводной файл с именем "
PPP_out.grd" в символьном grd-формате, готовый для
построения 2-мерных или 3-мерных рельефов частотно-временных диаграмм для статистики
в пакете Surfer. При этом временные метки соответствует значениям
ln( )LΔ
( 1) , 1,...,
s
hift L L W
TTk Tk N++−⋅Δ = ,
где индекс нумерует временные окна, - общее число временных окон. Таким образом,
временные метки соответствуют правым концам скользящих временных окон, взятых со
смещением
k
W
N
s
hift
T . Параметр
s
hift
T может быть полезен, например, в ситуации, когда времена
событий указываются в годах от начала 1900-го года, а на диаграмме полезно иметь метки в
общепринятых годах – тогда задается
s
hift
T =1900.
6. Примеры применения.
Во всех рассматриваемых ниже примерах, если они касаются анализа последовательности
землетрясений, данные взяты из источника:
http://neic.usgs.gov/neis/epic/ .
6.1. Периодические компоненты интенсивности сильнейших землетрясений для всего
мира. На рис.1. представлена последовательность сильнейших сейсмических событий для
глобального сейсмического процесса с начала эпохи инструментальной сейсмологии по 
                                                       13
                               ϕ                                  lg(Tmax ) − lg(Tmin )
                      T j = 10 j , ϕ j = lg(Tmin ) + ( j − 1) ⋅                         ,   j = 1,..., N p   (41)
                                                                       ( N p − 1)


   4. Задать режим вычислений – по всей выборке или в скользящем временном окне.
   5. Если задан режим оценки в скользящем окне, то необходимо ответить на
       дополнительные вопросы:
       5.1. Задать длину TL скользящего временного окна в размерном времени.

       5.2. Задать смещение ΔTL скользящих временных окон в размерном времени.

       5.3. Задать смещение Tshift временных меток в выводном grd-файле.

  Если задан режим оценки по всей выборке, то в текущей директории создается выводной
файл с именем "PPP_out.dat", который представляет собой таблицу из 3-х колонок длиной
N p строк. В первой колонке последовательно идут тестированные значения периодов T j , во

второй – разность Δ ln( L) между максимумами логарифмических функций правдоподобия, а
в 3-й – значение безразмерного параметра a, 0 < a ≤ 1 амплитуды гармонического колебания
интенсивности процесса в модели (31).
  Если задан режим оценки в скользящем временном окне то, в текущей директории
создается выводной файл с именем "PPP_out.grd" в символьном grd-формате, готовый для
построения 2-мерных или 3-мерных рельефов частотно-временных диаграмм для статистики
Δ ln( L) в пакете Surfer. При этом временные метки соответствует значениям


                                    Tshift + TL + (k − 1) ⋅ ΔTL , k = 1,..., NW ,


где индекс k нумерует временные окна, NW - общее число временных окон. Таким образом,
временные метки соответствуют правым концам скользящих временных окон, взятых со
смещением Tshift . Параметр Tshift может быть полезен, например, в ситуации, когда времена

событий указываются в годах от начала 1900-го года, а на диаграмме полезно иметь метки в
общепринятых годах – тогда задается Tshift =1900.

  6. Примеры применения.
  Во всех рассматриваемых ниже примерах, если они касаются анализа последовательности
землетрясений, данные взяты из источника: http://neic.usgs.gov/neis/epic/ .
  6.1. Периодические компоненты интенсивности сильнейших землетрясений для всего
мира. На рис.1. представлена последовательность сильнейших сейсмических событий для
глобального сейсмического процесса с начала эпохи инструментальной сейсмологии по

Страницы