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

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

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

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

@Kuzia domovenok · 03.04.2012, 00:52

Логарифмическая функция распределения
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);
}

@slavikg · 03.04.2012, 17:13 **[ТС]**

както слишком просто)

Добавлено через 4 минуты
Компилятор на всё ето говорит ерор

@Kuzia domovenok · 03.04.2012, 17:32

Ну наверное надо ещё добавить, какой еррор.

@slavikg · 03.04.2012, 17:36 **[ТС]**

1>c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(7): error C2601: getLrand: недопустимые локальные определения функций
1> c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(5): эта строка содержит "{", которая пока не имеет парной
1>c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(8): error C2065: RAND_MAX: необъявленный идентификатор
1>c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(8): error C3861: ln: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(8): error C3861: rand: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(8): error C3861: ln: идентификатор не найден

Добавлено через 1 минуту
И как ещо вывести на екран примерно 100 результатов?

@Kuzia domovenok · 03.04.2012, 18:16

Так ты это не тупо копируй в программу.
Я думал ты умеешь хоть что-то программировать.
Когда ты заговорил о числах с экспоненциальным законом распределения, я подумал: ну уж наверное он знает, что такое случайные числа и может написать в программе простейший генератор равномерно распределённых на интервале чисел.
А это оказывается твоя первая программа! Брось её и пиши Hello World! Приди к тому кто тебе дал задание и скажи "для меня это сложно, научи меня программировать сначала".
Ты хоть понимаешь, что если я напишу её целиком, у тебя умений от этого всё равно не прибавится?
Открой Шилдта или Страуструппа и читай!

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
29
30
31
32
33
34
35
//По теме: надо подключить библиотеки 
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>
//Далее пишем наше изобретение
#define LOFFSET ln(RAND_MAX)
double getLrand(double l){
   return (ln(rand()+1)-LOFFSET)/(-l);
}
 
//Затем надо написать главную функцию!
int main(){
  // В главной функции выводишь на экран текст
  printf("Enter size of random sample:");//"введите размер выборки"
  int N;//в этой переменной будем хранить размер выборки
  scanf("%d", &N);//считываем с клавиатуры десятичное целое число и запоминаем его в переменной N
  printf("Enter lambda of distribution");//"Введите лямбду"
  double lambda;//это параметр твоего логарифмического распределения
  scanf("%f", &lambda);//считываем с клавиатуры десятичную дробь и запоминаем её в lambda
  //мы не будем даже массивов использовать или условий проверять - сразу выводим выборку
 //правда для генерации псевдослучайных чисел надо сначало задать seed
  srand ( time(NULL) );
//Итак: цикл вывода
  while (N>0) {//начало цикла. повторения будут происходить пока оставшееся число повторений не 
                  //дойдёт до 0
       printf("%f ", getLrand(lambda));//генерируем элемент выборки и тут же выводим
       N=N-1;//цикл подходит к концу, нам его осталось повторять на один раз меньше, поэтому N-1
   }             ///конец цикла, повторяем всё с места, обозначенного "начало цикла"
 
  getch();//в конце работы программа останавливается и ждёт ввода с клавиатуры, чтобы
 //ты увидел результат перед выходом из программы
  return 0;//возврат  из программы.
}

@slavikg · 03.04.2012, 18:29 **[ТС]**

Очень болшое спасибо. Я уже умею немного, просто с генератором перый раз столкнулся.

@alexey31415 · 03.04.2012, 18:57

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

#define LOFFSET ln(RAND_MAX)
double getLrand(double l){
* *return (ln(rand()+1)-LOFFSET)/(-l)

RAND_MAX - это символическая константа или определяемое нами число?а то не сталкивался с RAND_MAX прежде

@slavikg · 03.04.2012, 19:56 **[ТС]**

а ln ето логарифм натуральный? А то мне ошыбку выдает

@Kuzia domovenok · 03.04.2012, 20:01

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

RAND_MAX - это символическая константа или определяемое нами число?а то не сталкивался с RAND_MAX прежде

из stdlib оттуда же, откуда и rand

Добавлено через 3 минуты

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

а ln ето логарифм натуральный? А то мне ошыбку выдает

Нет такой функции ln, надо писать log вместо ln. Это уж я лоханулся.

@slavikg · 03.04.2012, 20:09 **[ТС]**

если я делаю замену на log то мне выбивает
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(7): error C2668: log: неоднозначный вызов перегруженной функции
1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(575): может быть "long double log(long double)"
1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(527): или "float log(float)"
1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(120): или "double log(double)"
1> при попытке сопоставить список аргументов "(int)"
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(7): error C2668: log: неоднозначный вызов перегруженной функции
1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(575): может быть "long double log(long double)"
1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(527): или "float log(float)"
1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(120): или "double log(double)"
1> при попытке сопоставить список аргументов "(int)"
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(13): error C3861: printf: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(15): error C3861: scanf: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(16): error C3861: printf: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(18): error C3861: scanf: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(21): warning C4244: аргумент: преобразование "time_t" в "unsigned int", возможна потеря данных
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(25): error C3861: printf: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(29): error C3861: getch: идентификатор не найден

@panicwassano · 03.04.2012, 20:14

так надо подключить библиотеки
stdio.h
conio.h
math.h

@slavikg · 03.04.2012, 20:19 **[ТС]**

подключены

@Kuzia domovenok · 04.04.2012, 00:18

Вот такой вариант однозначно работает

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

Поставил float перед RAND_MAX. Вообще дефайны, макросы такая тема, в которых я ошибки ищу в последнюю очередь. Ведь компилятор ругается на совершенно другую строчку. А зря, иногда очень неочевидные вещи происходят.

@Nekto · 04.04.2012, 00:50

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

Вот такой вариант однозначно работает

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

Зачем несколько раз логарифм брать?

C++
1
-1.0/lambda+log(rand()*1.0/RAND_MAX);

@slavikg · 04.04.2012, 16:00 **[ТС]**

А почему мне всегда видает результат -0,0000 ?

@Nekto · 04.04.2012, 16:36

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

А почему мне всегда видает результат -0,0000 ?

потому что для double надо в считывании\вводе писать "%lf", а не "%f"

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;
}

@alexey31415 · 04.04.2012, 19:36

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

(float)RAND_MAX

а разве и так можно привести переменную к типу float?

@slavikg · 04.04.2012, 21:35 **[ТС]**

а какой вариант правилный тот который ты дал или тамтой что выше ? Ато результати разные

@Nekto · 04.04.2012, 21:36

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

а какой вариант правилный тот который ты дал или тамтой что выше ? Ато результати разные

последний. В том я формулу не ту задал.

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11680&d=1772460536 Одним из. . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .
SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .

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

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

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71
	03.04.2012, 17:13 [ТС]
	както слишком просто) Добавлено через 4 минуты Компилятор на всё ето говорит ерор 0

@Kuzia domovenok 4268 / 3327 / 926 Регистрация: 25.03.2012 Сообщений: 12,531 Записей в блоге: 1
	03.04.2012, 17:32
	Ну наверное надо ещё добавить, какой еррор. 1

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71
	03.04.2012, 17:36 [ТС]
	1>c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(7): error C2601: getLrand: недопустимые локальные определения функций 1> c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(5): эта строка содержит "{", которая пока не имеет парной 1>c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(8): error C2065: RAND_MAX: необъявленный идентификатор 1>c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(8): error C3861: ln: идентификатор не найден 1>c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(8): error C3861: rand: идентификатор не найден 1>c:\users\славiк\documents\visual studio 2010\projects\slavik15\slavik15\slavik15 .cpp(8): error C3861: ln: идентификатор не найден Добавлено через 1 минуту И как ещо вывести на екран примерно 100 результатов? 0

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71
	03.04.2012, 18:29 [ТС]
	Очень болшое спасибо. Я уже умею немного, просто с генератором перый раз столкнулся. 0

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71
	03.04.2012, 19:56 [ТС]
	а ln ето логарифм натуральный? А то мне ошыбку выдает 0

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71
	03.04.2012, 20:09 [ТС]
	если я делаю замену на log то мне выбивает 1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(7): error C2668: log: неоднозначный вызов перегруженной функции 1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(575): может быть "long double log(long double)" 1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(527): или "float log(float)" 1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(120): или "double log(double)" 1> при попытке сопоставить список аргументов "(int)" 1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(7): error C2668: log: неоднозначный вызов перегруженной функции 1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(575): может быть "long double log(long double)" 1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(527): или "float log(float)" 1> e:\program files (x86)\microsoft visual studio 10.0\vc\include\math.h(120): или "double log(double)" 1> при попытке сопоставить список аргументов "(int)" 1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(13): error C3861: printf: идентификатор не найден 1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(15): error C3861: scanf: идентификатор не найден 1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(16): error C3861: printf: идентификатор не найден 1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(18): error C3861: scanf: идентификатор не найден 1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(21): warning C4244: аргумент: преобразование "time_t" в "unsigned int", возможна потеря данных 1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(25): error C3861: printf: идентификатор не найден 1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.c pp(29): error C3861: getch: идентификатор не найден 0

@panicwassano 601 / 569 / 104 Регистрация: 07.11.2010 Сообщений: 2,004
	03.04.2012, 20:14
	так надо подключить библиотеки stdio.h conio.h math.h 0

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71
	03.04.2012, 20:19 [ТС]
	подключены 0

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71
	04.04.2012, 16:00 [ТС]
	А почему мне всегда видает результат -0,0000 ? 0

@slavikg 0 / 0 / 0 Регистрация: 03.04.2012 Сообщений: 71
	04.04.2012, 21:35 [ТС]
	а какой вариант правилный тот который ты дал или тамтой что выше ? Ато результати разные 0

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

Решение