Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса

Представлено аналитическую модель, описывающую симметричные логарифмически экспоненциальные распределения концентраций гидрохимических компонент. Рассчитаны погрешности оценивания вероятностей превышения ПДК по модели Гаусса, построены диаграммы погрешностей. Установлен диапазон значений эксцесса ра...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2011
Автор: Ковальчук, Л.А.
Формат: Стаття
Мова:Russian
Опубліковано: Український науково-дослідний гідрометеорологічний інститут МНС та НАН України 2011
Назва видання:Наукові праці Українського науково-дослідного гідрометеорологічного інституту
Теми:
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/51544
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса / Л.А. Ковальчук // Наукові праці Українського науково-дослідного гідрометеорологічного інституту: Зб. наук. пр. — 2007. — Вип. 256. — С. 308-317. — Бібліогр.: 1 назв. — рос.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-51544
record_format dspace
spelling irk-123456789-515442013-12-02T03:11:14Z Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса Ковальчук, Л.А. Моніторинг і стан довкілля Представлено аналитическую модель, описывающую симметричные логарифмически экспоненциальные распределения концентраций гидрохимических компонент. Рассчитаны погрешности оценивания вероятностей превышения ПДК по модели Гаусса, построены диаграммы погрешностей. Установлен диапазон значений эксцесса распределений, при которых вероятности превышения ПДК, оценённые по модели Гаусса с учётом погрешностей оценивания, приемлемы для решения практических задач. Представлено аналітичну модель, що описує симетричні логарифмічно експонентні розподіли концентрацій гідрохімічних компонент. Розраховані похибки оцінювання ймовірностей перевищення ГДК по моделі Гауса, побудовані діаграми похибок. Встановлено діапазон значень ексцесу розподілів, при яких імовірності перевищення ГДК, оцінені по моделі Гауса с урахуванням похибок оцінювання, прийнятні для рішення практичних задач. 2011 Article Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса / Л.А. Ковальчук // Наукові праці Українського науково-дослідного гідрометеорологічного інституту: Зб. наук. пр. — 2007. — Вип. 256. — С. 308-317. — Бібліогр.: 1 назв. — рос. XXXX-0054 http://dspace.nbuv.gov.ua/handle/123456789/51544 556.531:556.18 ru Наукові праці Українського науково-дослідного гідрометеорологічного інституту Український науково-дослідний гідрометеорологічний інститут МНС та НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language Russian
topic Моніторинг і стан довкілля
Моніторинг і стан довкілля
spellingShingle Моніторинг і стан довкілля
Моніторинг і стан довкілля
Ковальчук, Л.А.
Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса
Наукові праці Українського науково-дослідного гідрометеорологічного інституту
description Представлено аналитическую модель, описывающую симметричные логарифмически экспоненциальные распределения концентраций гидрохимических компонент. Рассчитаны погрешности оценивания вероятностей превышения ПДК по модели Гаусса, построены диаграммы погрешностей. Установлен диапазон значений эксцесса распределений, при которых вероятности превышения ПДК, оценённые по модели Гаусса с учётом погрешностей оценивания, приемлемы для решения практических задач.
format Article
author Ковальчук, Л.А.
author_facet Ковальчук, Л.А.
author_sort Ковальчук, Л.А.
title Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса
title_short Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса
title_full Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса
title_fullStr Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса
title_full_unstemmed Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса
title_sort погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели гаусса
publisher Український науково-дослідний гідрометеорологічний інститут МНС та НАН України
publishDate 2011
topic_facet Моніторинг і стан довкілля
url http://dspace.nbuv.gov.ua/handle/123456789/51544
citation_txt Погрешности интеграла вероятностей статистического распределения отличного от нормального, но вычисленного по модели Гаусса / Л.А. Ковальчук // Наукові праці Українського науково-дослідного гідрометеорологічного інституту: Зб. наук. пр. — 2007. — Вип. 256. — С. 308-317. — Бібліогр.: 1 назв. — рос.
series Наукові праці Українського науково-дослідного гідрометеорологічного інституту
work_keys_str_mv AT kovalʹčukla pogrešnostiintegralaveroâtnostejstatističeskogoraspredeleniâotličnogootnormalʹnogonovyčislennogopomodeligaussa
first_indexed 2025-07-04T13:37:14Z
last_indexed 2025-07-04T13:37:14Z
_version_ 1836723716827906048
fulltext УДК 556.531:556.18:311 Л.А. Ковальчук ПОГРЕШНОСТИ ИНТЕГРАЛА ВЕРОЯТНОСТЕЙ СТАТИСТИЧЕСКОГО РАСПРЕДЕЛЕНИЯ ОТЛИЧНОГО ОТ НОРМАЛЬНОГО, НО ВЫЧИСЛЕННОГО ПО МОДЕЛИ ГАУССА Представлено аналитическую модель, описывающую симметричные логарифмически экспоненциальные распределения концентраций гидрохимических компонент. Рассчитаны погрешности оценивания вероятностей превышения ПДК по модели Гаусса, построены диаграммы погрешностей. Установлен диапазон значений эксцесса распределений, при которых вероятности превышения ПДК, оценённые по модели Гаусса с учётом погрешностей оценивания, приемлемы для решения практических задач. Состояние проблемы Достоверность оценивания вероятностей превышения предельно допустимых концентраций (ПДК) гидрохимических компонент посредством интеграла вероятностей зависит от законов статистических распределений концентраций. При гидрохимических исследованиях часто принимается гипотеза о нормальном или логарифмически нормальном распределении концентраций, что не всегда подтверждается результатами статистического анализа. Процедуры идентификации законов статистических распределений требуют привлечения специального математического аппарата и связаны со значительными объёмами вычислений. В случаях, когда анализируемые распределения близки к нормальным распределениям, оценки статистических параметров гидрохимических компонент, в частности вероятности превышения ПДК, выполняются по правилам нормального распределения, но при этом не указываются возникающие погрешности оценивания. Цель настоящего исследования – установить допустимые пределы отклонений эксцесса от нормального значения, при которых погрешности интеграла вероятностей и вероятностей превышения ПДК, оценённых по нормальному закону с учётом погрешностей оценивания, приемлемы для решения практических задач. 308 Наук. праці УкрНДГМІ, 2007, Вип. 256 Материалы и методика В основу исследования положены материалы сезонных наблюдений Государственной гидрометеорологической службы Украины за последние три десятилетия на пятидесяти стационарных пунктах в основных речных бассейнах Днепра, Дуная, Днестра, Южного Буга, Западного Буга, Сиверского Донца, Приазовья и Крыма. Анализу были подвергнуты временные ряды сульфатов (SO4 2- ), хлоридов (Cl-), кальция (Са2+), магния (Mg2+), натрия (Na+), железа (Feобщ. ), меди (Cu ), цинка (Zn ), хрома (Cr6+ ), аммония (NH4 +), нитратов (NO2 -), нитритов (NO3 -), фосфатных ионов, кремния, фенолов, нефтепродуктов, синтетических поверхностно активных веществ (СПАВ). Установлено, что гистограммы статистических распределений логарифмированных значений концентраций зачастую имеют куполообразный вид. Методическая часть включала четыре этапа. На первом этапе, для идентификации законов статистических распределений концентраций, была привлечена аналитическая модель, описывающая широкий класс симметричных экспоненциальных распределений [1]: ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ − −= λσαλσ α α Хx Г xp цexp )/1(2 )( , (1) где )/3( )/1( α αλ Г Г = , Г(z) – гамма–функция: ( dttezГ zt 1 0 )( −∫ ∞ −= , при Re > 0), σ – стандартное отклонение, Хц – координата центра распределения, α – некоторая характеристика постоянная для анализируемого статистического распределения, которая однозначно определяет параметр формы экспоненциальных распределений – эксцесс (ε): [ ])/3(/)/5()/1( 2αααε ГГГ= , (2) При α < 1 модель описывает распределения, близкие по своим свойствам распределению Коши; α = 1 модель (1) представляет распределение Лапласа, имеющее более островершинный вид, чем у нормального распределения; α = 2 модель соответствует нормальному распределению; α > 2 модель идентифицирует распределения, близкие по своим свойствам к трапецеидальным распределениям; α → ∞ модель выражает равномерное распределение. 309 Приняв условие α = 2, Хц = 0, σ = 1, модель (1) была преобразована в модель нормального распределения: ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎣ ⎡ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ −= 22 1 2 exp)( xxp π , (3) посредством которой были рассчитаны значения плотности вероятностей распределения от 0 до 6. Для этого распределение центрировалось: Хц = 0 и представлялось к вычислениям в виде отклонений от центра через 0,1 от -3 до +3, как аналог ранжированной последовательности при построении гистограмм распределений. Значения плотности вероятностей интегрировались методом трапеций, в результате чего были получены значения функции распределения (интеграл вероятностей). На втором этапе был установлен диапазон значений эксцесса временных рядов гидрохимических компонент. Значения эксцесса были пересчитаны посредством выражения (2) в значения параметра α (табл. 1). Сохраняя условия: Хц = 0, σ = 1 и меняя значения α согласно табл. 1, посредством модели (1) были вычислены с дискретностью 0,1 значения кривых распределения, плотности вероятностей и интегралы вероятностей различных законов статистических распределений от -3 до +3. Таблица 1 Эксцессы распределений и соответствующие им значения α, использованные при численных экспериментах Эксцесс Параметр α Эксцесс Параметр α 1,900000 15,326920 3,761722 1,500000 2,000000 8,114994 4,000000 1,383831 2,631091 2,500000 5,000000 1,130523 2,787330 2,250000 6,000000 0,987010 3,000000 2,000030 7,000000 0,892513 3,302482 1,750000 8,000000 0,824588 На третьем этапе были вычислены разницы между значениями функции нормального распределения и значениями функций распределений, отличных от распределения Гаусса, когда α ≠ 2, но вычисленных по модели Гаусса. Полученные разницы были пересчитаны в относительные погрешности функции нормального распределения и представлены на диаграммах в процентах. 310 На четвёртом этапе для характерных временных рядов концентраций гидрохимических компонент посредством модели (1) были идентифицированы законы статистических распределений концентраций и вычислены вероятности превышения ПДК. Затем по модели Гаусса для тех же временных рядов были вычислены вероятности превышения ПДК, которые корректировались при помощи диаграмм и сопоставлялись с истинными значениями вероятностей превышения ПДК. Обсуждение результатов исследования При расчётах по модели Гаусса интегралов вероятностей пологих распределений со значениями эксцессов от 2,78 до 3,0, возникают как положительные, так и отрицательные погрешности: положительные – когда значения функции распределения Гаусса оказываются больше значений функций распределений идентифицированных статистических законов; отрицательные – когда функция распределения Гаусса меньше фактической функции распределения (рис. 1). В диапазоне вероятностей от 0 до 5% относительные погрешности расчётов стремительно уменьшаются с 25% до нуля; при вероятностях от 5 до 50% погрешности становятся отрицательными, не опускаясь ниже 5%, и остаются близкими к нулю при вероятностях более 50%. В случае распределений со значениями эксцессов от 3 до 3,3, более островершинных чем кривая Гаусса, погрешности интегралов вероятностей, рассчитанные по модели Гаусса, в диапазоне вероятностей 0-5% имеют отрицательный знак и стремительно уменьшаются от 30% до 0%. Затем, при вероятностях от 5 до 50%, погрешности приобретают положительные значения и по абсолютной величине не превышают 5%, а при вероятностях более 50% – колеблются около 0% (рис. 1). В диапазоне значений эксцесса 2,63-2,79 и 3,30-4,53 относительные погрешности интегралов вероятностей, вычисленных посредством модели Гаусса, достигают 50-80% при минимальных вероятностях. Тем не менее, сохраняются общие черты погрешностей: они максимальные при вероятностях от 0 до 5-6%, изменяют свой знак в диапазоне вероятностей 6-60%, и незначительны при вероятностях более 60% (рис. 2). 311 -30 -20 -10 0 10 20 30 0, 19 0, 82 2, 44 5, 94 12 ,4 22 ,5 36 ,2 51 ,9 67 ,3 80 ,1 89 ,3 94 ,9 97 ,9 99 ,2 99 ,7 значения функции нормального распределения (%) от но си те ль ны е по гр еш но ст и (% ) 2.79 2.91 3.11 3.3 Рис. 1. Погрешности интегралов вероятностей распределений со значениями эксцесса 2,79, 2,91, 3,11, 3,30, оценённых по модели Гаусса В диапазонах значений эксцесса 1,90-2,63 и 4,53-8,00 сохраняются общие черты относительных погрешностей интегралов вероятностей, вычисленных по модели Гаусса (рис. 3). Однако их значения существенно увеличиваются. При трапецеидальных распределениях с эксцессом около 2 и значениями интеграла вероятностей менее 4% относительные погрешности достигают +100%, опускаются до -40% в диапазоне вероятностей 10-50%, снова поднимаются до +10% при вероятностях 312 около 70-80%, после чего плавно уменьшаются до незначительных величин. -100 -80 -60 -40 -20 0 20 40 60 0, 19 0, 82 2, 44 5, 94 12 ,4 22 ,5 36 ,2 51 ,9 67 ,3 80 ,1 89 ,3 94 ,9 97 ,9 99 ,2 99 ,7 значения функции нормального распределения (%) от но си те ль ны е по гр еш но ст и (% ) 2.63 2.79 3.3 3.76 4 4.53 Рис. 2. Погрешности интегралов вероятностей распределений со значениями эксцесса 2,63, 2,79, 3,30, 3,76, 4,0 и 4,53, оценённых по модели Гаусса 313 -120 -100 -80 -60 -40 -20 0 20 40 60 80 100 120 0, 19 0, 82 2, 44 5, 94 12 ,4 22 ,5 36 ,2 51 ,9 67 ,3 80 ,1 89 ,3 94 ,9 97 ,9 99 ,2 99 ,7 значения функции нормального распределения (%) от но си те ль ны е по гр еш но ст и (% ) 1.9 2 2.63 4.53 5 6 7 8 Рис. 3. Погрешности интегралов вероятностей распределений со значениями эксцесса 1,90, 2,00, 2,63, 4,53, 5,00, 6,00, 7,00, 8,00, оценённых по модели Гаусса В диапазоне значений эксцесса 4,53-8,00 распределение погрешностей однородно: при малых вероятностях от 0 до 5% отрицательные значения погрешностей уменьшаются от 80-100% до 0%, в диапазоне вероятностей 5-50% погрешности увеличиваются до +20-+40%, при вероятностях более 50% погрешности сначала снижаются до -15%, а затем уменьшаются до 2-5%. 314 Таблица 2 Сопоставление вероятностей превышения ПДК фосфатных ионов, вычисленных по фактическим логарифмически экспоненциальным моделям и по моделям Гаусса с последующим учётом погрешностей оценивания Местоположение пункта наблюдений Вероятности превышения ПДК (%), вычислен. по Корректировка вероятн. превыш. ПДК, вычисленной по модели Гаусса (%) Бас- сейн Река Пункт Экс- цесс факти- ческой моде- ли моде- ли Гаус- са поп- равка исправ- ленное значе- ние По- грешн. коррек- тир-ки 1 2 3 4 5 6 7 8 9 Дунай Прут Яремча 3,86 10,81 11,14 -1,31 9,83 0,98 Западн. Буг Луга Вл.-Во- лынский 4,21 19,78 24,87 -3,93 20,94 -1,16 Закар- патье Лато- рица Подпо- лозье 2,19 33,22 26,82 3,48 30,30 -2,92 Днестр Зол. Липа Бережа- ны 4,68 83,74 81,03 4,86 85,89 -2,15 При- пять Ту- рья Ковель 5,26 5,91 6,89 1,3091 5,58 -0,33 Днепр Хо- рол Мирго- род 3,22 91,15 90,09 0,09 90,18 -0,97 Сив. Донец Уды Харьков 4,97 98,61 99,35 -0,49 98,85 0,24 Южн. Буг Ин- гул Кирово- град 3,18 93,13 92,93 0,09 92,84 -0,29 Приа- зовье Каль- чик Креме- нёвка 3,05 29,37 28,19 0,06 28,13 -1,24 Крым Сал- гир Пионер- ское 1,96 30,25 23,82 3,79 27,61 -2,64 Сопоставление вероятностей превышения ПДК фосфатных ионов в поверхностных водах Украины, вычисленных по фактическим логарифмически экспоненциальным моделям и по моделям Гаусса, с последующим учётом погрешностей оценивания показало, что наибольшие различия, достигающие 2,64% и 2,92% характерны для трапецеидальных распределений со значениями эксцессов 1,96 и 2,19 (табл. 2, рис. 4 а, б). 315 а) 0,00 0,05 0,10 0,15 0,20 0,25 0,30 0,35 -5, 52 -5, 25 -4, 97 -4, 70 -4, 42 -4, 15 -3, 87 -3, 60 -3, 32 -3, 05 -2, 78 -2, 50 -2, 23 -1, 95 -1, 68 -1, 40 ln (концентрации фосфатов) пл от но ст ь ве ро ят но ст и б) 0,00 0,10 0,20 0,30 0,40 -5 ,1 2 -4 ,9 1 -4 ,7 1 -4 ,5 0 -4 ,3 0 -4 ,1 0 -3 ,8 9 -3 ,6 9 -3 ,4 8 -3 ,2 8 -3 ,0 8 -2 ,8 7 -2 ,6 7 -2 ,4 6 -2 ,2 6 -2 ,0 6 ln (концентрации фосфатов) пл от но ст ь ве ро ят но ст и в) 0 0,2 0,4 0,6 0,8 1 1,2 1,4 -4, 61 -4, 42 -4, 24 -4, 05 -3, 87 -3, 68 -3, 50 -3, 31 -3, 13 -2, 94 -2, 76 -2, 57 -2, 39 -2, 20 -2, 02 -1, 83 ln (концентрации фосфатов) пл от но ст ь ве ро ят но ст и Рис. 4. Плотность вероятностей логарифмов концентрации фосфатных ионов: а) р. Латорица, c. Подполозье, эксцесс 2,19; б) р. Салгир, с. Пионерское, эксцесс 1,96; в) значения р(х) кривой распределения логарифмов концентрации фосфатных ионов р. Турья, г. Ковель, эксцесс 5,26 В тоже время обострённость кривой распределения со значением эксцесса 5,26 не всегда влечёт за собой значительные различия в оценках вероятностей превышения ПДК (табл. 2, рис. 4 в). 316 Выводы 1. Статистические распределения логарифмов концентраций гидрохимических компонент поверхностных вод Украины характеризуются значениями эксцесса от 1,90 до 8,00, что охватывает диапазон симметричных логарифмически экспоненциальных распределений. 2. Если не требуется высокая точность расчёта, то оценки вероятностей превышения ПДК гидрохимических компонент допустимо проводить согласно логарифмически нормальной модели с учётом погрешностей оценивания посредством предложенных диаграмм. 3. При необходимости высокоточного оценивания вероятностей превышения ПДК гидрохимических компонент рекомендуется использовать аналитическую модель, описывающую весь диапазон симметричных экспоненциальных распределений. * * Представлено аналітичну модель, що описує симетричні логарифмічно експонентні розподіли концентрацій гідрохімічних компонент. Розраховані похибки оцінювання ймовірностей перевищення ГДК по моделі Гауса, побудовані діаграми похибок. Встановлено діапазон значень ексцесу розподілів, при яких імовірності перевищення ГДК, оцінені по моделі Гауса с урахуванням похибок оцінювання, прийнятні для рішення практичних задач. * * 1. Новицкий П.В., Зограф И.А. Оценка погрешностей результатов измерений. – Л.: Энергоатомиздат, Ленингр. отд-ние, 1985. – 248 с. Украинский научно-исследовательский гидрометеорологический институт, Киев 317