
1025
.pdfОпределение 4. Основной (канонической) задачей ЛП называется задача, которая состоит в определении максимального значения функC ции (2.1) при выполнении условий (2.3) и (2.4), где k=0 и l=n.
Определение 5. Совокупность чисел Х = (х1, х2, ..., хn), удовлетвоC ряющих ограничениям задачи (2.2)–(2.4), называется допустимым реC шением (или планом).
Определение 6. План Х*, при котором целевая функция задачи (2.1) принимает свое максимальное (минимальное) значение, называC ется оптимальным.
Значение целевой функции при плане Х обозначим F(Х). СледоваC тельно, если Х* – оптимальный план, то для любого Х выполняется F(Х) F(Х*).
Переход между формами ЗЛП:от минимума к максимуму;
от ограничений неравенств к ограничениям равенств.
Если требуется найти минимум F, то можно перейти к поиску макC симума – F.
ОграничениеCнеравенство вида можно преобразовать в ограничеC ниеCравенство добавлением к его левой части дополнительной неотриC цательной переменной, а ограничение вида – в ограничениеC равенство вычитанием из его левой части дополнительной переменной. Число вводимых дополнительных неотрицательных переменных при преобразовании ограничений неравенств в ограничения равенства равC но числу преобразуемых неравенств. Вводимые дополнительные переC менные имеют определенный экономический смысл (объем неиспольC зуемого ресурса).
Если переменная хk не подчинена условию неотрицательности, то ее следует заменить двумя неотрицательными переменными uk и vk, приняв хк=uk –vk.
Векторная форма записи ЗЛП имеет вид: max F( Х* )=CX
при ограничениях A1x1+A2x2+…+Anxn=B, x 0,
где C=( c1 ,c2,,…,cn), Х = (х1, х2, ..., хn); CX – скалярное произведение векC торов C, X; Aj и B – векторCстолбцы.
Матричная форма записи ЗЛП:
max F(Х *)=CX
при условиях AX=B, X 0.
Здесь C=( c1, c2, …, cn) – векторCстрока; А=(aij) – матрица размерноC сти m n, столбцами которой являются векторCстолбцы Aj.
21
2.2. Математический аппарат
При изучении математических и прикладных моделей необходимы знания основных понятий высшей математики, матричной и векторной алгебры.
2.2.1. Действия над матрицами
Прямоугольная таблица чисел (действительных или комплексных)
a11 |
a12 |
a13 |
a1n |
|
|
|
A a21 |
a22 |
a23 |
a2n |
, |
(2.5) |
|
|
|
|
|
|
|
|
|
am2 |
am3 |
|
|
|
|
am1 |
amn |
|
|
состоящая из т строк и n столбцов, называется матрицей размера т n. Числа aij i 1,2, m, j 1,2,3, n , составляющие данную матрицу,
называются ее элементами. Здесь первый индекс i обозначает номер строки элемента, а второй j – номер его столбца.
Для матрицы (2.5) часто употребляется сокращенная запись
|
|
i 1,2, m, j 1,2,3, n . |
A aij |
Если число строк матрицы равно числу столбцов ( m n ), то матC рица называется квадратной порядка п. Если же m n , то матрица наC зывается прямоугольной. В частности, матрица типа 1 n называется
вектором строкой, а матрица типа m 1 – векто ром столбцом. Множество всех элементов квадратной матрицы, котоC рые лежат на отрезке, соединяющем левый верхний угол и правый нижний угол, называется главной диагональю, а на отрезке, соединяюC щем верхний правый угол с левым нижним, – побочной диагональю. Квадратная матрица вида:
1 |
0 |
0 |
|
0 |
|
|
|
|
2 |
|
|
|
|
A |
0 |
0 |
0 |
(2.6) |
||
|
|
|
||||
|
0 |
0 |
0 |
|
|
|
|
n |
|
||||
называется диагональной и |
обозначается |
кратко |
1, 2, n , где в |
скобках указаны элементы, находящиеся на главной диагонали.
22
В случае, если i 1 |
(i 1,2, m) , матрица (2.6) называется единич |
||||||
ной и обозначается обычно буквой E, т. е. |
|
|
|||||
|
|
1 |
0 |
0 |
|
0 |
|
E |
0 |
1 |
0 |
|
0 |
|
|
|
|
|
|
|
. |
||
|
|
|
|
||||
|
|
0 |
0 |
0 |
|
1 |
|
Матрица, все элементы которой равны нулю, называется нулевой и обозначается через 0.
С квадратной матрицей связано понятие определителя (детер минанта)
a11 |
a12 |
a13 |
a1n |
|
det A a21 |
a22 |
a23 |
a2n . |
|
|
|
|
|
|
|
an2 |
an3 |
|
|
an1 |
ann |
Не следует отождествлять эти два понятия: матрица представляет собой упорядоченную систему чисел, записанную в виде прямоугольC ной таблицы, а ее определитель det A есть число.
Равенство матриц
Две матрицы A aij и B bij считаются равными: A B , если они имеют одинаковое число строк и столбцов и соответствующие элементы их равны, т.е.
|
|
|
|
aij |
bij . |
|
|
|
|
|
|
Сумма и разность матриц |
|
||||||
Суммой двух матриц A |
|
|
и |
|
|
|
|||
aij |
B bij называется матрица |
||||||||
|
|
элементы которой cij |
равны суммам соответствующих элеC |
||||||
C cij , |
|||||||||
ментов aij |
и bij матриц A и B, т. е. cij |
aij bij . Таким образом, |
|||||||
|
|
a11 b11 |
a12 b12 |
a1n b1n |
|||||
|
|
A B a21 b21 |
a22 b22 |
a2n b2n . |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
am2 bm2 |
|
|
|
||
|
|
am1 bm1 |
amn bmn |
Аналогично определяется разность матриц:
23
a11 b11 |
a12 b12 |
a1n b1n |
||||
A B a21 b21 |
a22 b22 |
a2n b2n . |
||||
|
|
|
|
|
|
|
|
|
am2 bm2 |
|
|
|
|
am1 bm1 |
amn bmn |
|||||
Умножение матрицы на число |
|
|||||
|
|
|
|
|
|
|
Произведением матрицы A aij на число (или произведением
числа на матрицу A) называется матрица, элементы которой полуC чены умножением всех элементов матрицы A на число , т. е.
a11 |
a12 |
a13 |
a1n |
||
|
|
|
|
|
|
A A |
a21 |
a22 |
a23 |
a2n . |
|
|
|
|
|
|
|
|
|
am2 |
am3 |
|
|
am1 |
amn |
Умножение матриц
Пусть
a11 |
a12 |
a13 |
a1n |
|
||
A a21 |
a22 |
a23 |
a2n |
; |
||
|
|
|
|
|
|
|
|
am2 |
am3 |
|
|
|
|
am1 |
amn |
|
||||
b11 |
b12 |
b13 b1q |
|
|
||
b21 |
b22 |
b23 |
b2q |
|
|
|
B |
|
|
|
|
– |
|
|
|
|
|
|||
bp1 |
bp2 |
bp3 bpq |
|
|
матрицы размерностью m n и p q . Если число столбцов матрицы A равно числу строк матрицы B, т. е.
n p ,
то для этих матриц определена матрица C размерностью m q , назыC ваемая их произведением:
24
c |
c |
c |
|
|
11 |
12 |
13 |
c21 |
c22 |
c23 |
|
C |
|
|
|
|
|||
c |
c |
c |
|
|
m1 |
m2 |
m3 |
c1q
c2q
.
cmq
Из определения вытекает следующее правило умножения матриц:
чтобы получить элемент, стоящий в i й строке и j м столбце произве дения двух матриц, нужно элементы i й строки первой матрицы умно жить на соответствующие элементы j го столбца второй и полученные произведения сложить.
Произведение AB имеет смысл тогда и только тогда, когда матрица A содержит в строках столько элементов, сколько элементов имеется в столбцах матрицы B. В частности, можно перемножать квадратные матрицы лишь одинакового порядка.
Транспонированная матрица
Заменив в матрице
a11 |
a12 |
a13 |
a1n |
|
A a21 |
a22 |
a23 |
a2n |
|
|
|
|
|
|
|
am2 |
am3 |
|
|
am1 |
amn |
размерностью m n строки соответственно столбцами, получим так наC зываемую транспонированную матрицу:
a11 |
a12 |
a13 |
a1m |
|
A AT a21 |
a22 |
a23 |
a2m |
|
|
|
|
|
|
|
an2 |
an3 |
|
|
an1 |
anm |
размерностью m n . В частности, для вектораCстроки a a1,a2, ,an
транспонированной матрицей является векторCстолбец
a1 a a2 .
an
25

Если матрица A —квадратная, то
det A det A.
Матрица A aij называется симметрической, если она совпадает со своей транспонированной, т. е. если
A A .
Из равенства вытекает, что: 1) симметрическая матрица – кваC дратная (m=n) и 2) элементы ее, симметричные относительно главной диагонали, равны между собой, т. е.
aij a ji .
Обратная матрица
О п р е д е л е н и е 1 . Обратной матрицей по отношению к данной называется матрица, которая, будучи умноженной как справа, так и слева на данную матрицу, дает единичную матрицу.
Для матрицы A обозначим обратную ей матрицу через A 1 . Тогда по определению имеем:
AA 1 A 1A E ,
где E — единичная матрица.
Нахождение обратной матрицы для данной называется обращением данной матрицы.
Матрица называется невырожденной, если её столбцы линейно неC зависимы. Квадратная матрица обратима тогда и только тогда, когда она невырожденная. Квадратная матрица называется неособенной, если определитель ее отличен от нуля.
В противном случае матрица называется особенной, или сингу лярной.
Теорема. Всякая неособенная матрица имеет обратную матрицу.
Для данной матрицы A ее обратная матрица A 1 единственна. ОсоC бенная квадратная матрица обратной не имеет.
Основные свойства обратной матрицы.
1. Определитель обратной матрицы равен обратной величине опре делителя исходной матрицы.
det A 1 det1 A .
26
2. Обратная матрица произведения квадратных матриц равна про изведению обратных матриц сомножителей, взятому в обратном по рядке, т. е.
AB 1 B 1A 1 .
Спомощью обратной матрицы легко решаются матричные уравнеC
ния
AX B .
Действительно, если det A 0 , то
X A 1B
Вычисление определителей
Элементарные преобразования матрицы дают наиболее удобный способ вычисления определителя этой матрицы. Пусть, например,
|
a11 |
|
a12 |
a13 |
a1n |
|
|||
n |
a21 |
a22 |
a23 |
a2n |
|
||||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
an2 |
an3 |
|
|
|
|
an1 |
ann |
|
||||||
Предполагая, что a11 0 , будем иметь: |
|
|
|
||||||
|
|
1 a12 |
a13 |
a1n |
|
||||
|
|
a |
|
|
|
|
|
|
|
|
|
|
21 |
a22 |
a23 |
|
a2n |
|
|
|
|
|
|||||||
n a11 |
a |
|
|
|
|
|
|||
|
11 |
|
|
|
|
|
|||
|
|
|
|
|
|||||
|
|
a |
n1 |
an2 |
an3 |
|
ann |
|
|
|
|
|
|
|
|||||
|
|
|
|
||||||
|
|
a11 |
|
|
|
|
|
Отсюда, вычитая из элементов aij , принадлежащих jCму столбцу ( j 2 ), соответствующие элементы первого столбца, умноженные на aij , получим:
|
|
|
1 |
0 |
0 |
|
0 |
|
|
|
|
a |
1 |
1 |
|
1 |
|
|
|
|
|
|
21 |
a22 |
a23 |
|
a2n |
|
|
|
|
|
|
||||||
n |
a11 |
a |
|
|
|
|
|
a11 n 1 , |
|
|
11 |
|
|
|
|
|
|||
|
|
|
|
|
|
||||
|
|
a |
1 |
1 |
|
1 |
|
|
|
|
|
|
n1 |
an2 |
an3 |
|
ann |
|
|
|
|
|
|
||||||
|
|
a11 |
|
|
|
|
|
|
27

где
|
a221 |
a231 |
a241 |
a21n |
|||||
|
|
|
|
a331 |
a341 |
|
|
||
n 1 |
a321 |
a31n |
|||||||
|
|
|
|
|
|
|
|||
|
|
|
|
||||||
|
|
|
|
|
|
|
|
||
|
|
1 |
|
|
1 |
|
1 |
|
1 |
|
an2 |
an3 |
an4 |
ann |
|||||
и |
|
|
|
|
|
|
|
|
|
aij1 |
aij |
ai1a1 j |
|
i, j 2,3, ,n |
|||||
|
|||||||||
|
|
|
|
a |
ii |
|
|
|
|
|
|
|
|
|
|
|
|
||
К определителю n 1 применяем тот же прием. Если все элементы |
|||||||||
aiii 1 0 |
i 1,2, ,n 1 , |
то окончательно находим:
n a11a221 annn 1 .
Если в какомCнибудь промежуточном определителе n k левый верхний элемент akk1,k 1 0 , то следует переставить строки или столбC цы определителя n k так, чтобы нужный нам элемент был отличен от
нуля (это возможно всегда, если определитель 0). Конечно, при этом следует учесть изменение знака определителя n k . Можно дать
более общее правило. Пусть определитель n det ij преобразован так, что pq 1 ( pq — главный элемент), т. е.
11
i1n
p1n1
Тогда
1q 1 j
iq ij
1 rj
nq nj
1n
in
.
pn
nn
n 1 p q n 1,
28

1 |
есть определитель n 1 Cго порядка, получаюC |
где n 1 det ij |
щийся из n путем выбрасывания pCй строки и qCго столбца с послеC дующим преобразованием элементов по формуле
ij1 ij iq pj ,
т. е. каждый элемент aij1 определителя n 1 , равен соответствующему элементу определителя ij , уменьшенному на произведение его «проек ций» iq и pj на отброшенные столбец и строку исходного определи
теля.
Заметим, что число умножений и делений, нужных для вычисления определителя n го порядка, равно:
n 1 n2 n 3 . 3
2.2.2. Решение систем линейных уравнений
Способы решения систем линейных уравнений в основном делятся на две группы: 1) точные методы, представляющие собой конечные алC горитмы для вычислений корней системы (таковы, например, правило Крамера, метод Гаусса, метод главных элементов, метод квадратных корней и др.), и 2) итерационные методы, позволяющие получать корC ни системы с заданной точностью путем сходящихся бесконечных процессов (к числу их относятся метод итераций, метод Зейделя, ме тод релаксации и др.).
Решение систем с помощью обратной матрицы. Формулы Крамера
Пусть дана система n линейных уравнений с n неизвестными:
a11x1 a12 x2 a1n xn b1, |
|
|
||||||||||
a21x1 a22 x2 a2n xn b2 |
, |
(2.7) |
||||||||||
|
|
|
|
|
||||||||
|
|
|
|
|
||||||||
a |
x |
a |
n2 |
x |
a |
x |
b . |
|
||||
|
n1 1 |
|
2 |
|
|
nn |
n |
|
n |
|
||
Обозначим через |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a11 |
a12 |
a13 |
a1n |
|||||||
|
A |
a21 |
a22 |
a23 |
a2n |
|||||||
|
|
|
|
|
|
|
|
|||||
|
|
|
|
am2 |
am3 |
|
|
|
|
|||
|
|
am1 |
amn |
29
матрицу из коэффициентов системы (2.7), через
|
b1 |
|
|
b |
b |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
bn |
|
|
столбец ее свободных членов, через |
|
|
|
|
x |
|
|
|
|
1 |
|
x |
x2 |
||
|
|
|
|
|
|
|
|
|
xn |
|
столбец из неизвестных (и с к о м ы й |
в е к т о р ). Тогда система (2.7) |
кратко может быть записана в виде матричного уравнения |
|
Ax b |
(2.8) |
Совокупность чисел x1, x2, , xn (или, короче, вектор x), обращаюC |
щих систему (2.7) в тождество, называется решением этой системы, а сами числа xi —ее корнями.
Если матрица— неособенная, т. е.
a11 |
a12 |
a13 |
a1n |
|
|
det A a21 |
a22 |
a23 |
a2n |
0 , |
(2.9) |
|
|
|
|
|
|
an1 |
an2 |
an3 |
ann |
|
|
то система (2.7), или эквивалентное ей матричное уравнение (2.8), имеет единственное решение.
В самом деле, при условии det A 0 существует обратная матрица
A 1 . Умножая обе части уравнения (2.8) слева |
на матриC |
цу A 1 , получим: |
|
A 1Ax A 1b |
|
или |
|
x A 1b |
(2.10) |
Формула (2.10), очевидно, дает решение уравнения (2.8), причем так как каждое решение имеет вид (2.10), то решение единственно.
30