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

FA Арсеньев Функ.Ан

.pdf
Скачиваний:
52
Добавлен:
30.05.2015
Размер:
1.69 Mб
Скачать

Обыкновенные дифференциальные операторы.

Пусть

L

d

 

=

d

 

n

 

 

d

 

n 1

 

 

 

 

 

+ a(n 1)

 

 

+ + a0 ; t 2 R1:

(6.80)

dt

dt

 

dt

Дифференциальному оператору (6.80) мы поставим в соответствие по-

лином

L( ) = n + a(n 1) n 1 + + a0:

Напомним

Определение 6.3.2. Функцией Коши K(t) для оператора (6.80) называется решение уравнения

LK(t) = 0;

K(n 1)(0) = 1 ; K(n 2)(0) = = K(0) = 0:

Простое вычисление показывает, что функция Коши может быть вы- числена по формуле

K(t) = 2 i

I

L( ) d ;

1

 

exp( t)

j j=R

где радиус R больше модуля всех корней многочлена L( ).

Утверждение 6.3.1. Решение уравнения

d

L dt E(t) = (t)

дается формулой

E(t) = (t)K(t):

(6.81)

(6.82)

(6.83)

Это утверждение можно проверить непосредстеннным вычислением. Делается это так. По определению производной от распределения, равенство (6.83) означает, что

 

 

+1

 

d

dt =

8( 2 D(R)):

Z1 E(t)L

 

 

dt

+1

 

 

d

 

 

 

 

Z0

K(t)L

 

dt = (0):

 

dt

 

Для доказательства последнего равенства перебрасываем оператор L на K(t) и учитываем подстановку в нуле.

Можно заметить, что формула (6.83) фактически есть лишь другая запись формулы Дюамеля (3.275).

439

Волновое уравнение в размерности 1+1.

Найдем решение уравнения

 

@2

 

@2

E(x ; t) = (t) (x):

(6.84)

@t2

@x2

Беря от (6.84) преобразование Фурье по переменной x (преобразование

Фурье функции f(x) мы будем обозначать символом fb( )), мы получаем

уравнение

d2Eb( ; t) + 2Eb( ; t) = (t): dt2

Из формулы (6.84) следует, что

 

Eb( ; t) = (t)

sin( t)

; 2 R1:

 

 

 

 

 

Следовательно,

 

 

 

1

 

 

 

d :

 

E(x ; t) = 2( ) Z

exp(i x)

 

 

t

 

 

 

sin( t)

 

 

 

 

1

 

 

 

 

 

 

При вычислении этого интеграла полезно заметить, что

1

Z

exp( i( ; x)) (a jxj)dx = 2sin(a );

1

поэтому

 

 

 

sin( t)

 

1

 

 

 

F 1

 

=

 

(t jxj):

 

 

 

2

 

Мы получаем

 

 

 

 

 

 

 

 

 

Утверждение 6.3.2. Решение уравнения (6.84) дается формулой

 

 

E(x ; t) =

1

(t jxj):

(6.85)

 

 

 

 

2

Волновое уравнение в размерности 1+3.

 

Найдем решение уравнения

 

 

 

 

 

 

 

 

@2

x E(x ; t) = (t) (x) ; t 2 R1 ; x 2 R3:

(6.86)

@t2

440

Беря от (6.86) преобразование Фурье по переменной x, мы получаем уравнение

d2Eb( ; t) + 2Eb( ; t) = (t) ; t 2 R1 ; 2 R3: dt2

Следовательно,

Eb( ; t) = (t)sin(j jt) ; t 2 R1 ; 2 R3:

j j

Применяя формулу обращения преобразования Фурье,получаем:

8(

 

 

 

0):

 

(

 

) =

 

2

3

 

 

 

 

 

exp( (

 

))

 

jjjj

 

 

=

 

 

 

 

 

R!1 Z

 

 

t)

 

 

 

t >

 

 

 

E

x ; t

 

 

 

 

1

 

 

lim

 

 

 

i

x

 

sin(

d

 

 

 

 

 

 

 

 

 

 

 

 

 

j j R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

!1 Z

 

Z

 

 

 

 

 

j j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 Rlim

 

0

 

 

exp(ir x cos( )) sin( )d 1sin(rt)rdr =

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

@

0

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

1

R!1 Z0

2 sin( j

 

j) sin(

)

 

 

=

 

 

 

 

 

 

 

 

2

 

jxj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

 

 

r

x

 

 

 

rt

dr

 

 

 

 

 

 

 

 

 

 

1

 

2

 

1

R!1

sin((t

 

x )R)

 

 

sin((t +

x )R)

 

 

 

 

2

 

 

x

 

(t

jx j)

 

 

 

 

(t +

jx j)

 

 

=

 

 

 

 

 

1

 

 

 

 

j

j

lim

 

 

 

 

j

j

 

 

 

 

 

j j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(jxj t):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 jxj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, справедливо

Утверждение 6.3.3. Решение уравнения (6.86) есть функция

E(x ; t) = 1 (jxj t) ; x 2 R3:

4 jxj

Волновое уравнение в размерности 1+2.

Найдем решение уравнения

@2

@t2 x E(x ; t) = (t) (x) ; t 2 R1 ; x 2 R2:

Применяя формулу (6.87) к основной функции вида

(x1 ; x2 ; x3 ; t) = (x1 ; x2 ; t){(x3 ; 0 ; R ; 1) ; R ! 1;

(6.87)

(6.88)

441

где функция {(x3 ; 0 ; R ; 1) задается формулой (6.18), мы получаем, что решение уравнения есть функция

1p

 

 

 

 

t

Z

(

x2 + z2

t

 

 

 

 

 

 

 

E(x ; t) =

( )

 

p

 

 

)

dz =

 

 

 

 

4

 

 

 

 

 

 

 

x2 + z2

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

2

 

x2

2 t2

x 2

1

2

 

2 R

 

 

2 Z

 

 

 

 

(t)

(

t)

d =

 

(t

jxj)

; x = (x

; x

)

 

2:

 

 

jxj

p

 

 

 

 

p

j j

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, справедливо

Утверждение 6.3.4. Решение уравнения (6.88) есть функция

E(x ; t) =

(t jxj)

; x

2

R

2

:

 

 

 

 

2 pt2 jxj2

 

 

Уравнение теплопроводности.

Найдем решение уравнения

@

x

E(x ; t) = (t) (x) ; t 2 R1 ; x 2 Rd:

@t

Совершая преобразование Фурье по x, мы получаем уравнение

dEb( ; t) + 2Eb( ; t) = (t); dt

Eb( ; t) = (t) exp( 2t):

Совершая обратное преобразование Фурье, мы получаем

Утверждение 6.3.5. Решение уравнения (6.90) есть функция

E(x ; t) = (t)(4 t) d=2 exp( x2=4t):

Уравнение Лапласа.

Нам нужно найти решение уравнения

( x)E(x) = (x) ; x 2 Rd ; d 2:

(6.89)

(6.90)

(6.91)

(6.92)

(6.93)

442

Сначала приведем фомальные выкладки, аналогичные тем, которые делаются в аналогичных случаях физиками. Имеем:

2E

) = 1 ;

 

(d

3) : E(x) = F 1(

2)(x) =

 

b1(

1

 

8

2

 

1

 

 

1

 

 

j j

 

2

 

 

 

F (Z0

exp( t )dt)(x) =

Z0

F

(exp( t

)(x)dt =

1

 

 

 

 

 

 

 

1

 

1

 

d=2

 

1

 

d 2

Z0

(4 t) d=2 exp(jxj2=4t)dt =

 

 

 

 

 

 

 

 

 

(d=2 1);

4

 

 

 

jxj

 

 

 

 

 

 

1

4 (exp(jxj2=4t) exp(jyj2=4t))dt =

d = 2 : E(x) E(y) = Z0

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

21 (ln jxj ln jyj):

Для обоснования этих выкладок используем теорию полугрупп. В пространстве L2(Rd ; d ) функция

Tb(t) : ( ) 7!exp( 2t) ( )

(6.94)

есть полугруппа класса C0 (проверку этого утверждения мы оставляем читателю в качестве упражнения). Так как преобразование Фурье F есть взаимно однозначное унитарное (с точностью до множителя (2 ) d) ïðå- образование пространства L2(Rd ; dx) в пространство L2(Rd ; d ), функция

 

 

 

 

T (t) := F 1T (t)F

 

(6.95)

 

L

(R

 

; dx) åñòü

b

 

C0. В пространстве

в пространстве

2

 

d

 

 

полугруппа класса

 

 

L2(Rd ; dx) оператор T (t) есть интегральный оператор

 

 

T (t)

 

(x) = Z

G(t ; x ; y)

(y)dy;

 

(6.96)

 

G(t ; x ; y) = (4 t) d=2 exp( (x y)2=4t):

(6.97)

Инфинитезимальный оператор полугруппы T (t) в пространстве L2(Rd ; dx) есть оператор Лапласа:

x = F 1( 2)F:

Пусть g(x ; y) -интегральное ядро оператора ( x) 1. Из формулы (3.10.8)

443

и (6.96) следует, что

 

 

 

 

 

 

 

1

 

 

 

 

 

 

g(x ; y) = Z0

(4 t) d=2 exp( (x y)2=4t)dt =

 

1

 

1

 

d=2

 

 

1

 

 

 

d 2

 

 

 

 

 

 

 

 

 

 

(d=2 1) ; d 3:

(6.98)

4

 

 

x

 

y

 

 

 

 

 

 

 

 

 

j

 

j

 

 

Ïðè d = 2 необходима регуляризация интеграла:

g(x ; y) g(x ; z) =

 

 

 

 

 

 

 

 

 

 

 

1

lim

1exp(

t)

exp(

 

(x

 

y)2=4t) exp( (x z)2=4t)

dt =

 

 

 

 

 

 

 

 

 

 

 

4

!+0 Z0

 

 

 

 

 

 

 

 

 

t

 

1

(ln jx zj ln jx yj) ; x ; y ; z 2 R2:

(6.99)

 

 

 

2

Мы доказали

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Утверждение 6.3.6. Решение уравнения (6.93) дается формулой

 

 

E(x) =

8

4

 

1

 

jxj

 

d 2

(d=2 1) ; d 3;

(6.100)

 

 

 

 

 

1

 

 

1

d=2

1

 

 

 

 

 

 

 

<

 

 

 

 

ln x

; d = 2:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

j j

 

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

6.4Пространства Соболева.

6.4.1Преобразование Фурье-Планшереля.

Напомним определение преобразования Фурье функций из пространства Шварца:

Z

8( 2 S(Rd)) : F ( )( ) = exp( i(x ; )) (x)dx; (6.101)

формулу обратного преобразования Фурье:

Z

8( 2 S(Rd)) : (x) = (2 ) d exp(i(x ; ))F ( )( )d ; (6.102)

и равенство Парсеваля:

Z Z

8( 2 S(Rd)) : j (x)j2dx = (2 ) d jF ( )( )j2d : (6.103)

444

(6.106)

Интегрирование везде ведется по пространству Rd.

В дальнешем нам будет удобно использовать введенное ранее обозна- чение

b( ) := F ( )( ):

(6.104)

Ясно, что формулой (6.101) можно определить преобразование Фурье только для тех функций , для которых интеграл в (6.101) сходится.

Наша ближайшая цель состоит в том, чтобы распространить определение преобразования Фурье на все функции из L2(Rd ; dx).

Лемма 6.4.1. Пространство Шварца S(Rd) плотно в L2(Rd ; dx) по метрике пространства L2(Rd ; dx):

Cl(S(Rd)) = L2(Rd ; dx):

Доказательство. Пусть

L2(Rd ; dx) = Cl(S(Rd)) H0:

Докажем, что

H0 = 0:

Åñëè f0 2 H0, то функция f0 ортогональна любой функции из S(Rd). В частности, функция f0 ортогональна любой функции типа гриб , по-

этому функция f0 ортогональна характеристической функции любого параллеллипипеда K Rd:

Z

8K : f0(x)dx = 0:

K

Следовательно, функция f0 ортогональна любой ступенчатой функции и поэтому равна нулю почти всюду. Лемма доказана.

Следствие 6.4.1. Åñëè 2 L2(Rd ; dx), то существует такая последовательность n 2 S(Rd), ÷òî

k n j L2(Rd ; dx)k ! 0 ; n ! 1:

(6.105)

Если последовательность f ng удовлетворяет условию (6.105), то она фундаметальна в L2(Rd ; dx). Из равенства Парсеваля следует, что после-

довательность преобразований Фурье fbng фудаментальна в L2(Rd ; d ). Следовательно,

9b( ) : kb cn j L2(Rd ; d )k ! 0 ; n ! 1:

445

Определение 6.4.1. Если функции 2 L2(Rd ; dx) è b 2 L2(Rd ; d )

связаны равенствами (6.105)-(6.106), то мы говорим, что функция b 2

L2(Rd ; d ) есть преобразование Фурье-Планшереля функции 2 L2(Rd ; dx):

 

 

lim ( );

(6.107)

b(

 

) := n!1 bn

 

где последовательность f ng удовлетворяет условию (6.105).

Ясно, что преобразование Фурье-Планшереля определено на всем пространстве L2(Rd ; dx), оно преобразует пространство L2(Rd ; dx) â ñåáÿ,

для него справедлива (соответственно обобщенная) формула обращения и формула Парсеваля. В дальнешем мы не будем делать различия между заданным формулой (6.101) преобразованием Фурье и определенным формулой (6.107) преобразованием Фурье-Планшереля.

6.4.2Определение и основные свойства пространств Соболева.

Напомним, что заданная в Rd функция f(x) называется локально инте-

грируемой:

f 2 Lloc(Rd);

если для любого компакта K Rd справедливо включение:

f 2 L1(K):

Определение 6.4.2. Распределение f 2 S(Rd)? принадлежит простран- ству Соболева Hs(Rd), åñëè

1. Преобразование Фурье распределения f задается локально интегрируемой функцией:

9(f( ) 2 Lloc(Rd)) ; 8( 2 S(Rd)) :

 

F (b)( ) = ( ( )) =

( ) ( )

 

Z

b

(6.108)

f f F

f d :

2. Задающая преобразование Фурье функция удовлетворяет условию:

kf j Hs(Rd)k2 = Z

jf( )j2

(1 + j j2)sd < 1:

(6.109)

def

b

 

 

 

d

), то входящая

Если распределение f задается функцией f(x) 2 S(R

в (6.108) функция fb( ) есть преобразование Фурье функции f(x), íî â

446

общем случае в определении не требуется, чтобы функция fb( ) обяза-

тельно была бы преобразованием Фурье некоторой функции.

Входящй в определние пространства Соболева Hs(Rd) параметр s может принимать любое значение:

1 < s < 1:

Определение 6.4.3. Заданная в пространстве Rd функция f(x) принадлежит пространству Соболева Hs(Rd), если порожденное функцией f(x) распределение

Z

f : 7! f(x) (x)dx

принадлежит пространству Соболева Hs(Rd).

ßñíî, ÷òî

8(s 0) : Hs(Rd) L2(Rd ; d );

а так как преобразование Фурье переводит пространство L2(Rd) в себя, пространство Hs(Rd) ïðè s 0 может быть отождествлено с подпространством пространства L2(Rd).

Приведем примеры.

1. Пусть f = -функция. Преобразование Фурье распределения f задается функцией fb( ) 1. Имеем:

Z

(1 + j j2)s 1d < 1 ; åñëè s < d=2:

Следовательно,

(x) 2 Hs(Rd) ïðè s < d=2:

2. Пусть f(x) = (1 jxj) ; x 2 R1. Преобразование Фурье функции f(x) есть функция

 

1

exp( ix ) (1 jxj)dx = 2 :

 

Z

 

 

 

 

 

 

 

sin

 

1

 

 

 

 

 

 

 

Имеем:

 

 

 

 

 

 

 

 

Z

(1 + j j2)s4

sin

2

 

 

 

d < 1 ; åñëè s < 1=2:

 

Следовательно,

(1 jxj) 2 Hs(R1) ïðè s < 1=2:

447

Теорема 6.4.1. 1. Пространство Соболева Hs(Rd) есть гильбертово пространство со скалярным произведением

< f ; g >=

Z

f( ) g( )(1 + j j2)sd :

(6.110)

2. Распределения, задаваемые

 

b

b

f(x)

 

 

S(

R

)

H

 

(

R

 

)

 

функциями

 

 

 

 

2

 

 

d , плотны в

 

s

 

d

 

по метрике пространства Hs(Rd).

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство. Для доказательства первого утверждения нам до-

 

статочно доказать, что каждая функция f

 

 

) 2

L2

 

 

d ; (1 + 2)sd ) åñòü

 

 

 

 

 

(

 

 

 

(R. Для доказательстваj j

преобразование Фурье распределения из fb2 S(Rd)?

 

 

 

 

 

 

 

этого утверждения заметим, что, как следует из (6.108), распределение f 2 Hs(Rd) действует на основную функцию по правилу:

f( ) = Z

f( )F 1

( )( )d ;

(6.111)

 

b

 

 

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

8(f 2 Hs(Rd) ; 2 S(Rd)) :

Z 1=2

jfb( )j2(1 + j j2)sd

S

! 0 ; ! 0:

jf( )j

Z 1=2 jF 1( )( )j2(1 + j j2) sd

(6.112)

Из этого неравенства следуют два утверждения: во-первых, что каждая функция fb( ) 2 L2(Rd ; (1 + j j2)sd ) задает по правилу (6.111) линей-

ный непрерывный функционал на пространстве S(Rd), во-вторых, что

из сходимости последовательности распределений в метрике пространства Hs(Rd) следует сходимость в топологии пространства S(Rd)?.

Переходим к доказательству второго утверждения теоремы.

 

Пусть f( )

 

L2

(Rd ; (1 +

 

2)sd ) è fn 2 S(Rd) такая последователь-

ность, что b

2

k

f

n

f

j

Lj2(j

d ; (1 +

j

 

2)sd )

0:

 

Тогда

 

 

 

 

R

j

 

k !

 

 

F (fn)( ) =

Z

fn( ) ( )d ! Z

f( ) ( )d ;

 

à òàê êàê â ñèëó

 

 

 

 

 

 

 

b

 

 

 

 

b

 

 

 

 

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

 

 

 

 

 

 

lim F (fn)( ) = F ( lim fn)( );

 

 

 

 

 

 

 

n!1

 

 

 

 

n!1

 

 

 

то функция f( ) есть преобразование Фурье распределения lim f .

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n!1

n

448

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