Рішення систем лінійних рівнянь
Рішення систем лінійних рівнянь
Метод оберненої матриці: для системи з n рівнянь з n невідомими. за умови що визначник матриці не дорівнює нулю, єдине рішення можна представити у вигляді. Для того щоб вирішити систему лінійних рівнянь методом оберненої матриці, необхідно виконати наступні дії:
- сформувати матрицю коефіцієнтів і вектор вільних членів заданої системи;
- вирішити систему, представивши вектор невідомих як твір матриці, оберненої до матриці системи, і вектора вільних членів.
Дана система рівнянь:
Вирішуємо на MATLAB:
A = [1 -2 1; 2 -5 -1; -7 0 1];
x = inv (A) * b% Рішення системи x = A -1 b
Рішення системи лінійних рівнянь за допомогою методу Гаусса ґрунтується на тому, що від заданої системи, переходять до системи еквівалентної, яка вирішується простіше, ніж вихідна.
Метод Гаусса складається з двох етапів:
- Перший етап - це прямий хід, в результаті якого розширена матриця системи шляхом елементарних перетворень (перестановка рівнянь системи, множення рівнянь на число, відмінне від нуля, і складання рівнянь) приводиться до ступінчастому увазі.
- На другому етапі (зворотний хід) ступінчасту матрицю перетворять так, б в перших n шпальтах вийшла одинична матриця. Останній, n +1 стовпець цієї матриці містить рішення системи лінійних рівнянь.
Порядок вирішення завдання в MATLAB наступний:
- сформувати матрицю коефіцієнтів і вектор вільних членів заданої системи;
- сформувати розширену матрицю системи, об'єднавши і;
- використовуючи функцію rref, привести розширену матрицю до ступінчастого вигляду;
- знайти рішення системи, виділивши останній стовпець матриці, отриманої в попередньому пункті;
- виконати обчислення; якщо в результаті вийшов нульовий вектор, задача вирішена вірно.
A = [1 -2 1; 2 -5 -1; -7 0 1];
C = rref ([A b]); % Приведення розширеної матриці до трикутного вигляду
x = C (1: 3,4: 4)% Виділення останнього стовпчика з матриці