Двухпараметрические попеременно-треугольные и двуциклические методы решения сильно несимметричных систем линейных алгебраических уравнений

• Автор: Крукиер, Борис Львович
• Шифр специальности: 01.01.07
• Научная степень: Кандидатская
• Год защиты: 2006
• Место защиты: Москва
• Количество страниц: 163 с.

I. ИТЕРАЦИОННЫЕ МЕТОДЫ РЕШЕНИЯ СЛАУ

1.1. НЕКОТОРЫЕ СВЕДЕНИЯ ИЗ ТЕОРИИ МАТРИЦ И ТЕОРИИ ИТЕРАЦИОННЫХ МЕТОДОВ

I. 1.1. Типы матриц и их основные свойства

1.1.2. Сведения из теории матриц

1.1.3. Локализация спектра матриц

1.1.4. Лемма Келлога

1.1.5. Основные сведения из теории итерационных методов

1.1.6. Методы ускорения сходимости

1.2. КЛАССИЧЕСКИЕ ИТЕРАЦИОННЫЕ МЕТОДЫ (ОБЗОР)

1.2.1. Метод Якоби

1.2.2. Метод Гаусса-Зейделя

1.2.3. SOR (метод последовательной верхней релаксации)
1.2.3.1. Модифицированный SOR (modified successive overrelaxation) MSOR
1.2.3.2. Метод релаксации с ускорением (accelerated overrelaxation) - AOR
1.2.4. SSOR (симметричный метод SOR) и USSOR (Несимметричный SOR)
1.2.5. Треугольные методы
1.2.6. Попеременно-треугольные методы
1.2.7. LU - разложение
1.3. ИТЕРАЦИОННЫЕ МЕТОДЫ РЕШЕНИЯ СИЛЬНО НЕСИММЕТРИЧНЫХ СЛАУ
1.3.1. Вариационные методы
1.3.2. Кососимметричные итерационные методы (КМ)
1.3.2.1. Треугольные КМ (ТКМ)
1.3.2.2. Попеременно-треугольные КМ (ПТКМ)
1.3.2.3. Двуциклические треугольные КМ (ДТКМ)
1.3.2.4. Численное исследование на модельной задаче
1.3.3. Методы эрмитова и косоэрмитова разложения
II. ДВУХПАРАМЕТРИЧЕСКИЕ КОСОСИММЕТРИЧЕСКИЕ ТРЕУГОЛЬНЫЕ ИТЕРАЦИОННЫЕ МЕТОДЫ
II. 1. Двухпараметрический ДТКМ
II. 1.1. Условия сходимости метода ;
II. 1.2. Нахождение оптимального параметра метода
11.2. Двухпараметрический ПТКМ
11.2.1. Условия сходимости метода
11.2.2. Нахождение оптимального параметра метода
11.3. УСКОРЕНИЕ ТРЕУГОЛЬНЫХ КОСОСИММЕТРИЧНЫХ МЕТОДОВ
11.4. ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ НА МОДЕЛЬНОЙ ЗАДАЧЕ
III. ИСПОЛЬЗОВАНИЕ КОСОСИММЕТРИЧЕСКИХ ИТЕРАЦИОННЫХ МЕТОДОВ ДЛЯ ПЕРЕОБУСЛАВЛИВАНИЯ ВАРИАЦИОННЫХ МЕТОДОВ
111.1. МЕТОДЫ ПОДПРОСТРАНСТВА КРЫЛОВА
111.2. ПЕРЕОБУСЛАВЛИВАНИЕ
111.3. GMRES И ЕГО МОДИФИКАЦИИ
111.4. BlCG И ЕГО МОДИФИКАЦИИ
111.5. ПЕРЕОБУСЛАВЛИВАНИЕ GMRES И BlCG
111.6. ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ НА МОДЕЛЬНОЙ ЗАДАЧЕ
ЛИТЕРАТУРА

Теория итерационных методов является интенсивно развивающейся областью численного анализа и занимает важное место в вычислительной математике и механике.
Для решения задач математической физики широко используются методы дискретизации исходных дифференциальных или интегральных уравнений, краевых и начальных условий, которые позволяют преобразовать исходную непрерывную задачу в дискретную, т.е. перейти из бесконечномерного в конечномерное пространство, как правило, достаточно большой размерности. Далее, в этом конечномерном пространстве задачу преобразуют в систему линейных алгебраических уравнений, которую затем надо решить на ЭВМ. Такая технология решения сложных научно-технических задач, описываемых системами интегро-дифференциальных уравнений, краевых и начальных условий была разработана в начале 60-тых годов А.А.Самарским и была названа им вычислительным экспериментом. В данной работе особое внимание уделяется предпоследнему этапу технологии вычислительного эксперимента — решению системы линейных алгебраических уравнений. В соответствии с мировой статистикой 80% задач, решаемых на ЭВМ - это задачи нахождения решения системы линейных алгебраических уравнений (СЛАУ). В работе рассматриваются итерационные методы решения этой задачи, т.к. речь идет о СЛАУ содержащих сотни тысяч неизвестных и уравнений, а прямые методы их решений при таком размере СЛАУ не эффективны. Несмотря на то, что теория итерационных методов в достаточной степени разработана для достаточно большого класса матриц, остаются проблемы по созданию новых эффективных итерационных методов решения СЛАУ для матриц, обладающих достаточно специфическими свойствами. Одним из таких классов матриц являются сильно несимметричные матрицы, которые получаются, например, при центрально-разностной аппроксимации уравнения конвекции-диффузии с малым параметром при старшей производной.
В связи с этим актуальность работы обусловлена потребностью в эффективных методах решения такого класса СЛАУ.

Построение «быстрых» итерационных методов решения сильно несимметричных систем в данной работе основываются на включении в обращаемый оператор итерационного метода треугольной части лишь кососимметрической составляющей исходной матрицы.
Целью данной работы является разработка эффективных численных методов решения сильно несимметричных систем линейных алгебраических уравнений. В работе предложены двухпараметрические попеременно-треугольный и двуциклический методы для решения СЛАУ из рассматриваемого класса.
В соответствии с этими целями решен ряд задач:
• разработаны, теоретически обоснованы и численно проверены двухпараметрические попеременно-треугольный (ПТКМ) и двуциклический (ДТКМ) итерационные методы решения СЛАУ с сильно несимметричной матрицей;
• рассмотрены вопросы ускорения двухпараметрических ПТКМ и ДТКМ;
• предложено использование параметрических и безпараметрических ПТКМ и ДТКМ в качестве переобуславливателей для методов вариационного типа.
Научная новизна работы определяется полученными теоретическими результатами исследования:
• доказательством сходимости предложенных новых классов двухпараметрических попеременно - треугольных и двуциклических кососимметрических методов решения СЛАУ с сильно несимметричными матрицами;
• определением (в частных случаях) оптимальных параметров двухпараметрических ПТКМ и ДТКМ;
• исследованием возможности использования этих методов в качестве переобуславливателей для методов вариационного типа.
Разработанные итерационные методы вносят вклад в развитие численных методов решения систем линейных алгебраических уравнений специального вида. Вместе с тем, с помощью разработанных методов можно эффективно решать задачи типа «пограничного слоя» при конвективно-диффузионном переносе с преобладанием конвекции.