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

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 191

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

возрастает с увеличением исходного коэффициента вариации и при­ нятого отношения Cs/Cv', точно так же на величину смещения дейст­

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

ктически ощутимой в расчетах при малых объемах выборок

(п =

= 10-ь20), высоких коэффициентах вариации исходного ряда

(С„ =

= 1,0) и больших значениях внутрирядной корреляции (г = 0,3-^0,5). При использовании в расчетах значений выборочных коэффици­ ентов асимметрии (т. е. в случае трех свободно назначаемых пара­ метров) смещенность ординат кривой обеспеченности увеличивается вследствие смещенности коэффициентов асимметрии по сравнению

со случаем двух свободно назначаемых параметров (х, Cv).

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

оценка случайных колебаний выборочных ординат кривых обеспеченностей

В математической статистике [111] доказывается, что дисперсия cj2 выборочных оценок статистических параметров (у), в частности,

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

распределения х0 и центральных моментов р2, Цз, щ, в общей форме

может быть представлена в виде зависимости

av=

*ч‘3J + \ 2 2

з 2

2

 

 

, 2 2-j- 2Х|А2зjo№Г -a-]-2А,Х3а-а|1зГ

-[-

—1—2ai).,3—о

 

г — —I—2У,Уз

з

 

 

 

 

 

X

М-2ЛГЦз

 

 

 

 

(5.64)

!!.,

3

г

!j.2;j.3

4-2ХД.З

з

ц ,

г

ц2р.,

-4-2ХА,а

а

р ,

г

РзрЦ

I 1 1

X

X I , I

2 3

 

 

1

2 4

P.J

 

I 3 4

Из

 

 

.

 

<5®

 

 

 

.

 

(Эф

 

 

 

<5ф

 

 

 

 

 

Здесь /.1= —zr,

Х2 — — , л3=

-------, Хц—-----------частные произ-

 

 

 

дх

 

дцг

 

 

<5р3

 

 

 

 

дри

 

 

 

 

 

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

ческих моментов у=<р (х, ц2, Цз, ц4) ; о и г — среднее квадратическое

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

Опираясь на приведенное общее выражение для оценки диспер­ сии, можно получить записанные ранее формулы (5.48), (5.54) и (5.62). Общую схему использования уравнения (5.64) рассмотрим на примере вывода формулы стандартной ошибки ординат кривой

308


обеспеченности (вхт) применительно к биномиальному распреде­

лению с соотношением С«= 2СГ.

Врассматриваемом случае функция, связывающая величину хР

сопределяющими ее величинами, как известно, имеет вид

Х р = Х - \ - а х Ф = Х -{-У \> -2 ф ,

где Ф — табличное значение биномиального закона, зависящее от коэффициента асимметрии (Cs) и обеспеченности (р), которой от­ вечает величина хР>т. е. Ф= Ф(Св, р); ох — среднее квадратическое

отклонение ряда величин х.

Если CS = 2CV, величина Ф становится функцией выборочных оце­ нок второго центрального момента цг, поскольку Cs уже не высту­

пает в качестве самостоятельного параметра. Применительно к рас­

сматриваемому случаю Яз= Я4= 0 и уравнение

(5.64)

становится до­

статочно простым

 

 

 

 

 

 

о

дхр

dx„

 

Г- .

(5.65)

с +

дц2

—d—a-о

dp-2

х

111jrp-J

 

В свою очередь

 

 

 

 

 

 

д х г

 

 

дФ

 

(5.66)

Фг

2 Ур.2

Vv".2

^^2

 

 

 

Учитывая, что

 

2 _

|12

 

 

 

 

 

 

 

 

 

 

 

УГ ’

 

 

 

 

aL— 2|j.2(1 -f-3[i.2),

 

 

 

 

a- a

Г—

2(1-2

 

 

 

 

X

|Л2

п•

 

 

 

получаем

—5 = —)_|/|х2ф Л2-5^-(1+ 3(а2) +

р

п ~ ( 2 /(х2 '

I

п

 

 

+ 2

 

 

Л

2(4

(5.67)

 

2/(12

-]/(12Ф'

 

 

 

 

 

 

 

 

где Ф' = дФ_

 

 

 

 

 

 

 

d[i2

 

 

 

 

 

 

 

т-

 

,

дФ

1

дФ

 

Имея в виду, что ji2= crz и Ф = —------------- , можно получить

окончательно

 

 

 

дц2

о

оСа

 

 

 

 

 

 

 

 

- У

1 + 2 ( 4 - + а Ф , ) 2 ( 1 + 3 32) +

 

4 о ( А + а Ф , ) = - ^ - Т / Л ,

/-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.68)

3 0 9


где Л = 1+2^-^- + Ф8а ) (1 +3а2) + 4ог( “ТГ+фаа) ; Ф=

■нормированное отклонение от среднего ординат кривой

Со

обеспеченности; величины Ф приведены в таблице биномиального

Л

<?Ф

А

Cs, определяется как отно-

закона; Ф8= — -----производная Ф по

 

dCs

АФ

 

шение приращений

 

для тех обеспеченностей р, для которых

рассчитывается величина ох .

В том случае, если при построении кривой обеспеченности зна­ чение Cs определяется не по соотношению с величиной коэффици­

ента вариации, а непосредственно из выборочных оценок второго

и третьего моментов ( Cs~ %—), значение о*

выражается урав-

'

U-'2 '

Х-п

 

нением

 

 

 

% = 1^Г ^ Л + (6Ф^ + 4Фф^ + 34Ф"в,+

24Ф/

- 4Ф/ ) • (5.69)

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

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

(х, Сг>) и при фиксированном отношении Cs/Cv= 2. Уточнение, про­

изведенное Цинь Гу-анем, касалось учета корреляции между выбо­

рочными оценками х и а. Эта связь отсутствует для выборок из со­

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

/ 1+ 2 (4 + фл ) 2(1+ зс')+ *

+ 4 С „ ( 4 + ф,С „ )-4 Ф ,^ (1 + Ф С |;+ Ф 'Й ). (5.70)

Здесь выражение, стоящее под радикалом, представляет собой уточненное значение параметра А формулы (5.68).

Таким образом, зависимость (5.68) можно записать в виде

- ■ - к * * -

зю


Введя значение параметра А' (вместо Л) в уравнение (5.69), по­

лучаем уточненное выражение для а* » соответствующее случаю

использования трех свободно назначаемых параметров (х, Cv, Cs)

~ Y А +(бФ ?+4ФФ 5а + 34Ф?а2+24Ф 404 —4Ф^02 . (5.71)

РVn

Выражение (5.70) использовано для составления номограммы «гарантийной поправки», включенной в «Указания по определению расчетных ги­ дрологических характеристик» СН 435-72.

Эта номограмма представлена на рис. 5.13. При построении номограммы использо­ вано уравнение (5.70) в виде

У п

V А ' = f (Сv, Р).

В ранее существовавшем норматив­ ном документе для расчета максималь­ ных расходов воды (ГОСТ 3999/48) ис­ пользовалась номограмма, построенная на основании формулы (5.68). Е. Г. Блохинов предложил более простое по струк­ туре выражение для средней квадратиче­ ской ошибки ординат биномиальной кри­ вой распределения при CS = 2CV

1 +

Ci

{I OkDd p \2

 

 

 

 

Ы г ) • (5-72)

 

Здесь kP— табличное

значение

мо­

 

дульного коэффициента, соответствующее

 

рассматриваемой величине оХр.

 

5 10р%

Вывод формулы (5.72)

приведен в ра­

 

боте Крицкого и Менкеля

[79].

Рис. 5.13. Номограммы для

Имея в виду, что

ахР = CvХр kv,

определения

гарантийной

поправки.

формулы (5.72) можно получить выраже­ ние для коэффициента вариации выборочных ординат кривых обес­ печенности

р

__ ° гр

__

Cv

-1 [ .

,

1 + Су

dkp

(5.73}

dCv

 

Ьр

~

V n

V

^

2kl

 

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

3 1 1