
Lumf_p1&2(2010)
.pdfЗдесь Dn – произвольная постоянная. Найденным собственным значениям соответствуют решения уравнения для функции T :
Tn(t) = An cos |
n |
at + Bn sin |
n |
at |
(79) |
|
|
||||
|
l |
l |
|
Здесь An и Bn – произвольные постоянные. Таким образом, мы нашли частные решения исходного уравнения колебаний струны:
un(x; t) = Xn(t)Tn(t) |
|
|
(80) |
||||
или |
|
|
|
|
|
|
|
un(x; t) = An cos |
n |
at + Bn sin |
n |
at sin |
n |
x |
(81) |
l |
l |
l |
Очевидно, что сумма частных решений также будет удовлетворять исходному уравнению и граничным условиям:
1 |
lnat + Bn sin lnat sin lnx (82) |
|||||
u(x; t) = n=1 An cos |
||||||
X |
|
|
|
|
|
|
Неизвестные константы надо определить из начальных условий:
u(x; 0) = '(x); ut(x; 0) = (x): |
(83) |
41

Т.е.,
1 |
An sin nx = '(x) |
(84) |
||||||
X |
|
|
|
|
|
|
|
|
n=1 |
|
|
|
l |
|
|||
|
|
|
|
|
|
|
|
|
1 B na sin nx = (x) |
(85) |
|||||||
X |
n |
|
|
|
|
|
|
|
n=1 |
l |
|
|
|
l |
|
||
|
|
|
|
|
|
|
|
Формулы (84) и (85) представляют из себя разложения функций '(x) и (x) в ряд Фурье. Для нахождения неизвестных констант умножим левую и правую части уравнения (84) на sin lmx и проинтегрируем их по dx от 0 до l:
1 |
l |
|
|
|
|
|
l |
|
||
|
|
|
|
|
|
|
|
|
|
|
n=1 An Z |
sin |
n |
x sin |
m |
x dx = Z |
'(x) sin |
m |
x dx |
(86) |
|
l |
l |
l |
||||||||
X |
0 |
0 |
|
|
|
|
Для вычисления интеграла в левой части последнего равенства
42

воспользуемся тригонометрической формулой
sin sin = 12(cos( ) cos( + ))
l
Z
sin lnx sin lmx dx =
0
l
Z
= 12
0
cos (n l |
m)x dx 2 Z |
l |
l |
x dx = |
||||
cos |
||||||||
|
|
|
1 |
|
(n + m) |
|
||
|
|
|
|
|
0 |
|
|
|
|
1 |
|
l |
(n |
m) |
l |
|
1 |
l |
|
|
|
(n + m) |
l |
|
||||
= |
|
|
|
sin |
|
|
x 0 |
|
|
|
|
|
|
sin |
|
|
x 0 |
= |
|
2 |
(n m) |
|
l |
2 |
(n + m) |
l |
|
||||||||||||
|
|
|
|
|
|
= 0; |
если |
|
|
|
m = n: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
1 |
|
если |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
l; |
m = n: |
||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
43
Таким образом,
Z |
l |
|
l |
|
x dx = mn2 |
(87) |
|||||||||||||
sin lnx sin |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
l |
|
||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставляя (87) в (86), получаем |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
2 |
|
l |
|
|
m |
|
||||||||||||
|
|
|
|
|
|
||||||||||||||
|
Am = |
|
Z |
'(x) sin |
|
|
|
x dx |
(88) |
||||||||||
|
l |
|
l |
||||||||||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Аналогично для Bm получаем |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
l |
|
(x) sin |
|
l x dx |
(89) |
||||||
Bm = ma Z |
|
|
|||||||||||||||||
|
2 |
|
|
|
|
|
|
|
|
m |
|
||||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
44
|
Физическая интерпретация решения |
|
|
|
|
|
|
|
|
|
|||||||||||
|
Перепишем функцию un(x) в другом |
виде |
|
|
|
|
|||||||||||||||
u |
(x; t) = |
|
n |
n |
|
|
|
sin |
|
n |
x = |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
An cos l at + Bn sin |
|
|
at |
|
|||||||||||||||||
|
n |
|
l |
n |
|
l |
n |
||||||||||||||
|
|
|
|
|
|
|
= Cn sin |
|
|
x cos |
|
|
a(t + n) (90) |
||||||||
|
|
|
|
|
|
l |
|
|
|||||||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l |
||||||
|
|
|
|
|
|
n |
|
|
|
|
|
|
Bn |
||||||||
|
|
|
Cn = qAn2 + Bn2; |
|
|
arctg |
|||||||||||||||
|
|
|
|
|
a n = |
|
|
||||||||||||||
|
|
|
l |
|
An |
Таким образом, каждая определенная точка струны с координатой x0 колеблется по закону
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
u |
n |
(x ; t) = C |
n |
sin |
n |
|
x |
|
cos |
n |
a(t + ) |
(91) |
||||||||||
|
|
|
|
|
|||||||||||||||||||
|
|
0 |
|
|
|
|
|
l 0 |
|
|
|
l |
|
|
n |
|
|||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z (t) = Z |
n |
cos |
n |
a(t + |
n |
) |
|
(92) |
|||||||||||
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
n |
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
|||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
Zn = Cn sin |
n |
x |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
|
|
45

– амплитуда колебаний. Т.е. все точки струны колеблются в одинаковой фазе, но с разными амплитудами. Такое движение струны представляет из себя стоячую волну. Точки, у которых амплитуда колебаний равна нулю называются узлами стоячей волны, точки у которых амплитуда максимальная – пучности стоячей волны. Частоты колебаний всех точек струны одинаковы и равны
|
|
|
|
! = |
n |
a |
(93) |
|
|||
n |
l |
|
|
|
|
и носят название собственных частот колебаний струны.
Самая низкая частота (n = 1) или самый низкий тон называется
основным тоном струны: |
|
|
|
|
|
|
|
|
|
|
|
|
!1 = |
|
a; |
(94) |
|
|
|
||||
|
|
l |
|
|
остальные тона, соответствующие частотам, кратным !1, называются обертонами.
46

Вынужденные колебания струны
Метод разделения переменных позволяет решить задачу о вынужденных колебаниях струны, уравнение которых имеет вид:
a2uxx + f(x; t) = utt |
(95) |
|
Начальные и краевые условия: |
|
|
u(0; t) = 0; |
u(l; t) = 0 |
(96) |
u(x; 0) = '(x); |
ut(x; 0) = (x): |
(97) |
Будем искать решение в виде суммы двух функций u(x; t) = v(x; t) + w(x; t)
При этом функция v(x; t) будет решением однородного уравнения
a2vxx = vtt
с начальными и краевыми условиями
v(0; t) = 0; v(l; t) = 0 |
(98) |
47
v(x; 0) = '(x); vt(x; 0) = (x): |
(99) |
а функция w(x; t) должна удовлетворять неоднородному уравнению
a2wxx + f(x; t) = wtt |
(100) |
с нулевыми начальными и граничными условиями w(0; t) = w(l; t) = 0
w(x; 0) = wt(x; 0) = 0
Функция v(x; t) описывает свободные колебания струны, происходящие вследствие начального возмущения, w(x; t) – вынужденные колебания без начальных возмущений. Решение v(x; t) нам уже известно. w(x; t) будем искать в виде ряда по собственным функциям однородной задачи:
w(x; t) = |
1 |
(t) sin kx |
(101) |
||
|
X |
k |
|
|
|
|
k=1 |
l |
|
||
|
|
|
|||
|
|
|
|
|
48

Очевидно, что при таком выборе решения граничные условия удовлетворяются автоматически. Для того, чтобы удовлетворить начальным условиям, надо потребовать
k(0) = k0 (0) = 0
Перепишем уравнение (100) в виде
|
wtt a2wxx = f(x; t) |
|
(102) |
||||
и подставляя сюда (101) получаем |
|
|
|
||||
1 |
" 00(t) + |
|
2k2a2 |
k |
x = f(x; t) |
(103) |
|
X |
l2 |
k(t)#sin |
l |
||||
|
|
|
|
|
|
|
k=1
Разлагая функцию f(x; t) в ряд по той же системе функций, по-
лучим
1
f(x; t) = k(t) sin kx (104)
X
l
k=1
49
где
k(t) = l |
Z |
l |
(105) |
|||
f(x; t) sin lkx dx |
||||||
2 |
|
|
|
|
||
|
|
0 |
|
|
|
|
Подставляя (104) в (102) и приравнивая коэффициенты при одинаковых собственных функциях получим обыкновенные дифференциальные уравнения для нахождения неизвестных функций
k(t):
|
00(t) + |
2k2a2 |
|
|
(t) = (t) |
(106) |
|
k |
|||||
|
k |
l2 |
|
k |
|
|
|
|
|
|
|
|
Общее решение этого неоднородного уравнения представляется в виде сумму общего решения однородного уравнения и частного решения неоднородного уравнения:
|
(t) = A |
|
cos |
ka |
t + B |
|
sin |
ka |
t + |
но(t) |
|
(107) |
k |
|
k |
|
|||||||||
k |
|
|
l |
|
l |
k |
|
|
||||
|
|
|
|
|
|
|
|
|
50