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

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

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

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

Добавлен: 15.10.2024

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

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

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

ных последовательностей, действующих на входе и в цепи обрат­

ной

связи,

 

в

схеме устанавливается

динамическое равновесие.

Как

обычно,

условия

генерирования

выходных

символов

в последовательности р (z)

определяются соотношениями

zt =

1,

если

< Сч > ,

>

< ГСЧ )

I,

zt = 0, если

( Сч > sg < ГСЧ >

где

z, — значение

выходной

переменной

в

i-м такте работы устрой­

ства;

( С ч )

i,

( ГСЧ ) i содержимое счетчика и ГСЧ в i-м такте.

 

 

 

 

 

 

 

 

 

 

 

 

tt

 

 

 

Рис. 66. Блок-схема следящего сто­

 

 

 

 

 

 

хастического интегратора:

 

 

 

 

 

 

 

1 — входная

логика;

2 — реверсивный

 

 

 

 

 

 

счетчик; 3

— схема сравнения;

4 — гене- п(х)

 

 

 

 

 

 

 

ратор

случайных чисел

—■7

 

 

 

 

 

 

Если генератор случайных чисел вырабатывает I статистически

независимых

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

с р (1)

1

в

каждой

р (0) = —

из них,

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p(z) = M (z) =

2г-1

 

 

 

 

 

 

 

 

 

2

 

р < Сч > , 2

Р< ГСЧ >,

 

 

 

 

 

 

 

 

 

 

1=

1

.

1=0

 

 

 

 

 

 

 

 

 

 

 

2 1- 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

2

Р < Гч > /-К- = 2~1М ( С ч ) .

 

 

 

Таким образом, наличие статистической связи между после­ довательностями, образующимися на различных выходах двоич­ ного счетчика, не влияет на установку величины выходной вероят­ ности схемы р (z). Поэтому для математических ожиданий содер­ жимого счетчика можно записать следующее рекуррентное соот­ ношение

M ( C 4 ) i = M ( C 4 ) i_1 ( i ---- jL )+ p ,

i =

1, 2,

3,

. . . ,

где i — номер машинного такта;

р — р (а^)

=

р

(х2)

= .

. .= ,р(а;„),

п = 21.

 

 

 

 

 

 

Заметив, что М ( С ч ) х = р

и приняв М (

Сч ) 0= 0,

выпишем

несколько последовательных значений М ( С ч ) р .

М ( Сч ) 2 = р + р (1 — 4 " ) '

М ( С ч ) 3^ р + р ( Х - ± - ) т Р ( 1- 4 ) 2,

М ( С ч ) m = p ( i + a^ -a2 + . . . + ат'1),

151


где

Выражение в квадратных скобках представляет сумму геомет­ рической прогрессии, следовательно

М < Сч > т= пр (1 ат).

Сомножитель 1 — ат можно преобразовать так:

 

\—ап = 1 : ( * - 4 г ) г -

(4.3)

При больших п, выражение (4.3) стремится к 1 — е

" , а потому

М (С ч } т = п р ( 1 — е " )

 

или

 

p ( z ) = p { i — e п ) •

(4.4)

Таким образом, если входной процесс стационарный, то про­ цесс изменения содержимого в реверсивном счетчике подчиняется

о)

 

 

S)

 

-1 Р С

Р(Т)

 

- 1рс

№)

Г Д

Р(г)

L

+ i

 

&

-1РС

 

 

о

-

 

P(Z)

~ !PC

 

 

-1РС

 

 

 

 

 

 

Рис. 67. Состав блока входной логики следящего инте­ гратора: а — при однополярном однолинейном кодирова­ нии информации; б — при ОЛС кодировании

экспоненциальному закону с постоянной 21. Рассматриваемый случай соответствует однолинейному однополярному способу кодирования информации. Состав блока входной логики для него представлен на рис. 67, а.

Четыре логических элемента в этой схеме не допускают одно­ временной подачи сигналов + 1 , —1 на входы реверсивного счет­ чика. Заметим, что форма кодирования информации на входе и выходе устройства остается одинаковой, т. е.

М { Сч >

у

21

°'

Для перехода от машинной переменной X к истинной пере­ менной Х 0 структура блока входной логики должна быть изме-

152


йена (рис. 67, б). При совпадении единичных приращений после­ довательностей р (х) и р (z) вырабатывается приращение —2, уменьшающее < Сч ) на две единицы. Так реализуется уже упо­ минавшееся нами преобразование вида 2Х ’ — 1 (стр. 31). При этом содержимое счетчика отражает действительное значение истинной переменной Х 0.

Для определения дисперсии выходной величины z следящего интегратора с входной логикой, соответствующей рис. 67, а,

применим теорему о дисперсии суммы

(2.36)

т

 

 

D(z) = Z D ( z i) + 2 £

K tf.

i=0

i<i

Корреляционный момент K tj будем искать в виде

К ( x p iZ j) = Р (xiziz>) —p{x{) р (г,-) р (zj) =

= р (Xi) р (Zi/Xi) р (Zj/XiZi) р (xt) р (Zi) р (Zj),

где р (z, |Xi) — условная вероятность равенства z-t = 1 при усло­

вии, что в том же такте xt = 1; р (zj

|xtzj)

— условная вероятность

равенства Zj = 1 при условии xt

=

1 и zt =

1.

Так как состояние входа не

влияет

на

выходную величину

в одном и том же

такте,

то р (zt |xt) =

р (z,).

Кроме того, р

(Zj\XiZi)

= р (zj),

так

как

при условии xt = 1

и z{ = 1 состояние реверсивного счетчика не изменяется, что обеспечивается включением логической схемы на рис. 67, а.

Итак, К (XiZiZj) = 0 и, следовательно,

D{z) = m p (z )[i—p(z)\.

Воспользовавшись соотношением (1.26), для точности в опреде­ лении результата ep ^ 2 _<1+i) получим

т ^ 2 ,2 Ь -2 * а+1\

(4.5)

что определяет число тактов, необходимое для достижения задан­ ной точности в установившемся режиме. Для определения общего времени интегрирования с заданной точностью гр к величине т, получаемой по формуле (4.5), нужно добавить слагаемое

т 1 4- In ■ 1

1 — Р у ст

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

шейся

вероятности р уСтр в уравнение (4.4).

В

установившемся режиме работы реверсивного счетчика

(рис. 66) возможно его переполнение. Это происходит в такой си­

туации, когда на входах + 1

и —1 счетчика одновременно появля­

ются серии из s единиц и

s нулей. При этом PC сбрасывается

и процесс интегрирования

повторяется сначала. Вероятность

153


события П, заключающегося в переполнении счетчика при подаче на его вход стохастической последовательности с р (1) = р (х), равна

Р {Щ = И15 [1 —Р (^)]s = {Р (х) (я)]2}*-

Математическое ожидание количества переполнений за т так­ тов равно

М (П) = т {р (х) — [р (х)Г2}3.

(4.6)

Количество единиц (нулей), достаточных для опрокидывания счетчика в исходное состояние, составляет

s = 2 l [i —р(х)\.

Подставляя это значение в выражение (4.6) и разрешив урав­ нение относительно I, можем получить минимально необходимое количество разрядов Zmin, достаточное для исключения нежела­ тельных переполнений

ъ М Ж

1 m in ^ 3 , 3 2 l g 1 - Р ( Х )

l g ( * ) - р 2 ( * ) ] ■

Отклонение количества переполнений счетчика N (П ) от своего центра распределения для биномиального распределения можно найти, воспользовавшись интегральной теоремой Муавра — Лап­ ласа

Р р ( Д ) - 1 / ( Л ) К в } ^ Ф ( а 0),

(4.7)

где 8 — заданная точность, равная

 

 

е = а0 ]/ D (П) =

 

= а0 V т\р (х) —р 2 (ж)]2 tx р (*)] [1 р (х) + р 2 ( z) f l 11 р <х)1,

(4.8)

Ф (а 0) — интеграл

Лапласа.

 

Подставляя (4.8)

в уравнение (4.7), окончательно получим

P { \ N ( I I ) - M ( I I ) \ ^

s ; ос0 V т . (х) —2р2 (х) + 2р3(X) —pl (х)]2‘ [1_р W]} = Ф (а0).

Пример. Определим минимально допустимое количество раз­ рядов /тш при следующих исходных данных:

доверительная вероятность 0,95;

М (П ) ^ 1; т = 220; р (х) =

= 0,5.

1

Is 2-20

^min =S= 3,32 lg

0,5

lg2-«

сокруглением в большую сторону. Точность оценки е определим из (4.8)8

8 = 1,96 | /2 20 ( - ^ - )25'° ’6 ^ 2 - 8 .

1 5 4


Таким образом, с надежностью 0,95 возможно переполнение счетчика 1 ± 2~8 раз, то есть практически не более одного раза за 220 тактов.

22. Выполнение вычислительных операций

Рассмотрим некоторые операции, наиболее эффективно и просто реализуемые при помощи следящих стохастических интег­ раторов (СлСтИ). При изображении схем таких интеграторов будем пользоваться обозначением (рис. 68), заимствованным у аналоговых вычислительных машин. Знак «-Т» с цифрой, стоящий в левой части схемы, обозначает суммирующий вход интегратора и вес приращений по этому входу. Соответственно знаком «—»

обозначен вычитающий вход СлСтИ.

 

 

 

 

Алгебраическое сложение. Сложение и

 

 

 

вычитание

нескольких

переменных Х х,

 

 

 

Х 2, . . .

осуществляется наиболее просто:

 

 

 

путем суммирования приращений р (ж,-)j ,

 

 

 

р (Xi) 2, .

. ., поступающих на различные

 

 

 

входы интегратора. Действительно, при

 

 

 

подаче на суммирующие входы интегра­

Рис. 68.

Условное обоз­

тора на

рис. 68 двух стационарных после­

начение

следящего

сто­

довательностей Х х

и Х 2 с

математиче­

хастического интегра­

скими

ожиданиями

р {х^}1 = р (х .2)1

 

тора

 

 

 

 

= . . . =

р ( s J i =

Pi,

Р (si)s

= Р (*а)г

=

зависимости (4.4),

при

= . . . =

(рхп)2 =

р 2

и,

повторяя

вывод

нулевых начальных

условиях

получим1

 

 

 

 

 

 

Р (2) =

(Рх -гРг)

1 —е

п

) .

 

 

или для

установившегося режима

 

 

 

 

 

 

 

 

 

p (z )= p 1-r pi .

 

 

 

(4.9)

Если переменные Х х и Х 2 представлены в однолинейном несим­ метричном виде, то из (4.9) следует

Z = X 1 + X 2.

Это свойство интеграторов может быть с успехом использовано для вычисления сумм, рядов и т. п. Например, на рис. 69 пред­ ставлена схема, реализующая вычисление функции

Z --- sin X

при однолинейном однополярном несимметричном кодировании информации и представлении функции с помощью трех членов степенного ряда

sin X

X

 

Х5

1 !

3 !

5! '

 

155