Вариационное исчисление и методы оптимизации
Курс – 3, семестр — 5
- Часть 1. ВАРИАЦИОННОЕ ИСЧИСЛЕНИЕ
- Лекция № 3. Уравнение Эйлера для задачи Лагранжа
- 3.1. Задача Лагранжа
- 3.2. Уравнение Эйлера
- 3.3. Примеры
- 3.4. Частные случаи уравнения Эйлера
- 3.5. Задача падения тела
- 3.6. Принцип Ферма в геометрической оптике
- Выводы
- Задания на самостоятельную работу
- эйлера -лагранжа уравнение
- 📽️ Видео
Видео:Основы вариационного исчисления | уравнение Эйлера Лагранжа | 1Скачать
Часть 1. ВАРИАЦИОННОЕ ИСЧИСЛЕНИЕ
Видео:9. Метод вариации произвольной постоянной ( метод Лагранжа ). Линейные дифференциальные уравнения.Скачать
Лекция № 3. Уравнение Эйлера для задачи Лагранжа
Распространим полученные результаты на класс экстремальных задач, в которых требуется отыскать минимум не на числовом множестве, а на существенно более сложном объекте. Речь идет о минимизации не обычной функции, а функционала – отображения, определенного на множестве произвольной природы и принимающего значения из множества действительных чисел. Первым рассматриваемым примером является простейший вариант классической задачи Лагранжа, относящейся к вариационному исчислению. Конкретные реализации задачи Лагранжа рассматриваются в качестве семинарских занятий, а также заданий на самостоятельную работу. Имеются два частных случая, в которых имеется возможность достаточно легко понизить порядок уравнения Эйлера. В качестве практических приложений исследуются задача о падении тела под действием собственного веса и задача о распространении света в неоднородной среде.
3.1. Задача Лагранжа
где F – известная функция своих аргументов, а – неизвестная функция, удовлетворяющая граничным условиям
(3.1)
а и – заданные числа. Ставится следующая задача.
Задача 3.1. Найти функцию v, минимизирующую функционал I при выполнении
граничных условий (3.1).
Определение 3.1. Задача 3.1 называется задачей Лагранжа.
3.2. Уравнение Эйлера
Попытаемся свести рассматриваемую задачу к виду, пригодному для применения теоремы 2.1, т. е. к задаче 2.1. Предположим, что некоторая функция u является решением задачи Лагранжа. Определим функцию
где s – число, а h – достаточно гладкая функция, определенная на отрезке и удовлетворяющая однородным граничным условиям
(3.2)
Тогда величина наверняка удовлетворяет граничным условиям (3.2) (см. рис. 3.1).
Определение 3.1. Величина называется вариацией функции u.
Вопрос: Что дает введение вариации функции?
Очевидно, функционал I достигает своего минимума в точке u тогда и только тогда, когда число 0 является точкой минимума функции f.
Вывод: Задача минимизации функционала сводится
к задаче минимизации функции одной переменной.
В соответствии с теоремой 2.1 необходимым условием минимума (более точно, локального экстремума) функции в данной точке является равенство нулю его производной в этой точке, если таковая, конечно, существует. Будем полагать, что функция F дифференцируема по совокупности аргументов. Обозначим через и частные производные от функции F по второму и третьему аргументу. Тогда справедливо следующее утверждение.
Рис. 3.1. Вариация функции.
Лемма 3.1. Производная функции f в нуле равна
(3.3)
Доказательство. Определим величину
Пользуясь разложением в ряд Тейлора, находим значение
где при Отсюда следует равенство
После деления на s и перехода к пределу при получаем
Применяя интегрирование по частям с учетом граничных условий (3.2), будем иметь
В результате предшествующее равенство принимает следующий вид (3.3).
Определение 3.2. Производную функции f в нуле называют вариацией функционала I в точке u.
Вариацию функционала обозначают через Учитывая, что вариация функционала зависит также и от функции h, используют и более полное обозначение
Замечание 3.1. Здесь используется также термин производная по направлению.
Итак, из леммы 2.1 следует, что в том случае, когда u является решением задачи Лагранжа и существует вариация функционала I по любому направлению h, справедливо соотношение
(3.4)
для любых функций h, удовлетворяющих граничным условиям (1.6). Таким образом, справедливо следующее утверждение:
Теорема 3.1. Вариация функционала в точке его минимума обращается в нуль.
Замечание 3.2. Это утверждение служит обобщением теоремы Ферма, согласно которой производная функции в точке экстремума равна нулю.
Замечание 3.3. Практически все получаемые в дальнейшем утверждения вариационного
исчисления (см. лекции №№ 000) будут реализацией условия равенства нулю вариации функционала.
Для дальнейшего преобразования равенства (3.4) с целью практического использования этого результата используется следующее утверждение, называемое леммой Лагранжа – Эйлера или основной леммой вариационного исчисления.
Лемма 3.2. Если для некоторой функции g, непрерывной на отрезке справедливо равенство
(3.5)
для любой непрерывной функции h, то g тождественно равна нулю.
Доказательство. Рассмотрим произвольную точку и положительное число e, столь малое, чтобы выполнялось вложение
Подберем функцию h в равенстве (3.5) таким образом, чтобы она обращалась в нуль вне отрезка и удовлетворяла равенству (см. рис. 3.2). Тогда соотношение (3.5) принимает вид
Пользуясь теоремой о среднем, приходим к соотношению
где После деления на 2e и перехода к пределу при с учетом определения функции h установим равенство Отсюда в силу произвольности точки следуют утверждения леммы.
Рис. 3.2. Функция h в лемме 3.2.
Замечание 3.4. Мы выбрали такой (не самый простой) способ доказательства леммы, поскольку использованная здесь техника будет широко использоваться в дальнейшем для перехода от интегральных условий экстремума к поточечным. Классическая теорема о среднем является частным случаем теоремы Лебега.
Применяя лемму 3.2 к вытекающему из равенств (3.3), (3.4) соотношению
(3.6)
установим справедливость следующего утверждения.
Теорема 3.2. Если достаточно гладкая функция u является решением задачи Лагранжа, то она удовлетворяет уравнению Эйлера
(3.7)
Замечание 3.5. Мы не будем здесь и далее уточнять, что подразумевается под достаточной гладкостью рассматриваемых функций. На данном этапе исследования нас интересует исключительно принципиальная возможность сведения задачи Лагранжа к рассмотренной ранее задачи 3.1 нахождения экстремума функции.
Определение 3.3. Гладкое решение уравнения Эйлера называют экстремалью.
Замечание 3.7. Уравнение Эйлера реализует слабый минимум. Это означает, что функционал на его решении и не превосходит значения функционала на любой другой функции, достаточно близкой к и в норме пространства непрерывно дифференцируемых функций, т. е. гарантируется близость не только функций, но и их производных. Если же близость обеспечивается в смысле более широкого класса непрерывных функций, то получается понятие сильного минимума. Естественно, любой сильный минимум является слабым, но, вообще говоря, не наоборот, поскольку гарантировать близость в более узком классе функции легче, чем в более широком. В частности, изображенные на рис. 3.3 функции u, v и w достаточно близки как непрерывные функции. Однако они уже не будут таковыми, если оценивать не только близость самих функций, но и их производных. Мы не будем приводить условия сильного экстремума для задачи Лагранжа, поскольку при переходе к общим экстремальным задачам с произвольными функциональными пространствами понятия слабого и сильного экстремума не столь содержательны.
Рис. 3.3. Разная степень близость непрерывных функций.
Замечание 3.8. Как видно из формулы (3.6), вариация функционала оказывается линейной относительно функции h. Величина в формуле вариации функционала, умножаемая под интегралом на h (т. е. левая часть уравнения Эйлера), называется производной Гато функционала I в точке u.
Замечание 3.9. Вариационное исчисление предоставляет значительное количество условий экстремума в дополнение к уравнению Эйлера. Однако их изложение не входит в наши планы, если не считать рассматриваемое на заключительной стадии данной части курса (вариационного исчисления) условия Лежандра.
Вопрос: К какому типу уравнений относится уравнение Эйлера?
Очевидно, уравнение Эйлера является обыкновенным дифференциальным уравнением второго порядка, поскольку подынтегральная функция в задаче Лагранжа зависит от первой производной искомой функции, а второе слагаемое в левой части равенства (3.7) включает в себя дифференцирование относительно аргумента х.
Вопрос: Решение дифференциального уравнения второго порядка определяется с точностью до двух произвольных постоянных.
Как узнать их конкретные значения?
В дополнение к уравнению Эйлера имеются также краевые условия (3.1), которым
непременно должна удовлетворять искомая функция.
Вывод: Решение задачи Лагранжа удовлетворяет краевой задаче (3.1), (3.7).
Процесс практического решения конкретной задачи Лагранжа представлен в Табл. 3.1.
Табл. 3.1. Процесс практического решения конкретной задачи Лагранжа.
Задание конкретных значений подынтегральной
функции F, границ и граничных значений .
Приведение конкретной задачи
к стандартному виду.
Запись уравнения Эйлера.
Определение производных
подынтегральной функции и их подстановка в уравнение (3.7).
Решение общего решения уравнения Эйлера.
Определение решения
с точностью до двух констант.
Подстановка решения уравнения в краевые условия.
Определение неизвестных констант из условий (3.1).
Вычисление соответствующего значения функционала.
Определение значения интеграла для найденной экстремали.
Анализ полученных результатов.
Вообще говоря, не очевидно,
что получено решение задачи.
Воспользуемся описанной методикой для исследования некоторых частных случаев задачи Лагранжа.
3.3. Примеры
Для прояснения сути полученных результатов рассмотрим некоторые примеры.
Пример 3.1. Требуется минимизировать функционал
на множестве дважды непрерывно дифференцируемых функций, удовлетворяющих граничным условиям
В результате уравнение Эйлера принимает вид
Полученная краевая задача имеет решение которое и минимизирует рассматриваемый функционал. Действительно, данный функционал можно представить в следующем виде
Понятно, что его значения никак не могут быть меньше, чем -1, причем последнее может достигаться исключительно при
Пример 3.2. Рассматривается задача минимизации функционала
(3.8)
на множестве функций удовлетворяющих условиям
(3.9)
Мы имеем задачу Лагранжа с подынтегральной функцией
.
Соответствующее уравнение Эйлера имеет вид
Отсюда следует, что выражение в круглых скобках последнего равенства есть константа, т. е.
.
Таким образом, производная равна некоторой константе . Интегрируя равенство получаем
Для нахождения двух неизвестных констант имеем равенства (3.8). Получаем
В результате искомая функция определяется по формуле
Отметим, что функционал (3.8) выражает расстояние между точками, определяемыми равенствами (3.9), вдоль кривой Полученный результат оказывается вполне тривиальным: кратчайшим расстоянием между двумя точками оказывается прямая, соединяющая эти точки.
Замечание 3.10. Естественно, задачи максимизации рассмотренных функционалов приводят к тем же уравнениям Эйлера с теми же краевыми условиями. Однако решения соответствующих краевых задач (те же самые, что были получены выше), уже не будут решениями поставленных задач, доставляя минимума, а не максимум данным функционалам.
Вывод: Уравнение Эйлера является необходимым, но,
вообще говоря, не достаточным условием экстремума.
Полученный результат следует считать вполне естественным – ранее для существенно более простой задачи минимизации функции одной переменной общего вида мы также установили лишь необходимое условие экстремума.
Вопрос: Почему уравнение Эйлера оказалось лишь необходимым, но, вообще говоря, не достаточным условием экстремума?
Мы получили уравнение Эйлера как следствие равенства нулю производной соответствующей функции f в нуле. Однако этот результат сам является лишь необходимым условием экстремума. Таким образом, здесь уравнение Эйлера попросту наследует свойства условия Ферма.
3.4. Частные случаи уравнения Эйлера
Рассмотрим некоторые частные случаи уравнения Эйлера, когда последнее допускает более простой анализ. Пусть функция F не зависит от искомой функции и, т. е. Тогда уравнение Эйлера (3.7) принимает вид
Отсюда следует равенство
(3.10)
где – произвольная постоянная. Таким образом, если подынтегральное выражение в минимизируемом функционале не зависит от искомой функции и, то уравнение Эйлера легко сводится к дифференциальному уравнению (3.10) первого порядка. Именно с такой ситуацией мы столкнулись в последнем примере.
Второй специальный случай реализуется, когда функция F не зависит явным образом от х, т. е. Установим следующее соотношение
.
Выражение в правой части последнего равенства обращается в нуль в силу уравнения Эйлера. В результате приходим к соотношению
откуда следует равенство
(3.11)
где – произвольная постоянная. Таким образом, если подынтегральное выражение в минимизируемом функционале не зависит от независимой переменной х, то уравнение Эйлера также сводится к дифференциальному уравнению (3.11) первого порядка.
Вывод: Если подынтегральное выражение в задаче Лагранжа
не зависит явным образом от искомой функции или ее аргумента,
то уравнение Эйлера сводится к уравнению первого порядка.
Теперь рассмотрим некоторые физические примеры.
3.5. Задача падения тела
Рассмотрим процесс падения тела под действием собственного веса. В качестве функции состояния выбираем высоту тела у над землей. Определим его полную механическую энергию E(t) в произвольный момент времени t. Она складывается из кинетической энергии K(t) и потенциальной энергии U(t)
Потенциальная энергия представляет собой энергию сил тяготения и равна произведению веса тела Р на его высоту над землей
где m – масса тела, g – ускорение свободного падения, а знак «минус» обусловлен тем, что действие сил тяготения направлено в сторону, противоположную возрастанию координаты у. Кинетическая энергия пропорциональна квадрату скорости v движущегося тела.
В результате находим значение полной механической энергии падающего тела в произвольный момент времени
(3.12)
Рассмотрим некоторый интервал времени [t0 , t1] , на протяжении которого тело продолжает падать. При этом совсем не обязательно, чтобы момент времени t0 соответствовал началу, в t1 – концу падения тела. Предположим, что начальное и конечное положения тела известны и принимают некоторые значения х0 и х1 соответственно, т. е. справедливы равенства
Все траектории y = y(t) при tÎ[t0,t1] , удовлетворяющие условиям (3.13), будем называть допустимыми. Зададимся вопросом, какой из допустимых траекторий соответствует минимальные затраты энергии?
Если бы в процессе движения энергия тела не менялась и была бы равна некоторому значению Е*, то энергетические затраты за время от t0 до t1 были бы равны Е*(t1 – t0). Однако, как видно из формулы (3.12), по мере падения тела его энергия меняется, поскольку со временем происходит движение тела. Тогда для вычисления затрат энергии на данном интервале времени, соответствующей допустимой траектории х, следует проинтегрировать равенство (3.12) по времени. В результате находим величину
которая называется действием системы на интервале времени с допустимой траекторией y и выражает затраты энергии в процессе движения тела от точки х0 до точки х1 по траектории y = y(t).
Для решения полученной задачи Лагранжа воспользуемся уравнением Эйлера. В данном случае подынтегральная функция равна
Тогда уравнение Эйлера (3.7) имеет вид
Полученное соотношение соответствует второму закону Ньютона для прямолинейного движения для случая под действием силы тяготения. В частности, выражение есть вес тела. Итак, уравнением Эйлера для задачи минимизации действия системы оказывается классическое уравнение движения
Замечание 3.11. Используемый выше принцип наименьшего действия (эволюция системы осуществляется таким образом, чтобы затраты энергии в процессе движения были минимальны) является одним из глубочайших законов природы. Различные проявления этого закона будут рассматриваться и в последующих лекциях.
Замечание 3.12. То обстоятельство, что в задаче Лагранжа – основной задаче вариационного исчисления функционал зависит явным образом от искомой функции и ее первой производной в значительной степени связано с принципом наименьшего действия. Дело в том, что входящая в определение действия потенциальная энергия связана с координатой движущегося тела, а кинетическая энергия – с его скоростью, т. е. первой производной от координаты.
3.6. Принцип Ферма в геометрической оптике
В оптике известен принцип Ферма, согласно которому свет распространяется от одной точки к другой по такому пути, который соответствует минимальному времени на его преодоление. Отсюда, в частности, следует, что в однородной среде свет распространяется прямолинейно, поскольку скорость света в однородной среде постоянна. Попытаемся вывести из этого принципа законы преломления света в неоднородной среде. Пусть свет распространяется на плоскости по некоторой кривой , причем начальная и конечная точки известны:
(3.14)
Скорость движения света определяется по формуле
откуда следует равенство
Выше уже отмечалось, что пройденный путь по кривой характеризуется равенством
Тогда предшествующее равенство записывается в виде
Точка у нас соответствует начальному моменту времени а точка –конечному моменту времени Тогда, интегрируя предшествующее равенство от 0 до Т, получаем значение
(3.15)
Согласно принципу Ферма свет движется от одной заданной точке к другой таким образом, чтобы функционал (3.15), выражающий время движения, был наименьшим.
Предположим, что точка соответствует границе двух сред (см. рис. 3.4). В каждой из них скорость света постоянно. Таким образом, имеем следующее представление скорости света:
(3.16)
Рис. 3.4. Преломление света.
Отметим, что подынтегральное выражение для рассматриваемой вариационной задачи не зависит явным образом от независимой переменной х. Тогда уравнение Эйлера может быть сведено к уравнению первого порядка (3.11)
В рассматриваемом случае оно принимает вид
Отсюда следует равенство
(3.17)
.
Решая уравнения (3.17), можно найти искомый закон движения .
Отметим, однако, что из равенства (3.17) можно вывести законы преломления света. Действительно, как известно, производная есть тангенс угла наклона траектории, т. е. . Учитывая известную формулу
Тогда равенство (3.17) принимает вид
Отсюда следует соотношение
связывающее углы падения и преломления со скоростями света в рассматриваемых средах. Это и есть закон преломления света, называемый законом Снеллиуса.
Выводы
На основании полученных результатов можно сделать следующие выводы:
· Задача минимизации функционала может быть сведена к задаче минимизации функции одной переменной с помощью вариаций.
· Необходимым условием экстремума в задаче Лагранжа является уравнение
Эйлера.
· Уравнение Эйлера представляет собой обыкновенное дифференциальное уравнение второго порядка, которое решается совместно с заданными краевыми условиями.
· Решение уравнения Эйлера (т. е. экстремаль) может быть решением задачи
Лагранжа, но не обязано быть им.
· Задача Лагранжа возникает в механике и в оптике.
Задания на самостоятельную работу
Требуется найти функцию минимизирующую функционал
на множестве всех функций удовлетворяющих условиям
В Табл. 3.2 задаются значения параметров задачи для различных вариантов.
Табл. 3.2. Варианты параметров для самостоятельной работы.
Видео:Принцип наименьшего действия #2 - Уравнение Эйлера-ЛагранжаСкачать
эйлера -лагранжа уравнение
ЭЙЛЕРА -ЛАГРАНЖА УРАВНЕНИЕ -необходимое условие экстремума в задачах вариационного исчисления, полученное Л. Эйлером в 1744. Впоследствии, используя другой метод, это ур-ние вывел Ж. Лагранж (J. Lagrange) в 1759.
Пусть поставлена задача вариац. исчисления, состоящая в определении экстремума функционала
при известных условиях на концах
И пусть непрерывно дифференцируемая ф-ция x(t), есть решение задачи (1), (2). Тогда x(t)удовлетворяет Э.- Л. у.:
Ур-ние (3) можно записать в развёрнутом виде:
Гладкое решение ур-ния (3) [или (4)] наз. экстремалью. Если F xx =0 в точке (t, х), лежащей на экстремали, то в этой точке экстремаль имеет непрерывную 2-ю производную х. Экстремаль, во всех точках к-ройназ. неособенной. Для неособенной экстремали Э.- Л. у. можно записать в виде, разрешённом относительно 2-й производной х.
Решение вариац. задачи (1), (2) необязательно должно быть непрерывно дифференцируемым. В общем случае оптимальное решение x(t) может быть кусочно дифференцируемой ф-цией. Тогда в угл. точках х (t) должны выполняться необходимые условия Вейерштрасса — Эрдмана, обеспечивающие непрерывность при переходе через угл. точку выраженийа на отрезках между соседними угл. точками ф-ция должна удовлетворять
Э—Л. у. Кусочно гладкие линии, составленные из кусков экстремалей и удовлетворяющие в угл. точках условиям Вейерштрасса-Эрдмана, наз. ломаными экстремалями.
В общем случае дифференциальное Э.- Л. у. является ур-нием 2-го порядка и, следовательно, его общее решение зависит от двух произвольных постоянных
Эти произвольные постоянные можно определить из граничных условий (2):
Если рассматривается функционал, зависящий от неск. ф-ций,
то вместо одного Э.- Л. у. приходят к системе n Э—Л. у.:
Общее решение системы (7) зависит от 2n произвольных постоянных, к-рые определяются из заданных 2n граничных условий (для задачи с закреплёнными концами).
В случае вариац. задач с подвижными концами, в к-рых левый и правый концы экстремали могут смещаться по нек-рым заданным гиперповерхностям, недостающие граничные условия, позволяющие получить замкнутую систему соотношений типа (5), определяются с помощью необходимого условия трансверсальности. Для простейшей задачи типа (1), в к-рой точка
не фиксируется, а может принадлежать нек-рому множеству, условие трансверсальности записывается в виде
оно должно выполняться при любых значениях дифференциалов dt 1 , dx 1 , dt 2 , dx 2 , удовлетворяющих проварьиро-ванным граничным условиям. Если левый и правый концы экстремали могут смещаться вдоль заданных линий то в силу условий
и независимости вариаций dt 1 и dt 2 из (8) получают
Если уравнения линий, вдоль к-рых смещаются левый и правый концы экстремали, заданы в неявном виде то условие трансверсальности записывается так:
Если на один из концов экстремали не наложено никаких ограничений, то на этом конце в силу независимости соответствующих концевых вариаций dt и dx условие трансверсальности принимает вид
Для функционалов, содержащих производные высших порядков [а не только 1-го, как (1), (6)], необходимое условие, аналогичное Э—Л. у., записывается в виде диф-ференц. ур-ния Эйлера-Пуассона (см. [1 ]).
Для вариац. задач, в к-рых разыскивается экстремум функционалов, зависящих от ф-ций неск. переменных, аналогичное необходимое условие записывается в виде ур-ния Эйлера — Остроградского, представляющего собой дифференц. ур-ние с частными производными (см. [2]).
В случае вариац. задач на условный экстремум получение системы Э.- Л. у. связано с использованием множителей Лагранжа. Напр., для т. н. задачи Больца, в к-рой требуется найти экстремум функционала, зависящего от n ф-ций
при наличии дифференц. ограничений типа равенств и граничных условийс помощью множителей Лагранжаиз составляется ф-цияи Э—Л. у. записываются в виде
Т.о., оптимальное решение вариац. задачи (9) — (11) должно удовлетворять системе (12), причём первые т из этих ур-ний совпадают с заданными условиями связи (10). Используя дополнительно необходимое условие трансверсальности, получают замкнутую краевую задачу для определения решения вариац. задачи (9) — (II).
Помимо Э.- Л. у. и условий трансверсальности оптимальное решение вариац. задачи должно удовлетворять и др. необходимым условиям [условию Клебша (Лежанд-ра), условию Вейерштрасса и условию Якоби].
Лит.: 1) Ахиезер H. И., Лекции по вариационному исчислению, M., 1955; 2) Лаврентьев M. А., Люстерник Л. А., Курс вариационного исчисления, 2 изд., M.-Л., 1950.
📽️ Видео
Уравнения Лагранжа второго рода. Задача 1Скачать
Лагранжева механика I: уравнение Эйлера - ЛагранжаСкачать
№9. Элементы вариационного исчисления. Уравнения Лагранжа-Эйлера.Скачать
Решение системы дифференциальных уравнений методом ЭйлераСкачать
Основы вариационного исчисления | уравнение Эйлера Лагранжа | конкретные примеры | 1Скачать
#161. САМАЯ КРАСИВАЯ ФОРМУЛА В МАТЕМАТИКЕ — ФОРМУЛА ЭЙЛЕРА: e^(iπ)+1=0Скачать
Найти экстремаль функционалаСкачать
Метод множителей ЛагранжаСкачать
Условный экстремум и функция ЛагранжаСкачать
Основы вариационного исчисления | уравнение Эйлера Лагранжа | 5Скачать
Основы вариационного исчисления | многомерные вариационные задачи | уравнение Эйлера ЛагранжаСкачать
Уравнение Эйлера - bezbotvyСкачать
Уравнение ЛагранжаСкачать
Метод ЭйлераСкачать
Основы вариационного исчисления | уравнение Эйлера Лагранжа | конкретные примеры | 4Скачать
19. Метод вариации произвольных постоянных. Линейные неоднородные диф уравнения 2-го порядкаСкачать
Курс по ОДУ: Уравнения Клеро и Лагранжа | Занятие 8Скачать