Технофашист
225 / 213 / 11
Регистрация: 11.03.2009
Сообщений: 886
1

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

12.11.2010, 21:08. Показов 30930. Ответов 6
Метки нет (Все метки)

Какнить можно заставить rand() работать по аналогии в Паскале, т.е. выдавать числа от 0 до 1 (0,000000054, 0,743325235, 0,0052411 и подобно)???
Или может библиотеку подскажите, или может функция у кого есть. Буду благодарен

p.s. строго говоря, нужна случайная велечина с равномерным распределением на [0;1]
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
12.11.2010, 21:08
Ответы с готовыми решениями:

Генерация случайных чисел в диапазоне 0 - 1
Привет. Прошу помочь с рандомными числами. Нужны числа от 0 или 1

Генерация случайных чисел в заданном диапазоне.
Здесь сказано Почему?

Генерация случайных чисел на заданном диапазоне
Здравствуйте, прошу прощения, если такая тема уже проскакивала. Есть простой генератор чисел...

Рандом: генерация случайных чисел в заданном диапазоне
Вобщем такая задача: пользователь должен вводить верхнюю границу диапазона случайных чисел...

6
Эксперт С++
5052 / 3113 / 271
Регистрация: 11.11.2009
Сообщений: 7,045
12.11.2010, 21:35 2
srand((unsigned)time(NULL));
x = (double)rand() / RAND_MAX;

не забываем подключить stdlib.h и time.h
5
Эксперт С++
3064 / 1404 / 423
Регистрация: 19.01.2009
Сообщений: 3,816
12.11.2010, 21:39 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;
}
 
// ///////////////////////////////////////////////////////////////////////// //
1
Технофашист
225 / 213 / 11
Регистрация: 11.03.2009
Сообщений: 886
12.11.2010, 22:11  [ТС] 4
silent_1991, сори, но твой метод работает ужасно.

Oleg, спасибо за код, подходит. Но в Паскале разброс получше всётаки.
0
Эксперт С++
5052 / 3113 / 271
Регистрация: 11.11.2009
Сообщений: 7,045
12.11.2010, 22:25 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...
0
2831 / 1640 / 254
Регистрация: 03.12.2007
Сообщений: 4,222
12.11.2010, 23:01 6
Цитата Сообщение от darkAngel Посмотреть сообщение
заставить rand() работать по аналогии в Паскале
Цитата Сообщение от darkAngel Посмотреть сообщение
нужна случайная велечина с равномерным распределением на [0;1]
А в Паскале на [0; 1). Так что внимательнее - делить на RAND_MAX + 1 или на RAND_MAX.
0
Эксперт С++
2343 / 1716 / 148
Регистрация: 06.03.2009
Сообщений: 3,675
13.11.2010, 00:40 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/... 50591c333c
2
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
13.11.2010, 00:40
Помогаю со студенческими работами здесь

Генератор случайных вещественных чисел в пользовательском диапазоне
Добрый день, подскажите пожалуйста. Как написать правильный генератор вещественных чисел (в...

Генерация случайных чисел типа double во всем диапазоне
Нужна помощь в написании программы, генерирующей псевдослучайные числа типа double на всем...

Генерация вещественных случайных чисел с 3-мя знаками после запятой
Добрый день!хочу генерировать числа от -100 до 100 (вещественные) с тремя знаками после запятой...

Бинарные файлы. Написать функцию, которая записывает в бинарный файл N случайных вещественных чисел в указанном диапазоне
Добрый день! Подскажите пожалуйста! Генерирую рандомное кол-во вещественных чисел, и записываю...


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

Или воспользуйтесь поиском по форуму:
7
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2022, CyberForum.ru