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

Генерация случайных вещественных чисел в диапазоне 0..1 - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 72, средняя оценка - 4.88
darkAngel
Технофашист
206 / 187 / 3
Регистрация: 11.03.2009
Сообщений: 829
12.11.2010, 21:08     Генерация случайных вещественных чисел в диапазоне 0..1 #1
Какнить можно заставить rand() работать по аналогии в Паскале, т.е. выдавать числа от 0 до 1 (0,000000054, 0,743325235, 0,0052411 и подобно)???
Или может библиотеку подскажите, или может функция у кого есть. Буду благодарен

p.s. строго говоря, нужна случайная велечина с равномерным распределением на [0;1]
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
12.11.2010, 21:08     Генерация случайных вещественных чисел в диапазоне 0..1
Посмотрите здесь:

C++ Генерация случайных чисел
C++ Генерация случайных чисел типа double во всем диапазоне
Генерация случайных чисел в заданном диапазоне. C++
C++ Генерация случайных чисел на заданном диапазоне
C++ Генерация случайных чисел
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
silent_1991
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
12.11.2010, 21:35     Генерация случайных вещественных чисел в диапазоне 0..1 #2
srand((unsigned)time(NULL));
x = (double)rand() / RAND_MAX;

не забываем подключить stdlib.h и time.h
schdub
 Аватар для schdub
2902 / 1246 / 222
Регистрация: 19.01.2009
Сообщений: 3,210
Завершенные тесты: 1
12.11.2010, 21:39     Генерация случайных вещественных чисел в диапазоне 0..1 #3
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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
// ///////////////////////////////////////////////////////////////////////// //
//                                                                           //
// L'Ecuyer algorithm for uniform random generator with practically endless  // 
 
// period. Combining 2 sequences.                                            //
//                                                                           //
// ///////////////////////////////////////////////////////////////////////// //
 
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
 
// ///////////////////////////////////////////////////////////////////////// //
 
long dummy;
 
// ///////////////////////////////////////////////////////////////////////// //
 
void seed(long d)
 
{
 
    dummy = d;
 
}
 
// ///////////////////////////////////////////////////////////////////////// //
 
float rnd()
 
{
 
    const long IM1 = 2147483563;
 
    const long IM2 = 2147483399;
 
    const double AM = (1./IM1);
 
    const long IMM1 = (IM1 - 1);
 
    const int IA1 = 40014;
 
    const int IA2 = 40692;
 
    const long IQ1 = 53668;
 
    const long IQ2 = 52774;
 
    const int IR1 = 12211;
 
    const int IR2 = 3791;
 
    const int NTAB = 32;
 
    const int NWUP = 8;
 
    const double NDIV = (1 + IMM1 / NTAB);
 
    const double EPS = 1.2e-7;
 
    const double RNMX = (1.0 - EPS);
 
    
 
    static long dummy2 = 123456789;
 
    static long iy = 0;
 
    static long iv[32];
 
 
 
    long k = 0;
 
 
 
    // initialize the random sequence (first set of coefficients, the 
 
    // routine close to that in the function above
 
    if (dummy <= 0 || !iy)
 
    {
 
        // avoid negative or zero seed
 
        if (dummy < 0)
 
            dummy =- dummy;
 
        else
 
            if (!dummy)
 
                dummy = 1;
 
        dummy2 = dummy;
 
        // after NWUP warmups, initialize shuffle table
 
        for (int j = NTAB + NWUP - 1; j >= 0; --j)
 
        {
 
            k = dummy / IQ1;
 
            if ((dummy = IA1 * (dummy - k * IQ1) - IR1 * k) < 0) dummy += IM1;
 
            if (j < NTAB) iv[j] = dummy;
 
        }
 
        // first specimen from the table
 
        iy = iv[0];
 
    }
 
    // regular work: generate 2 sequences.
 
    k = dummy / IQ1;
 
    if ((dummy = IA1 * (dummy - k * IQ1) - IR1 * k) < 0) dummy+=IM1;
 
    k = dummy2 / IQ2;
 
    if ((dummy2 = IA2 * (dummy2 - k * IQ2) - IR2 * k) < 0) dummy2+=IM2;
 
    // shuffle output combining 2 sequences
 
    int j = 0;
 
    iy = iv[j = iy / NDIV] - dummy2;
 
    iv[j] = dummy;
 
    if(iy < 1) iy += IMM1;
 
    // return the result, as in the previous function
 
    float temp;
 
    return ((temp = AM * iy) > RNMX ? RNMX : temp);
 
}
 
// ///////////////////////////////////////////////////////////////////////// //
 
 
int main()
{
    seed(time(NULL));
    for (int i = 0; i < 100; ++i)
        printf("%f\n", rnd());
    return 0;
}
 
// ///////////////////////////////////////////////////////////////////////// //
darkAngel
Технофашист
206 / 187 / 3
Регистрация: 11.03.2009
Сообщений: 829
12.11.2010, 22:11  [ТС]     Генерация случайных вещественных чисел в диапазоне 0..1 #4
silent_1991, сори, но твой метод работает ужасно.

Oleg, спасибо за код, подходит. Но в Паскале разброс получше всётаки.
silent_1991
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
12.11.2010, 22:25     Генерация случайных вещественных чисел в диапазоне 0..1 #5
darkAngel, странно... Тысячи программистов им пользуются - и ничего, не жалуются...

Код
0.391034
0.552843
0.552721
0.340861
0.382183
0.233924
0.672872
0.109348
0.832606
0.187262
0.182318
0.622852
0.715262
0.62331
0.561632
0.604328
0.867641
0.510575
0.96231
0.927976
0.181249
0.132023
0.162847
0.875729
0.45732
0.643239
0.519333
0.411267
0.209449
0.788873
0.44087
0.878018
0.195196
0.433485
0.422132
0.0268563
0.957152
0.0631733
0.163427
0.780694
0.154942
0.914182
0.810297
0.641041
0.96997
0.884457
0.843043
0.143925
0.362987
0.884976
0.648305
0.270577
0.467086
0.786645
0.303323
0.219459
0.896237
0.61916
0.604328
0.446089
0.251228
0.177984
0.870479
0.376324
0.604511
0.695975
0.554216
0.700217
0.243934
0.735374
0.189673
0.651662
0.45848
0.957305
0.710318
0.181036
0.645589
0.923856
0.666524
0.426466
0.522233
0.590014
0.375439
0.390179
0.720786
0.641255
0.53328
0.226356
0.550798
0.332804
0.321818
0.0462355
0.208197
0.98764
0.540513
0.647481
0.212928
0.808283
0.794946
0.0730308
По-моему превосходный набор из сотни псевдослучайных чисел от 0 до 1...
Somebody
2770 / 1583 / 141
Регистрация: 03.12.2007
Сообщений: 4,139
Завершенные тесты: 1
12.11.2010, 23:01     Генерация случайных вещественных чисел в диапазоне 0..1 #6
Цитата Сообщение от darkAngel Посмотреть сообщение
заставить rand() работать по аналогии в Паскале
Цитата Сообщение от darkAngel Посмотреть сообщение
нужна случайная велечина с равномерным распределением на [0;1]
А в Паскале на [0; 1). Так что внимательнее - делить на RAND_MAX + 1 или на RAND_MAX.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
13.11.2010, 00:40     Генерация случайных вещественных чисел в диапазоне 0..1
Еще ссылки по теме:

Бинарные файлы. Написать функцию, которая записывает в бинарный файл N случайных вещественных чисел в указанном диапазоне C++
генерация рандомных чисел в диапазоне C++
Генерация случайных чисел в диапазоне 0 - 1 C++

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

Или воспользуйтесь поиском по форуму:
CyBOSSeR
Эксперт C++
 Аватар для CyBOSSeR
2293 / 1663 / 86
Регистрация: 06.03.2009
Сообщений: 3,675
13.11.2010, 00:40     Генерация случайных вещественных чисел в диапазоне 0..1 #7
darkAngel, решение с помощью boost::random:
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
#include <iostream>
#include <boost/random/mersenne_twister.hpp>
#include <boost/random/uniform_real.hpp>
 
int main() {
  boost::mt19937 generator;
  
  boost::uniform_real<> distribution(0.0, 1.0);
 
  const int count = 100;
  for(int i = 0; i < count; ++i)
     std::cout << distribution(generator) << std::endl;
     
  return 0;
}
Результат: http://liveworkspace.org/code/30fc55...13d750591c333c
Yandex
Объявления
13.11.2010, 00:40     Генерация случайных вещественных чисел в диапазоне 0..1
Ответ Создать тему
Опции темы

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