Последовательность расчета

1. Выбираются два многолетних ряда наблюдений, один из которых имеет меньшую продолжительность наблюдений (n), а другой большую – N. Этими рядами могут быть следующие:

- ряды сумм осадков за один и тот же месяц на двух разных станциях;

- ряды среднемесячных температур воздуха за один и тот же месяц на двух разных станциях;

- ряд температур и осадков на одной станции за один и тот же месяц;

- два ряда температур или осадков за соседние месяцы.

2. Определяются коэффициенты уравнения связи между коротким (Y) и продолжительным (X) рядами наблюдений за совместный период:

Y= k0+k1X ( 30 )

где:

k1 = (∑ (Yi - Ycp )(Xi - Xcp ))/(∑(Xi - Xcp )2 , ( 31 )

k0 = Ycp - k1Xcp. ( 32 )

3. Рассчитываются коэффициенты и параметры, входящие в неравенство (28):

_____________________

R = (∑ (Yi - Ycp )(Xi - Xcp ))/√( ∑(Yi - Ycp )2∑(Xi –Xcp )2) , ( 33 )

_____ _____

sk= sY √(1-R2)/( sX √(n-1)) ( 34 )

4. Оценивается эффективность уравнения (30) по условиям (28) и если оно эффективно, то на основе этого уравнения восстанавливаются значения ряда Y за период N-n.

5. В соответствии с условиями (28), из восстановленных данных исключаются те значения относительная погрешность которых превосходит предельную заданную.

6. В восстановленные значения вносится поправка на смещенность дисперсии по формуле (29).

7. Рассчитываются показатели эффективности восстановления:

- абсолютное количество восстановленных лет ∆n;

- относительное количество восстановленных лет Δn’(%) ;

- критерий Фишера для оценки однородности дисперсий восстановленных и наблюденных данных;

- критерий Стьюдента для оценки однородности средних восстановленных и наблюденных данных.

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

Выбраны многолетние ряды сумм осадков за июль на двух метеостанциях:

- короткий ряд наблюдений - метеостанция Мыс Каменный (1951 -1994 гг.) n=44 года;

- продолжительный ряд-аналог – метеостанция Новый Порт (1925 – 2004 гг.) N=80 лет.

Прежде всего, был определен коэффициент парной корреляции R между двумя рядами за совместный период, который был равен 0.75, относительная погрешность восстановления задана sε/Yp ≤40%, условие k/sk ³ 2 также было выполнено и k/sk =3.1.

Полученное уравнение имеет вид:

Y= 0.6382X+ 15.7868 ( 35 )

при стандартной погрешности расчета sε=19.5 мм и относительных погрешностей восстановления от 20% до 40%.

На основе полученного уравнения (35) были восстановлены данные наблюдений за 36 лет и из восстановленных данных выбраны те, которые удовлетворяли условиюsε/Yp ≤40%. Таких восстановленных значений оказалось 16 и они были пересчитаны по формуле (29), чтобы исключить систематическое преуменьшение дисперсии. Результатом восстановления является ряд, включающий измеренные и восстановленные значения, который приведен в Табл.9, а его график – на Рис.3. В этой же таблице для восстановленных данных приведена стандартная погрешность восстановления.

Таблица 9.

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

Год Значение (мм) Погрешность восст. (мм.)
58.6 19.5
50.5 19.5
51.4 19.5
51.4 19.5
53.2 19.5
54.1 19.5
91.6 19.5
96.1 19.5
48.8 19.5
77.3 19.5
69.3 19.5
91.6 19.5
60.4 19.5
63.0 0.0
52.0 0.0
87.0 0.0
20.0 0.0
48.0 0.0
77.0 0.0
9.0 0.0
88.0 0.0
11.0 0.0
27.0 0.0
13.0 0.0
87.0 0.0
32.0 0.0
29.0 0.0
63.0 0.0
107.0 0.0
39.0 0.0
93.0 0.0
88.0 0.0
46.0 0.0
49.0 0.0
50.0 0.0
12.0 0.0
8.0 0.0
9.0 0.0
10.0 0.0
18.0 0.0
62.0 0.0
48.0 0.0
73.0 0.0
4.0 0.0
58.0 0.0
44.0 0.0
36.0 0.0
18.0 0.0
54.0 0.0
41.0 0.0
25.0 0.0
13.0 0.0
50.0 0.0
50.0 0.0
17.0 0.0
4.0 0.0
51.0 0.0
57.7 19.5
85.4 19.5
72.0 19.5

Последовательность расчета - student2.ru

Рис.3. Восстановленный ряд осадков за июль на метеостанции Мыс Каменный.

Результаты оценки эффективности наблюденных и восстановленных данных приведены в Табл.10, из которой следуют выводы:

- восстановлены осадки июля за 16 лет (с 1929 по 1949 гг. и с 1999 по 2002 гг.);

- эффект восстановления составляет 36%;

- стандартная погрешность восстановления составляет 19.5 мм;

- отношение стандартных отклонений наблюденных и восстановленных рядов составляет 1.7, что дает статистику критерия Фишера 2.8;

- отношение средних значений наблюденных и восстановленных рядов составляет 0.64, что дает статистику критерия Стьюдента 3.19;

- гипотеза об однородности дисперсий и средних наблюденных и восстановленных рядов отклоняется, но их отклонение обусловлено небольшим объемом восстановленных данных (16 лет), который не включает полные циклы лет повышенных и пониженных осадков.

Таблица 10.

Результаты оценки эффективности восстановления ряда осадков за июль на метеостанции Мыс Каменный.

Код Наименование Кол-во факт. Кол-во восст. Эффект.восст.% Средн.погр. восст. Стд.откл.исх. Стд.откл.восст. Критерий Фишера Средн.исх. ряда Средн.восст. ряда Критерий Стьюдента
Мыс Каменный 19.5 27.8 16.7 2.79 42.8 66.8 -3.19

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