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

УМФ ЛК

.pdf
Скачиваний:
17
Добавлен:
17.05.2015
Размер:
1.57 Mб
Скачать

61

 

1

 

 

 

 

 

4 (x) , что и требовалось.

 

r

 

 

Определение (фундаментальное решение оператора Лапласа).

Фундаментальным решением оператора Лапласа в Rx3 называют функцию

(x, ) D': (x, ) (x ), x, R3 (8)

Замечание 1.

Из примера вытекает, что фундаментальным решением оператора Лапласа является

(x, )

1

1

 

 

(9)

4

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это вытекает из того, что

1

x

4 (x ) в силу примера.

Замечание 2 (решение уравнения Пуассона).

 

 

Из (8) и (9) вытекает решение уравнения Пуассона: u 4 (x),

x R3 ,

- фи-

нитная с носителем 0 - ограниченная область в R3. Его решение:

 

 

(10) u

 

 

 

1

 

 

( )d 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эта формула следует из того, что (x) (x ) ( )d 3 . Для строгого доказа-

 

 

 

 

 

 

 

 

0

 

 

тельства следует подставить (10) в уравнение Пуассона и воспользоваться результатом примера (формулой (1)).

Уравнение Гельмгольца.

 

 

 

 

 

(11)

k 2

u F(x), k R - уравнение Гельмгольца.

Рассмотрим функцию V (x,t) e i t u(x)

 

(13)

 

 

Эта функция удовлетворяет уравнению

2V

a2 V f (x) (12)

 

 

 

 

 

 

 

 

 

 

t 2

 

 

Подставим (13) в (12) и обозначим

a2

 

k 2 ,

 

f (x)

F (x) . Тогда мы получим (11).

 

 

 

 

 

 

 

 

 

2

 

 

 

2

Здесь k

a

 

 

2

- волновой вектор.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Фундаментальное решение - k (x, ):

(14) ( k 2 ) k (x, ) (x ).

Факт 1 (фундаментальные решения уравнения Гельмгольца).

Имеются два фундаментальных решения:

Уравнения математической физики, семестр 1. Глава 3 «Элементы теории обобщенных функций».

62

(15) (x, )

1

 

e ik

 

x

 

 

 

 

 

 

k

4

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При k=0 эта функция становится фундаментальным решением оператора Лапласа:

 

 

1

 

 

 

k 0

 

 

 

 

 

 

.

4

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача.

Доказать факт 1, воспользовавшись правилом Лейбница и формулой:

 

 

1

 

 

 

 

 

 

 

 

 

(x ) .

 

 

 

 

 

4

x

 

 

 

 

 

 

 

 

 

 

 

 

Проблема.

Какое из этих двух решений плохое?

Уравнения математической физики, семестр 1. Глава 3 «Элементы теории обобщенных функций».

63

Глава 4 «Метод ВКБ – Маслова для уравнений, содержащих малый параметр при производной»

ВКБ – Вентуель Крамерс Бриллиэн.

Метод ВКБ – для О.Д.У. Метод Маслова – для У.Ч.П.

§1. Метод регулярной теории возмущения (для О.Д.У.)

Постановка задачи.

Дано уравнение:

(1) y 2 (x, ) y 0.

0< <<1 и 2 (x, ) 0

2 (x) 1 (x) 2 2 (x) , 0, x - асимптотический ряд.

Если =0, то y

2 (x) y

0

0

(2)

 

0

0

 

 

 

 

 

 

 

Если 2 (x) const

2 y

0

C

ei 0 x C

e i 0 x .

0

 

0

 

 

 

 

 

Требуется изучить y(x, ) при 0, если все известно о y0(x).

Решение задачи.

Вид решения:

(3)y(x, ) y0 (x) y1 (x) 2 y2 (x)

Здесь y0, y1, y2,… - подлежат определению так, чтобы функция (3) приближенно удовлетворяло уравнению (1).

y y y 2 y - подставим это выражение в (1).

 

 

0

1

2

 

 

 

 

 

(x) y y 2 y 0

(1) y y

2 y

2 (x)

1

(x) 2

2

0

1

2

0

 

 

0

1

2

Соберем все коэффициенты при одинаковых степенях :

(4)

y

 

2 (x) y

 

y 2

(x) y

(x) y

 

O( 2 )

d 2

2 (x, ) y

 

y

 

0

0

0

 

 

2

0

 

 

0

 

 

 

 

 

 

 

1 0

 

 

 

1

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

Обозначим:

 

 

d 2

2

(x, ) L (i

 

 

 

, x, ) - по 0 оператор зависит регулярно.

 

 

 

dx2

 

x

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

( i).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx2

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используем принцип независимых коэффициентов в (4).

 

 

 

 

 

Если L1( y(x, ))=Ox( ),

 

2

(x) y0 0

 

 

 

 

 

 

 

 

 

то y0 0

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

Если L1( y(x, ))=Ox(

 

), то y1 0

 

 

(x) y1 y0 (x) 1 (x)

(6)

 

 

 

 

 

И т.д.

Уравнения математической физики, семестр 1. Глава 4 «Метод ВКБ – Маслова для урав- нений, содержащих малый параметр при производной».

64

§2. Метод сингулярной теории возмущения для О.Д.У. (метод ВКБ)

Постановка задачи.

 

2

y

2

 

 

 

 

 

 

2

 

2

 

 

 

 

(1)

 

 

(x) y 0,

1 L2

( i

 

, x)

i

 

 

 

 

(x) : L2

( i

 

, x) y(x, ) 0

 

 

x

 

 

x

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

lim y(x, ) y0 (x).

0

Проблема: как ведет себя y=y(x, ) при 0?

Замечание 1.

К (1) приводит уравнение из предыдущего параграфа, если 2 (x, ) 2 (x ). Действительно, введем новую переменную: x x . Тогда:

 

d

( 2 (x ))

 

d

2 (x ) O

( ) ( 0) , т.е. мы имеем уравнение (1) с медленно

 

 

 

 

dx

 

 

 

 

 

 

 

dx

x

 

 

 

 

 

 

 

 

 

 

 

меняющейся частотой. Действительно:

 

d 2

2

 

d 2

 

 

уравнение из предыдущего пункта принимает вид:

 

dx2

 

dx 2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

d

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

(x ) y(x , ) 0

, что есть уравнение (1).

 

 

dx 2

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание 2.

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

("2") y(x, ) y0 (x) y1 (x)

Подставим (2) в (1), соберем члены при k и приравняем их к нулю. Получим:

0 : 2 (x) y0 (x) 0,

1 : 2 (x) y1 (x) 0, y 0 нет такого анзатца.

Замечание 3.

Пусть 2 (x) const 2 . Тогда решение уравнения (1) легко находится:

 

 

i

x

 

e

i

x

 

 

 

 

 

 

 

 

y(x, ) C

e

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim y(x, ) , т.к. lim sin

a

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

S ( x)

 

 

i

S ( x)

 

Запишем решение в следующем виде: y(x, ) C

 

 

C

 

 

, где S (x) (x)

e

0

e

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

Идея решения.

Анзатц приближенного решения (1) будем искать в виде искаженной плоской волны, т.е. в следующем виде:

Уравнения математической физики, семестр 1. Глава 4 «Метод ВКБ – Маслова для урав- нений, содержащих малый параметр при производной».

 

 

 

 

 

 

 

 

 

 

 

i

S ( x) (x);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

y(x, ) e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь S(x), (x) - функции, которые подлежат определению.

 

 

Подставим (2) в (1). Для этого понадобится формула Лейбница:

 

 

 

 

d

 

i

S ( x)

 

 

 

 

 

i

 

 

 

 

i

S ( x)

 

i

 

S ( x)

 

d

 

i

S ( x)

S ( i )

 

 

 

 

 

 

 

 

S (x)e

 

 

 

 

 

 

i

 

 

e

 

 

(x)

i

 

(x) e

 

 

 

i

 

 

e

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

d 2

 

 

 

i

 

 

 

(3)

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S ( x)

d

 

 

 

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

e

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3)

i

S ( x) S (S (i ) ) (i )(S S (i ) )

 

 

 

 

e

 

 

 

 

 

i

S ( x) S 2 ( i ) 2S S ( i )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

65

(3)

 

 

 

 

 

 

 

2

 

 

i

 

i

 

2

d

 

 

 

 

e

 

 

S ( x) (x) e

 

S ( x) S 2 i 2S S 2

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание.

 

 

 

 

 

 

 

 

 

 

 

 

 

1

k

 

- волновое число. 0 k , т.е. в данной задаче рассматривается

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

коротковолновое приближение.

Подставим (2) в (1), используя (4), и соберем члены 0, 1, 2.

 

 

i

 

 

 

 

 

 

 

 

 

i

S ( x) 0 S 2

2 i 2S S 2 (x) (5)

L

e

 

S ( x)

 

 

e

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

e

 

 

 

S ( x) O( ) S 2 2 (x) 0

Пусть L

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S 2 2 (x) 0 (i)

 

 

Это уравнение на фазу S.

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

Пусть L

e

 

S ( x) O( 2 ) 2S S 0 (ii)

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Утверждение 1.

Если S(x) и (x) удовлетворяют (i) и (ii), то:

 

d 2

 

 

 

 

i

S ( x)

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

(x)

 

e

(x) O( 2 ).

 

 

 

 

 

2

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Утверждение 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

1) S(x), удовлетворяющее (i), имеет вид: S (x) (x)dx

 

 

 

 

 

 

 

 

 

x0

 

 

2) (x), удовлетворяющее (ii), имеет вид: (x)

 

C

 

 

 

 

 

 

S (x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6)

 

 

 

C

 

(7)

 

 

 

 

(x)

 

 

 

 

 

 

 

Уравнения математической физики, семестр 1. Глава 4 «Метод ВКБ – Маслова для урав- нений, содержащих малый параметр при производной».

66

Обоснование формулы (7).

 

 

 

 

 

1

 

 

 

 

 

 

x

Сделаем замену: (x)

 

 

 

 

f (x);

S (x)

S(x) (x)dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S (x)

 

 

 

 

x0

Пусть

 

S (x)

 

(x) 0,

 

 

S (x)

 

- гладкая функция.

 

 

 

 

Подставим в (ii):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

f

 

 

 

 

S (x)

 

 

 

d

 

 

 

 

 

 

 

 

 

S

 

 

 

 

0

2S

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S (x)

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

12

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

2

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S (x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S (x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

d

 

 

S (x)

 

 

 

S

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dfdx 0 f const , что и требовалось.

Теорема (о решении задачи).

Пусть (x)>0.

Тогда ВКБ-решение уравнения (1) имеет вид:

2S (x)

 

 

df

0

 

 

 

 

 

 

S (x)

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

i

x

 

 

 

i

x

 

 

 

( x)dx

 

 

 

( x)dx

 

 

 

 

 

 

 

 

y(x, )

C e x0

 

C e x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x)

 

 

(x)

Эта функция удовлетворяет уравнению (1) с точностью до O( 2 ) ( 0).

Эта теорема автоматически вытекает из утверждений 1 и 2.

§3. Метод ВКБ-Маслова для уравнения Шредингера

Постановка задачи.

Вспомним уравнение Шредингера и поставим для него задачу Коши:

 

 

 

 

h2

 

 

ih

 

t

 

2m V (x,t)

(1)

 

 

 

 

i

 

 

 

 

 

 

 

S0 ( x)

 

 

 

 

e h

(x)

(2)

 

 

 

 

 

 

t 0

 

 

0

 

 

 

 

 

 

 

Cчитается, что h<<1 (т.е. здесь роль играет h).

Замечание 1 (переход к классической механике).

Если h 0, то данная задача переходит в задачу классической механики.

(1')

H (x, p)

p2

V (x,t)

2m

 

 

 

 

 

 

 

 

p xV ,

 

 

 

 

 

 

 

 

 

(2')

 

 

p

 

m x xV (x,t) - уравнение Ньютона.

 

 

 

 

 

x

 

.

 

 

 

 

 

 

 

 

 

 

m

 

 

 

Уравнения математической физики, семестр 1. Глава 4 «Метод ВКБ – Маслова для урав- нений, содержащих малый параметр при производной».

67

Идея решения.

Решение уравнения (1) ищем в виде:

i S ( x,t )

(3) (x,t) e h (x,t)

Распишем вторую частную производную по x этого решения. Результат будет аналогичен полученному в предыдущем параграфе.

 

 

 

 

 

 

 

 

 

3

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h2 h2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x j

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

i

S ( x,t )

 

 

 

 

 

 

 

i

S ( x,t )

 

 

 

 

 

 

S

2

 

S

 

 

 

2

S

 

 

2

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

h

 

 

 

 

h

 

 

 

(x,t)

e h

 

 

( 1)

 

 

 

ih 2

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

, j 1,2,3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x j

 

 

 

x j x j

 

 

 

 

 

2

 

 

 

 

 

x j

2

 

 

 

 

x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

2

 

 

 

2

 

 

i

S

 

 

i

S

1

S

2

 

 

 

2 S

 

 

1

2

S

 

 

 

h

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

(4)

 

 

 

 

 

 

 

 

 

 

 

h

 

e

h

 

 

 

 

 

 

 

 

(ih)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

e

 

 

 

 

 

 

2m

 

x j

 

 

2m x j x j

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

2m x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2m x j

 

 

2m x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем частную производную по t функции (3):

 

 

i

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

S

S

 

 

e h

 

 

 

 

(4') ih

 

 

e h

 

(ih)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

t

 

 

t

 

 

 

 

 

 

 

 

 

 

Представим (1) в операторной форме:

 

ˆ

 

 

 

 

 

 

 

 

h2

 

 

 

 

(1) L

 

ih

 

 

 

 

 

V (x,t) 0

 

1

 

 

 

 

 

t

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставим (3) в (1), используя (4) и (4’) и соберем члены при степенях h. В результате получим:

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5)

 

 

ˆ

 

 

h

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

S

 

 

0 S

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

1

 

 

 

h2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

e h

 

h

 

 

t

 

 

 

 

S

V (x,t)

(ih)

 

 

t

 

 

 

S,

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

2m

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

S ( x,t )

 

 

 

 

 

 

 

 

 

 

 

 

S

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Пусть

 

ˆ

e

h

 

 

(x,t)

 

O(h) (h 0)

 

 

 

 

 

 

S

V (x,t) 0

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(i)

 

 

S

 

 

1

 

 

S 2

V (x,t) 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S(x,t) ?

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

S ( x,t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть

 

ˆ

e

h

 

 

(x,t)

 

O(h

2

) (h 0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ii)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 S, 1 S 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

m

 

 

 

 

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Утверждение 1 (о решении уравнения Шредингера).

Если S и - гладкие решения (i) и (ii), то функция (3) является приближенным с точностью до O(h2) решением уравнения (1), т.е.

Уравнения математической физики, семестр 1. Глава 4 «Метод ВКБ – Маслова для урав- нений, содержащих малый параметр при производной».

68

 

 

i

S

 

h2

i

S

 

 

 

L

 

 

O(h

 

) (h 0)

e

 

 

 

e

h

 

2

ˆ

 

h

 

 

 

 

 

 

 

1

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство очевидно в силу конструкции.

Решение задачи.

Из начальных условий для (1) следуют начальные условия (i’) и (ii’) для уравнений

(i) и (ii) соответственно:

S(x,t)

(x,t)

t 0

S0 (x)

(i )

t 0

0 (x)

(ii )

 

(i)

 

 

 

 

 

 

S

 

 

 

-

задача

Коши для уравнения Гамильтона Якоби

 

H (x, S,t) 0 , где

(i )

 

 

 

 

 

t

 

H (x, p,t)

1

 

 

2

V (x,t) .

 

 

 

 

p

 

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

Ее решение определяется системой Гамильтона-Якоби:

p

 

V

(x,t),

p

 

 

 

 

 

 

x

t 0

x

S(x

0

),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

p

 

 

x

 

 

.

 

 

 

 

 

 

 

 

 

t 0

0

 

 

 

 

x

.

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

(ii) линейное уравнение 1-ого порядка. Способ его решения подробно был рассмотрен в §1 и §3 главы 1.

Утверждение 2 (о переходе к уравнению неразрывности).

Замена (x,t) 2 (x,t) сводит задачу (ii)-(ii’) к следующей задаче:

(6)

 

 

 

 

 

 

 

 

div V (x,t) 0,

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 (x)

2

(x).

 

 

 

(7)

 

t 0

0

 

 

Уравнение (6) – уравнение неразрывности.

 

1

 

Здесь V (x,t)

 

S - потенциальное поле.

m

 

 

 

 

rotV (x,t) 0.

 

 

Замечание 2 (решение уравнения неразрывности).

 

 

0

(x

0

)

 

 

 

 

Решение этой задачи: (x,t)

 

 

 

 

 

x x (t) (8)

 

 

 

 

 

J (x0 ,t)

 

0 0

 

Здесь x0(x,t) определяется из уравнения: x=X(x0,t). А X и J определяются как решение задачи:

Уравнения математической физики, семестр 1. Глава 4 «Метод ВКБ – Маслова для урав- нений, содержащих малый параметр при производной».

69

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V (x,t)

 

 

 

 

 

S,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 0 x0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(9)

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

J (x

 

 

 

 

 

 

 

 

 

 

 

 

 

,t) det

 

 

 

 

 

(x ,t) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

x0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание 3 (о доказательстве утверждения 2).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оно доказывается на основе формулы: div V

V

, divV.

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

divV

div

 

S

 

 

 

 

 

 

 

S.

 

 

 

 

 

2m

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Утверждение 3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

(x

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x,t)

 

 

(x,t)

 

 

 

 

 

0

 

 

 

 

 

 

x (x,t) (8')

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J (x0 ,t)

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это автоматически следует из утверждения 2.

Замечание 4.

Итак, мы получили, что асимптотическое решение (1)-(2):

i

(x,t, h) e h S ( x,t ) (x,t) O(h2 ) (10)

Здесь функции S и определяются формулами (8’), (9) и решениями (i)-(i’) и (ii)-

(ii’).

X(x0,t), x0 R3, определяется из уравнения Ньютона: m x xV (x,t) .

Заключительное замечание ко всему курсу.

Спасибо всем, кто слушал (читал) мои лекции. Желаю всем благополучно сдать экзамены.

Профессор Белов В.В. 20.12.2000г.

Студент 3 курса Ярчук С.В. 3.02.2001г.

Уравнения математической физики, семестр 1. Глава 4 «Метод ВКБ – Маслова для урав- нений, содержащих малый параметр при производной».

70

Предметный указатель

А

алгоритм решения задачи Коши для квазилинейного уравнения, 14

Б

безопасность ядерного реактора, 34

В

векторное поле, 14 волновое уравнение колебаний, 8 волновой вектор, 61

Д

действие по Гамильтону, 22 дельта образная последовательность, 49

З

задача «ГИБДД», 12 задача для струны, 39

задача Коши для квазилинейного уравнения, 13

задача о трубе, 19 задача Штурма-Лиувилля, 32 закон Фурье, 30

И

инвариантность преобразования Фурье,

58

интеграл в смысле главного значения, 48 интегральная поверхность, 14

К

квазилинейное уравнение, 11 классическое преобразование Фурье, 55 компакт Шварца, 41 концентрация, 29 корректно-поставленная задача, 36 коэффициент диффузии, 29

коэффициент диффузии нейтронов, 35 коэффициент теплопроводности, 29 краевое условие III-го рода, 30 краевое условие II-го рода, 30

краевое условие I-ого рода, 30

Л

лагранжева поверхность, 28 лемма Гамильтона, 25 лемма Лебега, 55 лемма Лиувилля, 21

линейная замена переменных в обобщенных функциях, 52

линейное уравнение 1-ого порядка, 10 линейное уравнение бегущей волны, 9

М

матрица Якоби, 21 метод ВКБ, 64 метод Фурье, 32 мираж, 27

Н

начально-краевая или смешанная задача,

30

непрерывная среда, 8 нестационарное уравнение Гамильтона-

Якоби, 22 носитель, 40

носитель - функции, 47

О

обобщенная функция, 46 обобщенно регулярные функции, 46

обобщенное преобразование Фурье, 57 обобщенное уравнение Гамильтона-

Якоби, 27 обобщенные функции медленного роста,

58

общее линейное уравнение, 15

П

первообразная от - функции, 47 переменные Лагранжа, 12 переменные Эйлера, 12 переход к классической механике, 66 плотность вещества, 29

плотность источников вещества, 29 плотность тепловых источников, 29 пористость среды, 29 порядок уравнения, 8

Уравнения математической физики, семестр 1. Глава 4 «Метод ВКБ – Маслова для урав- нений, содержащих малый параметр при производной».

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