* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки
Работа с пакетом LinearAlgebra и алгоритмами NAG
627
8.3.3. Методы решения систем линейных уравнений средствами пакета LinearAlgebra
Конечной целью большинства матричных операций является решение систем ли нейных уравнений. Для этого пакет LinearAlgebra предлагает ряд методов и средств их реализации. Основными методами решения являются следующие: • обращение матрицы коэффициентов уравнений и решение вида X=A 1*B; • применение метода LU – декомпозиции (method='LU'); • применение метода QR – декомпозиции (method='QR'); • применение метода декомпозиции Холесски (method='Cholesky'); • метод обратной подстановки (method='subs'). Решение с применением обращения матрицы коэффициентов левой части сис темы уравнений A уже не раз рассматривалось и вполне очевидно. В связи с этим отметим особенности решения систем линейных уравнений другими методами. Любопытно отметить, что указание метода может быть сделано и без его заключе ния в одинарные кавычки.
8.3.4. Решение системы линейных уравнений методом LU декомпозиции
Зададим матрицу A левой части системы уравнений и вектор свободных членов B:
> restart; with(LinearAlgebra): UseHardwareFloats := false: > A:=<<4|.24|-.08>,<.09|3|-.15>,<.04|-.08|4>>; B:=<8,9,20>;