%Тестирование функции решения систем линейных алгебраических уравнений
SLAE
%Пример 1
% Матрица коэффициентов при неизвестных
A = [ 1 -3 6 -5 0; 4 2 1 10 2; 2 0 -9 1 6 ]
% Матрица свободных членов
B = [ 3; 5; 7 ]
% –== 1 ==–
disp(‘- – = = 1 = = – -‘);
disp(‘Стандартное решение посредствам системы MatLab X = A\B’);
X = A\B;
disp(‘X = ‘);
disp(X);
disp(‘Невязка Eps = ‘)
disp(B – A*X);
disp(‘Евклидова норма N = ‘)
disp(norm(X, ‘fro’));
% –== 2==–
disp(‘- – = = 2 = = – -‘);
disp(‘Решение MatLab c первоначальной диагонализацией по методу Гауса’);
% Приведение расширенной матрицы A|B к диагональному виду
[m, n] = size(A);
A = rref([A, B]);
B = A(:, n + 1);
A = A(:, 1 : n);
X = A\B
disp(‘Невязка Eps = ‘);
disp(B – A*X);
disp(‘Евклидова норма N = ‘);
disp(norm(X, ‘fro’));
% –== 3 ==–
disp(‘- – = = 3 = = – -‘);
disp(‘Решение системы функцией SLAE’);
% Повторный ввод параметров
A = [ 1 -3 6 -5 0; 4 2 1 10 2; 2 0 -9 1 6 ];
B = [ 3; 5; 7 ];
[X, N3] = SLAE(A, B);
disp(‘Невязка Eps = ‘);
disp(B – A*X);
Нашли опечатку? Выделите и нажмите CTRL+Enter