- •Содержание
- •Предисловие
- •Глава I. Нелинейное программирование: экстремумы функций нескольких переменных
- •§1. Классические методы оптимизации
- •1.1. Безусловный экстремум функции одной переменной
- •1.2. Условный экстремум функции одной переменной
- •1.3. Безусловный экстремум функций двух переменных
- •2.1. Необходимые и достаточные условия
- •1.4. Упражнения.
- •§2. Безусловный экстремум функций нескольких переменных
- •§3. Условный экстремум функций нескольких переменных
- •3.1. Основные определения и факты
- •3.2. Условный экстремум при ограничениях типа равенств
- •3.3. Условный экстремум при ограничениях типа неравенств
- •3.4. Условный экстремум при смешанных ограничениях
- •3.5. Понятие о методах штрафных функций
- •Глава II. Численные методы нелинейного программирования
- •§1. Общие положения
- •1.1. Постановка проблемы
- •1.2. Общие принципы.
- •§2. Методы нулевого порядка одномерной минимизации
- •2.1. Общие положения
- •2.2. Метод равномерного поиска
- •2.3. Метод деления интервала пополам
- •§3. Методы первого и второго порядка
- •3.1. Метод градиентного спуска с постоянным шагом
- •3.2. Метод Ньютона
- •§4. Элементы численных методов задачи условной оптимизации выпуклого программирования
- •4.2. Постановка задачи выпуклого программирования.
- •4.3. Численные методы на основе метода штрафных функций
- •4.4. Метод проекции градиента.
- •Приложение 1. Варианты индивидуальных заданий
- •Задание нп-1
- •Задание нп-2
- •Задание нп-3
- •Задание нп-4
- •Задание нп-5
- •Задание нп-6
- •Задание нп-7
- •Задание нп-8
3.4. Условный экстремум при смешанных ограничениях
В этом случае утверждения, 3.1.1 3.1.4 остаются без изменения.
Таким образом, для того, чтобы исследовать функцию на условный экстремум, достаточно:
Составить обобщённую функцию Лагранжа:
L(X, 0, )=0f(X)+
Составить систему:
(3.4.1)
Решить систему (3.4.1) для двух случаев:
а) 0=0;
б) 0=1.
В результате находятся условно-стационарные точки X*.
Для условно стационарных точек X*, полученных в пункте 3, проверить достаточные условия экстремума. Для этого:
1) определить число l ограничений-равенств и активных ограничений-неравенств в точке X*;
2) если l=n и >0 для всех iIa, то в точке X* локальный минимум. Если l=n и <0 для всех iIa, то в точке X* локальный максимум. Если l<n или соответствующие множители Лагранжа не удовлетворяют достаточным условиям первого порядка, то проверить достаточные условия второго порядка. Для этого:
1) записать выражение для второго дифференциала классической функции Лагранжа в точке (X*, *):
d2L(X*, *)= ;
2) записать условия, накладываемые на первые дифференциалы активных ограничений:
(3.4.2)
3) исследовать знак второго дифференциала функции Лагранжа для ненулевых dx, удовлетворяемых системе (3.4.2). Если d2L(X*, *)>0 (и 0), то в точке X* условный локальный минимум. Если d2L(X*, *)<0 (и 0), то в точке X* условный локальный максимум.
Если достаточные условия первого и второго порядка не выполняются, то проверяется выполнение необходимых условий второго порядка. Если они не выполняются, то в точке X* нет условного экстремума, а если выполняются, то требуется дополнительное исследование.
5. Вычислить значения функции в точках условного экстремума.
Пример I. Исследовать на условный экстремум функцию:
f(X)=
extr
Решение. Имеем
1(Х)=
,
2(Х)=
Далее действуем по вышеприведённой схеме:
1. Составим функцию Лагранжа:
L(X,
0,
)=0(
)+1(
)+2(
).
2. Составим систему (3.4.1). Имеем
=20x1+21x1+22x1=2x1(0+1+2),
=20(x24)+21(x2+4)+22x2.
Поэтому
3. Решим систему для двух случаев:
a) 0=0. Система принимает вид
а.1) x1=0. Тогда из третьего уравнения имеем x2=0, а второе и последнее дают 1=2=0, то есть 0=1=2=0, что невозможно.
а.2) x1≠0. Тогда 1+2=0 и 1=2≠0 и второе уравнение даёт 1(x2+4x2)=0, то есть 41=0, что невозможно. (Здесь тот аргумент, что 1 и 2 имеют разные знаки, не уместен, так как 1 соответствует ограничению-равенству, а 2 ограничению неравенству.)
Таким образом, условно-стационарных точек при 0=0 нет.
б) 0=1. Система принимает вид
б.1)
x1=0.
Тогда из третьего уравнения имеем x2=0,
подставляя которое во второе, получаем
4+41=0,
то есть 1=1.
Последнее уравнение системы даёт 2=0.
Таким образом, (
)=(0,
0; 1, 0)
условно-стационарная точка.
б.2) x1≠0. Тогда 1=12. Предположим, 2≠0. Тогда имеет место система
откуда
.
Решаем последнее уравнение: х2=
.
Подставляя во второе, получаем х1=
.
Подставляя найденное значение х2
во второе уравнение исходной системы,
получаем 712=9,
и приходим к системе
решением которой является 1=1, 2=2 (как и выше, знаки 1 и 2 могут быть различными). Это можно найти и по-другому: преобразовав второе уравнение, получаем
(так
как
),
откуда 1=1, 2=2.
Таким образом,
(
)=(
;
;
1 2),
(
)=(
;
;
1 2)
ещё две условно-стационарные точки.
Если предположить, что 2=0, то 1=1. Тогда второе уравнение системы даёт х24(х2+4)=0, то есть 8=0 противоречие.
Таким образом, условно-стационарных точек три:
( )=(0; 0; 1; 0), ( )=( ; ; 1 2),
( )=( ; ; 1 2),
при этом в точке ( ) имеем 2=0, то есть 20 и 20, в точке выполняются необходимые условия как для минимума, так и для максимума. В точках ( ) и ( ) имеем 2=20, и в этих точках выполняются необходимые условия для максимума.
4. Проверим достаточные условия экстремума в найденных точках.
Имеем
= =2(0+1+2), = =0.
Поэтому d2L=2(0+1+2)(d +d ).
а) В точке X0=(0, 0) ограничение не является активным, так как (X0)=4<0. Поэтому достаточные условия первого порядка не выполняются. Проверим достаточные условия второго порядка. В точке ( ) имеем 0=1=1, 2=0. Поэтому d2L(А)=4(d +d ). Имеем d1=2x1dx1+2(x2+4)dx2, d1(А)=8dx2=0 dx2=0. Поэтому d2L(А)=4d >0 при dx1≠0 (так как ограничение 2(Х)0 пассивно, то Ja= и условие d2(Х)=0 не рассматривается). Следовательно, в точке X1 условный локальный минимум.
б) В точках ( ) и ( ) ограничение-неравенство активно, так как 2(Х*)=0. Поэтому суммарное число активных ограничений-неравенств в точках и и ограничений равенств равно числу переменных. Так как 2=2, то в точках выполняются достаточные условия максимума первого порядка и оно является точкой локального максимума.
5. Вычислим значение функции в точках условного экстремума:
fmin(X)=f(X1)=f(0, 0)=16,
fmax(X)=f(
)=f(
)=
=24.
Ответ:
X1=(0,
0)
точка регулярного условного локального
минимума, fmin(X)=16;
=
и
=
точки регулярного условного локального
максимума, fmax(X)=24.
