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

Содержание

1.

Нахождение экстремумов

3

2.

Построение графика

6

3.

Численный метод нахождения экстремума

7

4.

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

8

 

4.1.

Метод конфигураций (Хука-Дживса) . . . . . . . . . . . . . . . . . . . . . . . . .

8

 

4.2.

Метод деформируемого многогранника (Нелдера Мида) . . . . . . . . . . . . . .

8

 

4.3.

Метод наискорейшего спуска . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

8

 

4.4.

Метод сопряженных направлений и его модификации . . . . . . . . . . . . . . .

9

 

4.5.

Метод Ньютона и его модификации . . . . . . . . . . . . . . . . . . . . . . . . . .

9

 

4.6.

Метод дробления шага . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

9

A Скрипт MATLAB для построения графика

10

 

1.1.

Описание функции . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

10

 

1.2.

Текст скрипта . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

10

B

Реализация алгоритма Ньютона

11

C

Тест реализации алгоритма Ньютона

21

1.Нахождение экстремумов

Дана

функция

 

x2

x12 2 + 100 (1 x1)2

(1)

Определение оптимума функции.

находим первые частные производные функции

fx01 = 4x31 4x1x2 + 200x1 200 fx02 = 2x2 2x21

приравниваем их к нулю, чтобы найти стационарные точки

(2x2

 

2x12

+ 200x1 200

= 0

4x13

4x1x2

= 0;

 

 

 

 

 

• сокращаем первое уравнение на 4, второе на 2

(x2

x12

= 0

x13

x1x2

+ 50x1 50 = 0;

• из второго уравнения видно, что x2 = x21, заменяем в первом уравнении x31 x31 + 50x1 50 = 0;

50x1 = 50; x1 = 1; x2 = 1

стационарная точка найдена: (1; 1);

находим Гессиан функции

fx001x1 = 12x21 4x2 + 200; fx001x2 = 4x1;

fx002x2 = 2; fx001x1 (1; 1) = 208;

fx001x2 (1; 1) = 4; fx002x2 (1; 1) = 2;

 

 

 

Hf(1; 1) =

208

4

 

4

2

определяем знак угловых миноров: 1 > 0, 2 > 0 ) гессиан положительно определён ) стационарная точка является точкой минимума;

значение функции в точке минимума: f(1; 1) = 0.

3

Определение вогнутости и выпуклости функции.

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

@x12 > 0;

@x22 > 0;

@x12 @x22

x1x2!

2

> 0

@2f

 

@2f

 

@2f @2f

@2f

 

 

 

 

 

 

 

 

 

 

 

 

• составляем систему неравенств и решаем её

(2 (12x12

 

 

4x2

+ 200)

 

16x12 > 0

12x12 4x2 + 200

 

> 0;

(3x12

 

 

 

 

 

 

 

 

x2

+ 50

 

2x12

> 0

3x12

x2

+ 50

 

 

> 0;

(x2

 

 

 

 

 

 

 

 

< x12

+ 50

 

 

 

x2

< 3x12 + 50;

 

 

 

очевидно, выпуклое вниз подмножество функции: x2 < x21 + 50;

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

@x12 < 0;

@x22 < 0;

@x12 @x22

x1x2!

2

> 0

@2f

 

@2f

 

@2f @2f

@2f

 

 

 

 

 

 

 

 

 

 

 

 

• составляем систему неравенств и решаем её

(2 (12x12

 

 

4x2

+ 200)

 

16x12 > 0

12x12 4x2 + 200

 

< 0;

(3x12

 

 

 

 

 

 

 

 

x2

+ 50

 

2x12

> 0

3x12

x2

+ 50

 

 

< 0;

(x2

 

 

 

 

 

 

 

 

> x12

+ 50

 

 

 

x2

< 3x12 + 50;

 

 

 

• очевидно, функция не имеет выпуклого вверх подмножества.

4

Определение овражности.

• Находим собственные числа гессиана в точке 1 1

Hu = u;

 

 

 

u (H E) = 0;

 

H E = 0;

0

 

4

2

 

208

4

 

 

0 = 0;

208

4

= 0;

4

2

 

(208 ) (2 ) = 0;

2 210 + 416 = 0; p

210 2102 4 416

1;2

=

 

 

 

 

 

;

 

 

2

 

 

 

 

 

 

 

 

 

1;2

= 105 103

 

 

 

• находим число обусловленности в точке 1 1

 

 

max

 

208

 

cond H =

 

=

 

= 104 > 10

min

2

• очевидно, cond H > 10 функция в окрестности точки 1 1 является овражной.

5

2.Построение графика

Для построения графика функции был выбран математический пакет MATLAB. Получившийся результат на рисунке 1. Красной точкой на рисунке обозначен минимум функции 1 1 . В приложении A описан скрипт для генерации такого графика.

Рис. 1: График функции 1 с точкой минимума

6