- •Нахождение экстремумов
- •Построение графика
- •Численный метод нахождения экстремума
- •Теоретические вопросы
- •Метод деформируемого многогранника (Нелдера Мида)
- •Метод наискорейшего спуска
- •Метод сопряженных направлений и его модификации
- •Метод Ньютона и его модификации
- •Метод дробления шага
- •Скрипт MATLAB для построения графика
- •Описание функции
- •Текст скрипта
- •Реализация алгоритма Ньютона
- •Тест реализации алгоритма Ньютона
Содержание
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