Пример расчета и оформления результатов

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

Пример расчета и оформления результатов - student2.ru

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

Прежде всего, был проведен анализ однородности эмпирического распределения осадков по статистическим критериям Диксона и Смирнова - Граббса. Для этой цели ряд наблюдений был ранжирован и полученные значения приведены в табл.4. Эмпирическая обеспеченность ранжированного ряда (Р%) определялась по формуле:

Пример расчета и оформления результатов - student2.ru , ( 24 )

где m - порядковый номер членов ряда метеорологической характеристики, расположенного в убывающем порядке; n - общее число членов ряда.

Таблица 4.

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

Обеспеченность, P % Значение ранжированных осадков (мм) Год
1.370 57.0
2.740 52.0
4.110 46.0
5.479 44.0
6.849 44.0
8.219 40.0
9.589 38.0
10.959 37.0
12.329 37.0
13.699 35.0
15.068 35.0
16.438 33.0
17.808 31.0
19.178 29.0
20.548 29.0
21.918 29.0
23.288 28.0
24.658 26.0
26.027 25.0
27.397 23.0
28.767 22.0
30.137 21.0
31.507 20.0
32.877 20.0
34.247 20.0
35.616 20.0
36.986 20.0
38.356 19.0
39.726 19.0
41.096 19.0
42.466 18.0
43.836 18.0
45.205 17.0
46.575 17.0
47.945 17.0
49.315 17.0
50.685 16.0
52.055 16.0
53.425 15.0
54.795 15.0
56.164 15.0
57.534 13.0
58.904 13.0
60.274 13.0
61.644 13.0
63.014 13.0
64.384 13.0
65.753 13.0
67.123 13.0
68.493 12.0
69.863 12.0
71.233 12.0
72.603 10.0
73.973 10.0
75.342 10.0
76.712 9.00
78.082 9.00
79.452 9.00
80.822 9.00
82.192 8.00
83.562 8.00
84.932 8.00
86.301 7.00
87.671 7.00
89.041 7.00
90.411 7.00
91.781 6.00
93.151 6.00
94.521 4.00
95.890 4.00
97.260 3.00
98.630 1.00

Эмпирическое распределение осадков за ноябрь по метеостанции Hовый Порт приведено на рис.2.

Пример расчета и оформления результатов - student2.ru .

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

По ряду наблюдений сумм осадков за ноябрь на метеостанции Новый Порт за период 1924-2003 гг. были получены следующие значения параметров: среднее ряда = 19.2 мм, Cv=0.6424, Cs = 1.0804, Cs/Cv = 1.6818, r(1)= 0.3667. При объеме выборки 72 года (для автокорреляции n=71 год) и числе степеней свободы 69 (n-2) автокорреляция, как следует из формулы (18) и табл.1, является статистически значимой величиной. По формулам (1)-(12) были вычислены расчетные значения статистик Диксона и Смирнова-Граббса, которые приведены в табл.5. Критические значения были получены из таблиц приложения при имеющемся объеме ряда наблюдений (72 года), вычисленным значения коэффициентов асимметрии и автокорреляции и принятом уровне значимости 5%. Полученные критические значения статистик критериев также приведены в табл.5.

Таблица 5

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

Экстремум Критерий Расчетн.знач. Критич.знач. Уровень значимости расч. Вывод
max Диксон 1 0.0893 0.3268 11.0000 однороден
max Диксон 2 0.0926 0.3340 11.0000 однороден
max Диксон 3 0.2037 0.4102 11.0000 однороден
max Диксон 4 0.2075 0.4068 11.0000 однороден
max Диксон 5 0.1964 0.4041 11.0000 однороден
min Диксон 1 0.0357 0.0629 11.0000 однороден
min Диксон 2 0.0392 0.0729 11.0000 однороден
min Диксон 3 0.0588 0.0913 11.0000 Однороден
min Диксон 4 0.0667 0.1006 11.0000 Однороден
min Диксон 5 0.0536 0.0824 11.0000 Однороден
max Cмирнов-Граббс 3.0692 4.3010 11.0000 Однороден
min Cмирнов-Граббс 1.4754 1.8608 11.0000 Однороден

В табл.5 приведены названия критериев, расчетные значения статистик, их критические значения и уровень значимости, соответствующий расчетным значениям, который также определяется интерполяцией из таблиц критических значений Приложения. В данном случае при проверке экстремумов уровень значимости составляет более 10% и поэтому в таблице приведено значение 11%, свидетельствующее о том, что гипотеза об однородности принимается с высокой степенью достоверности. В остальных случаях, если уровень значимости, соответствующий расчетному значению статистки, находится в интервале от 1% до 10% - приводится его численное значение. Это позволяет дифференцировать вывод о принятии или отклонении нулевой гипотезы об однородности. Так, если уровень значимости близок к 5% (находится в диапазоне от 4% до 5%), то гипотеза принимается, но с учетом этого обстоятельства и вывод помещается в скобки: (+).

Для оценки стационарности дисперсий и средних значений многолетний ряд наблюдений вначале был разбит на две равные части или подвыборки: первая с 1924 по 1964 гг., вторая с 1965 по 2003 гг. Для каждой подвыборки были рассчитаны средние значения и дисперсии, которые соответственно равны: YcpI.=13.7 мм, σ 2I. =115.98 мм2 и YcpII.=24.6 мм, σ 2II. =130.75 мм2. По полученным дисперсиям и средним значениям были определены расчетные значения статистик Фишера и Стьюдента (формулы (19) и (22)), которые приведены в табл.6. На основе вычисленных значений новых степеней свободы (формулы (20) и (21), таблицы 2 и 3) и таблиц приложения были определены критические значения статистик Фишера и Стьюдента, которые также даны в табл.6. Из результатов табл.6 следует, что расчетная статистика критерия Фишера меньше критической, поэтому вывод о стационарности дисперсий принимается. Вместе с тем, расчетное значение статистики Стьюдента значительно превышает его критическое значение и гипотеза о стационарности средних значений двух частей временного ряда осадков отклоняется с высокой степенью достоверности.

Таблица 6.

Оценка стационарности средний значений и дисперсий ряда наблюдений сумм осадков за ноябрь на метеостанции Новый Порт.

Критерий Расчетное зн. Критич.зн. Уровень значимости расч. Вывод
Критерий Фишера 1.1274 1.9141 11.0000 однороден
Критерий Стъюдента 4.1116 2.8580 0.9000 неоднороден

Вторым способом временной ряд осадков, приведенный на рис.1, был разделен по дате разрыва данных наблюдений, которые закончились в 1955 г. и были возобновлены в 1961 г. Поэтому ряд наблюдений был разбит на следующие две части: 1924 -1955 гг. и 1961 – 2003 г. Для каждой подвыборки были рассчитаны средние значения и дисперсии, которые соответственно равны: YcpI.=10.7 мм, σ 2I. =25.6 мм2 и YcpII.=26.0 мм, σ 2II. =149.7 мм2. Как следует из сравнения параметров, полученных первым и вторым способами, во втором случае как существенно отличаются средние значения (в 2.6 раза), так уже и дисперсии (в 5.8 раза). Полученные расчетные и критические значения статистик критериев Фишера и Стьюдента при втором способе разбиения ряда приведены в Табл.7 и свидетельствуют о существенной нестационарности как средних значений, так и дисперсий.

Таблица 7.

Оценка стационарности средний значений и дисперсий ряда наблюдений сумм осадков за ноябрь на метеостанции Новый Порт при разбиении ряда по датам производства наблюдений.

Критерий Расчетное зн. Критич.зн. Уровень значимости расч. Вывод
Критерий Фишера 5.9457 1.9494 0.9000 неоднороден
Критерий Стъюдента 7.2174 2.8727 0.9000 неоднороден

Оценка стационарности осадков были проведена и для остальных месяцев года и результаты анализа помимо индивидуальных таблиц расчетов для каждого месяца (Табл.5- 7) были представлены в виде одной обобщенной таблицы 8. В таблице 8 знаком «+» представлен вывод о принятии гипотезы однородности и стационарности, а знаком «-» - об отклонении. Если расчетное значение статистики критерия близко к критическому, но превышает его, то в таблице также приводится уровень значимости, соответствующий расчетному значению критерия. В этом случае вывод о принятии или отклонении гипотезы однородности и стационарности является сомнительным и заключается в скобки.

Таблица 8.

Результаты оценки стационарности месячных сумм осадков по критериям

Диксона (Д) , Фишера (Ф) и Стьюдента (Ст) на метеостанции Новый Порт

Д Ф Ст Д Ф Ст Д Ф Ст
январь февраль Март
+ + - + + + + - -
апрель Май Июнь
+ (-)3.0 - + - + + + -
июль август Сентябрь
+ + + + + + + - +
октябрь Ноябрь Декабрь
+ + - + + - + (+)4.4 -

Анализ табл.8 позволяет сделать выводы, что все распределения сумм осадков по месяцам являются однородными. Установлена нестационарность дисперсий и особенно средних в холодные месяцы года. Полученная нестационарность вполне может быть обусловлена сменой приборов для регистрации осадков как раз в 1960е годы, причем предшествующий осадкомер не позволял нейтрализовать ветровое выдувание осадков, что особенно характерно для твердых осадков и для станций, находятся на морском побережье в зоне больших ветров.

Наши рекомендации