Файл: Рождественский, А. В. Статистические методы в гидрологии.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 |
(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 |