Метод Коробова. 
Получение случайных чисел на ЭВМ

Псевдослучайные числа вычисляются по формуле.

(6).

где; ;; i=1,2,…, n;

m — простое число вида (- тоже простое число);

n — число псевдослучайных точек.

Заданиe: Получить десять псевдослучайных точек, равномерно распределенных на интервале (0,1), методами Неймана, Лемера и Коробова. Необходимые для вычислений значения параметров указаны в вариантах задания (таблица).

Выполнение:

А) Метод середины квадрата (Неймана) Листинг:

случайный величина монте карло.

#include.

#include.

using namespace std;

double F (double x).

{ int k=1;

int celoe= pow (x, 2)*pow (10, 3*k);

double drobnoe = (celoe * pow (10, -2*k));

if (drobnoe>0).

{.

drobnoe-=floor (drobnoe);

}.

return drobnoe;

}.

int main ().

{ //Русский язык.

setlocale (LC_ALL," russian");

double x0; //начальное число.

int const nn=20; //максимальное количество элементов в массиве.

int n;

double x[nn+1]; //массив цифр

int y[nn+1]; //массив хранения полученных чисел.

cout<<" Генерация случайных чисел с помощью метода середины квадрата (Неймана)" <

cout<<" Сколько псевдослучайных точек вы бы хотели видеть на интервале (0,1)?" <<" n = «;

cin>>n;

cout<<" Введите начальное число x0″ <

cin>>x0;

x[0]=x0;

x[1]=F (x[0]);

for (int i=1;i.

{.

x[i]=F (x[i-1]);

} cout<

for (int i=1;i.

{.

cout<<

}.

return 0;

}.

Результат работы программы:

Б) Метод сравнений (Лемера) Листинг:

#include.

#include.

using namespace std;

int main ().

{ setlocale (LC_ALL," Russian");

cout<<" Генерация случайных чисел методом сравнений (Лемера)." <

int const n = 10;

double m0 = 1.6;

double g = pow (5,5);

double M = pow (2,15);

double y[n+1];

double m[n+1];

m[0]=m0;

for (int i = 1; i.

{.

m[i]=(g*m[i-1])/M;

m[i]-=floor (m[i]); //получим остаток (дробь).

// cout<<

y[i]=g*m[i-1]/M;

y[i]-=floor (y[i]);

cout<<

}.

return 0;

}.

Результат работы программы:

В) Метод Коробова Листинг:

#include.

#include.

using namespace std;

int main ().

{ setlocale (LC_ALL," Russian");

cout<<" Генерация случайных чисел методом Коробова." ;

int const m1 = 19;

int const n = 10;

int h0=1;

double m = 2*m1+1;

double r[n];

double y[n];

double h[n];

h[0]=h0;

r[1]=m/2*h[0];

for (int i = 0; i.

{.

r[i+1] = m/2*h[i];

h[i+1] = r[i+1] - m*((int)(r[i+1]/m));

y[i+1] = r[i+1]/m;

y[i+1]-=floor (y[i+1]);

cout<

}.

return 0;

}.

Результат работы программы:

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