Вероятностное устройство для решения конечно-разностных уравнений

 

ВЕРОЯТНОСТНОЕ УСТРОЙСТВО ДЛЯ РЕШЕНИЯ КОНЕЧНО-РАЗНОСТНЫХ УРАВНЕНИЙ , содержащее блок оперативной памяти, первые вход и выход которого являются соответственно информационным входом и информационным выходом устройства, информационный регистр, первые вход и выход которого соединены соответственно с вторыми выходом и входом блока оператив ной памяти, блок постоянной памяти, вход которого является входом задания закона распределения устройства , блок усреднения вход которого соединен с вторым выходом информационного регистра, блок вычисления требуемого числа испытаний, первый вход которого является входом установки степени точностиустройства , а второй вход соединен с первым выходом блока усреднения, управляемый вероятностный

(19) (11) СОЮЗ СОВЕТСКИХ

СОЦИАЛИСТИЧЕСКИХ

РЕСПУБЛИК

3(5)) G 06 F 15/36.

ГОСУДАРСТВЕННЫЙ КОМИТЕТ СССР

flO ДЕЛАМ ИЗОБРЕТЕНИЙ И ОТКРЫТИЙ (21) 331 f887/18-24 (22) 16.07.81 (46) 30.05.84. Бюл. N 20 (72) В ° А. Песошин (71) Казанский ордена Трудового

Красного Знамени авиационный институт им. A.Н. Туполева (53) 681.3(088.8) (56) 1. Авторское свидетельство СССР

Ф 328475, кл. С 06 G 7/40, 1973.

2. Авторское свидетельство СССР

)) 369580, кл. G 06 G 7/48, 1974.

3. Авторское свидетельство СССР ,Р 362304, кл. G 06 F 15/36, 1973.

4. Авторское свидетельство СССР

М 477418, кл. G 06 F 15/32, 1976

{прототип) . (54)(57) ВЕРОЯТНОСТНОЕ УСТРОЙСТВО

ДЛЯ РЕШЕНИЯ КОНЕЧНО-РАЗНОСТНЫХ УРАВНЕНИЙ,. содержащее блок оперативной памяти, первые вход и выход которо-. го являются соответственно информационным входом и информационным выходом устройства, информационный регистр, первые вход и выход которого соединены соответственно с вторыми выходом и входом блока оперативной памяти, блок постоянной памяти, вход которого является входом задания закона распределения устройства, блок усреднения вход котоУ рого соединен с вторым выходом информационного регистра, блок вычисления требуемого числа испытаний, первый вход которого является входом установки степени точности устройства, а второй вход соединен с первым выходом блока усреднения, управляемый вероятностный (1- n ) -полюсник, первый вход которого является входом размерности исследуемой области устройства, счетный регистр, первый вход которого является входом установки начального адреса устройства, первый элемент И, первый вход которого соединен с первым вйходом счетного регистра, первый регистр, первый вход которого соединен с выходом первого элемента И, а выход соединен с третьим входом блока оперативной памяти, первый, элемент:ИЛИ, выход которого соединен с вторым входом первого элемен- та И, первый триггер, единичный вход которого соединен с первым входом первого элемента ИЛИ и является пусковым входом устройства, первый переключатель, выход которого соединен с входом сброса первого триггера, генератор тактовых импульсов, второй элемент И, первый вход которого соединен с единичным выходом первого триггера, второй соединен с выходом генератора тактовых импульсов, а выход соединен с вторым входом управляемого вероятностйого (1- n )-полюсника и четвертым входом блока оперативной памяти, второй переключатель, первый вход которого. соединен с выходом блока постоянной памяти, второй вход соединен с третьим выходом информационного регистра, а выход соединен с третьим входом управляемого вероятностного (1- )-полюсника, шифратор, выход которого соединен с вторым входом информационного регистра", первый дешифратор, вход которого соединен с четвертьв выходом информационного регистра, первый управляемый счетчик, первый выход которого соединен с, третьим входом блока вычисления требуемого числа испытаний, а второй выход соединен с входом шифратора, второй управляемый счетчик, управляющий вход которого является входом установки числа шагов устройства, а счетный вход соединен с выходом второго элемента И, третий переключатель, вход которого соединен с выходом второго управляемого счетчика, второй элемент ИЛИ, первый вход

1095190

20 которого соединен с первым выходом управляемого вероятностного (1- и )—

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

ИЛИ, причем выход третьего управляемого счетчика соединен с нулевым входом первого триггера и является выходом окончания счета устройства, управляющий вход соединен с выходом блока вычисления требуемого числа испытаний, а счетный вход соединен с вторым выходом счетного регистра и нулевым входом второго триггера, единичный вход которого соединен с первым входом шестого элемента

ИЛИ и подключен к пусковому входу устройства, единичный выход второго триггера соединен с первым входом

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

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

- шестого элемента И, а нулевой выход соединен с первыми входами четвертого и пятого элементов И, вторые входы которых соединены соответственно с выходом второго элемента И и выходом второго дешифратора, вход второго дешифратора соединен с первым выходом четвертого управляемого счетчика, второй выход которого соединен с первыми входами четвертого и пятого элементов ИЛИ, счетный вход соединен с выходом четвертого элемента И, а управляющий вход подключен к входу размерности исследуемой области устройства, второй выход управляемого вероятностного (1- n )-полюсника соединен с вторым входом шестого элемента И, выход которого соединен с первым входом третьего элемента ИЛИ, второй вход третьего элемента ИЛИ соединен с выходом пятого элемента И, а выход соединен с вторым входом первого регистра, вход второго регистра является входом установки числа узлов исследуемой области устройства, а выход соединен с первым входом седьмого элемента И, второй вход которого соединен с выходом шестого элемента ИЛИ, а выход соединен с управляющим входом первого управляемо,го счетчика, вторые входы четвертого, пятого и шестого элементов ИЛИ подключены к второму выходу первого управляемого счетчика, выход пятого элемента ИЛИ соединен с вторым входом третьего элемента И, выход четвертого элемента ИЛИ соединен с вторым входом счетного регистра.

В внутренних узлов сеточной области Г11, С21 и Г31.

Недостатком этих устройств является невысокое быстродействие.

Наиболее близким к предлагаемому. по технической сущности является вероятностное устройство для решения конечно-разностных уравнений, содержащее линию задержки, выход которой подключен через первую схему ИЛИ к второму входу первой схемы И, триггер начало-останова работы, вход которого соединен с выходом регистра стартового адреса, другой выход которого соединен с первым входом первой схемы И, выход последней подключен к первому входу регистра блужданий, второй вход которого подключен к первому выходу управляемого вероятностного (1 в n ) 1095190

Цель изобретения — повышение быстродействия.

Поставленная цель достигается тем,ба что вероятностное устройство для решения конечно-разностных уравнений, содержащее блок оперативной памяти, 1 первые вход и выход которого являются соответственно информационным 5 цолюсника, а выход — к оперативноу запоминающему устройству, первый вход управляемого вероятностного, (1-n )-полюсника соединен с вторым переключателем, второй выход подключен к первому входу второй схе мы ИЛИ, .первый вход второй схемы И подключен к триггеру начало-останова работы, второй вход - к генератору тактовых импульсов, регистр числа, выходы которого с эон ве« роятностей перехода, начальных и граничных условий, признаков узла и всего регистра числа в целом подключены соответственно .к второму пе- . реключателю, первому входу сумки- 15 рующего и усредняющего блока к блоку анализа признаков узла и к оперативному запоминающему устройству, а выходы зон начальных и граничнщс условий, признаков узла и всего ре- 2р гистра в целом соединены соответственно с выходом третьей схемы И, выходом блока изменения признака узла и оперативным запоминающим уст-ройством, суммирующий и усредняющий блок, первый и второй выходы которого подключены соответственно к первому входу третьей схемы И и первому входу процессора, второй вход которого соединен с первым выходом первого. управляемого счетчикаделителя, первый вход которого подключен к выходу процессора, а второй выход через первый переключатель к триггеру начало-останова, а также к входу регистра стартового адреса, постоянное запоминающее устройство, выход которого подключен к второму переключателю, второй управляющий счетчик-делитель, выход которого через третий переключатель подключен 4а к третьему входу второй схемы ИЛИ, второй вход которой подключен к выхо-. ду блока анализа признаков узла f43.

Однако в известном устройстве для решения указанных задач моделируются 45 случайные блуждания из стартового узла до попадания на границу сеточной области. Для. получения решения задачи по всей области с необходимой точностью d и заданной надежностью 5а вычисления (Ъ проводят требуемое число испытаний .для каждого узла области (на практике H q =10 -10 ), что существенно увеличивает объем вычислений для получения приближен- 55 ного решения задачи по всей области и снижает быстродействие устройства. входом и информационным выходом устройства, инФормационный регистр, первые вход и выход которого соединены соответственно с вторыми выходом и входом блока оперативной памяти, блок постоянной памяти, вход которого является входом задания закона распределения устройства, блок ус- реднения, вход которого соединен с вторым выходом информационного регистра, блок вычисления требуемого числа испытаний,.первый вход которо- го является входом установки степе ни точности устройства, а второй вход соединен с первым выходом блока усреднения, управляемый вероятностный

{1- n )-полюсник, первый вход которо, го является входом размерности исследуемой области устройства, счетный регистр, первый вход которого явля-" ется входом установки начального адреса устройства, первый элемент И, первый вход которого соединен с первым выходом счетного регистра, первый регистр, первый вход которого

-соединен с выходом первого элемента И, а выход соединен с третьим входом блока оперативной памяти, первый элемент ИЛИ, выход которого соединен с вторым входом первого элемента И, первый триггер, единичный вход которого соединен с первым входом первого элемента ИЛИ и является

Ьусковым входом устройства, первый переключатель, выход которого соединен с входом сброса первого триггера, генератор тактовых импульсов, второй элемент И, первый вход которого соединен с единичным выходом первого триггера, второй соединен с выходом генератора тактовых импульсов, а выход соединен с вторым входом управляемого вероятностного (1- )-полюсника и четвертым входом блока оперативной памяти, второй переключатель, первый вход которого: соединен с выходом блока постоянной памяти, второй вход соединен с третьим выходом информационного ре1 гистра, а выход соединен с третьим входом управляемого вероятностного {1- h )-полюсника, шифратор, выход которого соединен с вторым входом информационного регистра", первый дешифратор, вход которого соединен с четвертым выходом информационного регистра, первый управляемый счетчик, первый выход которого соединен с . третьим входом блока вычисления требуемого числа испытаний, а второй выход соединен с входом шифратора, второй управляемый счетчик, управляющий вход которого является входом установки числа шагов устройства, а счетный вход соединен с выходом второго элемента И, третий переключатель, вход которого соединен с выходом второго управляемого счет1095190

6 чика, второй элемент ИЛИ, первый вход которого соединен с первым выходом управляемого вероятностного (1-И )-полюсника, второй вход соединен с выходом первого дешифратора, третий вход соединен с выходом

5 третьего переключателя, а выход соединен со счетным входом первого управляемого счетчика и вторым входом первого элемента ИЛИ, элемент задержки, вход которого соединен с входом пер10 вого переключателя и вторым входом счетного регистра, а выход соединен с третьим входом первого элемента ИЛИ третий элемент И, первый вход которого соединен с вторым выходом блока 15 усреднения, а выход соединен с третьим входом информационного регистра, дополнительно содержит третий. и четвертый управляемые счетчики, второй дешифратор, второй регистр, 20 второй триггер, четвертый, пятый, шестой и седьмой элементы И, третий ,четвертый, пятый и шестой элементы

:.ИЛИ, причем выход третьего управляемого счетчика соединен с нуле- 25 вым входом первого триггера и является выходом окончания счета устройства, управляющий вход соединен с выходом блока вычисления требуемого числа испытаний, а счетный вход З0 соединен с вторым выходом счетного регистра и нулевым входом второго триггера,. единичный вход которого соединен с первым входом шестого элемента ИЛИ и подключен к пусковому входу устройства, единичный выход второго триггера соединен с первым входом шестого элемента И, а нулевой выход соединен с первыми входами четвертого и пятого элементов И, вторые входы которых соединены соответственно с выходом второго элемента И и выходом второго дешифратора> вход второго дешифратора соединен с первым выходом четвертого управляемого счетчика, второй выход которого соединен с первыми входами четвертого и пятого элементов ИЛИ, счетный вход соединен с выходом четвертого,элемента И, а управляющий вход подключен к входу размерности 50 исследуемой области устройства, второй выход управляемого вероятностного (1- n )-полюсника соединен с вторым входом шестого элемента И, выход которого соединен с первым 55 входом третьего элемента ИЛИ, второй вход третьего элемента ИЛИ соединен с выходом пятого элемента И, а выход соединен с вторым входом первого регистра, вход второго регистра явля- 60 ется входом установки числа узлов исследуемой области устройства, а выход соединен с первым входом седьмого элемента И, второй вход которого соелинен с выходом шестого эле мента ИЛИ, а выход соединен с управляющим входом первого управляемого счетчика, вторые входы четвертого, пятого и шестого элементов ИЛИ подключены к второму выходу первого управляемого счетчика, выход пятого элемента ИЛИ соединен с вторым входом третьего элемента И, выход четвертого элемента ИЛИ соединен с вторым входом счетного регистра.

На фиг. 1 приведена структурная схема устройства; на фиг. 2 — экспериментальные графики зависимости средней ошибки решения задачи ДирихI ле для уравнения Лапласа от числа итераций; на фиг.3 — структурная схема блока вычисления требуемого числа испытаний на фиг. 4 — структурная схема дисперсиометра.

Устройство имеет пусковой вход

1, первый триггер 2, первый переключатель 3, выход 4 окончания счета, генератор 5 тактовых импульсов, второй элемент 6 И, третий управляемый счетчик 7, четвертый элемент 8 И, четвертый управляемый счетчик 9, второй дешифратор 10, второй триггер 11, пятый элемент 12 И, вход

13 установки начального адреса, счетный регистр 14, первый элемент

15 И, первый регистр 16, третий элемент 17 ИЛИ, шестой элемент 18 И, вход 19 размерности исследуемой области, управляемый вероятностный (1-К ) -полюсник 20, первый элемеНт

21 ИЛИ, элемент задержки 22, вход

23 установки степени точности, вход

24 задания закона распределения, блок 25 постоянной памяти, второй переключатель 26, блок 27 усреднения, блок 28 вычисления требуемого числа испытаний, четвертый элемент

29 ИЛИ, информационный вход 30, информационный выход 31, блок 32 оперативной памяти, информационный регистр 33, третий элемент 34 И, пятый элемент 35 ИЛИ. шестой элемент 36 ИЛИ, шифратор 37, первый дешифратор 38, второй элемент 39 ИЛИ первый управляемый счетчик 40, вход 41 установки числа узлов исследуемой области, второй регистр 42, седьмой элемент 43 И. второй управляемый счечик 44, вход 45 установки числа шагов, третий переключатель 46.

Вход элемента 22 задержки соединен с входом первого переключателя

3 и с вторым входом счетного регистра 14, единичный вход первого триггера 2 соединен с пусковым входом 1 устройства и с первым входом первого элемента 21 ИЛИ, нулевой вход первого триггера 2 соединен с выходом

4.окончания счета устройства, первый вход счетного регистра 14 подключен к входу 13 установки начального ад10 95190

65 реса устройства, а первый выход соединен с первым входом первого элемента 15 И, первый вход управляемого вероятностного (1- и )-полюсника 20 соединен с входом 19 размерности исследуемой области устройства, а второй вход соединен с выходом второго элемента б И, с четвертым входом блока 32 оперативной памяти и со счетным входом второго управляемого счетчика 44, первый вход блока 28 вычисления требуемого числа испытаний соединен с входом 23 установки степени точности устройства, а третий вход соединен с первым выходом первого управляемого счетчика 40, вход блока 25 постоянной памяти соединен с входом 24 задания закона распределения, устройства, а выход » с первым входом второго переключателя

26, первый вход блока 32 оперативной памяти соединен с информационным входом 30 устройства, а первый выход с информационным вь)ходом 31 устройства, управляющий вход второго управляемого счетчика 44 соединен с входом 45 установки числа шагов устройства, а выход соединен с входом третьего переключателя 46, единич-

;ный выход первого триггера 2 соединен с первым входом второго элемен; та 6 И, второй вход которого подключен к выходу генератора 5 тактовых импульсов, выход первого переключателя 3 соединен с входом с6роса первого триггера 2, выход первОго элемента 15 И соединен с первым входом первого регистра 1б,,а второй вход подключен к выходу первого элемента 21 ИЛИ, выход первого регистра 16 подключен к третьему входу блока 32 оперативной памяти, вторые вход и выход которого соединены соответственно с первыми выходом и входом информационного регистра

33, первый выход управляемого вероятностного (1 в n )-полюсника 20 соединен с первым входом второго элемента 39 ИЛИ, третий вход подключен к выходу второго переключателя 26, выход элемента 22 задержки соединен с третьим входом первого элемента

21 ИЛИ, второй вход которого под; ключен к выходу второго элемента !

39 ИЛИ и к счетному входу первого управляемого счетчика 40, второй вход блока 28 соединен с первым выходом блока 27 усреднения, второй выход которого соединен с первым .входом третьего элемента 34 И, второй вход второго переключателя 26 соединен с третьим выходом информационного регистра 33, второй выход которого подключен к входу блока 27 усреднения, выход шифратора 37 подключен к второму входу информационного регистра 33, четвертый выход которого соединен с входом первого дешифрато ра 38, выход третьего элемента 34 И соединен с третьим входом информационного регистра 33, выход первого дешифратора 38 подключен к второму входу второго. элемента 39 ИЛИ, третий вход которого соединен с выходом третьего переключателя 46, второй выход первого управляемого счетчика

40 соединен с входом шифратора 37, выход третьего управляемого счетчи10 ка 7 подключен к нулевому входу первого триггера 2, а управляющий вход соединен с вторым выходом счетного регистра 14 и с нулевым входом второго триггера 11, единичный вход

15 второго триггера 11 соединен с единичным входом первого триггера 2 и с первым входом шестого элемента 36 ИЛИ,,единичный выход второго триггера 11 соединен с первым входом шестого элемента 18 И, второй вход которого подключен к адресному выходу управляемого вероятностного

{1- a )-полюсника 20, нулевой выход второго триггера 11 соединен с первым входом четвертого элемента 8 И и с первым входом пятого элемента 12 И, вход второго дешифратора 10 соединен с первым выходом четвертого управляемого счетчика 9, второй выход

30 которого дключен к перв " входам четвертого 29 и пятого 35 элементов ИЛИ, выход четвертого элемента

8 И соединен со счетным входом четвертого управляемого счетчика 9, управляющий вход которого подключен к первому входу управляемого вероятностного (1 в д )-полюсника 20, второй вход четвертого элемента 8 И соединен с выходом второго элемента б И, вторые входы четвертого 29, 40 пятого 35 и шестого 36 элементов ИЛИ подключены к второму выходу первого управляемого счетчика 40, выход пятого элемента 35 ИЛИ соединен с

,вторым входом третьего элемента 34 И, вход второго регистра 42 соединен с входом 41 установки числа. узлов исследуемой области,. а выход соединен с первым входом седьмого элемента 43 И, выход шестого элемента 36

ИЛИ соединей с- вторым входом седь- мого элемента 43 И, выход которого подключен к управляющему входу первого управляемого счетчика 40.

Н аз наче ние элементов устройства следующее. На пусковой вход 1 по-. дается сигнал "Начало счета", с помощью первого триггера 2 осуществляется пуск и останов устройства первый переключатель 3 используется при решении задачи для одного узла или для нескольких узлов сеточной области, с выхода 4 окончания счета поступает сигнал "Окончания счета", генератор 5 тактовых импульсов задает временную сетку работы устройст1095190

10 ва, третий управляемый счетчик 7 осуществляет подсчет числа итераций на второй ступени вычислительной процедуры, четвертый управляемый счетчик 9 обеспечивает подсчет числа тактов обхода соседних узлов и исследуемого узла сеточной области, второй дешифратор 10 формирует сигналы изменения координат узлов области в соответствии с заданным направлением их обхода; второй триггер 11 управляет .режимами первой и второй ступени вычислительной процедуры, на вход 13 установки начального адреса поступает сигнал выборки начального адреса, счетный регистр 14 15 фиксирует адрес исследуемого узла области, первый регистр 16 состоит из реверсивных счетчиков координат и обеспечивает возможность решения задачи в одно-, двух- и трехмерных областях, на вход 19 размерности исследуемой области подается сигнал установки размерности исследуемой области путем переключения выходов управляемого вероятностного (1- и )-полюсника 20 и путем изменения коэффициента деления четвертого управляемого счетчика 9, управляемый вероятностный (1 в A )-полюсник 20 предназначен для вероят- 3О ностного преобразования информации и реализации закона случайного блуждания частицы по узлам сеточной области, на вход 23 установки степени точности подается сигнал для установ ки необходимой точности решения d с,заданной надежностью вычислений р, на вход 24 задания закона распреде ления блуждания частицы по узлам сеточной области, в блоке 25 постоянной памяти хранятся значения ве-. роятностей переходов из исследуемого узла в соседние узлы, блок 27 усреднения предназначен для получения первой оценки решения задачи (среднего арифметического "веса" по траекториям45 начинающим .с узла, для которого находится решение) на первой ступени вычислительной процедуры, второй оценки .(среднего арифметического "весов" соседних узлов) и результат 5р решения среднего арифметического (пер ;вой и второй оценок решения задачи) на второй ступени вычислительной процедуры, блок 28 вычисления требуемого числа испытаний осуществляет вычисление дисперсии "веса", определение требуемого числа испытаний и числа итерации, значение которого заносится в третий управляемый счетчик 7, на информационный вход 30 поступают сигналы значений законов вероятностных переходов, начальных и граничных условий, признаков граничных или внутреннйх узлов, с информационного выхода 3 1 выдаются результаты решения, информационный регистр 33 предназначен для органиэации обмена анализируемой и результирующей информации; шифратор 37 предназначен для выработки кода, " 11" соответствующего признаку граничного

I узла и маркировки этим кодом узлов, для которых найдено решение; первый дешифратор 38 осуществляет дешифрирование содержимого последних двух разрядов информационного регистра

33, фиксирует окончание построения одной траектории, когда последняя попадает на границу; первый управляемый счетчик 40 обеспечивает построение точно и траекторий из данного узла на вход 41 установки числа узлов исследуемой области второй регистр 42 фиксирует число узлов исследуемой области; второй управляемый счетчик 44 осуществляет подсчет числа шагов траектории блуж дания; третий переключатель включается при решении параболических уравнений; на вход 45 установки числа шагов подается сигнал для устанонки числа шагов (путем изменения коэффициента деления второго управляемого счетчика).

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

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

Число испытаний на первой ступени М определяется числом узлов ис- следуемой области n., что обосновывается экспериментально на примере решения задачи Цирихле для уравнения Лапласа с числом узлов прямоугольной сеточной области р =100 и различном числе испытаний (фиг.2, где для графиков а -Mn =50, в — Nq

-100, с — Мн --- 200). Как видно из

1095190

10

15 фиг. 3 график в имеет более быструю сходимость средней ошибки решения к заданной за меньшее число итераций на второй ступени вычислительной процедуры, чем графики а и с. Число итераций на второй ступени определяется из экспериментального соотноШения =м пр ь /н.,где требуемое число испытаний laws вычисляется с .помощью процессора по первым 50-100 испытаниям с заданной точностью 8 и надежностью вычислений

При решении сеточных уравнений устройство работает следующим образом.

После занесения исходных данных в блок 32 оперативной памяти по входам 30 установки, по входу 23 необходимой точности решения с заданной вероятностью в блок 28 вычисления требуемого числа испытаний, записи по входу 41 установки числа узлов исследуемой области Н в регистр 42, по входу 13 устанавливается начальный адрес в счетный регистр 14. По начальному импульсу по входу 1 триггеры 2 и 11 переходят в единичное состояние, координаты начального узла с помощью элементов 21 ИЛИ и 15 И переносятся в регистр 16, значение из регистра 30

42 с помощью элементов 36 ИЛИ и 43 И заносится в счетчик 40. Тактовые импульсы от генератора 5 начинают поступать через открытый элемент б И на блоки устройства. Координаты на- 35 чального узла из регистра 16 поступа« ют в блок 32 оперативной памяти, где возбуждают соответствующую чис; ловую линейку. Из блока 32 в регистр

33 списываются все данные, относящиеся к данному узлу. При этом в блок 27 усреднения поступает начальное значение функции в данном узле (при решении параболических уравнений). Одновременно в управляемый вероятностный (1- н )-полюсник 20 списываются значения вероятностей перехода из узла в соседние узлы.

В соответствии с этими вероятностями начинается построение случайной траектории — на одном из выходов 50 (1, — )-полюсника появляются импульс, определяющий направление перемещения траектории. Таким образом, в следующий тактовый момент времени в регистре 16 с помощью элементов 18 И и 17 ИЛИ будет новый адрес, — адрес узла соседнего с начальным. Далее процесс повторяется. При каждом переходе траектории в новый узел дешифратор 38 анализирует последние 60 два разряда регистра 33. Данная траектория заканчивается, если дешифратором 38 зафиксирован признак граничного узла (код 11). Траектория может заканчиваться не только на 65 границе, но и в том случае, если она попадает в условное поглощающее состояние непосредственно с внутреннего узла. При решении параболических уравнений (переключатель

46 включен) траектория может закончиться после К шагов еще до попадания на границу. Во всех перечислен-

1ных случаях в счетчик 40 через эле мент 39 ИЛИ поступает импульс, уве; личивающий его содержимое на единицу. Одновременно накопившаяся к концу траектории сумма в блоке 27 поступает в блок 28 вычисления требуемого числа испытаний. Из счетного регистра 14 в регистр 16 снова переносятся координаты начального узла и подобно вышеописанному строится новая случайная траектория. После первых () = 50-100 испытаний блок 28 вычисления требуемого числа испытаний вычисляет в соответствии с полученными оценками М () Ц"(P,,заданными значениями ь, 8 требуемое число испытаний N ws и число итераций (e +1). Значение числа (т +1) поступает в счетчик 7. После переполнения счетчика 40 вычисление функции для узла заканчивается. При этом из блока 27 полученное значение функции записывается через элемент 36 И в соответствующие разряды регистра 33, одновременно последние два разряда регистра 33 маркируются единичками с помощью шифратора 37, и числовая информация из регистра 33 переписывается в блок 32 оперативной памяти по адресу начального узла. Кроме этого, в счетный регистр 14 с помощью элемента 29 ИЛИ прибавляется единица и в нем устанавливается адрес нового узла, соседнего с начальным узлом, а в счетчике 40 из регистра 42 снова заносится значение числа rl . По задержанному элементом

22 задержки импульсу переполнения счетчика 40 с помощью элементов 21 ИЛИ и 15 И новый стартовый адрес из счетного регистра 14 переписывается в регистр 16. После того, как первая оценка решения вероятностным методом будет получена для всех узлов. сеточной области, по сигналу переполнения с выхода счетного ре|гистра 14 из счетчика 7 вычитается единица, триггер 11. устанавливается в нулевое состояние и адрес начального узла сеточной области переносится в регистр 16. При этом элемент 18 И,обеспечивающий построение . случайных траекторий, запирается, а элементы 8 и 12 И, обеспечивающие выполнение итерационного алгоритма, открываются. Таким образом, устройство переходит к реализации второй ступени вычислительной процедуры.

Рассмотрим работу устройства в указанном режиме на примере решения

14

60 плоской задачи. В соответствии с состояниями счетчика 9, подсчитывающего такты обхода соседних узлов и исследуемого узла дешифратор 10 формирует сигналы йзменения координат узлов области в соответствии с за- . данным направлением их обхода, поступающие через элементы 12 И и 17 ИЛИ на реверсивные счетчики координат регистра 16. Например, узлы, соседние к исследуемому узлу Р(yIq;),мож- 10 но обходить,по следующему направлению: в первом такте в блок 27 из блока .32 через регистр 33 поступает значение оценки решения,(ф,) в узле Р ()(ь Я;),во Втором такте -(Ма) в 15 узле Й (; Ч-,, в третьем такте ,(йз в узле (>с,ф, в четвертом . такте,- (© в узле .Р,(х;,Щ ), после чего в блоке 27 вычисляется среднее арифметическое значение Я,=Щ+,М4Ф в пятом такте †(4) в узле Я (д1 ф, после чего в блоке 27 вычисляется среднее арифметическое значение „

= ф(Щ кд), которое на шестом такте по сигналу переполнения счетчика 9 через элементы 35 ИЛИ и 34 И заносится в соответствующие разряды регистра 33. Затем возбуждается адрес узла Р(„, г)в блоке 32 и числовая информация из регистра 33 переписывается в блок 32 оперативной памяти, по указанному адресу. Одновременно в счетный регистр 14 прибавляется единица и начинается аналогично описанному вьппе вычислительная процедуРа ДлЯ Узла сосеДнего с Узлом (1 г,ф, Подобные операции вычисляются для всех узлов сеточной области N число раз, после чего по сигналу переполнения счетчика 7 триггер 2 устанавливается в нулевое состояние, соот- 4О ветствующее окончанию счета, а ре зультирующая информация из блока 32 оперативной памяти выводится по выходу 31 на печать или в универсаль-. ную ЭВИ. 45

Блок 28 вычисления требуемого числа испытаний включает как показано на фиг. 3, дисперсиометр 47 и умножитель 48. Требуемое число испытаний Йщрвь вычисляется по формуле йщр 9 ффОценка дисперсии и находится с помощью дисперсиометра после построения первых 50-100 траекторий случайного блуждания иэ стартового узла, которая затем умножается 55 в умножителе 48 на коэффициент » 8 надежности и точности вычислений, задаваемый по внешнему входу 23 вероятностного устройства. Для вычисления числа итераций г по формуле

= Мщ у/п=Р+/ . // I в блоке 28 исключена операция усреднения значения D ,по д, так как и выбирается равным ближайшей меньшей степени 2 и процесс деления D» íà h состоит 65 в перенесении запятой в значении j3» на m разрядов, начиная с младшего, в сторону старшего разряда.

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

49, логический блок 50,, блок 51,управления, элементы 52, 53, 54 и 55 И элемент 56 ИЛИ, сумматоры 57 и 58.

Приемный регистр 49 служит для приема двоичного числа последовательности значений "весов" первых 50-100 ,траекторий, поступающих иэ блока 27

,усреднения вероятностного устройст-. ва в параллельном коде, выходные шины регистра 49 соединены с одним входом элемента 53 И и логическим блоком 50, состоящим из схем совпадения на два входа, выход логического блока 50 соединен с входом элемента 52 И, выходы элементов

52, 53 и 54 И соединены с соответствующими входами сумматора 57 через элемент 56 ИЛИ. Выходные шины старших разрядов сумматора 57 через элемент 55 И соединены с входами младших разрядов сумматора 58, инверсные выходы сумматора 58 через элемент 54 И и элемент 56 ИЛИ сое-. динены с входами сумматора 57. Вторые входы элементов 52, 53, 54 55 И соединены с выходами блока 51 управления представляющего собой . регистр опроса, на вход которого подается из счетчика 40 сигнал окончания приема числа в приемный регистр 49.

Число A (зйачение "веса" траекто-, рии), представленное двоичным кодом

A (а„,, а,J, подается из блока 27 усреднения параллельным кодом на приемный регистр 49. По окончании приема числа А, т.е. в тот момент, когда все разряды приемного регистра установятся в соответствующие состояния а,, на вход блока 51 управления подается сигнал окончания приема числа. Блок 51 управления вырабаты- вает управляющий сигнал разрешения

УИ, на элемент 53 И. Выходные сигналы от регистра 49 через элемент

53 И и элемент 56 ИЛИ подаются на сумматор 57 таким образом, что сигнал j -ro разряда регистра 49 подается на 2 j -й разряд сумматора 57, т.е. после окончания действия импульса разрешения от блока 51 в сумматоре 57 будет установлено число(а„о,а Я...о,а„,Оса и блок 51 управления вырабатывает второй сигнал разрешения УИ, который подается на элемент И 52. На второй вход элемента И 52 подаются сигналы с выходом логического блока 50, состоящего иэ С элементов И, гдеС„ г, число сочетаний. На каждую из этих схем совпадений подаются сигналы

1095190

20 фиг,1 с выходов двух разрядов регистра 49 сигналы с выходов блока 50 через элемент 52 И и элемент 56 ИЛИ поступают на входы сумматора 57 таким образом, что на вход К -ro разря.да сумматора 57 подаются сигналы с выхода той схемы совпадения логического блока 50, на вход которой поданы сигналы с выходов j -го и (-го разрядов регистра 49, причем

k =- +Ь 1. Укаэанная операция равносильна введению поправки Еа О для определения квадрата числа A:

Таким образом, по окончании УИ ! в сумматоре 57 получено число !! А

Ъ

Затем блок 51 управления вырабатывает сигнал УИ, подаваемый на вход элемента 54 H. Ha второй вход эле-. мента 54 И подаются сигналы с инверсных выходов сумматора 58. Под действием сигнала УИ сигналы с элемента 54 И через элемент 56 ИЛИ подаются на одноименные входы сум матора 57, при этом осуществляется операция вычитания и по окончании

УИ в сумматоре 57 находится число

А - Dz- где Зj;! — значение оценки дисперсии для предыдущего отсчета, хранимое в сумматоре 58 ° По оконЧании УИ> блок 51 управления вырабатывает сигнал УИ, подаваемый на элемент 55 И. Сигналы с выходов младших разрядов сумматора 57 пода« ются на входы младших разрядов сумматора 58 со смещением на т разря- дов, а на входы старших разрядов сумматора 58 подаются сигналы со э.накового разряда сумматора 57. Таким образом, по окончании УИ в сумматоре 58 установлено число соответствующее несмещенной оценке дисперсии. Приемный регистр может вновь принимать следующее число по окончании действия УИ2 .

Решение конечно-раэностных уравнений во всех узлах сеточной области с помощью. предлагаемого устройства, реализук!щего двухступенчатую вероятностно-детерминированную вычислительную процедуру дает более чем на порядок сокращение объема вычислительных эатрат по сравнению с выбранным прототипом. !

1095190

Составитель О. Майоров Редактор В. Ковтун Техред A.A÷ Корректор С Шекмар

Заказ 3599/31 Тираж 699 Подписное

ВНИИПИ Государственного комитета СССР по делам изобретений и открытий

113035, Москва, Ж-35, Раушская наб., д. 4/5

Филиал ППП "Патент", г. Ужгород, ул. Проектная, 4

Вероятностное устройство для решения конечно-разностных уравнений Вероятностное устройство для решения конечно-разностных уравнений Вероятностное устройство для решения конечно-разностных уравнений Вероятностное устройство для решения конечно-разностных уравнений Вероятностное устройство для решения конечно-разностных уравнений Вероятностное устройство для решения конечно-разностных уравнений Вероятностное устройство для решения конечно-разностных уравнений Вероятностное устройство для решения конечно-разностных уравнений Вероятностное устройство для решения конечно-разностных уравнений Вероятностное устройство для решения конечно-разностных уравнений Вероятностное устройство для решения конечно-разностных уравнений 

 

Похожие патенты:

Изобретение относится к компьютерному проектированию и компьютерному дизайну, и в частности к системе и способу улучшенного параметрического геометрического моделирования
Наверх