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

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

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

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

Добавлен: 15.10.2024

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

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

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

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

Однако этот вариант нормализации случайной переменной х

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

При выполнении нормализации по второму способу (т. е. через кривую обеспеченности, принятую в качестве аппроксимирующей эмпирический ряд) полностью сохраняются свойственные исход­ ному эмпирическому ряду флуктуации. Например, довольно часто величины эмпирического ряда в зоне малых обеспеченностей не со­ гласуются с общим направлением аналитической кривой обеспе­ ченности (так называемые иногда «отскакивающие точки»). Такое рассеяние, свойственное исходному ряду, сохраняется и при нор­ мализации по рассматриваемому второму способу.

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

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

с параметрами cr= 1 и х = 0; для осуществления такого преобразо­ вания в нижней части графика (рис. 2.17) необходимо использовать нормальную кривую распределения с этими параметрами.

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

В том случае, когда закон распределения рассматриваемой ста­ тистической совокупности (кривая а на рис. 2.17) и закон распре­

деления, к которому преобразуется случайная переменная (кри­ вая б на рис. 2.17), табулированы, преобразование может быть осу­

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

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

125


При нормализации непрерывной случайной величины х

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

вформе, удобной для решения гидрологических задач, приведена

вработе [9].

Описанный прием нормализации исходного ряда непосредст­ венно через эмпирические обеспеченности использовал Г. А. Алек­ сеев [9] при исследовании корреляционных связей характеристик гидрометеорологического режима. В связи с этим он пересчитал таблицу нормального закона распределения в таблицу, позволяю­ щую определить по заданным величинам эмпирической обеспечен­ ности ординаты кривой нормального закона [9].

II. Преобразование нормального закона распределения в асим­ метричный закон, простирающийся от нуля до бесконечности, осно­ вывается на использовании общих соотношений, определяющих аналитическое преобразование одного закона распределения в дру­ гой. Задача ставится следующим образом. Задан закон распреде­ ления случайной величины у в виде Р (у); случайная величина х связана с у функциональной зависимостью x = f(y). Нужно опре­ делить закон распределения Р (х) случайной величины х.

Будем считать, что обратная функция у=<р{х) однозначна (рис. 2.18). Вероятность того, что случайная переменная у\ попадет

в интервал (у, y + dy), должна

быть в точности равна

(вследствие

однозначной связи между х н у )

вероятности того, что значение xt

заключено в соответствующем промежутке х, x + dx, т.

е.

P ( x < x l < x - { - d x ) = P ( y < y i < y + d y ) ,

 

или

 

 

P (x)dx= P * (у) dy.

(2.107)

Аналогичным образом (в силу принятой однозначной связи между х и у) в интегральном выражении вероятность непревышения (или превышения) данного значения исходной переменной х и преобразованной у будет одинакова. В этом заключается важная

особенность рассматриваемого преобразования как равновероятно­

стного.

Из равенства (2.107) следует, что

Р(х) = Р*[у(х) 1

(2.108)

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

Введя вместо у функцию у=ц> (х), получаем

(2.109)

126


Таким образом, выражение (2.109) определяет закон распреде­ ления случайной величины х по известному закону распределения непрерывной случайной переменной у в зависимости от вида функ­

ции преобразования х = / (у).

Если функция у=кр(х) является двузначной (рис. 2.18), т. е. каждому значению х соответствует два значения у, формулу (2.109) нужно применить дважды — к каждой ветви кривой у =<pi (лг) и У = фг (х)

Я(х) = Я* (у) [^(х)] dx + Я*(У)ЫА>]

Рис. 2.18. Иллюстрация однозначной (я) и двузначной (б) функций.

Числовые характеристики (моменты) закона распределения случайной величины х можно определить, не устанавливая выра­ жение Р (у)\ это можно сделать проще. В самом деле, математиче­

ское ожидание величины х равно

 

с о

 

 

х =

J х P(x)dx,

 

 

— с о

 

 

но для однозначной функции, как показано выше,

 

Я (х) d x = Я* (у) dy,

 

а * = / (У); следовательно,

 

 

 

 

с о

 

 

■*=jf ( y ) P * ( y ) d y .

 

 

с о

 

 

Аналогично для начальных (т)

и центральных (ц)

моментов

в общем виде получим:

 

 

 

СО

 

СО

(2.110)

mk= jx kP ( x ) d x =

j-■[ / {y)]kP* (y)dy,

—co

—oo

 

oo

 

co

(2.111)

H== J (x —x)* P ( x ) d x =

j[f(y) —x]kP*{y)dy.

127


В частности, для дисперсии имеем

со

 

4 = j [ f ( y ) - x ] 2P*{y) dy.

(2.112)

Опираясь на полученные уравнения, покажем, что осуществлен­

ное в § 8 преобразование переменной по соотношению

2 = ---------

не приводит к нарушению нормального закона. Учтем,

Ож

что z явля­

ется однозначной функцией х, а пределы изменения как х, так и 2 от —оо до оо.

Выразим х через функцию преобразования, т. е. получим обрат­

ную функцию

X = Z a x -\-Х ,

dx

 

 

 

 

 

откуда —~ = a A-.

 

 

 

 

 

Подставляя в (2.108)

эти соотношения и уравнение (2.98), полу­

чаем:

 

 

 

 

 

 

( г а + ж - * )

_ _ f l

 

 

7

s e~

2"

<2-|13)

Уравнение (2.113) показывает, что

выполненное

преобразова­

ние величины х, распределенной

по нормальному закону,

не изме­

няет нормальности распределения и для случая х = 0 и а^=1.

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

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

Исходное уравнение нормального закона для плотности вероят­ ности величины у, как известно, может быть записано в форме

 

(у - у )2

Р( У)= Су К2тс

(2.114)

Чтобы приспособить этот закон распределения для описания асимметричных статистических совокупностей переменных, изме­ няющихся в пределах от 0 до оо, преобразуем г/= 1пх; это равно­ сильно тому, что принята обратная функция x = f(y) в виде х = еу.

В таком случае

rfy__ 1 dx х

128


Учитывая это, распределение (2.114) преобразуем в логарифми- чески-нормальное распределение

 

( i n

— I n Л1)2

 

Р(х)

___ 1__ е

2

(2.115)

I n X

Таким образом,

X V2n

 

распределение

логарифмически-нормальное

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

Поскольку отрицательные величины не имеют логарифмов, то лога­ рифмически-нормальное распределение описывает статистическую совокупность величин, изменяющихся в пределе от 0 до оо.

Логарифмически-нормальное распределение соответствует та­ кой статистической модели, в которой предполагается, что форми­ рование случайной величины происходит в результате умножения

большого числа влияющих на эту величину факторов,

каждый

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

возникает

Уместно напомнить, что нормальное распределение

в результате сложения равнозначно

влияющих на рассматривае­

мую случайную величину факторов

(или ошибок измерения).

Рассмотрим зависимости, устанавливающиеся между

парамет­

рами логарифмически-нормальной кривой распределения исходной величины х и нормальной кривой величины у = \пх. Для этого зави­

симость (2.115) запишем в виде

 

1

 

(у-у)2

 

Р(х)

е

2°у

(2.116)

 

еу ау У 2к

 

 

 

Применяя рассмотренный ранее метод моментов

к логарифми-

чески-нормальному закону распределения случайной переменной х,

записанной в виде (2.116), получаем для начальных моментов

о о

о о

 

aA= j

x kP (х) dx = J ek-yP(x)dx,

(2.117)

о

о

 

или, имея в виду, что dx = ev dy, получаем

 

 

 

(2.118)

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

степени,

выделяя

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

 

 

2йуа2 — (у — у)2

(у — у)2 — 2£з2 (у — у) +

26уа2 +

k2^

ц

(2.119)

9 Зак . № 88

129