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

1017

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

6

 

5

 

 

 

 

x2

 

 

 

1 x

 

 

 

x 1

5

0

5

x

 

 

 

 

 

5

 

6

 

 

 

 

6

x x x1

6

 

Рис 1. График исследуемой функции

 

Порядок выполнения работы

В процессе проведения лабораторной работы студент должен выполнить следующие задания:

1.Ознакомиться с теоретическим материалом.

2.Исследовать функцию, заданную в табл. 4 в соответствии с номером варианта, методами дифференциального исчисления и построить график.

3.Составить отчет.

4.Защитить работу.

 

 

 

Таблица 4

№ вари-

f(x)

№ вари-

f(x)

анта

анта

 

 

1

2

3

4

1

 

y

x3

 

11

y

 

x3

 

 

x2

1

3

x2

 

 

 

 

 

 

 

 

 

 

 

 

x

3

 

 

 

 

 

 

 

 

4

 

 

 

2

y

 

 

 

 

 

12

y

x

 

 

 

2 x 1 2

 

 

 

 

 

 

x3 1

3

y

 

 

1

 

 

 

13

y

 

 

2x 1

x 2

x 2 1

 

x 1 2

 

11

Окончание табл. 4

1

 

 

 

 

2

 

 

 

 

 

3

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

y

 

x2

1

14

y

 

4x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2

1

 

x3

1

 

 

 

 

 

 

 

5

y

 

 

 

x

 

15

y

 

 

x3

 

1

x2

 

x 1

 

 

 

 

 

 

6

y

 

x2

1

16

y

 

2 4x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2

1

1

4x2

 

 

 

 

 

 

7

y

 

 

 

x

 

17

y

x 2 x 1

 

 

 

 

x

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 1

8

y

 

x 2

 

18

y

 

 

4x 2

 

 

x

1

3

x 2

 

 

 

 

 

 

9

y

 

 

x2

 

19

y

 

x2

3x 3

 

 

x

2

1

 

 

 

x 1

 

 

 

 

 

 

 

 

 

10

y

 

 

x3

 

20

y

 

x 2

4x 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2

1

 

 

 

x 4

 

 

 

 

 

 

 

 

 

21

y

 

 

2

 

 

 

 

 

 

26

y

 

 

8 x 1

 

 

 

 

 

 

 

 

 

 

 

x 1 2

 

 

x2 2x

 

22

y

 

12x

 

27

y

 

 

 

4

 

 

 

 

 

 

 

9 x2

 

x 2 2x 3

 

 

 

 

 

 

 

23

y

 

4 x3

 

28

y

 

x3

4

 

 

 

x2

 

 

x2

 

 

 

 

 

 

 

 

 

24

y

 

 

 

x 2

 

29

y

 

x 2

6x 9

 

 

x

1 2

 

 

x 1 2

 

 

 

 

 

 

 

 

25

y

12 3x2

30

y

 

3x4 1

 

x2 12

 

 

 

 

x3

 

 

 

 

 

 

 

 

 

 

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

1.Из каких этапов состоит алгоритм исследования функции?

2.Как определить интервалы монотонности функции?

3.Что такое точка перегиба и как ее вычислить?

4.Как определить интервалы выпуклости – вогнутости функции?

5.Какназываетсяфункция, котораянеявляетсячетнойинеявляетсянечетной?

6.Графиккакойфункциисимметриченотносительноначалакоординат? Приведите примеры.

7.Какая асимптота является частным случаем наклонной асимптоты?

12

Лабораторная работа 2. Исследование функции градиентным методом

Целью проведения студентом лабораторной работы является:

– формирование умения и навыков градиентного исследования функций.

Теоретические основы

Рассмотрим некоторую точку X 0 из области определения функции и вектор D d1 , d2 , , dn , задающий направление в пространстве Rn :

D 1, di 0 .

Совокупность точек X X 0 D (τ – некоторый одномерный параметр) образует прямую, если , и луч, если 0 .

Градиентом функции x1 , x2 , , xn называют вектор частных производных

 

 

 

 

 

 

 

x

 

 

X

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

x2

 

(1.1)

X

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

x

 

 

 

 

 

 

 

 

 

n

 

 

 

 

Обозначение: grad (X) или X

, – оператор набла.

 

Вектор X , противоположный вектору градиента,

называется ан-

тиградиентом функции.

Точка, в которой градиентфункции f x(k ) 0 , называется стационар-

ной.

Модуль градиента показывает максимальную скорость изменения функциивокрестностиM0, тоестьчастотулинийуровня. Модульградиента равен наибольшей производной по направлению в данной точке.

Теоремы о функции, имеющей производную по заданному направлению

Пусть D d1 , d2 , , dn – некоторое заданное направление. ВеличинаT X 0 D T (X ), D – скалярное произведение, показывающее

мгновенную скорость изменения функции Φ вдоль направления D и равное производной от функции Φ в направлении D:

T x0 D lim

X 0

D X 0

.

(1.2)

 

 

0

 

 

 

13

Теорема 1.

Производная по направлению D равна произведению градиента на вектор D:

D T X 0 D .

Теорема 2.

Если T X 0 D 0, то существует такая величина σ, при которойX 0 D X 0 для любого 0, .

Теорема 3.

Градиент X 0 показывает направление наискорейшего возрастания

функции Φ в точке X 0 . Теорема 4.

В точке максимума X* функции (X) градиент равен нулю:

X * 0 .

Пример 2.

Задана функция x1, x2 x1 0,5x22 и вектор D 1; 0,5 .

1.Определить единичное направление e , совпадающее с направлением вектора D.

2.Из точки X 0 1;1 провести прямую в направлении D и написать

уравнение этой прямой.

3.Вычислить градиент функции X 0 .

4.Определить производную функции Φ по направлению e : e X 0 .

5.Найти точки, в которых X 0, т.е. точки, подозрительные на

экстремум. Решение.

1. Определим единичный вектор:

 

 

 

 

 

 

1;

0,5

 

1

 

0,5

 

 

 

e

 

D

 

 

 

0,8944;

0,4472 .

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

12

0,52

 

1,116

 

1,116

 

 

 

2. Для построения прямой в направлении вектора D найдем еще одну точку на искомой прямой, например точку X 1 X 0 D, т.е. X 1 x10 d1; x20 d2 , изапишемуравнениепрямой, проходящейчерездве точки в виде:

 

 

 

 

 

x

x0

 

x

2

x0

 

 

 

 

 

1

 

1

 

2

.

 

 

 

 

x1

x0

x1

 

 

 

 

 

 

 

x0

 

 

 

 

1

 

1

 

 

2

2

 

Для точек X 0

1;1 и X 1

2;

1,5 уравнение будет иметь вид:

x1 1

 

x2 1

или

x 0,5x

0,5.

 

 

 

 

 

 

 

 

2 1

1,5 1

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14

3. Градиент функции в точке X x1 , x2 находим по формуле (1.1):

 

 

 

 

 

 

 

 

 

 

 

x

X

1

 

 

 

 

1

 

 

x

.

 

 

X

 

 

 

 

 

 

x

X

2

 

 

 

 

 

 

 

 

2

 

 

 

 

В точке X 0 1,1

1

 

 

 

 

 

X 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

4.

 

1

 

 

 

 

 

Производную функции Φ по направлению e вычислим, исходя из

утверждения теоремы 1:

 

 

 

 

 

 

 

e T X 0 e ,

 

где T X 0 1; 1 ; e 0,8944;

0,4472 .

 

 

Вычислим скалярное произведение двух векторов (как сумму произве-

дений соответствующих координат):

 

 

 

 

e T X 0 e 1; 1 0,8944;

0,4472 0,8944 0,4472 1,3416.

5.

Определим точки, подозрительные на экстремум:

 

 

X

1

 

 

 

 

 

 

0 .

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

Видим, что X не обращается в 0 ни в одной точке (в том числе и

при x2

0 ); следовательно, для данной функции не существует точек, подо-

зрительных на экстремум. Область изменения функции , .

Порядок выполнения работы

В процессе проведения лабораторной работы студент должен выполнить следующие задания:

1.Ознакомиться с теоретическим материалом.

2.Для функции x1, x2 x1 10m x22 выполнить следующее2:

определить ed d – единичное направление, совпадающее с направле-

нием d n, m ;

10

из точки x0 x10 , x20 , где x10 x20 0,01, провести прямую в направле-

 

нии d (написать уравнение этой прямой);

x0 x0

, x0

;

 

 

 

вычислить градиент функции в точке

 

 

 

 

1

2

x0 x0

 

;

 

определить производную по направлению e в точке

, x0

 

 

 

 

1

2

 

2 Здесь и далее n – номер группы, m – номер по списку в журнале.

15

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

0;

найти T d .

3.Составить отчет.

4.Защитить работу.

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

1.Что такое градиент функции, антиградиент?

2.Чему перпендикулярен градиент функции?

3.Чему равен модуль градиента?

2.Как определить производную по направлению D?

3.Как найти точку экстремума с помощью градиента функции?

16

Раздел 2. ЭЛЕМЕНТЫ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ

Лабораторная работа 3.

Решение задачи линейного программирования

Целью проведения студентом лабораторной работы является:

формирование умения и навыков математической формулировки задачи линейного программирования (ЗЛП);

формированиеуменияинавыковпостроениядопустимойобластиЗЛП

илинии уровня;

формирование умения и навыков решения ЗЛП графическим и аналитическим методами.

Теоретические основы

Говорят, что в пространстве Rn поставлена задача линейного программирования, если определена линейная целевая функция f (x), подлежащая мини-

мизации, а также задана система линейных ограниченийна компоненты векторов X x1, x2 , , xn в виде равенств и неравенств, определяющая некоторую

область Rn :

 

n

 

f (x) c j x j min ;

 

j 1

X

 

 

n

 

i I;

aij x j ai ,

j 1

 

 

n

 

 

bkj xj bk , k K;

j 1

 

 

x

0,

l L.

l

 

 

 

 

 

(2.1)

(2.2)

Здесь I, K, L – заданные целочисленные множества индексов. Канонической формой записи ЗЛП называют задачу

 

 

n

 

i 1, , m;

n

 

aij xj ai ,

f (x) c j x j max;

j 1

 

(2.3)

j 1

X

 

0,

j 1, ,n.

 

 

xj

То есть в канонической форме ЗЛП все ограничения системы (2.2) (кроме неравенств, выражающих условие неотрицательности переменных) представлены в виде равенств и требуется найти максимум целевой функции.

17

Приведение задачи к стандартной или канонической форме легко осу-

ществить, если использовать следующие положения.

 

1)

Задача минимизации может быть преобразована в задачу максими-

зации

путем

смены

знака

у

критерия

оптимальности:

min f (x) max f (x) и наоборот, max f

(x) min f (x) .

 

X

X

 

 

X

X

 

2) Ограничения типа равенств и неравенств можно привести к необходимой форме, используя следующие приемы:

а) ограничение типа можно привести к ограничению путем умножения обеих частей неравенства на –1;

б) ограничениетипа= приводитсякдвумограничениямтипанеравенстви , если потребовать их одновременного выполнения. Например, равенство x1 x2 1 эквивалентно системе двух неравенств

x1 x2 1;x1 x2 1;

в) ограничение типа неравенства можно всегда привести к ограничению типа равенства путем добавления еще одной переменной xn 1 в левую часть

неравенства. Так, например, ограничение ai1 x1 ai2 x2 ... ain xn bi

можно привести к ограничению типа равенства добавлением положитель-

ной переменной xn 1 : ai1 x1 ai 2 x2 ... ain xn xn 1 bi , xn 1 0. При этом нужно помнить, что размерность задачи (число переменных хi) станет

на единицу больше, а неизвестная xn 1 войдет в качестве слагаемого в целевую функцию f(x) с коэффициентом cn 1 , равным нулю ( cn 1 0 );

г) если не задано ограничение вида xi 0 на некоторую переменную, а оно необходимо для канонической формы, то вместо одной переменной xi можно рассмотреть две переменные xi и xi , такие, что xi 0 и xi 0 , а

xi xi xi .

Используя перечисленные приемы, произвольную задачу линейного программирования можно привести к основной канонической форме: мак-

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

симизировать целевую функцию f (x) c j

x j

по переменным x1 ,..., xn

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

при выполнении ограничений:

 

 

 

 

 

 

 

 

 

 

a11 x1 a12 x2

... a1n xn b1;

 

a

21

x a

 

x

... a

2n

x

b ;

 

1

22

 

2

 

 

 

n

 

2

 

................................................

 

(2.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

m1

x a

m2

x

2

... a

mn

x

n

b ;

 

1

 

 

 

 

 

 

m

x

j

0,

j 1,...,n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

18

В матричной форме задача записывается следующим образом:

(C, X ) max;

A X T B; (2.5)

X 0.

Всякое решение системы (2.4) называется допустимым решением. Допустимое решение, при котором линейная функция f(x) принимает оптимальное (в данном случае максимальное) значение, называется оптималь-

ным решением или оптимальным планом.

Пример 3.

Привести к канонической форме задачу минимизации функции f 3x1 4x2 при ограничениях:

2x1 x2 5;

x1 x2 3;x1 0.

Чтобы привести задачу к канонической форме, выполним следующее. 1) Задачу минимизации функции f заменим задачей максимизации пу-

тем смены знака:

minx1 , x2 3x1 4x2 maxx1 , x2 3x1 4x2 .

2)Ограничение типа неравенства 2x1 x2 5 приведем к ограничению типа равенства путем добавления переменной x3 0 : 2x1 x2 x3 5 .

3)Вместо переменной x2, для которой не задано ограничение x2 0 , введем в рассмотрение две переменные x2 0 и x2 0 : x2 x2 x2 .

Для удобства записи задачи переименуем все переменные в соответствии с порядком их вхождения в целевую функцию:

x1 x1 , x2 x2 , x3 x2 , x4 x3 .

В канонической форме ЗЛП примет вид:

3x1 4x2 4x3 0x4

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

2x1 x2 x3 x4 5;x1 x2 x3 3;

x1 0, x2 0, x3 0, x4 0.

max

x1 ,x2 ,x3 ,x4

Рассмотрим некоторые типичные задачи линейного программирования.

19

Задача о раскрое

Пусть имеется n вариантов раскроя листа жести на заготовки m типов. Заготовок i-го типа требуется в количестве bi (i 1, .., m) . При j-м способе

раскроя j 1, , n получаются заготовки i-го вида в количестве aij , а отходы (по весу, по площади и др.) равны c j . Обозначим как x j количество листов, раскраиваемых по варианту j. Вектор X (x1 , x2 ,..., xn ) представляет

собой общий план раскроя. Тогда задача по отысканию оптимального плана раскроя может быть записана следующим образом.

Определить

вектор

X,

 

минимизирующий

функцию отходов

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (X ) c j x j

приограничениях, направленныхнавыполнениепланавы-

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пуска заготовок требуемых типов:

 

 

 

 

 

 

 

 

 

 

a

 

x

a

x

2

... a

x

n

b

;

 

11

1

 

12

 

 

 

1n

 

 

1

 

 

..............................................

 

 

 

 

 

x

a

 

x

 

 

... a

 

x

 

 

b ;

 

a

i1

i2

2

in

n

 

 

1

 

 

 

 

i

 

 

..............................................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

am2 x2 ... amn xn bm;

 

am1 x1

 

x

j

0,

j 1, ..., n.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При постановке конкретной задачи удобно сначала все исходные дан-

ные свести в таблицу (см. пример 4).

 

 

 

 

 

 

Пример 4.

Имеется четыре варианта раскроя j 1, 2, 3, 4 на заготовки двух типов A, B (i 1, 2 ) при заданных потребностях b1 180 , b2 480 и образующихся отходах c1 12, c2 5, c3 3, c4 0 .

Количество заготовок i-го типа при использовании j-го способа произ-

 

4

3

2

1

 

водства aij задается матрицей

 

 

 

 

.

 

0

4

6

10

 

 

 

Для решения задачи минимизации отходов удобно применять следующую таблицу (табл. 5).

 

 

 

 

 

 

Таблица 5

Заготовка

 

Варианты раскроя (j)

 

Потребность

 

(i)

1

2

3

 

4

bi

 

A

4

3

2

 

1

180

 

B

0

4

6

 

10

480

 

Отходы cj

12

5

3

 

0

 

 

20

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