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

Lektsii_zubova_1

.pdf
Скачиваний:
6
Добавлен:
03.06.2015
Размер:
303.94 Кб
Скачать

Задача Коши для уравнения малых колебаний струны.

 

utt

a2uxx = 0

 

 

 

 

t=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(y, x)

 

= u0(x)

l < x < l

( )

t=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ut(t, x)

 

= u1(x)

 

 

 

 

 

 

(dx)2 a2(dt)2 = 0 (dx + adt)(dx adt) = 0

{

ξ = x + at η = x at

yˆξη = 0, в силу линейности замены от вторых производных добавок нет

uˆ(ξ, η) = f (ξ) + g(η)

u(x, y) = f (x + at) + g(x at) - сумма двух волн. Даламбер в 1747 году открыл этот факт. Задача Коши будет "хорошей".

u(t, x) =0= f (x) + g(x) = u0(x),

 

 

x

< l

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

x-at=0

 

 

< l

 

 

 

 

 

 

 

 

 

 

 

ut(t, x)tt=0= a f (x)

ag(x) = u1(x|),|

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

| |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g(x)]=

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

, ãäå

 

 

 

 

 

 

 

 

 

[ f (x)

 

u1(x)

 

 

 

 

f (x)

 

 

g(x) =

U1

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U1(x) = l

u1(ξ)dξ + C, òî åñòü:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

f (x) =

1

u0(x) +

1

 

U1(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

x

< l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x) =

 

 

 

u

(x)

 

 

 

 

 

 

U

(x)

 

| |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

u0(x + at) + u0(x at) U1(x + at) U1(x at)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

( ,

x

) =

 

(x+at)+

 

U (x+at)+

 

u (x at)

 

 

 

U (x at)

 

 

+

 

 

 

u t

 

2

u

0

 

 

 

 

 

 

 

2a 1

 

 

 

 

 

 

2 0 − −2a 1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x+at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

+

 

)

U1

(

 

 

at

) =

 

 

 

 

 

 

(ξ)dξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U1

 

x

 

at

 

 

x

 

 

 

 

(

xat

u1

 

 

 

x+at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u0(x + at)

+

u0

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x, t) =

 

 

 

x at

 

+ xat

u1(ξ)dξ

 

 

 

( )

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это формула Даламбера, которую Эйлер доказал в 1748 году.

t

Найд¼м область, в которой решение определено однозначно.

Ýòî ðîìá Q, ограниченный неравенствами |x + at| < l; |x at| < l

x

"можно заглянуть в прошлое"

Эта область будет максимальной, так как например при x > l можно взять любое непрерывное продолжение u0(x), и получить решение в боль-

шем ромбе.

мы только что доказали следующую теорему:

Теорема 2.1 Пусть u0(x) C2((−l, l)), u1(x) C1((−l, l)), тогда Задача Коши (*) имеет в Q единствен-

ное решение u(t, x) C2(Q), оно называется классическим, и представимо формулой Даламбера (**)

Покажем, что решение непрерывно зависит от u0 è u1

 

utt1

a2uxx1 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

utt2 a2uxx2 = 0

 

 

 

u1

 

t=0

= u01

 

 

|

|

< l

 

 

 

(1)

 

 

 

 

 

u2

 

t=0

= u02

| |

< l

(2)

 

 

(x), |x|

 

 

 

 

 

 

 

 

 

(x), |x|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u1

 

 

= u1

(x), x < l

 

 

 

 

 

 

 

 

 

 

u2

 

= u2

(x), x < l

 

 

|

t

 

t=0

 

1

 

 

|

 

 

 

|

 

 

1

 

|

 

 

 

 

 

t

t=0

1

 

 

 

 

u

0

 

 

0

(x)

6

0

;

1

(x)

(x)

6

δ

1

|

x

|

 

 

 

 

 

 

 

 

 

(x)

 

u

 

 

δ

 

u

u

 

 

 

 

< l

 

 

 

 

 

v0(x) = u01(x) − u02(x);

 

v1(x) = u11(x) − u12(x)

 

 

 

 

 

 

 

 

11

u1(t, x) − u2(t, x) = v(t, x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vtt a2vxx = 0

 

(x, t) Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t=0 = v0(x)

 

 

|

|

 

 

 

|

 

 

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

|x|

 

< l

 

|v0(x)|

6 δ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vt

 

 

 

= v1(x)

 

 

x < l

 

 

v1(x) 6 δ1

 

 

 

 

 

 

 

 

v(t, x) 6 |

 

 

 

 

 

|

 

|

 

 

 

 

 

|

 

+

 

 

 

x+at

v (ξ)dξ 6 δ

 

+ sup v (ξ)

 

 

 

 

|

6

 

 

(

 

+

 

+

 

(

 

 

 

1

 

 

 

 

 

+

+

 

 

 

 

 

 

v0

x

at

)

 

v0

x

at)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

at x

at

 

 

|

 

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xat

|

1

 

 

|

 

0

 

 

 

|

1 |

|

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2a

 

 

 

|ξ|<l

 

 

 

2a

 

 

 

 

6 δ0 +

1

δ12at = δ0

+ δ1t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Åñëè l < +∞, òî ýòî 6 δ0 + δ1

 

→ 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если ограничено время, то это 6 δ0 + δ1T → 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача (*) корректна.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

 

{

L(x, D)u(x) = f (x),

 

x Ω Rn

 

 

 

 

 

 

Линейная задача.

 

 

 

 

 

 

 

 

Bj(x, D)u(x) = gi(x),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x Γ (j = 0, n − 1)

 

 

 

 

 

 

 

 

 

 

L(x, D) - линейный оператор порядка p.

Γ Ω - поверхность

Пусть

1)Линейные нормированные пространства F(Ω) è H(Ω) â Ω

2)ËÍÏ G0(Γ), . . . , Gn−1(Γ) в Γ такие, что f (x) F(Ω) è gj(x) Gj(Γ) (j = 0, n − 1) Краевая задача (*)

n−1

имеет единственное решение u(t, x) H(Ω), и справедлива оценка u H(Ω) 6 C· f F(Ω) + Cj gj Gj(Γ)

j=0

C, Cj - универсальные константы.

В этом случае говорят, что краевая задача (*) является корректной в системе пространств

F(Ω), G0(Γ), . . . , Gn−1(Γ), H(Ω)

В задаче Коши было u0(x) C2(. . .); u1(x) C1(. . . ); решение было C2(. . .) по норме C.

В одной системе пространств задача может быть корректной, а в другой - нет.

Пример Адамара

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

 

uxx + uyy = 0

 

(x, y) Ω R2

 

 

 

y=0

 

 

 

 

 

 

u

 

(x) = e− √n cos(nx) → 0

 

y=0

= u0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

uy

 

 

= u1(x)

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Несложно проверить, что решениями будет

(равномерно, и все производные тоже)

u(x, y) = en cos(nx)· cosh(ny)

Рассмотрим точку (x , y ); x = 0, y > 0, тогда:

n)

→ +∞, ïðè n → +∞

u(x , y ) = e

n cosh(ny ) > 2 e− √neny =

2 eny

(1 − y

 

 

 

1

 

 

1

 

1

 

 

 

То есть ни в какой разумной норме эта задача корректной не будет .

Обобщ¼нное решение задачи Коши

До сих пор мы рассматривали задачи Коши с гладкими начальными данными. А интересно было бы обобщить наш результат на более общий случай, например как вед¼т себя гитарная стуна, при начальной деформации, отпущенная без начальной скорости.

Попробуем аппроксимировать функциями класса C2

1)Определение Пространство L2((−l, l)) - линейное нормированное пространство функций u(x), èí-

l

тегрируемых по Риману( в несобственном смысле) и таких, что u L2((−l,l)) =

u(x) 2dx < +∞

 

l

12

(u, v) =

l

u(x)v(x)dx - скалярное произведение.

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

Пусть l - конечное число.

 

 

 

 

 

 

 

2)Пусть u0k (x) C2([−l, l]);

u1k (x) C1([−l, l])

u0(x)

 

uttk a2uxxk = 0, (t, x) Q

 

 

 

 

 

 

 

 

 

 

 

 

uk

 

t=0

= u0k

 

 

uk(t, x) C2(Q)

 

 

 

 

(x), x [−l, l]

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

uk

 

 

 

= uk

(x), x [ l, l]

 

 

 

 

 

 

 

t=0

 

 

 

 

 

 

 

t

1

 

 

 

 

 

 

 

 

 

Доказательство непрерывности в замыкании проходит анало-

 

 

гично доказательству, которое мы привели ранее.

 

 

3) а) Пусть u0(x) C([−l, l]) и последовательность u0k (x)

 

 

C2([−l, l]) такая, что

 

 

 

 

 

 

 

 

u0k (x) u0(x) C([−l,l]) = xmax[ l,l] |u0k (x) u0(x)| → 0 ïðè k → ∞

 

 

 

 

 

 

 

 

 

последовательность u1k (x) C1([−l, l]) такая, что

 

 

б) Пусть u1(x) L2([−l, l]);

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

u1k (x) − u1(x) L2([−l,l]) = l |u1k (x) − u1(x)|dx → 0 ïðè k → ∞

 

 

Определение: Говорят, что u(t, x) C(Q) явлвяется обобщенным (сильным) решением задачи Коши:

 

 

 

 

 

 

 

 

u t=0 = u0

(x) C([−l,

l])

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

 

utt

t=0

a2uxx

=

0; (t, x)

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

åñëè

 

k

 

 

2

 

 

k

 

 

 

 

 

1

 

 

 

 

 

 

таких, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ut

 

 

= u1(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2([

 

l, l])

 

 

 

 

 

 

 

u0

(x)

 

C

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

([ l, l]);

ua x

)

 

 

 

 

([

l l

])

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

u0k (x) − u0(x) C([−l,l]) → 0 ïðè k → ∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u1k (x) − u1(x) L2([−l,l]) → 0 ïðè k → ∞

 

 

 

 

 

 

 

 

 

 

k( ,

 

 

 

 

выполняется условие

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max

( ,

)

)

 

0 ïðè

 

(t, x)

u(t, x)|C(Q)

=

| →

k → ∞

 

 

 

 

 

 

 

|u

(t,x)

Q

|u t x

 

u t

x

 

Докажем следующую теорему:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 2.2: u0(x) C([−l, l]);

 

 

u1(x) L2([−l, l]) единственное обобщ¼нное решение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

задачи Коши (*), и оно представимо формулой Даламбера.

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

x+at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)u1(x) L2([l, l]) xat

u1(y)dy (t, x) Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|u1(y)|· 1· dy 6

 

 

 

 

 

·

 

 

 

 

 

 

 

 

 

x+at

u1(y)dy

6

x+at |u1(y)|dy 6

 

 

 

l

 

l u1(y)2dy

 

l dy 6 u1 L2([−l,l])·

 

 

 

 

 

 

 

 

 

 

 

2l

 

 

 

 

 

 

 

 

 

xat

 

 

 

l

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

xat

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x+at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u0(x

+ at) + u0(x

 

at)

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)u(t, x) =

 

 

 

 

 

 

 

+

 

 

 

 

xat

u1(y)dy C(Q)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k( , )

( , ) =

u0k (x + at)

+ u0k (x at) + 1

 

x+at

k ( )

 

 

 

u0(x + at) − u0(x at)

+ 1 x+at

 

( )

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|u t x u t x |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xat

u1 y dy

 

 

 

 

 

 

 

 

 

xat

u1

y dy

 

 

 

 

 

2

 

 

 

 

 

 

 

2a

 

 

 

 

 

2

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u0k (x + at) − u0(x + at)

 

 

 

u0k (x

at) − u0(x at)

 

 

1

 

x+at

k

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xat

|u1(y) − u1

(y)|dy 6

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

u0k (y) − u0(y) C([−l,l])(

 

 

+

 

) +

 

u1k (y) − u1(y) L2([−l,l]) 2l → 0 ïðè k → ∞

 

 

 

 

 

 

 

 

 

 

2

 

2

2a

 

 

 

 

 

Получили равномерную сходимость. То есть формула Даламбера действительно да¼т обобщ¼нное решение задачи Коши. Более того оно будет единственным, так как все последовательности сходятся кu(t, x)

13

Смешанная задача Коши для полубесконечной струны с 1 закрепл¼нным концом.

u t=0

= u0

(x)

ut

t=0

= u1(x),

x > 0 (начальные условия)

( )

 

utt a2uxx = 0;

t

> 0, x > 0

 

 

u = 0, t > 0

граничные(краевые условия)

 

 

x=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u0(x) C2([0, +∞)) u1(x) C1([0, +∞))

классическое решение u(t, x) C1(t > 0, x > 0) ∩ C2(t > 0, x > 0)

Общее решение нашего уравнения:

u(t, x) = f (x + at) + g(x at) Подставим данные Коши:

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{

u

 

 

0 = f (x) + g(x) = u0(x), x > 0

 

 

 

 

utt=t=0 = a f (x)

 

g(x) = ua(x), x > 0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

U1(x) =

u1(ξ)dξ + C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

x > 0

 

 

 

 

 

 

 

 

 

 

 

(x)

 

 

 

(x)

 

 

 

 

 

 

 

 

1

 

+

1

 

 

 

 

 

 

f (x) = 2 u0

2aU1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

x=0

 

f

 

 

 

 

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

u0

(x)

(x)

 

 

 

 

 

g(x) =

 

 

U1

 

= (at)

+ g( at) = 0

t > 0

 

 

 

 

 

 

 

t

x-at=0

 

 

x

 

x+at=0

В штрихованной области решение однозначно определено начальными условиями. Используя граничные условия, найд¼м решение в оставшейся области.

Пусть η = −at, η 6 0

 

g(η) = − f (−η)

 

 

 

f (x) =

1

u0

(x) +

1

U1

(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2a

 

2

 

10

 

 

 

2a

11

 

следует из g(η) =

f ( η)

 

 

 

g(x)

=

g x

 

 

 

 

 

 

 

 

 

 

 

 

 

(

) =

1

u

 

(x)

 

 

1

U

(x),

x > 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2a

 

 

 

 

 

 

 

 

 

g(x) =

 

u

(

x)

 

U (

− −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x), x < 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x) C2([0, +∞)) g(x) C2([0, +∞)) g(x) C2((−∞, 0])

Сошь¼м решения:

 

 

 

 

 

 

g(+0) = g(−0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g(+0) = g(−0)

 

 

 

g(x) будет

 

C2((−∞, ∞))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g′′(+0) = g′′( 0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)g(+0) = g(−0)

1

u0(+0) −

1

U1

(+0) =

1

u0

(+0) −

1

 

U1

(+0) u0(0) = 0

2

2a

 

2

 

2a

2)g(+0) = g( 0)

 

1

u

(+0)

 

1

u

 

(+0) =

 

1

u(+0) +

1

u

(+0)

 

u

(0) = 0

 

 

 

 

 

 

 

2 0

 

 

 

2a

1

 

 

2 0

 

 

2a

1

 

 

 

1

 

3)g′′(+0) = g′′( 0)

 

 

 

1

u′′(+0)

 

1

u(+0) =

 

 

 

1

u′′(+0)

 

 

 

1

u(+0)

 

u′′(0) = 0

 

 

2

2a

 

 

 

2a

 

 

0

 

 

 

1

 

 

 

2

0

 

 

1

 

0

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

 

u0

(

x

+

at

2 u0

(

x at

)

+ 2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(t, x) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u0

(x + at) + u0

(x at)

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x+at

u1(ξ)dξ

xat

x+at

u1(ξ)dξ

atx

{

x + at > 0

x at > 0

{

( )

x + at > 0

x at 6 0

14

Вообще говоря, надо добавить t > 0, но оказыватеся решение представимо в таком же виде и при t < 0 - "можно заглянуть в прошлое"

Докажем следующую теорему:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 2.3 Пусть в смешанной задаче (*) функции u0(x) è u1(x) удовлетворяют:

 

 

 

 

 

 

а) условию гладкости: u0(x) C2([0, +∞));

u1(x) C1([0, +∞))

 

 

 

 

 

 

 

 

 

 

 

б) условию согласованности: u0(0) = u1(0) = u0′′(0) = 0

 

 

 

 

 

 

 

 

 

 

 

 

Тогда смешанная задача (*) имеет единственное решение(классическое) u(t, x) C2(t > 0, x > 0) è îíî

 

 

 

представимо формулой (**)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Метод продолжений

 

 

 

 

 

 

 

 

 

 

 

Попробуем свести задачу (*) к задаче Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где в качестве uˆ0 è uˆ1 выберем:

utt a2uxx = 0, t

> 0, x R1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{

 

u

 

= uˆ0(x); ut

 

t=0

= uˆ1(x), x R1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t=0

 

 

 

 

 

 

 

 

 

 

 

 

uˆ0(x) = {

u (x), x > 0

 

 

 

 

 

 

 

 

 

 

 

uˆ1(x) = {

u (x), x > 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C2(R1)

 

 

 

 

 

 

 

 

 

0u0(−x), x < 0

C2(R1)

 

 

 

 

1u1(−x), x < 0

 

 

 

 

 

 

По формуле Даламбера имеем:

 

 

x+at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

(

x

+

) ˆ

(

 

)

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

uˆ(t, x) =

u0

 

 

at u0

 

x at

 

+

 

 

 

 

xat

uˆ1(ξ)dξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x at < 0, x + at > 0

 

 

 

 

 

1

 

x+at

 

 

 

 

 

 

1

 

0

 

 

u0(x + at)

 

 

1

x+at

 

 

 

 

ˆ( ,

x

) =

u0(x + at) − u0(at x)

+

 

0

 

u1

(ξ)

d

ξ

 

 

 

( ξ)

d

ξ =

u0(at x)

+

atx

u1

(ξ)

d

ξ

 

2a

 

2a

 

2

2a

u t

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

atx u1

 

 

 

 

 

Таким образом, мы свели смешанную задачу к задаче Коши, и тем самым доказали теорему. Пример (отражение волны от закрепл¼нного конца) Воспользуемся методом продолжений. Для это-

го построим продолжение начального возмущения, симметричное относительно (t,x) = (0,0). Каждая из начальных волн состот из двух с вдвое меньшей амплитудой, которые побегут в разные стороны.

t

x

t

x

t

x

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача Коши для волнового уравнения в R2,3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегралы, зависящие от параметров

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема: Пусть f (x, y), x

 

 

x Rn, y

 

y Rn непрерывна в

 

x ×

 

y. Тогда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)J(y) = x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x, y)dx C(

y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) Åñëè

 

 

 

(x, y)

C(Ωx ×

y), òî J(y) имеет непрерывные производные

 

,

y y,

 

 

 

 

 

 

 

yk

yk

 

 

 

 

 

 

 

 

 

J(y)

= x

f (x, y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ïðè÷¼ì:

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yk

 

 

yk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим функцию, зависящую от параметра:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ug(t, x, τ) =

1

 

|ξ−"x|=at

g(ξ, τ)dSξ, ãäå:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

aπa2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a > 0, t > 0 x = (x1, x2, x3) R3, τ > 0, ξ = (ξ1, ξ2, ξ3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Лемма 3.1

Пусть:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)g(ξ, τ) C({ξ R3, τ > 0})

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)Dξα g(ξ, τ) C({ξ R3, τ > 0}) α, |α| 6 p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тогда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)Dtα,xug(t, x, τ) C({t > 0 x R3, τ > 0}) α, |α| 6 p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) lim ug(t, x, τ) = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t→+0

 

 

 

 

 

 

 

 

 

 

 

ug(t, x, τ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3)Åñëè p > 1, òî lim

 

= g(x, τ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t→+0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ãäå C(

) :

 

 

C(Ω), допускающие непрерывное продолжение на границу.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1) Зафиксируем (t, x)

пусть η =

ξ − x

ξ = x + atη

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|ξ − x| = at; |η| = 1;

 

 

dSξ = a2t2dSη

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ug(t, x, τ) =

 

 

a2t2

 

|"η|=1

 

g(x + atη, τ)dSη = t· Jg(t, x, τ), ãäå:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Jg(t, x, τ) =

 

1

|"η|=1

g(x + atη, τ)dSη

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g(ξ, τ) C({ξ R3, τ > 0})

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g˜(t, x, η, τ) = g(x + atη, τ) C({t > 0,

x R3,

τ > 0, |η| = 1}), как суперпозиция непрерывных.

 

 

 

Jg(t, x, τ) C({t > 0, x R3,

 

 

τ > 0}), по теореме, сформулированной в начале параграфа.

 

 

 

 

 

 

Dtα,x Jg(t, x, τ) C({t > 0,

x R3,

τ > 0})

α, |α| 6 p по той же теореме.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)lim

( ,

x

, τ) = lim[

t·

( ,

x

, τ)] = lim

lim

( ,

x

, τ) = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t→0 ug t

 

 

 

 

 

 

t→0

Jg

 

t

 

 

 

 

 

t→0 t·

t→0

Jg t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

так как последний предел при фиксированных x, τ. Можно продолжить нул¼м в 0.

 

 

 

 

 

 

3)lim

ug(t, x, τ)

 

= lim

 

 

 

 

 

( ,

 

, τ)

= lim

( ,

 

, τ) + lim

Jg

=

 

 

(0,

 

 

, τ) + 0

 

 

 

=

1

 

 

 

(

 

, τ)

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|"η|=1

 

 

 

t→0

 

 

 

t

 

 

t→0

t

[t· Jg t

x

 

 

]

 

t→0 Jg t x

 

t→0 t· t

 

Jg

 

x

 

 

 

 

· C

 

g

x

 

dSτ

 

1

|"η|=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g(x, τ)

 

dSη = g(x, τ)

 

 

 

 

 

 

 

 

 

 

 

Доказательство завершено .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

utt a2(ux1x1

+ ux2x2 + ux2x2 ) = 0,

t > 0, x R3

 

простая задача Коши

Теорема 3.1

Пусть в задаче

Êîøè

 

 

3

( ) −

 

 

{ u

t=0 = 0; ut

t=0 = u1(x), x R3

 

 

 

 

 

 

 

 

 

u1(x)

 

C(R )

 

 

 

 

 

 

|x"ξ|=at

 

 

 

 

 

Тогда u(t, x) =

 

1

u1(ξ)dSξ, ïðè÷¼ì:

 

 

 

 

a2t

 

 

 

1)u(t, x) C2({t > 0, x R3})

2)u(t, x) - классическое решение задачи Коши (*)

Это формула Пуассона - Кирхгофа, открытая в 1818 году.

Она является хорошим кандидатом на классическое решение, так как граничные условия выпол-

нены в силу леммы 3.1

"

u(t, x) =

t

u1(x + atη)dSη C2({t > 0, x R3})

 

 

|η|=1

Дифференцируем по t, x :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xu(t, x) =

 

 

t

|"η|=1 ξu1(x +ξatη)dSη

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

| {z }

 

t

 

 

 

3

u1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|"η|=1

 

 

 

 

 

 

 

|"η|=1

 

 

 

 

 

 

 

 

 

ut(t, x) =

 

 

 

 

 

 

 

 

 

 

 

∂ξk (x + atη)· ηk· a· dSη

 

 

 

u1(x + atη)dSη +

k=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

k

=

ξk xk

=

ξk xk

= η

, ãäå

n - вектор нормали к сфере

 

 

 

 

 

|

ξ

x

|

 

 

 

 

at

k

 

 

→−

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(t, x)

1

 

 

3

u1(ξ)

 

 

u(t, x)

1

 

 

u1(ξ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|ξ−"x|=at

 

 

 

 

 

 

 

 

 

 

|ξ−"x|=at

 

 

 

ut(t, x) =

 

 

 

t +

at

∂ξ nk(ξ)Sξ = t +

 

→−

Sξ =

 

 

 

 

 

k=1

at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

n

at

x

1 + I

at at

Используем формулу Остроградского-Гаусса:

$

$

"

 

ξu1(ξ)dξ =

div( u1(ξ))dξ =

|x−ξ|<at

|x−ξ|<at

|x−ξ|=at

$ at "

I = ∆ξu1(ξ)dξ = dρ ∆ξu1(ξ)dSξ

" u1(ξ)

(n(ξ), u1(ξ))dSξ = ∂→− dSξ n

 

|ξ−x|=at

It = a

" ξu1(ξ)dSξ

|x−ξ|<at

 

 

u(t, x)

 

0

 

|x−ξ|=ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|x−ξ|=at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

u

 

u

 

 

I

 

 

 

I

 

 

1 u

I

 

u

 

 

 

I

 

I

 

I

utt(t, x) =

 

[

 

 

+

 

 

] =

t

 

 

 

 

 

+

 

t

=

 

 

(

 

 

+

 

)

 

 

 

 

 

 

+

t

=

t

t

t

 

at

t

t2

 

at2

at

t

t

at

t2

 

 

at2

 

at

at

Отсюда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

a2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

utt a2xu =

 

 

|x"ξ|=at

ξu1(ξ)dSξ

 

|"η|=1

ξu1(x + atη)dSη =

 

 

 

 

 

 

 

 

at

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

a2t

ξu1(x + atη)dSξ(

 

1

 

) = 0

 

 

 

 

 

 

 

=

 

 

|x"ξ|=at

ξu1(ξ)dSξ

 

 

|x"ξ|=at

 

 

 

 

 

 

 

 

 

at

a2t2

 

 

 

 

Получаем дополнительную информация, что utt = 0 ïðè t = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

utt a2(ux1x1 + ux2x2 + ux2x2 ) = 0, t > 0, x R3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

 

{ u t=0 = u0(x); ut

t=0 = 0, x R3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u0(x) C3(R3)

17

{

vtt a2(vx1x1

+ vx2x2 + vx2x2 ) = 0,

t > 0,

 

x R3

решение класса

 

3

v t=0 = 0; vt

t=0 = u0(x), x R3

 

 

 

 

 

( ) −

 

C

 

1

 

 

 

 

3

 

 

 

 

3

 

 

 

 

 

 

 

|x"ξ|=at

u0(ξ)dSξ

 

C

{

 

 

R

 

}

)

 

 

 

 

a2t

 

 

 

 

 

v(t, x) =

 

 

( t > 0, x

 

 

 

 

 

 

Докажем, что решением ( ) будет u(t, x) = vt(t, x) C2({t > 0, x R3}) vtt a2xv = 0

0 = vttt a2(∆xv)t = vttt a2x(vt) = utt a2xu = 0

u t=0 = vt t=0 = u0(x)

ut t=0 = vtt t=0 = 0

 

То есть, мы только что доказали следующую теорему:

 

Теорема

3 3

 

 

 

 

 

 

u(t, x) = t[

a2t

|x"ξ|=at u0(ξ)dSξ], t > 0, x R

 

 

 

3.2 Функция u(t,x):

 

 

 

 

1

 

3,

ãäå u0(x) C (R ) является классическим решением задаче Коши ( )

 

 

 

 

 

utt a2

(ux1x1 + ux2x2 + ux2x2 ) = f (t, x), t > 0, x R3

 

 

Для решения этой задачи

Коши, используем метод Дюамеля

 

 

 

 

{

u

 

 

 

0; ut

 

 

 

 

 

( )

 

 

 

t=0 =

t=0 = 0, x R3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим следующее семейство задач:

 

 

 

 

 

 

{

W

 

(t, x, τ)

 

a2xW(t, x, τ) = 0, t > τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x R3; f − класса C2 ïî x

 

 

 

Wttt= 0;

Wt t= f (τ, x),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решения для каждого

W(x, t, τ) =

1

a2(t − τ)

τ следует из формулы Пуассона-Кирхгофа:

"

f (τ, ξ)dSξ

 

 

 

 

 

 

 

 

|ξ−x|=a(t−τ)

 

 

 

 

 

W(t, x, τ) C2(t, τ, x) - непрерывно дифференцируема

 

 

 

Wf (x, t, τ) =

1

|ξ−"x|=at

f (τ, ξ)dSξ

 

 

 

 

 

a2t

 

 

 

 

 

Пусть Dxα f (t, x) C({t > 0, x R3}) α, |α| 6 2

 

Только при таком предположении, мы сможем корректно решить задачу.

 

 

 

Ïî

лемме 3.1

имеем:

 

 

 

 

 

 

Dtα,xWf (t, x, τ) C({t > 0, x R3, τ > 0}) α, |α| 6 2

 

 

 

W(t, x, τ) = Wf (t − τ, x, τ) Dtα,xW(t, x, τ) C({t > τ, x R3, τ > 0}) α, |α| 6 2

 

 

Определим:

 

 

 

 

 

 

 

 

 

 

u(t, x) = 0 t W(t, x, τ)dτ Утверждается, что е¼ можно выбрать в качестве решения (***)

 

 

1)Dtα,xu(t, x) C({t > 0, x R3}) α, |α| 6 2

 

 

 

 

 

2)∆xu(t, x) = 0 t xW(t, x, τ)dτ

 

 

 

 

 

3)ut(t, x) = W(t, x, t) + 0 t Wt(t, x, τ)dτ, íî W t= 0: первый член = 0

 

 

 

 

 

 

 

 

 

t

 

 

t

 

 

 

 

 

 

 

o

 

 

 

2

2

 

 

 

 

 

 

Wtt(t, x, τ)dτ = f (t, x)

0

a xW(t, x, τ)dτ = f (t, x) + a xu(t, x)

utt(t, x) = Wt(t, x, t) +

+

 

 

То есть это действительно решение. Проверим начальные данные:

 

u

 

0

 

 

 

 

 

 

 

 

t=0 =

0 0

W(t, x, τ)dτ = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ut

0

 

 

 

 

 

 

 

 

 

 

Wt(t, x, τ)dτ = 0

 

 

 

 

 

 

 

 

 

t=0 =

 

 

 

 

То есть это действительно классическое решение задачи Коши (***).

Его можно интерпретировать как набор точечных источников, которые потом суммируются.

18

Сформулируем только что доказанный факт в теорему:

Теорема 3.3 Пусть Dαx f (t, x) C({t > 0, x R3}) α, |α| 6 2

Тогда u(t, x) = 0

t

(

a2(t − τ)

"

f (τ, ξ)dSξ)dτ

 

 

 

1

 

 

|ξ−x|=a(t−τ)

1) Dαt,xu(t, x) C({t > 0, x R3}) α, |α| 6 2

2) Является классическим решение задачи Коши (***)

 

 

 

{

utt a2xu = f (t, x), t > 0, x R3

 

 

 

 

 

 

 

= u1(x), x R3;

 

 

 

u t=0

= u0(x); ut t=0

 

 

 

 

 

 

 

 

 

Теорема 3.4 Пусть в задаче Коши (****):

 

3

3

 

 

 

 

 

 

 

à)u0(x) C2

(R3)

 

 

 

 

 

 

 

á)u1(αx) C (R

 

3

})

α, |α| 6 2

 

â) Dt,xu(t, x) C({t > 0, x

R

 

Тогда функция:

"

u0(ξ)dSξ] + 1a2t

" u1(ξ)dSξ +

u(t, x) = t[1a2t

 

 

 

 

 

 

 

 

 

( )

a2

(

ξ x

 

 

)dξ C (R )

1

$

f

t

|ξ−x|

, ξ

 

 

 

a

 

2 3

 

 

|

 

|

 

 

 

|ξ−x|=at

|ξ−x|=at

|ξ−x|<at

-является классическим решением задачи Коши (****)

Докажем. Используем принцип суперпозиции.

{

utt0

a2xu0

= 0

 

 

{

utt1

a2x10 = 0

 

{

utt2

a2xu2 = f (t, x)

u0

 

 

 

= 0

u1

 

 

= 0; ut1

 

= u1(x)

u2

 

 

= 0; ut2

 

= 0

t=0 = u0(x); ut0

t=0

t=0

t=0

t=0

t=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(t, x) = u0(t, x) + u1(t, x) + u2(t, x)

Насч¼т единственности пока ничего сказать не можем. Осталось привести u2(x, t) к нужному виду:

 

 

 

 

 

 

 

ρ

 

 

 

 

dρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(пусть a(t − τ) = ρ;

τ = t

 

; dτ

= −

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

1

 

 

|ξ−x"|=a(t−τ)

 

 

 

 

 

 

0 dρ

"x|=ρ

 

(

 

 

ρ

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u2(t, x) =

 

(a2(t

 

 

 

 

 

 

at a2ρ

 

 

 

 

 

0

 

 

τ)

 

 

f (τ, ξ)dSξ)dτ = −

 

 

 

f t a , ξ dSξ =

=

 

 

 

 

 

[

 

(

 

 

 

 

 

)

dSξ dρ =

 

 

 

(

 

 

 

 

)

dξ

 

 

 

 

 

 

 

 

 

 

ξ x

 

 

 

 

 

ξ x

 

 

 

a2

0

 

 

 

a2

 

 

 

 

 

 

 

 

1

 

at

 

 

 

f t

|ξ−a x|

, ξ

 

 

1

 

 

f t |ξ−a x| , ξ

 

 

 

 

 

 

 

 

 

 

 

 

"x|=ρ

| − |

 

 

 

 

 

$x|<at

 

| − |

 

 

 

 

 

Это так называемый запаздывающий потенциал. Принцип Гюйгенса

f(t,x)=0

t = R0

0 a

x

19

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

Пусть где-то произош¼л взрыв (то есть f = 0). Проследим за распространением волны. Через некоторое маленькое время t будет тишина, через некоторое время t0

фронта, потом будет сбор информации по сфере и ещ¼ через некоторое времязаднего. То есть передний и задние фронты ч¼тко выражены.

Возмущение, локализованное в пространстве приводит к действию, локализованному по времени.

Решение задачи Коши для волнового уравнения в R2

 

 

 

 

 

 

 

 

 

 

 

utt a2(ux1x1 + ux2x2 ) = f (t, x1, x2), t > 0, x = (x1, x2) R2

 

 

 

 

Используем метод спуска

èç

R

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{ u

 

 

 

 

 

 

 

 

 

t=0

= u1

(x), x R2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t=0

= u0

(x); ut

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

utt a2(ux1x1 + ux2x2

 

+ ux3x3 ) = f (t, x1, x2), t > 0, x = (x1, x2, x3) R3

 

 

 

 

Используем результат

 

теоремы 3.4,

 

тогда в н¼м

ξ = (ξ1

, ξ2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{

u

 

 

 

 

 

 

(x1, x2); ut

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t=0 = u0

 

t=0 = u1(x1, x2),

 

 

x R2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

|ξ−"x|=at

u11, ξ2)dSξ = V(t, x).

А будет ли оно подходить в качестве решения? Смотрим на

 

 

 

 

a2t

Заметим, что сфера разбивается на две полусферы Satè Sat+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sat± ξ3 = x3 ±

 

 

 

 

 

 

 

 

 

 

= x3 ±

 

 

 

 

 

x| < at ξ= (ξ1, ξ2)

a2t2 − (ξ1 x1)2 − (ξ2 x2)2

a2t2 − |ξx|2

 

"Sat

 

"Sat

V(t, x) =

 

1

 

 

"Sat

u11, ξ2)dSξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

=

 

a2t

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вспомним факт из математического анализа:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S : ξ3 = f 1, ξ2) - поверхность, тогда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

"S

u(ξ)dSξ = "D

u[ξ1, ξ2, F1, ξ2)] √

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + Fξ21

+ Fξ22 dξ1dξ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

То есть, в данном

случае F1, ξ2) = x3 +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2t2 − |ξ

x|2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1 +

1 x1)2 + (ξ2

x2)2

 

 

 

 

 

 

 

 

 

at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + 2 + 2

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( , ) =

1

Fξ2

 

(ξ , ξ )

 

 

 

 

 

ξ

 

 

 

1

2

 

 

 

 

 

ξ

ξ =

1

 

 

2

)

 

1

 

 

ξ

 

Fξ1

 

 

 

 

 

a2t2

− |

 

x

 

 

 

 

 

 

 

2

t

2

− |ξ

 

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

 

 

 

 

 

 

 

 

a

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

V t a

 

 

"x|<at

u1 1 2

 

 

 

 

 

 

 

d 1d 2

 

 

 

 

"x|<at u1

 

 

 

d

 

 

 

 

a

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

a2t2

 

− |ξx|2

 

 

 

 

a2t2 − |ξx|2

 

 

 

Сформулируем теорему:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 3.5: Пусть в задаче Коши

{

utt a2(ux1x1 + ux2x2 ) = f (t, x), t > 0, x R2 u t=0 = u0(x); ut t=0 = u1(x), x R2

 

 

 

3

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

à) u0(x) C2

(R2)

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

á)

u1α x

)

C

(R )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

â) Dx f (t, x) C({t > 0, x R }) α, |α| 6 p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

|ξ−"x|<at

 

 

 

(ξ)

 

 

 

 

1

|ξ−"x|<at

 

 

(ξ)

 

 

 

t

 

1

|ξ−x"|<a(t−τ)

 

f (τ, ξ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− |

ξ

|

2 dξ]+

− |

ξ

|

 

 

[

a2t2

− |

ξ

|

2 dξ]dτ

u(t, x) = t[a

 

a2t2u0

 

 

x

a

 

a2t2u1

 

 

x

2 dξ+ 0

 

a

 

 

 

x

1)u(t, x) C2({t > 0, x R2})

2)является классическим решением задачи Коши.

 R2 принципа Гюйгенса не будет, при t > T0 решение всегда будет не нулевым. Это можно интер-

претировать двумя способами:

1) в координатах (x1, x2) интегрирование вед¼тся по кругу, который при достаточно больших значениях времени всегда будет пересекать область возмущения.

20

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