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

Исследование операций и методы оптимизации в примерах. Практикум

.pdf
Скачиваний:
0
Добавлен:
09.08.2026
Размер:
470 Кб
Скачать

МИНИСТЕРСТВО НАУКИ И ВЫСШЕГО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ «ДОНСКОЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»

Кафедра «Вычислительные системы и информационная безопасность»

ИССЛЕДОВАНИЕ ОПЕРАЦИЙ И МЕТОДЫ ОПТИМИЗАЦИИ В ПРИМЕРАХ

ПРАКТИКУМ

Ростов-на-Дону ДГТУ

2019

УДК 681. 5

Составители: А.Ю. Полуян, С.Б. Петренкова

Исследование операций и методы оптимизации в примерах: практикум. – Ростов-на-Дону : Донской гос. техн.

ун-т, 2019. – 28 с.

Описаны методы решения некоторых задач с использованием оптимизационных методов.

Предназначены для обучающихся направлений подготовки 10.03.01 Информационная безопасность, 09.03.02 Информационные системы и технологии, 10.05.02 Информационная безопасность телекоммуникационных систем.

УДК 681. 5

Печатается по решению редакционно-издательского совета Донского государственного технического университета

Научный редактор канд. техн. наук, доцент А.И. Зотов

Ответственный за выпуск зав. кафедрой «Вычислительные системы и информационная безопасность» д-р техн. наук, профессор В.А. Фатхи

_____________________________________________________

В печать 08.02.2019 г.

Формат 60×84/16. Объем 0,8 усл. п. л. Тираж 50 экз. Заказ № 152.

__________________________________________________

Издательский центр ДГТУ Адрес университета и полиграфического предприятия:

344000, г. Ростов-на-Дону, пл. Гагарина, 1

© Донской государственный технический университет, 2019

2

1. БЕЗУСЛОВНЫЙ ЭКСТРЕМУМ ФУНКЦИИ ДВУХ ПЕРЕМЕННЫХ

Цель работы: Исследование стационарных точек функции двух переменных, определяемых аналитическим методом.

Задание: Для функции f(x,y) и уравнения связи g(x,y)=0 найти аналитическим методом:

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

2)определить вид каждой из стационарных точек (максимум, минимум, нет экстремума), воспользвавшись достаточными условиями экстремума.

При несоблюдении требований правила проверки достаточного условия (B2 – АС = 0, разные знаки у А и С и т.д.) указать, по какой причине заключение

охарактере соответствующей стационарной точки сделать нельзя.

При выполнении преобразований делать максимальные упрощения с использованием правильных дробей, а все конечные численные значения координат и значений функций указывать в десятичной форме с округлением дробной части до трёх знаков (например, вместо 5/3 записать 1,66(6), а вместо 22 значение 2,628);

3) Результаты выполнения задания оформить в виде отчёта о лабораторной работе.

ПРИМЕР:

f(x,y)=2x3+xy2+5x2+y2

Определение координат стационарных точек

 

 

 

 

f

2

2

 

 

 

 

 

 

 

6x

y 10x 0,

(1)

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

2xy 2y 0.

(2)

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

Так как из (2) следует, что 2y(x 1) 0, то

 

 

 

 

 

 

 

 

либо x1 1, либо y1

0.

Подставляя оба эти значения в (1), получаем

 

6x2 y2 10x

 

 

6 ( 1)2 y2 10( 1) y2 4 0,

 

 

 

x1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

y2 2;

y3 2;

 

6x2 y2 10 x

 

 

6x2 10x x(6x 10) 0,

 

 

 

y1 0

 

5

 

 

 

 

 

 

 

 

 

 

при этом x2 0,

x3

1,667.

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

Проверка:

6x2 10 x

 

 

x2 0

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6x

2

10 x

 

 

5

6(

5

)

2

10(

5

)

50

 

50

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

3

 

3

3

3

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, имеем 4 стационарные точки (результаты указаны в

десятичной форме с округлением дробной части до трёх знаков):

(x1,y2) ( 1,2),

 

(x1,y3) ( 1, 2),

 

 

 

 

 

 

(x2,y1 ) (0,0),

 

(x3 ,y1 ) ( 1,667;0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определение вида стационарных точек

A

2f

 

(6x2 y2 10x)

12x 10,

 

 

 

 

 

 

x2

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B (

 

 

2 f

 

 

)

(6x2 y2

 

 

10 x )

2y,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

C (

 

2 f

 

)

(2xy 2y)

2x 2;

 

 

 

 

 

 

 

y2

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A 12x 10 10, B 2y 0, C 2x 2 2. Здесь

 

 

 

 

в точке (0,0):

 

 

B2 AC 0 10 2 20 0; значит (0,0 ) является точкой экстремума, а так как А и С положительны, то точкой минимума; значение f(0,0) = 2x3+xy2+5x2+y2

= 0; в точке (–1,667;0):

A 12x 10 12 ( 5) 10 10, B 2y 0,

3

C 2x 2 2 ( 5 ) 2 1,333.

3

Так как

B2 AC 0 ( 10) (1,333) 13,33 0,

то рассматриваемая точка является точкой экстремума, а так как А и С отрицательны, то точкой максимума;

Значение функции:

f( 1,67;0)=2x

3

2

 

2

 

2

 

5

 

3

 

5

 

 

2

 

 

+xy

+5x

 

+y

 

=2(

 

)

 

+5(

 

 

)

 

= 4,630

 

 

 

3

 

3

 

 

в точке (-1,2):

 

A 12x 10 2,

 

B 2y 4,

C 2x 2 0. Так как

B2 AC 16 ( 2) 0 16 0, в рассматриваемой точке нет экстремума;

значение f(1,2) = 2x3+xy2+5x2+y2 = 3;

 

 

 

 

 

 

в точке (1,2): A 12x 10 2,

B 2y 4,

 

 

C 2x 2 2 2 0. Так

как B2 AC 16 ( 2) 0 16 0, в рассматриваемой точке нет экстремума;

значение f(1, 2) = 2x3+xy2+5x2+y2 = 3.

4

2. УСЛОВНЫЙ ЭКСТРЕМУМ ФУНКЦИИ ДВУХ ПЕРЕМЕННЫХ

Цель работы: исследование особенностей определения условного экстремума функции двух переменных методом множителей Лагранжа.

Задание: для функции f(x,y) и уравнения связи g(x,y)=0:

1) найти аналитически методом множителей Лагранжа координаты и вид условных экстремумов функции f(x,y) для указанного уравнения связи g(x,y)=0, а также значения функции в точках условного экстремума.

При выполнении преобразований делать максимальные упрощения с использованием правильных дробей, доводя решаемое уравнение до записи с целочисленными коэффициентами. Конечные численные значения указать в десятичной форме с округлением дробной части до трёх знаков (например, вместо 5/3 записать 1,66(6), а вместо 22 значение 2,628).

2) результаты выполнения задания оформить в виде отчёта о лабораторной работе (без приводимых далее дополнительных текстовых пояснений).

ПРИМЕР:

f (x,y)=2x3+xy2+5x2+y2 ; g(x,y)=3x+4y – 1=0.

1. Составляем вспомогательную функцию:

Ф(х,у)= f(x,y)+λ g(x,y)= 2x3+xy2+5x2+y2+λ(3x+4y1).

2. Приравниваем частные производные Ф (х, у) по x и по y нулю:

 

Ф

6x2

y2

10x 3 0;

Ф

2xy 2y 4 0.

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

2.1. Исключаем λ:

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как

(2xy 2y)

, получаем уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

6x2 y2

10x 3

(2xy 2y)

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Полученное уравнение совместно с уравнением связи

 

 

 

 

Получаем из уравнения связи

g(x,y)=3x+4y–1=0 значение x

1 4y

и

 

 

 

 

подставляем в полученное уравнение:

 

 

 

3

 

 

 

 

 

 

 

 

 

1 4y

 

 

 

 

 

 

1 4y

 

 

(2

1 4y

y 2y)

 

 

 

 

 

2

 

2

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6(

 

 

)

 

y

 

10(

 

) 3

 

 

 

 

 

0.

 

3

 

 

3

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

Сделать максимальные упрощения с использованием, для сохранения точности, правильных дробей и получением квадратного уравнения с целочисленными коэффициентами:

6

 

 

 

 

 

 

 

2

 

 

 

2

 

10 40y 3

 

 

 

 

 

 

 

 

2y 8y2

 

 

 

 

 

 

 

(1

 

4y)

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

 

 

 

 

2y) 0,

9

 

 

 

 

3

 

 

 

 

 

 

3

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(×6)

 

2

(1 8y

16y

2

) y

2

 

 

10

 

 

 

 

40y

 

 

y

2y

2

 

3y

0,

3

 

 

3

 

 

 

 

3

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 32y 64y2 6y2 20 80y 3y 12y2 9y 0,

 

82y2 124y 24 0 или

 

 

 

41y2

62y 12 0.

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y4,5

 

 

62

 

( 62)2

4 41 12

 

62 43,313

.

 

 

 

 

 

 

 

 

 

 

 

 

 

2 41

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

82

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y4

105,313

1,284; x4

 

1 4y

 

1 4 1,284

1,379;

 

 

 

 

 

 

 

 

 

 

 

 

 

82

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

y5

18,69

0,228;

 

 

x5

1 4y

 

 

 

1 4 0,228

0,029.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

82

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

4. Определяем вид экстремумов

Рассматриваемая функция f(x,y)=2x3+xy2+5x2+y2 третьего порядка, то есть имеет один минимум и один максимум. С учётом этого обстоятельства можно избежать аналитических преобразований и определить вид экстремумов в точках условного экстремума, сравнив значения функции в этих точках:

f(x4,y4) = 2 ( 1,379)3+( 1,379) 1,2842+5 ( 1,379)2+1,2842= 3,638; f(x5,y5) =2 0,0293 +0,029 0,2882+5 0,0292+0,2882= 0,058.

Таким образом:

вточке (x4,y4) максимум,

вточке (x5,y5) минимум.

6

3. ЭЛЕМЕНТАРНАЯ ЗАДАЧА ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ Цель работы: закрепление навыков аналитического решения задач

оптимизации с использованием метода множителей Лагранжа

Задание: для заданного функционала при указанных граничных условиях: 1) Найти экстремаль (аналитическое выражение) с проверкой частного решения дифференциального уравнения, полученного при указанных

граничных условиях.

Примечания: результаты решения линейного однородного дифференциального уравнения (в частности, выражения с корнями), а также все окончательные результаты представлять десятичными числами с точностью до трёх знаков; проверить решения дифференциального уравнения (отклонение от 0 должно быть не более 0,01 – менее 1%).

2) Результаты выполнения задания оформить в виде отчёта о лабораторной работе.

ПРИМЕР:

2

y 1 1, y 1 2 .

( x2 y2 y'2 )dx ,

1

 

 

 

 

Алгоритм выполнения задания:

 

dFy

 

1. Cоставим уравнение Эйлера (Fy

 

0 )

 

 

 

 

dx

1.1. Исходя из F = x2 + y2 + y 2, получим Fy и Fy

Fy (x2 y2 y 2 )y 2y,

Fy (x2 y2 y 2 )y 2y .

1.2. Запишем выражение для полной производной dFy , учитывая, что dx

F 2y :

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

dFy (x,y,y )

 

 

Fy

 

Fy

 

 

Fy

y

 

 

dx

 

 

 

x

 

y

y

 

 

 

y

 

0 0 2y

 

 

 

 

 

 

 

 

 

 

 

 

 

2y .

 

 

 

 

 

 

 

 

7

1.3. Подставим

полученные

выражения

для Fy 2y

и

dFy

2y в

dx

 

 

dFy

 

 

 

 

 

 

уравнение Эйлера (Fy

 

 

0 ) :

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

2y 2y 0, или после упрощения

y y 0.

 

 

 

 

2. Решить уравнение Эйлера.

 

 

 

 

 

2.1. Анализ уравнения.

 

 

 

 

 

 

Полученное уравнение

y y 0 является

линейным

дифференциальным

второго порядка, однородным (без правой части), с постоянными коэффициентами

– [4], с. 589.

Решение такого уравнения ищется с использованием корней характеристического уравнения.

В общем случае для дифференциального уравнения y by cy = 0

характеристическое уравнение имеет вид r2 br c 0

([4], с.с. 587–588), а

общее решение дифференциального уравнения у = С1 er1x + С2 er2x.

В данном случае дифференциальное уравнение

характеристическое

уравнение корни характеристического уравнения общее решение дифференциального уравнения y y 0:

r2 1=0;

r1 = 1; r2 = –1,

у = С1 ex + С2 e x .

2.3. Проверка подстановкой С1 и С2 в (2) и (3):

С1·e 1 + С2·e =0,701/2,718+0,273·2,718 =0,258+0,742=1,

С1·e+С2·e 1 =0,701·2,718+0,273/2,718=1,905+0,1= 2.005,

отклонение менее 1 %.

2.4.Аналитическое выражение, описывающее экстремаль:

у= С1 ex + С2 e x = 0,701ex + 0,273e x .

8

4.ЭКСТРЕМАЛЬФУНКЦИОНАЛА,ЗАВИСЯЩЕГООТ ДВУХФУНКЦИЙ Цель работы: закрепление навыков решения задачи нахождения

экстремалей функционала, зависящего от двух функций.

Задание: для нахождения экстремалей заданного функционала, зависящего от двух функций, при указанных граничных условиях:

1)составить систему уравнений Эйлера.

2)в порядке решения системы дифференциальных уравнений Эйлера получить систему линейных уравнений.

ПРИМЕР

J y,z

2

2

 

2

 

y 2 1;

z 2 0;

(3y'

9z'

7 yz)dx;

 

 

 

0;

 

2.

 

2

 

 

 

 

y 2

z 2

 

 

 

 

 

 

 

 

 

Алгоритм выполнения задания:

1.Составим систему уравнений Эйлера при заданных исходных данных.

1.1.Общий вид системы уравнений

 

 

 

dFy

F

y

 

 

0,

 

 

 

dx

 

 

 

dFz

0.

1.2.НаходимчастныепроизводныеFy,Fy ,Fy и Fz :Fz dx

Fy (3y'2 9z'2 7 yz)y 7z,

Fy (3y'2 9z'2 7 yz)y 6y ,

Fz (3y'2 9z'2 7 yz )z 7 y,

Fz (3y'2 9z'2 7 yz )z 18z .

1.3 Вычислим полные производные

1.4.1. Вспомним выражения для полных производных:

dFy ( x,y,y )

 

Fy

 

Fy

y

Fy

 

y ,

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

y

 

dx

 

 

 

 

 

 

 

 

 

dFz ( x,z,z )

Fz

 

Fz

z

Fz

z

.

 

 

 

 

 

 

 

 

 

 

 

x

 

z

z

dx

 

 

 

 

 

 

 

 

 

1.4.2. Вычислим частные производные (с учётом того, что, как уже показано в п.1.1, Fy 6y ,Fz 18z ):

9

F

 

 

 

F

 

 

 

F

 

y

0,

 

y

0

 

y

6,

x

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

y

 

F

 

F

 

 

 

 

F

 

z

 

 

0,

 

z

0,

 

z

18.

x

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

z

 

1.4.3. Значения полных производных:

dFy

dF

 

 

 

 

6y ,

z

 

18z

.

 

dx

 

dx

 

 

 

1.5. Подставим все полученные значения в систему уравнений Эйлера общего вида, при этом получим искомую систему уравнение Эйлера:

 

7z 6y 0,

 

 

 

 

7 y 18z 0.

 

 

 

 

 

 

 

 

2. Решение системы уравнений Эйлера

 

 

 

 

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

7z 6y 0,

 

 

 

Из первого уравнения полученной системы

 

 

следует, что

 

 

7 y 18z

0,

 

 

 

z 76 y . После подстановки во второе уравнение системы получаем

7 y 18 d2 (6 y ) 0 , dx2 7

или после упрощений y( IV ) 10849 y 0.

Выполняя аналогичные преобразования относительно другой функции,

получаем

y

18

z , а после подстановки 7z 6

d2

(

18

z ) 0,

7

dx2

7

 

 

 

 

 

и z( IV ) 49 z 0.

108

Таким образом, приходим к эквивалентной системе уравнений

y( IV )

z( IV )

10849 y 0,

49 z 0.

108

10

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