Аналитические и имитационные модели. Финаев В.И - 76 стр.

UptoLike

76
- NZ - заданное число тактов моделирования;
- m - число интервалов разбиения отрезка [0,1],
происходит обнуление необходимых для работы
идентификаторов и счетчиков.
В блоке 2 происходит наращивание тактов
моделирования. В блоке 3 датчиком случайных чисел
генерируется число Х. Затем в блоках 4 – 6 определяется
индекс I. В блоке 7 датчиком случайных чисел
генерируется второе число Х (тем самым образована
пара
(Х
1
,Х
2
)). В блоках 8 – 10 определяется индекс J.
В блоке 11 в счетчиках K[I,J] осуществляется подсчет
частот попадания случайных пар (Х
k
,Х
р
) в
соответствующие деления квадратной таблицы. В блоке 12
осуществляются проверки выполнения числа опытов. Если
датчик псевдослучайных, квазиравномерно
распределенных чисел выполнил генерацию заданного
числа опытов NZ/2, то в блоке 13 происходит вычисление
случайной величины
χ
2
. В подпрограмме WIWOD (см.
блок 14) на экран дисплея выводятся значения случайной
величины
χ
2
, счетчиков K[I,J], значения эмпирических
частостей
*
ij
P
=2n
ij
/N=2K[I,J]/NZ. Генерируемая выборка в
данном методе используется неэффективно. Можно пары
образовать в следующем виде (Х
1
,Х
2
),(Х
2
,Х
3
),(Х
3
,Х
4
),... .
Этот метод образования пар более эффективен, т.к.
полнее использует выборку чисел, но из-за зависимостей
пар случайная величина
χ
2
определится по формуле
=−
==
∑∑
22
mm
2
mNmN
2
χ
(n ) (n )
ij i
2
NNm
m
i,j 1 i 1
], (4.9)
где
=
=
m
nn
ii
j
j
1
.