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

Lumf_p1&2(2010)

.pdf
Скачиваний:
38
Добавлен:
21.05.2015
Размер:
1.03 Mб
Скачать

Теперь мы должны потребовать, чтобы решение (59) удовлетворяло начальным условиям:

u(x; 0) = f1(x) + f2(x) = '(x)

(60)

ut(x; 0) = af10(x) af20(x) = (x)

(61)

Проинтегрируем (61):

 

 

 

 

x

(z)dz + C

 

f1(x) f2(x) = a Z

 

1

 

 

 

 

 

x0

 

 

В результате получаем систему для нахождения f1 и f2:

 

f1(x) + f2(x) = '(x)

 

(62)

 

 

x

(z)dz + C

(63)

f1(x) f2(x) = a Z

1

 

 

 

 

 

x0

 

 

31

Складывая и вычитая, находим:

 

 

f1(x) = 2'(x) +

 

x

(z)dz +

2

 

 

 

(64)

 

 

2a Z

 

 

 

 

 

1

 

 

 

1

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

f2(x) = 2'(x)

 

x

(z)dz

2

 

 

 

(65)

 

 

2a Z

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

Подставляем найденные f1 и f2 в (59):

 

 

 

 

 

 

 

 

u(x; t) =

1

('(x + at) + '(x at))+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

+

1

 

2 x+at

(z)dz

 

x at

(z)dz3

(66)

 

 

2a

 

 

Z

 

 

 

 

 

 

 

 

Z

 

 

7

 

 

 

 

 

 

 

 

 

 

6 x0

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

5

 

32

u(x; t) = 2

('(x + at) + '(x at)) +

2a

x+at

(z)dz

(67)

Z

1

 

1

 

 

 

 

 

 

x at

 

 

 

Формула (67) – формула Даламбера. Она была получена в предположении существования решения рассматриваемой задачи. Любое решение задачи Коши для бесконечной струны дается формулой Даламбера, что доказывает единственность решения. Сам метод вывода формулы Даламбера доказывает существование решения.

Полученное решение с физической точки зрения представляет собой процесс распространения начального отклонения и начальной скорости. Функция f(x at) представляет собой неизменный профиль f(x), перемещающийся в положительном направлении оси x со скоростью a распространяющаяся или бегущая волна; функция f(x+at) – волна, бегущая в отрицательном направлении оси x. Таким образом, общее решение задачи Коши для бесконечной струны представляет собой суперпозицию двух волн, одна из

33

которых распространяется направо со скоростью a, другая налево с той же скоростью.

Для исследования решения (67) удобно ввести плоскость состояний или фазовую плоскость (x; t). Рассмотрим фиксированную точку M (x0; t0) и проведем через нее характеристики x

at = C1 = x0 at0 и x+at = C2 = x0 + at0. Очевидно, что эти характеристики

пересекут ось x в точках x1 = x0 at0 и x2 = x0 +at0. Найдем значение функ-

ции u(x; t) в точке M:

u(x0; t0) = f1(x0 at0) + f2(x0 + at0) = f1(x1) + f2(x2) (68)

Т.о., отклонение струны в точке M определяется начальным отклонением в вершинах характеристического треугольника PQM

34

и значением начальной скорости на стороне PQ:

(69)

u(M) = 2

('(P ) + '(Q)) +

2a Z

(z) dz

1

 

1

 

 

 

P Q

3.4 Метод разделения переменных

Метод разделения переменных носит также название метода Фурье и является наиболее распространенным методом решения уравнений с частными производными. Рассмотрим его на примере струны с закрепленными концами. Уравнение колебаний

utt = a2uxx

(70)

Граничные условия

 

 

u(0; t) = 0;

u(l; t) = 0

(71)

Начальные условия

 

 

u(x; 0) = '(x);

ut(x; 0) = (x):

(72)

35

Будем искать решение в виде произведения функции зависящей только от x и только от t:

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

(73)

Подставляя (73) в (70) получаем

 

X00T =

1

T 00X

 

 

 

 

a2

 

Разделим левую и правую часть нашего равенства на произведение XT :

X00

=

1

 

T 00

(74)

X

a2 T

 

 

В (74) левая часть является функцией только x, правая часть – только t, причем оно должно выполняться во всей области значений переменных. Это возможно только в том случае если правая и левая часть равны некой константе:

X00

 

1 T 00

(75)

 

=

 

 

 

=

X

a2

T

36

В результате получаем ОДУ для нахождения неизвестных функций X и T :

X00 + X = 0

(76)

T 00 + a2 T = 0

(77)

Из граничных условий

u(0; t) = X(0)T (t) = 0 ) X(0) = 0 u(l; t) = X(l)T (t) = 0 ) X(l) = 0

Таким образом для нахождения функции X(x) мы получили задачу на собственные функции и собственные значения (задачу Штурма-Лиувилля):

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

X00 + X = 0

(78)

X(0) = X(l) = 0

37

а также соответствующие им решения – собственные функции. Рассмотрим возможные значения параметра .

1. < 0 В этом случае общее решение уравнения (78) ищем в виде:

X = Cex

Тогда:

X0 = C ex

X00 = C 2ex

Подставляем в (78):

C 2ex + Cex = 0

Отсюда

2 + p= 0

=

И в результате общее решение имеет вид

p p

X(x) = C1e x + C2e x

38

Из граничных условий

X(0) = C1 + C2 = 0

p p

X(l) = C1e l + C2e l = C1e + C2e = 0

Из первого уравнения находим C1 = C2, подставляем во второе

C1(e e ) = 0

Отсюда получаем C1 = 0, тогда и C2 = 0.

Таким образом, мы показали, что при < 0 задача не имеет нетривиальных решений.

2.= 0. В этом случае тоже не возникает нетривиальных решений.

Упражнение. Доказать, что при = 0 рассматриваемая задача не имеет нетривиальных решений.

3.> 0. В этом случае общее решение имеет вид

p p

X(x) = D1 cos x + D2 sin x

39

Из граничных условий находим

X(0) = D1

= 0

p

 

 

 

X(l) = D2 sin

l = 0

Отсюда

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

sin l = 0

 

p

 

=

n

 

 

 

 

l

 

 

 

 

 

 

 

 

где n любое целое число.

Таким образом нетривиальные решения нашей задачи возможны лишь при значениях

n = n 2 l

Таким образом, мы нашли собственные значения, им будут соответствовать собственные функции

n Xn(x) = Dn sin l x:

40

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