Файл: Рождественский, А. В. Статистические методы в гидрологии.pdf

выборки п = 10 величина отрицательного смещения составляет 25%; при тех же исходных статистических параметрах, но при п = 25 ве­

личина смещения получена равной 12%. Указанные величины сме­ щения соответствуют биномиальному распределению, принимае­ мому за исходное. Связь между математическим ожиданием (средним значением) выборочных значений коэффициентов вариации

(С„выб) и несмещенной оценкой коэффициента вариации (С„) для рядов, не обладающих внутрирядной корреляцией, выражается фор­ мулой Е. Г. Блохинова [18]

Наличие внутрирядных связей увеличивает отрицательное сме­ щение выборочных коэффициентов вариации. Это обстоятельство учтено в формуле Крицкого—Менкеля [78]

Здесь С' виб— среднее значение выборочных коэффициентов ва­ риации для рядов объемом п членов, обладающих внутрирядным

коэффициентом корреляции (г); Cv Выб — среднее значение выбо­

рочных коэффициентов вариации для рядов, не имеющих внутри­ рядной связи.

Подставляя в зависимость (5.52) равенство (5.51), получаем

Зависимость (5.53) соответствует равенству CS = 2CV в ряду ис­

ходных величин. Из уравнения (5.53) можно получить выражение для несмещенной оценки коэффициента вариации (Cv) в зависи­

мости от среднего значения выборочных оценок коэффициентов вариации (С ввыб), числа членов совокупности п и коэффициента

внутрирядной корреляции г. Однако такое решение громоздко и не­

удобно. Поэтому более целесообразно применять графическую за­ висимость Cv = f (С' б; п; г). Пример такой зависимости при г =

= 0,3 представлен на рис. 5.7.

Использование метода статистических испытаний позволило оце­ нить возможность применения формулы (5.53) при С3ф2Су. Ока­

залось, что при С „^0,5 оценки несмещенных значений коэффициен­ тов вариации, полученные по зависимостям, вытекающим из урав­ нения (5.53), хорошо соответствуют оценкам, получаемым методом статистических испытаний при. всех обычно используемых в гидро­ логических расчетах отношениях Cs/Cv. При коэффициентах вариа­ ции, приближающихся к 1,0, и отношении Cs/Cv> 2 величины не-

2 9 4

смещенных значений С„, получаемые по зависимостям (рис. 5.7), вытекающим из уравнения (5.53), оказываются ниже, чем по ме­ тоду статистических испытаний. Из этого следует, что величина сме­ щения выборочных значений коэффициентов вариации при указан­ ных значениях параметров исходного распределения уравнением (5.53) учитывается не полностью. Более точный учет смещенности

Рис. 5.7. График зависимости Cv= f(C„Bbl6) при различных объемах

выборки (тг=10; п = 25; /г=50), построенной по формуле (6.53) при

Ca=2Cv и /•(,,+1=0,3.

выборочных оценок коэффициентов вариации возможен на основа­ нии использования зависимостей Cv = f ( С 'выб, г, п, Cs/Cv), полу­

ченных методом статистических испытаний.

Примеры таких зависимостей, соответствующих биномиальному распределению и распределению Крицкого—Менкеля, при Cs/Cv=

= 3 и 4 и при г= 0 и г= 0,3 представлены на рис. 5.8. Эти примеры показывают, что выборочные оценки коэффициентов вариации, полученные с применением кривой Крицкого—Менкеля, имеют

295

меньшее отрицательное смещение, чем оценки, соответствующие би­ номиальному распределению. Следует, однако, отметить, что это различие является достаточно существенным лишь при малых объ­ емах выборок и при больших значениях отношения С«/С„; эти раз­ личия увеличиваются также с возрастанием внутрирядной связанно­ сти. При истинных значениях С„=1,0, CS = 4CV, г = 0,3 и п=10 ве­

личина смещения достигает при биномиальном распределении 27%, а при использовании распределения Крицкого—Менкеля 19%. С уменьшением C„, Cs/Cv и увеличением r u n отрицательное сме­

щение выборочных коэффициентов вариации уменьшается. Рассмотрим случайное рассеивание выборочных коэффициентов

вариации (0^ ,). Общее представление о случайных колебаниях выборочных значений коэффициентов вариации, определенных для

генеральной совокупности с параметрами х=1,0, С„ = 0,3, С«= 0,6, дают кривые обеспеченности, приведенные на рис. 5.9. Для построе­ ния этих кривых использовано 600 выборок, включающих 25 членов каждая. Выборочные значения коэффициентов вариации опреде­ лены методом моментов при различных коэффициентах корреляции между смежными членами (г = 0,0; 0,5; 0,7; 0,9). На рис. 5.9 видно, что рассеяние выборочных коэффициентов вариации увеличивается с возрастанием коэффициента корреляции между смежными чле­ нами ряда; при этом также возрастает положительная асимметрия закона распределения. На рис. 5.9, наряду с эмпирическими пред­ ставлены и аналитические кривые обеспеченности, соответствующие биномиальному закону распределения. Аналитические кривые по­

строены по параметрам (С®, Cv , Cs ) смоделированных рядов.

Хорошее совпадение эмпирических и аналитических кривых обеспе­ ченностей дает основание при обработке рядов, полученных мето­ дом статистических испытаний, использовать не всю совокупность величин, образующих статистический ряд коэффициентов вариации, а лишь статистические параметры этих рядов. Такой вывод был сде­ лан при рассмотрении выборочных оценок средних.

Для рядов, не обладающих внутрирядной связью, среднее ква­ дратическое отклонение выборочных коэффициентов вариации (Ocv), полученных методом статистических моментов, может быть

определено по формуле [18]

(5.54)

которая рекомендуется в «Указаниях по определению расчетных гидрологических характеристик» СН 435-72.

Впоследствии Е. Г. Блохинов [18] ввел поправку, установленную эмпирическим путем,

(5.55)

298

Рис. 5.9. Кривые обеспеченности выборочных коэффициентов вариации по данным моделируемых

методом Монте-Карло рядов (дс=1; Си=0,3; Cs= 0,6) с различными коэффициентами корреляции между смежными членами ряда.

/ — г - 0; 2 — г - 0 . 5 ; 3 — г - 0 . 7 ; 4 — г = 0,9.

где п —-объем выборки; С„ — коэффициент вариации величин ис­

ходного ряда.

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

в «Указаниях по определению расчетных гидрологических харак­ теристик» СН 435-72. Формулы (5.54) — (5.56) соответствуют усло­ вию, когда внутрирядные связи отсутствуют, и получены примени­ тельно к соотношению CS = 2C„.

А. Ш. Резниковский [37] для выборочных Cv, полученных ме­

тодом моментов, опираясь на метод статистических испытаний, оце­ пил влияние внутрирядной корреляции для случая CS=2CV на ве­ личину ocv. Введя в зависимость (5.55) эмпирическую поправку,

Резниковский представил ее в виде

На основании указанного выше обширного материала, получен­ ного в ГГИ методом статистических испытаний, была произведена оценка формулы (5.57) в отношении возможности ее применения при соотношениях Cs/Cv ^=2 и при использовании в качестве исход­

ных законов распределения кривых Крицкого—Менкеля и бино­ миального закона. Указанное сопоставление выполнялось в отно­ шении коэффициентов вариации выборочных коэффициентов вариа­ ции, в связи с чем формула (5.57) применялась в виде

Результаты сопоставления показывают, что при отсутствии вну­ трирядной корреляции в исходном ряду (г = 0) и при CS = 2C'„ фор­ мула (5.58) показала хорошее согласование с данными статистиче­ ского моделирования. При исходных соотношениях, уклоняющихся от CS = 2CV, по формуле (5.58) получаются меньшие значения оце­

нок, чем оценки, полученные методом статистических испытаний. Это различие особенно существенно при малых п и больших значе­ ниях коэффициента вариации исходного ряда (Cv), отношения Cs/Cv и коэффициента внутрирядной корреляции (г). Для этих условий

эмпирическая поправка Резниковского, введенная в зависимость (5.55), должна быть уточнена. Впредь до завершения такого иссле­ дования оценку случайного рассеивания выборочных коэффициен­ тов вариации для соотношения CS¥=2CV и при наличии внутриряд­

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

300

Коэффициенты вариации выборочных коэффициентов вариации (Cv ), полученные по выборкам, смоделированным с использова-

нием биномиального закона и распределения Крицкого—Менкеля при С8ф2Съ, различаются несущественно. Ряды выборочных коэф­

фициентов вариации имеют положительную асимметрию. Эта асим­ метрия увеличивается с возрастанием отношения Cs/Cv, Cv и г и

уменьшается с увеличением объема выборки. Отношение коэффи­ циента асимметрии к коэффициенту вариации выборочных коэффи­ циентов вариации (Csc /С„с ) увеличивается с возрастанием ис-

ходного коэффициента вариации.

оценка выборочных коэффициентов асимметрии

В качестве иллюстрации случайных рассеиваний выборочных коэффициентов асимметрии на рис. 5.10 представлены кривые обес­ печенности этого параметра. При построении указанных кривых обеспеченности использовано 600 значений Cs, полученных для вы­ борок, включающих каждая 25 членов; общий объем смоделиро­ ванной последовательности jV = 15 000; исходные параметры, приня­

тые при моделировании, следующие: х=1, С„ = 0,3, Cs= 0,6; коэф­ фициенты корреляции между смежными членами ряда принимались равными г= 0,0; 0,5; 0,7; 0,9.

Как и при исследовании смоделированных статистических сово­ купностей других параметров (среднего значения, коэффициента ва­

риации), были рассчитаны стандартные параметры (Cs, С» , Cs )

рассматриваемых совокупностей. На основании этих параметров построены биномиальные кривые обеспеченности (рис. 5.10), кото­ рые вполне удовлетворительно соответствуют распределению выбо­ рочных Cs.

На рис. 5.10 видно, что распределение выборочных коэффициен­ тов асимметрии обладает большим случайным рассеиванием, чем

ряды выборочных значений х и Cv. Это распределение характери­

зуется существенным отрицательным смещением, которое увеличи­ вается с увеличением коэффициента корреляции между смежными членами ряда. Отрицательное смещение выборочных коэффициентов асимметрии при отсутствии внутрирядных связей может быть уст­ ранено путем расчета несмещенной оценки (С8цесм) по одной из сле­ дующих формул:

(формула, применяемая американскими гидрологами),

301