Моделирование случайных чисел с экспоненциальным законом распределения

@slavikg · Регистрация: 03.04.2012

Студворк — интернет-сервис помощи студентам

Программирования случайных чисел с экспоненциальным законом распределения. Товарищи помогите пожалуйста.

@slavikg · 04.04.2012, 22:14 **[ТС]**

можеш кинуть а то я непонял какой ето последній?)

@Nekto · 04.04.2012, 22:32

Сообщение от slavikg

можеш кинуть а то я непонял какой ето последній?)

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>
 
double getLrand(double l)
{
return -log((double)(rand())/RAND_MAX)/l;
}
 
int main()
{
printf("Enter size of random sample:");
int N;
scanf("%d", &N);
printf("Enter lambda of distribution");
double lambda;
scanf("%lf", &lambda);
srand ( time(NULL) );
while (N>0)
{
 printf("%lf ", getLrand(lambda));
 N=N-1;
}            
getch();
return 0;
}

@slavikg · 04.04.2012, 23:08 **[ТС]**

спасибо

@Twi7ter · 18.12.2013, 23:38

Сообщение от Kuzia domovenok

Логарифмическая функция распределения
F=1-exp(-lx)
Тебе нужна обратная ей функция. То есть по оси F откладываем случайные числа от 0 до 1 с равномерным распределением, а соответствующие им значения на оси x будут случайными, распределёнными по нужному нам закону.

Функция распределения F=1-exp(-lx)
Обратная ей будет x=ln(1-F)/(-l)
В итоге получаем
F=rand()/RAND_MAX;//нецелочисленное деление
x=ln(1-F)/(-l)
1-F распределена точно так же, как F, поэтому
x=ln(F)/-l=ln( rand()/RAND_MAX )/-l=( ln(rand()) - ln(RANDMAX))/(-l)

Итог:

C
1
2
3
4
#define LOFFSET ln(RAND_MAX)
double getLrand(double l){
   return (ln(rand()+1)-LOFFSET)/(-l);
}

ДО сего поста не мог понять как числа генерить с заданным законом распределения,огромное спасибо)

@Sergey7777 · 21.11.2014, 22:57

А что такое lambda of distribution? Т.е. double lambda;? У наж же ведь экспоненциальный закон распределения? А не какая-то там лямбда.

@Kuzia domovenok · 22.11.2014, 10:47

Сообщение от Sergey7777

У наж же ведь экспоненциальный закон распределения?

а теперь напиши мне, какова плотность экспоненциального распределения? Или все формулы из головы вылетели?

@Sergey7777 · 22.11.2014, 15:53

Вот те формула L*exp^(L*x) при x>0, где L-лямбда, но в программе нет такой формулы! Там только такая есть F=1-exp(-lx).

Добавлено через 6 минут
Но F=1-exp(-lx) - это логарифмическая функция распределения, а у нас говорится про экспоненциальный закон распределения.

@Kuzia domovenok · 22.11.2014, 17:38

Sergey7777, https://en.wikipedia.org/wiki/... stribution

Новые блоги и статьи Все статьи Все блоги /
Транслитерация текста с кириллицы на латиницу (вариант №2) Maks 13.11.2026 Обновленная версия ранее созданной обработки. Добавлено: - Проверка исходного текста на кириллицу, через обращение к стандартной функции "ТолькоКириллицаВСтроке" из общего модуля. . .	Часы с термометром и праздничным календарем yariko 12.11.2026 Сегодня мы рассмотрим как можно сделать отличные многофункциональные часы. Итак, поставим конкретную задачу: разработать часы с календарем, двумя термодатчиками, выводом поздравительных сообщений. . .	Транслитерация текста с кириллицы на латиницу (вариант №1) Maks 11.11.2026 Решил облегчить жизнь себе и коллегам админам при создании электронных почтовых ящиков для пользователей, в результате написал обработку "Транслитерация текста" с кириллицы на латиницу, а в ней. . .	Книги по геймдизайну, левел-дизайну, игровым механикам и проектированию игрового опыта 8Observer8 08.11.2026 Содержание блога \| English Version 1. 2003 - 06 - Designing Virtual Worlds (1st Edition) Автор: Richard A. Bartle \| Объем: 768 стр. Ссылки: Amazon . . .
Родители токсичные абьюзеры по отношению к своим детям ругают их оскорбляют не любят zorxor 08.11.2026 Жить в атмосфере постоянных оскорблений, критики и холодности со стороны родителей — это тяжелое и болезненное испытание, которое оставляет глубокие эмоциональные раны. Токсичное и абьюзивное. . .	Как только я захожу к себе в свою комнату я сразу же с себя снимаю свои трусы. Ответ от Джемени ИИ. zorxor 08.11.2026 Джемени, как только я захожу к себе в свою комнату я сразу же с себя снимаю свои трусы. Многие люди, возвращаясь домой, сразу переодеваются в максимально свободную одежду или вовсе предпочитают. . .	Паттерн очереди действий/команд (Action / Command Queue). Система очереди действий и кат-сцен с SDL3 Callbacks, cglm и OpenGL 8Observer8 07.11.2026 Содержание блога \| English Version Руководство для начинающих по реализации паттерна очереди действий/ команд (Action / Command Queue) на C++ (C++20) с использованием официального SDL3 Callbacks. . .	Проверка полного ФИО в типовом справочнике Maks 07.11.2026 Решение реализовано в конфигурации КА2, на примере справочника "КонтактныеЛицаПартнеров". Задача: Запрет инициалов. Не допускаеть использование сокращений вида "И. И. ", "И. И. ", "И И" и. . .

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71

	Моделирование случайных чисел с экспоненциальным законом распределения 03.04.2012, 00:11. Показов 38152. Ответов 27 Метки нет (Все метки) Программирования случайных чисел с экспоненциальным законом распределения. Товарищи помогите пожалуйста. 0

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71
	04.04.2012, 22:14 [ТС]
	можеш кинуть а то я непонял какой ето последній?) 0

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71
	04.04.2012, 23:08 [ТС]
	спасибо 0

@Sergey7777 7 / 7 / 0 Регистрация: 03.12.2013 Сообщений: 44
	21.11.2014, 22:57
	А что такое lambda of distribution? Т.е. double lambda;? У наж же ведь экспоненциальный закон распределения? А не какая-то там лямбда. 0

@Sergey7777 7 / 7 / 0 Регистрация: 03.12.2013 Сообщений: 44
	22.11.2014, 15:53
	Вот те формула Lexp^(Lx) при x>0, где L-лямбда, но в программе нет такой формулы! Там только такая есть F=1-exp(-lx). Добавлено через 6 минут Но F=1-exp(-lx) - это логарифмическая функция распределения, а у нас говорится про экспоненциальный закон распределения. 0

@Kuzia domovenok 4268 / 3328 / 926 Регистрация: 25.03.2012 Сообщений: 12,545 Записей в блоге: 1
	22.11.2014, 17:38
	Sergey7777, https://en.wikipedia.org/wiki/... stribution 0