Форум программистов, компьютерный форум CyberForum.ru

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

Восстановить пароль Регистрация
 
 
Рейтинг: Рейтинг темы: голосов - 73, средняя оценка - 4.82
slavikg
0 / 0 / 0
Регистрация: 03.04.2012
Сообщений: 71
03.04.2012, 00:11     Моделирование случайных чисел с экспоненциальным законом распределения #1
Программирования случайных чисел с экспоненциальным законом распределения. Товарищи помогите пожалуйста.
Лучшие ответы (1)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
03.04.2012, 00:11     Моделирование случайных чисел с экспоненциальным законом распределения
Посмотрите здесь:

Сгенерировать последовательность 60 случайных чисел с экспоненциальным законом распределения с параметром C++
ГСЧ для равномерного распределения случайных чисел на заданном интервале C++
Написать датчик случайных чисел с таблично заданной функцией распределения C++
[C++] Генератор случайных чисел... закон распределения задан C++
C++ Моделирование случайных величин
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Kuzia domovenok
 Аватар для Kuzia domovenok
1883 / 1738 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
03.04.2012, 00:52     Моделирование случайных чисел с экспоненциальным законом распределения #2
Сообщение было отмечено автором темы, экспертом или модератором как ответ
Логарифмическая функция распределения
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
0 / 0 / 0
Регистрация: 03.04.2012
Сообщений: 71
03.04.2012, 17:13  [ТС]     Моделирование случайных чисел с экспоненциальным законом распределения #3
както слишком просто)

Добавлено через 4 минуты
Компилятор на всё ето говорит ерор
Kuzia domovenok
 Аватар для Kuzia domovenok
1883 / 1738 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
03.04.2012, 17:32     Моделирование случайных чисел с экспоненциальным законом распределения #4
Ну наверное надо ещё добавить, какой еррор.
slavikg
0 / 0 / 0
Регистрация: 03.04.2012
Сообщений: 71
03.04.2012, 17:36  [ТС]     Моделирование случайных чисел с экспоненциальным законом распределения #5
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
 Аватар для Kuzia domovenok
1883 / 1738 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
03.04.2012, 18:16     Моделирование случайных чисел с экспоненциальным законом распределения #6
Так ты это не тупо копируй в программу.
Я думал ты умеешь хоть что-то программировать.
Когда ты заговорил о числах с экспоненциальным законом распределения, я подумал: ну уж наверное он знает, что такое случайные числа и может написать в программе простейший генератор равномерно распределённых на интервале чисел.
А это оказывается твоя первая программа! Брось её и пиши 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
0 / 0 / 0
Регистрация: 03.04.2012
Сообщений: 71
03.04.2012, 18:29  [ТС]     Моделирование случайных чисел с экспоненциальным законом распределения #7
Очень болшое спасибо. Я уже умею немного, просто с генератором перый раз столкнулся.
alexey31415
 Аватар для alexey31415
59 / 59 / 3
Регистрация: 16.05.2010
Сообщений: 632
03.04.2012, 18:57     Моделирование случайных чисел с экспоненциальным законом распределения #8
Цитата Сообщение от Kuzia domovenok Посмотреть сообщение
#define LOFFSET ln(RAND_MAX)
double getLrand(double l){
* *return (ln(rand()+1)-LOFFSET)/(-l)
RAND_MAX - это символическая константа или определяемое нами число?а то не сталкивался с RAND_MAX прежде
slavikg
0 / 0 / 0
Регистрация: 03.04.2012
Сообщений: 71
03.04.2012, 19:56  [ТС]     Моделирование случайных чисел с экспоненциальным законом распределения #9
а ln ето логарифм натуральный? А то мне ошыбку выдает
Kuzia domovenok
 Аватар для Kuzia domovenok
1883 / 1738 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
03.04.2012, 20:01     Моделирование случайных чисел с экспоненциальным законом распределения #10
Цитата Сообщение от alexey31415 Посмотреть сообщение
RAND_MAX - это символическая константа или определяемое нами число?а то не сталкивался с RAND_MAX прежде
из stdlib оттуда же, откуда и rand

Добавлено через 3 минуты
Цитата Сообщение от slavikg Посмотреть сообщение
а ln ето логарифм натуральный? А то мне ошыбку выдает
Нет такой функции ln, надо писать log вместо ln. Это уж я лоханулся.
slavikg
0 / 0 / 0
Регистрация: 03.04.2012
Сообщений: 71
03.04.2012, 20:09  [ТС]     Моделирование случайных чисел с экспоненциальным законом распределения #11
если я делаю замену на log то мне выбивает
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.cpp(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.cpp(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.cpp(13): error C3861: printf: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.cpp(15): error C3861: scanf: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.cpp(16): error C3861: printf: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.cpp(18): error C3861: scanf: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.cpp(21): warning C4244: аргумент: преобразование "time_t" в "unsigned int", возможна потеря данных
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.cpp(25): error C3861: printf: идентификатор не найден
1>c:\users\славiк\documents\visual studio 2010\projects\slavi16\slavi16\slavik16.cpp(29): error C3861: getch: идентификатор не найден
panicwassano
591 / 559 / 20
Регистрация: 07.11.2010
Сообщений: 2,004
03.04.2012, 20:14     Моделирование случайных чисел с экспоненциальным законом распределения #12
так надо подключить библиотеки
stdio.h
conio.h
math.h
slavikg
0 / 0 / 0
Регистрация: 03.04.2012
Сообщений: 71
03.04.2012, 20:19  [ТС]     Моделирование случайных чисел с экспоненциальным законом распределения #13
подключены
Kuzia domovenok
 Аватар для Kuzia domovenok
1883 / 1738 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
04.04.2012, 00:18     Моделирование случайных чисел с экспоненциальным законом распределения #14
Вот такой вариант однозначно работает
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
342 / 287 / 10
Регистрация: 23.03.2012
Сообщений: 838
04.04.2012, 00:50     Моделирование случайных чисел с экспоненциальным законом распределения #15
Цитата Сообщение от 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);
http://upload.wikimedia.org/wikipedi...7bc3280b86.png
slavikg
0 / 0 / 0
Регистрация: 03.04.2012
Сообщений: 71
04.04.2012, 16:00  [ТС]     Моделирование случайных чисел с экспоненциальным законом распределения #16
А почему мне всегда видает результат -0,0000 ?
Nekto
342 / 287 / 10
Регистрация: 23.03.2012
Сообщений: 838
04.04.2012, 16:36     Моделирование случайных чисел с экспоненциальным законом распределения #17
Цитата Сообщение от 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
 Аватар для alexey31415
59 / 59 / 3
Регистрация: 16.05.2010
Сообщений: 632
04.04.2012, 19:36     Моделирование случайных чисел с экспоненциальным законом распределения #18
Цитата Сообщение от Nekto Посмотреть сообщение
(float)RAND_MAX
а разве и так можно привести переменную к типу float?
slavikg
0 / 0 / 0
Регистрация: 03.04.2012
Сообщений: 71
04.04.2012, 21:35  [ТС]     Моделирование случайных чисел с экспоненциальным законом распределения #19
а какой вариант правилный тот который ты дал или тамтой что выше ? Ато результати разные
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
04.04.2012, 21:36     Моделирование случайных чисел с экспоненциальным законом распределения
Еще ссылки по теме:

Моделирование случайной величины с заданным законом распределения C++
C++ Написание своего генератора случайных чисел с различными плотностями распределения
Как подобрать последовательность 100 случайных чисел с экспоненциальным законом на C++ C++

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
Nekto
342 / 287 / 10
Регистрация: 23.03.2012
Сообщений: 838
04.04.2012, 21:36     Моделирование случайных чисел с экспоненциальным законом распределения #20
Цитата Сообщение от slavikg Посмотреть сообщение
а какой вариант правилный тот который ты дал или тамтой что выше ? Ато результати разные
последний. В том я формулу не ту задал.
Yandex
Объявления
04.04.2012, 21:36     Моделирование случайных чисел с экспоненциальным законом распределения
Ответ Создать тему
Опции темы

Текущее время: 20:22. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru