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

Прикладная математика. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

 

Введя дополнительные переменные

и

и переходя от

 

к

 

, где

 

сформулируем эту,

задачу следующим

образом.

min

,

 

 

max

 

 

 

 

, , , ,

,

удовле-

Найти неотрицательные значения переменных

творяющие системе уравнений:

 

 

 

428;

10;

(1.62)

324,

для которых функция

3

2

(1.63)

принимает наименьшее значение.

Для возможности применения симплекс-метода необходимо выразить первые т неизвестных через остальные. Здесь выгоднее последние неизвестные выразить через первые переменные , . Перепишем систему (1.62) в виде

 

28

 

4 ;

 

(1.64)

 

10

 

;

 

 

Базисное решение

системы (1.62) имеет вид

 

 

24

3

.

 

 

Очевидно, оно0,

,

28,

10,

24.

(1.65)

0

с учетом

является допустимым; при этом

0коэффици-

(1.65). Согласно симплекс-методу из всех отрицательных

ентов функции (1.63) выбираем наиболее отрицательный, т.е. наибольший по модулю. У нас это коэффициент при ; его берем в ка-

честве

. Поэтому имеем случай 3. Сделаем свободную перемен-

ную

ɤбазисной. Чтобы установить, какую базисную переменную

нужно при этом сделать свободной, определим ведущий коэффициент системы (1.64). Так как

min | | ,| | ,| | | | ,

то ведущим коэффициентом является коэффициент –3 в последнем уравнении системы (1.64) и, следовательно, базисную переменную

надо сделать свободной.

Разрешая последнее уравнение системы (1.64) относительно и исключая затем эту переменную из остальных уравнений системы (1.64) и целевой функции (1.63), находим:

31

 

 

8

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

;

 

 

(1.66)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

;

 

 

 

(1.67)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Новое допустимое

базисного решения системы (1.66) имеет вид

 

24

 

 

.

 

 

 

(1.68)

8,

0

,

20,

2, .

 

0.

 

 

 

Новое значение целевой функции

 

 

 

 

 

 

 

Поскольку значение

целевой функции уменьшилось с нуля до

–24,

то решение (1.68)

24

 

 

 

лучше решения (1.65).

Далее процесс решения повторяем. Имеем случай 3. Переменную х2 переводим в разряд базисных, так как при в (1.66) все коэффициенты отрицательные. Составляем ряд для определения ведущего коэффициента. Он определяется из условия

 

 

 

 

min

 

 

,

 

 

,

 

 

 

 

 

 

 

.

 

 

при пере-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ведущим коэффициентом является коэффициент

 

 

 

менной

в третьем уравнении системы (1.66). На основании этого

коэффициента нужно сделать базисную переменную

свободной.

Разрешая третье уравнение системы (1.66) относительно

и исклю-

чая

 

из остальных уравнений системы (1.66) и целевой функции

(1.67), получаем

3

 

 

1,5

 

 

0,5

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.69)

 

 

 

 

 

 

7

 

 

0,5

 

 

0,5

 

;

 

 

 

 

Так как все

 

 

9

 

 

5,5

 

 

1,5

 

;

 

 

(1.70)

 

 

 

27

1,5

 

0,5

 

.

положительные

 

 

 

 

коэффициенты при переменных

 

 

 

(

γ

), то имеем случай 1.

Очевидно, допустимое,

базисное реше-

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ние системы (1.69) имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

(1.71)

 

 

 

 

7,

 

3

,

 

9,

. 0,

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

Новое значение функции

 

 

 

 

 

 

Так как значение целевой

27

функции уменьшилось (с –24 до –27), то решение (1.71) лучше решения (1.68). Процесс решения снова повторяем. Так как имеем случай 1,

32

то решение (1.71) является оптимальным для задачи (1.62), (1.63). Таким образом, допустимое оптимальное решение исходной задачи – это 7, 3, аmax 27.

1.6. Выбор допустимого базисного решения

Мы решили систему уравнений (1.41) относительно т неизвестных с помощью симплекс-метода при выполнении условий (1.46). В результате получили решение в виде (1.44). Теперь предположим,

что некоторая переменная

входит только в первое уравнение сис-

темы (1.41) со знаком плюс, переменная

входит во второе урав-

нение этой системы уравнений со знаком плюс, и так далее и, нако-

нец, некоторая переменная

входит в m-е уравнение системы

(1.41) со знаком плюс.

 

 

 

.

Разрешаем систему (1.41) относительно переменных

 

Поскольку коэффициенты при переменных были

положительными,

 

, ,…,

 

то сразу находим неотрицательное базисное решение. Именно этот случай имел место в предыдущем примере.

Один или несколько коэффициентов «b» в системе уравнений

(1.41) могут быть отрицательными, например 0. Тогда надо умножить все соответствующее уравнение в (1.41) на (–1), а рассуждения повторить. Для выбора допустимого базисного решения посту-

пают

следующим

образом. Вводят дополнительные переменные

 

 

 

 

и при помощи рассмотренного симплекс-метода

находят такие неотрицательные значения переменных

, ,…, ,

ƺ

,…,ƺ

, удовлетворяющие уравнениям

 

ƺ

,…,ƺ

 

,

1,2,…,

,

(1.72)

для которых ƺ

 

 

 

 

 

вспомогательная целевая функция

 

 

принимает

 

Φ

∑ ƺ

, ƺ

0

 

(1.73)

 

 

 

 

 

 

наименьшее значение. В этом случае базисное решение

системы (1.72) выглядит так:

 

0.

 

 

 

 

 

,

,…,

,

 

 

Применяя шаг за шагом симплекс-метод Данцига, можно получить оптимальное допустимое решение. Минимальное значение функции (1.73) в точности равно нулю, когда каждая из ƺ в точности равна нулю, т.е.

min ∑ ƺ

0.

33

Тогда все базисные переменные

перешли в разряд свободных

или независимых переменных и,

следовательноƺ

, их надо полагать

равными нулю. Поэтому оптимальное допустимое решение имеет вид

 

ƺ

 

β′

α′

,

1,2,…, ,

а все

0

, т.е. это означает, что

β′

0 .

 

 

 

Рассмотрим такой случай: минимум функции (1.73) может быть больше нуля. Покажем, что при этом допустимого базисного решения не существует. Доказательство проведем от противного. Пусть у системы (1.41) существует допустимое базисное решение, полученное с помощью минимизации функции (1.73). Сличим между собой системы уравнений (1.41) и (1.72). Смысл формулы (1.72): вычитание

из bk

левых частей (1.41) и получение

. Очевидно, что система

(1.41)

совпадает с системой (1.72). Следовательноƺ

, допустимое ба-

зисное решение системы (1.41) будет допустимым базисным реше-

нием (1.72) при

 

 

. Но при

 

 

0

имеем, что это допусти-

мое решение

обеспечивает достижение

 

 

 

. И оно же является

 

ƺ

0

 

 

ƺ

 

 

0

 

 

 

допустимым для (1.72), при этом

 

minΦ.

 

 

 

 

 

 

 

 

 

 

 

 

 

Имеем противоречие.

Следовательно, в этом случае

допустимого базисного решения у сис-

 

minΦ

0

 

 

minΦ

 

темы (1.41) не существует. Это очевидно, т.к. если

, то это

означает, что некоторые ƺ

0 .

 

 

 

 

 

 

0

Контрольные вопросы

1.Как исходя из неравенства (1.22) определить необходимое число итераций для получения решения с заданной точностью?

2.Сколько допустимых решений получается при решении задачи линейного программирования в типичных случаях?

3.Найдите неотрицательные значения переменных x1 и x2, удовлетворяющие неравенствам

 

9

4

15 ;

 

3

27 ;

 

 

5

20 ,

длякоторыхфункция

 

 

принимаетнаименьшеезначение.

4. Найдите значения

переменных x

4

2

1 и x2, удовлетворяющие нера-

венствам

34

 

 

 

4

 

28;

 

 

 

 

 

3

 

10 ;

0 ;

 

0 ,

длякоторыхфункция

 

 

24;

 

 

 

 

принимаетнаибольшеезначение.

5. Решите симплекс-

методом следующие задачи:

 

2

 

2

 

 

 

 

а)

,

→ max

б)

,

 

3

→ max

при

 

 

 

2

16 ;

при

3 5

2 4

3 ;

9 ;

 

 

2 ;

 

 

0 ,

1,2,3 .2

 

5 ;

 

,0 .

6.Симплекс-методом решите задачу

2→ min

при ограничениях

3 2 12;

28 ;

2 3 6 ; , 0 .

7. Решите симплекс-методом задачу

4→ max

при ограничениях

4 8 0;

22 0 ;

2

10

0;

,

0 .

2

0 ;

8.Приведите примеры действий над матрицами.

9.Как находится обратная матрица?

10.Сформулируйте постановки двух основных задач линейной алгебры.

11.Каковы недостатки и преимущества известных Вам методов решения основных задач линейной алгебры? Сравните эти методы.

12.В чем состоит идея решения системы алгебраических уравнений методом обратной матрицы?

13.Охарактеризуйте и приведите формулы Крамера.

14.Опишите и объясните метод Гаусса.

15.В чем заключается метод простой итерации для линейных алгебраических систем?

16.Приведите достаточное условие сходимости метода простой итерации.

35

17.Что такое норма матрицы? Как она вычисляется?

18.Приведите решение второй основной задачи линейной алгебры.

19.Сформулируйте постановку основной задачи линейного программирования.

20.Как осуществляется сведение простых задач линейного программирования к каноническому виду?

21.Первый случай симплекс-метода Данцига, его смысл и основные формулы.

22.Идея и основные формулы второго случая симплекс-метода Данцига.

23.Третий случай симплекс-метода Данцига, его смысл и вывод основных формул.

24.Придумайте и решите конкретную задачу линейного программирования.

25.Дайте обоснование выбора допустимого базисного решения при b2 < 0 в системе уравнений (1.41).

26.Может ли целевая функция (1.73) иметь минимальное значение, большее нуля?

27.Почему минимум функции (1.73) должен быть равен нулю? Докажите это.

36

2. ЧИСЛЕННЫЕ И ПРИБЛИЖЕННЫЕ МЕТОДЫ ПРИКЛАДНОЙ МАТЕМАТИКИ

2.1. Решение трансцендентных уравнений

Уравнение вида

 

 

произвольная функция0,от

 

(2.1)

где

, называется трансцендент-

ным уравнением.

sin

10

. Как частный слу-

Это, например, уравнение

чай,

это

будет уравнение

0

 

 

 

алгебраическое.

Предположим, что

дважды непрерывно дифференцируемая функция.

Чтобы найти решение уравнения (2.1), надо найти такую постоянную ƺ, которая будучи подставлена в (2.1), обратит его в равенство ƺ 0. Другими словами, это значит найти корни уравнения (2.1). Так как точные значения корней уравнения (2.1) в общем случае найти невозможно, то их находят приближенно, т.е. будут погрешности. Приближенное нахождение корней уравнения (2.1) обычно осу-

ществляется в два этапа.

1-й этап. Нахождение грубых значений корней.

Обозначим через

отрезок, в котором имеется только один

корень ƺ уравнения (2.1),,

т.е.

 

. Предположим, что этот ко-

рень вещественный и простой, т.е.ƺ b

 

. Любое число из данно-

го отрезка можно взять за корень. Чемƺ

меньше отрезок, тем точнее

 

0

корень. Приближенное решение записывается в виде

Грубые, с 1–2 десятичнымиƺ ƺ

знаками.

, значения действительных

корней уравнения (2.1) можно найти графически. С этой целью стро-

ят график функции

и абсциссы точек пересечения кривой

с осью

принимают за приближенные значения корней

уравнения (2.1). При этом отрезок

берется как можно меньше;

важно, чтобы только корень лежал на,этом отрезке (рис. 2.1).

Здесь предполагали, что отрезок

выбирается так, чтобы на

нем находился точно один корень. Если,

будут два корня, метод рабо-

тать не будет, т.е. корни должны быть однократными.

37

Рис. 2.1. Графическое решение уравнения (2.1)

В некоторых конкретных случаях уравнение (2.1) выгоднее представить в виде

 

 

 

 

 

 

 

 

 

(2.2)

 

 

 

 

выбирать так, чтобы их графики были в

причем эти функции надоφ

Ψ

,

 

построении проще функции

.

 

 

 

 

Пример. Решим уравнение

1

 

 

0.

 

Представим его в виде

lg

 

 

 

Окончательно через

соотношение (2.2) получим

 

 

 

 

 

lg

1.

 

где

 

 

 

(рисlg. 2.2)

 

,

 

 

 

 

 

.

 

 

 

 

 

 

 

2φ-й этап.lg ; Ψ

 

 

 

 

 

 

 

 

 

Уточнение грубых значений корней.

 

Так как корни ищем приближенно, то возникает вопрос: а как же

далек приближенный корень от истинного значения

? Другими сло-

вами, насколько велика погрешность

 

?

ƺ

ƺ

на

 

,

приближенное

значение|

корня

Если

 

 

, то

 

называется невязкой для . Очевидно, для корня невязка рав-

на нулю. Поскольку

не совпадает с ,

то

в каком-то смысле

близко к

. По величине невязки нельзяƺ

судить о точности найденно-

го приближенногоƺ

корня, так как часто говорят, что если невязка мала,

токореньнайденточно. Покажемэтонапримере.

38

(х),

 

(х)

lgх

 

Рис. 2.2. Иллюстрация графического метода

Пример. Имеем уравнение

2 1 0. Тогда корень

ƺ1 точный корень. Запишем уравнение в виде

Возьмем в

10

 

 

2

 

10

10

.

 

 

.

 

 

 

 

 

 

0

100

 

 

 

 

качестве приближенного решения:

 

 

 

 

 

В этом случае невязка составит

 

 

 

 

 

 

Но это 10

 

10

2

 

10 100

10

10

.

от

 

 

от-

 

 

100

ƺ

1

мизерная величина, а в то же время

 

 

 

 

личается очень сильно.

 

 

 

 

 

 

 

 

 

2.2. Теорема об оценке погрешности корня через невязку

Пусть функция

 

достаточно

гладкая, т.е. имеет столько

производных, сколько надо. Рассмотрим величину

ƺ

. По

теореме Лагранжа о среднем имеем

 

 

,

 

где

.

ƺ

= ( ƺ

)

 

 

(2.3)

Переходяƺ

в (2.3) к модулям и учитывая очевидное соотношение

| ƺ | 0, получим

, последнее| |ƺ

||

|

.

 

 

 

|

 

 

 

Поскольку

|

равенство позволяет написать

 

причем

|

min|

||ƺ

 

|,

 

 

 

.

 

 

 

 

 

 

 

39

Откуда для погрешности ƺ

можно дать следующую

оценку:

(2.4)

 

min

 

|

|

|

,

 

 

 

где

|

.

 

 

 

 

 

 

 

Первая производная|

 

на отрезке

,

в нуль не обращается.

Поэтому имеет место приведенная оценка

(2.4). Однако эта оценка

непрактичная, т.е. корень будет ближе к истинному, чем дает оценка (2.4). Это связано с тем обстоятельством, что при оценке погрешности в (2.4) делим на минимум производной функции . Для этого случая введено понятие плохо обусловленных уравнений. Если

на отрезке

очень мала, то уравнение называется плохо

обусловленным. Если,

уравнение является плохо обусловленным, то

сходимость решения плохая, а раз так, то и метод ведет себя плохо. А это означает, что ошибка велика, так как делим невязку на малую величину.

2.3. Метод итераций для уравнения

Предположим, что на отрезке

 

находится единственный ко-

рень и нам удалось записать уравнение,

(2.1) в виде

 

 

 

 

 

 

 

(2.5)

где

 

достаточно гладкая

функция, т.е. у нее существуют произ-

φ

 

φ

,

 

 

 

 

 

 

 

 

водные на отрезке

.

 

 

 

 

Решить уравнение (2.5)ƺ

– это значит найти его корень:

 

 

 

 

=

 

 

(2.6)

Соотношение (2.6) есть

точное

равенство.

ƺ

 

φ ƺ .

 

Возьмем приближенное значение корня

,, .

Организуем процесс решения в виде цепочки итераций:

φ ,

φ ,…,

.

(2.7)

Если этот процесс сходится, то он сходится к точному решению задачи (2.1). Действительно, если процесс (2.7) сходится, то

lim

ƺ.

→ ∞

40

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