Материал: Основы высшей математики. Алейников С.М., Горяйнов В.В

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

k =

2 + (2)2 4 2

= 2 + −4 = 2 + −1 4 = 2 + 4 1 = 2 + 2i =1+i .

 

1

2

2

2

2

2

 

Как известно, если квадратный трехчлен с действительными коэффициентами имеет комплексный корень k1 =α +iβ , то он имеет и сопряженный ко-

рень k2 =α iβ . Поэтому второй корень характеристического уравнения мож-

но не вычислять, а написать сразу

k2 =1i .

Корни характеристического уравнения получились комплексные. Значит, по формуле (4.42) получаем общее решение данного ЛОДУ

y = eαx (C1 cos βx +C2 sin βx)= ex (C1 cos x +C2 sin x).

Таким образом, нахождение общего решения ЛОДУ второго порядка с постоянными коэффициентами (4.37) сводится к нахождению корней характеристического уравнения (4.39) и использованию формул (4.40), (4.41), (4.42) общего решения уравнения (не прибегая к вычислению интегралов).

4.4.4. Исследование дифференциального уравнения свободных колебаний

Ранее, в п. 4.2, было выведено уравнение свободных колебаний (4.3):

y′′(t )+ 2λy(t )+ω2 y(t )= 0 .

ДУ (4.3) является линейным однородным дифференциальным уравнени-

ем второго порядка. Найдем решение этого уравнения.

 

I. Предположим, что отсутствует сила трения, т.е. μ = 0

λ = 0 . Та-

ким образом, уравнение (4.3) принимает вид

 

y′′(t )+ω2 y(t )= 0 .

(4.43)

Для решения ДУ (4.43) составим характеристическое уравнение

p2 +ω2 = 0 .

Его корнями являются комплексные числа p1 =iω и p2 = −iω . Следовательно, применяя формулу (4.42), получим решение ДУ (4.43)

141

y =C1 cosωt +C2 sinωt ,

(4.44)

где C1 и C2 – произвольные постоянные.

Запишем функцию (4.44) в более компактной форме. Для этого введем новые постоянные A и ϕ0 , а C1 и C2 будут связаны с A и ϕ0 соотношениями:

C1 = Asinϕ0 , C2 = Acosϕ0 .

(4.45)

Тогда постоянные A и ϕ0 через C1 и C2 определяются следующим обра-

зом:

 

 

 

= arctg C1 .

 

A = C2 +C2

и ϕ

0

 

1

2

 

C2

 

 

 

 

 

 

Подставляя значения C1 и C2

из (4.45) в уравнение (4.44), получим

 

y = Asinϕ0 cosωt + Acosϕ0 sinωt

 

или

 

 

 

 

 

y = Asin (ωt +ϕ0 ).

(4.46)

Колебания, описываемые уравнением (4.46), называются гармоническими. Интегральными кривыми в этом случае являются синусоиды. График функции (4.46) изображен на рис. 4.5.

х

ϕ0

А

 

0

t

T = 2ωπ

Рис. 4.5. График функции (4.46)

Периодом колебаний называется промежуток времени Т, за который ар-

гумент синуса изменяется на 2π (в данном случае T = 2ωπ ).

142

Частотой колебаний называется число колебаний за время 2π (в нашем случае частота равна ω).

Амплитудой колебаний называется величина наибольшего отклонения от положения равновесия (в данном случае амплитудой является величина А).

Число ϕ0 называется начальной фазой.

II. Теперь рассмотрим решение ДУ (4.3) для случая, когда сила трения присутствует, т.е. μ 0 λ 0 . Характеристическое уравнение для (4.3)

имеет вид

p2 + 2λ p +ω2 = 0 .

(4.47)

Корнями уравнения (4.47) являются числа

p1,2 = −λ ± λ2 ω2 = −λ ± λ ,

где λ < λ .

В зависимости от соотношения между λ и ω могут возникнуть три слу-

чая:

1. При λ =ω имеем p1 = p2 = −λ ( p1,2 < 0 ). Таким образом, применив формулу (4.41), запишем решение ДУ (4.3) в виде

y = eλt (C +C

t ).

(4.48)

1

2

 

 

Из уравнения (4.48) следует, что отклонение у при любых начальных условиях асимптотически стремится к нулю, если t →∞. В этом случае колебаний груза не будет, так как сила сопротивления велика по сравнению с коэффициентом упругости пружины k и в зависимости от начальных условий может наблюдаться один из видов апериодического движения (рис. 4.6).

2. При λ >ω имеем p1 = −λ + λ и p2 = −λ λ ( p1,2 < 0 ). Значит, используя формулу (4.40), получим

y =C e(λ+λ)t +C e(λλ)t .

(4.49)

1

2

 

Здесь отклонение у также как и в первом случае будет асимптотически стремится к нулю при t → ∞, но, в отличие от предыдущего случая, гораздо быстрее.

143

у

0

t

а)

у

0

t

б)

t

в)

Рис. 4.6. Виды апериодического движения

3. При λ <ω имеем комплексные корни p = −λ ±i

ω2 λ2 = −λ ±ωi .

 

 

1,2

 

Тогда, используя формулу (4.42), получим

 

 

 

y = eλt (C cosωt +C

2

sinωt )

(4.50)

1

 

 

или

 

 

 

y = Aeλt sin (ωt +ϕ0 ).

(4.51)

Здесь в качестве амплитуды нужно рассматривать величину Aeλt , зави-

сящую от времени. Так как « λ < 0 », то амплитуда Aeλt будет стремиться к нулю при t → ∞, т.е. в этом случае мы имеем дело с затухающими колеба-

144

ниями, период которых равен T =

2π

~

 

 

2

2

 

(ω

=

ω

 

λ

– частота затухающих ко-

~

 

 

ω

 

 

 

 

 

 

лебаний). График затухающих колебаний изображен на рис. 4.7.

у

Амплитуда

0

t

y = Aeλt sin (ωt +ϕ0 )

Рис. 4.7. Затухающие колебания

145

Источник: https://studfile.net/preview/16565699/