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

1

§7 Метод Фурье. Понятие корректности задач математической физики

 

Обоснование метода Фурье в простейших случаях.

 

 

>

utt = uxx; (t; x) G = f(t; x); t > 0; 0 < x < 1g;

 

 

 

 

= '0(x); u2t t=0 = '1(x); x

 

 

 

 

 

 

 

 

<

 

t=0

 

- =

 

x; 0 < x < 1 ;

 

 

 

8 u

2

f

(1)

>

j

 

j

 

 

 

 

g

>

 

 

 

 

 

 

 

 

 

 

 

:

 

 

= ujx=1 = 0; t > 0:

 

 

 

 

 

 

 

 

 

 

> ujx=0

 

 

 

 

 

 

 

 

Физический смысл задачи (1).

 

 

 

 

 

 

 

 

Класс решения задачи (1): u(t; x) 2 C2(G), u непрерывна в

G

, ut

непрерывна

в

G

.

 

 

 

 

 

 

 

 

 

 

 

Мы не будем, пользуясь идеями x6, заменять задачу (1) на задачу для системы

Ut = AU:

Применим метод разделения переменных к исходной задаче (1). С этой целью, будем искать у волнового уравнения

utt = uxx

частные, нетривиальные решения вида

u(t; x) = T (t)v(x);

удовлетворяющие краевым условиям

ujx=0 = ujx=1 = 0;

TT00 (t) = vv00 (x) = const = ¡a2 !

T 00 + a2T = 0; v00 + a2v = 0;

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

v(0) = v(1) = 0:

Итак, для нахождения функции v(x) мы вновь имеем задачу на собственные значения (спектральную задачу) того же вида, что и в x6:

(

v00 + a2v = 0; x 2 -; v(0) = v(1) = 0:

Нетривиальные решения этой задачи:

ak = k¼; vk(x) = sin k¼x; k = 1; 2; ::: :

Из уравнения

T 00 + a2T = 0

находим:

Tk(t) = Ak cos k¼t + Bk sin k¼t;

Ak, Bk - пока произвольные постоянные. Итак, функции

uk(t; x) = sin k¼xfAk cos k¼t + Bk sin k¼tg; k = 1; 2; :::

удовлетворяют уравнению и граничным условиям при любых Ak, Bk.

 

 

2

Попытаемся определить эти постоянные так, чтобы бесконечный ряд

 

 

1

 

 

Xk

(2)

u(t; x) =

uk(t; x)

 

=1

 

давал решение задачи (1). При t = 0:

ujt=0 = '0(x) = P1 Ak sin k¼x;

k=1

utjt=0 = '1(x) = P1 k¼Bk sin k¼x:

k=1

Пусть '0(x); '1(x) 2 C04(-). Тогда

Ak = 2 R1 '0(x) sin k¼xdx;

0

R1

Bk = 2 0 '1(x) sin k¼xdx;

причем

jAkj ·

const

; jBkj ·

const

:

 

 

 

k4

k5

Задача.

Покажите, что ряд (2) можно дифференцировать почленно два раза по t и x. Покажите, что u(t; x) ! '0(x), ut(t; x) ! '1(x) при t ! +0 равномерно по x.

Единственность:

uI ; uII ; u = uI ¡ uII ;

utt = uxx ! 2ututt = 2utuxx !

 

! (ut2)t = 2(utux)x ¡ (ux2)t !

 

! (ut2 + ux2)t = 2(utux)x ! 1R0

dx :

 

= 0 ! J(t) = J(0); J(t) = R0

1

 

dt

(ut2 + ux2)dx:

dJ(t)

 

 

 

Поскольку J(0) = 0, то ut = ux ´ 0 ! u ´ 0 в G. Непрерывная зависимость решения от начальных данных.

Разница между задачей (1) из x6 и x7: гладкость решения (2) можно повысить, только повышая гладкость начальных данных.

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

4x;yu = uxx + uyy; (x; y) 2 - ½ R2:

Определение. Непрерывные решения уравнения 4x;yu = 0 называются гармоническими функциями.

Задача Дирихле для уравнения Лапласа:

(

 

4x;yu = 0; (x; y) 2 -;

(3)

 

uj@- = f(x; y)j@-:

 

Физический смысл: задача о прогибе мембраны.

 

Пусть - = f(x; y); x2 + y2 < R2g - круг.

 

4x;yu = 0 !

( r

@

(r

@

) + @'@2u2 = 0;

 

 

@r

@r

(30)

 

ujr=R = f('):

 

3

y

R

r

0

x

Частные решения уравнения Лапласа

u(r; ') = A(r)B(');

r(rA0)0 = ¡B00 = a2 = const;

A B

B(')½= c1 sin a' + c2 cos a';

A(r) = c3ra + c4r¡a; (a =6 0); c5 + c6 ln r; (a = 0);

c1;2;3;4;5;6 - некоторые постоянные; B(') - периодическая с периодом 2¼, a = n, n = 0; 1; ::: .

Функции r¡n, n = 1; 2; 3; ::: имеют особенности при r = 0, поэтому непрерывные внутри круга частные решения уравнения Лапласа имеют вид:

un(r; ') = rn(bn sin n' + an cos n'); n = 1; 2; ::: :

 

При n = 0 (a = 0):

 

 

 

 

 

 

 

 

 

a0

 

 

 

 

 

 

 

 

 

 

u0(r; ') = const =

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

an, n = 0; 1; :::, bn, n = 1; 2; ::: - некоторые постоянные.

 

 

 

 

Линейная комбинация частных решений тоже решение

 

 

 

 

 

a0

N

 

 

rn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

u(r; ') =

2

+

 

Rn

(an cos n' + bn sin n'):

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

Пусть a0, an, bn, n = 1; 2; ::: - ограниченные постоянные

 

 

 

 

 

 

ja0j; janj; jbnj < M; M < 1:

 

 

 

 

Тогда ряд

 

a0

1

 

 

rn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

(4)

u(r; ') =

 

 

+

 

 

 

(an cos n' + bn sin n')

 

 

 

 

2

n=1

 

 

Rn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

является при r < R решением уравнения Лапласа.

 

 

 

 

В самом деле, ряд (4) можно представить так:

 

 

¾:

 

u(r; ') =

20

+ n=1

Re½

)Rn n ¡

n

)

 

 

 

 

 

a

1

 

 

 

 

(x + iy n(a

ib

 

 

Но ряд

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

1

 

an ¡ ibn zn = w(z); z = x + iy

 

a0 +

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

n=1

 

 

Rn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

является степенным рядом с радиусом сходимости не меньшим, чем R. Отсюда следует, что при r < R: w(z) - аналитическая функция, а u(r; ') = Re w(z)- гармоническая функция (ТФКП).

Если предположить равномерную сходимость ряда вплоть до границы круга r = R, то для граничных значений u(r; ') = f(') мы будем иметь представление рядом Фурье

f(') = a20 + X1 (an cos n' + bn sin n'):

n=1

Коэффициенты Фурье вычисляются по формулам

 

 

2R¼

 

>

1

2¼

 

>

 

 

 

8 an = ¼1

0

f(') cos n'd'; n = 0; 1; ::: ;

>

 

R0

 

>

 

 

:

 

 

< bn = ¼

f(') sin n'd'; n = 1; 2; ::: :

Задача.

Пусть f(') - периодическая функция с периодом 2¼, причем f Тогда

janj ·

const

; jbnj ·

const

:

 

 

n2

n2

(5)

(6)

2 C2([0; 2¼]).

Заметим, что эти неравенства обеспечивают сходимость рядов (4), (5).

Мы показали, что ряд (4) с коэффициентами (6) является в круге r · R решением задачи Дирихле для уравнения Лапласа (3):

(

4x;yu = 0;

ujR = f(');

если f(') достаточно гладкая.

Представим решение задачи (3) в другой форме. Преобразуем ряд (4), воспользовавшись формулами (6):

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

 

1 rn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(r; ') =

 

 

 

+

 

 

 

 

 

(an cos n' + bn sin n') =

 

 

 

 

 

 

2

 

 

 

R

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2¼

 

 

 

 

 

 

 

 

 

P1 1 rn

2¼

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

R0

f(Ã)+ ¼

 

 

 

 

 

 

 

f(Ã) cos n(Ã ¡ ')=

 

2¼

=1 Rn

 

 

 

 

 

2¼

 

 

 

 

 

 

 

nP

 

 

 

 

 

R0n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1 r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

f(Ã)½1 + 2 n=1

 

 

 

cos n(Ã ¡ ')¾=

 

 

 

 

 

 

2¼ 0

Rn

(7)

 

 

 

 

 

 

 

 

R

(

Ã

¡

'

=

¾; z

 

 

r

e;

 

z

j

<

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P= R

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2¼

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

f(Ã)½¡1 + 2Re n=0 zn¾=

 

 

 

 

 

 

 

 

 

 

2¼ 0

 

 

 

 

1

 

2¼

 

 

 

 

 

 

 

R

 

 

 

 

1

 

 

 

 

 

 

 

 

P

2¼

 

2

2

)

=

 

 

f(Ã)½¡1 + 2Re

 

 

 

 

 

 

 

 

 

1

 

 

 

 

f(Ã)(R

¡r

2¼ 0

1¡z ¾=

2¼ 0

 

R2+r2¡2Rr cos(á')

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

Формула (7) называется интегралом Пуассона. Она дает решение задачи (3), если f(') - достаточно гладкая.

Определение. Задача называется корректной (по Адамару), если она разрешима при любых начальных данных, принадлежащих к некоторому классу, имеет единственное решение и это решение непрерывно зависит от начальных данных.

Задача называется некорректной, если она разрешима не при любых начальных данных, либо если она имеет неединственное решение, либо если нет непрерывной зависимости решений от начальных данных.

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

5

Примеры некорректных задач:

а) задача Коши с данными на характеристике; б) задача (6) из x6 при t < 0;

и) задача Коши 8

> uxx + uyy = 0;

>

<

> ujy=0 = '0(x);

> u j = ' (x):

: y y=0 1

Соседние файлы в папке 2010-10-28