Интересное распределение случайных величин
В аппроксимации я, как говорится, собаку съел. Но тут крепко споткнулся. Никакая формула безупречно не совпадала с давдцатью девятью точками экспериментов! То макушка оказывалась другой, то хвост не так сближался с осью абсцисс. Стал усложнять формулы, то есть добавлять дополнительные параметры, и это не помогало. Неделю провозился, забросил даже русфорус ру. В конце концов начал наобум писать разные математические структуры и оптимизировать параметры методом Монте-Карло. И неожиданно одна из них дала сумму квадратов отклонений на три порядка меньше, чем оптимальное нормальное распределение! Оказалось - всего лишь однопараметрическое выражение в точности соответствовало всем линиям изгиба. Конкретно, такая аппроксимация плотности вероятности [Рис. b)]:
Тут же звоню Диме (то есть моему работодателю) и диктую формулу. Он попросил взять интеграл правой ветви колокола. Но не просто численно, а алгебраически. То есть вместо числа в кубе желательно иметь "a" в кубе. Теперь другая головная боль - интеграл. В принципе Вольфрам его брал, но уж очень корявый ответ. Некрасивый какой-то. Плюнул на Вольфрам и, внимательно следуя рекомендациям хороших учебников, стал методично разворачивать решение. В итоге получил красоту неописуемую! Можете посмотреть сами [Рис. c)]:
Верхний знак дает правую ветвь колокола, и нижний знак - левую ветвь.
Такое не стыдно предлагать как очередной табличный интеграл! Сколько студентов будет избавлено от мучений и нервных переживаний! Я тут же вписал его в Википедию. Однако прошла пара дней, и этот замечательный интеграл изъяли. Мол, нужен известный литературный источник. Надежный такой, общепризнанный. Ну как с такими недалекими разговаривать? Хочется, как лучше, а они делают полное обнуление.
Но это я так, отвлекся немного. Дальше было вот что. Нужно же взять определенный интеграл. Потеть пришлось всю ночь, но совершенно не зря! На Рис. d) виден невероятно компактный результат! И после подстановки с большой точностью вычислилась ожидаемая (в теории функции плотности вероятности) единица! Конечно приближенная единица, но очень-очень близкая. Конкретно 1.0000004. И тут я (опять же по просьбе Димы) уточнил сам параметр "а", при котором площадь под кривой окажется абсолютно точной единицей:
a=2*sqrt(pi/(3*sqrt(3)))
Наглядно формула эта - на Рис. f).
Абсолютизация параметра "а" настолько покорила коллегу, что тот обещал мне выслать ящик коньяка. Ко дню космонавтики и дню моего рождения.
Догадайтесь с трех раз - получу ли я столь бесценный приз?
24 января 2021 г.
Свидетельство о публикации №221012401405