Метод вращений решения СЛАУ

МИНИСТЕРСТВО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

ГОСУДАРСТВЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО И ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ

«МАГНИТОГОРСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ им. Г. И. НОСОВА»

Кафедра вычислительной техники и прикладной математики

КУРСОВАЯ РАБОТА

по дисциплине: «Вычислительная математика»

на тему: «Метод вращений решения СЛАУ»

Исполнитель: Сысоев Н.В,, студент 2 курса, АВ-09-1.

Руководитель: Филиппов Е. Г.

Магнитогорск, 2011.

Содержание

Введение

1 Теоретический обзор

1.1 Прямые методы

1.2 Метод Гаусса

1.2.1 Описание метода

1.2.2 Сходимость метода простой итерации

1.2.3 Апостериорная оценка погрешности

1.2.4 Пример

1.3 Метод вращений линейных систем

1.3.1 Описание метода

1.3.2 Контроль точности и уточнение приближенного решения в рамках прямого метода

1.3.3 Апостериорная оценка погрешности

1.3.4 Пример

1.4 Метод релаксации

1.4.1 Пример

2 Практическая часть

2.1 Таблица идентификаторов

2.2 Листинг программы

2.3 Пример

2.4 Сравнительная таблица

Заключение

Библиографический список

Введение

Как утверждается в книге известного американского математика Валяха, 75% всех расчетных математических задач приходится на решение СЛАУ. Это не удивительно, так как математические модели тех или иных явлений или процессов либо сразу строятся как линейные алгебраические, либо сводятся к таковым посредством дискретизации и/или линеаризации. Поэтому трудно переоценить роль, которую играет выбор эффективного способа решения СЛАУ. Современная вычислительная математика располагает большим арсеналом методов, а математическое обеспечение ЭВМ – многими пакетами прикладных программ, позволяющих решать различные возникающие на практике линейные системы. Чтобы ориентироваться среди методов и программ и в нужный момент сделать оптимальный выбор нужно разбираться в основе построений методов и алгоритмов, учитывающих специфику постановок задач, знать их сильные и слабые стороны и границы применимости.

1 Теоретический обзор

1.1 Прямые методы

математический модель итерация погрешность

Все методы решения линейных алгебраических задач можно разбить на два класса: прямые и итерационные. Прямые методы – это такие методы, которые приводят к решению за конечное число арифметических операций. Если операции реализуются точно, то и решение также будет точным (в связи с чем к классу прямых методов применяют название точные методы). Итерационные методы – это методы в которых точное решение может быть получено лишь в результате бесконечного повторения единообразных действий.

Эффективность способов решения системы

или

иначе, векторно-матричных уравнений Ах=f, где f=(f1, f2, …,fn)T – вектор свободных членов и

х=( х1, х2, …,хn)T – вектор неизвестных, а – вещественная n×n-матрица коэффициентов данной системы, во многом зависит от структуры и свойств матрицы А : размера, обусловленности, симметричности, заполненности и др.

Так размерность системы (т.е число n) является главным фактором, заставляющим вычислителей отвернуться от весьма привлекательных в теоретическом плане и приемлемых на практике при небольших n формул Крамера.

1.2 Метод Гаусса

1.2.1 Описание метода

Рассмотрим один из самых распространенных методов решения СЛАУ – метод Гаусса. Этот метод (который называют также методом последовательного исключения неизвестных) известен в различных вариантах уже более 2000 лет.

Вычисления с помощью метода Гаусса состоят из двух основных этапов, называемых прямым ходом и обратным ходом. Прямой ход метода Гаусса заключается в последовательном исключении неизвестных из системы (1):

для преобразования её к эквивалентной системе с верхней треугольной матрицей. Вычисления значений неизвестных производят на этапе обратного хода.

1.2.2 Алгоритм.

1.2.3 Апостериорная оценка погрешности.

1.2.4 Пример

1.3 Метод вращений линейных систем

1.3.1 Описание метода.

Как и в методе Гаусса, цель прямого хода преобразований в этом методе – приведение системы к треугольному виду последовательным обнулением поддиагональных элементов сначала первого столбца, затем второго и т.д.

Пусть с1 и s1 – некоторые отличные от нуля числа. Умножим первое уравнение исходной системы (1) на с1, второе на s1 и сложим их; полученным уравнением заменим первое уравнение системы. Затем первое уравнение исходной системы умножаем на –s1, второе – на c1 и результатом их сложения заменим второе уравнение. Таким образом, первые два уравнения (1) заменяются уравнениями

Отсюда .

Эти числа можно интерпретировать как косинус и синус некоторого угла (отсюда название метод вращений, каждый шаг такого преобразования можно рассматривать как вращение расширенной матрицы системы в плоскости обнуляемого индекса).

В результате преобразований получим систему

где

Далее первое уравнение системы заменяется новым, полученным сложением результатов умножения первого и третьего уравнений соотведлтственно на

а третье – уравнением, полученным при сложении результатов умножения тех же уравнений соответственно на –s2 и с2.

Внимание, отключите Adblock

Вы посетили наш сайт со включенным блокировщиком рекламы!
Ссылка для скачивания станет доступной сразу после отключения Adblock!

Скачать