Материалы сайта
Это интересно
Исследование систем линейных уравнений неполного ранга
%Тестирование функции решения систем линейных алгебраических уравнений 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);