Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
АДУС.pdf
Скачиваний:
119
Добавлен:
21.03.2016
Размер:
2.37 Mб
Скачать

Перейти к оглавлению на странице: 256

то есть f первый интеграл системы (2.18). Что и требовалось.

Построение интеграла в стационарном случае

Если H = H(q, p), то, как мы знаем, H(q, p) первый интеграл системы (2.18) (интеграл энергии, см. (2.2)). Предположим, что ϕ(q, p, t) еще один интеграл этих уравнений. Тогда по теореме Пуассона функция

 

 

f = {ϕ, H}

(2.26)

тоже будет интегралом этой системы.

 

Из равенства

 

 

 

 

 

∂ϕ

+ {ϕ, H} = 0,

(2.27)

 

 

 

∂ t

следует, что

 

∂ϕ

 

 

 

f = −

(2.28)

 

 

 

.

 

 

∂ t

Таким образом, мы доказали следующий результат:

Теорема 2.2. Если H = H(q, p), ϕ(q, p, t) первый интеграл

системы (2.18), а функции mϕ при m = 1, . . . , k непрерывно

∂ tm

дифференцируемы по всем своим аргументам, то они также будут интегралами этой системы (или постоянными).

Упражнение 2.2. Примените теорему Пуассона в задаче о движении точки в центральном поле.

§3. Метод Якоби решения канонических уравнений

Полный интеграл

Равенство

F(x, z, p) = 0,

(3.1)

где x = (x1, . . . , xn), p = (p1, . . . , pn), pj = ∂z/∂xj , j [1 : n], рассмотрим как уравнение относительно неизвестной вещественнозначной функции z от n независимых вещественных аргументов x1, . . . , xn . Его называют уравнением в частных производных первого порядка.

162

Перейти к оглавлению на странице: 256

Пусть D область в Rn . Любая функция z аргумента x, которая обращает равенство (3.1) в тождество при x D, называется

интегралом или частным интегралом уравнения (3.1) в области D. Полным интегралом уравнения (3.1) в области D называют такое его n-параметрическое семейство интегралов

 

z = z(x, a), x D,

a = (a1, . . . , an) A Rn,

(3.2)

что z, zx0

j C1(D × A), j [1 : n], а из равенств

 

 

 

z = z(x, a), p1 = zx0

1 (x, a), . . . , pn = zx0

n (x, a)

(3.3)

исключением параметров a1, . . . , an можно получить уравнение (3.1).

Отметим, что из этого определения не следует, что, во-первых, полный интеграл определяется уравнением (3.1) однозначно, и, во вторых, что из полного интеграла может быть получен любой частный интеграл в D. Если функция F не зависит явно от z , то полный интеграл может быть рассмотрен в виде z(x, a1, . . . , an−1)+an . Если функция F не зависит явно от неизвестной z и переменных x1, . . . , xk , k < n, то можно положить

 

 

k

 

 

z =

Xj

(3.4)

 

ajxj+ζ(xk+1, . . . , xn),

 

 

=1

 

тогда

 

 

 

pj = aj, j [1 : k], pk+j = ∂ζ/∂ xk+j, j [1 : n − k],

(3.5)

и уравнение (3.1) преобразуется к уравнению вида

 

Φ(xk+1, . . . , xn, a1, . . . , ak, ∂ζ/∂xk+1, . . . , ∂ζ/∂xn) = 0,

(3.6)

имеющему полный интеграл вида

 

 

ζ(xk+1, . . . , xn, a1, . . . , an−1) + an,

(3.7)

откуда следует, что

 

 

 

k

 

 

z =

Xj

 

(3.8)

ajxj+ζ(xk+1, . . . , xn, a1, . . . , an−1) + an

 

=1

 

 

163

Перейти к оглавлению на странице: 256

полный интеграл уравнения (3.1).

Уравнение Гамильтона-Якоби

Систему обыкновенных дифференциальных уравнений в симметричной форме

dx1/P1 = . . . = dxn/Pn = −dp1/X1 = . . . = −dpn/Xn,

(3.9)

где Pj = ∂F/∂pj, Xj = ∂F/∂xj, j [1 : n], называют уравнениями характеристик для уравнения в частных производных первого порядка вида:

 

 

 

F(x, p) = 0.

 

 

 

(3.10)

Записав канонические уравнения в симметричной форме

dq1

= . . . =

dqn

= −

dp1

= . . . = −

dpn

=

 

dt

 

 

 

 

 

 

,

∂H/∂ p1

∂H/∂ pn

∂H/∂ q1

∂H/∂ qn

1

 

 

 

 

 

 

 

 

(3.11)

можно видеть, что это система уравнений характеристик для уравнения в частных производных вида (3.10):

∂S

+ H(q1, . . . , qn, ∂S/∂ q1, . . . , ∂S/∂ qn, t) = 0.

(3.12)

∂ t

 

 

Равенство (3.12) называют уравнением Гамильтона-Якоби. Неизвестную функцию S аргументов q1, . . . , qn и t называют главной функцией Гамильтона. Так как левая часть уравнения (3.12) не зависит от S , то его полный интеграл можно записать в виде:

S = S(t, q1, . . . , qn, a1, . . . , an) + an+1.

(3.13)

Метод Якоби

Метод Якоби решения канонических уравнений (1.10) основан на предположении, что известен некоторый полный интеграл (3.13) уравнения (3.12).

Теорема 3.1. (Якоби) Если D R2n+1 область, а

S(t, q, a) C2(D) полный интеграл вида (3.13) уравнения (3.12), удовлетворяющий условию

 

2S

n

 

 

det

i,k=1

6= 0, (t, q, a) D,

(3.14)

∂ qi∂ ak

164

Перейти к оглавлению на странице: 256

то существует такая область D0 D × Rn , что равенства

p =

∂S

, b =

∂S

, (t, q, a, b) D0

(3.15)

 

 

 

∂ q

∂ a

представляют собой 2n независимых интегралов канонических уравнений (1.10), где b = (b1, . . . , bn), как и a = (a1, . . . , an), произвольные постоянные.

Доказательство. Если

(t0, q0, a0) D, p0

= ∂q

(t,q,a)=(t0,q0,a0) , b0 =

∂ a (t,q,a)=(t0,q0,a0) ,

 

 

∂S

∂S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то по теореме о неявных функциях, из условия (3.14) следует, что в некоторой окрестности точки (t0, q0, a0, b0) равенства (3.15) однозначно определяют дважды непрерывно дифференцируемые функции p, q аргументов t, a, b. Нам остается доказать, что эти p, q , рассматриваемые как функции времени, удовлетворяют каноническим уравнениям.

Подставим удовлетворяющую равенствам (3.15) функцию q во вторую группу этих равенств и полученные тождества bi = ∂S/∂ ai , i [1 : n] продифференцируем по t:

2S

 

n

2S

 

d qk

 

 

X

 

 

 

+

 

 

 

= 0, i [1 : n].

(3.16)

∂ t ∂ ai

k=1

∂ qk∂ ai

 

dt

 

 

 

 

 

 

 

Сдругой стороны, подставляя полный интеграл S (см. (3.13))

вуравнение Гамильтона-Якоби (3.12) и дифференцируя полученное тождество по ai , приходим к равенствам:

2S

 

n

∂H

 

2S

 

 

+

X

 

= 0, i [1 : n].

(3.17)

 

 

 

 

∂ t ∂ ai

k=1

∂(∂S/∂ qk)

 

∂ qk ∂ ai

 

 

 

 

 

 

 

Вычитая равенства (3.17) из соответствующих равенств (3.16) и учитывая равенства (3.15), получаем, что

n

2S

dq

 

∂H

= 0, i [1 : n],

 

 

(3.18)

k=1 ∂ qk ∂ ai

dtk

∂ pk

X

 

 

 

 

 

 

 

165

Перейти к оглавлению на странице: 256

а из этих равенств и условия (3.14) следует, что функции p, q , удовлетворяющие равенствам (3.15) удовлетворяют первой группе канонических уравнений. Подставим теперь удовлетворяющую равенствам (3.15) функцию q в первую группу этих равенств и полученные тождества pi = ∂S/∂ qi , i [1 : n] продифференцируем по t:

dpi

 

2S

n

2S dqk

 

 

X

 

 

=

 

+

 

 

 

, i [1 : n].

(3.19)

dt

∂ qi∂ t

k=1

∂ qi ∂ qk

 

dt

 

 

 

 

 

 

 

 

 

С другой стороны, подставляя полный интеграл S (см.(3.13)) в уравнение Гамильтона-Якоби (3.12) и дифференцируя полученное тождество по qi , приходим к равенствам:

2S ∂H

 

n

∂H

 

2S

 

 

X

 

 

 

+

 

+

 

 

 

= 0, i [1 : n]. (3.20)

∂ qi∂ t

∂ qi

k=1

∂(∂S/∂ qk)

 

∂ qk∂ qi

 

 

 

 

 

 

 

 

Вычитая равенства (3.20) из соответствующих равенств (3.19) и учитывая уже доказанные равенства первой группы канонических уравнений, получаем, что функции p, q (удовлетворяющие равенствам (3.15)) удовлетворяют второй группе канонических уравнений.

Что и требовалось.

166