Вычислительные методы линейной алгебры. Горбаченко В.И - 35 стр.

UptoLike

ВУЗ: 

ПГУ | Пенза

Составители: 

Рубрика: 

35
x =
-21
15
-7
10
Невязка
r =
0
0
0
0
Пример sam2_03.m демонстрирует попытку решения СЛАУ с выро-
жденной матрицей
% Метод Гаусса решения СЛАУ
% Главная программа
% Попытка решения СЛАУ с вырожденной матрицей
clear;
clc;
% Пример матрицы и правой части
disp('Матрица решаемой системы')
% Вырожденная матрица
A=[1 2 3
4 5 6
7 8 9]
disp('Вектор правой части')
b=[ 1
2
3]
[x,fl]=gauss(A,b);% Вызов метода Гаусса
if fl
disp('Решение')
x
disp('Невязка');
r=b-A*x
end
В результате получаем:
Матрица решаемой системы
A =
1 2 3
4 5 6
7 8 9
Вектор правой части
b =
1
2
3
Матрица A вырождена
Проверим, что вырожденность матрицы не обнаруживается сравнением
направляющего элемента с нулем (из-за погрешностей вычислений направ-
x =
   -21
    15
    -7
    10
Невязка
r =
     0
     0
     0
     0
     Пример sam2_03.m демонстрирует попытку решения СЛАУ с выро-
жденной матрицей
% Метод Гаусса решения СЛАУ
% Главная программа
% Попытка решения СЛАУ с вырожденной матрицей
clear;
clc;
% Пример матрицы и правой части
disp('Матрица решаемой системы')
% Вырожденная матрица
A=[1 2 3
   4 5 6
   7 8 9]
disp('Вектор правой части')
b=[ 1
    2
    3]
[x,fl]=gauss(A,b);% Вызов метода Гаусса
if fl
    disp('Решение')
    x
    disp('Невязка');
    r=b-A*x
end
В результате получаем:
Матрица решаемой системы
A =
     1     2     3
     4     5     6
     7     8     9
Вектор правой части
b =
     1
     2
     3
Матрица A вырождена
    Проверим, что вырожденность матрицы не обнаруживается сравнением
направляющего элемента с нулем (из-за погрешностей вычислений направ-

                                                                   35

Страницы