Экспоненциальное распределение
May. 10th, 2012 07:29 pmРешил похвастаться получившимися результатами. Я не только реализовал это, но и ещё сделал полную проверку. Нет, конечно я лукавлю, по критерию Пирсона или
я не проверял, но на глаз всё хорошо.

Гистограмма распределения и проверочная кривая экспоненты
На картинке приведена гистограмма распределения для 1000 000 сгенерированных таким образом случайных чисел.
Выражаю благодарность товарищу
eddy_em за наводку на свой пост http://eddy-em.livejournal.com/2573.html и функцию drand48(). В общем-то за основу и была взята програма, там уже приляпал функцию построения гистограммы и проверки.
Свой быдлокод приводить не буду, скажу, что получить экспоненциальное распределение очень просто:
rand_number=-lambda*(log(drand48()));
Где lambda - это интенсивность потока (отсылаю к википедии http://ru.wikipedia.org/wiki/%D0%AD%D0%BA%D1%81%D0%BF%D0%BE%D0%BD%D0%B5%D0%BD%D1%86%D0%B8%D0%B0%D0%BB%D1%8C%D0%BD%D0%BE%D0%B5_%D1%80%D0%B0%D1%81%D0%BF%D1%80%D0%B5%D0%B4%D0%B5%D0%BB%D0%B5%D0%BD%D0%B8%D0%B5)
я не проверял, но на глаз всё хорошо.
Гистограмма распределения и проверочная кривая экспоненты
На картинке приведена гистограмма распределения для 1000 000 сгенерированных таким образом случайных чисел.
Выражаю благодарность товарищу
Свой быдлокод приводить не буду, скажу, что получить экспоненциальное распределение очень просто:
rand_number=-lambda*(log(drand48()));
Где lambda - это интенсивность потока (отсылаю к википедии http://ru.wikipedia.org/wiki/%D0%AD%D0%BA%D1%81%D0%BF%D0%BE%D0%BD%D0%B5%D0%BD%D1%86%D0%B8%D0%B0%D0%BB%D1%8C%D0%BD%D0%BE%D0%B5_%D1%80%D0%B0%D1%81%D0%BF%D1%80%D0%B5%D0%B4%D0%B5%D0%BB%D0%B5%D0%BD%D0%B8%D0%B5)