WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Pages:     | 1 |   ...   | 4 | 5 || 7 | 8 |

В результате решения задачи безусловной минимизации получаем x x(3) = (3,00; - 2,00), f f (x(3)) = -2,50.

Ответ: x (3,00; - 2,00), f -2,50.

Задачи 1. Найти матрицу A-1, обратную матрице A вида 2 1 A = 4 3 1.

- 2 3 2. Решить методом Ньютона задачу безусловной минимизации 2 f (x) = x1 + 2x2 - 2x1 + x2 - 5 min при = 0,4, x(0) = (0, 2).

3. Решить методом Ньютона задачу безусловной минимизации 3 f (x) = x1 + x2 - x1x2 + x1 - 3x2 - 4 min при = 0,2, x(0) = (2, 2).

8. МЕТОД АППРОКСИМИРУЮЩЕГО ПРОГРАММИРОВАНИЯ Метод аппроксимирующего программирования (МАП) относится к численным методам решения задач условной оптимизации. Наиболее простой и хорошо изученной задачей условной оптимизации является задача линейного программирования (ЛП).

МАП является одним из методов решения задач нелинейного программирования. В данном случае исходная задача нелинейного программирования преобразуется в последовательность задач ЛП с помощью процедур линеаризации. Рассмотрим задачу условной минимизации вида f (x) min, n x X = {x R+ : gi (x) 0, i = 1, r}, где f (x), gi (x) - произвольные нелинейные функции.

МАП является итерационным методом. На k-й, k = 1,2,..., итерации определяется точка x(k ) - k-е приближение к точке минимума x ; при этом исходной точкой для k-й итерации является точка x(k-1). В результате, исходя из заданной начальной точки x(0), находится последовательность точек x(1), x(2),..., сходящаяся при определенных условиях к решению x* исходной задачи нелинейного программирования. Отметим, что в качестве начальной точки x(0) выбирается некоторая точка из допустимого множества X, т.е. x(0) X ; точка x(0) выбирается произвольно, принадлежность к X определяется проверкой выполнения ограничений для этой точки.

На k-й, k=1,2,…, итерации в окрестности точки x(k-1) осуществляется линейная аппроксимация (линеаризация) задачи нелинейного программирования, т.е. каждая из нелинейных функций исходной задачи заменяется двумя первыми членами в разложении в ряд Тейлора. В результате получается задача ЛП ~ f (x) = f (x(k -1) ) + f (x(k -1)),(x - x(k -1)) min, ~ gi (x) = gi (x(k -1) ) + gi (x(k -1) ),(x - x(k -1) ) 0, i = 1, r, n x R+.

Затем находится решение x0 задачи ЛП, после чего определяется точка x(k ) по известным x(k-1) и x0. Эта точка должна удовлетворять условиям x(k ) X, (8.1) < f (x(k )) f (x(k -1)).

Существует много способов определения x(k ). На практическом занятии x(k ) определяется из соотношения x(k ) = x(k -1) + k (x0 - x(k -1)), 0 k 1, k = 1,2,..., где k - параметр (скаляр), определяющий длину шага из точки x(k-1) в направлении точки x0.

Очевидно, что при k = 1 выполняется x(k ) = x0, при k = 0 x(k ) = x(k -1).

Величина k выбирается так, чтобы выполнялись условия (8.1). Процесс выбора шага, удовлетворяющего данным условиям, во многом аналогичен соответствующему процессу в случае градиентного метода с дроблением шага. Выбирается константа 0 < < 1 (в данном случае =1). На k-й, k = 1,2,..., итерации проверяется выполнение условий (8.1) при k=1, т.е. для x(k) = x0. Если они не выполняются, то производится дробление шага, т.е. полагается k =, и вновь проверяется выполнение условий (8.1). Процесс дробления, т.е. умножения текущего значения k на, продолжается до тех пор, пока условия (8.1) не окажутся выполненными.

Алгоритм решения задачи условной минимизации методом аппроксимирующего программирования заключается в следующем.

1. Задаются, 1, ; выбирается x(0) ; полагается k = 1.

2. Осуществляется линеаризация исходной задачи в окрестности точки x(k -1). В результате получается задача ЛП.

3. Находится решение x0 задачи ЛП.

4. Полагается k = 1.

5. Вычисляется x(k ) = x(k -1) + k (x0 - x(k -1)).

6. Проверяются условия выбора x(k) :

gi (x(k )) 0, i = 1, r;

f (x(k )) < f (x(k -1)).

Если они выполняются, то осуществляется переход к п.7.

Если условия не выполняются, то полагается k = k и осуществляется переход к п.5.

7. Проверяются условия окончания решения исходной задачи:

f (x(k )) - f (x(k -1)) 1, f (x(k -1)) x(jk ) - x(jk -1), j = 1, n.

x(jk -1) Если они выполняются, то полагается x* x(k ), * f f (x(k )) и вычисления завершаются.

Если условия не выполняются, то полагается k = k +1 и осуществляется переход к п.2.

Пример. Решить методом аппроксимирующего программирования задачу условной минимизации f (x) = 4x1 - x2 -12 min, 2 x1 + x2 25, 2 10x1 - x1 +10x2 - x2 34, x1 0, x2 при = 0,7, 1 = 0,1, = 0,3.

Решение.

Преобразуем ограничения исходной задачи к виду gi (x) 0 :

2 g1(x) = x1 + x2 - 25 0, 2 g2(x) = x1 -10x1 + x2 -10x2 + 34 0.

Находим f (x), g1(x), g2 (x) :

f f = 4, = -2x2 f (x) = (4,-2x2 );

x1 xg1 g = 2x1, = 2x2 g1(x) = (2x1,2x2 );

x1 xg2 g = 2x1 -10, = 2x2 -10 g2 (x) = (2x1 -10,2x2 -10).

x1 xВыберем x(0) = (2, 4).

Проверяем принадлежность точки x(0) к допустимой области X:

g1(x(0) ) = 22 + 42 - 25 = -5 < 0, g2 (x(0) ) = 22 -10 2 + 42 -10 4 + 34 = -6 < 0, ( ( x10) > 0, x20) > 0.

Поскольку ограничения выполняются, то точка x(0) = (2, 4) является допустимой, т.е. x(0) X.

Первый этап (первая итерация) Осуществляем линеаризацию исходной задачи в окрестности точки x(0) :

f (x(0)) = 4 2 - 42 -12 = 8 -16 -12 = -20, f (x(0) ) = (4, - 2 4) = (4, - 8), ~ f (x) = -20 + (4, -8), (x1 - 2, x2 - 4) = = -20 + 4(x1 - 2) - 8(x2 - 4) = 4x1 - 8x2 + 4;

g1(x(0)) = (2 2, 2 4) = (4, 8), ~ g1(x) = -5 + (4, 8), (x1 - 2, x2 - 4) = = -5 + 4(x1 - 2) + 8(x2 - 4) = 4x1 + 8x2 - 45;

g2(x(0)) = (2 2 -10, 2 4 -10) = (-6, - 2), ~ g2(x) = -6 + (-6, -2), (x1 - 2, x2 - 4) = = -6 - 6(x1 - 2) - 2(x2 - 4) = -6x1 - 2x2 +14.

Составляем задачу ЛП:

~ f (x) = min, ~ g1(x) 0;

~ g2 (x) 0;

x1 0, x2 0.

~ ~ ~ Подставляем f (x), g1(x), g2 (x) :

~ f (x) = 4x1 - 8x2 + 4 min, 4x1 + 8x2 45, (1) 6x1 + 2x2 14, (2) x1 0, x2 0.

Решаем задачу ЛП графическим методом (рис.8.1):

4x1 + 8x2 = 45 : x1 = 0 x2 = 5,625, x2 = 0 x1 = 11,25;

6x1 + 2x2 = 14 : x1 = 0 x2 = 7, x2 = 0 x1 = 2,33;

f (x) = (4,-8), ~ где f (x) - градиент целевой функции задачи ЛП.

Из рис. 8.1 следует, что задача ЛП имеет решение x0.

Точка x является решением системы уравнений 4x1 + 8x2 = 45, + 2x2 = 14.

6xxxx1 2 3 4 5 6 7 8 9 10 11 (1) (2) f x ( ) Рис. 8.Находим x0 :

4x1 + 8x2 = - x1 = = 0,55;

24x1 + 8x2 = 56 - 20x1 = -4 0,55 + 8x2 = 45 8x2 = 45 - 2,2 = 42,8 x2 = 5,35;

x0 = (0,55; 5,35).

Полагаем 1 = 1. Вычисляем x(1):

x(1) = x(0) + 1(x0 - x(0) ) = x0 = (0,55; 5,35).

Проверяем условия выбора x(1) :

g1(x(1) ) = 0,552 + 5,352 - 25 = 3,9 > 0.

Поскольку условия не выполняются, то полагаем 1 = 1 = 0,7. Вычисляем x(1):

x(1) = (2, 4) + 0,7(0,55 - 2; 5,35 - 4) =.

= (2, 4) + (-1,015; 0,945) = (0,985; 4,945).

Проверяем условия выбора x(1) :

g1(x(1) ) = 0,9852 + 4,9452 - 25 = 0,423 > 0.

Поскольку условия не выполняются, то полагаем 1 = 1 = 0,49. Вычисляем x(1):

x(1) = (2, 4) + 0,49(-1,45; 1,35) = (2, 4) + (-0,71; 0,66) = (1,29; 4,66).

Проверяем условия выбора x(1) :

g1(x(1) ) = 1,292 + 4,662 - 25 = -1,62 < 0, g2 (x(1) ) = 1,292 -10 1,29 + 4,662 -10 4,66 + 34 = -2,12 < 0, ( ( x11) > 0, x21) > 0, f (x(1) ) = 4 1,29 - 4,662 -12 = -28,56 < f (x(0) ) = -20.

Поскольку условия выполняются, то x(1) = (1,29; 4,66).

Проверяем условия окончания решения исходной задачи f (x(1) ) - f (x(0) ) - 28,6 + = = 0,428 > 1 = 0,1.

- f (x(0) ) Поскольку условия не выполняются, то выполняем второй этап.

Второй этап (вторая итерация) Осуществляем линеаризацию исходной задачи в окрестности точки x(1) :

f (x(1) ) = (4; -2 4,66) = (4; -9,32), ~ f (x) = -28,56 + (4; -9,32),(x1 -1,29; x2 - 4,66) = = -28,56 + 4(x1 -1,29) - 9,32(x2 - 4,66) == 4x1 - 9,32x2 + 9,71;

g1(x(1) ) = (2 1,29; 2 4,66) = (2,58; 9,32), ~ g1(x) = -1,62 + (2,58; 9,32),(x1 -1,29; x2 - 4,66) = = -1,62 + 2,58(x1 -1,29) + 9,32(x2 - 4,66) = 2,58x1 + 9,32x2 - 48,4;

g2(x(1)) = (21,29 -10; 2 4,66 -10) = (-7,42; - 0,68), ~ g2(x) = -2,12 + (-7,42; -0,68),(x1 -1,29; x2 - 4,66) = = -7,42x1 - 0,68x2 +10,62.

Составляем задачу ЛП:

~ f (x) = min, ~ g1(x) 0;

~ g2 (x) 0;

x1 0, x2 0.

~ ~ ~ Подставляем f (x), g1(x), g2(x) :

~ f (x) = 4x1 - 9,32x2 + 9,71 min, 2,58x1 + 9,32x2 48,4, (1) 7,42x1 + 0,68x2 10,6, (2) x1 0, x2 0.

Решаем задачу ЛП графическим методом (рис. 8.2):

xxx2 4 6 8 10 12 14 16 (1) (2) f x ( ) Рис. 8.2,58x1 + 9,32x2 = 48,4 : x1 = 0 x2 = 5,2, x2 = 0 x1 = 18,8;

7,42x1 + 0,68x2 = 10,6 : x1 = 0 x2 = 15,6, x2 = 0 x1 = 1,43;

~ f (x) = (4; -9,32).

Из рис. 8.2 следует, что задача ЛП имеет решение x0.

Точка x0 является решением системы уравнений 2,58x1 + 9,32x2 = 48,4, + 0,68x2 = 10,6.

7,42xНаходим x0 :

2,58x1 + 9,32x2 = 48,44,- x2 = = 4,92;

2,58x1 + 0,236x2 = 3,7 9,9,084x2 = 44,2,58x1 + 9,32 4,92 = 48,4 2,58x1 = 2,546 x1 = 0,987;

x0 = (0,987; 4,92).

Полагаем 2 = 1. Вычисляем x(2):

x(2) = x(1) + 2(x0 - x(1)) = x0 = (0,987; 4,92).

Проверяем условия выбора x(2) :

g1(x(2)) = 0,9872 + 4,922 - 25 = 0,181 > 0.

Поскольку условия не выполняются, то полагаем 2 = 2 = 0,7. Вычисляем x(2):

x(2) = (1,29; 4,66) + 0,7(0,987 -1,2; 4,92 - 4,66) = = (1,29; 4,66) (-0,21; 0,18) (1,08; 4,84).

+ = Проверяем условия выбора x(2) :

g1(x(2)) = 1,082 + 4,842 - 25 = -0,41 < 0, g2(x(2)) = 1,082 -10 1,08 + 4,842 -10 4,84 + 34 = -0,61 < 0, ( ( x12) > 0, x22) > 0, f (x(2)) = 4 1,08 - 4,842 -12 = -31,11 < f (x(1)) = -28,56.

Поскольку условия выполняются, то x(2) = (1,08; 4,84).

Проверяем условия окончания решения исходной задачи f (x(2)) - f (x(1)) - 31,11+ 28,= = 0,089 < 1 = 0,1, - 28,f (x(1)) ( ( x12) - x11) 1,08 -1,= = 0,163 < 2 = 0,3, ( 1,x11) ( ( x22) - x21) 4,84 - 4,= = < 2 = 0,039 0,3.

( 4,x21) Поскольку условия выполняются, то полагаем * x* x(2) = (1,08; 4,84), f f (x(2)) = -31,1 и вычисления завершаются.

* Ответ: x* (1,08; 4,84), f -31,1.

Задачи 1. Дана задача условной минимизации 2 f (x)= x1 + x2 -16x1 -10x2 min, x1 - 6x1 + 4x2 11, exp(x1 - 3)- x1x2 + 3x2 1, x1 0, x2 0.

Определить методом аппроксимирующего программирования точку x(1) при x(0) = (4, 3), = 0,25.

2. Решить методом аппроксимирующего программирования задачу условной максимизации f (x) = x1 + x2 max, 2x1 - x2 1, 0,8x1 + 2x2 9, x1 0, x2 при= 0,7, 1 = 0,1,2 = 0,2, x(0) = (3, 0).

9. МЕТОД ШТРАФНЫХ ФУНКЦИЙ Метод штрафных функций относится к численным методам решения задач условной оптимизации. В данном случае исходная задача условной оптимизации преобразуется в последовательность задач безусловной оптимизации путем введения штрафных функций. Рассмотрим задачу условной минимизации вида f (x) min, x X = {x Rn : gi (x) 0, i = 1, r}.

На ее основе строится задача безусловной минимизации P(x, R) = f (x) + (R, g(x)) min, x Rn, где P(x, R) - расширенная функция, ( R, g ( x )) - штрафная функция, R - штрафной параметр.

Задача условной минимизации f (x) заменяется последовательностью задач безусловной минимизации P(x, Rt -1) при t = 1,2,.... При этом, исходя из заданной начальной точки x[0], находится последовательность точек x[1], x[2],..., сходящаяся при x определенных условиях к решению исходной задачи. При минимизации расширенной функции P(x, Rt -1), t = 1,2,..., исходной (начальной) точкой является x[t-1], а решение задачи безусловной минимизации P(x, Rt -1) определяет точку x[t].

Методы штрафных функций разделяются на методы внутренней точки и методы внешней точки. Метод штрафных функций называется методом внутренней точки (внешней точки), если все точки последовательности x[t ], t = 0,1,2,..., являются допустимыми (недопустимыми). Вид метода (внутренней или внешней точки) определяет вид штрафной функции и правило, по которому производится пересчет штрафного параметра после решения очередной задачи безусловной минимизации.

Для методов внутренней точки штрафные функции должны обладать следующими свойствами:

- на большей части допустимого множества X внутренние штрафные функции близки к нулю;

- при приближении изнутри к границе допустимого множества X внутренние штрафные функции достаточно быстро возрастают.

В качестве внутренней штрафной функции часто используются логарифмическая штрафная функция r (Rt, g(x)) = -Rt ln(-gi (x)), i=а также обратная штрафная функция r (Rt, g(x)) = -Rt gi (x) i =Внутренние штрафные функции имеют смысл только внутри допустимого множества X, в связи с этим необходимо проверять соблюдение ограничений при решении задач безусловной оптимизации.

Для того чтобы обеспечить сходимость последовательности точек x[t] к точке x*, в качестве последовательности Rt, t = 0,1,2,..., для методов внутренней точки следует выбирать монотонно убывающую сходящуюся к нулю последовательность положительных чисел, т.е. Rt +0 при t. При этом P(x, Rt ) f (x), x[t ] x*. Для вычисления Rt используется рекуррентное соотношение Rt = Rt -1 c, t = 1,2,…, где R0 > 0 (часто R0 = 1), c > 1 (часто c = 10 ).

Для методов внешней точки штрафные функции должны обладать следующими свойствами:

- во всех точках допустимого множества X внешние штрафные функции равны нулю;

- при выходе за пределы допустимого множества X внешние штрафные функции становятся положительными и достаточно быстро возрастают.

В качестве внешней штрафной функции часто используется штрафная функция типа квадрата "срезки" r (Rt, g(x)) = Rt gi (x).

i =Здесь gi (x) - "срезка" функции gi(x), определяемая следующим образом:

gi (x), если gi (x) 0, gi (x) = 0, если g (x) < 0.

i Для того чтобы обеспечить сходимость последовательности точек x[t] к точке x*, в качестве последовательности Rt, t = 0,1,2,..., для методов внешней точки следует выбирать монотонно возрастающую последовательность положительных чисел, т.е. Rt при t. Для вычисления Rt используется рекуррентное сооотношение Rt = Rt -1 c, t = 1,2,…, где R0 > 0 (часто R0 = 1), c > 1 (часто c = 10 ).

Метод штрафных функций позволяет в простых случаях явно (аналитически) решить задачу условной оптимизации. Рассмотрим в качестве иллюстрации аналитического решения следующий пример.

Пример. Дана задача условной минимизации f (x)= x min, x 2.

Легко видеть, что решением данной задачи является точка x*=2, при этом f=2.

Рассмотрим аналитическое решение задачи: а) методом внешней точки, б) методом внутренней точки (логарифмическая штрафная функция), в) методом внутренней точки (обратная штрафная функция).

Решение. Преобразуем ограничение исходной задачи к виду g(x)0:

g(x)= 2 - x 0.

а) Метод внешней точки.

Штрафная функция типа квадрата «срезки» имеет вид (R, g(x))= R 2 - x.

Pages:     | 1 |   ...   | 4 | 5 || 7 | 8 |



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.