«Статическое моделирование Лабораторная работа №1»

Задание

1. Построить генератор псевдослучайных чисел, равномерно

распределенных на интервале (0,1), на основе следующих соотношений(генератор Лемера):

Здесь:

− некоторые вспомогательные величины(реальные!);

γ − псевдослучайные числа, равномерно распределенные на

интервале (0,1);

знак(mod 4) означает, что берется остаток от деления на 4

2. Построить гистограмму распределения чисел γ на интервале (0,1) при .

3. На основе генератора из п.1 получить 1000 случайных цифр по формуле

4. Вычислить частоты каждой цифры. Вычислить значение

Исходный код

string str = "37542048056489474296248052403720636104020842268953196450930323209025601595334764990190252909376707153831131165886767439712807999708015736147640323665398951168776606574717340727685036697361706581339885";

double[] ga = new double[10010];

double[] U= new double[10010];

double[] E = new double[10010];

double[] C = new double[10010];

double[] p = new double[10010];

double[] r = new double[10010];

int k; int k2; int k3;

double q;

static double Factorial(double n)

{

if (n == 0)

return 1;

else

return n * Factorial(n - 1);

}

private void button1_Click(object sender, EventArgs e)

{

dataGridView1.Rows.Clear();

k=Convert.ToInt32(textBox1.Text);

U[0]=3.14;

U[1] = (Math.Pow(5, 17)) * (Math.Pow(2, -42));

textBox2.Text = Convert.ToString(U[1]);

for (int i = 2; i

Random rnd = new Random();

private void button4_Click(object sender, EventArgs e)

{

dataGridView3.Rows.Clear();

dataGridView4.Rows.Clear();

double n= Convert.ToInt16(textBox5.Text);

double p1 = Convert.ToDouble(textBox6.Text);

double la = Convert.ToDouble(textBox7.Text);

k3 = Convert.ToInt32(textBox4.Text);

double M; double P; double CC; int k;

for (int i = 0; i < k3; i++)

{

M = rnd.Next(99);

M = M / 100;

dataGridView3.Rows.Add();

dataGridView3[0, i].Value = i + 1;

dataGridView3[1, i].Value = M;

k = 0;

CC = (Factorial(n)) / (Factorial(k) * Factorial(n - k));

q = 1 - p1;

p[k] = CC * (Math.Pow(p1, k)) * (Math.Pow(q, n - k));

P = p[k];

if ((M - P) >= 0)

while (M >= 0)

{

CC = (Factorial(n)) / (Factorial(k) * Factorial(n - k));

p[k] = CC * (Math.Pow(p1, k)) * (Math.Pow(q, n - k));

r[k] = ((n - k) / (k + 1)) * (p1 / (1 - p1));

P = P * r[k];

M = M - P;

k++;

}

dataGridView3[2, i].Value = k;