Файл: Яковлев, В. В. Стохастические вычислительные машины.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 + 360g2g3— 60g2g-4 — 75g%+ 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