Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

1025

.pdf
Скачиваний:
0
Добавлен:
16.06.2024
Размер:
1.06 Mб
Скачать

Определение 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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]