Введем понятие системы единичных векторов n-мерного векторного пространства. Если положить, что j-я компонента вектора Ej равна единице, а остальные равны нулю, то при j =1, 2,..., n получим систему
(*)
которая называется системой единичных векторов n-мерного векторного пространства.
Система единичных векторов (*) образует одинизбазисов n-мерного векторного пространства; компоненты же любого n-мерного вектора можно считать координатами этого вектора в единичном базисе.
Докажем это утверждение. Покажем, что система единичных векторов линейно независима. Составим соотношение
k1E1+k2E2+...+knEn =0.
Подставляя вместо E1, E2,...,En их выражения, умноженные соответственно на kj, и складывая, получим
(k1,k2,…,kn) = (0; 0;...; 0),
откуда k1=0, k2=0,…,kn =0. Таким образом, система векторов (*) линейно независима.
С другой стороны, всякий n-мерный вектор А = (a1, a2,...,an) представляет собой линейную комбинацию векторов этой системы:
A= a1E1+a2E2+...+anEn = (a1, a2,...,an)
Пример, вектор А = (3;-2; 4; -5) можно записать как линейную комбинацию единичных векторов E1= (1; 0; 0; 0); E2 = (0; 1;0;0); Ез = (0; 0; 1; 0); Е4 =(0; 0; 0; 1) с коэффициентами, которые равны компонентам вектора А:
А = ЗЕ1 - 2Е2 + 4Ез - 5Е4 = (3; - 2; 4; - 5).
Вычислим ранг матрицы, используя теоремы 1.3. Назовем элементарными следующие преобразования матрицы:
- перемену мест двух строк или столбцов;
- умножение строки или столбца на произвольное, отличное от нуля число;
- прибавление (вычитание) к одной строке (столбцу) кратного другой строки (столбца).
Эти преобразования не меняют ранг матрицы. Матрицы, получаемые в процессе преобразования, эквивалентны исходной. В результате элементарных преобразований векторы-строки (столбцы) матрицы можно преобразовать в единичные или в нулевые векторы. Число единичных векторов в преобразованной матрице определяет ее ранг.
2.5. Решение системы линейных уравнений методом Жордана – Гаусса
Рассмотрим систему из т линейных уравнений с n неизвестными:
Обозначим через А = (aij) матрицу системы, через X =(xj) - матрицу-столбец, состоящую из неизвестных, через В = ( bi ) - матрицу-столбец, состоящую из свободных членов; тогда системуможно записать в виде матричного уравнения
AX = B.
При решении системы линейных уравнений воспользуемся численным методом, который позволяет с помощью элементарных преобразований за конечное число шагов найти решение (если оно существует) и при необходимости получить обратную матрицу. Этот метод называется методом полного исключения неизвестных или методом Жордана — Гаусса.
Неизвестное хк называется разрешенным, если какое-нибудь уравнение содержит хк с коэффициентом 1, а в остальных уравнениях хк не содержится.
Суть метода состоит в том, что, рассмотрев первое уравнение, а в нем неизвестное, например, х1 с коэффициентом а11, отличным от нуля (он в дальнейшем называется разрешающим элементом ), и, разделив первое уравнение на этот коэффициент, с помощью первого уравнения исключают это неизвестное из всех уравнений, кроме первого. Для этого преобразуют систему, умножая обе части 1-го уравнения на а21 и вычитая из соответствующих частей 2-го уравнения. Затем обе части 1-го уравнения умножают на а31 и вычитают из соответствующих частей 3-го уравнения и т.д. Если коэффициент при х 1 в первом уравнении равен нулю,то можно выбрать любое уравнение, в котором коэффициент при х1 отличен от нуля, и с помощью этого уравнения исключить х1 из всех остальных уравнений.
Выбрав во втором уравнении неизвестное с коэффициентом, отличным от нуля, и разделив на него второе уравнение, с помощью второго уравнения исключают другое неизвестное из всех уравнений, кроме второго, и т. д., т. е. с помощью одного уравнения производят полное исключение одного неизвестного. Процесс продолжают до тех пор, пока не будут использованы все уравнения. При этом возможны следующие случаи.
1. В процессе исключений левая часть i- го уравнения системы обращается в нуль, а правая часть равна некоторому числу, отличному от нуля, т. е. имеет место равенство 0=bi¹0. Это означает, что система несовместна, т.е. не имеет решений.
2. Левая и правая части i-го уравнения обращаютсяв нуль. Это означает, что i-е уравнение является линейной комбинацией остальных, ему удовлетворяет любое найденное решение системы, поэтомуономожет быть отброшено.
3. После того, как все уравнения использованыдля исключения неизвестных, либо будет получено решение,либо доказано, что система несовместна.
Рассмотрим применение метода Жордана — Гауссана примере системы линейных уравнений, в которой количество уравнений меньше количества неизвестных.
Пример 1. Решить системууравнений
или найти коэффициенты разложения вектора В по векторам A1, А2, А3, А4:
A1x1 +А2x2 + А3x3 + А4x4 =B.
Решение. Запишем коэффициенты при неизвестных и свободные члены в таблицу, к элементам которой и применим метод полных исключений.
В первом уравнении коэффициент при x1 отличен от нуля: a11=5¹0; выбираем его за разрешающий элемент. Разделим первое уравнение на 5. Исключим x1 из всех уравнений, кроме первого.
А1 | А2 | А3 | А4 | А0 |
-2 | -2 | |||
2/5 -14/5 14/5 | 3/5 19/5 1/5 | 3/5 4/5 1/5 | 1/5 18/5 -13/5 | |
8/7 -19/14 | 5/7 -2/7 | 5/7 -9/7 | ||
3/7 3/56 1/4 | 3/7 -53/56 1/4 |
Выделены разрешающие элементы уравнений, с помощью которых производим полные исключения неизвестных. Окончательно система принимает вид
Откуда находим общее решение: х1 = 3/7-3/7x4,
x2 =- 53/56-З/5б х4,
x3 = 1/4-1/4x4.
Общее решение системы – это равносильная разрешенная система, в которой разрешенные неизвестные выражены через свободные.
Векторы A1, А2, А3 линейно независимые и составляют базис в трехмерном пространстве; соответствующие им неизвестные x1, х2, х3 называются базисными, а х4, соответствующее вектору A4 - свободным (ему можно придавать произвольные значения).
Минор, состоящий из коэффициентов при базисных неизвестных в линейно независимых уравнениях, называется базисным минором.
Любые m переменных системы линейных уравнений с n переменными (m<n) называются базисными, если определитель (базисный минор) матрицы коэффициентов при них отличен от нуля, остальные n-m переменных называются свободными (небазисными).
Базисными могут быть разные группы из n переменных. Максимально возможное число групп равно числу сочетаний .
Если в общем решении системы свободные неизвестные приравнять нулю, то такое решение называется базисным.
Таким образом, базисное решение рассматриваемой системы таково: х1 = 3/7, х2= - 53/56, x3 = 1/4, x4 = 0.
Для базисного решения получаем разложение вектора B по векторам базиса A1, А2, А3
3/7A1 - 53/56А2 + 1/4А3 = B.
Базисное решение, содержащее точно т отличных от нуля неизвестных, где т — количество линейно независимых уравнений системы, называется невырожденным.
Если хотя бы одно из базисных неизвестных равно нулю, то базисное решение называется вырожденным.
Решение X=(x1,x2,…,xn) называется допустимым, если оно содержит лишь неотрицательные компоненты, т.е. xj³0 для любых j=1,2,…,n. В противном случае решение называется недопустимым.
В ЗЛП особый интерес представляют допустимые базисные решения, называемые опорными планами.
На методе Жордана-Гаусса основаны все известные варианты симплексных методов и, в частности, симплексный метод, который рассматривается ниже.
Контрольные вопросы