- •МЕТОДЫ ОБРАБОТКИ ГИДРОЛОГИЧЕСКОЙ ИНФОРМАЦИИ
- •Оглавление
- •I. Основные понятия теории вероятностей, используемые в гидрологии
- •1. Обоснование применения вероятностных методов в инженерной гидрологии
- •2. Случайные события и случайные величины в гидрологии
- •3. Числовые характеристики случайных величин
- •4. Равномерное и нормальное распределения случайных величин
- •5. Распределение Пирсона III типа (непрерывное биномиальное распределение)
- •6. Распределение Крицкого-Менкеля (трехпараметрическое гамма-распределение)
- •7. Основные задачи математической статистики
- •8. Точечные и интервальные оценки параметров распределения
- •9. Статистическая гипотеза. Критерий статистической гипотезы
- •10. Случайные процессы в гидрологии
- •II. Статистические методы в инженерной гидрологии
- •11. Построение кривых обеспеченности по эмпирическим данным
- •12. Оценка параметров распределения по эмпирическим данным
- •13. Оценка погрешностей выборочных параметров распределения
- •14. Проверка однородности гидрологических рядов
- •15. Общие критерии согласия
- •16. Критерии случайности
- •17. Корреляционный анализ гидрологических и климатических рядов
- •18. Уравнение линейной регресии. Метод наименьших квадратов
- •19. Оценка точности уравнения линейной регрессии для двух переменных
- •20. Уравнения линейной регрессии при нескольких аргументах
- •21. Основные характеристики случайных процессов
- •22. Моделирование искусственных гидрологических рядов
- •23. Понятие о динамико-стохастических моделях многолетних колебаний речного стока
- •24. Вводные положения и термины основного нормативного документа
- •25. Определение расчетных гидрологических характеристик при наличии данных гидрометрических наблюдений
- •26. Учет неоднородности исходных данных и выдающихся значений
- •27. Расчетное внутригодовое распределение месячного стока
- •29. Наивысшие уровни воды
- •30. Определение расчетных гидрологических характеристик при недостаточности данных гидрометрических наблюдений
- •31. Учет данных кратковременных наблюдений
- •32. Расчетные расходы и уровни воды гидротехнических сооружений
- •33. Строительная климатология (СП 131.13330.2012)
- •Список литературы
- •Приложение А
Функцией распределения вероятностей F(x) СВ X в точке x называется вероятность того, что в результате опыта СВ примет значение, меньшее чем x, т.е. F(x) = p{X < x}. Если существует производная F(x), то она называется
плотностью распределения вероятностей СВ: f (x) = F (x) . Вероятность попа- |
|
′ |
|
дания случайной величины на заданный отрезок [α, β] равна: |
|
β |
|
p{α ≤ x ≤ β}= F(β) − F(α) = ∫ f (x) dx . |
(2.1) |
α
В иженерной гидрологии часто используют понятие обеспеченности СВ (вероятности превышения). Обеспеченность гидрологической величины Р(х) – это вероятность того, что рассматриваемое значение х гидрологической величины Х может быть превышено среди совокупности всех возможных ее значений
P(x) = p{X > x}=1− F(x) . |
(2.2) |
Обычно обеспеченность выражают не в долях единицы, как вероятностьСВ, а в процентах. Например, если среднегодовой расход имеет обеспеченность 1 %, это означает, что такой расход будет превышен в среднем один раз в 100лет.
Значения обеспеченности определяются для генетически однородных гидрологических характеристик: объемов весенних половодий, объемов дождевых паводков теплого периода, объемов годового стока, максимальных расходов половодий и максимальных расходов дождевых паводков, минимальных расходов межени и т.д.
3. Числовые характеристики случайных величин
Определение математического ожидания и дисперсии дискретной СВ:
n |
n |
|
mx = ∑ pi xi , Dx = ∑ pi (xi −mx )2 . |
(3.1) |
|
i=1 |
i=1 |
|
Математического ожидание (МО) СВ является главной характеристикой положения, дисперсия – рассеивания СВ.
Основные числовые характеристики положения непрерывной СВ – математическое ожидание mx (МО), медиана (Ме), мода (Мо) определяются так:
|
+∞ |
|
mx = |
∫x f (x) dx , |
(3.2) |
|
−∞ |
|
Me |
+∞ |
|
∫ f (x) dx = ∫ f (x) dx , |
(3.3) |
|
−∞ |
Me |
|
max (f (x))= f (Mo). |
(3.4) |
|
Числовые характеристики рассеивания непрерывной СВ: дисперсия, сред- |
||
нее квадратическое отклонение (СКО), коэффифиент вариации: |
|
|
Dx = µ2 = +∞∫(x −mx )2 f (x) dx , |
(3.5) |
|
|
−∞ |
|
11
σx = Dx , Cv = σx .
mx
Центральный момент s-го порядка
µS = M [(X −mx )S ].
(3.6)
(3.7)
Коэффициент асимметрии характеризует степень симметричности рассеяния относительно МО (см. рис. 3.1):
|
µ3 |
|
M [(X −mx )3 |
] |
|
1 |
+∞ |
3 |
|
|
|||
Cs = |
|
|
= |
|
|
= |
|
|
∫ |
(x −mx ) |
f (x) dx . |
(3.8) |
|
σ 3 |
σ 3 |
|
σ 3 |
||||||||||
|
|
x |
|
x |
|
|
|
x |
−∞ |
|
|
|
|
При симметричном распределении (Cs=0) Mo = Me = mx = a (рис. 3.1).
Рис. 3.1. Соотношение между МО, Ме и Мо для распределения с отрицательной Cs < 0 (а) и положительной Cs > 0 (б) асимметрией
Как уже отмечалось, каждая случайная величина исчерпывающим образом характеризуется своим законом распределения. Аналитическим выражением закона распределения является функция распределения (дифференциальная или интегральная), которую записывают в виде параметрического выражения, где в качестве параметров обычно используют числовые характеристики случайной величины (МО, СКО и т.д.). При такой записи для каждого значения СВ соответствующее значение функции распределения однозначно определяется параметрами аналитического выражения (параметрами распределения). В гидрологической практике, как правило, рассматриваются случайные величины, функции распределения которых зависят от небольшого числа параметров – обычно двух или трех.
12
В одних ситуациях закон распределения можно вывести из теоретических соображений, в других – он неизвестен, и аналитическое выражение функции распределения является лишь аппроксимацией истинного распределения. В последнем случае можно предложить несколько аналитических выражений. Исследователь вынужден производить выбор подходящего, опираясь на практический опыт и некоторые априорные соображения. Ниже рассматриваются аналитические функции распределения, наиболее часто используемые в практике гидрологических расчетов.
4. Равномерное и нормальное распределения случайных величин
Пусть возможные значения СВ X лежат в интервале [а, b] и нет ос нований для того, чтобы отдать предпочтение какому-либо из этих значений. При таких условиях СВ X подчиняется закону равномерной плотности (рис. 4.1). Это распределение называют также равномерным или прямоугольным.
а |
б |
Рис. 4.1. Плотность (а) и функция (б) равномерного распределения СВ
Закон равномерной плотности определяется двумя параметрами: началома и концом b интервала изменения СВ X. Дифференциальная и интегральная функции распределения определяются равенствами:
0 |
при x ≤ a; |
|
|
при a <x ≤b; |
(4.1) |
f (x) = 1/(b −a) |
||
|
при x >b. |
|
0 |
|
0 |
при x ≤ a; |
|
|
при a <x ≤ b; |
(4.2) |
F(x) = (x − a) /(b − a) |
||
|
при x > b. |
|
1 |
|
Подставляя (4.1) в (3.2) и (3.5) получим МО, дисперсию и СКО, соответственно:
mx = |
a +b |
; |
Dx = |
(b −a)2 |
; σx = |
(b −a) |
. |
(4.3) |
|||
2 |
12 |
|
|
|
|||||||
4 |
3 |
||||||||||
|
|
|
|
|
|
||||||
13
Медиана совпадает с МО. Закон равномерной плотности моды не имеет, так как все значения плотности вероятности равны между собой. Закон равномерной плотности играет важную роль при моделировании искусственных гидрологических рядов, поскольку значения любой случайной величины можно получить преобразованием значений СВ, равномерно распределенной на [0, 1].
В природе и различных областях человеческой деятельности весьма распространены СВ, представляющие собой сумму большого числа независимых или слабо зависимых СВ, дисперсии которых малы по сравнению с дисперсией всей суммы. Как следует из центральной предельной теоремы теории вероятностей, распределение таких СВ при весьма общих дополнительных условиях хорошо аппроксимируется нормальным распределением. Этим объясняется весьма широкое распространение последнего. Нормальное распределение применяется и в тех случаях, когда истинный закон распределения известен, но вычисления по нему затруднены, а аппроксимация его нормальным законом не приводит к большим ошибкам.
Плотность вероятности нормального распределения
f (x) = |
|
1 |
|
|
−(x −mx ) |
2 |
|
|
|
|
exp |
. |
(4.4) |
||||
|
|
|
|
|||||
σx |
|
2π |
2σx2 |
|||||
|
|
|
|
|
|
|||
Нормальное распределение является двухпараметрическим, т.е. зависит от двух параметров: математического ожидания mx и среднеквадратического отклонения σx. График плотности вероятности нормального распределения представлен на рис. 4.2. Как видно на рис. 4.2, нормальное распределение является симметричным и, следовательно, для него коэффициент асимметрии равен нулю (СS = 0), а мода, медиана и МО совпадают. Область возможных значений СВ, подчиняющейся нормальному распределению – вся числовая ось (−∞;+∞) .
Интегральная функция нормального распределения имеет вид
F(x) = |
|
1 |
|
x |
|
−(z −m2 |
x ) |
2 |
dz . |
|
|
|
|
∫ |
exp |
|
(4.5) |
||||||
|
|
|
|
||||||||
σx |
|
2π |
|
||||||||
|
|
|
−∞ |
|
2σx |
|
|
|
|
||
Интеграл, стоящий в правой части выражения (4.5), нельзя выразить через элементарные функции, поэтому ординаты нормального закона распределения обычно представляют в виде таблиц. Чтобы не публиковать множество таблиц, для различных сочетаний тх и σх вместо СВ X рассматривают нормированную СВ t, значения которой связаны со значениями X формулой
t = (x −mx ) /σx . |
(4.6) |
В настоящее время нет необходимости пользоваться таблицами, так как интеграл в (4.5) может быть легко найден в любом современном математическом пакете, например в среде Mathcad:
f (x) := dnorm(x,m,σ ) F(x) := pnorm(x,m,σ )
14
