Математическое обеспечение информационных технологий. Непрерывные системы. Сольницев Р.И - 55 стр.

UptoLike

55
I
k
(i–1)
= (I
k
)
(i–1)
– b
k
(I
k
)
(i)
, k =1, ..., n. (2.16)
В результате вычислений формируются новые значения вектора пра-
вых частей, в которых теперь хранятся элементы Y. Решение уравнения
(2.11) проводится с помощью простой операции деления элементов век-
тора Y на элементы вектора D. В результате снова преобразуется вектор
правых частей
i
i
i
Y
z
d
=
. (2.17)
Последний этап называется “обратной подстановкой”, на котором
вычисляются значения U:
1
.
n
jj kj k
k j
UZ bU
=+
=−
(2.18)
Для хранения и обработки информации о системе по описанному
методу используется списочная структура, включающая че-
тыре списка. Первый из них назовем DIAG, размерностью n, в кото-
ром хранятся диагональные элементы матрицы G в порядке возрас-
тания номеров строк и столбцов. В списке ELEM содержатся эле-
менты матрицы G, находящиеся под диагональю. Размерность его
равна количеству ненулевых элементов под диагональю матрицы. Сле-
дующий список STROK содержит номера строк соответствующих эле-
ментов из ELEM и имеет такую же размерность. Последний список
ANS (адреса начала столбцов) размерностью n содержит номера по-
зиций в ELEM, с которых начинаются элементы данного столбца.
Для того чтобы исключить из рассмотрения нулевые элементы, из
указанных списков формируется специальная структура из двух спис-
ков. В первом из них (ADR) на каждый ненулевой элемент матрицы
B отводится по три элемента списка, где первый содержит адрес эле-
мента в массиве ELEM, второй – номер строки, а третий – адрес
следующего элемента списка. Второй список LS содержит количе-
ство ненулевых элементов в соответствующем столбце. Такая струк-
тура дает возможность легко провести вычисления адресов новых
ненулевых элементов (ННЭ) и добавление их в соответствующие
столбцы матрицы. Далее проводится обратное преобразование до-
полненных списков в исходные, но с учетом позиций ННЭ матрицы.
В ELEM в этих позициях появятся нули.