Основные свойства обратного итерационного алгоритма решения систем линейных уравнений с положительно определенными матрицами
Иванов В.Н.

Рассматривается задача решения уравнений движения механических систем относительно ускорений при их численном интегрировании. Задача сводится к решению положительно определенных систем линейных алгебраических уравнений с медленно меняющимися коэффициентами. Предложена новая модификация метода переменной метрики Пауэлла-Бройдена, основанного на симметричной формуле ранга один пересчета матрицы, обратной к матрице системы. Получены условия локальной и глобальной сходимости алгоритма в приложении к поставленной задаче и обсуждаются его основные свойства. Доказывается, что в случае точной арифметики метод сходится за конечное число итераций, которое не превосходит ранг матрицы возмущений линейной системы. На примерах интегрирования уравнений движения конкретных механических систем показана сравнительная эффективность метода. Работа выполнена при частичной финансовой поддержке РФФИ (код проекта 11-01-96024-р_урал_а).

Ключевые слова: системы линейных алгебраических уравнений, итерационные методы, методы переменной метрики, механические системы, уравнения движения, численное интегрирование

Название статьи, аннотация и ключевые слова на английском языке

Иванов В.Н., доцент, e-mail: precol@psu.ru - Пермский государственный университет, механико-математический факультет, ул. Букирева, 15, 614990, г. Пермь