Мегаобучалка Главная | О нас | Обратная связь

Метод Гаусса с выбором главного элемента





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

Анализ алгоритма решения показывает, что точность определяется значениями коэффициентов, расположенных на главной диагонали , которые называют ведущимиэлементами. На каждом шаге прямого хода предполагалось, что .Если это не так, то в качестве ведущего элемента можно использовать любой другой ненулевой коэффициент поменяв местами строки или столбцы матрицы. Однако если он окажется малым то после деления на этот элемент и последующего вычитания возникают большие погрешности округления. Чтобы избежать этого, на каждом этапе строки и столбцы переставляют так, чтобы на главной диагонали оказался наибольший по модулю элемент. Поскольку одновременный перебор строк и столбцов увеличивает количество операций и усложняет алгоритм, то обычно перебирают только строки. Найденный элемент, наибольший по модулю называют главным. В этом состоит суть метода Гаусса с выбором главного элемента.

Метод имеет малые погрешности округления при условии если матрица системы хорошо обусловлена. Матрица называется плохо обусловленной, если малые изменения её элементов приводят к существенным изменениям элементов обратной матрицы . Произведение ведущих элементов равно определителю системы:

.

Откуда следует, что если система имеет решение - , то всегда можно найти ведущие элементы отличные от нуля и, следовательно, найти решение методом Гаусса с выбором главного элемента.

Метод Гаусса с выбором главного элемента надежен, прост и наиболее выгоден для систем ЛАУ с плотно заполненной матрицей.

Численные методы решения нелинейных и трансцендентных уравнений

Метод простой итерации

Система уравнений приведена к виду: X = G(X)

 

Если она задана в виде: , то преобразуем её, добавляя X к левой и правой частям.

X = F(X)+X.

 

В этом случае решение ищется по следующей итерационной формуле:



 

Xk = G(Xk-1).

Итерации прекращаются, когда ║Xk-Xk-1 ║≤ ε1,

ε1 – заданная погрешность решения.

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

 

│∂G/∂Х│< 1.

 

Метод простой итерации является методом первого порядка, т.е.

 

║Хк – Х* ≤ с║Хк-1 – Х*║, где с – постоянная величина.

 

Метод Ньютона

Система нелинейных уравнений может быть записана виде:

F(X) = 0

 

f1(x) = 0

f2(x) = 0

………..

fn(x) = 0

 

Х = (х1; х2; …; хn)

Мы линеаризуем эту функцию в точке Хк-1 , т.е. приведем к линейному виду.

Для этого раскладываем в ряд Тейлора.

 

 

,

где Xk= Xk-1+ ∆Xk .

 

∂f1/∂x1; ∂f1/∂x2; ….∂f1/ ∂xn

∂f2/∂x1; ∂f2/∂x2;….∂f2/∂xn;

∂F/∂X = ………………………………… – матрица первых

производных.

∂fn/∂x1; ∂fn/∂x2;…..∂fn/∂xn

 

Итерационный процесс заканчивается, когда ║∆Xk ║≤ ε1

 

║ F(Xk ) ║≤ ε2,

где норма вектора отождествляется с его длиной:

 
 


А║ = √ А12 + А22 + …. +Аn2,

А= (А1, А2, …Аn)

 

Недостатком метода является то, что не всегда имеет место сходимость решения. Если имеется сходимость, то она носит квадратичный характер, т.е.

 

║Хк – Хк-12≤║Хк – Х*║,

 

где Х* - точное решение.

 

Для линейной системы уравнений решение находится за одну итерацию.

 





Читайте также:


Рекомендуемые страницы:


Читайте также:
Личность ребенка как объект и субъект в образовательной технологии: В настоящее время в России идет становление новой системы образования, ориентированного на вхождение...
Модели организации как закрытой, открытой, частично открытой системы: Закрытая система имеет жесткие фиксированные границы, ее действия относительно независимы...
Почему люди поддаются рекламе?: Только не надо искать ответы в качестве или количестве рекламы...

©2015-2020 megaobuchalka.ru Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. (1190)

Почему 1285321 студент выбрали МегаОбучалку...

Система поиска информации

Мобильная версия сайта

Удобная навигация

Нет шокирующей рекламы



(0.005 сек.)