УМК (1013374), страница 27

Файл №1013374 УМК (Учебно-методический комплекс) 27 страницаУМК (1013374) страница 272017-06-17СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 27)

СМЕШАННЫЕ ОГРАНИЧЕНИЯПостановка задачиДаны дважды непрерывно дифференцируемые целевая функция f (x) = f ( x1,…, xn )и функции ограничений типа равенств и неравенств: g j ( x) = 0 , j = 1, … , m ; g j ( x) ≤ 0 ,j = m + 1, … , p , определяющие множество допустимых решений X .Требуется исследовать функцию f ( x) на экстремум, т.е. определить точки x ∗ ∈ Xее локальных минимумов и максимумов на множестве X :f ( x ∗ ) = min f ( x ) ;x∈X⎧⎪ g j ( x) = 0, j = 1,… , m; m < nгде X = ⎨ xg j ( x) ≤ 0, j = m + 1,… , p⎪⎩f ( x ∗ ) = max f ( x) ,x∈X⎫⎪⎬.⎪⎭Алгоритм решения задачиШаг 1. Составить обобщенную функцию Лагранжа:pL ( x, λ 0 , λ ) = λ 0 f ( x ) + ∑ λ j g j ( x ) .j =1Шаг 2. Записать необходимые условия минимума (максимума) первого порядка:∂ L( x ∗ , λ ∗0 , λ ∗ )= 0,i = 1,… , n ;а)∂ xiб) g j ( x ∗ ) = 0 , j = 1, … , m ; g j ( x ∗ ) ≤ 0 , j = m + 1, … , p ;в) λ∗j ≥ 0 , j = m + 1, … , p (для минимума),λ∗j ≤ 0 , j = m + 1, … , p (для максимума);г) λ ∗j g j ( x ∗ ) = 0 , j = m + 1, … , p .Шаг 3.

Решить систему для двух случаев:Первый случай: λ∗0 = 0 .Второй случай: λ∗0 ≠ 0заменитьλ∗jλ∗0(приэтом поделить условия «а», «в», «г» наλ∗0 ина λ∗j ).В результате найти условно-стационарные точки x ∗ , выделив из них полученныепри λ∗0 ≠ 0 (они могут быть регулярными точками экстремума). В каждом из двух160случаев следует начинать с рассмотрения 2 p − m вариантов удовлетворения условия «г»дополняющей нежесткости.Шаг 4. Для выделенных на шаге 3 точек проверить выполнение достаточныхусловий экстремума первого или второго порядка.Для проверки достаточных условий первого порядка следует:а) определить число l ограничений-равенств и активных ограничений-неравенств;б) если l = n и λ∗j > 0 для всех j ∈ J a , т.е. для всех активных ограниченийнеравенств, то в точке x ∗ – локальный минимум.

Если l = n и λ∗j < 0 для всехj ∈ J a , то в точке x ∗ – локальный максимум. Если l < n или соответствующиемножители Лагранжа не удовлетворяют достаточным условиям первогопорядка, проверить достаточные условия второго порядка.Для проверки достаточных условий второго порядка следует:а) записать выражение для второго дифференциала классической функцииnnЛагранжа в точке ( x ∗ , λ ∗ ) : d 2 L( x ∗ , λ ∗ ) = ∑ ∑i =1 j =1∂ 2 L( x ∗ , λ ∗ )dx i dx j ;∂x i ∂x jб) записать условия, накладываемые на первые дифференциалы ограниченийравенств и активных в точке x ∗ ограничений-неравенств:n ∂ g (x∗ )j∗dg j ( x ) = ∑dx i = 0 , j = 1, … , m и j ∈ J a ; λ∗j > 0 ( λ∗j < 0 );∂ xii =1n∂ g j (x∗ )i =1∂ xidg j ( x ) = ∑∗dx i ≤ 0 , j ∈ J a , λ∗j = 0 ;в) исследовать знак второго дифференциала функции Лагранжа для ненулевых dx ,удовлетворяющих системе из п.б.

Если d 2 L( x ∗ , λ ∗ ) > 0 , то в точке x ∗ –условныйлокальный минимум. Если d 2 L( x ∗ , λ ∗ ) < 0 , то в точке x ∗ –условный локальный максимум. Если достаточные условия экстремума невыполняются, следует проверить выполнение необходимых условий второгопорядка (см. утверждение 3.10 – лекция 2), следуя аналогичной процедуре. Еслиони выполняются, то требуется дополнительное исследование, а если нет, то вточке x ∗ нет условного экстремума.Шаг 5. Вычислить значения функции в точках условного экстремума.Пример 3.

Найти условный экстремум в задачеf ( x) = x1 − x 22 → extr ,g1 ( x) = x1 − x 2 − 1 = 0 ,g 2 ( x) = x12 + x 22 − 5 ≤ 0 .† 1. Составим обобщенную функцию Лагранжа:L ( x, λ 0 , λ ) = λ 0 ( x1 − x 22 ) + λ 1 ( x1 − x 2 − 1) + λ 2 ( x12 + x 22 − 5) .2. Выпишем необходимые условия минимума и максимума первого порядка:∂ L (x , λ 0 , λ )∂ L (x , λ 0 , λ )а)= λ 0 + λ1 + 2λ 2 x1 = 0 ,= −2λ 0 x 2 − λ1 + 2λ 2 x 2 = 0 ;∂ x1∂ x2161б) x1 − x 2 − 1 = 0 , x12 + x 22 − 5 ≤ 0 ;в) λ 2 ≥ 0 (для минимума), λ 2 ≤ 0 (для максимума);г) λ 2 ( x12 + x 22 − 5) = 0 .3. Решим систему для двух случаев.Первый случай: λ 0 = 0 , тогда условие «а» имеет видλ1 + 2λ 2 x1 = 0,− λ1 + 2λ 2 x 2 = 0 .Рассмотрим два варианта удовлетворения условия «г»:1) λ 2 = 0 , тогда λ1 = 0 и не удовлетворяется условие утверждения 3.8;2) λ 2 ≠ 0 , тогда система уравненийx12 + x 22 − 5 = 0,x1 − x 2 − 1 = 0удовлетворяется в двух точках: x1 = 2, x 2 = 1; x1 = −1, x 2 = −2 .

Складывая двауравнения в условии «а», получаем 2λ 2 ( x1 + x 2 ) = 0 . Так как λ 2 ≠ 0 , то x1 = − x 2 , что неудовлетворяется в обеих найденных точках.Второй случай: λ 0 ≠ 0 . Поделив уравнения приведенной в п.2 системы на λ 0 иλλзаменив 1 на λ1 , 2 на λ 2 , запишем условие «a» в видеλ0λ0∂ L (x , λ )∂ L (x , λ )= 1 + λ1 + 2λ 2 x1 = 0 ,= −2 x 2 − λ1 + 2λ 2 x 2 = 0 .∂ x1∂ x2Остальные условия сохраняют вид.Рассмотрим два варианта удовлетворения условия «г»:1) λ 2 = 0 , тогда1 + λ1 = 0,−2 x 2 − λ1 = 0,x1 − x 2 − 1 = 0 .31Отсюда имеем λ1 = −1, x 2 = , x1 = .

В результате получили условно-стационарную2231точку А: x1∗ = , x 2∗ = , λ∗1 = −1, λ∗2 = 0 , в которой удовлетворяется необходимое22условие и минимума, и максимума;2) λ 2 ≠ 0 , тогдаx12 + x 22 − 5 = 0,x1 − x 2 − 1 = 0,1 + λ1 + 2λ 2 x1 = 0,−2 x 2 − λ 1 + 2 λ 2 x 2 = 0 .Получаем еще две условно-стационарные точки:1552В : x1∗ = 2, x 2∗ = 1, λ∗1 = − , λ∗2 = > 0 ; С : x1∗ = −1, x 2∗ = −2, λ∗1 = , λ∗2 = > 0 ,3636в которых удовлетворяются необходимые условия минимума.4. Проверим выполнение достаточных условий экстремума первого порядка.162В точке А ограничение-неравенство не является активным, поэтому l = 1 < n = 2 иусловия не выполняются. В точках В и С ограничение-неравенство активное, поэтомуl = n = 2 . В обеих точках λ∗2 > 0 , поэтому в них достигается условный локальныйминимум.Проверим достаточные условия экстремума второго порядка из методическихсоображений (в точке А это требуется обязательно).В точке А ограничение-неравенство не является активным:d 2 L ( A ) = 2λ∗2 dx12 + 2λ∗2 − 2 dx 22 = −2dx 22 ,()dg 1 ( A ) = dx1 − dx 2 = 0 .Отсюда имеем: dx1 = dx 2 и d 2 L ( A ) = −2dx12 < 0 при dx1 ≠ 0 .

Поэтому в точке А –локальный условный максимум.В точках B и C ограничение-неравенство активно.15d 2 L (B ) = dx12 − dx 22 ,33Исследуем точку В : dg1 (B ) = dx1 − dx 2 = 0,dg 2 (B ) = 2 x1∗ dx1 + 2 x 2∗ dx 2 = 4dx1 + 2dx 2 = 0 .Отсюда следует, что dx1 = dx 2 = 0 и d 2 L (B ) ≡ 0 , поэтому требуется дополнительноеисследование.51d 2 L ( C ) = dx12 − dx 22 ,33dg1 ( C ) = dx1 − dx 2 = 0,Исследуем точку С:dg 2 ( C ) = −2dx1 − 4dx 2 = 0.Отсюда имеем: dx1 = dx 2 = 0 и d 2 L (q ) ≡ 0 .

Требуется дополнительное исследование.Из рис. 4 следует, что в точке В – условный локальный минимум, а в точке С –условный глобальный минимум.55. Значения функции в точках экстремума: f ( A ) = , f (B ) = 1, f (C ) = −5 . „4x2g1 ( x) = 0f ( x) =1g 2 ( x) = 0B12−11−554A2x1− 1X−2CРис. 4163f ( x) = 1f ( x) = 0f ( x) = −5Занятие 3.

ЧИСЛЕННЫЕ МЕТОДЫ ПОИСКА БЕЗУСЛОВНОГОЭКСТРЕМУМА. МЕТОДЫ ПЕРВОГО ПОРЯДКАПостановка задачиПусть дана функция f ( x) , ограниченная снизу на множестве R n и имеющаянепрерывные частные производные во всех его точках.Требуется найти локальный минимум функции f ( x) на множестве допустимыхрешений X = R n , т.е. найти такую точку x ∗ ∈ R n , чтоf ( x ∗ ) = minn f ( x) .x∈ RА. МЕТОД ГРАДИЕНТНОГО СПУСКА С ПОСТОЯННЫМ ШАГОМАлгоритмШаг 1. Задать x 0 , 0 < ε < 1 , ε1 > 0 , ε 2 > 0 , M – предельное число итераций.

НайтиT⎛ ∂ f (x )∂ f (x ) ⎞⎟ .,...,градиент функции в произвольной точке ∇f ( x ) = ⎜⎜∂ x n ⎟⎠⎝ ∂ x1Шаг 2. Положить k = 0 .( )Шаг 3. Вычислить ∇f x k .( )Шаг 4. Проверить выполнение критерия окончания ∇f x k< ε1 :а) если критерий выполнен, расчет закончен, x ∗ = x k ;б) если критерий не выполнен, то перейти к шагу 5.Шаг 5. Проверить выполнение неравенства k ≥ M :а) если неравенство выполнено, то расчет окончен: x ∗ = x k ;б) если нет, то перейти к шагу 6.Шаг 6. Задать величину шага t k .( )Шаг 7. Вычислить x k +1 = x k − t k ∇f x k .Шаг 8.

Проверить выполнение условия() ( )f x k +1 − f x k < 0(или() ( )( )f x k +1 − f x k < − ε ∇f x k2):а) если условие выполнено, то перейти к шагу 9;tб) если условие не выполнено, положить t k = k и перейти к шагу 7.2Шаг 9. Проверить выполнение условийx k +1 − x k < ε 2 ,() ( )f x k +1 − f x k164< ε2 :а) если оба условия выполнены при текущем значении k и k = k − 1 , то расчетокончен, x ∗ = x k +1 ;б) если хотя бы одно из условий не выполнено, положить k = k + 1 и перейти кшагу 3.Геометрическая интерпретация метода для n = 2 приведена на рис.

1.C1 > C 2 > C 3x2f ( x ) = C1x0− ∇f ( x 0 )x1x∗x2f (x ) = C3f (x ) = C 20x1Рис. 1Пример 1. Найти локальный минимум функцииf ( x ) = 2 x12 + x1x 2 + x 2 2 .† I. Определим точку x k , в которой выполнен по крайней мере один из критериевокончания расчетов.1. Зададим x 0 , ε1 , ε 2 , M : x 0 = (0,5;1) , ε1 = 0,1 ; ε 2 = 0,15 ; M = 10 . Найдем граTдиент функции в произвольной точке ∇f ( x ) = (4 x1 + x 2 ; x1 + 2 x 2 )T .2. Положим k = 0 .( ) ( )∇f (x ) : ∇f (x )3 0 . Вычислим ∇f x 0 : ∇f x 0 = (3; 2,5) .4 0 . Вычислим00T= 3,9 > 0,1 .

Перейдем к шагу 5.5 0 . Проверим условие k ≥ M : k = 0 < 10 = M . Перейдем к шагу 6.6 0 . Зададим t 0 = 0,5 .7 0 . Вычислим x 1 : x 1 = (0,5;1) − 0,5(3; 2,5)T( )( )( )= ( −1;− 0,25) ; f x 1 = 2,31 .TT( ) ( )80 . Сравним f x 1 с f x 0 = 2 . Имеем f x 1 > f x 0 . Вывод: условие() ( )f x k +1 < f x k для k = 0 не выполняется. Зададим t 0 = 0,25 , перейдем к повторениюшагов 7, 8.7 01 . Вычислим x 1 : x 1 = (0,5;1) − 0,25(3; 2,5)TT165( )= ( − 0,25; 0,375) ; f x 1 = 0,171 .T( )( )( ) ( )− x и f (x ) − f (x ) := 0,976 > 0,15 ;f (x ) − f (x ) = 1,829 > 0,15 .801 .

Характеристики

Тип файла
PDF-файл
Размер
2,34 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Почему делать на заказ в разы дороже, чем купить готовую учебную работу на СтудИзбе? Наши учебные работы продаются каждый год, тогда как большинство заказов выполняются с нуля. Найдите подходящий учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6417
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее