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

Климатология лабы / Климатология лаб1 готовая - копия

.docx
Скачиваний:
15
Добавлен:
15.03.2019
Размер:
115.45 Кб
Скачать

Министерство образования и науки Российской Федерации

Федеральное государственное бюджетное образовательное учреждение высшего образования

Российский государственный гидрометеорологический университет

(РГГМУ)

Кафедра МКОА

Лабораторная работа

по дисциплине «Климатология»

по теме: «Оценка однородности и стационарности»

Выполнил: ст. гр. ПМ-Б15-1-3

Скачков И.В.

Проверил: к.г.н., доцент каф. МКОА

Шадурский А.Е.

Санкт-Петербург 2017

Задачи работы:

1. Расчет основных параметров рядов наблюдений (среднее значение, дисперсия, асимметрия, автокорреляция) и их эмпирических распределений (ранжирование ряда и расчет эмпирической обеспеченности).

2. Оценка случайных погрешностей и статистической значимости параметров распределений.

3. Оценка однородности интегральных эмпирических распределений по статистическим критериям Диксона и Смирнова – Граббса на резко отклоняющиеся экстремумы.

4. Оценка стационарности средних значений и дисперсий двух частей временного ряда (при разбиении ряда на две равные части и при произвольном разбиении ряда) по критериям Стьюдента и Фишера.

Исходные данные

Многолетние ряды наблюдения среднемесячных температур и осадков в п. Илирней, Чукотского А.О.(67,3 с.ш 168 в.д) (25248).

Обработка результатов для осадков.

Рассчитанные значения эмпирической обеспеченности и ряд ранжированых осадков и температур за ноябрь на метеостанции Илирней.

Хронологический график сумм осадков и температур за январь по метеостанции Илирней.

Pm,%

значение ранжир.осадков

Pm,%

значение ранжир.темп

1,92

57

1,92

-13,9

3,84

53

3,84

-19,7

5,76

48

5,76

-23,9

7,69

46

7,69

-24,2

9,61

43

9,61

-24,4

11,53

36

11,53

-25,6

13,46

33

13,46

-27

15,38

31

15,38

-27

17,30

30

17,30

-27,4

19,23

27

19,23

-27,6

21,15

22

21,15

-28,1

23,07

21

23,07

-28,2

25

21

25

-28,6

26,92

20

26,92

-30,4

28,84

20

28,84

-30,6

30,76

17

30,76

-31,6

32,69

17

32,69

-31,6

34,61

17

34,61

-31,7

36,53

16

36,53

-31,9

38,46

14

38,46

-32,1

40,38

14

40,38

-32,2

42,30

14

42,30

-32,2

44,23

14

44,23

-32,3

46,15

14

46,15

-32,4

48,07

13

48,07

-32,9

50

12

50

-33

51,92

12

51,92

-33,4

53,84

12

53,84

-33,8

55,76

12

55,76

-33,8

57,69

11

57,69

-34,1

59,61

11

59,61

-34,3

61,53

11

61,53

-34,3

63,46

11

63,46

-34,4

65,38

11

65,38

-35

67,30

10

67,30

-35,1

69,23

10

69,23

-36,1

71,15

9

71,15

-36,1

73,07

8

73,07

-36,2

75

8

75

-36,2

76,92

8

76,92

-36,9

78,84

8

78,84

-37,2

80,76

7

80,76

-37,8

82,69

7

82,69

-38,1

84,61

6

84,61

-38,3

86,53

6

86,53

-38,9

88,46

6

88,46

-39,1

90,38

3

90,38

-40,3

92,30

3

92,30

-41,4

94,23

2

94,23

-42,3

96,15

2

96,15

-43,7

98,07

2

98,07

-44

Обработка результатов для осадков.

Рис.1. Хронологический график сумм осадков за январь по метеостанции Илирней.

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

Из графика можно сделать вывод что резких экстремумов нет.

По ряду наблюдений сумм осадков за январь по метеостанции Илирней за период с 1944-2004 гг. были получены следующие значения параметров: Cs=0,34299717 r(1)= 0,17628357. По формулам (1)-(12) были вычислены расчетные значения статистик Диксона и Смирного-Граббса, которые приведены в табл 1. Критические значения были получены из таблиц приложения при имеющемся объеме ряда наблюдений (51 год), вычисленным значениям коэффициентов асимметрии и автокорреляции и принятом уровне значимости 5 %. Полученные критические значения статистик критериев.

Таблица1.Результаты оценки однородности эмпирического распределения сумм осадков за январь на метеостанции Илирней по критериям Диксона и Смирного-Граббса.

Экстремум

Критерий

Расчетное значение

Критическое значение

Вывод

max

Диксон1

0,0727

0.31

Однороден

max

Диксон 2

0,0727

0.31

Однороден

max

Диксон 3

0,0636

0.30

Однороден

max

Диксон 4

0,0636

0.30

Однороден

max

Диксон5

0,0636

0.30

Однороден

min

Диксон1

0,0181

0.04

Однороден

min

Диксон 2

0,00196

0.16

Однороден

min

Диксон 3

0,00196

0.16

Однороден

min

Диксон4

0,00196

0.16

Однороден

min

Диксон5

0,00196

0.16

Однороден

max

Смирнов-Граббс

7,2111

8,59

Однороден

max

Смирнов-Граббс

2,6993

3.19

Однороден

Вывод: расчётные значения критериев Диксона и Смирнова-Граббса ниже критических значений, при уровне значимости 5%, из этого делается вывод, что ряд наблюдений однороден.

Для оценки стационарности дисперсий и средних значений многолетний ряд наблюдений вначале был разбит на две части: первая – с 1945 по 2002г., вторая – с 1970 по 2001 г. Для каждой подвыборки были рассчитаны средние значения и дисперсии, которые соответственно равны: Yср1. = 20,4 мм, Дисп.1= 257,36 мм^2 и Yср2. = 13,92мм, Дисп.2= 77,91 мм^2. По полученным дисперсиям и средним значениям были определены расчетные значения статистик Фишера и Стьюдента (формулы (19) и (22)), которые приведены в таблице ,на основе вычисленных значений новых степеней свободы(формулы (20) и (21),таблица 2 и 3 из приложения) и таблиц приложения были определены критические значения статистик Фишера и Стьюдента ,которые также даны в табл.3.

Таблица 2. Анализ параметров ряда наблюдений для осадков.

r1

t

t кр.

Вывод

0,17

1,25

3,38

Нулевая гипотеза отклоняется

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

Критерий

Расчетное значение

Критическое значение

Вывод

Критерий Фишера

1,74

2,45

однороден

Критерий Стьюдента

2,11

0,77

неоднороден

Обработка результатов для температур.

Рис 1. Хронологическое распределение температур за январь на метеостанции Илирней.

Рис 2.Эмпирическое распределение температур за январь на по метеостанции Илирней.

Из графиков можно сделать вывод что резких экстремумов нет.

Таблица1. Результаты оценки однородности эмпирического распределения сумм температур за январь на метеостанции Илирней по критериям Диксона и Смирного-Граббса.

Экстремум

Критерий

Расчетное значение

Критическое значение

Вывод

max

Диксон1

0,19269

0.34

Однороден

max

Диксон 2

0,19463

0.33

Однороден

max

Диксон 3

0,33557

0.44

Однороден

max

Диксон 4

0,35211

0.45

Однороден

max

Диксон5

0,33222

0.44

Однороден

min

Диксон1

0,00996

0.01

Однороден

min

Диксон 2

0,00996

0.01

Однороден

min

Диксон 3

0,33557

0.44

Однороден

min

Диксон4

0,33557

0.44

Однороден

min

Диксон5

0,33557

0.44

Однороден

max

Смирнов-Граббс

7,2111

8,54

Однороден

min

Смирнов-Граббс

4,2911

6.15

Однороден

Вывод: расчётные значения критериев Диксона и Смирнова-Граббса ниже критических значений, при уровне значимости 5%, из этого делается вывод, что ряд наблюдений однороден.

Для оценки стационарности дисперсий и средних значений многолетний ряд наблюдений вначале был разбит на две части: первая – с 1945 по 2002г., вторая – с 1970 по 2001 г. Для каждой подвыборки были рассчитаны средние значения и дисперсии, которые соответственно равны: Yср1. = 20,4 мм, Дисп.1= 257,36 мм^2 и Yср2. = 13,92мм, Дисп.2= 77,91 мм^2. По полученным дисперсиям и средним значениям были определены расчетные значения статистик Фишера и Стьюдента (формулы (19) и (22)), которые приведены в таблице ,на основе вычисленных значений новых степеней свободы(формулы (20) и (21),таблица 2 и 3 из приложения) и таблиц приложения были определены критические значения статистик Фишера и Стьюдента ,которые также даны в табл.3.

Таблица 2. Анализ параметров ряда наблюдений для температур за январь на метеостанции Илирней..

r1

t

t кр.

Вывод

0,85

2,11

3,19

Нулевая гипотеза отклоняется

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

Критерий

Расчетное значение

Критическое значение

Вывод

Критерий Фишера

1,34

1,22

однороден

Критерий Стьюдента

1,17

0,69

неоднороден

Вывод: расчётные значения критериев Фишера и Стьюдента ниже критических значений, при уровне значимости 5%, из этого делается вывод, что ряд наблюдений однороден.

Вывод:

Цели лабораторной работы достигнуты полностью. Ряды наблюдений сумм осадков в температуры на метеостанции Илирней были исследованы статистическими методами. Критерии Диксона и Смирнова-Граббса выявляют резко отклоняющиеся экстремумы, которые появляются при аномально высоких отклонениях в ряде. За весь исследуемый период таких значений не наблюдалось. За лабораторную работу были рассчитаны средние дисперсии и следовательно были рассчитаны коэффициенты Стьюдента и Фишера. В конечном результате:

Ряд

Gn

G1

Диксон

Фишер

Стьюдент

Температура

+

+

+

+

+

Осадки

+

+

+

+

+