Файл: Яковлев, В. В. Стохастические вычислительные машины.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