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

курсовая работа

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

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

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

Отсюда .

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

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

где

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

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

где

Выполнив преобразование m-1 раз, придем к системе

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

Далее по этому же алгоритму преобразуется матрица

и т.д.

В результате m-1 этапов прямого хода система будет приведена к треугольному виду.

Нахождение неизвестных не отличается от обратного хода метода Гаусса.

Треугольная, точнее, трапециевидная структура последней системы позволяет последовательно одно за другим вычислять значения неизвестных, начиная с последнего:

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

Прямые методы часто приводят к точному решению СЛАУ при точном выполнении предусматриваемых соответствующими алгоритмами арифметических операций (без округлений).

Реальные же вычисления базируются на арифметике машинных (т.е. усеченных до определенного количества разрядов) чисел. Как отражается на результате решения системы подмена арифметики действительных чисел машинной арифметикой, зависит от самой решаемой системы, параметров применяемого компьютера и системы представления данных, способов реализации алгоритмов. В любом случае, практически вместо точного решения СЛАУ прямой метод дает приближенное решение*) (обозначим его х(0)). Подставив х(0) в выражение о:=f-Ax, называемое невязкой, по малости полученного вектора значения о(0)=f-Ax(0) можно с осторожностью судить о близости найденого решения x(0) к точному решению x. Если, напимер,

|| о(0)|| - недостаточно малая величина, то следует искать вектор-поправку p такой, что x(0)+р=х есть точное решение системы

т.е. А(х(0)+р)=f.

Последнее равносильно векторно матричному уравнению

Ар = о(0).

Таким образом, нахождение поправки сводится к решению такой же системы, как и

,

где в качестве вектора свободных членов должен быть взят вектор невязок.

Делись добром ;)