RU2012049C1 - Device for solution of system of linear algebraic equations - Google Patents
Device for solution of system of linear algebraic equations Download PDFInfo
- Publication number
- RU2012049C1 RU2012049C1 SU4915869A RU2012049C1 RU 2012049 C1 RU2012049 C1 RU 2012049C1 SU 4915869 A SU4915869 A SU 4915869A RU 2012049 C1 RU2012049 C1 RU 2012049C1
- Authority
- RU
- Russia
- Prior art keywords
- input
- output
- elements
- inputs
- group
- Prior art date
Links
Images
Landscapes
- Complex Calculations (AREA)
Abstract
Description
Изобретение относится к вычислительной технике и может быть использовано в высокопроизводительных специализированных вычислительных машинах и устройствах обработки сигналов для решения систем линейных алгебраических уравнений. The invention relates to computer technology and can be used in high-performance specialized computers and signal processing devices for solving systems of linear algebraic equations.
Целью изобретения является повышение точности вычисления за счет выбора ведущего элемента по столбцам матрицы. The aim of the invention is to increase the accuracy of the calculation by selecting a leading element in the columns of the matrix.
На фиг. 1 представлена структурная схема устройства; на фиг. 2 - структурная схема устройства для случая n= 3 и m= 1; на фиг. 3 - схема вычислительного модуля второго типа; на фиг. 4 - схема вычислительного модуля первого типа; на фиг 5 - схема коммутатора. In FIG. 1 shows a block diagram of a device; in FIG. 2 is a block diagram of a device for the case n = 3 and m = 1; in FIG. 3 is a diagram of a computing module of the second type; in FIG. 4 is a diagram of a computing module of the first type; Fig 5 is a diagram of a switch.
Устройство для решения систем линейных алгебраических уравнений (фиг. 1) содержит информационный вход 1, настроеные входы 2-8, синхровход 9, вычислительный модуль 10 второго типа, вычислительные модули 11i (i= ) первого типа и выход 12.A device for solving systems of linear algebraic equations (Fig. 1) contains
Вычислительный модуль 10 второго типа (фиг. 3) содержит информационный вход 13, первый 14 и второй 15 настроечные входы, синхровход 16, первый узел 17 задержки, регистр 18, узел 19 сравнения, второй узел 20 задержки, триггеры 21-24, группы 25-30 элементов И, элементы И 31-39, группы 40 и 41 элементов ИЛИ, элементы ИЛИ 42-45, элементы НЕ 46, 47, 481 и 482 и выход 49.
Вычислительный модуль 11 первого типа (фиг. 4) содержит информационный вход 50, первый 51, второй 52, третий 53, четвертый 54 и пятый 55 настроечные входы, синхровход 56, четвертый 57 и первый 58 узлы задержки, регистр 59, второй узел 60 задержки, регистр 61, третий узел 62 задержки, делитель 63, умножитель 64, вычитатель 65, коммутатор 66, триггеры 67-73, пятый 74, шестой 75, седьмой 76 узел задержки, группы 77-88 элементов И, элементы И 89-101, группы 102-106 элементов ИЛИ, элементы ИЛИ 107 и 108, элементы НЕ 109-113, элемент 114, информационный выход 115, первый 116, второй 117, третий 118 и четвертый 119 настроечные выходы.
Коммутатор 66 (фиг. 5) содержит первый 120 и второй 121 входы, узел 122 сравнения, группы элементов И 123-126, группы элементов ИЛИ 127 и 128, элемент НЕ 129, первый 130, второй 131 и третий 132 выходы. The switch 66 (Fig. 5) contains the first 120 and second 121 inputs,
Вычислительный модуль 10 второго типа обладает возможностью реализации следующих функций:
Aj+1= , где bj= aj-1, если ( αj-1, βj-1 )= (1, 1)V
( αj-1, βj-1, γj-1 )= (1, 0, 1)V
( αj-1, βj-1, γj-n-2 )= (0, 0, 1),
γj= ,
где aj - значение на информационном входе на j-м такте вычислительного модуля;
αjиβj - значения соответственно на первом и втором настроечных входах на j-м такте вычислительного модуля;
Aj+1 - значение на выходе на (j+1)-м такте вычислительного модуля.The
A j + 1 = , where b j = a j-1 , if (α j-1 , β j-1 ) = (1, 1) V
(α j-1 , β j-1 , γ j-1 ) = (1, 0, 1) V
(α j-1 , β j-1 , γ jn-2 ) = (0, 0, 1),
γ j = ,
where a j is the value at the information input on the j-th clock of the computing module;
α j and β j are the values, respectively, at the first and second tuning inputs on the j-th clock of the computing module;
A j + 1 is the value at the output at the (j + 1) -th clock cycle of the computing module.
Вычислительный модуль первого типа обладаети возможностью реализации функций
Uj+2n+2= Uj ;
Vj+2n+2= Vj;
Wj+2n+2= ωj;
( где Uj, Vj, ωj и εj - значение соответственно на первом, втором, третьем и четвертом настроечных входах на j-м такте вычислительного модуля, Uj, Vj, Wj и Ej - значения соответственно на первом, втором, третьем и четвертом настроечных выходах на j-м такте вычислительного модуля,
где fj= dj-1-bj-1 cj-1, если α5 j-1= (0, 1, 0, 0)V
α6 j-1 = (0, 0, 1, 0) V α7 j-1 = (0, 0, 0, 0, 1),
bj-1= dj-2/dj-n-4, если α2 j-2= (1, 1, 0),
, и 0
,
Pj = , где
τ
τ
τ
τ
τ
τ
τ
Aj+1 - значение на информационном выходе на (j+1)-м такте вычислительного модуля.The computing module of the first type has the ability to implement functions
U j + 2n + 2 = U j ;
V j + 2n + 2 = V j ;
W j + 2n + 2 = ω j ;
( where U j , V j , ω j and ε j are the values at the first, second, third and fourth tuning inputs, respectively, at the jth clock of the computing module, U j , V j , W j and E j are the values at the first, respectively second, third and fourth tuning outputs on the j-th clock of the computing module,
where f j = d j-1 -b j-1 c j-1 if α 5 j-1 = (0, 1, 0, 0) V
α 6 j-1 = (0, 0, 1, 0) V α 7 j-1 = (0, 0, 0, 0, 0, 1),
b j-1 = d j-2 / d jn-4 if α 2 j-2 = (1, 1, 0),
, and 0
,
P j = , Where
A j + 1 is the value at the information output at the (j + 1) -th cycle of the computing module.
Рассмотрим работу вычислительных модулей второго и первого типа. Consider the work of computing modules of the second and first type.
Вычислительный модуль 10 второго типа работает в четырех режимах, которые задаются комбинацией управляющих сигналов αиβ , подаваемых соответственно на входы 14 и1 5. The
В первом режиме ( α, β )= (1, 1), на выходе элемента И 35 формируется единичный сигнал, который открывает элемент И 32. При этом через элементы И 32 и ИЛИ 42 тактовый импульс подается на синхровод регистра 18, что обеспечивает запись в регистр 18 элемента а, подаваемого на вход 13. Узел 17 задержки содержит n+1 последовательно соединенных регистров. Информация из регистров 17i(i= ) записывается в регистры 17i+1. Триггер 21 устанавливается в единичное состояние, единичный сигнал с выхода которого через элемент ИЛИ 45 открывает элементы И группы 29. Содержимое регистра 17n через элементы И группы 29 и ИЛИ группы 40 выдается на выход 49 вычислительного модуля 10. В триггерах 20j (j= ) происходит циклическая перезапись информации (элемент И 39 открыт).In the first mode (α, β) = (1, 1), a single signal is generated at the output of the And 35 element, which opens the And 32 element. At the same time, the clock pulse is fed through the And 32 and OR 42 elements to the
Во втором режиме ( α, β )= (1, 0), на выходе элемента И 36 формируетя единичный сигнал, который открывает элементы И 31, 34 и подается на входы элементов И группы 26. Входной элемент а, подаваемый на вход 13, сравнивается по модулю с содержимым регистра 18 в узле 19 сравнения. Если < Рег. 18 > ≥ a, то на выходе узла 19 сравнения формируется нулевой сигнал, который открывает элементы И группы 25.In the second mode (α, β) = (1, 0), a single signal is generated at the output of the And 36 element, which opens the And 31, 34 elements and is fed to the inputs of the And elements of
Входной элемент а через элементы И группы 25 и или группы 41 записывается в регистр 171. Содержимое регистров 17i (i= ) записывается в регистры 17i+1. Кроме того, нулевой сигнал с выхода узла 19 сравнения через элементы И 34 и ИЛИ 43 записывается в триггер 201. Содержимое триггеров 20i (i= ) записывается в триггеры 20i+1.The input element a through the AND elements of
Если < Рег. 18 > < a, то на выходе узла 19 сравнения формируется единичный сигнал, который подается на входы элементов И 31 и 34 и элементов И группы 26. Содержимое регистра 18 через элементы И группы 26 и ИЛИ группы 41 записывается в регистр 171. Входное данное а записывается в регистр 18 (тактовый импульс через элементы И 31 и ИЛИ 42 подается на синхровход регистра 18). Единичный сигнал с выхода узла 19 сравнения через элементы И 34 и ИЛИ 43 записывается в триггер 201. Триггер 22 устанавливается в единичное состояние, что обеспечивает выдачу содержимого регистра 17n+1 на выход 49.If <Reg. 18> < a , then at the output of the
В третьем режиме ( α, β )= (0, 0). На выходе элемента И 38 формируется единичный сигнал, который подается на вход элемента И 33 и на входы элементов И групп 27 и 30. Если триггер 20n+1 находится в единичном состоянии, то элемент И 33 открывается и разрешается запись входного данного а в регистр 18. Кроме того, открываются элементы И группы 27 и через элементы И группы 27 и ИЛИ группы 41 содержимое регистра 18 записывается в регистр 171. Информация из регистров 17i(i= ) записывается в регистры 17i+1. Если триггер 20n+1 находится в нулевом состоянии, то на выходе элемента НЕ 482 формируется единичный сигнал, который открывает элементы И группы 30. Входное данное записывается в регистр 171 через элементы И групп 30 и ИЛИ группы 41. На выход 49 вычислительного модуля 10 выдается содержимое регистра 17n+1.In the third mode, (α, β) = (0, 0). At the output of the And 38 element, a single signal is generated, which is fed to the input of the And 33 element and to the inputs of the And elements of
В четвертом режиме ( α, β )= (0, 1). При этом на выходе элемента И 37 формируется единичный сигнал, который через элемент ИЛИ 44 открывает элемент И 39, что обеспечивает циклическую запись информации в триггерах 20i (i= ). Информация из регистров 17i (i= ) записывается в регистры 17i+1. Триггер 23 устанавливается в единичное состояние и открываются элементы И группы 28. Содержимое регистра 18 выдается через элементы И группы 28 и ИЛИ группы 40 на выход 49 вычислительного модуля.In the fourth mode (α, β) = (0, 1). At the same time, at the output of the And 37 element, a single signal is formed, which opens the And 39 element through the
Вычислительный модуль 11 первого типа работает в семи режимах, которые задаются комбинацией внешних управляющих сигналов u, v, ω, ε и внутреннего управляющего сигнала φ , который формируется только в седьмом режиме работы в результате сравнения чисел по модулю на выходе 132 узла сравнения в коммутаторе 66 (фиг. 5). Соответствующие значения управляющих сигналов uj, vj, ωj, εj и φj на j-м такте определяют управляющие сигналы αi j+1= τi j= (uj, vj, ωj, εj, φj), (i= 1,7), где i - номер режима работы, которые подаются на соответствующие логические элементы вычислительного модуля 11. Причем управляющий сигнал τi j обеспечивает запись данных в регистр или триггер на j-м такте, а управляющий сигнал αi j+1 - на (j+1)-м такте. На вход 55 постоянно подается единичный сигнал.The
Управляющие сигналы u, v, ω и ε подаются соответственно на входы 51, 52, 53 и 54 и выдаются соответственно на выходы 116, 117, 118 и 119. Причем сигналы u, v и ωзадерживаются на 2n+2 тактов соответственно триггерами 74i, 75i, 76i (i= ) соответствующих узлов 74, 75 и 76 задержки, а сигнл εзадерживается на n тактов триггерами 57i (i= ) узла 57 задержки.The control signals u, v, ω, and ε are supplied to the
Первый режим работы вычислительного модуля 11 задается управляющими сигналами τ1 j= (uj, vj, ωj)= (1, 1, 1) (значение управляющего сигнала ε на режим работы не влияет). На выходе элемента И 95 формируется единичный сигнал τ1, на выходе элемента НЕ 111 - нулевой сигнал , которые подаются соответственно на входы элементов И групп 87, 88 и элемента И 94. При этом входной элемент aj, подаваемый на вход 50, записывается в регистр 601 узла 60 задержки на n+1 тактов через элементы И группы 87 и ИЛИ группы 106. Информация из регистров 60i (i= ) записывается в регистры 60i+1. Аналогисчно информация записывается и в триггеры 57i (i= ) и 61i (i= ).The first operating mode of the
Второй режим работы задается управляющими сигналами α
Третий режим работы определяется управляющими сигналами τ3 j= α3 j+1= (uj, vj, ωj)= (1,0,0). Уп-равляющие сигналы τ3 jи α3 j+1 формируются соответственно на выходах элемента И 97 и триггера 68. Сигнал τ3 jподается на входы элементов И групп 78, 79 и на вход элемента И 91. Сигнал α3 j+1 подается на вход элементов И группы 86, а сигнал α3 j+1 - на вход элементов И группы 80. При этом в регистр 581 записывается содержимое регистра 61 через элементы И 78 и ИЛИ группы 102. Содержимое регистра 59 через элементы И группы 86 и ИЛИ группы 104 выдается на (j+1)-м такте на выход 115. Единичный сигнал, подаваемый на вход 55, через элементы И 91 и ИЛИ 107 устанавливает триггер 571 в единичное состояние. Информация из триггера 57i (i= ) записывается в триггер 57i+1. Аналогично информация циклически переписывается в регистрах 60i(i= ) и триггерах 62i (i= ).The third mode of operation is determined by the control signals τ 3 j = α 3 j + 1 = (u j , v j , ω j ) = (1,0,0). The control signals τ 3 j and α 3 j + 1 are formed respectively at the outputs of the element And 97 and
В четвертом режиме работы подаются управляющие сигналы τ4 j= α4 j+1= (uj, εj)= (0,1) (сигналы ujиωjна режим работы не влияют). На выходах элемента И 98 и триггера 69 формируются соответственно управляющие сигналы τ4 jиα4 j+1. Управляющий сигнал τ4 j подается на входы элементов И группы 77, управляющий сигнал α4 j+1 - на входы элементов И группы 81. При этом входной элемент aj записывается в регистр 581. На выходе умножителя 64 формируется произведение <Рег. 601>˙<Рег. 61>, которое подается на вход вычитателя 65. На второй вход последнего подается содержимое регистра 581 и на его выходе формируется значение <Рег. 581>-<Рег. 601>˙<<Рег. 61>, которое на (j+1)-м такте через элементы И группы 81 и ИЛИ группы 103 записывается в регистр 582. Информация из регистров 58i (i= ) записывается в регистры 58i+1. Содержимое регистра 58n через элементы И группы 80 и ИЛИ группы 104 ( α3= 1) выдается на выход 115. Аналогично информация записывается в регистры 60i (i= ) и триггеры 57i (i= ) и 62i (i= ).In the fourth operating mode, control signals τ 4 j = α 4 j + 1 = (u j , ε j ) = (0,1) are given (the signals u j and ω j do not affect the operating mode). The outputs of the element And 98 and the trigger 69 are formed, respectively, the control signals τ 4 j and α 4 j + 1 . The control signal τ 4 j is supplied to the inputs of the AND elements of group 77, the control signal α 4 j + 1 is fed to the inputs of the elements AND of
Пятый режим работы задается управляющими сигналами τ5 j= α5 j+!= (uj, vj, ωj, εj)= (0,1,0,0, ). На выходах элемента И 99 и триггера 70 формируется соответственно сигналы τ5 j и α5 j+1. Управляющий сигнал τ5 j подается на входы элементов И группы 77 При этом в регистр 581записывается входной элемент aj, на выходе вычитателя 65 формируется значение <Рег. 581>-<Рег. 601>˙<Рег. 61>, которое на (j+1)-м такте по управляющему сигналу α5 j+1 через элементы И группы 83 и ИЛИ группы 106 записывается в регистр 59 (элемент И 90 открыт и тактовый импульс подается на синхровход регистра 59). На выход 115 выдается содержимое регистра 58n ( = 1). Аналогично осуществляется запись информации в регистры 58i (i= ), 60i (i= ) и триггеры 57i (i= ) и 62i(i= ).The fifth operating mode is set by control signals τ 5 j = α 5 j +! = (u j , v j , ω j , ε j ) = (0,1,0,0,). At the outputs of the element And 99 and
В шестом режиме работы подаются управляющие сигналы τ6 j= α6 j+1= (uj, vj, ωj, εj)= (0, 0, 1, 0). На выходе элемента И 100 формируется сигнал τ6 j, на выходе триггера 71 - сигнал α6 j+1. Управляющий сигнал α5 j+1 подается на входы элементов И группы 84 и на вход элемента И 90. Сигнал подается на входы элементов И группы 77. Сигнал τ6 j подается на входы элементов И групп 82, 85 и на входы элементов И 90 и 93. При этом в регистр 581 записывается входной элемент aj, на выходе вычитателя 65 формируется значение A= <Рег. 581>-<Рег. 601>˙<Рег. 61> которое подается на первый вход коммутатора 66, на второй вход которого подается содержимое B регистра 59. Если B ≥ A , то на первом выходе коммутатора 66 выдается значение A, на втором выходе - значение B и на третьем выходе - сигнал φj= 0. На (j+1)-м такте по сигналу α6 j+1 значение A записывается через элементы И группы 82 и ИЛИ группы 103 в регистр 582, значение B записывается через элементы И группы 85 и группы ИЛИ 105 в регистр 59 и сигнал φj = 0 - через элементы И 93 и ИЛИ 108 в триггер 621. Если B < A , то на первом выходе коммутатора 66 выдается значение B, на втором выходе - значение A и на третьем выходе - сигнал φj= 1. На (j+1)-м такте по сигналу α6 j+1значение B записывается в регистр 582, значение A - в регистр 59 и сигнал φ = 1 - в триггер 621. Содержимое 58n выдается на выход 115 ( = 1). Аналогично информация записывается в регистры 58i (i= ), 60i(i= ) и триггеры 57i (i= ), 62i (i= ).In the sixth mode of operation, control signals τ 6 j = α 6 j + 1 = (u j , v j , ω j , ε j ) = (0, 0, 1, 0) are given. The signal τ 6 j is formed at the output of the And 100 element, and the signal α 6 j + 1 is generated at the output of the trigger 71. The control signal α 5 j + 1 is fed to the inputs of the AND elements of group 84 and to the input of the AND 90 element. The signal is fed to the inputs of the elements AND of group 77. The signal τ 6 j is fed to the inputs of the AND elements of
Седьмой режим работы определяется управляющими сигналами τ7 j= α7 j+1= (uj, vj, ωj, εj, φj)= = (0,0,0,0,1). Сигнал τ7 j формируется на выходе элемента И 101, сигнал α7 j+1 - на выходе триггера 72, сигнал α7 j+1- на выходе триггера 73. Сигнал τ7 j подается на входы элементов И группы 77, сигнал α7 j+1 - на входы элементов И групп 83 и 84 и на вход элемента И 90, сигнал α7 j+1- на входы элементов И группы 81. При этом в регистр 581 записывается элемент ai, на выходе вычитателя 65 формируется значение A= <Рег. 581>-<Рег. 601>˙<Рег. 61> содержимое регистров 58i (i= ) записывается в регистры 58i+1, содержимое регистра 58n выдается на выход 115 ( = 1). Если триггер 61n+1 в нулевом состоянии ( φ 0), то = 1 и значение A записывается в регистр 582 на (j+1)-м такте через элементы И группы 81 и ИЛИ группы 103. Если триггер 61n+1 в единичном состоянии, то α7 j+1= 1, значение A записывается в регистры 59 на (j+1)-м такте через элементы И группы 84 и ИЛИ группы 105, содержимое B регистра 59 через элементы И группы 83 и ИЛИ группы 103 записывается в регистр 582.The seventh mode of operation is determined by the control signals τ 7 j = α 7 j + 1 = (u j , v j , ω j , ε j , φ j ) = (0,0,0,0,1). The signal τ 7 j is generated at the output of the And 101 element, the signal α 7 j + 1 is at the output of the
Коммутатор 66 (фиг. 5) работает следующим образом. На входы 119 и 120 подаются соответственно значения A и B. Если B ≥ A , то на выходе узла 122 сравнения формируется единичный сигнал, на выходе элемента НЕ 129 - нулевой сигнал ( φ= 0), элементы И групп 125 и 126 закрыты, а групп 123 и 124 открыты, значение A через элементы И группы 123 и ИЛИ группы 127 подается на выход 130, а значение B через элементы И группы 124 и ИЛИ группы 128 - на выход 131. Если B < A , то на выходе узла 122 сравнения формируется нулевой сигнал, элементы И групп 123 и 124 закрыты, а элементы И групп 126 и ИЛИ группы 128 подается на выход 131, а значение B через элементы И группы 125 и ИЛИ группы 127 на выход 130 и φ = 1.The switch 66 (Fig. 5) operates as follows.
В основу работы устройства для решения систем линейных алгебраических уравнений вида A1 ˙X= A2, где A1= (aij), i, j= , A2= (aij), i= , j= , X= (xij), i= , j= , положен метод Гаусса-Жордана с частичным выбором ведущего элемента по столбцам матрицы в виде следующих рекуррентных соотношений:
b1j= a1j, j= , ,
i= :
если b ≥ a , то a(o) ij= aij, j= , иначе a(o) ij= b1j,b1j= aij, j= ; a(o) 1j= b1j, j=
k = :
a
a
b
i =
c
иначе a
a
a
a
a
k = :
a
a
a
a
a
xij= a
Рассмотрим работу устройства.The basis of the device for solving systems of linear algebraic equations of the form A 1 ˙ X = A 2 , where A 1 = (a ij ), i, j = , A 2 = (a ij ), i = , j = , X = (x ij ), i = , j = , put the Gauss-Jordan method with a partial selection of the leading element in the columns of the matrix in the form of the following recurrence relations:
b 1j = a 1j , j = ,,
i = :
if b ≥ a , then a (o) ij = a ij , j = otherwise a (o) ij = b 1j, b 1j = a ij , j = ; a (o) 1j = b 1j , j =
k = :
a
a
b
i =
c
otherwise a
a
a
a
a
k = :
a
a
a
a
a
x ij = a
Consider the operation of the device.
На вход 1 устройства подаются элементы aij в моменты времени taij= i+(n+1)j-n-1, i= , j= . На входы 5 и 6 подаются соответственно управляющие сигналы αиβ , принимающие значение 0 или 1. Первый режим работы вычислительного модуля 10 задается комбинацией сигналов τ1= (α, β)= (1,1), второй режим - τ2(1, 0), третий режим - τ3= (0, 0) и четвертый режим - τ4= (0, 1). Организация ввода управляющих сигналов τi t (i= ), где t - номер такта, следующая:
(вторая группа повторяется n+m-2 раз, т. е. всего подается одна первая группа и n+m-1 вторых групп).Elements a ij are fed to input 1 of the device at time ta ij = i + (n + 1) jn-1, i = , j = . The control signals α and β, taking the
(the second group is repeated n + m-2 times, i.e., only one first group and n + m-1 second groups are served).
На входы 2, 3 и 4 подаются соответственно управляющие сигналы u, v и ω , принимающие значение 0 или 1. На входы 2, 3 и 4 подаются следующие комбинации сигналов u, v и ω : δ1= (u, v, ω )= (1, 1, 1), δ2= (1, 1, 0), δ3= (0, 1, 0), δ4= (0, 0, 1), δ5 = (1, 0, 0) и δ6 = (0, 0, 0). Организация ввода управляющих сигналов δi τ (i= ) следующая:
. . . (третья группа повторяется n+m-3 раз, т. е. всего подается одна первая группа, одна вторая группа и n+m-2 третьих групп). На вход 7 постоянно подается нулевой сигнал, на вход 8 - единичный сигнал.Control inputs u, v, and ω, which take the
. . . (the third group is repeated n + m-3 times, i.e., one first group, one second group and n + m-2 third groups are served in total). At
Элементы xij формируются на выходе 12 устройства в моменты времени
txij= 2n(n+1)+i+j(n+1)-1.Elements xij are formed at the output of the
t xij = 2n (n + 1) + i + j (n + 1) -1.
Последний элемент xnm формируется на (2n2+m(n+1)+3n-1)-м такте. Период ввода элементов aij очередной задачи решения систем линейных алгебраических уравнений равен (n+1)(n+m) тактом.The last element x nm is formed at the (2n 2 + m (n + 1) + 3n-1) -th beat. The input period of elements a ij of the next task of solving systems of linear algebraic equations is equal to (n + 1) (n + m) tact.
Рассмотрим работу устройства для случая n= 3 и m= 1 (фиг. 2). Пусть требуется решисть систему линейных алгебраических уравнений вида
Входные и выходные данные, состояния регистров и триггеров вычислительных модулей 10, 111, 112 и 113 приведены в табл. 1-4.Consider the operation of the device for the case n = 3 and m = 1 (Fig. 2). Let a system of linear algebraic equations of the form
The input and output data, the status of the registers and triggers of the
Вчислительный модуль 10 обеспечивает выбор ведущего элемента по столбцам матрицы, т. е. перемещает строки так, чтобы строка, содержащая максимальный элемент первого столбца, была первой. С первого по третий такты выбирается ведущий элемент (амакс= а31= 2), который записывается в регистр 18, и на третьем такте в триггер 201 записывается значение γ= 1. На последующих тактах в соответствии с приведенным алгоритмом работы вычислительного модуля 10 происходит перемещение местами элементов первой и третьей строк исходной матрицы. На выходе 49 формируются элементы исходной матрицы в моменты времени
(i= , j= ).The calculating
(i = , j = )
В вычислительном модуле 111 с пятого по седьмой такты в регистры 60i (i= 1, n+1) записываются элементы a(o) i1 (i= ), которые на последующих тактах циклически переписываются. На девятом такте формируется элемент a(1) 42= a(1) 12= a(o) 12/a(o) 11= 3/2, который записывается на девятом такте в регистр 61. На десятом такте формируется элемент a(1) 22= a(o) 22-a(o) 21 ˙a(1) 12= 2-0 ˙3/2= 2, который записывается на одиннадцатом такте в регистр 59. На одиннадцатом такте формируется элемент a(1) 32= a(o) 32-a(o) 31˙ a(1) 12= 2-1 ˙3/2= 1/2, который записывается в регистр 582 на двенадцатом такте, и в триггер записывается сигнал φ = 0, так как a<a. На тринадцатом такте в регистр 582 записывается элемент а(1) 42= 3/2. Элементы а(1) 22, а(1) 32 и а(1) 42 выдаются на выход 115 соответственно на двенадцатом, тринадцатом и четырнадцатом тактах. На тринадцатом такте формируется элемент а(1) 43= а(1) 13= а(о) 13/a(o) 11= 0/2= 0, который на четырнадцатом такте записывается в регистр 61. На четырнадцатом такте формируется элемент a(1) 23= а(о) 23-а(о) 21 ˙а(1) 13= 1-0˙ 0= 1, который на пятнадцатом такте записывается в регистр 59. На пятнадцатом такте формируется элемент а(1) 33= а(о) 33-а(о) 31 ˙а(о) 13= 0-1 ˙0= 0, который на шестнадцатом такте записывается в регистр 582, так как триггер 624 находится в нулевом состоянии ( φ= 0). На семнадцатом такте элемент а(1) 43= 0 записывается в регистр 581. Элементы а(1) 23, а(1) 33 и а(1) 43выдаются на выход 115 соответственно на шестнадцатом, семнадцатом и восемнадцатом тактах. Аналогично формируются элемет a(1) 44= a(1) 14= a(o) 14/a(o) 11= 4/2= 2 (записывается в регистр 61 на восемнадцатом такте), элемент a(1) 24= a(o) 24-a(o) 21 ˙a(1) 14= 5-0 ˙2= 5 (записывается в регистр 582 на двадцатом такте), элемент a(1) 34= a(o) 34-a(o) 31 ˙a(1) 14= 3-1 ˙2= 1 (записывается в регистр 582 на двадцать первом такте). Элемент а(1) 44= 2 записывается в регистр 581 на двадцатом такте. Элементы а(1) 24, а(1) 34 и а(1) 44 выдаются на выход 115 соответственно на двадцатом, двадцать первом и двадцать втором тактах.The
В вычислительном модуле 112 с тринадцатого по пятнадцатый такты записываются в регистры 60i (i= ) элементы a(1) i2 (i= ), которые на последующих тактах циклически переписываются. На семнадцатом такте формируется элемент а(2) 53= а(2) 23= а(1) 23/а(1) 22= 1/2, который на восемнадцатом такте записывается в регистр 61. На восемнадцатом такте формируется элемент а(2) 33= а(1) 33-а(1) 32˙ а(2) 23= -1/4, который записывается в регистр 59 на девятнадцатом такте. Элемент а(2) 43= а(1) 43-а(1) 42 ˙а(2) 23= 0-3/2 ˙1/2= -3/4 формируется на девятнадцатом такте, который записывается в регистр 582 на двадцатом такте. Также на двадцатом такте в регистр 581 записывается элемент а(2) 53= 1/2. Элементы а(2) 33, а(2) 43 и а(2) 53 задаются на выход 115 соответственно на двадцатом, двадцать первом и двадцать втором тактах. На двадцать первом такте формируется элемент а(2) 54= а(2) 24-а(1) 24˙а(1) 22= 5/2, который записывается на двадцать втором такте в регистр 61. На двадцать втором такте формируется элемент а(2) 34= а(1) 34-а(1) 32˙ а(2) 24= 1-1/2˙ 5/2= -1/4 , который записывается в регистр 59 на двадцать третьем такте. Элемент а(2) 44= а(1) 44-а(1) 42˙ а(2) 24= -7/4 формируется на двадцать третьем такте и на двадцать четвертом такте записывается в регистр 582. Элемент а(2) 54 записывается на двадцать четвертом такте в регистр 581. Элементы а(2) 34, а(2) 44 и а(2) 54 выдаются на выход 115 соответственно на двадцать четвертом, двадцать пятом и двадцать шестом тактах. Состояние триггера 624 (значение φ ) на режим работы вычислительного модуля 112 не влияет, так как комбинация управляющих сигналов (u, v, ω, ε )= (0, 0, 0, 0) на входах 51, 52, 53 и 54 соответственно отсутствует.In the
В вычислительном модуле 113 двадцать первого по двадцать третий такты в регистры 60i (i= ) записываются элементы а(2) 33, а(2) 43 и а(2) 53, которые на последующих тактах циклически переписываются. На двадцать пятом такте формируется элемент а(3) 64= а(3) 34= а(2) 34/а(2) 33= 1, который записывается в регистр 61 на двадцать шестом такте. Элемент а(3) 44= а(2) 44-а(2) 43˙ а(3) 34= -7/4+3/4 ˙1= -1 формируется на двадцать шестом такте и записывается в регистр 582 на двадцать сельмом такте. На двадцать седьмом такте формируется элемент a(3) 54= a(2) 54-a(2) 53˙ a(2) 34= 5/2-1/2˙ 1= 2, который записывается на двадцать восьмом такте в регистр 582. Состояние триггера 624 (значение φ ) не влияет на режим работы вычислительного модуля 113 (отсутствует комбинация входных сигналов (u, v, ω, ε )= = ( 0, 0, 0, 0). Элементы х1= а(3) 44= -1, х2= а(3) 54= 2 и х3= а(3) 64= 1 выдаются на выход 115 соответственно на двадцать восьмом, двадцать девятом и традцатом тактах.In the
Таким образом, предлагаемое устройство по сравнению с прототипом обладает более высокой точностью вычисления за счет возможности выбора ведущего элемента. В прототипе не выбирается ведущий элемент и возможны случаи давления на нуль или числа, близкие по величине к нулю. Это приводит к неработоспособности прототипа. Thus, the proposed device in comparison with the prototype has a higher calculation accuracy due to the ability to select a leading element. In the prototype, the leading element is not selected and there may be cases of pressure on zero or numbers close in magnitude to zero. This leads to the inoperability of the prototype.
Claims (2)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SU4915869 RU2012049C1 (en) | 1991-03-31 | 1991-03-31 | Device for solution of system of linear algebraic equations |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SU4915869 RU2012049C1 (en) | 1991-03-31 | 1991-03-31 | Device for solution of system of linear algebraic equations |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2012049C1 true RU2012049C1 (en) | 1994-04-30 |
Family
ID=21563096
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
SU4915869 RU2012049C1 (en) | 1991-03-31 | 1991-03-31 | Device for solution of system of linear algebraic equations |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2012049C1 (en) |
-
1991
- 1991-03-31 RU SU4915869 patent/RU2012049C1/en active
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kung et al. | A highly concurrent algorithm and pipeleined architecture for solving Toeplitz systems | |
RU2012049C1 (en) | Device for solution of system of linear algebraic equations | |
Zeman et al. | A high-speed microprogrammable digital signal processor employing distributed arithmetic | |
RU2037199C1 (en) | Device for inverting n x n matrices | |
SU1444760A1 (en) | Device for squaring a sequential series of numbers | |
RU2037197C1 (en) | Device for solving systems of linear algebraic equations | |
RU2022339C1 (en) | Multiplier | |
RU2012050C1 (en) | Device for computation of eigenvalues of (n x n) matrix | |
SU1667055A1 (en) | Device for modulo m multiplication | |
RU2324219C1 (en) | Logic computing device | |
SU1119025A1 (en) | Device for implementing fast fourier transform of sequence with zero elements | |
SU1508233A1 (en) | Device for computing fast fourier transform | |
SU1621033A1 (en) | Device for with check for multiplying numbers | |
SU1716536A1 (en) | Device for multiplying matrices | |
RU2051412C1 (en) | Device for solving systems of linear algebraic equations | |
RU1777155C (en) | Device for lv-resolution of matrixes | |
SU1170462A1 (en) | Fast fourier transformer | |
SU1633529A1 (en) | Device for majority sampling of asynchronous signals | |
SU256367A1 (en) | ACCUMULATING PARAMEL ACTION | |
RU2049354C1 (en) | Device implementing method of least squares | |
RU1778762C (en) | Matrix inversion device | |
RU1807499C (en) | Matrix multiplier | |
RU2012048C1 (en) | Device for computation of two-dimensional discrete fourier transform | |
SU1705836A1 (en) | Matrix multiplier | |
SU1262519A1 (en) | Device for logical processing of information |