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

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

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

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

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

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

@Nekto · 04.04.2012, 22:32

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

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

C++

#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

#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

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

@Kuzia domovenok 4064 / 3318 / 924 Регистрация: 25.03.2012 Сообщений: 12,493 Записей в блоге: 1
	22.11.2014, 10:47	26
	Сообщение от Sergey7777 У наж же ведь экспоненциальный закон распределения? а теперь напиши мне, какова плотность экспоненциального распределения? Или все формулы из головы вылетели? 0

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

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

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

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

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