Лабораторная № 6. Численное решение систем линейных алгебраических уравнений

Сентябрь 12, 2022

Главная
Математика
Лабораторная № 6. Численное решение систем линейных алгебраических уравнений

Содержание

2. СЛАУ n-го порядка В векторной форме:
3. Здесь индексы α, β, ..., ω пробегают все возможные n! перестановок номеров 1, 2, ..., n;
4. Метод исключения Гаусса
5. Идея: последовательное исключение неизвестных, приводящее исходную систему к треугольному виду, в котором все коэффициенты ниже главной
6. Получим новую систему, в которой первое уравнение осталось неизменным, а остальные больше не содержат член с
7. При помощи второго уравнения исключаются коэффициенты при х2 из третьего, четвертого и последующих уравнений. Продолжая этот
8. На каждом шаге новые значения коэффициентов определяются через значения на предыдущем шаге согласно где m –
9. Обратный ход метода Гаусса состоит в последовательном вычислении xn, xn–1, ..., x1 по алгоритму Во время
10. Если элемент на главной диагонали мал, то эта строка умножается на большие числа, что приводит к
11. Метод прогонки
12. Модификация метода Гаусса, применяемая к системам с матрицей трехдиагонального типа (часто встречаются при решении краевых задач
13. В развернутом виде: b1x1 + c1x2 = d1; a2x1 + b2x2 + c2x3 = d2; a3x2
14. На этапе прямого хода каждое неизвестное xi выражается через xi+1 xi = Ai ⋅ xi+1 +
15. Расчет начинается с вычисления коэффициентов
16. На этапе обратного хода из последнего уравнения системы при i = n–1 Далее посредством прогоночных коэффициентов
17. Задание 1. Решить СЛАУ методом исключения Гаусса. Привести результат решения, а также вид матрицы системы и
18. Кому какая система? Для Гаусса
20. Скачать презентацию

Слайд 2

СЛАУ n-го порядка
В векторной форме:

Слайд 3

Здесь индексы α, β, ..., ω пробегают все возможные n! перестановок

номеров 1, 2, ..., n; k – число инверсий в данной перестановке.

Определитель (детерминант) матрицы А n-го порядка:

Слайд 4

Метод исключения Гаусса

Слайд 5

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

котором все коэффициенты ниже главной диагонали равны нулю. Процесс приведения матрицы коэффициентов к треугольному виду называется прямым ходом метода Гаусса.
Возьмем первое уравнение системы и вычтем его из второго, предварительно умножив на такое число, чтобы уничтожился коэффициент при х1. Затем таким же образом вычтем первое уравнение из третьего, четвертого и т.д.

Слайд 6

Получим новую систему, в которой первое уравнение осталось неизменным, а остальные

больше не содержат член с х1, т.е. исключили все коэффициенты первого столбца матрицы, лежащие ниже главной диагонали:

Слайд 7

При помощи второго уравнения исключаются коэффициенты при х2 из третьего, четвертого

и последующих уравнений. Продолжая этот процесс, можно исключить из матрицы все коэффициенты, лежащие ниже главной диагонали. Исключение неизвестных повторяется до тех пор, пока в левой части последнего n-го уравнения не останется одно неизвестное хn

Слайд 8

На каждом шаге новые значения коэффициентов определяются через значения на предыдущем

шаге согласно
где m – номер уравнения, из которого исключается xk; k – номер неизвестного, которое исключается из оставшихся (n–k) уравнений (номер столбца, из которого исключаются элементы); l – номер столбца исходной матрицы.

Слайд 9

Обратный ход метода Гаусса состоит в последовательном вычислении xn, xn–1, ...,

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

Слайд 10

Если элемент на главной диагонали мал, то эта строка умножается на

большие числа, что приводит к значительным ошибкам округления при вычитаниях. Чтобы избежать этого, каждый цикл всегда начинают с перестановки строк. Среди элементов столбца находят главный, т. е. наибольший по модулю в k-м столбце, и перестановкой строк переводят его на главную диагональ, после чего делают исключения.

Слайд 11

Метод прогонки

Слайд 12

Модификация метода Гаусса, применяемая к системам с матрицей трехдиагонального типа (часто

встречаются при решении краевых задач для дифференциальных уравнений второго порядка). Каноническая форма:
aixi–1 + bixi + cixi+1 = di;
i= 1…n; a1 = cn = 0

Слайд 13

В развернутом виде:
b1x1 + c1x2 = d1;
a2x1 + b2x2 + c2x3 =

d2;
a3x2 + b3x3 + c3x4 = d3;
. . .
an–1xn–2 + bn–1xn–1 + cn–1xn = dn–1;
anxn–1 + bnxn = dn .
При этом, как правило, все коэффициенты bi ≠ 0.

Слайд 14

На этапе прямого хода каждое неизвестное xi выражается через xi+1
xi =

Ai ⋅ xi+1 + Bi для i = 1,2, ..., n–1,
посредством прогоночных коэффициентов Ai и Bi, которые имеют следующий вид
где еi = аi ⋅ Аi–1 + bi (i=2,3, ..., n–1).

Слайд 15

Расчет начинается с вычисления коэффициентов

Слайд 16

На этапе обратного хода из последнего уравнения системы при i =

n–1
Далее посредством прогоночных коэффициентов последовательно вычисляем xn–1, xn–2, ..., x1.
При условии | bi | ≥ | ai | + | ci |, деление на «0» исключается, и система имеет единственное решение.

Слайд 17