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

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

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

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

Добавлен: 15.10.2024

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

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

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

gig’s = 3r\ —4r0rx \rl — 3r|r2 + 3r0r2 —r0r3+ rxr3,

gig4= M — 5/y4+ r%— 6rxr2+ 6r0r2—4 iy 3- 4г/, —r4r4+ V i-

g\ = 4/4 — 4r/! + r%—Arxr2 2r0r2+ rf,

gigs = 6?-i — 5r0rx 4- — 9ггг2+ 4r0r2—r0r:i4- 2rxr3+ 3r\ r2r3,

gigi = 8r| — б/у'2-f — 167-^2 + Ir{)r2— 4r0r3-

87/3 + 6/'?2—4л2Л4—

2rxrx+ r0r4 -f r2r4,

 

g\ = 9r\ — 6/y4 -r ri — 18rxr2+ 6r0r2 — 2r0r3 +

6rxr3 + 9rl — 6л2г3 -f r\,

g3g4= 12/4 — ? V i ^0— 30rxr2 + 9r0r2 — 5r0r3+ 16/y3 + 18rl —

— 18л/з + 4rl — 3r/4 + r0r4 + Зл/4 —r3r4,

g\ = 16/4 — 8r0r4 + r%—48rxr2 + 12r0r2 —8r/3 + 32/y3 + 36r| —

—48л/3 + 16г| — 8г/ 4 + 2r0r4 + 12/y4 — 8r3r4 + r|.

3. Коэффициенты полинома (3.77) представляем в виде алгеб раической суммы нескольких произведений gngt, каждое из кото­ рых будет входить в эту сумму с некоторым весом wni. Такую замену производим в следующем порядке.

Выбираем любой коэффициент полинома (3.77), записав его в форме

W00r 0 Г W’o/ol + •••+ Wlirl + W12r12 + W13r13 “Г •••-f WnO'ni. (3.78)

Примем подобную же форму для записи произведений gngt

gngi = W^r\ +

w[n2l)r12+ . . . +

W%iyrnl.

 

 

 

Среди слагаемых

w’nirni

коэффициента отмечаем то из них,

которое имеет наибольшее и и г.

 

 

 

 

 

Умножаем gngt на

w„i.

Тогда искомая

величина

wni

равна

= ^nf J IffnU

если i

Ф п (gn_ xgL, если i = п),

на

такую

Умножим gng i - x,

величину wni - i {wn -u при

i = п), чтобы

 

 

 

 

 

wnl-lrni-l = ^пфЯ-Фт-! +

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

— wniw (ni)

 

 

 

 

 

 

w n i-1 =

n i- l

 

 

 

 

 

 

( r n ' - l )

 

 

 

 

 

 

 

n i-

 

 

 

 

 

Умножаем g„g/_ 2,

если

i — 1 ф n (gnxgi,

если i — 1 =

n),

на величину wni - s { wn-\i

ПРИ * — 1 =

n),

определяемую

из

условия

 

 

 

 

 

 

 

w'ni-2^ni-2 = ^ n iW ^ - V n i - i + 1 т - i "

140


откуда

 

(т)

(nt-1)

«V

u'ni-4— WniWni- 2

W ni-lW ni- 2

>{ni-2)

 

 

ni-

 

Этот процесс будем продолжать до совпадения соотношений

(3.78) и

wng\ -г w12g12+ wl3g13+ . . . + wnigni.

Аналогичные операции применяем для всех остальных коэф­ фициентов полинома (3.77).

Например, для коэффициента при р 2 полинома (3.77) запишем

Wmrl + « V 01 + Ш02Г02 + Wlirl + WW\2 + W22rl-

Используя приведенный алгоритм для определения весовых коэф­ фициентов wnl, получаем: ш22 = 4, ш12= 20~^Ц —= —36,ши =

—62 — 4 - 4 — (—36) 2

= -------------------- ---------------------

= - 6.

Таким образом, коэффициент при р 2 теперь может быть записан более компактно: 4g\ — 3Qgig3 — 6gf .

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

При I = 2

п 2 K ( z tz t') = 2 - 2 (2 h s) [ p g \ - p 2 (g21 + 2 g 1g 2) + p

2 ( 2 g 1g 2 + g D - p igl} -

t<v

 

 

 

При I =

3

 

 

п 2

K (z lzr) = 2-2 ^

s)\p3g\-p2(3g21 +

l2glg2- g l ) +

t<t'

 

 

+ р 3 (i2g,g2-4- 6gxg3 - lOgl -

2g2g3) —pi ((‘Ш з + H g\+ Sg2g3 —ga) +

 

+ P5 (10g2g3 + g%) —p e2g%].

 

При l =

4

 

 

n 2

K (ztzr) = 2-2(i+s)\p6g\-p2(6g\ + 36glg2- 4 g l) +

t< r

 

 

 

+ p3(36g±g2+ 3 6 ^ 3 + 46g| — I6g2g3+ g%) —pi ( З 6 ^ 3 + i2glgi +

A 50g| + 76g2g3

—8g2gt — 13g| + 2g3gi ) + p5 (12^ 4 + 92g2g3+

+ 20g-2g4 + 25g| —

10g-3g4 + g\) p 6 (28g2g4+ 39g2 + t O g ^ - g l) +

 

! - P 1 (2 2 g 3g ^ - g \ ) - p * 3 g \ \ .

141


При 1 = 5

п 2 К (ztzr ) = 2~2 (5; s) [plOgt - р* (I0g\ - 8С)glg2 - Щ \) + t< t’

+ Р3(80glg2+ I20glg3+ HOgi — (>0g2g3+ 5gl) —p4 (i20g1g3+ 80^^4 ~ -f 150gt + 360g2g360g2g-475g%+ 20g3g4 g\) +

+ p5 (80glg3-j- 20gig5+ 420^2^з + 200£2g4 —20g2g5+ 195gi —

— 120gr3gr4 + 1 0 g3g5+ 1 6g24 - 2gigb) -

p6 (20gl g b 260g2^4 + 40g2g5 +

 

+ 275g| + 180^3^4 — 30g3g5

36g4

+ t2g4g&gl) g-

 

~r P 1 ( 3 0 g 2g 5 + 3 2 0 g 3g 4 + З 0 ^ 3^5 -f- 36 ^ 1 — 12g'4g,5 4 - g l)

-

P8 (70g3g5+ 89g| +

12g4p5 - g l) + p 9 (38g4g5+

gl) —pi°4g§].

На рис. 62 показаны зависимости дисперсии

(3.71) и суммар­

ной

дисперсии (3.72)

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

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

nD(z) пл(ъ )

Рис. 62. Зависимости дисперсии (2) и суммарной

дисперсии (1) последовательности на выходе ФПВВ, воспроизводящего функцию sin X :

выходе ФПВВ, вычисленные с применением полученных формул для случая воспроизведения функции sin X . Из рисунка видно, что при 1 = 5 длина выходной последовательности должна быть увеличена по крайней мере в 1,5 раза в сравнении с процессом

1 4 2

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

1. Преобразование входных переменных. В этих случаях логи­ ческий преобразователь ФПВВ реализует одну из ФАЛ вида (3.65), а последовательности р (а г), р (а2), . . ., р (a t) — регу­ лярны.

2.Вычисление элементарных функций, имеющих разложение

вряд Маклорена. Примеры решения таких задач были нами уже рассмотрены.

3.Решение дифференциальных уравнений с помощью степен­ ных рядов.

4.

Решение краевых задач

математической

физики [10].

5.

Вычисление интегралов,

зависящих от

параметра.

6.Интерполирование функций.

7.Деление Х /У , основанное на использовании формулы

1 + ( 1 - Х ) + ( ! - * ) * + ( ! - Х ) » + . . . =

±

 

х


Г л а в а IV

СТОХАСТИЧЕСКИЕ ИНТЕГРАТОРЫ

Эту главу мы посвятим в основном изучению решающих бло­ ков, содержащих элементы памяти — триггеры. В СтВМ такие блоки чаще всего применяются для построения счетчиков и реги­ стров с обратной связью, выполняющих операции интегрирова­ ния, воспроизведения некоторых целых и дробных рациональных функций, итерационных и других процессов. Таким образом, класс аналитических функций, реализуемых автоматами с па­ мятью, существенно расширяется. Вместе с тем анализ таких схем в сравнении со схемами, рассмотренными в предыдущей главе (там триггеры используются только в качестве элемента задержки) оказывается достаточно сложным. Наличие памяти приводит к зависимости между состояниями автомата в соседних тактах даже в тех случаях, когда входные сигналы для разных тактов независимы.

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

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

19. Принципы стохастического интегрирования

Стохастические интеграторы составляют особый класс цифро­ вых интеграторов, на основе которых строятся ЦДА (цифровые дифференциальные анализаторы) — наиболее распространенные цифровые машины специального назначения [34, 45]. В методике подготовки задач для решения на «классических» и стохастиче­ ских интеграторах можно найти много общих черт.

Как и в детерминистском случае, при выполнении стохастиче­ ского интегрирования истинная непрерывная функция ф (Х 0)

1 4 4

обычно

аппроксимируется каким-либо известным многочленом

3 (Х 0),

над которыми производятся необходимые аналитические

операции. Построение приближающего многочлена осущест­ вляется при помощи информации об исходной функции ф (Х 0) в виде конечного множества ее квантованных значений в узлах интерполяции. В результате замены функции ф (Х 0) другой функ­ цией $ (Х 0), над которой можно выполнить операции интегриро­ вания с помощью дискретных устройств, появляется погрешность многочленной аппроксимации, сопровождающаяся погрешностью квантования.

В цифровых стохастических интеграторах обычно применяют те же формулы численного интегрирования: по прямоугольникам, трапециям или параболам, в основе которых лежат соответст­ венно ступенчатая, линейная или квадратичная аппроксимация исходной функции [33].

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

исущественные недостатки. К ним относятся большая погреш­ ность и ярко выраженная обратная зависимость между точностью

искоростью вычислений.

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

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

1. Интеграторы на основе структур комбинационного типа. Для таких структур характерно отсутствие элементов памяти. Сведения об исходных зависимостях, например, о подынтеграль­

ной функции

ф (Х 0),

хранятся

в ПЗУ стохастической машины

в виде выборочных значений $

(Х п), где п = 1, 2, . . ., N.

Если воспользоваться преобразованием независимой пере­

менной X в

виде

[93]

 

 

 

v

2га —1

то дискретную функцию <^7 (Х п) можно передать в форме совокуп­ ности стохастических последовательностей $ (Х г), $ (Х 2), . . .

. . ., ^7 (XN), заключенной в пучке из N линий.

10 В . В . Яковлев

т


Применив формулу численного интегрирования по прямо­ угольникам, составим сумму

 

 

 

5 „ - 5 0+

i ^

( X r),

 

 

 

(4.1)

 

 

 

 

 

 

Г = 1

 

 

 

 

 

где S 0 — начальное значение суммы.

 

 

 

 

 

В

свою очередь, эта сумма близка интегралу

 

 

 

 

 

 

 

 

Хп

 

 

 

 

 

 

 

 

 

 

+ l

J { X ) d X .

 

 

 

(4.2)

 

 

 

 

 

о

 

 

 

 

 

 

Логическая структура, воспроизводящая линейную зависи­

мость,

подобную

(4.1), представлена на рис.

63.

Как

известно

 

 

 

 

 

 

 

из (2.31) суммирование г пе-

*1_ £

 

£

Х/у

£

 

 

ременных

в

стохастических

*»_

 

к'

•* •• kN

 

 

 

машинах

требует

использо­

 

 

 

 

 

 

 

вания г несовпадающих коэф­

 

 

L 1

1

 

1

 

фициентов

к х+ к 2+

&3+ . . .

 

 

 

—. . . . -г

 

kikj =

0 (i=kj),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а потому формула суммиро­

 

 

si

 

 

 

Sn

вания (4.1) должна быть пре­

 

 

 

 

 

образована:

П

 

 

Рис.

63. Интегрирующая матрица

 

 

 

 

 

 

s n’ = *V t - 2 К З ( X r) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r= i

 

 

Если

 

k x k.> = . . .== к,

то

для схемы на рис. 63 окончательно

получаем

 

 

 

 

 

 

 

 

 

 

 

Sn = s 0-4-к 2

3

(Xr) ^

s 0+ k f V

(X) dX.

 

 

 

 

 

r =1

 

 

у

 

 

 

 

Если выходные последовательности, отображающие перемен­

ные

Sn, должны

быть преобразованы в цифровую форму, то в

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

Воспользовавшись формулой (1.26), определим ошибку пре­ образования е = 1,5Г-2 , где Т — целое число, представляющее время выборки.

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

Увеличивая число логических элементов и изменяя топологию связей в структуре на рис. 63, можно реализовать преобразования более сложного вида, в том числе преобразования по нескольким аргументам. Однако устройства подобного типа требуют значи­ тельного объема П ЗУ для хранения исходных данных, а также большого количества управляющих случайных последователь­ ностей для преобразования этих данных в вероятностную форму.

146