Файл: Яковлев, В. В. Стохастические вычислительные машины.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 95

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

угольником, вершинами которого, кроме занятых клеток, могут быть и нули главной диагонали. Однако необходимо напомнить, что при синтезе такого ГПСЧ особого внимания требует вопрос, касающийся интервалов сдвига между двоичными последователь­ ностями, генерируемыми на выходах полусумматоров.

Рассмотрим пример, иллюстрирующий метод синтеза парал­

лельного ГПСЧ на основе двух регистров сдвига.

 

 

 

Пример. Пусть требуется реализовать операцию

А~В,

где А,

В — входные

переменные,

заданные

8-разрядными двоичными

 

 

 

числами.

Для

 

преобразования

этих

 

 

 

переменных

в стохастические перемен­

 

 

 

ные, необходимо

иметь

две

последо­

 

 

 

вательности

 

псевдослучайных

чисел

 

 

 

X

= (ж1?

ж2,

. . ., xs)

и

У =

{ух,

 

 

 

 

• •» Ув)-

 

 

 

 

 

 

 

 

 

 

Положим, что операция возведения

 

 

 

в квадрат выполняется с помощью

 

 

 

конъюнктора

с

элементом

задержки

 

 

 

входной случайной последовательности

 

 

 

на один такт. Тогда для точного вы­

Рис. 128. Таблица выходов

полнения операции необходимо, чтобы

псевдослучайные

числа

Xk,

 

Хк+ 1

параллельного ГПСЧ с од­

и

Yk были

статистически независимы.

ним

регистром сдвига

Очевидно, для этого нужно взять число

 

 

 

 

 

 

разрядов в регистрах сдвига

т +

п ^

3-8 = 24.

С другой

стороны, чтобы

обеспечить

погрешность

вычисления

не хуже

е =

2~а+1) =

 

2-9,

необходимо

генериро­

вать

последовательность псевдослучайных чисел

длиной

п ^

^ 22/ = 2й . Но для этого требуется, чтобы суммарное число разрядов в регистрах сдвига было т + п Ss 16.

Желая получить наиболее простое устройство, пойдем на компромиссное решение, выбрав в качестве генераторов исход­ ных последовательностей максимальной длины 11-разрядный и 7-разрядный регистры сдвига с одним полусумматором в цепи

обратной связи (фх (х)

= 1 -f х + х7,

ср2 (х) = 1 4-

х +

х11,

т + п = 18). Очевидно,

такое решение будет приемлемым,

если

максимальная ошибка операции, вызванная линейной

зависи­

мостью между псевдослучайными числами

Xk, Хк+ 1, Y k,

не пре­

высит величины в = 2~9.

 

 

 

 

Перейдем к выбору связей от разрядов регистров к полу­ сумматорам, представляющим разряды ГПСЧ. Для этого по­

строим таблицу выходов и заполним в ней

16

клеток

соответ­

ственно

суммарному количеству

разрядов

в

числах

X = (хг,

я 2, . .

я8)

и У =

(уг, у2, . . .,

у8). Поскольку т -f

п — 1 =

=

17 >

16, то выходы в таблице можно разместить таким обра­

зом, чтобы

числа

Xk, Хк+ 1 и Y k

были попарно независимыми.

Пример

такого размещения показан на рис. 129,

где

знаками

О .

*>

обозначены выходы,

соответствующие

числам Xk,

290


X ft + 1, Y k. Нетрудно видеть, что на любой совокупности клеток с двумя различными знаками не удается построить замкнутый многоугольник.

Общее количество клеток, занятых в таблице, равно 23, в то время как максимальное число линейно независимых выходов

тп — 1 = 17. Поэтому число независимых линейных соот­

ношений, которыми

связаны в совокупности

числа X ft, X fe + 1

и Yk, будет равно 6.

 

 

 

 

 

1 1

 

 

ж

1

 

 

 

ж

2

'

Ж

 

ж

о

о 3

1

о

*

о

4

ж

*

о

 

5

□!o

О

 

 

 

6

 

 

 

 

 

O :

 

 

 

 

7

 

 

 

 

 

 

1 2

3 Ч 5

6 7

8 9

10 11

Рис. 129. Пример синтеза параллельного ГПСЧ на основе семи- и одиннадцатиразрядных регист­ ров сдвига:

клетки, обозначенные символами О. *, □> соответ­ ствуют разрядам двоичных чисел Х^, Х^+1, Y ^

Произведем нумерацию выходов ГПСЧ таким образом, чтобы любая совокупность линейно зависимых выходов, соответству­ ющая в таблице вершинам замкнутого многоугольника, включала

в

себя разряды чисел Xk, X ft + 1, Y k с большими номерами /,

g,

h, т. е. чтобы получить величину min (/ +

g + h) ^ l + 1 = 9.

 

В окончательном виде таблица выходов,

определяющая топо­

логию связей для всех 16 полусумматоров и их распределение

согласно номерам

разрядов двоичных чисел X =

(хг,

х 2, . .

., xs)

и Y = (уг, у 2, .

. ., уя), представлена на рис.

130.

Для

этого

случая все шесть независимых линейных соотношений имеют вид!

хг ф

Х3 0

х6 ф х\ 0

^6 0

Уа= о,

f + g + h =

18,

Хз ®

х 2®

x s

0

x t ф

у2 ®

ys =

0

,

f + g + h = 15,

Хь е

х'з0

х'ъ

0

х'й 0

уъ 0

у 2 =

0

,

/ + £ + А =18,

х-г ®

хъ0

х'в ® х'ч 0

у2® уь = 0,

f + g +

h =

19,

хъ ® x's ® х\ ® х'з ® у2 ф г/4 = 0,

 

f + g +

h =

17,

х$ Ф Xj ф у2 ф у3= 0, f-\-g + A = 18.

19*

291


Суммируя в различных комбинациях эти уравнения, можно получить все возможные соотношения, связывающие чисела X*,

X ft + 1, Y k. При

этом оказывается, что минимальная сумма наи­

больших номеров f , g , h

у

разрядов

чисел Х А=

X, X ft + 1 =

X ' и

Y k =

Y , входящих в эти соотношения, равна min (/ + g -f- h) = 15.

Это

значит, что максимальная ошибка операции А 2В ,

вызванная

влиянием

линейной

зависимости,

не

превысит

величины |е |«=*

2“14 <

2-9.

 

 

 

 

 

 

 

 

 

 

 

 

Определим теперь интервалы сдвига s;7 между последова­

тельностями, генерируемыми

в разрядах

ГПСЧ.

Для

этого

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

ж ш

 

 

 

 

 

 

 

 

 

0

ж

 

*

 

©

 

 

 

 

 

 

0

*

 

*

 

©

 

©

 

 

 

 

 

*

 

*

 

©

 

©

 

 

и

 

 

 

 

 

ш ©

 

©

 

 

 

 

 

 

 

 

 

 

 

©

 

 

 

 

 

ш

 

 

 

 

 

 

 

 

f

2

3

9

5

6

7

8

9

10

11

 

 

 

 

Рис. 130. Таблица выходов,

определяющая струк­

 

 

 

 

туру синтезируемого ГПСЧ (обозначения те же, что

 

 

 

 

 

 

 

и на

рис.

133)

 

 

 

 

 

решим уравнение (7.24), подставляя в правую часть все возмож­

ные

значения

h =

i — /,

в

результате

получим:

 

 

 

 

h - 1 0 - 9 - 8 - 7 - 6 - 5 - 4 - 3 —2 - 1

0

1

2 3

4

5

6

У

84

101

118

8

25

42

59

76

93

111

0

17

34

51

68

85

102

Анализируя полученную таблицу и учитывая выражение (7.22), обнаружим, что минимальный интервал сдвига между последо­ вательностями будет равен min (s/7) я» 8 N = 16376 за исклю­ чением двух пар разрядов (h = —10 и 5, h = —9 и 6), для кото­ рых min (Sji) ^ N = 2047.

Рассмотренный пример показывает, что табличный метод анализа линейной зависимости между разрядами ГПСЧ позволяет сравнительно просто, не прибегая к услугам ЦВМ, синтезиро­ вать параллельный генератор на основе двух регистров сдвига. Благодаря простоте и наглядности метода удалось сократить

2 9 2


число разрядов в регистрах сдвига с 24 до 18. При этом ошибка вычислений, возникающая вследствие такого сокращения остается

пренебрежимо малой по сравнению с едоп =

2~9. Период после­

довательности генерируемых

псевдослучайных чисел

Xk и

Y k

Т = MN = 245969 >

216, а

интервалы

корреляции

т =

sjt

между числами Xk и

Xk + т, Y k к Y k + х,

X k

и Y k+T

слишком

велики, чтобы существенно влиять на точность вычислений.

 

44.

Сравнение генераторов

случайных

 

 

 

 

и

псевдослучайных

чисел

 

 

 

 

 

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

Чтобы ответить на этот вопрос, необходимо сравнить возмож­ ности генераторов случайных и псевдослучайных чисел с точки зрения быстродействия, качества генерируемых последователь­ ностей, а также сложности их реализации. В связи с этим прежде всего определим требования, предъявляемые к параметрам гене­ раторов в СтВМ. Как известно, для получения высоких характе­ ристик СтВМ (точности и быстродействия) быстродействие ГСП должно быть максимально приближено к быстродействию логи­ ческих элементов, используемых для построения операционной части машин, т. е. должно лежать в диапазоне от единиц до со­ тен МГц. Требования к качеству случайных последовательностей, т. е. к степени соответствия их статистических характеристик характеристикам гипотетической последовательности, можно опре­ делить исходя из необходимой точности вычислений.

Будем характеризовать степень приближения распределения чисел в реальной последовательности к равномерному распре­ делению отклонением вероятности появления единиц в разрядах

числа от 0,5

е = р

(1) — 0,5.

Оценим допустимую величину

этого отклонения еДоп,

исходя из условия обеспечения необходи­

мой точности

преобразования в

случайную последовательность

двоичных чисел А = (а1, а 2, . . ., at).

Так как погрешность преобразования еи = р (и) — А 2~1 зави­ сит не только от степени равновероятности 0 и 1 в разрядах слу­ чайного числа, но также и от логической схемы преобразователя «код — вероятность», для конкретности произведем эту оценку, как и раньше, для преобразователя, построенного на основе схемы сравнения (см. стр. 78). Причем ограничимся рассмотре­ нием частных случаев, поскольку для произвольных значений ех, е2, . . ., еi определение погрешности еи представляет большие трудности.

Как известно, в случае преобразования кода А = (1, 0, . . ., 0) выход преобразователя и есть функция только одной переменной

293