Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология
Образование Политология Производство Психология Стандартизация Технологии


Определение коэффициентов вариации и ассиметрии при наличии длительного ряда наблюдений



 

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

При наличии длительных наблюдений Указаниями СН 371-67 предусматриваются два метода определения в зависимости от изменчивости годового стока.

Если изменчивость годового стока невелика и характеризуется коэффициентом вариации , рекомендуется следующая формула

 

(3.1)

 

где kiмодульный коэффициент стока каждого года; n – число лет наблюдений (число членов статистического ряда).

При n< 30 формула (3.1) используется в виде

 

(3.1')

 

Формула (3.1) представляет собой выражение второго момента площади кривой распределения относительно центральной ординаты, а метод определения и по формулам (3.1 ) и (3.6) называется методом моментов.

Относительная средняя квадратическая ошибка коэффициента вариации, вычисленная по этим формулам, равна

 

(3.2)

 

Формула (4.2) предложена Е. Г. Блохиновым и рекомендована Указаниями СН 371-67.

Для нормального распределения параметра или близкого к нему более обоснованной формулой определения является:

(3.3)

 

Значения в зависимости от п и вычисленные по формулам (3.2) и (3.3), приводятся в таблицах 3.1 и 3.2.

 

Таблица 3.1

Относительные средние квадратические ошибки (%) коэффициента вариации, вычисленные по формуле (3.2)

 

  Число лет наблюдений n
0, 20 9, 3 8, 6 8, 1 7, 6 7, 2
0, 30 9, 6 8, 8 8, 3 7, 8 7, 4
0, 40 9, 8 9, 1 8, 5 8, 0 7, 6

 

Примечание. При промежуточных значениях Cv и п определяется по интерполяции.

 

Таблица 3.2

Относительные средние квадратические ошибки (%) коэффициента вариации, вычисленные по формуле (3.3)

 

n
0, 20 0, 40 0, 60 0, 80 1, 00
23, 2 25 7 29, 4 33, 7 38, 6
16, 5 18, 2 20, 7 23, 9 27, 4
13, 4 14, 8 16, 9 19, 5 22, 3
10, 4 11, 5 13, 1 15, 1 17, 3
7, 4 8, 2 9, 3 10, 7 12, 3

 

Из этих таблиц видно, что для определения < 0, 50 с точ-ностью ±10% необходимо иметь ряд наблюдений продолжительностью 50—60 лет, а вычисления по формулам (3.2) и (3.3) дают малую разницу, которая практически не сказывается на конечных результатах расчетов.

При большой изменчивости годового стока ( > 0, 50), что характерно для рек засушливых районов, коэффициент вариации определяется методом наибольшего правдоподобия. По этому методу значение Cv устанавливается в зависимости от параметра l, который представляет собой среднее значение логарифмов модульных коэффициентов lg k и вычисляется по равенству

 

(3.4)

Имея значение l, по таблице 3.3 легко установить величину Cv. Относительная средняя квадратическая ошибка Cv, установленного с помощью параметра l, приближенно равна

 

(3.5)

Более точные значения , вычисленные методом наибольшего правдоподобия, приводятся в таблице 3.4.

Значения коэффициентов вариации, определенные методом моментов и методом наибольшего правдоподобия при Сv > 0, 5, различаются не более чем на 2–3% (таблица 3.4). В то же время метод наибольшего правдоподобия отличается сложностью и громоздкостью вычислений. Поэтому можно считать, что метод наибольшего правдоподобия для определения Cv не имеет преимущества перед методом моментов и при массовых расчетах Cv лучше пользоваться последним, т.е. формулой (3.1), который оправдал себя при многочисленном практическом применении.

Коэффициенты вариации, установленные по формуле (4.1) и по табл. 3.4 с помощью параметра l, принимаются в качестве расчетных, когда их средняя квадратическая ошибка, вычисленная по формуле (3.2), не превышает следующие пределы:

 

Интервалы значений Cv 0, 2—0, 3 0, 4—0, 8 больше 0, 8

Допустимая средняя

квадратическая ошибка, % 20 10 5

 

Таблица 3.3

Относительная средняя квадратическая ошибка (%) коэффициента изменчивости

 

Cv Число лет наблюдений п
0 50 21, 5 15, 2 12, 4 10, 8 9, 6 8, 8 8, 1 7, 6 7, 2 6, 8
0, 75 20, 5 14, 5 11, 8 10, 3 9, 2 8, 4 7, 7 7, 2 6, 9 6, 5
1, 00 19, 4 13, 7 11, 2 9, 7 8, 7 7, 9 7, 3 6, 8 6, 5. 6, 1
1, 25 18, 2 12, 8 10, 5 9, 1 8, 1 7, 5 6, 9 6, 4 6, 1 5, 8
1, 50 16, 9 12, 0 9, 8 8, 5 7, 6 6, 9 6, 4 6, 0 5, 7 5, 4
1, 75 15, 8 11, 2 9, 1 7, 9 7, 0 6, 4 6, 0 5, 6 5, 3 5, 0
2, 00 14, 6 10, 4 8, 4 7, 3 6, 6 6, 0 5, 5 5, 2 4, 9 4, 6

Примечание. При промежуточных значениях Cv и п средняя квадра. тическая ошибка определяется по интерполяции.

 

Таблица 3.4

Средние значения Cv, определенные методами моментов и наибольшего правдоподобия

 

 

 

Истинное значение Cv Объем выборки и Значение Cv
по методу моментов по методу наибольшего правдоподобия
0, 25   0, 50   1, 00 0, 25 0, 24 0, 49 0, 50 0, 97 0, 98 0, 25 0, 25 0, 49 0, 50 0, 99 1, 00

 

Максимальные ошибки Cv значительно превышают среднее квадратическос значение. Они возможны тогда, когда в короткий ряд наблюдений, по которому определяются Cv, входит очень многоводный или маловодный год, повторяемость которого значительно реже чем один раз в п лет. Однако вероятность больших ошибок Cv, в два-три раза превышающих ее среднюю квадратическую величину, очень мала.

 

Коэффициент асимметрии Cs характеризует несимметричность ряда величин стока относительно его среднего значения. Это менее устойчивый параметр кривой распределения или обеспеченности и для надежного его определения требуется ряд наблюдений над стоком более 100—150 лет. По имеющимся рядам наблюдений можно установить лишь приближенное значение коэффициента асимметрии. Для этой цели используется формула третьего момента

(3.6)

 

Относительная средняя квадратическая ошибка Cs зависит от коэффициента вариации Сv и числа лет наблюдений п. Ее значение при Cs = 2Cv с учетом асимметричного распределения годовых величин стока может быть определена по формуле С. Н. Крицкого и М. Ф. Менкеля. Так при обычных значениях CV=0, 20—0, 40 средняя квадратичная ошибка при n=20—30 составляет около 14—16% от значения коэффициента, а при n = 50 уменьшается до 10—12%. Относительная среднеквадратичная ошибка определения коэффициента Cs вычисляется по формуле

 

(3.7)

 

Величины этих ошибок коэффициентов Cv и Cs, определенных по формулам (3.2-3.3) и (3.7), даны в таблицах 3.2–3.3 и 3.6.

Для наиболее распространенных величин годового стока рек лесной и лесостепной зон Cv=0, 20÷ 0, 50. Средняя ошибка его опре­деления при n=20-30 лет находится в пределах 15, 5-19, 5%.

Для надежного подсчета коэффициента Сs необходимо иметь, длительный ряд (порядка 100 и более), что практически бывает крайне редко. Коэффициент обладaeт более значительной устойчивостью во времени и для надежного его подсчета достаточен более короткий ряд наблюдении. Поэтому обычно величину Сs принимают кратно значению Сv т.е.

Сs=а Сv (3.8)

 

 

Таблица 3.5

Относительные средние квадратические ошибки определения коэффициента асимметрии в % при Сs =2Cv

 

 

п Коэффициент вариации Cv
0, 10 0, 20 0, 40 0, 80 1, 00

 

Величина а для различных гидрологических характеристик принимается различной. Так, для биноминальной асимметричной кривой пределом этой величины могут быть от а=2 до ; Кминотношение наинизшего расхода к среднему данного ряда). Для среднегодового стока равенство коэффициентов CS=Cv справедливо только при малых величинах Cv. При больших величинах Cv следует применять , что дает, как правило, величину CS> 2CV. Для других гидрологических величин соотношение между Cv и CS, будет указано ниже.

При наличии ряда наблюдений порядка 20 и более построение кривой обеспеченности выполняется с помощью таблицы вспомогательных величин, образец которой приведен в таблице. 3.6.

В графе 12 подсчитывается обеспеченность точек (модульных коэффициентов) эмпирической кривой р % определяется по формуле

 

(3.9)

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

 

Таблица 3.6

Вспомогательные величины для построения кривой обеспеченности

 

№ пп Год Q м3/сек в убывающ ем порядке (k -1) (k -1)2 (k -1)3 р, %  
  +   ­     +     ­  
годы Q, м3/сек    
 
                   
                   
                     
                     
Сумма        

 

На график (рис. 3.1) наносят точки ряда наблюдений (полые кружки) и получают кривую обеспеченности. Пользуясь табл. 3.6.можно построить на этом же графике теоретическую кривую обеспеченности. Для этого из графы 5 надо получить величину откуда можно получить

 

 

Для контроля правильности подсчета сумма модульных коэффициентов должна удовлетворять равенству , а в графах 7, 8 . Используя сумму , из графы 9 можно получить значение значение Cv по формуле (3.1-3.1'). В первом приближении применяем CS=2CV.

Для получения теоретических точек кривой обеспеченности следует пользоваться таблицей Фостера – Рыбкина, в которой приведены отклонения ординат кривой обеспеченности от середины Ф при Хср = 1 и Cv=1 (см. табл. 3.9).

Пользуясь табл. 3.7, можно подсчитать значения модульных коэффициентов по формуле

(3.10)

 

Затем по формуле (4.5) определить значения расходов разной обеспеченности:

(3.11)

р, % 0, 1 99, 9
Сs=0, 57
Фр                          
Фр Cv                          
kp = Фр Cv+1                          
Qp = Q0 kp                          
 
kp                          
Qp = Q0 kp                          

Результаты заносятся в табл. 3.7.

 

Таблица 3.7

 

Полученные значения Qр – наносят на график, приведенный на рис. 3.1 (на рисунке значения Qр. обозначены сплошными кружками). Если проведенная по этим кружкам кривая обеспеченности (пунктирная) не совпадает с кривой, построенной по наблюденным данным, то это обозначает, что принятое значение Cs для построения теоретической кривой не соответствует значению Cs соответствующего фактическим данным. Это значение необходимо изменять до тех пор, пока теоретическая кривая не совпадет с кривой натурных наблюдений.

Биноминальная кривая обеспеченности хорошо согласуется с натурными данными. Однако она имеет ряд недостатков, одним из которых является то, что при соотношении Cs < 2 Cv встречающемся на практике при расчетах годового стока рек засушливых районов, ординаты кривой обеспеченности в нижней части имеют отрицательные значения. Для исключения этого недостатка С. Н. Крицкий и М. Ф. Менкель предложили новые кривые обеспеченности [35]. По этим кривым при любом соотношении Cs и Cv нулевая ордината получается только при обеспеченности в 100%. При CS=2CV кривая обеспеченности совпадает с биномальной ассиметричной кривой. Ассиметричные кривые на клетчатке вероятности спрямляются не полностью и имеют выпуклость тем большую, чем больше коэффициент асимметрии. При положительной асимметрии выпуклость обращена вниз, при отрицательной – вверх.

Характеристика связи между обеспеченностью и повторяемостью расходов приведена в табл. 3.8.

Применение методов теории вероятности и математической статистики в решении гидрологических задач получило широкое распространение как у нас, так и за границей. Однако следует отметить несовершенство и условность этих методов. Кривые обеспеченности, построенные теоретическим путем, необходимо каждый раз анализировать и проверять, используя фактические наблюдения.

 

 

Рис. 3.1. Кривая обеспеченности средних годовых расходов на клетчатке вероятности

 

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

Таблица 3.8

Характеристики лет разной обеспеченности

 

Обеспечен­ность, % Повторяется 1 раз в n лет Характеристика года
0, 1 Катастрофически многоводный
Очень многоводный
Многоводный
Средней многоводности
Средней многоводности
Медианный
Средней маловодности
Средней маловодности
Маловодный
Маловодный
Очень маловодный
99, 9 Катастрофически маловодный

Таблица 3.9

Отклонение ординат биноминальной асимметрической кривой обеспеченности от середины (от 1, 0) при Cv=1

CS Обеспеченность р % CS
0, 01 0, 1 99, 9
0, 00 3, 72 3, 09 2, 33 1, 88 1, 64 1, 28 0, 84 0, 67 0, 52 0, 25 0, 00 -0, 25 -0, 52 -0, 67 -0, 84 -1, 28 -1, 64 -1, 88 -2, 33 -3, 09 0, 00
0, 05 3, 83 3, 16 2, 36 1, 90 1, 65 1, 28 0, 84 0, 66 0, 52 0, 24 -0, 01 -0, 26 -0, 52 -0, 68 -0, 84 -1, 28 -1, 62 -1, 86 -2, 29 -3, 02 0, 05
0, 10 3, 94 3, 23 2, 40 1, 92 1, 67 1, 29 0, 84 0, 66 0, 51 0, 24 -0, 02 -0, 27 -0, 53 -0, 68 -0, 85 -1, 27 -1, 61 -1, 81 -2, 25 -2, 95 0, 10
0, 15 4, 05 3, 31 2, 44 1, 68 1, 30 0, 84 0, 66 0, 50 0, 23 -0, 02 -0, 28 -0, 54 -0, 68 -0, 85 -1, 26 -1, 60 -1, 82 -2, 22 -2, 88 0, 15
0, 20 4, 16 3, 38 2, 47 1, 96 1, 70 1, 30 0, 83 0, 65 0, 22 -0, 03 -0, 28 -0, 55 -0, 69 -0, 85 -1, 26 -1, 58 -1, 79 -2, 18 -2, 81 0, 20
0, 25 4, 27 3, 45 2, 50 1, 98 1, 71 1, 30 0, 82 0, 64 0, 49 0, 21 -0, 04 -0, 29 -0, 56 -0, 70 -0, 85 -1, 25 -1, 56 -1, 77 -2, 14 -2, 74 0, 25
0, 30 4, 38 3, 52 2, 54 2, 00 1, 72 1, 31 0, 82 0, 64 0, 48 0, 20 -0, 05 -0, 30 -0, 56 -0, 70 -0, 85 -1, 24 -1, 55 -1, 75 -2, 10 -2, 67 0, 30
0, 35 4, 50 3, 59 2, 58 2, 02 1, 73 1, 32 0, 82 0, 64 0, 48 0, 20 -0, 06 -0, 30 0, 56 -0, 70 -0, 85 -1, 24 -1, 53 -1, 72 -2, 06 -2, 60 0, 35
0, 40 4, 61 3, 66 2, 61 2, 04 1, 75 1, 32 0, 82 0, 63 0, 47 0, 19 -0, 07 -0, 31 -0, 57 -0, 71 -0, 85 -1, 23 -1, 52 -1, 70 -2, 03 -2, 54 0, 40
0, 45 4, 72 3, 74 2, 64 2, 06 1, 76 1, 32 0, 82 0, 62 0, 46 0, 18 -0, 08 -0, 32 -0, 58 -0, 71 -0, 85 -1, 22 -1, 51 -1, 68 -2, 00 -2, 47 0, 45
0, 50 4, 83 3, 81 2, 68 2, 08 1, 77 1, 32 0, 81 0, 62 0, 46 0, 17 -0, 08 -0, 33 -0, 58 -0, 71 -0, 85 -1, 22 -1, 49 -1, 66 -1, 99 -2, 40 0, 50
0, 55 4, 94 3, 88 2, 72 2, 10 1, 78 1, 32 0, 80 0, 62 0, 45 0, 16 -0, 09 -0, 34 -0, 58 -0, 72 -0, 85 -1, 21 -1, 47 -1, 64 -1, 92 -2, 32 0, 55
0, 60 5, 05 3, 96 2, 75 2, 12 1, 80 1, 33 0, 80 0, 61 0, 44 0, 16 -0, 10 -0, 34 -0, 59 0, 72 -0, 85 -1, 20 -1, 45 -1, 61 -1, 88 -2, 27 0, 60
0, 65 5, 16 4, 03 2, 78 2, 14 1, 81 1, 33 0, 80 0, 60 0, 44 0, 15 -0, 11 -0, 35 -0, 60 -0, 72 -0, 85 -1, 44 1, 19- 1, 59 -1, 84 -2, 20 0, 65
0, 70 5, 28 4, 10 2, 82 2, 15 1, 82 1, 33 0, 79 0, 59 0, 43 0, 14 -0, 12 -0, 36 -0, 60 -0, 72 -0, 85 -1, 18 -1, 42 -1, 57 -1, 81 -2, 14 0, 70
0, 75 5, 39 4, 17 2, 86 2, 16 1, 83 1, 34 0, 78 0, 58 0, 42 0, 13 -0, 12 -0, 36 -0, 60 -0, 72 -0, 86 -1, 18 -1, 40 -1, 54 -1, 78 -2, 08 0, 75
0, 80 5, 50 4, 24 2, 89 2, 18 1, 84 1, 34 0, 78 0, 58 0, 41 0, 12 -0, 13 -0, 37 -0, 60 -0, 73 -0, 86 -1, 17 -1, 38 -1, 52 -1, 74 -2, 02 0, 80
0, 85 5, 62 4, 31 2, 92 2, 20 1, 85 1, 34 0, 78 0, 58 0, 40 0, 12 -0, 14 -0, 38 -0, 60 -0, 73 -0, 86 -1, 16 -1, 35 -, 49 -1, 70 -1, 96 0, 85
0, 90 5, 73 4, 38 2, 96 2, 22 1, 86 1, 34 0, 77 0, 57 0, 40 0, 11 -0, 15 -0, 38 -0, 61 -0, 73 -0, 85 -1, 15 -1, 35 -1, 47 -1, 66 -1, 90 0, 90
0, 95 5, 84 4, 46 2, 99 2, 24 1, 87 1, 34 0, 76 0, 56 0, 39 0, 10 -0, 16 -0, 38 -0, 62 -0, 73 -0, 85 -1, 14 -1, 34 -1, 44 -1, 62 -1, 84 0, 95
1, 00 5, 96 4, 53 3, 02 2, 25 1, 88 1, 34 0, 76 0, 55 0, 38 0, 09 -0, 16 -0, 39 -0, 62 -0, 73 -0, 85 -1, 13 -1, 32 -1, 42 -1, 59 -1, 79 1, 00
1, 05 6, 07 4, 60 3, 06 2, 26 1, 88 1, 84 0, 75 0, 54 0, 37 0, 08 -0, 17 -0, 40 -0, 62 -0, 74 -0, 85 -1, 12 -1, 30 -1, 40 -1, 56 -1, 74 1, 05
1, 10 6, 18 4, 67 3, 09 2, 28 1, 89 1, 34 0, 74 0, 54 0, 36 0, 07 -0, 18 -0, 41 -0, 62 -0, 74 -0, 85 -1, 10 -1, 28 -1, 38 -1, 52 -1, 68 1, 10
1, 15 6, 30 4, 74 3, 12 2, 30 1, 90 1, 34 0, 74 0, 53 0, 36 0, 06 -0, 18 -0, 42 -0, 62 -0, 74 -0, 84 -1, 09 -1, 26 -1, 36 -1, 48 -1, 63 1, 15
1, 20 6, 41 4, 81 3, 15 2, 31 1, 91 1, 34 0, 73 0, 52 0, 35 0, 05 -0, 19 -0, 42 -0, 63 -0, 74 -0, 84 -1, 08 -1, 24 -1, 33 -1, 45 -1, 58 1, 20
1, 25 6, 52 4, 88 3, 18 2, 32 1, 92 1, 34 0, 72 0, 52 0, 34 0, 04 -0, 20 -0, 42 -0, 63 -0, 74 -0, 84 -1, 07 -1, 22 -1, 30 -1, 42 -1, 53 1, 25
1, 30 6, 64 4, 95 3, 21 2, 34 1, 92 1, 34 0, 72 0, 51 0, 33 0, 04 -0, 21 -0, 43 -0, 63 -0, 74 -0, 84 -1, 06 -1, 20 -1, 28 -1, 38 -1, 48 1, 30
1, 35 6, 76 5, 02 3, 24 2, 36 1, 93 1, 34 0, 72 0, 50 0, 32 0, 03 -0, 22 -0, 44 -0, 64 -0, 74 -0, 84 -1, 05 -1, 18 -1, 26 -1, 35 -1, 44 1, 35
1, 40 6, 87 5, 09 3, 27 2, 37 1, 94 1, 34 0, 71 0, 49 0, 31 0, 02 -0, 22 -0, 44 -0, 64 -0, 73 -0, 83 -1, 04 -1, 17 -1, 23 -1, 32 -1, 39 1, 40
1, 45 6, 98 5, 16 3, 30 2, 38 1, 94 1, 34 0, 70 0, 48 0, 30 0, 01 -0, 23 -0, 44 -0, 64 -0, 73 -0, 82 -1, 03 -1, 15 -1, 21 -1, 29 -1, 35 1, 45
1, 50 7, 09 5, 23 3, 33 2, 39 1, 95 1, 33 0, 69 0, 47 0, 30 0, 00 -0, 24 -0, 45 -0, 64 -0, 73 -0, 82 -1, 02 -1, 13 -1, 19 -1, 26 -1, 31 1, 50
1, 60 7, 31 5, 37 3, 39 2, 42 1, 96 1, 33 0, 68 1, 46 0, 28 -0, 02 -0, 25 -0, 46 -0, 64 -0, 73 -0, 81 -0, 99 -1, 10 -1, 14 -1, 20 -1, 24 1, 60

Окончание 3.9


Поделиться:



Популярное:

  1. A. Притяжения и отталкивания, силы отталкивания больше на малых расстояниях, чем силы притяжения. Б. Притяжения и отталкивания, силы отталкивания меньше на малых расстояниях, чем силы притяжения.
  2. Adjective and adverb. Имя прилагательное и наречие. Степени сравнения.
  3. D. Правоспособность иностранцев. - Ограничения в отношении землевладения. - Двоякий смысл своего и чужого в немецкой терминологии. - Приобретение прав гражданства русскими подданными в Финляндии
  4. D. ПРЕИМУЩЕСТВА ПРИСОЕДИНЕНИЯ К ГААГСКОМУ СОГЛАШЕНИЮ
  5. F70.99 Умственная отсталость легкой степени без указаний на нарушение поведения, обусловленная неуточненными причинами
  6. F71.98 Умственная отсталость умеренная без указаний на нарушение поведения, обусловленная другими уточненными причинами
  7. I Использование заемных средств в работе предприятия
  8. I. Методические принципы физического воспитания (сознательность, активность, наглядность, доступность, систематичность)
  9. I. О НОВОПРИБЫВШИХ ГРАЖДАНАХ.
  10. I. Предприятия крупного рогатого скота
  11. I. Придаточные, которые присоединяются непосредственно к главному предложению, могут быть однородными и неоднородными.
  12. I. СИЛЬНЫЕ СТОРОНЫ ПРЕДПРИЯТИЯ


Последнее изменение этой страницы: 2016-03-16; Просмотров: 2161; Нарушение авторского права страницы


lektsia.com 2007 - 2024 год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! (0.054 с.)
Главная | Случайная страница | Обратная связь