Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Введение в теорию солитонов

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
874 Кб
Скачать

 

 

Уравнение ГЛМ

 

51

Функции Φ(λ, x) и a(λ) аналитичны по λ в верхней полуплоскости

Im λ > 0, причем имеет место оценка

 

 

 

Φ(λ, x)

eiλy e−iλ(x−y) = O λ1e−iλ(x−y)

,

λ → ∞.

 

 

 

a(λ)

Деформируем контур интегрирования в левой части (9.9) с вещественной оси на параллельную ей прямую −∞ < Re λ < +, Im λ > 0. Поскольку |e−iλ(x−y)| = eImλ(x−y) 0 при (x−y) < 0, Im λ → +, то в силу последней оценки экспонента пропадет и интеграл в левой части занулится.

Уравнение (9.9) принимает вид

 

+

 

 

 

0 =−∞ eiλy #

Θ(λ, x)e−iλx+r(λ)Θ(λ, x)$ dλ.

(9.10)

Будем искать его решение в виде

 

 

 

 

 

(9.11)

 

Θ(λ, x) = e−iλx+ K(x, y)e−iλy dy,

 

 

 

 

x

 

где K – новая неизвестная вещественная функция. Это представление справедливо в силу интегрального уравнения (7.12), выведенного выше при рассмотрении прямой задачи рассеяния:

 

Θ(λ, x) = eiλx+

sin λ(x−y)

u(y)Θ(λ, y)dy.

 

λ

 

x

Подставляя эту функцию в уравнение (9.10) и заменяя переменную интегрирования y на ξ, получим

0 = +eiλy

K(x, ξ)eiλξ + r(λ)

eiλx +

K(x, ξ)eiλξ dξ dλ.

−∞

x

 

x

 

(9.12) Первый двойной интеграл здесь – суть последовательно примененные прямое и обратное преобразования Фурье. Точнее, вводя функцию

Kˆ (x, y) =

0,

x > y,

 

K(x, y), x < y,

52

 

Лекция 9

 

имеем

+

 

 

 

+

 

 

 

eiλy

Kˆ (x, ξ)e−iλξ dξ dλ = 2πKˆ (x, y).

−∞

−∞

 

 

 

Уравнение (9.12) теперь перепишется при x < y так:

 

+

 

+

0 = 2πK(x, y) +

r(λ)e(x+y)+ r(λ)eiλy

K(x, ξ)eiλξ dξdλ.

−∞

 

 

−∞

x

Введем обозначение

 

 

 

 

 

def

1

+

 

 

 

 

 

F (z) =

r(λ)eiλz

 

 

2π

 

 

 

 

 

−∞

и в уравнении переставим порядок интегрирования в последнем двойном интеграле:

 

+

 

eiλy r(λ)

K(x, ξ)eiλξ dξdλ =

 

−∞

x

+

=

K(x, ξ) r(λ)e(ξ+y)dλdξ = K(x, ξ)2πF (ξ+y)dξ.

x

−∞

x

Окончательно получается уравнение Гельфанда–Левитана–Марченко

K(x, y)+F (x+y)+ K(x, ξ)F (ξ+y)= 0, x < y.

(9.13)

x

2. Формула обращения

Нам осталось сделать последний шаг в схеме метода обратной задачи рассеяния – перейти от решения K(x, y) уравнения ГЛМ к искомому решению уравнения КдФ u(x, t). Этот переход дается модификацией формулы обращения (9.6), выведенной в начале этой лекции.

Перепишем выражение (9.6) для потенциала u в терминах функции

Θ = χe−iλx:

 

 

 

 

 

 

 

 

u(x) =

lim 2(χ(λ, x)

1) =

lim 2(Θ(λ, x)e−iλx

1).

 

 

 

∂x λ

 

∂x λ

 

 

 

→∞

 

 

 

→∞

 

 

 

 

 

 

 

 

 

 

 

 

 

Формула обращения

 

 

 

 

 

 

 

 

 

53

Подставим сюда выражение (9.11) для Θ(λ, x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Θ(λ, x) = eiλx+

K(x, ξ)eiλξ dξ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что дает

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x) = ∂x λlim

2

K(x, ξ)e

(ξ

x)

=

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

→∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂x

λ→∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 x

K(x, ξ)d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

lim

 

 

 

 

 

e(ξ

 

 

x)

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2 ∂x λlim

 

 

K(x, x)

 

 

 

 

Kξ (x, ξ)

=

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

→∞

 

 

x

 

(ξ

x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

2

 

 

 

 

 

 

 

 

 

 

 

Kξ (x, x) +

 

 

 

 

 

 

∂xK(x, x) + 2 ∂x λlim

 

 

e

 

 

 

Kξξ (x, ξ)dξ .

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

(ξ

 

x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь последнее слагаемое обращается в нуль, поскольку интеграл в силу гладкости K(x, ξ) имеет порядок O(λ−N ) при λ → ∞, N – любое положительное число (N раз проинтегрировать по частям).

Таким образом, окончательно получим следующую формулу, которая называется формулой обращения:

u(x) = 2

 

∂x K(x, x).

(9.14)

Подводя итог всему сказанному, сформулируем стандартный алгоритм решения обратной задачи рассеяния для уравнения КдФ.

По данным рассеяния a(λ) = a(λ, 0) и b(λ, t) = b(λ, 0)e83t находим коэффициент отражения

r(λ, t) = b(λ, t) = b(λ, 0) e83t. a(λ, t) a(λ, 0)

54

Лекция 9

Вычисляем ядро уравнения ГЛМ

 

F (z, t) =

 

1

 

+

 

 

r(λ, t)eiλz dλ.

 

2π

 

 

 

 

 

 

−∞

Решаем уравнение ГЛМ (9.13)

K(x, y, t)+F (x+y, t)+ K(x, ξ, t)F (ξ+y, t)= 0, x < y.

x

Находим решение уравнения КдФ по формуле обращения (9.14)

u(x, t) = 2 ∂xK(x, x, t).

Лекция 10

1. Дискретный спектр уравнения Шредингера

При выводе уравнения Гельфанда–Левитана–Марченко (9.13)

K(x, y) + F (x + y) + K(x, ξ)F (ξ + y)= 0

x

было существенно использовано условие a(λ) = 0, Imλ > 0. Это условие необходимо, чтобы сдвинуть контур интегрирования по λ для интеграла от левой части равенства

 

Φ(x, λ)

 

 

 

 

 

 

 

 

= Ψ(x, λ) + r(λ)Ψ(x, λ).

(10.1)

 

 

 

a(λ)

 

 

В силу аналитичности функции

Φ(x, λ)

 

можно сдвинуть контур интегри-

a(λ)

 

 

 

 

 

 

 

рования на Im(λ) вверх, а затем Im(λ) устремить в бесконечность. Если a(λ) имеет конечное число нулей в точках n, λn > 0, n = 1, 2, . . . , N , то справедливо разложение Тейлора

a(λ) = (λ − iλn)a (n) + O(λ − iλn)2, λ → iλn.

Предполагая a (n) = 0, имеем по теореме о вычетах

 

 

 

a(λ)

n=1

λ=n

a(λ)

 

=2πi n=1

a (n)

Φ(x, λ)eiλy

N

 

Φ(x, λ)e−iλy

N

Φ(x, iλn)e−λny

−∞

 

 

 

 

 

 

 

 

 

 

 

= 2πi

 

Res

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

(10.2)

В дальнейшем мы воспользуемся этой формулой для обобщения уравнения Гельфанда–Левитана–Марченко.

Введем важное определение дискретного спектра уравнения Шредин-

гера

 

Φxx + (λ2 + u(x))Φ = 0.

(10.3)

56

Лекция 10

Оказывается, что в классе быстроубывающих потенциалов

(1 + |x|)|u(x)|dx < ∞

−∞

уравнение (10.3), кроме знакомых нам осциллирующих решений

 

 

 

eiλx,

x → −∞,

Φ

 

 

(λ)eiλx +

b

(λ)e−iλx,

x → +∞,

 

a

имеет экспоненциально убывающие по x решения Φ(x, λ) при некоторых (комплексных) значениях λ. Такие решения соответствуют разрешенным состояниям квантовомеханической системы, описываемой уравнением (10.3), которое, как уже отмечалось выше, является одномерным стационарным уравнением Шредингера.

Точное определение дискретного спектра состоит в следующем. ОПРЕДЕЛЕНИЕ. Значение λ C принадлежит дискретному спектру

уравнения (10.3) с потенциалом u(x), если Φ(x, λ) L2(R), то есть

|Φ(x, λ)|2 dx < const.

−∞

Эти значения λ называются собственными значениями, а соответствующие решения – собственными функциями дифференциального оператора Шредингера d2/dx2 + u(x).

1.1. Свойства дискретного спектра

1.Точки дискретного спектра совпадают с нулями коэффициента про-

хождения, то есть λn - собственное значение тогда и только тогда, когда a(λn) = 0.

Для доказательства воспользуемся выражением a(λ) через вронскиан (см. Свойство 3 данных рассеяния)

W , Ψ) a(λ) = ,

2

откуда при a(λ) = 0, Im λ > 0 получим W (Φ(x, λn ), Ψ(x, λn)) = 0. Это означает, что Ψ и Φ линейно зависимы, то есть существует константа bn такая, что

Φ(x, λn) = bnΨ(x, λn).

 

 

 

Дискретный спектр

 

 

 

57

Оценивая эти функции при больших значениях x, имеем

 

 

 

 

 

|Φ(x, λn)| e−iλnx

= eImλnx,Im

Imλn > 0, x → −∞,

+

 

,

|

Φ(x, λn) = bn

||

Ψ(x, λn)

bn e− λnx, Imλn > 0, x

| |

| | |

 

 

 

откуда следует, по определению, что Ψ и Φ есть собственные функции,

а

λn – собственное значение. Обратно, из условия убывания

Ψ(x, λn)

на

бесконечности по x и постоянства вронскиана, следует

равенство

W (Φ(x, λn ), Ψ(x, λn)) = 0, что эквивалентно a(λ) = 0.

2. Теорема Марченко. Если u(x) – гладкий потенциал и

+

(1+ | x |) | u(x) | dx < const,

−∞

тогда множество собственных значений n} не более чем конечно.

В частности, если интеграл меньше ln(2 + 3), то дискретного спектра нет, то есть a(λ) нигде не обращается в нуль.

Доказательство этой теоремы достаточно сложно и выходит за рамки данного курса.

3.При вещественном u(x) все собственные значения лежат на мнимой оси:

λn = iκ n, κ n > 0, n = 1, 2, . . . , N.

Для доказательства воспользуемся тем, что если Φ(x, λn) – собственная

функция, то Φ(x, −λn) – тоже собственная функция,

Φxx + ((−λn)2 + u(x))Φ = 0,

поскольку u(x) = u(x). Граничные значения Φ(x, λn ) и Φ(x, −λn) также

совпадают:

Φ(x, λn) enx, x → −∞,

 

 

 

 

n)

 

 

 

nx = enx, x

 

Φ(x,

λ

e−iλ

.

 

 

 

 

 

 

 

→ −∞

Таким образом, эти функции совпадают и, следовательно, совпадают их

собственные значения λn = −λn, а это может быть, если они лежат на мнимой оси.

58

Лекция 10

4.Для различных собственных значений соответствующие собственные функции ортогональны:

+

Φn, Φm = Φ(x, λn)Φ(x, λm)dx = 0, λn = λm.

−∞

Для доказательства определим Φn = Φ(x, λn), Φm = Φ(x, λm), λn = = iκ n, λm = iκ m. Функции Φn и Φm удовлетворяют уравнению (10.3)

n)xx + (κ n2

+ u(x))Φn = 0,

(10.4)

m)xx + (κ m2

+ u(x))Φm = 0.

 

Умножим (10.4) на Φm и

проинтегрируем

обе части на

интервале

(−∞, +).

 

 

 

 

Из свойства 3 следует, что собственные функции на бесконечности

экспоненциально убывают

 

 

 

 

Φn eκ n|x|,

Φm eκ m|x|,

x → ±∞,

 

поэтому все интегралы ниже корректно определены. Интегрируя по частям, имеем

+

0 = Φm(x) #n)xx + (κ 2n + u(x))Φn(x)$ dx =

−∞

+

=#Φm(x)d n)x + (κ 2n + u(x))ΦmΦndx$ =

−∞

 

 

 

+

 

 

+

 

2

 

= Φmn)x −∞ +

#

n)xm)x + (κ n + u(x))ΦmΦn

dx =

−∞

 

$

 

 

 

 

 

 

 

 

+

=#(Φm)xdn) + (κ 2n + u(x))ΦmΦndx$ =

−∞

 

 

+

+

$dx =

2

= (Φm)xΦn −∞ +

Φn #m)xx + (κ n + u(x))Φm

 

 

 

 

 

 

−∞

Дискретный спектр

59

+

=Φn(x)[(Φm)xx + (κ 2n + u(x))Φm + (κ 2n + κ 2mm]dx =

−∞

= (κ n2 + κ m2 )

+

ΦnΦmdx.

 

−∞

Поскольку по условию κ n = κ m, последний интеграл равен нулю, что и требовалось.

5. Скалярный квадрат собственной функции равен

Φn(x, iκ n) 2 = Φn, Φn = ibna (iκ n),

где

 

Φn eκ nx,

 

x → −∞,

 

 

 

 

Φn → bneκ nx, x → +∞.

Для доказательства продифференцируем уравнение (10.3) по λ:

 

2

 

0 =

 

n)xx + (λ

 

+ u(x))Φn =

∂λ

 

 

= (Φ)xx + (λ2 + u(x))Φ+ 2λΦn.

Умножим последнее уравнение на Φn и проинтегрируем по x на ин-

тервале (−∞, +). Интегрируя по частям, имеем

 

+

+

 

2λ−∞ Φn2 dx =−∞ Φn #(Φ)xx + (λ2 + u(x))ΦnΦ)$dx =

 

+

 

 

+

 

 

= Φn)x −∞ +

n)x)x + (λ2 + u(x))ΦnΦ) dx =

 

−∞

#

$

 

 

 

 

 

 

 

 

 

+

 

= bnκ na (iκ n) +

#n)xd) + (λ2 + u(x))ΦnΦ)$dx =

 

= bnκ na (iκ n) n)x) +

+

 

 

 

 

−∞

60

Лекция 10

+

+#n)xxΦ+ (λ2 + u(x))ΦnΦ)$dx =

−∞

+

=Φ#n)xx + (λ2 + u(x))Φn$dx + 2κ nbna (iκ n) =

−∞

= 2κ nbna (iκ n).

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]