Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лабораторная работа 2. 6 вар. Трофимов-2019.pdf
Скачиваний:
247
Добавлен:
19.02.2020
Размер:
172.35 Кб
Скачать

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