Індивідуальний комплект навчально-методичних матеріалів з дисципліни "Надійність та технічна діагностика електрообладнання", страница 6

При X0=0,5 получается множество случайных чисел Х={0,152545; 0,954155; 0,084227, ...}. Программа занимает малый объем памяти и ее можно без редактирования фрагментарно вставлять в другие программы.

Для получения генератора k случайных чисел с равновероятными исходами интервал (0...1) делится на k равных частей и организуется вывод на индикацию появления номера соответствующего подинтервала.

Программа 9/34. Генератор псевдослучайных событий с равномерными исходами.

0

1

2

3

4

5

6

7

8

9

0

ИП9

Fp

+

Fex

П9

КИП9

XY

ИП9

-

П9

1

ИП8

X

1

+

П7

КИП7

ИП7

С/П

БП

00

Инструкция: k=P8; x0=Р9; В/О; С/П; РХ = y1; С/П; РХ = у2, ..., С/П; РХ = уn. В этой программе использована операция взятия целого от числа путем следующей подпрограммы

0

1

2

3

4

5

6

7

8

9

0

П7

КИП7

ИП7

С/П

БП

00

Инструкция: х = РХ; В/О; С/П; Int(х)= Р(х). При статистических испытаниях сложных систем возникает потребность в получении случайной перестановки К чисел, имитирующих случайную последовательность отказов элементов. С этой целью можно использовать алгоритм «перемешивания» массива.

Программа 10/34. Перемешивание массива чисел.

0

1

2

3

4

5

6

7

8

9

0

П0

ИПД

Fp

+

Fex

ПД

КИПД

XY

ИПД

-

1

ПД

ИПС

Х

1

+

ПА

ИП0

ПВ

КИПВ

КИПА

2

КПВ

XY

КПА

FLO

01

С/П

Инструкция: RАV=РД; x1=Pl, ...; xk = PK; K=PC; k = РХ; В/О; С/П; PI = x1, ..., PK = xk.

2.2. Генераторы псевдослучайных чисел с заданным законом распределения.

Эти генераторы используются для получения случайного времени безотказной работы элемента получают принимая случайное число RAV за вероятность безотказной работы RAV=R(t), откуда определяют t через обратную функцию R-1(t).

Например, при экспоненциальном законе отказов

RAV = ехр (-lt); t=-1/l ln (RAV).

Программа 11/34. Генерация псевдослучайных чисел с экспоненциальным законом распределения.

0

1

2

3

4

5

6

7

8

9

0

ИП9

Fp

+

Fex

П9

КИП9

XY

ИП9

-

П9

1

Fln

ИПА

¸

/-/

С/П

БП

00

Инструкция: RAV=P9; l=РА; В/О; С/П; t1 = РХ; С/П; t2 = РХ, ... .

Этот подход целесообразно использовать при простом виде обратной функции от R(t).

2.3. Статистические оценки числовых характеристик случайных величин.

           Поскольку всякие статистические оценки связаны с обработкой большого числа данных, очень важно использовать эффективные алгоритмы вычисления, некоторые из них приведены в данном разделе.

Оценка средней арифметической (математическое ожидание) выполняется согласно выражению

где ti—значение времени безотказной работы і-го элемента;

N—число испытуемых элементов.

Стандартная ошибка при этом составит , где D- дисперсия.

Оценка дисперсии осуществляется по формуле

Программа 12/34. Проверка генератора псевдослучайных чисел.

0

1

2

3

4

5

6

7

8

9

0

ИП9

Fp

+

Fex

П9

КИП9

XY

ИП9

-

П9

1

П1

ИПА

+

ПА

ИП1

Fx2

ИПВ

=

ПВ

ИП2

2

1

+

П2

FLO

00

ИПВ

ИПА

Fx2

ИП2

¸

3

-

ИП2

1

-

¸

С/П

ИПА

ИП2

¸

С/П