При 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 |
¸ |
С/П |
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.