Lumf_p1&2(2010)
.pdfТеперь мы должны потребовать, чтобы решение (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