Постановка задачи
Данная задача чаще всего формулируется так:
Найти экстремали функционала
J[y]=∫x0x1L(x,y,y′)dx
y(x0)=y0,y(x1)=y1
Шаблон решения
Вычисление производных. Уравнение Эйлера
Уравнение Эйлера для данного функционала:
ddx∂∂y′L(x,y,y′)=∂∂yL(x,y,y′)
Вычисляем производные:
∂∂yL(x,y,y′),
∂∂y′L(x,y,y′),
ddx∂∂y′L(x,y,y′)
Выписываем уравнение Эйлера и приводим его к максимально простому виду.
Решение уравнения
Решаем полученное уравнение и выписываем общее решение:
y=y(x,C1,C2).
Определение констант
Используя граничные условия, получаем систему для определения постоянных:
{y(x0,C1,C2)=y0y(x1,C1,C2)=y1
Решая данную систему, получаем значения констант
C1=C1(x0,x1,y0,y1),
C2=C2(x0,x1,y0,y1).
Уравнение экстремали
Подставляем найденные в предыдущем пункте значения констант в общее решение и получаем уравнение экстремали функционала
y∗(x)=y(x,C1(x0,x1,y0,y1),C2(x0,x1,y0,y1)).
Если подинтегральная функция L(x,y,y′) не зависит от первого аргумента, то уравнение Эйлера имеет первый интеграл
y′∂L(y,y′)∂y′−L(y,y′)=C1
Данное выражение определяет дифференциальное уравнение первого порядка, общее решение которого (зависящее от некоторой постоянной C2) совпадает с решением уравнения Эйлера (См. пример 2 ниже).
Примеры
Пример 1
Найти экстремали функционала
J[y]=∫0π/2(4ycosx+y′2−y2)dx
y(0)=0,y(π/2)=1.
Решение
- Уравнение Эйлера для данного функционала имеет вид:
ddx∂∂y′(4ycosx+y′2−y2=∂∂y(4ycosx+y′2−y2)
Вычислим производные:
∂∂y(4ycosx+y′2−y2)=4cosx−2y,∂∂y′(4ycosx+y′2−y2)=2y′,
ddx∂∂y′(4ycosx+y′2−y2)=ddx(2y′)=2y′′.
Таким образом, уравнение может быть записано как
2y′′=4cosx−2y⇔y′′+y=2cosx
- Уравнение Эйлера является неоднородным ЛОДУ второго порядка с постоянными коэффициентами.
Общее решение соответствующего однородного уравнения
y0=C1cosx+C2sinx.
Решение неоднородного уравнения ищем в виде
y1=Axcosx+Bxsinx.
Подстановка в уравнение Эйлера дает
y′′+y=−2Asinx+2Bcosx=2cosx.
Следовательно, A=0 и B=1.
Таким образом, общее решение уравнения Эйлера
y=C1cosx+C2sinx+xsinx.
- Используя граничные условия, получаем систему для определения постоянных
{C1cos0+C2sin0+0sin0=0C1cosπ/2+C2sinπ/2+(π/2)sinπ/2=1
Решая данную систему, получим
{C1=0C2+π/2=1⇔{C1=0C2=1−π/2
- Подставляем найденные в предыдущем пункте значения констант в общее решение и получаем уравнение экстремали функционала J
y∗(x)=(x+1−π2)sinx.
Пример 2
Найти экстремали функционала
J[y]=∫021+y′2ydx,
y(0)=0,y(2)=4.
Решение
1, 2. Подинтегральная функция L(x,y,y′)=1+y′2y не зависит от переменной x, поэтому уравнение Эйлера для данного функционала имеет первый интеграл
y′∂∂y′1+y′2y−1+y′2y=C,
где C – некоторая постоянная.
Вычислив частную производную и выполнив элементарные преобразования, получим уравнение
y′=C12y2−1,
где C1=1/C.
Отсюда следует, что
dx=ydyC12−y2.
Интегрируя, получим
x=∫ydyC12−y2=−12∫d(C12−y2)C12−y2=−C12−y2+C2.
Таким образом, общее решение уравнения Эйлера
y(x)=C12−(x−C2)2.
- Используя граничные условия, получаем систему для определения постоянных
{0=C12−(0−C2)24=C12−(2−C2)2⇔{C1=±C216=4C2−4⇔{C1=±5C2=5
- Подставляем найденные в предыдущем пункте значения констант в общее решение и получаем уравнение экстремали функционала J
y∗(x)=10x−x2.
Экстремаль является дугой окружности радиуса 5 с центром в точке (5,0).
Пример 3
Важно отметить, что существуют функционалы, у которых нет экстремалей, удовлетворяющих данным граничным условиям.
Найти экстремали функционала
J[y]=∫01(xy+y2−2y2y′)dx,
y(0)=y0,y(1)=y1.
Решение
Уравнение Эйлера для данного функционала имеет вид
ddx∂∂y′(xy+y2−2y2y′)=∂∂y(xy+y2−2y2y′)
Вычислим производные:
∂∂y(xy+y2−2y2y′)=x+2y−4yy′,∂∂y′(xy+y2−2y2y′)=−2y2,
ddx∂∂y′(xy+y2−2y2y′)=ddx(−2y2)=−4yy′.
Таким образом, уравнение Эйлера может быть записано как
x+2y−4yy′=−4yy′⇔y=−x2
и по сути не является дифференциальным.
Исходный функционал имеет только одну экстремаль y=−x2. Поставленная задача имеет решение только если y0=0 и y1=−12.
Вам нужно срочно заказать статью по математике для публикации? Обратитесь за помощью к нашим экспертам!
Комментарии