Устройство для решения интегральных уравнений фредгольма второго рода

 

Изобретение относится к вычислительной технике и предназначено для решения уравнений Фредгольма второго рода. Цель изобретения - расширение класса решаемых задач за счет увеличения допустимого диапазона изменения численного параметра уравнения. Устройство содержит генераторы 1 функций, блоки 2 сдвига, группы 3 элементов ИЛИ, сумматоры 4 невязок, кодирующие элементы 5 с приоритетом, регистры 6 старших разрядов, дешифраторы 7, сумматоры 8 искомой функции, триггеры 9, элементы И 10, элементы НЕ 11, сумматоры 12 адреса, элементы ИЛИ 13, реверсивный счетчик 14 численного параметра, первый блок 15 сравнения, сумматор 16 нормы невязок, коммутатор 17 знаков невязок, коммутатор 18 старших разрядов невязок, счетчик 19 "результативных" интераций, второй блок 20 сравнения. 2 ил.

СОЮЗ СОВЕТСКИХ

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

РЕСПУБЛИК (19) Я 0 (ll ) (51)5 Г 06 F 7/64

ОПИСАНИЕ ИЗОБРЕТЕНИЯ

К .А ВТОРСКОМУ СВИДЕТЕЛЬСТВУ

Риг. 7

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

ПО ИЗОБРЕТЕНИЯМ И ОТКРЫТИЯМ

ПРИ CHHT СССР (21) 4615790/24 (22) 05.12.88 (46) 30.12.90. Бюл, Р 48 (71) Институт кибернетики им. В.М.Глушкова (72) В .П.Боюн, Л.Г.Козлов и В.Г.Тракай (53) 681.32(088.8) (56) Авторское свидетельство СССР

У 687452, кл. G 06 Г 7/64, 1976.

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

У 1108444, кл, С 06 F 7/64, 1982.

2 (54) УСТРОЙСТВО ДЛЯ РГШЕНИЯ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ ФРЕДГОЛЬИА ВТОРОГО

РОДА (57) Изобретение относится к вычислительной технике и предназначено для решения уравнений Фредгольма второго рода. Пель изобретения — расширение класса решаемых задач за счет увеличения допустимого диапазона изменения численного параметра уравне- ° нпя. Устройство содержит генераторы функций, блоки 2 сдвига, группы 3

161 7438

Д )(„(Х,) = 4„ pl< (Х;) $ 8п (Я„((Х

ll если /К„(Х;)/С

И к p„q е ли .Е/Е (((Х )!

i 1

/И„q если (И„= (Л е., еес(ее (Eagq/E„, (K;) /I/ элементов ИЛИ, сумматоры 4 невязок, кодирующие элементы 5 с приоритетом, регистры 6 старших разрядов дешифраторы 7, сумматоры 8 искомой функции, триггеры 9, элементы И 10, элементы

НЕ 11, сумматоры 12 адреса, элементы

ИЛИ 13, реверсивный счетчик 14 чис-!

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

Цель изобретения — расширение класса решаемых задач за счет увеличения допустимого диапазона изготовления численного параметра уравнения.

На фиг. 1 представлена блок-схема предлагаемого устройства, на фиг. 2 схема. блока управления.

Предлагаемое устройство содержит 25 и генераторов 1 функций, и блоков 2 сдвига, и групп элементов ИЛИ 3, и сумматоров 4 невяэок, и кодирующих элементов 5 с приоритетом, и регистров 6 старших разрядов невязок, п дешифраторов 7, и сумматоров 8 исковой функции, и триггеров 9, и элементов И 10, и элементов НЕ 11, и сумматоров 12 адреса, и элементов

ИЛИ 13, реверсивный счетчик 14 численного параметра, первый блок 15 сравнения, сумматор 16 нормы невязок, коммутатор 17 знаков невязок,коммутатор

18 старших разрядов невяэок, счетчик

19 результативных итераций, второй 40 блок 20 сравнения, регистра 21 нормы невязок, блок 22 управления, вход

23 задания численного параметра устройства, вход 24 задания точности устройства, и входов 25 задания íà- 45 ленного параметра, первый блок 15 сравнения, сумматор 16 нормы невязок, коммутатор 7 знаков невязок, к()мму татор 18 старших разрядов невязок, счетчик 19 результативных итераций, второй блок 20 сравнения. 2 ил. чального значения искомой функции устройства, и входов 26 задания начальных значений невязок устройства, вход 27 задания начального значения нормы невязок устройства, вход 28 тактовых импульсов устройства, вход

29 задания порядка устройства, вход

30 начальной установки устройства, вход 31 запуска устройства и и выходов 32 устройства.

Блок 22 управления содержит первый 33 и второй 34 триггеры, элемент

ИЛИ 35, счетчик 36, узел 37 сравнения, первый 38, второй 39 и третий

40 элементы И, первый 41, второй 42, третий 43 и четвертый 44 элементы задержки, первый 45, второй 46, третий 47, четвертый 48 и пятый 49 выходы блока управления, первый 50 вход блока управления, второй, третий, четвертый и пятый входы блока 22 уп- равления являются соответственно входом 28 тактовых импульсов, входом

29 задания порядка, входом 30 начальной установки и входом 31 запуска устройства.

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

161 7 1

1,2,. ° .,п, 1,2,..., основание системь: счисления, принятой в устрой5 стве; значение невяэки на (1с-1) -й итерации. где (U

Ч

Е0

E,(X;) = 6 „,(X, )- ЬУ„(Х;)— и

1.Ь К(Х ; S ) gy (8 где К(Х S ) — ядро интегрального урав)

15 нения — известная функция

L — численный параметр уравнения1

h — шаг интегрирования.

На первой итерации, как и в прототипе, приращение искомой функции в 1-й точке равно по величине старшему разряду (обозначим его адрес через

1 ) невязки, полученной на предыду- 25 щей итерации в той же точке, а знак приращения определяется знаком этой невязки. Затем алгоритм вычислительного процесса изменяется. На каждой итерации вычисляется норма невяэок — 3Q сумма модулей невяэок, которая сравнивается с нормой невяэок, полученной на предыдущей итерации. В зависимости от результата сравнения величина приращения искомой функции на

35 следующей итерации либо тоже равна величине старшего разряда невязки (выполненная итерация считается результативной, она приблизила значение искомой функции к конечному резуль- 40 тату — решению интегрального уравнения), либо старший разряд невязки уменьшен на один разряд (значение ис" комой функции на текущей итерации оказалось дальше от решения, чем на предыдущей, или не приблизилось к нему — итерация нерезультативная), Ес. ли t последовательных итераций оказались результативными, то для уско, рения вычислительного процесса старшиР 5р разряд невязкп на следующей итерации увеличивается на один разряд (оптимальное значение t для разных классов задач находится в пределах

4-8 и задается разрядностью счетчи55 ка 19 результативных итераций) . .В любом случае приращение функции всегда остается величиной, кратной основанию q системы счисления, при38 6 пятой в устройстве, поэтому при вы-. числении невязки вместо операции умножения значений ядра на приращение функции выполняется операция сдвига значений ядра на соответствующее число разрядов, учитывая, что шаг интегрирования h выбирается таким образом, чтобы произведение Lh = q, где m— целое положительное число, причем

m (p (р — разрядность представления чисел в устройстве). Вместо операции умножения выражения, стоящего под знаком суммы в формуле невязки, на

Lh-,производится сдвиг его Hà m разрядов вправо — в сторону уменьшения.

Итерационный процесс заканчивается, если выпополняется условие

И

/Я (Х )/(r где z — заданная точность решение интегрального уравнения.

Предлагаемое устройство работает следующим образом.

Перед началом работы начальное приближение функции У (Х;) = 0 по входам 25 заносится в сумматоры 8 искомой функции, а соответствующее ему значение невяэок Я (Х ) — Е(X;) (r(Х;) — правая часть интегрального уравнения — известная функция) по входам 26 — в сумматоры

4 невязок. В зависимости от класса решаемых задач счетчик 19 результативных итераций принимается двух-, трех и т.д. -разрядным. Таким образом задается величин, t. На вход 23 задания численного параметра устройства подается значение m, определяемое из соотношения I h = q, на вход

24 задания точности устройства значение заданной точности r на вход 28 тактовых импульсов устройства — тактовые импульсы, на вход 29 задания порядка устройства — значение порядка и. На вход 27 задания начальной нормы невязок устройства падается значение F. = g )Я (Х,)) .

1е

Затем на вход 30 начальной установки устройства подается сигнал "Начальная установка", а на вход 31 запуска устройства — сигнал Запуск", после чего начинается работа устройства.

При выполнении очередной k-й итерации по сигналу блока 22 управления кодирующие элементы 5 с приорите1617438 том выделяют адреса старших разрядов

1; соответствующих невязок, вычисленных на (k-1)-й итерации, которые заносятся для хранения в регистры 6

5 старших разрядов невязок, а с них поступают на группу входов коммутатора

18 старших разрядов невязок, триггеры 9 выделяют знаки тех же невязок, которые поступают на группу входов коммутатора 17 знаков невязок и на информационные входы элементов И 10.

Затем бло 22 управления выдает последовательно и сигналов на входы генераторов 1 функции и управляющие входы коммутаторов 17 и 18 знаков невязок и старших разрядов невязок. При этом адреса старших разрядов 1 всех, невязок последовательно поступают на первые входы сумматоров 12 сдвига, на вторые входы которых поступает величина, хранящаяся в реверсивном счетчике 14. На первой итерации— это значение m на последующих— уточненное в соответствии с алгорит- 25 мом значение тп, = ш 11,. Полученная в сумматорах 12 сумма т„ + 1„ подает, ся на управляющие входы блоков 2 сдвига. Генераторы 1 функции выдают последовательно по словам, параллель- 30 но по разрядам значения ядер K(X, S ), которые через блоки 2 сдвига, где йро исходит сдвиг К(Х; S ) Hà m К+ 1; разрядов вправо, и группу элементов ИЛИ

3 поступают на первые информационные, входы сумматоров 4 невязок, где они складываются или вычитаются в зависимости от знака невязок f „, (Х ) которые последовательно поступают с. выхода коммутатора 17 знаков невязо ., 40 через элементы ИЛИ 13 и элементы HE „ (знак меняется на противоположный) на вторые информационные входы сумматоров 4 невязок. Таким образом, в сумматорах 4 невязок вычисляются ве-, 45 личины и

1.1+x(x,s ) gv„(s>)

1-I которые складываются с величинами не вязок „,(Х, ), хранящимися в тех же сумматорах. После этого блок 22 управления выдает сигнал на де.шифраторы 7, элементы И 10 и сумматор 16 нормы невяэок. С выхода дешифраторов

7 величины старших разрядов невяэок

F„,(Х, ) поступают в сумматоры 8 искомой функции, туда же поступают зна ки невязок, т.е. в сумматоры 8 по/ даются приращения функции gv„(X;)

В сумматорах 8 приращения фнукции складываются со значениями функции, полученными на предыдущей итерации.

Таким образом реализуется основная формула итерационного процесса: к(Х ) = к „(", ) + hV„(X, ).

Величины приращений функции с выхода дешифраторов 6 поступают такжс через группу элементов ИЛИ 3 на сумматоры 4 невяэок, знаки этих приращений с выхода блоков И 10 тоже поступают на сумматоры 4 невязок через элементы ИЛИ 13 и элементы НЕ 11.

Этим достигается получение на выходах сумматоров 4 невязок значений невязок Як(Х, ), которые поступают на сумматор 16 нормы невязки, где происходит вычисление нормы невязоа

Е „ = 5 (Я„(X;)I. Эта величина срав1 1 нивается в блоке 15 сравнения с величиной r, которая характеризует точность решения интегрального уравнения. Если Е (r, процесс счета заканчивается. С выхода блока 15 сигнал поступает в блок 22 управления, который останавливает работу устрой- ства. При Е к > г сигнал с выхода блока 15 сравнения не поступает.

Вычисленная норма Е невязок поступает также в блок 20 сравнения норм невязок, куда подается норма невязок, полученных на предыдущей итерации, хранящаяся в регистре 21 нормы невязок. По сигналу с блока 22 управления в блоке 20 сравнения вычисляется разность F. E к, . Если она окажется отрицательной, то итерация считается результативной. Сигнал с выхода блока 20 поступает в счетчик 19

I результативных итераций. При поступлении t последовательных сигналов происходит переполнение счетчика 19; сигнал с его выхода поступает на реверсивный счетчик 14, где хранится величина m . При этом вычитается единица из ш„, что влечет за собой увеличение приращения функции на один разряд на следующей итерации. Если разность Е к — F.к < оказывается поло жительной или равной нулю, то сигнал с другого выхода блока 20 поступает в счетчик 19 результативных итераций— происходит сброс его в ноль, и на реверсивный счетчик 14. При этом добавляется единица к его содержимому, 16174 что приводит к уменьшению m и приращению функции на один разряд на следующей итерации. После этого по сигналу блока 22 управления в регистр 21 нормы невязок заносится очередное значение Е.„ и »а этом заканчивается очередная итерация. После выполнения k-й итерации в сумма. торах 8 искомой функции содержатся значения функции

7 „(Х, ), которые подаются на выходы 32 устройства, а в сумматорах 4 невязок — значения соответствующих невязок f (Х, ). Последующие итерации выполняются аналогично. 15

В блоке 22 управления сигнал "Начальная -установка" по входу 30 устанавливает первый триггер 33 в состояние, открывающее первый 38 и второй

39 элемент И, подготавливая схему к 20 работе. Сигнал "Запуск" (начало первой итерации) по входу 31 через элемент ИЛИ 35 и второй элемент И 39 поступает по выходу 47 на управляющие входы кодирующих элементов 5 с 25 приоритетом и триггеров 9. Через первый элемент 41 задержки задержанный

Ф сигнал поступает на второй триггер 34, который открывает третий элемент И 40.

При этом тактовые импульсы с входа

28 устройства через третий элемент

И 40 поступают по выходу 45 на входы генераторов 1 функций и управляющие входы коммутаторов знаков 17 невязок и старших разрядов 18 невязок до тех пор, пока узел 37 сравнения, который

35 сравнивает величину порядка, поступающую с входа 29 задания порядка устройства, с количеством поступивших на счетчик 36 импульсов, не выда- 40 ет сигнал на второй триггер 34, который закрывает третий элемент И 40, прекращая подачу импульсов. Сигнал с узла 37 сравнения по выходу 46 поступает на управляющие входы дешифраторов 7, элементов И 10, сумматора 16 нормы невязки. Задержанный сигнал через открытый первый элемент

И 38 и второй элемент 42 задержки по выходу 49 поступает на управляющий вход блока 20 сравнения. Через третий элемент 43 задержки задержанный сигнал по выходу 48 поступает на управляющий вход регистра 21 нормы невязки. Через четвертый элемент 44

55 задержки сигнал поступает на элемент

ИЛИ 35, начиная следующую итерацию и т.д. После достижения заданной точности решения интегрального уравне38 1О ния блок 24 сравнения выдает сигнал, поступающий через вход 49 на первый триггер 33, который запирает первый

38 и второй 39 элементы И, останавливая работу блока управления и всего устройства.

Формула изобретения

Устройство для решения интегральных уравнений Фредгольма второго рода, содержащее п генераторов функций (n — число точек, в которых определяется искомая функция), п блоков сдвига, Il сумматоров невязок, ll триггеров, коммутатор знака невязок, блок управления, и групп элементов

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

161 соответствуюших элементов ИЛИ и входами знаков соответствующих сумматоров искомой функции, выходы коммутатора старших разрядов невязок соединены с первыми входами сумматоров адреса, выходы которых соединены через соответствующие элементы НЕ со вторыми информацио«..арми входами сумматоров невязок. выходы сумматора нормы невяэоч соединены с первыми входами первого блока сравнения, вторые входы которого соединены с входами задания точности устройства, выходы сумматоров искомой функции соединены с выходами устройства, вход таковых импульсов устройства, вход задания порядка устройства, вход начальной установки устройства и вход запуска устройства соединены с одноименными входами блока управления, выход rtepвого блока сравнения соединен с вхо дом останова блока управления, второй выход которого соединен с входами синхронизации триггеров и управляющими входами кодирующих элементов с приоритетом, третий выход блока управления соединен с вторыми входами элементов И и управляющими входами дешифраторов и сумматора нормы невязок, о т л и ч а ю щ е е с я тем, что, с целью расширения класса решаемых задач за с- ет увеличения допустимого диапазона изменения численно7438 12 го параметра уравнения, в него введены счетчик результативных итераций, второй блок сравнения, регистр нормы невяэки и реверсивный счетчик численного параметра, причем входы задания начального значения нормы невязки устройства соединены с установочными входами регистра нормы невязки, выходы которого соединены с входами первой группы второго блока сравнения, первый выход которого соединен с входом сброса счетчика результативных итераций и вычитающим входом реверсивного счетчика численного параметра, выходы сумматора нормы невязок соединены с информационными входами регистра нормы невяэки и входами второй группы второго блока сравнения, 20 второй выход которого соединен со счетным входом счетчика результативных итераций, выход переполнения которого соединен с суммирующим входом реверсивного счетчика численного

25 параметра, четвертый и пятый выходы блока управления соединены с входами синхронизации второго блока сравнения и регистра нормы невязки соответственно, входы задания начального значения, 30 численного параметра устройства соединены с установочными входами реверсивного счетчика численного параметра, выходы которого соединен с вторыми входами сумматоров адресов.

28 H 70,71