umf
.pdfУравнения математической физики
А. И. Шафаревич
осень 2010 года
Глава 1
Уравнения в частных производных первого порядка
Общая постановка задачи решения уравнения в частных производных 1-го порядка такова: имеется пространство Rn и его элементы x Rn. Требуется найти скалярную функцию u(x), удовлетворяющую уравнению:
∂u |
= 0. |
L x, u(x), ∂x |
Вся развиваемая здесь теория будет локальной, то есть справедливой в окрестности (для простоты, шаровой окрестности) какой-либо точки x.
Несколько слов об обозначениях и терминологии. У пространств и многообразий при необходимости верхним индексом указывается размерность. Если F скалярная функция от x Rn (и, возможно, каких-то
ещё переменных), то |
∂F |
или Fx это вектор из Rn, составленный из производных |
∂F . Различий между |
|
∂x |
|
∂xj |
векторами и ковекторами не делается. Чтобы различить объединение двух n-мерных векторов ξ и η в один
Pn
2n-мерный (ξ, η) и их скалярное произведение j=1 ξj ηj , последнее обозначается как (ξ · η). Итак, рассмотрим простейший случай: линейное однородное уравнение.
1.1Линейные однородные уравнения
Определение 1.1. Линейное однородное уравнение в частных производных это уравнение вида
n |
|
|
X |
∂u = 0 |
(1.1) |
vj (x) |
∂xj
j=1
Здесь vj (x) заданные функции.
В данном разделе мы подходим к функциям ¾по-физически¿: забудем про особые точки и тому подобные неприятности. Считаем, что все функции бесконечно дифференцируемы: vj C∞.
Итак, vj образуют векторное поле v(x) = (v1(x), . . . , vn(x)) Rn.
Определение 1.2. Поле v называется характеристическим полем для (1.1), а его траектории решения системы x˙ = v(x) называются характеристиками.
Точкой здесь и везде далее обозначается дифференцирование по параметру t, который можно рассматривать как координату на одномерной траектории.
Утверждение 1.1. Функция u(x) является решением системы (1.1) тогда и только тогда, когда u постоянна вдоль траекторий поля v. Это условие можно записать как u˙ (x(t)) = 0. Такие функции u(x) называются первыми интегралами уравнений характеристики x˙ = v(x).
Доказательство. Рассмотрим траекторию (характеристику), задаваемую параметрически: x = x(t). Значение функции u на этой траектории будет выражаться как u(x(t)). Покажем, что u постоянна вдоль траектории:
|
n |
|
|
|
|
n |
|
|
d |
X |
∂u dxj |
|
X |
∂u |
vj = 0. |
||
|
u(x(t)) = |
∂xj dt |
= |
|
|
|||
dt |
|
j=1 |
∂xj |
|||||
|
j=1 |
|
|
|
|
|
|
2
1.1. Линейные однородные уравнения |
3 |
Замечание. Возьмём точку P Rn, такую, что v(P ) =6 0. В окрестности такой точки траектории поля v ведут себя как семейство параллельных прямых, и по теореме о выпрямлении векторного поля мы можем сделать такую замену координат, что они и превратятся в семейство параллельных прямых: в этих координатах будет v = (1, 0, . . . , 0). Если эти координаты z1, . . . , zn, то u = F (z2, . . . , zn) как раз даст произвольное решение уравнения ( ) (F произвольная функция). Правда, придётся выразить координаты z через координаты x.
А если всё-таки получать решение в исходных координатах x? Тогда нам по-прежнему требуются функции, постоянные вдоль траекторий. Любая их комбинация также будет постоянной вдоль траекторий, а значит, являться решением уравнения (по доказанному утверждению). А функция, постоянная вдоль траектории, это, собственно, первый интеграл. Так что если f1(x), . . . , fn−1(x) функционально независимые первые интегралы уравнения характеристики (x˙ = v), то есть f˙k(x(t)) = 0, а F произвольная функция, то u = F (f1(x), . . . , fn−1(x)) будет являться общим решением уравнения ( ).
Как видно, общее решение зависит от одной неизвестной функции. Здесь видна аналогия с обыкновенными дифференциальными уравнениями первого порядка, в общее решение которых всегда входила одна неизвестная константа. И бороться с этой неопределённостью можно так же, как и в случае ОДУ: путём задания начальных условий. К уравнению вида ( ) необходимо добавить условие, что на некотором многообразии функция имеет определённый вид.
Простейший пример.
(
∂u = 0 ∂x1
u|x1=0 = ϕ(x2, . . . , xn)
Здесь ϕ это какая-то заданная функция, x1 = 0 то самое многообразие, на котором задано начальное условие, а u, соответственно, искомая функция. Как видно из первого уравнения, функции-коэффициенты vj все равны 0, кроме v1 = 1. Запишем характеристическую систему:
x˙ 1 = 1
x˙ 2 = 0
.
..
x˙ n = 0
При таком раскладе у этой системы первыми интегралами (и притом независимыми) будут как раз x2, . . . , xn. (Ведь первый интеграл это как раз такая комбинация неизвестных, производная от которой равна нулю. А в системе ровно это и написано.) Тогда u = F (x2, . . . , xn). На многообразии x1 = 0 должно быть u = ϕ(x2 . . . , xn). Значит, попросту F = ϕ, и всюду u = ϕ(x2, . . . , xn).
Теперь определение.
Рассмотрим многообразие (поверхность) M Rn, dim M = n − 1.
Определение. Задачей Коши для ( ) называется задача отыскания функции u(x), такой, что u(x) удовлетворяет ( ), и u(x)|M = u0 заданная гладкая функция.
Всегда ли у задачи Коши существует решение и как его найти? Воспользуемся геометрическими соображениями. Выпустим из каждой точки многообразия M траекторию поля v (характеристику). Так как u вдоль траекторий постоянно, то на каждой траектории значение u полагается равным значению u0 в той точке на M, откуда эта траектория выпущена. Собственно, осталось параметризовать M как-нибудь, и добавить ещё одну координату вдоль траектории (всего получится n координат).
Помешать нам проделать все эти операции может только одно обстоятельство: если где-то траектория будет касаться M. Тогда решения не будет u должна быть постоянна вдоль траектории, а на M задана u0: кто знает, как она там себя ведёт, никаких гарантий постоянства уж точно нет.
Определение. Точка P M называется нехарактеристической, если v(P ) / TP M, то есть если проходящий через неё отрезок траектории не лежит в касательном пространстве (а, наоборот, торчит куда-то в сторону).
Замечание. Точки из некоторой окрестности любой нехарактеристической точки P тоже нехарактеристические. (Скажем, в P у v(P ) есть какая-то ненулевая нормальная компонента, и тогда в силу непрерывности она и в окрестности ненулевая, согласно теореме об устойчивости знака непрерывной функции, доказывавшейся в курсе матанализа.)
Теорема. Пусть P нехарактеристическая точка. Тогда в окрестности точки P задача Коши для ( ) имеет единственное решение.
Доказательство. Итак, параметризуем поверхность M, то есть введём на M координаты: уравнение M будет выглядеть как x = r(y). (y = (y1, . . . , yn−1) как раз и будут координаты на M.) И добавим ещё одну координату t вдоль траекторий, которые мы выпускаем из каждой точки M в окружающее пространство Rn. (Мы это можем сделать, так как все рассматриваемые нами точки из окрестности точки P замечательным
4 Глава 1. Уравнения в частных производных первого порядка
образом нехарактеристические.) Для траекторий имеем систему:
x˙ = v(x) x(0) = r(y)
Таким образом, здесь есть неизвестная n-мерная функция x(t) и заданные зависимости r(y), v(x), а также параметры y1, y2, . . . , yn−1. Это система обыкновенных дифференциальных уравнений (при этом все функции, как мы уговорились, гладкие). В окрестности точки P у неё должно быть бесконечно дифференцируемое решение x = x(t) = X(t, y). При t = 0 (в частности, в точке P ) оно существует гарантированно (и равно r(y)), а чтобы доказать, что оно и при каких-то ещё t существует (то есть в окрестности точки P ), снова воспользуемся теоремой о спрямлении векторного поля; из неё следует, что для этого должен быть не равен нулю
якобиан отображения X, то есть должны быть линейно независимы векторы |
∂X |
, |
∂X , . . . , |
∂X |
. Но позволь- |
|||
∂t |
|
|||||||
|
|
|
|
∂y1 |
|
∂yn−1 |
||
те, на этом самом начальном многообразии t = 0 это просто векторы v, |
∂r |
, . . . , |
∂r |
. Последние образуют |
||||
|
∂yn−1 |
|||||||
|
∂y1 |
|
|
|
|
|
в касательной (к этому многообразию) плоскости TP M базис, а v в TP M как раз не лежит (поскольку P нехарактеристическая точка). Значит, они действительно линейно независимы, и решение указанной системы действительно существует в окрестности точки P . Заодно по теореме об обратной функции мы можем выра-
зить y = y(x). Однако нам-то нужна функция u(x). Что ж, нам известно, что она постоянна на траекториях, |
|
поэтому: |
( u t=0 = u0(y) |
|
|
|
∂u = 0 |
|
∂t |
|
Решение сразу видно: раз при t = 0 u = u0(y), и u с изменением t не меняется, то и при любых t имеем u = u0(y) = u0(y(x)).
Откуда взялась эта u0(y)? Нам были известны значения функции u на многообразии M, а y это координаты на этом самом многообразии. Проблема состояла ровно в том, чтобы найти связь между x и y.
Итог.
1) Общее решение. Пишем характеристическую систему уравнений x˙ = v, ищем независимые первые интегралы f1(x), . . . , fn−1(x), и тогда u = F (f1(x), . . . , fn−1(x)) общее решение (F произвольная функция).
2) Задача Коши. Решаем систему x˙ = v, x(0) = r(y) с y как параметрами, находим x, выражаем y(x), и тогда решение задачи Коши: u(x) = u0(y(x)).
У нас теперь все теоремы будут примерно об этом же что в окрестности нехарактеристической точки задача Коши имеет единственное решение. Только каждый раз будет своё определение задачи Коши и нехарактеристической точки.
1.2Квазилинейные уравнения
Здесь терминология и обозначения очень похожи на терминологию и обозначения из предыдущего раздела, однако путать их не следует.
Определение. Квазилинейное уравнение в частных производных это уравнение вида
n |
|
|
|
X |
∂u |
|
|
aj (x, u) |
∂xj |
= b(x, u) |
( ) |
j=1 |
|
|
|
Здесь aj (x, u), b(x, u) заданные функции.
Чем здесь ситуация хуже, чем в случае линейных уравнений? Мы уже не можем сказать, что aj это векторное поле в Rn, поскольку к n координатам добавилась ещё одна, соответствующая последнему аргументу (назовём его w, чтобы показать, что в пространстве Rn(x,w+1 ) у нас n + 1 независимая координата, и
w вовсе не обязательно имеет какое-то отношение к неизвестной функции u(x)).
Зато если v(x, w) = (a1(x, w), . . . , an(x, w), b(x, w)), то v векторное поле в Rn(x,w+1 ), и каждой системе вида ( ) однозначное соответствует поле v.
Характеристиками, соответственно, называются траектории этого поля, то есть решения системы:
w˙ = b(x, w) |
, или dt (x, w) = v(x, w) |
|
x˙ = a(x, w) |
|
d |
Здесь уже нельзя сказать, что u постоянна вдоль траекторий. Поэтому делаем так. Пусть u(x) произвольная гладкая функция, а многообразие (поверхность) M её график в Rn(x,w+1 ) (то есть M задаётся уравнением w = u(x)).
Утверждение. u(x) удовлетворяет ( ) тогда и только тогда, когда P M v(P ) TP M.
1.2. Квазилинейные уравнения |
5 |
Доказательство. Итак, уравнение этой поверхности: Φ(x, w) = u(x)−w = 0. Утверждение, что некоторый вектор лежит в касательном пространстве к поверхности, равносильно тому, что форма dΦ зануляет этот вектор.
|
n |
|
|
|
n |
|
|
|
X |
∂u |
|
|
X |
∂u |
|
ξ T M dΦ(ξ) = 0 |
|
|
ξxj |
− ξw = 0 |
|
|
aj − b = 0 |
j=1 |
∂xj |
j=1 |
∂xj |
||||
|
|
|
|
|
|
Последнее равенство и есть уравнение ( ).
Отступление. Инвариантные поверхности векторных полей.
Пусть у нас есть какое-то пространство Rm и векторы z из него. Далее, пусть имеется векторное поле v(z) и k-мерная поверхность M.
Определение. Поверхность M инвариантна относительно векторного поля v, если P M v(P ) TP M. Как видно, предыдущее утверждение можно переформулировать так: u(x) удовлетворяет ( ) тогда и только тогда, когда поверхность, задаваемая уравнением w = u(x), инвариантна относительно поля v(x, w).
Так вот. Что можно сказать про инвариантные поверхности?
Утверждение. Пусть M инвариантная поверхность. Тогда вместе с любой своей точкой M содержит и отрезок траектории поля v, проходящий через эту точку.
Доказательство. M задаётся уравнениями:
Φ1(z) = 0
· · ·
Φm−k(z) = 0
Кроме того, так как M инвариантна относительно v, то все dΦj (v) = 0. Рассмотрим траектории z(t). Это, по определению, решения уравнения z˙ = v(z). Посмотрим, как на траектории изменяются функции Φj (z(t)).
dtd Φj (z(t)) = dΦj (dz/dt) = dΦj (v) = 0
Итак, вдоль траекторий значение Φj не изменяется. В точке P все Φj равны нулю (поскольку P M). Значит, и дальше на траектории значения Φj равны нулю, а стало быть, точки траектории лежат в M.
Теперь. Пусть имеется Rm, поверхность k−1 Rm и гладкое векторное поле v(z). Определение. Точка P нехарактеристическая, если v(P ) / TP .
Соответственно, точки из в некоторой окрестности точки P тоже нехарактеристические. Утверждение. Пусть P нехарактеристическая точка . Тогда в окрестности P существует единствен-
ная k-мерная поверхность M такая, что M и M инвариантна относительна v.
Доказательство. Введём на координаты y = (y1, . . . , yk−1). Рассмотрим z(t) траектории поля v. Выпустим эти траектории из каждой точки многообразия . При этом пересечению траектории с будет соответствовать координата t = 0, а y будут выступать в качестве параметров (то есть конкретный набор игреков задаёт конкретную траекторию z(t)). Тогда уравнения, которыми задаются траектории, будут такие:
z˙ = v(z)
z(0) = r(y)
Итак, у этой системы есть решения z = z(t) = Z(t, y). Из тех же соображений, что использовались при доказательстве аналогичного утверждения из предыдущего раздела, получаем, что якобиан отображения Z ненулевой, и оно является диффеоморфизмом (и оно само, и обратное к нему отображение являются непрерывно дифференцируемыми).
Собственно, вот оно, многообразие M: оно задаётся параметрами (t, y), и Z(t, y) даёт точки многообразия M. Первое из уравнений системы означает инвариантность, второе что лежит в M. Размерность M на единицу больше размерности , поскольку для задания M используется на 1 параметр больше.
Конец отступления. Вернёмся к квазилинейным уравнениям.
Итак, чтобы найти решение ( ), воспользуемся доказанным утверждением и будем искать поверхность Mn Rn(x,w+1 ), инвариантную относительно поля v. Тогда, если выразить последнюю ((n + 1)-ую) координату (w) через остальные (x), то полученная связь w = u(x) и будет искомой функцией u(x).
Определение. Задача Коши для ( ) это задача отыскания функции u(x) такой, что u(x) удовлетворяет ( ), и u|γ = u0, где γ заданная (n − 1)-мерная поверхность, u0 заданная функция.
Определение. Задача Коши называется нехарактеристической в точке P γ, если a(P, u0(P )) = (a1(P, u0(P )), . . . , an(
TP γ.
В этом определении нехарактеричности присутствует не только поверхность γ и функции aj , но и функция
u0.
Теорема. Пусть задача Коши нехарактеристична в точке P поверхности γ. Тогда в окрестности P у этой задачи Коши существует единственное решение.
6 |
Глава 1. Уравнения в частных производных первого порядка |
Доказательство. Введём на γ координаты y = (y1, . . . , yn−1). Соответственно, в пространстве Rn поверхность γ будет задаваться уравнением x = r(y), где r какая-то заданная функция. Теперь построим в Rn(x,w+1 ) график функции u0 размерности n − 1. Это многообразие задаётся уравнениями x = r(y), w = u0(y). А теперь из каждой точки выпускаем траектории векторного поля v. Сейчас мы это будем описывать в формулах.
Убедимся, что v не касается . Рассмотрим базис в касательной плоскости к γ: |
∂r |
и к : ( |
∂r |
, |
∂u0 ). |
∂yj |
|
||||
|
|
∂yj |
∂yj |
Мы знаем, что a не касается γ, так как наша задача Коши нехарактеристическая. Иными словами, a не раскладывается по базису в касательной плоскости к γ. Но v = (a, b), так что v тем более не раскладывается по базису в касательной плоскости к . Значит, v не касается и мы можем выпустить из каждой точки траектории поля v и получить многообразие M размерности на 1 больше (то есть n).
Осталось доказать, что уравнение M можно представить в виде w = u(x), то есть однозначно выразить последнюю координату через первые. Иными словами, предлагается использовать x как координаты на M.
M задаётся в координатах следующей системой уравнений:
x˙ = a(x, w)
w˙ = b(x, w)
x(0) = r(y)
w(0) = u0(y)
Решения этой системы: x = x(t) = X(t, y), w = w(t) = U(t, y). Сейчас в качестве координат на M используются (t, y), а мы хотим x; для этого докажем, что можно выразить y = y(x), t = t(x). Это так, если отобра-
жение X является диффеоморфизмом, а это верно, когда его якобиан ненулевой, то есть |
∂X , . . . , |
|
∂x |
, |
∂x |
||
|
|
∂t |
|||||
|
∂y1 |
∂yn−1 |
|||||
линейно независимы. Но в любой точке поверхности γ эту систему можно представить как |
∂r |
, . . . , |
∂r |
, a; |
|||
|
∂yn−1 |
||||||
|
|
∂y1 |
|
|
здесь первые векторы образуют базис в касательном пространстве, а последний в нём не лежит, а значит, все они линейно независимы (и в окрестности тоже). Получается, что мы можем составить функцию
w= u(x) = U(t(x), y(x)), что и требовалось.
Ещё раз.
1)Общее решение. Пишем характеристическую систему уравнений x˙ = a, w˙ = b, ищем независимые первые интегралы f1(x, w), . . . , fn(x, w), и тогда F (f1(x, u), . . . , fn(x, u)) = 0 общее решение (F произвольная функция). Да, неявное.
2) Задача Коши. Решаем систему x˙ = a, w˙ |
= b, x(0) = r(y), w(0) = u0(y) с y как параметрами, находим |
x и w, выражаем y(x) и t(x), и тогда решение |
задачи Коши: u(x) = w(t(x), y(x)). |
1.3 Уравнение Гамильтона–Якоби
Определение. Уравнением Гамильтона–Якоби называется уравнение вида
H(x, ∂x∂u ) = 0 |
( ) |
Здесь H заданная функция. |
|
Рассмотрим пространство R2n. Его элементы обозначим как (x, p), то есть первые n координат, объединённые в n-мерный вектор, будут соответствовать x из уравнения и составлять пространство решений, а вторые n координат объединяются в n-мерный вектор p.
Если теперь имеется функция u(x), то мы можем построить поверхность Λn R2n 1-график этой функции, заданный уравнением p = ∂u∂x (так как p n-мерный вектор, то это фактически n уравнений). Координаты точек этой поверхности, таким образом, будут (x, ∂u∂x ). В этом определении уравнение ( ) никак
не участвует, есть только пространство R2n |
|
и функция u(x). |
|
|
||||||
|
|
|
|
(x,p) |
|
|
|
|
|
|
Найдём стандартный базис в касательном пространстве TP Λ к точке P Λ. |
∂xi∂xn ) |
|||||||||
ei = ∂xi |
= ( ∂xi , |
∂xi |
|
∂x ) = (0, . . . , 1, . . . , 0, ∂xi∂x1 , . . . , |
||||||
|
∂(x, ∂u∂x ) |
|
∂x |
∂ |
|
∂u |
|
∂ |
2u |
∂2u |
Значения всех частных производных здесь берутся в точке P . |
|
|
||||||||
Соответственно, вектор ξ TP Λ можно разложить по этому базису: |
|
|
||||||||
|
|
|
|
|
|
n |
n |
|
|
|
|
|
|
|
|
|
X |
X |
|
|
|
|
|
|
|
ξ = ξj ej = ξxj ej |
|
|
||||
|
|
|
|
|
|
j=1 |
j=1 |
|
|
Последнее равенство верно из-за того, что ξx, равное (ξx1 , . . . , ξxn ), оказывается также равным (ξ1, . . . , ξn). В таком случае оказывается, что для любого ξ Λ
1.3. Уравнение Гамильтона–Якоби |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
7 |
|||
|
|
|
|
|
|
n |
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
|
X |
|
∂2u |
|
|
|
|
|
|
|
||
|
|
|
ξpi = ξxj ejpi |
= ξxj |
∂xi∂xj |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
j=1 |
|
|
j=1 |
|
|
|
|
|
|
|
|
|
|
|
Напомним, что ξi это i-я координата в разложении вектора ξ по стандартному базису в TP Λ, ξx и ξp |
|||||||||||||||||||
это соответственно первые n координат и вторые n координат 2n-мерного вектора ξ, ξxi |
его i-я координата |
||||||||||||||||||
из первых n. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Теперь на касательном пространстве T |
P |
R2n |
для каждой точки P |
|
R2n |
введём билинейную форму |
|||||||||||||
n |
|
|
|
|
|
|
|
|
|
|
|
|
(x,p) |
|
|
|
|||
ω = Pj=1 dpj dxj . Теперь пусть ξ, η TP Λ. Сосчитаем ω(ξ, η). |
|
∂xi∂xj |
|
|
|
||||||||||||||
ω(ξ, η) = |
(ξpi ηxi −ξxi ηpi ) = |
|
ξxj |
|
∂xi∂xj ηxi − ξxi |
|
ηxj |
= |
∂xi∂xj (ξxj ηxi −ξxi ηxj ) = 0 |
||||||||||
|
n |
n |
n |
|
|
|
|
|
|
|
n |
|
|
|
|
|
n |
|
|
|
X |
X X |
|
|
∂2u |
|
|
X |
|
|
∂2u |
|
X |
∂2u |
|||||
|
i=1 |
i=1 |
j=1 |
|
|
|
|
|
|
j=1 |
|
|
|
|
|
i,j=1 |
|
|
Итак, мы доказали, что ω(ξ, η) = 0 при любых ξ, η TP Λ. Иными словами, ω|Λ = 0. (Более точно было бы
написать ω|TP Λ = 0 в любой точке P Λ.) При этом, напомним, поверхность Λn задаётся уравнением p = |
∂x∂u , |
||
а форма ω = |
n |
задана в любом касательном пространстве TP R(2x,pn ) (и может быть ограничена |
|
j=1 dpj dxj |
P
на любое его подпространство).
Заметим ещё, что для того, чтобы Λ была элементарной поверхностью и все проделанные нами операции (начиная с нахождения стандартного базиса в Λ) имели смысл, необходимо, чтобы параметризация Λ переменными x была диффеоморфизмом, то есть между точками поверхности Λ и пространством векторов x существовало взаимно однозначное, в обе стороны гладкое соответствие. Иначе это можно сформулировать так: πx : Λn → Rnx проекция Λ на Rnx должна быть диффеоморфизмом. Это условие действительно выполняется благодаря тому, что Λ задана как график функции.
Вернёмся к уравнению ( ), функции H и пространству R2(x,pn ). Уравнение H(x, p) = 0 задаёт некоторую (2n − 1)-мерную поверхность . Это утверждение не касается функции u и потому не связано само по себе с уравнением ( ). Однако если предположить, что функция u удовлетворяет уравнению ( ) (то есть H(x, ∂u∂x ) =
0) и точка P = (x, p) Λ (то есть p = ∂u∂x ), то из этих предположений вытекает, что P = (x, p) , то есть H(x, p) = 0. Следовательно, если u решение ( ), то Λ .
Всё вышесказанное можно кратко сформулировать в виде утверждения.
Утверждение. Пусть функция u(x) гладкое решение уравнения ( ). Тогда поверхность Λn R2(x,pn ),
заданная уравнением p = ∂u , удовлетворяет свойствам:
Pn ∂x
1) ω|Λ = 0, где ω = j=1 dpj dxj .
2)Λn 2n−1, где поверхность задана уравнением H(x, p) = 0.
3)πx : Λn → Rnx диффеоморфизм.
Получается, что задача поиска функции u(x) сводится к задаче поиска соответствующей поверхности Λ, причём найти её можно именно из того соображения, что ω|Λ = 0.
Рассмотрим касательное пространство TP R2n к точке P из R2(x,pn ) (тоже 2n-мерное). Компоненты вектора из этого касательного пространства тоже будем разбивать на две группы по n компонент: ξ = (ξx, ξp) =
(ξ |
, . . . , ξ |
xn |
, ξ |
p1 |
, . . . , ξ |
pn |
). Тогда форма ω |
= |
n dp |
j |
|
dx будет применяться к векторам ξ, |
η из T |
P |
R2n |
x1 |
|
|
|
|
|
Pj=1 |
|
j |
|
|
|||||
следующим образом: |
|
n |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ω(ξ, η) = |
ξpj ηxj − ξxj ηpj = (ξp · ηx) − (ξx · ηp) |
|
|
|
j=1
Здесь ( · ) обычное скалярное произведение n-мерных векторов. Таким образом, в координатах внешняя билинейная форма ω записывается так: ω = 0 −E , где E единичная матрица.
E0
Отступление. Линейная алгебра в симплектических пространствах.
Определение. Линейное симплектическое пространство это линейное пространство L, в котором задана кососимметричная билинейная невырожденная форма ω.
Кососимметричность: для любых ξ, η L выполнено ω(ξ, η) = −ω(η, ξ)
Невырожденность: ω(ξ, η) = 0 для любого η из L тогда и только тогда, когда ξ = 0. Это утверждение равносильно тому, что матрица Aω формы ω в любом базисе невырождена (det Aω =6 0).
Определение. Форма ω называется симплектической формой. w(ξ, η) называется кососкалярным про-
изведением векторов ξ и η. |
|
|
|
|
|
Отметим, что в симплектическом пространстве привычное определение длины (как |
|
|
для про- |
||
|
(ξ · ξ) |
||||
странств с обычным скалярным произведением) теряет смысл, так как ω(ξ, ξ) = 0 для |
любого ξ |
|
L. |
||
|
p |
|
|
8 |
Глава 1. Уравнения в частных производных первого порядка |
|||||
|
Пример. В качестве пространства L возьмём R2n, а форму ω зададим так: ω(ξ, η) = |
n |
(ξi+nηi −ξiηi+n). |
|||
|
i=1 |
|||||
|
|
. , ξ |
, ξ |
, . . . , ξ |
|
), то |
Отметим, что если группировать координаты векторов по n штук: ξ = (ξx, ξp) = (ξx1 , . .P xn |
p1 |
|
pn |
|
ω(ξ, η) записывается как ω(ξ, η) = (ξp · ηx) − (ξx · ηp).
Вчастности, при n = 1 мы получим пространство R2, и тогда ω(ξ, η) = ξ1η2 − ξ2η1 ориентированная
площадь параллелограмма, натянутого на вектора ξ и η.
Утверждение. Размерность всякого симплектического линейного пространства чётная. Доказательство. Рассмотрим линейное пространство L с симплектической формой ω. Пусть A матрица
формы ω (в некотором базисе). Тогда AT = −A, и det A 6= 0. Пусть dim L = m. Заметим, что det A = det AT = det(−A) = (−1)m det A. Отсюда (−1)m = 1, и m чётно.
Определение. Базис e1, . . . , en, u1, . . . , un в L называется симплектическим, если для любых i и j вы-
полняется ω(ei, ej) = ω(ui, uj ) = 0, ω(ei, uj) = δij . Матрица кососкалярного произведения в таком базисе |
||
|
В |
|
выглядит так: A = 0 −E |
||
|
E |
0 . |
Утверждение. |
|
любом симплектическом пространстве существует симплектический базис. |
Доказательство. Возьмём e1 любой ненулевой вектор из симплектического пространства L, dim L = 2n. Тогда существует такой вектор v1, что ω(e1, v1) =6 0. В таком случае, если u1 = v1/ω(e1, v1), то ω(e1, u1) = 1. Рассмотрим теперь пространство L1 = {ξ L: ω(ξ, e1) = ω(ξ, u1) = 0} так называемое косоортогональное дополнение к подпространству he1, u1i, натянутому на векторы e1 и u1. (Косоортогональное дополнение будет определено ниже.) Размерность этого пространства на 2 меньше размерности пространства L: dim L1 = 2n−2, так как ω(ξ, e1) = 0 и ω(ξ, u1) = 0 это два независимых уравнения на 2n координат.
Покажем, что ω|L1 невырождена. Предположим, что это не так; тогда ξ L1 : ω(ξ, η) = 0 η L1. Но
тогда по определению L1 вектор ξ лежит в линейной оболочке he1, u1i векторов e1 и u1, то есть представляет собой их линейную комбинацию. Но раз существует такой ненулевой вектор ξ, лежащий одновременно в he1, u1i и его косоортогональном дополнении, то ω также вырождена и на he1, u1i. Однако в этом подпространстве матрица ω равна Aω|he1,u1i = . Эта матрица невырождена, значит, невырождена и сама форма ω на указанном подпространстве противоречие.
Итак, форма ω, ограниченная на подпространство L1, невырождена, и he1, u1i не лежит в L1. Тогда мы можем проделать в L1 все те же операции, что проделали в L. Так мы получим ещё одну пару векторов e2, u2 и подпространство L2 размерностью 2n − 4 (при этом, так как L2 подпространство L1, то he1, u1i не лежит и в L2); и так далее до тех пор, пока не дойдём до пространства Ln−1 размерностью 2 и последней пары векторов en, un. Они будут обладать заявленным свойствами, так как векторы из одной и той же пары будут давать в кососкалярном произведении единицу по построению этой пары векторов, а векторы из разных пар будут давать в кососкалярном произведении ноль по построению подпространств, в котором они находятся.
Определение. Косоортогональное дополнение MÜ к подпространству M L это множество всех векторов из L, которые дают ноль при кососкалярном умножении на любой вектор из M, то есть MÜ = {ξ L: ω(ξ, η) = 0 η M}.
Определение косоортогонального дополнения MÜ подпространства M симплектического пространства L, очевидно, аналогично определению ортогонального дополнения N подпространства N евклидова пространства H. Евклидово пространство разлагается в прямую сумму подпространства и его ортогонального дополнения: H = N N , откуда следует dim N = dim H − dim N и N ∩ N = {0}. Соблюдаются ли эти свойства для косоортогонального дополнения?
Утверждение. Пусть L 2n-мерное симплектическое пространство, M его k-мерное подпространство. Тогда dim MÜ = 2n − k.
Доказательство. Выберем набор векторов v1, . . . , v2n так, чтобы это был базис в L, и при этом набор v1, . . . , vk был бы также базисом в M. Тогда векторы ξ множества MÜ задаются системой уравнений ω(ξ, vj ) = 0, j = 1, . . . , k. Эта система линейно независима, так как строки её матрицы являются также строками матрицы формы ω, а она невырождена. Таким образом, мы имеем k независимых уравнений на 2n координат вектора ξ, следовательно, размерность решения будет 2n − k. (Линейность кососкалярного произведения гарантирует, что MÜ также будет подпространством.)
Таким образом, для размерности косоортогонального дополнения выполняется то же соотношение, что и для обычного ортогонального дополнения. Однако разлагается ли симплектическое пространство в прямую сумму M MÜ?
Это не так. Пример: пусть dim M = 1, то есть M = hξi все вектора подпространства M отличаются двух векторов из M:
лишь на множитель, η = aξ η M. Вычислим кососкалярное произведение любых |
Ü |
, то есть M |
|
M |
Ü |
. |
|||||||||
ω(η , η ) = a |
a |
ω(ξ, ξ) = 0. Таким образом, любой вектор из M также принадлежит и M |
|
|
|
||||||||||
1 2 |
1 |
2 |
|
L = M M |
Ü |
|
|
|
|
|
|
|
|
||
В этом случае разложение |
|
невозможно. |
Ü |
|
|
|
|
|
|
|
|||||
Задача. Пусть M L, и M0 = ker ω|M . Доказать, что M ∩ M |
|
= M0. |
|
|
|
|
|
|
|||||||
Определение. Пусть M подпространство L. Тогда оно: |
|
|
|
|
|
|
|
|
|||||||
изотропно, если M MÜ;Ü |
|
|
|
|
|
|
|
|
|
|
|
||||
коизотропно, если M MÜ |
. |
; |
|
|
|
|
|
|
|
|
|
|
|||
лагранжево, если M = M |
|
|
|
|
|
|
|
|
|
|
|
1.3. Уравнение Гамильтона–Якоби |
9 |
Приведём несколько примеров.
Пример изотропного пространства уже приведён: это любое такое M, что dim M = 1.
Если dim M = 2n − 1, то M коизотропно. Действительно, ведь если MÜ 6 M, то неизбежно M MÜ = L, а это не так.
Далее, пусть (e1, . . . , en, u1, . . . , un) симплектический базис в L. Тогда he1, . . . , eki изотропно для любого k < n (это следует из определения симплектического базиса). В то же время hek+1, . . . , en, u1, . . . , ukiÜ =
hu1, . . . , uki, так что hek+1, . . . , en, u1, . . . , uki коизотропно. Наконец, he1, . . . , eni лагранжево. Пусть dim L = 2n, dim M = k, dim MÜ = k′. k + k′ = 2n. Тогда:
если M изотропно, то k 6 n 6 k′; если M коизотропно, то k > n > k′; если M лагранжево, то k = n = k′.
Замечание. M изотропно тогда и только тогда, когда ω|M = 0 и dim M < n. M лагранжево тогда и только тогда, когда ω|M = 0 и dim M = n.
Конец отступления. Вернёмся к поиску поверхности Λ и решения уравнения Гамильтона–Якоби.
Касательное пространство TP R2n к любой точке P из R(2x,pn ) с формой ω = |
n |
как раз явля- |
|||||
j=1 dpj dxj |
|||||||
ется симплектическим. Далее, если P |
2n−1 |
, то благодаря (2n − 1)- |
мерности (и, соответственно, T |
P |
) |
||
|
P |
|
|
|
подпространство TP будет коизотропным. Кроме того, так как ω|Λ = 0, то TP Λ лагранжево. Определение. Поверхность, касательное пространство к которой в каждой точке лагранжево, сама на-
зывается лагранжевой.
Это новое определение, так как предыдущее определение лагранжевости относилось к симплектическим подпространствам и применялось в наших геометрических рассуждениях к касательным плоскостям, а не к поверхностям.
|
|
Замечание. Пусть N подпространство в 2n-мерномÜсимплектическом пространстве, dim N = 2n − 1, |
|||||||||
и |
M |
N, dim M = n, M лагранжево. Рассмотрим N . Так как dim N = 2n |
− |
1, то N коизотропно, и |
|||||||
Ü |
|
|
Ü |
Ü |
M |
Ü |
|
|
|||
N |
|
N; кроме того, dim N |
|
= 1. Так как M N, то N |
|
= M (M лагранжево). |
|||||
|
|
|
Ü |
TP Λ. |
|
|
|||||
|
|
В нашем случае N это TP , M это TP Λ, так что TP |
|
|
|
Сейчас мы привлечём к нашим рассуждениям симплектическую дифференциальную геометрию.
Итак, пространство TP 2n−1 коизотропно, то есть lP = (TP )Ü TP , при этом одномерная прямая lP задана в каждой точке поверхности ; все векторы из lP отличаются только на множитель, и любой из них задаёт одно и то же направление для каждой точки P . Поэтому, а также благодаря тому, что lP TP , lP можно рассматривать как поле на .
Определение. Поле lP характеристическое поле поверхности . Кривые в R2(x,pn ), касающиеся lP в точке P характеристики.
Характеристическое оно в том же смысле, в каком назывались характеристическими поля в предыдущих разделах: при сдвиге вдоль них сохраняются нужные свойства. Правда, докажем мы это, как обычно, потом. Но как получить это поле не геометрически, а в координатах?
Ищем такое поле, чтобы ω, будучи применённой к вектору этого поля и к любому вектору пространства
TP , давала ноль. |
|
|
, то в этом случае dH(ξ) = 0. Покажем, |
|||||
Итак, пусть ξ TP . Так как задана уравнением H(x, p) = 02n |
|
|
|
|
|
|||
что в этом случае существует такое векторное поле v, заданное в R(x,p), что ω(ξ, v) = 0. |
|
|
|
|
||||
Для этого просто попробуем устроить ω(ξ, v) = dH(ξ) = 0. |
|
|
|
|
|
|
||
Но при этом |
· ξx) + ( ∂p |
· ξp) = j=1 |
∂xj ξxj + |
∂pj ξpj |
|
|
|
|
dH(ξ) = ( ∂x |
|
|
|
|
||||
|
|
n |
|
|
|
|
|
|
∂H |
∂H |
X |
∂H |
∂H |
|
|
|
|
|
|
|
|
|
||||
и в то же время |
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
ω(ξ, v) = (ξp · vx) − (ξx · vp) = |
(ξpj vxj − ξxj vpj ) |
|
|
|
|
|||
|
|
j=1 |
|
|
|
|
|
|
|
|
|
∂H |
∂H |
|
∂H |
|
|
Как видно, условие ω(ξ, v) = dH(ξ) = 0 выполняется, если v = ∂p , − |
∂x , то есть если vxj |
= |
∂pj |
, vpj |
= |
−∂H . Такое векторное поле v называется гамильтоновым.
∂xj
Итак, мы построили характеристическое поле это гамильтоново поле v. Формально говоря, мы доказали, что в каждой точке P вектор v(P ) лежит в одномерном пространстве lP , и lP = hvi; соответственно, в дальнейшем v и lP можно использовать почти как синонимы.
Продолжаем поиск поверхности Λn и операции с поверхностью 2n−1, а также симплектической формой
ω.
Утверждение. Пусть Λ лагранжева поверхность, Λ . Тогда Λ инвариантна относительно характеристик, то есть P Λ lP TP Λ.
Доказательство. Λ TP Λ TP TP Ü TP ΛÜ, но при этом TP Ü и есть lP .
10 Глава 1. Уравнения в частных производных первого порядка
Инвариантность относительно некоторого известного векторного поля это уже весьма примечательное свойство. Если теперь задать поверхность размерности на 1 меньшей, чем Λ, такую, чтобы она обязательно содержалась в Λ, то Λ будет задана уже однозначно. (См. последнее утверждение о инвариантных поверхностях векторных полей в предыдущем разделе.) Назовём эту поверхность Λ0.
Геометрически это можно представить так, как будто мы выпускаем из каждой точки поверхности Λn0 −1 траектории характеристического поля и получаем поверхность Λn.
Для того, чтобы Λ могла содержать Λ0 и удовлетворять требуемым условиям, тем же условиям должна
удовлетворять и Λ0: Λn0 −1 2n−1, ω|Λ0 = 0, проекция πx : Λ0 → γ диффеоморфизм (γ кривая в пространстве Rnx ).
Определение. Лагранжева задача Коши это задача нахождения n-мерной поверхности Λ, такой, что Λ0 Λ, ω|Λ = 0. Здесь Λ0 заданная (n − 1)-мерная поверхность, Λn0 −1 2n−1, ω|Λ0 = 0.
Определение. Лагранжева задача Коши называется нехарактеристической в точке P Λ0, если v / TP Λ0.
Теорема. Пусть некоторая лагранжева задача Коши нехарактеристична в точке P0. Тогда в в окрестности точки P0 у этой лагранжевой задачи Коши существует единственное решение.
Доказательство. Мы уже знаем, что существует единственная поверхность Λ, которая инвариантна относительно характеристик и содержит Λ0, и, таким образом, может быть решением. Докажем, что она действительно является решением.
Для этого покажем, что ω|Λ0 со сдвигом = 0. В этом поможет следующее утверждение. Утверждение. Поток любого гамильтонова поля сохраняет форму ω.
Доказательство. Все 2n координат будем рассматривать вместе: (x, p) = z. Траектории поля v задаются уравнениями z˙ = v(z), z(0) = z0.
Собственно, поток векторного поля это гладкое отображение gt : z0 → z(t) сдвиг вдоль траектории. Его дифференциал действует в касательных пространствах: dgt : TzR2n → TgtzR2n. Итак, возьмём любые два вектора ξ и η из TzR2n и покажем, что отображение dgt не меняет результата формы ω, применённой к ним.
Пусть ω(dgt(ξ), dgt(η)) = f(t); покажем, что производная этой функции по t всюду равна нулю. Достаточно показать это для точки t = 0, так как начало отсчёта времени можно поместить в любую удобную нам точку, поскольку H не зависит от t. Впрочем, все последующие рассуждения можно проделать и без этого соображения, просто тогда вместо ξ и η нужно будет использовать dgt0 (ξ) и dgt0 (η) соответственно (и
показать, что утверждение верно при всех t0). |
t=0 dgt(ξ), dgt |
(η) t=0 |
|
dgt(ξ) |
t=0, |
dt t=0 dgt(η) |
|
|||||||||||||||||
|
dt t=0 = dt t=0 ω(dgt(ξ), dgt(η)) = ω dt |
+ ω |
|
|||||||||||||||||||||
|
df |
|
d |
|
|
|
|
|
|
|
|
d |
|
|
|
|
|
|
|
|
|
d |
|
|
|
|
|
|
|
t = 0 оператор g |
|
(и |
|
|
|
просто тождественным оператором, |
то есть |
||||||||||||
Заметим, что при |
|
|
|
|
t |
|
|
dgt) становится |
|
|
|
|
|
|
|
|
|
|
||||||
dgt(η) t=0 = η. |
|
|
|
df |
|
= ω |
|
d |
|
|
dgt(ξ), η + ω ξ, |
|
d |
|
dgt(η) |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
t=0 |
|
dt |
t=0 |
|
|
|
|
dt |
t=0 |
|
|
|
|
|
Действие преобразования dgt на векторы из касательного пространства TP R2n выглядит как домножение на матрицу Якоби этого преобразования:
∂z(t) dgt(ξ) = ∂z0 ξ
В то же время z˙ = v. При этом v явно не зависит от времени, а только от x и p (то есть от компонент z).
|
|
|
|
d |
|
|
∂v |
|
|
|
|
|
dt |
t=0 dgt(ξ) = |
∂z |
ξ |
|
|
dt |
t=0 = ω( ∂z ξ, η) + ω(ξ, ∂z η) |
||||||
|
df |
|
|
|
|
∂v |
|
∂v |
|
|
|
|
ξx |
|
|
Hp |
|
Наконец, вспомним, что z = (x, p), ξ = |
ξp , и поле v |
гамильтоново, поэтому v = −Hx (напомним, |
нижние индексы x и p при H обозначает производные по n соответствующим компонентам, а нижние индексы
xи p у 2n-мерных векторов первые или вторые n компонент). Тогда
∂v = |
Hpx |
Hpp |
|
∂z |
−Hxx |
−Hxp |
ξp |
|
|
||
∂z ξ = −Hxxξx − Hxp |
|||
∂v |
Hpxξx + Hppξp |
ω( ∂v∂z ξ, η) = ((−Hxxξx − Hxpξp) · ηx) − ((Hpxξx + Hppξp) · ηp)