Речь пойдет о конкретной задаче из классического "Сборника задач по аналитической механике"
(Е. С. Пятницкий, Н. М. Трухан, Ю. И. Ханукаев, Г. Н. Яковенко; под ред. Е. С. Пятницкого, 4-е изд., Москва, МФТИ, 2018).
Формально она относится к вариационному исчислению, но по сути затрагивает одну из ключевых проблем аналитической механики - обратную задачу вариационного исчисления.
В задачнике формулируется вопрос, который можно переформулировать так:
Всякое ли скалярное (одномерное) дифференциальное уравнение второго порядка эквивалентно некоторому уравнению Лагранжа?
Эквивалентность понимается в строгом смысле: если функция (q(t)) является решением исходного уравнения второго порядка, то она же является решением и соответствующего уравнения Эйлера-Лагранжа, и наоборот; и это должно выполняться в некоторой окрестности рассматриваемой точки.
Официальный ответ, приведенный в задачнике, выглядит, мягко говоря, пугающе: он формален, практически не объясняет, откуда берется Лагранжиан, и почему он вообще должен существовать. Более того, при первом знакомстве не вполне ясно, существует ли в принципе функция (L), удовлетворяющая предъявленным требованиям, или это лишь красивая формулировка без реальной содержательной опоры.
Ниже изложим по сути авторское решение, то есть тот подход, который был предложен самим автором задачи (опирающимся, в частности, на идеи, изложенные в книге O. Bolza "Lectures on the Calculus of Variations", University of Chicago Press, 1904).
---
Формулировка проблемы в современном виде
Рассмотрим обыкновенное дифференциальное уравнение второго порядка для одной функции (q(t)):
[
ddot q = F(t, q, dot q),
]
где (F) - гладкая (достаточно дифференцируемая) функция своих аргументов в некоторой области пространства ((t, q, dot q)).
Вопрос:
можно ли для любого такого уравнения (локально, в малой окрестности некоторой точки ((t_0, q_0, dot q_0))) подобрать Лагранжиан (L(t, q, dot q)), такой что соответствующее уравнение Эйлера-Лагранжа
[
frac{d}{dt}frac{partial L}{partial dot q} - frac{partial L}{partial q} = 0
]
эквивалентно исходному уравнению (ddot q = F(t,q,dot q))?
Иначе говоря, спрашивается, всякое ли "разумное" уравнение второго порядка можно представить как уравнение Лагранжа с некоторой подходящей функцией (L). Подчеркнем: речь идет именно о локальном существовании (L) - в малой окрестности выбранной точки, а не на всей бесконечной области.
Ответ, который дается в задаче (и который, к сожалению, в задачнике оформлен крайне скупо и неинформативно), звучит так:
> Да. В малой окрестности любой точки фазового пространства ((t, q, dot q)) существует функция Лагранжа (L(t, q, dot q)), для которой уравнение Эйлера-Лагранжа эквивалентно заданному скалярному уравнению второго порядка.
Дальше - уже техника.
---
Вывод уравнения на Лагранжиан
Запишем уравнение Эйлера-Лагранжа в развернутом виде. Для функции (L(t,q,v)), где (v) мы временно пишем вместо (dot q), имеем:
[
frac{d}{dt}frac{partial L}{partial v} - frac{partial L}{partial q} = 0.
]
Подставим в него выражение для (ddot q) из исходного уравнения движения:
(ddot q = F(t, q, v)), где (v = dot q).
Тогда, выполняя дифференцирование по времени с учетом зависимости (q=q(t)), (v = dot q(t)), получаем (нижние индексы - частные производные):
[
frac{d}{dt} L_v = L_{vt} + L_{vq},dot q + L_{vv},ddot q
= L_{vt} + L_{vq},v + L_{vv} F(t,q,v).
]
Условие Эйлера-Лагранжа принимает вид:
[
L_{vt} + L_{vq},v + L_{vv} F(t,q,v) - L_q = 0.
tag{1}
]
Это уравнение в частных производных для неизвестной функции (L(t,q,v)). Цель - найти (L), удовлетворяющую (1), хотя бы локально.
Здесь впервые проявляется тонкость: (1) - нелинейное уравнение в частных производных, да еще и со смешанными производными. В общем виде выписать его решение в явном виде почти никогда нельзя, поэтому приходится действовать опосредованно.
---
Переход к вспомогательной функции
Чтобы упростить задачу, введем обозначение
[
P(t, q, v) = L_{vv}(t, q, v).
]
То есть (P) - вторая частная производная Лагранжиана по скорости. Именно эта величина часто интерпретируется в механике как "обобщенная масса" или коэффициент при квадратичном члене по скорости.
Теперь продифференцируем уравнение (1) по (v). Получим уравнение уже не для (L), а для (P). При аккуратном дифференцировании (1) по (v) возникает уравнение вида:
[
P_t + v,P_q + F(t,q,v),P_v + F_v(t,q,v),P = 0.
tag{2}
]
Это - линейное уравнение в частных производных первого порядка для функции (P). Его структура типична для транспортных и кинетических уравнений, и к нему вполне применим метод характеристик.
Чтобы уравнение (2) имело единственное решение, нужно задать подходящее начальное условие. Вводится условие вида
[
P(t_0, q, v) = P_0(q,v),
tag{3}
]
где (P_0(q,v)) - некоторая гладкая функция, определенная в окрестности точки ((q_0, v_0)); при этом ((t_0, q_0, v_0)) - та самая точка, в окрестности которой мы хотим построить Лагранжиан.
В задачнике начальное условие выбирается специальным образом, чтобы впоследствии обеспечить правильный вид решения; здесь важна лишь идея: начальное условие (3) задает данные на некоторой "поверхности" (например, на гиперплоскости (t = t_0)).
---
Что дает метод характеристик
Уравнение (2) - линейное уравнение первого порядка. Для него можно рассмотреть характеристическую систему:
[
begin{cases}
displaystyle
frac{dt}{ds} = 1,\[4pt]
displaystyle
frac{dq}{ds} = v,\[4pt]
displaystyle
frac{dv}{ds} = F(t,q,v),\[4pt]
displaystyle
frac{dP}{ds} = - F_v(t,q,v), P.
end{cases}
]
Первые три уравнения - это просто динамика исходной системы, записанная в расширенном пространстве ((t,q,v)). Четвертое описывает, как вдоль траектории в этом пространстве изменяется функция (P).
Теорема о существовании и единственности решений задачи Коши для системы обыкновенных дифференциальных уравнений гарантирует, что при гладкости (F) задача Коши для системы характеристик (с заданием (t,q,v) и (P) при (s=0)) имеет единственное решение, по крайней мере, в достаточно малой окрестности точки ((t_0, q_0, v_0, P_0)).
Отсюда следует, что и задача Коши для уравнения (2)-(3) имеет единственное локальное решение (P(t,q,v)). Явной формулы для (P) мы, как правило, не получим, однако существование и гладкость решения нам уже обеспечены.
---
Восстановление Лагранжиана по (P)
Зная (P = L_{vv}), можно восстановить сам Лагранжиан (L) двукратным интегрированием по переменной (v):
[
L(t,q,v) = int^v int^{tilde v} P(t,q,hat v), dhat v, dtilde v + a(t,q),v + b(t,q),
]
где (a(t,q)) и (b(t,q)) - произвольные гладкие функции, возникающие как "константы" интегрирования (поскольку мы интегрируем по (v)).
Однако тут скрыта тонкость: уравнение (1) мы перед этим дифференцировали по (v). Поэтому, построив (L) только из решения уравнения (2), мы обеспечиваем выполнение не исходного уравнения (1), а некоторого ослабленного варианта:
[
L_{vt} + v,L_{vq} + F L_{vv} - L_q = g(t,q),
tag{4}
]
где (g(t,q)) - некоторая функция, не зависящая от (v). Она появляется ровно потому, что при интегрировании по (v) "теряется" информация о слагаемых, зависящих только от ((t,q)).
То есть: построенный таким образом Лагранжиан не обязательно удовлетворяет исходному уравнению (1), а лишь уравнению (4) - с добавочным членом (g(t,q)).
---
Устранение лишней функции (g(t,q))
Ключевое наблюдение состоит в том, что мы обладаем свободой модифицировать Лагранжиан добавлением полной производной по времени от некоторой функции (phi(t,q)). Это классический факт аналитической механики:
Если заменить (L) на
[
tilde L(t,q,v) = L(t,q,v) + frac{d}{dt}phi(t,q) = L(t,q,v) + phi_t(t,q) + phi_q(t,q),v,
]
то уравнения Эйлера-Лагранжа не изменятся: обе функции (L) и (tilde L) порождают одни и те же уравнения движения.
Используя эту свободу, можно подобрать функцию (phi(t,q)) так, чтобы компенсировать лишний член (g(t,q)) в (4). При добавлении (phi_t + phi_q v) к Лагранжиану соответствующее уравнение для (tilde L) изменится: в левой части появятся дополнительные слагаемые, которые можно настроить так, чтобы "поглотить" (g(t,q)).
Иными словами, существует такая функция (phi(t,q)), что для (tilde L = L + frac{d}{dt}phi) уравнение (4) превратится точно в исходное (1) без добавочного члена. Тогда (tilde L) уже будет удовлетворять правильному уравнению Лагранжа, эквивалентному исходному дифференциальному уравнению второго порядка.
---
Локальность результата и роль начальных условий
Очень важно, что все рассуждения ведутся в рамках локального анализа.
Мы не утверждаем, что найденный Лагранжиан существует на всем множестве значений ((t,q,v)). Теорема о существовании решения задачи Коши для характеристической системы гарантирует решение лишь в достаточно малой окрестности исходной точки. Ровно это и требуется в задаче: показать, что в любой точке можно "развернуть" подходящий Лагранжиан.
Начальное условие для (P), то есть выбор функции (P_0(q,v)) в (3), также играет концептуальную роль. Фактически мы можем подобрать его так, чтобы обеспечить, например, ненулевую и положительную вторую производную по скорости (что соответствует "правильной" кинетической энергии с положительной "массой"). Это подчеркивает, что множество возможных Лагранжианов для данного уравнения движения вообще говоря не единственно, и дополнительными условиями (симметриями, знаками, физическим смыслом) можно выделять более "естественные" представители.
---
Что на самом деле говорит задача о механике
С математической точки зрения результат можно сформулировать так:
- Для любого достаточно гладкого скалярного ОДУ второго порядка (ddot q = F(t,q,dot q)) в окрестности произвольной точки существует Лагранжиан, порождающий это уравнение как уравнение Эйлера-Лагранжа.
С точки зрения физики и аналитической механики этот факт означает, что одномерное движение с заданным законом (ddot q = F(t,q,dot q)) всегда можно интерпретировать как движение по принципу наименьшего действия, по крайней мере локально.
Это довольно сильное утверждение: оно показывает, что вариационный подход в одномерном случае обладает максимальной общностью. Если вам дано любое разумное уравнение движения точки на прямой, его можно "восстановить" как результат некоторого вариационного принципа с соответствующим функционалом действия.
Разумеется, в многомерном случае (несколько степеней свободы) ситуация становится существенно сложнее. Там появляются глубокие интегрируемые условия (условия Гельмгольца и их обобщения), и далеко не всякая система уравнений второго порядка является лагранжевой. Но в одной степени свободы локально препятствий нет - в этом и заключается суть обсуждаемой задачи.
---
Почему ответ в задачнике кажется "ужасающим"
Причина, по которой официальный ответ в задачнике производит столь неприятное впечатление, в том, что он:
1. Практически не объясняет, зачем вводятся новые функции и какие за ними стоят идеи.
2. Сводится к быстрой формальной манипуляции: "продифференцируем, обозначим, применим метод характеристик, интегрируем дважды" - без комментария о смысле каждого шага.
3. Оставляет без пояснений ключевые моменты - например, появление функции (g(t,q)) и способ ее устранения добавлением полной производной.
В результате у читателя создается впечатление, что это не живое доказательство, а какой-то технический фокус. Между тем внутренняя логика решения довольно стройна:
- уравнение Эйлера-Лагранжа переписывается как уравнение на (L);
- после дифференцирования по скорости оно упрощается до линейного уравнения на (L_{vv});
- метод характеристик гарантирует локальное существование и единственность решения;
- двукратное интегрирование по скорости восстанавливает Лагранжиан;
- свобода добавления полной производной уничтожает лишний член, вызванный интегрированием.
---
Еще несколько замечаний по существу
1. Гладкость решения.
Предполагается, что функция (F(t,q,v)) достаточно гладкая (например, непрерывно дифференцируема нужное число раз). Тогда метод характеристик и теорема о существовании решений ОДУ обеспечивают гладкость найденной функции (P), а затем и (L).
2. Неявность формулы.
В практически всех нетривиальных примерах явное выражение для (L(t,q,v)) через элементарные функции получить невозможно. Но для теоретического результата это не важно: ключевым является факт существования и локальной единственности (с точностью до добавления полной производной и выбора начальных данных).
3. Свобода выбора Лагранжиана.
Помимо добавления полной производной, остается свобода выбора начальных условий для (P). Разные выборы могут привести к разным Лагранжианам, описывающим одно и то же уравнение движения. Это напоминает калибровочную свободу: динамика определяется не самим (L), а лишь его классом с точностью до "безопасных" преобразований.
4. Примерное построение.
В принципе, для конкретного уравнения (ddot q = F(t,q,dot q)) можно наметить алгоритм приближенного построения (L):
- численно решать характеристическую систему для (P),
- численно интегрировать (P) по (v),
- затем подстроить добавочную функцию (phi), добившись выполнения исходного уравнения Эйлера-Лагранжа.
На практике это громоздко, но концептуально осуществимо.
5. Связь с классическими учебниками.
Описанная схема восходит к работам начала XX века, в частности к лекциям Болцы по вариационному исчислению. Уже там обсуждалась обратная задача: по заданному дифференциальному уравнению восстановить вариационный функционал. Современный вид рассуждений мало отличается по сути, но опирается на стандартные теоремы о системах ОДУ и уравнениях в частных производных.
---
Итог
Задача из задачника по аналитической механике, на первый взгляд оформленная сухо и пугающе, на самом деле иллюстрирует важное и красивое утверждение:
- Любое достаточно гладкое одномерное уравнение второго порядка локально допускает лагранжево представление.
- Для его конструктивного получения нужно:
- переписать уравнение Эйлера-Лагранжа как уравнение (1) на (L);
- дифференцировать его по скорости и решить линейное уравнение (2) для (P = L_{vv}) методом характеристик;
- дважды проинтегрировать по скорости;
- скорректировать результат добавлением полной производной, устранив лишнюю функцию (g(t,q)).
Тем самым демонстрируется, что в одномерной механике любой закон движения можно рассматривать как следствие принципа наименьшего действия - возможно, с довольно экзотическим, но все-таки существующим Лагранжианом.
