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

1

§10 Гладкие решения задачи Коши для уравнения теплопроводности

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

 

 

 

 

 

 

ut = 4xu; (t; x) 2 G;

(1)

G = f(t; x); tP

0

 

 

½ ujt=0 = '(x); x 2 Rn:

 

2 g.

 

 

 

 

 

 

 

n

2

 

 

 

 

 

 

 

 

 

Здесь 4x = k=1

@

- оператор Лапласа,

 

@xk2

 

>

; x

Rn

 

 

 

 

 

 

 

Кроме того, введем следующие обозначения:

 

G

T = f(t; x); 0 < t < T

<

1

; x

2

Rn

g,

 

 

n

 

 

 

 

G±= f(t; x); t ¸ 0; x 2 R g,

n

g.

 

 

 

G = f(t; x); t ¸ ± > 0; x 2 R

 

 

 

В этом параграфе мы докажем следующую теорему

Теорема.

Для любой функции '(x) 2 L2(Rn) решение задачи (1) в области G дается

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

 

 

 

RZ

 

 

 

 

1

 

n

 

x¡yj2

 

 

 

2

 

 

 

u(t; x) = µ

4¼t

n

e¡j

4t '(y)dy; u = K'; K : ' ! u

(2)

при этом:

 

 

 

 

 

а) функция u(t; x) 2 C1(G),

 

n

),

 

б) при любом t ¸ 0: функция u(t; x) 2 L2(R

 

в) функция u(t; x) удовлетворяет начальному условию в следующем смысле:

½(u ¡ ') ! 0 при t ! +0 (в L2(Rn)):

Построенное решение является единственным в данном классе функций и непрерывно зависит от начальных данных.

Замечание к формулировке теоремы.

1) Можно показать, что теорема справедлива и тогда, когда '(x) 2 L2;k(Rn), k < 0, где (см. x8) L2;k(Rn) - весовое L2 пространство с нормой

 

 

 

 

jj'jjL2;k(Rn) = µ n (1 + jxj2)kj'(x)j2dx

1

 

 

 

 

 

2

; k < 0:

 

 

 

 

 

 

 

 

 

RZ

 

 

Решение задачи (1) вновь задается формулой (2).

 

 

Пункты: б) 8t ¸ 0: u(t; x) 2 L2;k(Rn),

 

 

в) jj

u

 

'

n

0

при

t

+0

 

 

¡

jjL2;k(R )

 

 

 

 

! n

)

 

!n

.

 

 

2) Если '(x) 2 C(R

L2(R

) - непрерывная функция, то в формулировке

теоремы пункт в)

можно переписать так:

 

 

 

 

T

 

 

 

 

 

в) функция u(t; x) 2 C(G) и удовлетворяет начальному условию в обычном

смысле:

lim u(t; x) = '(x0):

t!+0 x!x0

В этом легко убедиться, если в формуле (2) сделать замену

 

 

 

z =

y ¡ x

;

 

 

 

 

 

 

 

 

 

 

 

2pt

 

 

 

 

 

 

 

 

 

 

при этом формула (2) перепишется так:

 

 

 

 

u(t; x) = µ¼

n

Z

e¡jzj

'(x + 2zpt)dz:

(20)

 

 

1

 

2

 

2

 

 

 

 

 

 

 

 

 

 

Rn

2

Переходя к пределу под знаком интеграла (почему?) получаем:

x! x0

µ¼

n

Zn

lim u(t; x) =

1

 

2

e¡jzj2 '(x0)dz = '(x0);

 

 

 

t +0

 

 

 

 

!

 

 

 

R

поскольку 1 = ¡¼1 ¢n2 R e¡jzj2 dz.

Rn

Доказательство теоремы:

Оно состоит из нескольких частей.

I) получение априорной оценки и формулы (2). Запишем сначала формальное решение задачи (1) (см. x6):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

tk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xk

 

 

 

 

 

4xk'(x):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(t; x) = et4x '(x) =

=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее, возьмем '(») 2 C01(Rn), тогда

 

'(») = '(x)

2kJ

. Поэтому

 

F

 

 

 

P

 

 

k

 

 

R

 

 

 

2¼i(x;»)

 

 

 

 

 

 

 

 

 

b

 

P

 

 

R

2¼i(x;»)

 

 

 

 

 

1btk

 

µRn e

 

 

 

 

 

 

 

 

 

 

1

tk

 

 

 

 

 

 

u(t; x) = k=0 k! 4x

 

e

 

 

 

'(»)¶ = k=0 k! Rn

4x(e

 

 

 

 

 

 

)'(»)=

 

 

 

R P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

2

 

 

 

R

 

n

 

¼2 » 2

 

 

 

 

 

 

 

b

 

 

 

1 tk( 4¼2 »

j

2)k

 

 

 

 

 

 

 

 

'(»)= Rn e

 

 

 

e¡

 

2

 

 

2

'(»)d»:

 

 

 

= Rn µk=0

 

¡

 

k! j

 

 

 

 

 

 

 

 

 

x;»)

 

 

 

 

 

j j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2¼i(b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2¼i(x;») 4¼

t »

 

 

 

b

 

Поскольку (см. формулу (4) из x9 F e¡ajxj

 

= (¼a ) 2 e¡

aj

j

 

, a =

1

)

 

 

 

 

4t

 

 

 

 

 

F ·µ

 

4¼t

n

e¡j

4jt

¸ = e¡4¼

 

tj»j

 

 

 

 

 

для 8t > 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

 

 

x 2

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

h(») = f(»)g(») = e¡4¼2tj»j2 '(»);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

n

 

jxj2

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h(x)

 

 

функций f(x) = (

1

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 e¡

4t , g(x) = '(x)

 

 

 

 

 

 

 

где

 

- свертка

 

 

 

 

 

b

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

4¼t

 

 

 

 

 

µ

4¼t

 

Zn

e¡j

, то

 

 

(2)

 

 

u(t; x) = Zn

e2¼i(x;»)h(»)= h(x) =

n

4t '(y)dy;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

 

 

 

 

x¡yj2

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

при этом

 

 

 

 

F u(t; x) = u(t; ») = e¡4¼2tj»j2 '(»):

 

 

 

 

 

 

 

 

 

 

(20)

 

 

 

8 ¸

 

 

x

 

 

 

 

 

 

 

 

Ясно, что для

0: u(t; x)

2 Jb

 

(см. (2 ) и

 

 

 

 

 

 

 

b

 

 

 

9). В силу равенства

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

лемму из

 

 

Парсеваля (3) из x9 имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

n ju(t; x)j2dx =

 

n ju(t; »)j2=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

RR

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

RRn

 

8¼2Rt » 2

jb

 

2

 

 

 

 

·

 

 

 

 

 

 

 

2

 

 

 

 

 

RRn

j

 

 

 

 

j

2

 

 

 

 

 

 

 

 

 

 

 

 

j j

 

j

 

 

 

 

 

 

RRn jb

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= e¡

 

 

 

 

 

'(»)

 

 

 

 

 

 

'b(») =

 

 

'(x) dx

 

 

 

 

 

 

jju(t)jjL2

2(Rn) · jj'jjL2

2(Rn)

 

 

 

для 8t > 0:

 

 

 

 

 

 

 

(3)

Из (3) получаем также

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z0T dtRZn

ju(t; x)j2dx = Z0T jjujjL2 2(Rn)dt = jjujjL2 2(GT ) · T jj'jjL2 2(Rn):

(30)

(30) и является искомой априорной оценкой, которую мы получили при условии, что '(x) 2 J ('b(») 2 C01(Rn)). Решение задачи (1) в этом случае дается формулой (2).

3

II) доказательство пунктов а),б),в). Пусть теперь '(x) 2 L2(Rn), тогда существует функция 'b(») = F '(x) 2 L2(Rn) и существует последовательность f'bk(»)g, k = 1; 2; : : :, 'bk(») 2 C01(Rn) такая, что

 

 

jj'bk ¡ 'bjjL2(Rn) ! 0; k ! 1:

 

Записывая (30) для 8k:

 

 

 

jjukjjL2 2(GT ) · T jj'kjjL2 2(Rn);

(300)

где 'k(x) =

 

'bk(»), причем

 

F

 

jj'k ¡ 'jj2L2(Rn) = jj'k ¡ 'jj2L2(Rn) ! 0; k ! 1:

 

 

 

 

RZ

 

 

 

1b

n

b

 

jx¡yj2

 

2

 

uk(t; x) = µ

 

 

 

e¡

4t 'k(y)dy;

4¼t

 

n

и переходя к пределу при k ! 1 в (300) мы получаем, что последовательность fuk(t; x)g такая, что

jjuk ¡ ujj2L2(GT ) ! 0; k ! 1;

 

1

 

n

 

jx¡yj2

 

где u(t; x) =

 

 

2

e¡

4t

'(y)dy L2(GT ), т.е. функция u(t; x) определена

¼t

почти всюду¡в4 G¢T . RRn

 

 

 

2

Из (3) следует, что для 8t ¸ 0: u(t; x) 2 L2(Rn).

Заметим, что так определенная функция иногда называется обобщенным решением задачи (1) при '(x) 2 L2(Rn). Покажем, что на самом деле формула

(2) переводит функцию '(x)

 

2 L2(Rn)

в функцию u(t; x)

2 C1(G). Снова

доказываем это свойство сначала для функций из плотного множества.

Пусть '(») 2 C01(Rn), '(x) =

F

'(») 2 J ; ®0, ® - любые, D0 =

@

 

 

 

 

 

 

@t

 

 

 

 

 

b

D

®0

D

®

u

 

 

±

) =

sup

 

sup

 

D

®0

D

®

u =

 

sup

D

®0

D

®

u

 

 

(R )

 

 

 

0

x

 

 

 

 

0

 

x

 

0

 

x

 

 

 

 

 

jj

 

 

 

 

 

jjB

 

 

 

t¸±>0 x2R

j

 

 

 

 

 

 

 

 

j

 

 

 

 

jj

 

 

 

 

 

 

 

 

 

jjB

 

 

n

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t¸±>0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

по теореме вложения

 

· C(n) sup jjD0®0 Dx®ujjW2[ n2 ]+1(Rn) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t¸±>0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D®0 D®+¯u

 

2dx

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C n

sup

 

 

 

 

 

 

[ n ]+1 j

t; x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

(

) t¸±>0µZn ¯

 

 

0

 

 

x

 

 

(

 

)j

 

 

 

 

 

=

 

 

 

 

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2¼i»)®+¯ 2

 

F D0®0 u(t; x)

2

 

2

по рав-ву Парсеваля

= C(n) sup

 

 

 

 

¯

 

[ n2 ]+1 j

j

 

=

 

 

 

 

 

 

 

 

 

 

 

t¸±>0µRn

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R j

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2¼i»)®+¯

2(

 

4¼2

 

 

2)2®0 e¡8¼2tj»j2

 

 

 

 

 

 

2

2

 

= C(n) sup

 

 

 

¯

 

[ n2 ]+1 j

 

»

j

'(»)

 

 

·

 

 

t¸±>0µRn

 

 

 

 

 

 

 

 

 

 

j

 

¡

 

 

 

j

j

 

 

 

 

 

 

 

 

 

 

b

j

 

 

 

 

 

 

 

 

 

 

R j

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

· C1(n; ±; ®; ®0)jj'jjL2(Rn):

Если '(x) 2 L2(R ), '(e) =

(n ), тоe9f ke(

)g,

k( ) 2

01(

), f k( )g,

 

 

jjujjBm(G±) · C1jj'jjL2(Rn); C1

= C1

(n; m; ±); 8m:

 

 

 

n

»

F ' x

' » ' » C Rn

' x

'k(x) = F 'k(») J , 'k !

' в L2(R ), k ! 1.

 

 

 

 

 

 

b

 

 

b

 

b

 

 

Запишем (4)2для 'k(»b):

 

 

 

 

 

jjD0®0 Dx®ukjjB(G±) · C1(n; ±; ®; ®0)jj'kjjL2(Rn);

4

 

 

uk

 

Bm(G±) · C1jj'kjjL2(Rn)

где

jj

 

jj

 

 

 

n

e

 

 

 

uk(t; x) = µ

1

2

RZ

e¡j

x¡yj2

 

 

 

 

4t 'k(y)dy;

 

4¼t

 

n

понятно, что fuk(t; x)g фундаментальна в Bm(G±) для любого ±, m, т.е. предельная функция u(t; x) 2 C1(G) и удовлетворяет уравнению ut = 4xu в G (т.к. uk(t; x) - решения этого уравнения, то возможен предельный переход в этом уравнении при k ! 1).

Легко убедиться, что построенная функция удовлетворяет начальному условию в следующем смысле

 

 

 

 

 

 

 

 

 

 

lim

 

u t; x

) ¡

'

x

)jjL2

(R

 

 

 

:

 

 

 

 

 

 

 

 

 

 

t +0 jj

(

 

 

(

 

 

) = 0

 

 

 

 

 

 

 

 

 

 

!

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В самом деле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(t; x)

'(x)

L2

(Rn)

=

 

u(t; »)

 

'(»)

 

L2 2(Rn) = [e¡4¼

jj

 

 

¡

4¼2t »jj2

2

 

1]

2

 

 

jj

 

2

 

¡

b

 

jj

4¼2t »

2

 

RRn

2

 

'(»

=

[e¡

j j

 

 

 

 

'(»)b +

[e¡

 

j j

 

 

1]

 

 

j

R

 

 

¡

 

 

j

b

 

j

 

 

j

jR

 

 

 

 

 

 

¡

 

j

b

 

 

»j·D

 

 

 

 

 

 

 

 

 

 

 

 

» >D

 

 

 

 

 

 

 

 

 

 

 

(5)

2tj»j2 ¡ 1]2j'b(»)j2=

)j2d» · "22 + "22 = "2;

8" > 0; 9± > 0; D0 > 0; t < ±; D > D0:

Теорема единственности и непрерывной зависимости решений уравнения ut = 4xu от начальной функции '(x) следует из оценки (30).

Задача.

Найти решение задачи

ut = 4xu + f(t; x);

(t; x) 2 G;

( )

½ ujt=0 = '(x); x 2 Rn

¤

с помощью принципа Дюамеля (см. x8).

Пусть '(x) 2 C(Rn), f(t; x) 2 C1(G). Тогда задача (¤) имеет классическое решение, непрерывное в G.

Замечание.

Почему задача (1) не рассматривается при t < 0?

½ut = uxx; t < 0; x 2 R1 ujt=0 = '(x)

Существует достаточно гладкое решение u(t; x).

un = u + e¡pne¡n2t sin nx,

 

 

 

 

 

 

 

 

= e¡p

 

¡n2t

 

 

 

j,

j

u

n ¡

u

j

n

j

sin nx

 

 

 

 

¼

 

 

 

 

j sin nxj = 1, x = 2 , n - нечетное.

u

njt=0

= '

(x) = '(x) + e¡p

n

sin nx

 

 

 

 

n

 

 

 

 

 

p

 

 

 

 

 

junjt=0 ¡ 'n(x)j · e¡

 

n ! 0, n ! 1.

Соседние файлы в папке Lecture_10