Файл: Яковлев, В. В. Стохастические вычислительные машины.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 15.10.2024
Просмотров: 93
Скачиваний: 0
На рис. 26 кривая т = 1 представляет зависимость nD (zt) —
— р 2 (1 — р 2) в функции р — математического ожидания входной переменной квадратора, а кривая т = 2 — зависимость диспер сии, вычисленной по формуле (2.38). Первая функция достигает
максимума D x = 0,25 при р = 1/0,5; значение максимума второй «функции D 2 = 0,46 и определено при р = 0,728. Отношение
Рис. 26. Зависимость дисперсии частоты появле ния единиц в выходной последовательности дли ной п от р при возведении переменной в т-ю степень
дисперсий D JD X= 1,84. Следовательно, относительное увеличе ние длины случайной двоичной последовательности Бернулли равно 1,84, то есть п' = 1,84ге, что, очевидно, приведет к необхо димости повышения разрядности выходного накопителя СтВМ.
На рис. 27 представлена схема с двумя стохастическими умно жителями, позволяющая вычислить третью степень от входной переменной. Вероятность появления импульса на выходе схемы равна р 3, а дисперсия p s (1 — р 3). Вероятность совместного на ступления Событий Zi = 1 и zUx = 1
Р — Р fa ) Р (х{ |xt) р fa .j) Р (ж£_г |х{_х) р (xt_2) р (х,+1) = р4,
•54
а вероятность |
Р (ztzi+2) = р ъ. |
Ясно, |
что Р (ztz,) |
= |
р 6 для всех: |
||
) > 2. |
|
|
|
|
|
|
|
Корреляционная матрица в этом случае имеет |
вид |
||||||
р3(1—р 8) pi ( i —p 2) |
|
p 5q |
О |
|
О |
||
|
ps ( i —p3) |
pi ( i —p2) p bq ■■■ |
о |
||||
|
|
|
|
|
|
|
P3( l —P3) . |
Теперь, используя |
уравнение (2.36), получим |
|
|
||||
|
D ( т |
) = 1 Г (1 + |
2Р + |
2р2 - 5^3)- |
|
|
|
Рис. 27. Схема для возведения сто |
|
|
|
|
|
||
хастической переменной |
в третью |
|
|
|
|
|
|
степень |
|
|
|
|
|
|
|
|
|
|
|
|
т< |
|
т, |
На рис. 26 эта |
зависимость в |
функции р представлена кривой: |
|||||
т = 3. Видно, |
что |
отношение |
D sJDx я» 2,7. |
Следовательно, |
во столько же раз должна быть увеличена длина входной последо вательности, если мы хотим получить на выходе схемы значение ошибки не превышающее то, которое образуется при интегриро вании последовательности Бернулли.
Наконец, для произвольной степени т Ф 1 и re |
> |
т , а тг > |
тс ,. |
||
можно записать |
[74] |
|
|
|
|
|
|
т-1 |
|
|
|
nD |
1 + |
2 ^ р 1- { 2 т - \ ) р т |
(2.39) |
||
|
|
=i |
|
|
|
Значения максимума Dm этой |
функции зависят |
от |
т (табл. |
4). |
В третьей строке таблицы даны значения р т — второй координаты максимума.
|
|
|
|
|
|
|
|
Т а б л и ц а 4 |
|
Предельные |
значения дисперсии процесса на выходе стохастического |
||||||||
|
|
|
|
умножителя |
|
|
|
|
|
т |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
Dm |
0,46 |
0,67 |
0,9 |
1,12 |
1,34 |
1,56 |
1,78 |
2,01 |
2,23 |
Рт |
0,73 |
0,81 |
0,86 |
0,88 |
0,9 |
0,92 |
0,93 |
0,934 |
0,94 |
5 а
Отношение n'Jn получим из таблицы, умножив соответствующее Dm на величину D j1 = 4, т. е.
п" = inD m. |
(2.40) |
В качестве элементов задержки входящей последовательности можно использовать триггеры и тогда практическая схема для воз ведения переменной в т-ю степень будет состоять из т — 1-разряд ного регистра сдвига и такого же количества двухвходовых конъюнкторов, включенных цепочкой, как показано на рис. 28. Число вентилей И, используемых в схеме, в общем является избыточ ным и, следовательно, может быть сокращено. Хорошей иллюстра цией этому служит схема на рис. 27.
Рг
Рис. 28. Схема с одним регистром сдвига Рг для возведения переменной в произвольную степень т
Выпишем значения выходной функции схемы в последовательдше моменты времени с условием т х = т 2 = 1:
zi = XiXi.jXi.yXi.^
Z;+l x i+ ix ix ix i-y,
zn— XnXn^yXn_yXn_2-
Исключая тавтологии и переходя к вероятностям, при условии идеальности входящей последовательности получим
Р (zd = Р (xt) р (хи1) р (х£_?)
или, подставляя |
р (х£) = р (х^у) = р (х£_ 2) |
— • • • = p t окон |
чательно р (д,) = |
р 3. Заметим, что эта схема |
реализует операцию |
возведения переменной в третью степень и в том случае, если вто рой элемент задержки подключается к выходу первого конъюнкчора (эта связь показана на рисунке штрихами).
Если тА= 1, а т2 = 2, то
Z; = х (Xi ^\Х£.о'Г/-к-
Z;+1 — Xi+yX£Xi „уХ/_2,
Zn ---х пх п-\х п-2х п-Ь
и
Таким образом, схема на рис. 27 содержит на один вентиль меньше,, чем схема на рис. 28, когда они обе работают в режиме возведения переменной в четвертую степень.
В связи с этим представляет интерес разработка алгоритмов умножения многих переменных или возведения переменной в це лую степень на стохастических умножителях, приводящих к мини мальным затратам оборудования.
Запишем выражения этих интересующих нас операций:
* 01Х < А з ' ■' -^-0mi |
|
|
||||
Х0Х0Х0 . . . X 0= Xf , |
|
(2.41} |
||||
X di vd2 |
X т |
|
|
|||
|
0 1 Л 02 |
|
|
|||
|
|
|
• • ■Л- от> |
|
|
|
где все показатели и т — вещественные и целые числа. |
Х о1 яв |
|||||
Машинными отображениями истинных |
переменных |
|||||
ляются соответствующие вероятности р (х,-). |
Если область |
машин |
||||
ных переменных р (хг) |
£ (0, |
1) |
совпадает с диапазоном изменения |
|||
истинных переменных |
X oi, |
то соотношения |
|
|
||
Р Ы Р (х2) р (х а) . . . р (хт), |
| |
|
||||
Р (х) р {х )р {х ) |
. . . |
р (х) = рт(х), |
(2.42) |
|||
Ы |
pd‘ (х2) |
. . , |
pdm(xm) |
J |
|
фактически определяют алгоритмы вычисления функций системы (2.41). Следовательно, первое уравнение системы (2.42) может быть реализовано на некотором количестве N 0 г-входовых венти
лей И. Если г = |
т, то используется один такой элемент. |
В силу того, |
что стохастические элементы с раздельными вхо |
дами не имеют памяти, первая задача (умножение т различных переменных) решается за один шаг вычислений.
При этом количество r-входовых элементов И, требуемое для
вычисления произведения т различных переменных, |
равно [73] |
Г т — 1 ~ |
(2.43) |
= L j— 1 J 7 |
где символ [а] обозначает число а взятое по избытку до ближай шего целого.
Доказательство можно провести индукцией по г. При г = 2 каждый элемент выполняет одну операцию умножения. Поэтому общее число элементов определяется числом операций, необходи мых для умножения т переменных, т. е. N 0 = т — 1. При г = 3 каждый элемент выполняет две операции умножения, т. е. потре буется вдвое меньшее число элементов, и т. д. Ясно, что при г > 2 для некоторых т не все N 0r входов логической схемы окажутся занятыми. Для определения количества избыточных входов L 0 .
удобно |
воспользоваться выводом признаков делимости чисел |
т — 1 |
на г — 1 [9]. |
5 7 ;
Из |
|
|
|
to- 1 |
= &s10s-1 + &s.110s”2+ |
. . . + ^ - 1 , |
|
замечая, что 10 ~ |
О (mod 2), имеем |
|
|
|
т — 1 = J Ъ1— 1 1 (mod 2). |
|
|
■Следовательно, Ъг — 1 = t + 2q, где |
q — целое, а 0 ^ t |
2, |
откуда вытекает, что при перемножении нечетного числа перемен ных на трехвходовых элементах используются все входы; если т — четное, то не используется один вход.
Замечая, что 10 |
1 (mod 3), имеем |
т,— 1 = |
bs -f- bs_t + . . . -f Ъ1 — 1 (mod3). |
Рис, 29. Зависимости количества |
двухвходовых |
||
(— ) и |
трехвходовых (—) стохастических умножи |
||
телей, |
требуемого для умножения т различных пе |
||
ременных (1) и возведения переменной в натураль |
|||
|
ную степень т (2) |
|
|
Представим правую часть сравнения в виде |
|||
|
fes + frs- i + |
. ••+&i — 1 = |
3 ? !-f-£lt |
где q t — целое, |
а 0 ^ |
<( 3. Тогда L 0 — |
Следовательно, при перемножении т переменных на четырех входовых элементах все входы будут использованы, если сумма цифр, изображающих т, минус единица кратна трем. Например,
если т = 23, то 3 qx + jtx = 3 -7 + 2, откуда |
L Q= 2. |
Рассуждая аналогично, можно определить |
величину L 0 для |
иных г.
На рис. 29 показаны зависимости N 0 от т при использовании двухвходовых и трехвходовых стохастических умножителей в ре жиме перемножения т переменных.
Теперь перейдем к рассмотрению умножителей, содержащих на входах развязывающие триггеры (рис. 25, 27).
5 8
Представим т в системе счисления с основанием г
m = bs+1rs + bsrs~1-{- . . . + & ! •
Тогда минимальное число r-входовых стохастических умно жителей, с помощью которых можно вычислить р т{х), равно
|
S+1 |
bi~- s ( r — D — 1 |
|
|
2 |
|
|
#0 = |
;=1 |
0. |
(2.44). |
|
Действительно, из разложения т по степеням г вытекает, что-
сомножитель вида rs |
реализуется на s г-входовых элементах, со |
множитель вида г*-1 |
— на s — 1 элементах и так далее. Иными |
словами, все целые |
степени г разложения вычисляются на S'- |
r-входовых элементах. |
|
Теперь остается задача о суммировании |
|
|
и — bs+i + bs + . . . + &i |
чисел, среди которых могут быть |
и равные: fcs + 1 чисел равных rs,. |
bs чисел равных rs-1 и т. д. С |
рассматриваемой точки зрения |
задача сложения и чисел на стохастических умножителях ничем |
|
не отличается от аналогичной задачи умножения и чисел. Поэтому,, |
применяя формулу (2.43), получим
что и дает искомый результат (2.44).
На рис. 29 представлены зависимости iV0 от т, вычисленные* по формуле (2.44) для случаев использования двухвходовых и трехвходовых умножителей. Отметим, что для некоторых т (4, 8, 16 и т. д.) затраты оборудования одинаковы.
Определим разницу в оборудовании при вычислениях по фор мулам (2.43) и (2.44)
s+1
т —^ bi— s (r—1)
Разложив т по степеням г, получим
ANo = j~ |
fe»+i 9 s- l ) |
+ M r s- 1- i ) + |
• • • + Ы г - 1 ) - « ( г - ч ) |
j ^ (2.45) |
|||
после чего исследуем числитель этой дроби. |
|
|
|||||
При s = |
0 AiV0 = |
0, |
т. е. формулы (2.43) и (2.44) |
дают одина |
|||
ковый |
результат. |
0 лишь в том случае, когда b2= |
|
|
|||
При s = |
1 ДУ0= |
1. |
Если же- |
||||
Ь2 > 1 , |
то |
A7V0 > 0 . |
Учитывая |
независимость ДУ 0 |
от Ь1г |
5 3