- •Классификация
- •Классификация задач оптимизации
- •Классификация методов оптимизации
- •Геометрические и физические свойства производных
- •Критерии оптимальности
- •Линии уровня
- •Построение линий уровня для функции первого порядка
- •Построение линии уровня для функции второго порядка
- •Числовые характеристики симметричной квадратной матрицы
- •Определение квадратной матрицы и её свойства.
- •Практическое задание №1
- •Практическое задание №2
- •Формула Тейлора для функции одной и нескольких переменных
- •Квадратичная аппроксимация функции
- •Нахождение локальных экстремумов
7.Нахождение локальных экстремумов
Дана функция 15. Воспользуемся критериями оптимальности, чтобы найти экстремумы.
1.находим производные функции
fx0 = e x2+xy y2 + (x + y 6)( 2x + y)e x2+xy y2
fy0 = e x2+xy y2 + (x + y 6)(x 2y)e x2+xy y2
2.приравниваем производные к нулю и собираем систему уравнений
(
e x2 |
+xy y2 |
+ (x + y |
|
6)( |
2x + y)e x2+xy y2 |
= 0 |
||||
2 |
+xy y |
2 |
+ (x + y |
|
|
2 |
+xy y |
2 |
= 0 |
|
e x |
|
|
6)(x |
|
2y)e x |
|
||||
|
|
|
|
|
|
|
|
|
3. делим оба уравнения на e x2+xy y2 , т. к. он не равен 0
(
1 |
+ (x + y 6)( 2x + y) |
= 0 |
1 |
+ (x + y 6)(x 2y) |
= 0 |
4.вычитаем из первого уравнения второе
(x + y 6)( 2x + y x + 2y) = 0
5.собираем воедино, умножаем
3x2 + 3y2 + 18x 18y = 0
6.делим всё на 3
x2 y2 6x + 6y = 0
7.собираем из этого квадраты разниц
(x2 (2 3)x + 9) 9 (y2 (2 3)y + 9) + 9 = 0 ) (x 3)2 (y 3)2 = 0
8.из этого следует логический вывод, что x = y:
(x 3)2 = (y 3)2
9.упрощаем уравнение из 4, подставив в него вместо y x
1 + (x + x 6)( 2x + x) = 0 ) 2x2 6x 1 = 0
|
|
|
p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10. |
корни получившегося уравнения: |
3 11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
2 |
|
|
|
|
|
|
|
|
!, |
|
|
|
|
|
|
|
! |
|
|
|
|
3 + p |
|
|
; |
3 + p |
|
3 p |
|
|
; |
3 p |
|
|
|||||
11. |
собственно, стационарные точки: |
11 |
11 |
11 |
|
11 |
|
|||||||||||||
|
2 |
|
|
|
|
2 |
|
2 |
|
2 |
|
18
12.для того, чтобы проверить экстремум, нужно найти гессианы функции в обоих точках fxx00 = ( 2x + y)e x2+xy y2 + ( 4x y + 12)e x2+xy y2 +
+( 2x2 xy + y2 + 12x 6y)( 2x + y)e x2+xy y2 =
=e x2+xy y2 ( 2x + y 4x y + 12+
+4x3 24x2 3xy2 + 24xy + y3 6y2) =
=e x2+xy y2 (4x3 24x2 3xy2 + 24xy 6x + y3 6y2 + 12);
fxy00 = (x 2y)e x2+xy y2 + ( x + 2y 6)e x2+xy y2 +
+( 2x2 xy + y2 + 12x 6y)(x 2y)e x2+xy y2 =
=e x2+xy y2 (x 2y x + 2y 6
2x3 + 3x2y + 12x2 + 3xy2 30xy 2y3 + 12y2) =
=e x2+xy y2 ( 2x3 + 3x2y + 12x2 + 3xy2 30xy 2y3 + 12y2 6);
fyy00 = (x 2y)e x2+xy y2 + ( x 4y + 12)e x2+xy y2 +
+(x2 xy 2y2 6x + 12y)(x 2y)e x2+xy y2 =
=e x2+xy y2 (x 2y x 4y + 12+
+x3 3x2y 6x2 + 24xy + 4y3 24y2) =
=e x2+xy y2 (x3 3x2y 6x2 + 24xy + 4y3 24y2 6y + 12):
13.так как x = y, сократим уравнения для упрощения подсчётов
fxx00 = 2e x2 (x3 3x2 3x + 6); fxy00 = 2e x2 (x3 3x2 3);
fyy00 = 2e x2 (x3 3x2 3x + 6) fxx00 :
14.подставим первую найденную стационарную точку в уравнения
fxx00 |
3 +2p11; |
3 +2p11 |
= 2e ((3+p11)=2)2 |
0 |
|
3 +2p11 |
|
3 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
! |
|
|
|
|
|
|
|
2@ |
|
|
|
|
|
|
|
|
! |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
3 +2p11 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
3 |
3 |
|
3 +2p11 |
+ 61 |
= |
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
! |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
p |
|
|
|
p |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
63 + |
19 11 |
|
|
|
9 |
11 |
9 + 3 |
11 |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
= 2e (10+3p11)=2 |
|
|
|
|
|
|
|
|
|
30 + |
|
|
|
|
|
|
+ 6! = |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
2 |
|
|
|
2 |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
= e (10+3p11)=2 |
|
|
|
|
|
! 1:76505 10 4; |
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
19
fxy00 |
3 +2p11; |
3 +2p11 |
= 2e ((3+p11)=2)2 |
0 |
|
3 +2p11 |
3 |
|
3 |
3 +2p11 |
|
2 |
|
31 |
= |
|||||||||||||||||||||
|
|
|
|
|
|
! |
|
|
|
|
@ |
|
|
|
|
|
|
|
|
! |
|
|
|
|
|
|
|
|
! |
|
|
A |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
= 2e (10+3p11)=2 |
|
|
|
4 |
p |
|
|
|
|
|
|
p |
|
3! |
= |
|
||||||||||||||
|
|
|
|
|
|
63 + |
|
|
|
|
|
|
2 |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
19 |
|
11 |
|
30 + 9 |
11 |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
p |
|
|
|
|
|
|
11 |
|
9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
= e (10+3 11)=2 |
|
|
|
|
|
|
|
1:32287 10 4 |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|||||||||||||||||||
15. формируем, собственно, сам гессиан |
1:32287 |
1:76505 |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
Hf = 10 4 |
|
1:76505 |
1:32287 |
|
|
|
|
|
|
|
|
16. определяем угловые миноры 1 < 0 и 2 > 0 ) этот гессиан отрицательно определён |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
3 + p |
|
; |
3 + p |
|
! является точкой максимума (см. достаточное условие экс- |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
) точка |
11 |
11 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
2 |
|
2 |
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
тремальности). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
17. подставляем вторую точку в уравнения из 13 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
fxx00 |
3 2p11; |
3 2p11 |
= 2e ((3 p11)=2)2 |
|
0 |
|
|
3 2p11 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
! |
|
|
|
|
|
|
|
|
|
|
2@ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
! |
|
|
+ 61 = |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
3 2p11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
3 |
|
3 2p11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
! |
|
|
A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
= 2e (10 3p11)=2 |
|
|
|
63 |
4 |
p |
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
p |
|
|
|
= |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
2 |
|
|
+ 6! |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
19 |
|
|
11 |
|
30 |
|
9 |
11 |
|
9 |
|
|
3 |
11 |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
= e (10 3p11)=2 |
|
9 + |
|
p |
|
|
|
|
|
|
12:475; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
3 2p11; |
3 2p11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
3 2p11 |
|
|
|
|
31 |
|
|
|
||||||||||||||||||||
fxy00 |
= 2e ((3 p11)=2)2 |
|
0 |
|
|
3 2p11 |
3 |
|
3 |
|
|
2 |
|
|
= |
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
! |
|
|
|
|
|
|
|
|
|
|
|
@ |
|
|
|
|
|
|
|
|
|
|
|
|
! |
|
|
|
|
|
|
|
|
|
|
|
|
|
! |
|
|
A |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
= 2e (10 3p11)=2 |
|
|
|
63 |
4 |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
3! |
= |
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
19 |
|
|
11 |
|
|
30 |
|
9 |
11 |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
= e (10+3p11)=2 |
|
|
|
p |
|
|
|
|
! |
|
|
|
2:86683 10 4 |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
11 + 9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
18. формируем, собственно, сам гессиан |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12:475 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
Hf = |
2:86683 10 4 |
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12:475 |
|
|
|
|
|
|
|
|
|
|
|
2:86683 |
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
20
19. определяем угловые миноры 1 > 0 и 2 > 0 ) этот гессиан положительно опреде- |
||
лён ) точка |
3 2p11; |
3 2p11! является точкой минимума (см. достаточное условие |
экстремальности).
Проверка найденного решения. Для проверки верности нахождения экстремумов воспользуемся графиками функции вблизи стационарных точек.
Рис. 7: График функции вблизи точки минимума
Рис. 8: График функции вблизи точки максимума
Из графиков выше видно, что найденные стационарные точки являются экстремумами функции.
Теперь проверим экстремумы численно. Для начала вычислим значение функции в ближайших точках
pp !
3 |
|
11 3 |
|
11 |
||||
f |
|
2 |
|
; |
|
2 |
|
6:1602798 |
21
|
3 p |
|
|
|
|
3 p |
|
|
|
|
|
|
! |
|
|
f |
11 |
+ 10 4; |
11 |
|
+ 10 4 |
6:1602797 |
|||||||||
2 |
|
|
|
2 |
|
|
|||||||||
|
3 p |
|
|
|
|
3 p |
|
|
|
! |
|
||||
f |
11 |
10 4; |
11 |
|
10 4 |
6:1602797 |
|||||||||
2 |
|
|
2 |
|
Вывод. Значение функции в ближайших к экстремуму точках больше, чем значение минимума ) точка минимума найдена верно.
Проверим максимум функции, определив значение функции в ближайших точках
|
|
p |
|
|
|
|
p |
|
|
|
! |
|
|
|
|
|
|
|
3 |
+ 11 3 + |
11 |
|
|
|
|
|
|||||||||
f |
|
|
|
; |
|
|
|
|
|
|
1:4739567 10 5 |
||||||
|
2 |
|
|
2 |
|
|
|
|
|||||||||
|
|
+ p |
|
|
|
3 + p |
|
|
|
|
! 1:4739564 10 5 |
||||||
f |
3 |
11 |
|
+ 10 4; |
11 |
+ 10 4 |
|||||||||||
|
2 |
|
|
2 |
|
|
|||||||||||
|
|
+ p |
|
|
|
|
3 + p |
|
|
! 1:4739564 10 5 |
|||||||
f |
3 |
11 |
|
10 4; |
11 |
10 4 |
|||||||||||
|
2 |
|
|
2 |
|
Вывод. Значение функции в ближайших к экстремуму точках меньше, чем значение максимума ) точка максимума найдена верно.
22