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

Lumf_p1&2(2010)

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

Здесь 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

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