Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
[ Грикуров, Аленицын ] Обобщённые функции в матфизике.pdf
Скачиваний:
146
Добавлен:
22.08.2013
Размер:
658 Кб
Скачать
(4.4)

4.2ОБОБЩЕННАЯ ЗАДАЧА КОШИ

4.2.1Постановка и решение обобщенной задачи Коши для волнового уравнения

 

состоит в отыскании функции

 

удовлетворяющей при

3

Задачи Коши

u(~x, t),

~x R , t > 0

 

2

 

 

уравнению utt − a

 

u = f(~x, t) è ïðè t = 0 — начальным условиям u|t=0 = u0(~x), ut|t=0 =

u1(~x). В случае обобщенных решений говорить об их значенияхïðè t

= 0 не имеет

смысла, поэтому постановка обобщенной задачи Коши нуждается в переформулировке (ñì. также параграф 1.7.3 ).

Будем рассматривать задачу при всех, в том числе отрицательных, временах и под решением понимать функцию u(~x, t), удовлетворяющую волновому уравнению при t > 0 (ò.å. допускающую два дифференцирования по t ïðè t > 0) и продолженную нулем на область t < 0. Аналогично будем понимать правую часть f(~x, t). Вычислим вторую производную по времени от u(~x, t):

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

(utt, φ) = (u, φtt) =

d~x

 

 

dt uφtt

=

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

dt

 

u + f(~x, t) φ =

 

=

d~x

 

 

 

u0

φt

t=0 + u1

φ

 

t=0 +

a2

(4.3)

Z

 

 

 

 

|

 

 

|

 

Z

 

 

 

 

= a

2

 

 

 

 

 

 

 

 

 

0

 

 

u + f(~x, t) + u0 · δ0(t) + u1 · δ(t), φ .

 

 

Таким образом, мы приходим к следующему определению:

Определение 4.2. Обобщенной задачей Коши называется уравнение

 

2

u(~x, t) = f(~x, t) + u0 · δ0(t) + u1 · δ(t) , f(~x, t)|t<0 ≡ 0 ,

 

∂t2 − a2

(4.4)

понимаемое в смысле обобщенных функций.

Утверждение 4.4. Задача (4.4) имеет единственное решение в классе S0; это решение непрерывно зависит от данных задачи.

Решение задачи дается формулой u = E F, в которой E – фундаментальное решение волнового оператора, F = f(~x, t) + u0 · δ0(t) + u1 · δ(t) – данные задачи. Свертка существует по следующим причинам: 1) носитель фундаментального решения по переменным ~x ограничен (ñì. формулы раздела 4.1.1 ); 2) правая часть F полуфинитна по t. В силу непрерывности свертки (как непрерывного функционала) E Fn → E F, åñëè Fn → F, ò.å. E F непрерывно зависит от F. В частности, ïðè F = 0 имеем E F = 0, что и означает единственность.

4.2.2Запаздывающие потенциалы

Рассмотрим более подробно структуру полученного решения :

u = E F = E f(~x, t) + E (u0(~x) · δ0(t)) + E (u1(~x) · δ(t)) .

42

f(x − ξ, t − τ)dτdξ+
θ(τ)θ(aτ − |ξ|)u0(x − ξ)δ0(t − τ)dτdξ+
θ(τ)θ(aτ − |ξ|)u1(x − ξ)δ(t − τ)dτdξ =

Каждое слагаемое в правой части этой формулы носит название запаздывающего потенциала. Из определения 4.4 видно, что начальные данные можно интерпретировать как

источники, действующие мгновенно. Соответствующие им запаздывающие потенциалы могут быть несколько упрощены.

Действительно, по свойствам свертки

g(~x, t) t δ(t) = g(~x, t) , g(~x, t) t δ0(t) =

∂g(~x, t)

t δ(t) =

∂g(~x, t)

,

 

 

 

∂t

∂t

ãäå t означает свертку по переменной t. Поэтому мы можем переписать решение, вычислив свертку по t:

 

E

 

∂t ~x

0

 

E ~x

1

 

 

4.5)

u =

 

f(~x, t) +

∂E

 

u

(~x) +

 

u

(~x) .

(

 

 

 

 

 

Формула (4.5) и является окончательным решением задачи Коши в общем случае . Ниже мы применим эту формулу в различных размерностях к данным задачи , являющимися (кусочно) непрерывными функциями.

4.2.3Формула Даламбера

Рассмотрим решение задачи Коши в пространстве R1, считая данные задачи локально- интегрируемыми функциями. Не будем пользоваться здесь формулой (4.5), а вычислим свертку непосредственно. С учетом f(x, t)|t<0 получаем

u = E f(x, t) + E (u0(x) · δ0(t)) + E (u1(x) · δ(t)) =

Z Z

= 21a θ(τ)θ(aτ − |ξ|)f(x − ξ, t − τ)dτdξ+

Z Z

+ 21a

Z Z

+ 21a

t

Z Z

= θ(t) 2a

 

 

 

 

0 −aτ

 

 

 

 

 

∞ aτ

 

 

 

 

 

+ 2a

∞ aτ

u0(x − ξ)δ0(t − τ)dτdξ + 2a

u1(x − ξ)δ(t − τ)dτdξ+

 

Z Z

Z Z

 

 

θ(t)

 

 

 

θ(t)

 

 

 

 

 

 

 

 

 

0 −aτ

 

 

 

 

 

0 −aτ

 

 

 

 

= 2a

t x+a(t−τ)

at

u0(x − ξ)dξ +

2a

at

u1(x − ξ)dξ =

Z Z

f(ξ, τ)dτdξ +

2a dt Z

Z

 

θ(t)

 

 

 

 

θ(t) d

 

 

 

 

 

θ(t)

 

 

0 x−a(t−τ)

tx+a(t−τ

= θ(t) Z

1

Z

2a

 

 

 

−at

)

f(ξ, τ)dτdξ + u0(x + at) + u0(x at) 2

−at

x+at

Z

+ 1 u1(ξ)dξ .

2a

0 x−a(t−τ) x−at

Данная выкладка не является строгой ; ее следует уточнить, используя прием, примененный при выводе формулы (4.3).

43

Полученный результат носит название формулы Даламбера . Она дает обобщенное решение задачи Коши. В случае достаточно гладких данных задачи (а именно, f(x, t) C1, u0(x) C2, u1(x) C1) этой формулой выражаются классические решения .

Замечание 4.3. (принцип Дюамеля)

 

 

 

t

.

Запишем первое слагаемое в квадратных скобках как

R

v(x, t; τ)dτ, ãäå v(x, t; τ) =

 

 

x+a(t−τ)

 

 

 

 

0

 

=.

1

x−aR(t−τ) f(ξ, τ)dξ. Функция v(x, t; τ) обладает следующими свойствами:

2a

1) ∂t22 − a2 v(x, t; τ) = 0 τ;

2)v(x, t; τ)|t=τ = 0;

3)vt(x, t; τ)|t=τ = f(x, t);

Данные свойства означают, что неоднородность f(x, t), входящая в правую часть волнового уравнения, может быть учтена посредством вспомогательной задачи Коши 1)—3). То же самое имеет место и в старших размерностях .

Важнейшие из физических следствий формулы Даламбера следующие . 1) Принцип суперпозиции – вклад в решение от различных источников суммируется . 2) Решение u(x, t) при f(x, t) ≡ 0 зависит от комбинаций x + at и x − at; в частности, разрывы решения распространяются по линиям x ± at = const, т.е. вдоль характеристик. 3) В случае данных

задачи вида f(x, t) ≡ 0, u1(x) ≡ 0 и функции u0(x), имеющей ограниченный носитель, любая точка пространства при достаточно больших временах оказывается вне носителя

решения (говорят, что отсутствует диффузия волн; это явление имеет место только в нечетных размерностях). Эти свойства присущи решениям любых гиперболических уравнений.

4.2.4Формулы Кирхгофа и Пуассона

Применим формулу (4.5) â R3. Напомним, что если δSρ (x) – δ-функция, сосредоточенная на сфере радиуса ρ с центром в начале координат (для краткости пишем x вместо ~x, т.д.),

то интегралы вида

R

 

 

 

 

 

R

 

 

с центром в точке

 

R

 

 

 

δSρ (x)φ(x)dx понимаются как

φ(x)dSx.

Соответственно,

δSρ (x −

сделанного

 

R

 

 

 

:

 

 

 

 

 

 

 

 

 

 

x0)φ(x)dx =

φ(x)dSx,

ãäå

Sρ(x0) –

сфера радиуса

ρ

 

 

 

 

x0.

С учетом

 

 

 

 

 

 

 

 

 

 

 

 

Sρ(x0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

замечания получаем

 

 

 

 

 

 

 

 

 

 

 

 

E f(x, t) = θ(t) Z

t

 

 

 

 

 

 

 

 

 

Z

f(ξ, τ)dSξ ;

v(x, t; τ) , v(x, t; τ) = 4πa2(t τ)

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

Sa(t−τ)(x)

 

 

 

(видим, что и здесь выполняется принцип Дюамеля);

 

 

 

 

 

 

 

 

 

 

 

 

E x

u1(x) = θ(t)4πa2t

Z

u1(ξ)dSξ;

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

Sat(x)

Это свойство непосредственно вытекает из того факта , что v(x, t; τ ) фактически зависит от x + at и

x − at.

†Можно показать, что и здесь функция v(x, t; τ ) удовлетворяет принципу Дюамеля , т.е. обладает свойствами 1)–3) замечания 4.3.

44

Sat(x)

∂E

 

x u0(x) = θ(t)

 

1

 

Z

u0(ξ)dSξ .

∂t

∂t

2

t

 

 

4πa

 

 

 

 

 

 

 

Sat(x)

 

 

 

 

 

 

 

 

 

Таким образом, мы приходим к известной формуле Кирхгофа, дающей решение обобщенной задачи Коши в R3:

 

Z

t

 

 

 

 

 

 

 

 

 

u(x, t) = θ(t)

v(x, t; τ)dτ +

tMSat(x)[u0]

+

tMSat(x)[u1]

,

(4.6)

 

 

∂t

 

 

0

 

 

 

.

1

 

 

 

 

 

 

ãäå v(x, t; τ) = (t − τ)MSa(t−τ)(x)[f|t−τ ], MSat(x)[g] =

 

 

R

g(x)dSx среднее значение

4πa2t2

функции g(x) на сфере Sat(x). Нетрудно понять, что решение, даваемое формулой (4.6), не приводит к диффузии волн: если в начальный момент времени в точке x возмущение отсутствовало, то оно исчезнет при достаточно больших временах .

Из формулы Кирхгофа можно вывести формулу Пуассона , решающую аналогичную

R

задачу в R2. С этой целью заметим, что интегрирование по сфере gdS в случае g =

Sρ

g(x, y) сводится к интегрированию по кругу:

Z

g(x, y) dS(x,y,z) = 2

2

Z2 2

g(x, y)cos(~n, ~z)

= 2ρ

2

Z2 2

g(x, y)

ρ2

x2

y2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

Sρ

x +y 6ρ

 

 

x +y 6ρ

 

p

 

− −

 

 

В результате такой редукции мы получаем формулу , аналогичную (4.6):

u(x, y, t) = θ(t)

 

v(x, y, t; τ)dτ +

∂t t Cat(x,y)[u0]

+

t Cat(x,y)

 

 

t

 

 

 

 

 

Z e

n Mf

o

n Mf

 

 

o

[u1] ,

0

в которой, однако, усреднение устроено по-другому, а именно

Mat(ξ,η)[g] = 2πat

Z

 

(at)2

(x

 

ξ)2

(y

 

η)2 ξ,η .

1

 

 

 

g(ξ, η)

 

 

 

 

f

 

Cat(x,y)

p

 

 

 

 

ãäå Cat(x, y) – круг радиуса at с центром в точке (x, y).

Решение, описываемое формулой Пуассона, демонстрирует эффект диффузии волн: возмущение, порожд¼нное любыми начальными данными, пройдя точку (x, y), не исчезает в этой точке за конечное время.

4.2.5Задачи Коши для уравнения теплопроводности

Рассуждая точно так же, как и в начале раздела 4.2.1 , мы приходим к следующей постановке обобщенной задачи Коши для уравнения теплопроводности :

Определение 4.3. Обобщенной задачей Коши для уравнения теплопроводности называется (понимаемое в смысле обобщенных функций) уравнение

− a2

u(~x, t) = f(~x, t) + u0(~x) · δ(t) , f(~x, t)|t<0 ≡ 0 ,

∂t

45