Файл: Яковлев, В. В. Стохастические вычислительные машины.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 15.10.2024
Просмотров: 97
Скачиваний: 0
а суммарную дисперсию частоты появления единиц на т-м выходе соединения
nD |
) - р (1 - р ) т [1 - р (1 - |
р)т (1 + |
2т)]. |
(2.53) |
Если исследуется только один элемент |
(т = |
1), то |
|
|
и |
nD (zu) = p ( i — р)(1 —р + р 2) |
|
(2.54) |
|
|
|
|
|
|
|
nD (± ) = р (1 - р) (1 - Зр + Зр2). |
(2 .5 5 ) |
Обратим внимание, что функция (2.55) дает в четыре раза мень шую дисперсию, чем определенную по (2.52) для элемента «эквиваленция».
Обе зависимости (2.54) и (2.55) в функции р представлены на рис. 31.
Вычитая из (2.55) (2.54) и записав результат
( n - l ) K Zi ( 1 ) = - р 2 ( 1 - р ) 2,
обнаружим, что в данном случае имеет место отрицательная кор реляция. По этой причине происходит уменьшение вариаций на выходе логического элемента «запрет» (или «импликация»). С ростом т это свойство проявляется сильнее, так как функция (2.53) в пределе
lira р (1 —р)т [1 —р (1 — р)т (1 + 2т )] = О
т-> со
для любых значений р.
10. Точность и скорость выполнения операций
Рассмотрим типичное соединение блоков СтВМ (рис. 33), пред назначенное для выполнения какой-либо математической операции над т машинными переменнымир (х^),р (х2), . . ,,р (хт). При этом, если истинные переменные Х 01, Х 02, . . ., Х 0т представлены в виде двоичных кодов, то используются т входных преобразователей
код — вероятность П К В , |
осуществляющих кодирование инфор |
мации по одной из принятых схем (ДЛС или ОЛС). |
|
Машинные переменные |
поступают на вход логического пре |
образователя Л П — дискретного устройства, являющегося ориен тированным т,1 полюсником без обратных связей. Блок Л П содержит логические элементы И, ИЛИ, Не и т. п., образующие однотактную схему, и линии задержки. К выходу этого блока под ключается интегратор для определения среднестатистического значения реализации выходной последовательности ЛП с матема тическим ожиданием р (z).
Можно отметить следующие основные факторы, влияющие на точность вычислений в структуре на рис. 33.
5 В . В. Яковлев |
65 |
1. Неидеальность передаточной функции преобразователей П И В. Ранее уже отмечалось (стр. 18), что линейное преобразова ние реализуется лишь в том случае, если распределение состояний Х ;- регистра Рг2 преобразователя (рис. 6) равномерное.
Генерирование этих состояний обычно происходит так, что на
вход каждого i-го разряда /-разрядного регистра |
сдвига Рг2 |
|||||
подаются |
случайные |
последовательности |
р (уд |
от |
источников |
|
случайного |
процесса. |
|
|
|
|
|
Следовательно, вероятность записи единицы в первый разряд |
||||||
регистра равна р (уд , |
во второй р |
(у2) и т. д. |
Как и ранее, предпо |
|||
лагаем, что последовательности р |
(у д , р (у2), . . ., |
р |
(у,) незави |
|||
симы и стационарны. |
Если событие уг понимать в смысле появле |
ния значения 1 , а у,- — как противоположное событие, то р (у,) =
Рис. 33. Последовательность соединения основ ных стохастических блоков
=q(yt) и плотности вероятностей этих случайных величин задаем
ввиде
/ (Уд = Р (Уд 6 (У1- 1 ) . / (уд = Р {уд s (уд ■
Таким образом, могут быть вычислены вероятности p f событий Ху, заключающихся в появлении конкретного состояния Ху
(табл. 5).
Продемонстрируем правила составления этой таблицы на примере последней строки.
Имеем X 2i = у гу2. . •*//• Плотность вероятности /независимых случайных величин равна произведению их одномерных плотно
стей вероятности, т. е.
i
/ (Х 2/) = / (уг, у2, . . . . уд = П Р (Уд 6 (Vi — 1).
£=1
откуда вероятность сложного события р;- определяется как
I |
1+ 80 |
р2/= ИтП |
j р (Уд б (у,- — 1) dyu |
Ео->-о £=1 1-е0
66
|
|
|
|
Т а б л и ц а 5 |
|
В ер о я тн о сти |
о б р а зо ва н и я случ ай н ы х |
двоичн ы х наборов |
|
Состояния Ху |
|
|
|
|
|
Двоичное |
|
Плотность вероятности |
Вероятности |
|
|
f (ХР |
состояний |
|
х / |
представ |
|
р/ |
|
ление |
|
|
||
|
|
|
1 |
|
X i |
000. . . 0 |
|
П р (У{) ь (уР |
Я (Ур Я (Уг) ■■ - Я (Ур |
|
|
|
1=1 |
|
|
|
г - i |
|
|
Х 2 |
0 0 0 ... 1 |
П Р (Ур б (Ур Р (УР б (У; — 1) |
Я (Ур Я (Уг) ■ ■ - Р(Ур |
|
|
|
i =1 |
|
|
|
|
г - з |
г |
|
Х„ |
о о о .. . и И Р (ур б (г/£) 1 1 р (ур б (г/у — 1) я (уР ■•-Р(У1-рР(ур |
|||
|
|
1=1 |
1=1-2 |
|
|
|
|
г |
|
X,/ |
1 1 1 .. .1 |
|
П Р (Ур б (yi 1) |
Р (Ур Р (Уг) ■ ■ •Р (Ур |
|
|
|
i=1 |
|
или, используя свойства дельта-функции, окончательно запишем
P2i = Р (Ур Р (Уг) ■••Р(Ур-
Подставляя все полученные значения р;- в формулу (1.22), получим
р (ж) = 2 р /.
/=1
Очевидно, что, если р (г/,) = 0,5, то р;- — 2~l = const, и вновь подтверждается уже известное р;- = s2-/. Если р (уг) Ф 0,5, то вероятности состояний из множества Х у будут отличаться. Напри мер, для I = 3 при условии р (уР = р (у 2) = р (у3) = р обра зуются четыре типа разновероятных состояний с вероятностями:
Pi = ?8. Pi = Pe = P i= P 2Q,
Рь = Р2= Р з= Р Я 2, Ра = Р3-
При этом, чем больше отклоняется р от значения 0,5, тем силь нее выражены искажения при линейном преобразовании. Таким образом, к датчикам случайных чисел должны быть предъявлены довольно жесткие требования по стабилизации выходной вероят ности с номиналом р (1) = 0,5.
Г>* |
6 7 |
2. Наличие корреляции и автокорреляции в последовательно стях на выходе всех П К В .
Мы уже отмечали, что при условии взаимной независимости входных последовательностей математическое ожидание последова тельности на выходе произвольной комбинационной схемы опре деляется лишь математическими ожиданиями входных стохастиче ских процессов. Если это условие не выполняется, то на точность выполнения большинства операций влияет не только точность уста новки математических ожиданий, но и величина взаимной корре ляционной связи между входными последовательностями.
Мерой статистической зависимости между входными перемен ными Х[ и Xj может служит значение нормированной взаимной корреляционной функции
Р ( X j X j ) — p ( X j ) р ( X j )
Р ( * i ) Р (ж/) Р ( x j ) Р ( X j )
Тогда, предположив, что точность установки математических ожиданий значительно выше необходимой точности выполнения операции, можно поставить вопрос о допустимой величине значе ния корреляционной функции rt;-. Так, при выполнении операции умножения двух переменных, используя (2.7), запишем
p f a = P fa) Р {Xj) + о- fa) a (xj) г{]-(0), |
(2.56) |
где |
. |
a fa) = V p fa ) p fa), а (x}) = V p (x^ p (x;).
Поскольку необходимо выполнить условие
p(z) — p fa ) p (x^ = a (xt) a (x}) rlf (0) s? едоп>
то |
необходимо, чтобы |
|
|
|
|
rij (°) |
£доп |
|
|
|
a ( x i ) о |
( x j ) ’ |
|
|
|
Правая часть этого неравенства |
максимальна, |
когда р (xL) = |
|
= Р fa) = Р (xj) = Р fa ) = 0,5. |
|
|
|
|
В |
этом случае |
|
|
|
|
П,- (0) ^ 4едоп |
(2.57) |
Например, если истинные переменные представлены /-разряд ными двоичными кодами и соответствующие им стохастические переменные установлены с предельной абсолютной погрешностью, равной половине единицы младшего разряда* то
гц(0) 4 * 2 - «н> = 2-0-i).
Так что при I = 10 Гц (0) ^ 0,002.
68
Соотношение (2.57), таким образом, определяет допустимую степень статистической зависимости между разрядами регистров Рг2 преобразователей IIK B i и ПКВ-Г
Аналогично можно определить допустимую величину внутриразрядной корреляции (или нормированной автокорреляционной функции гхх%(т)). При возведении в квадрат переменной х резуль
тат определяется выражением
Р (xxj = р2(X) + О2(;г) гхх^(т),
где х и хх —■события, заключающиеся в совместном появлении единиц в последовательности в моменты времени t и t -f- т. Эта зависимость практически не отличается от выражения (2.56), и, следовательно, повторяя подобные рассуждения, окончательно получим
Г^ ( т ) < 2- « - 1>.
3.Конечность времени интегрирования. Если интегрируется идеальный процесс без последействия (ЛП не содержит линий за держки), то для определения состоятельного объема выборки, при котором частота появлений события z отличается от его ве роятности менее чем на величину ех при заданной доверительной вероятности р л, может быть использована ранее полученная зави симость (1.25, а). Подставляя в эту формулу значение
находим
«оР(z) ч (z) ei
Эта функция максимальна, когда р (z) = q (z) — 0,5, следова тельно, необходимый объем выборки (или время интегрирования) должен быть
2
п ■
или в логарифмической форме
lgra = 2 (lg - Y ----- |
lg e ^ . |
(2.58) |
На рис. 34 представлены некоторые зависимости, рассчитанные по этой формуле для ряда значений доверительной вероятности рд. Из графиков, в частности, видно что для достижения точности 10% при доверии р л = 0,95 достаточна выборка длиной примерно 100 тактов, для точности 1 % уже требуется 10 000 тактов, а полу чить е! 0,1 % , можно, проведя приблизительно 1000000 испы таний.
69