УМК
.PDFЗамечание. Если дифференциальное уравнение не содержит явно независимой переменной x , искомой функции y и ее первых (k −1) производных,
т.е. уравнение имеет вид |
F(y(k ), y(k +1),K, y(n ))= 0 , то |
его порядок |
можно |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
понизить на |
(k +1) |
|
единиц, |
применяя сначала подстановку y(k ) = z(x), |
а затем |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
z′ = p(y). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Примеры решения задач |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
ПРИМЕР 2.20. Решить уравнение y ′′= |
|
y′ |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
Решение. |
Это уравнение второго порядка, не содержащее явно x . Делаем |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
замену |
y′ = p(y), |
|
y ′′= p |
dp |
|
и подставляем в уравнение. |
Получим |
|
|
p |
dp |
= |
|
|
p |
|
, |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
dy |
|
|
dy |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
y |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
dp = y− |
1 |
|
|
|
|
|
|
|
|
|
|
p = ∫y− |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy = 2 |
|
+ C . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y′, |
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||
2 dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||
или |
2 |
y |
|
Найденное p |
заменим |
на |
|
тогда |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
y′ = 2 y + C1 . |
|
Разделяем |
|
переменные |
|
|
|
|
|
= dx . |
|
|
|
Интегрируем |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
2 |
y + C1 |
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
dy |
|
|
|
|
y = t 2 |
|
|
|
|
|
|
|
|
2 t dt |
|
|
|
|
|
|
|
|
2 t + C − C |
|
|
|
|
|
C |
|
|
|
|
2 dt |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
x = ∫ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
= ∫ |
|
|
|
|
|
|
|
|
= ∫ |
|
|
1 |
|
|
1 |
dt =∫dt − |
1 |
|
∫ |
|
|
|
|
|
|
|
|
|
. |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2t + C1 |
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
2 y + C1 |
|
dy = 2 t dt |
|
|
2 t + C1 |
|
|
|
2 t + C1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
x = t − |
C1 |
|
|
|
2t + C |
|
+ C |
|
|
|
|
|
Т.к. t = |
|
|
|
|
|
|
|
|
|
x = |
|
|
− 0,5 C ln |
|
2 |
|
|
+ C |
|
|
+ C |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
ln |
|
|
2 |
. |
|
|
|
|
y , |
то |
|
|
y |
|
|
y |
|
2 |
|
– |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
общий интеграл. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
ПРИМЕР 2.21. Найти решение задачи Коши для дифференциального |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
уравненияy3 y ′′= −1 с начальными условиями y(0) =1, y′(0) =1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
Решение. Данное уравнение является уравнением, не содержащим явно |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
независимую |
переменную. |
Сделаем |
|
|
замену y′ = p(y) . |
Тогда, |
учитывая, |
|
|
что |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
p′′ = p′ p , получим уравнение первого порядка y3p′p = −1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
Принимая во внимание, что p′ = |
dp |
|
, получим уравнение первого порядка с |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
разделяющимися переменными |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
dp |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
|
|
|
|
|
|
|
|
dy |
|
|
|
|
p2 |
|
|
|
|
|
|
|
p2 |
|
|
|
|
|
y |
−3+1 |
|
|
|
|
|
C |
|
|
||||||||||||||||||||||||||
|
y3 |
|
|
|
|
p |
= −1, |
|
|
pdp = − |
|
|
|
|
, |
|
∫pdp = ∫ |
− |
|
|
|
|
|
|
, |
|
|
|
|
|
|
= −∫ y−3dy , |
|
|
= − |
|
|
|
|
|
|
|
+ |
|
1 |
, |
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
y3 |
|
|
|
|
y3 |
|
|
2 |
|
|
− |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
3 +1 2 |
|
|
|||||||||||||||||||||||||||||||
|
p2 |
= |
1 |
+ |
C1 |
, p2 = |
1 |
+ C . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
2 |
|
|
|
|
2y2 2 |
|
|
|
|
y2 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Необходимо найти частное решение данного дифференциального уравнения при заданных начальных условиях.
Определим произвольную постоянную С1 учитывая, что y′(0) = 1 и y(0) = 1, то в силу замены y′ = p(y) получим p(1) = 1. Подставим y = 1 и p = 1 в полученное
промежуточное решение p2 = |
1 |
+ C : 1 = |
1 |
+ C , |
C = 0 . |
|||||||
|
|
|||||||||||
|
|
|
y2 |
1 |
1 |
1 |
1 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|||
Таким образом, p2 = |
1 |
или |
p = ± |
1 |
|
. Так как y′ = p(y) , то y′ = ± |
1 |
. |
||||
|
y2 |
|
|
|
y |
|
|
|
|
y |
Получили уравнение первой степени с разделяющимися переменными.
dy |
= ± |
1 |
, |
ydy = ±dx , ∫ ydy = ±∫ dx , |
y2 |
= ±x + |
C2 |
, y2 = ±2x + С2 . |
dx |
y |
|
|
|||||
|
|
2 |
2 |
|
Таким образом, решение примет вид: y = ±C2 ± 2x .
Найдем произвольную постоянную C2 , подставив x = 0 и y = 1 1 = ±С2 ± 0 , C2 = 1.
Окончательно получаем y2 = ±2x +1
или y = ±1 ± 2x – частное решение
2.10 ЛИНЕЙНЫЕ УРАВНЕНИЯ ВЫСШИХ ПОРЯДКОВ
Линейным дифференциальным уравнением с переменными коэффициентами называется уравнение вида
y(n ) + a1(x)y(n −1) + a 2 (x)y(n −2) +K + a n (x)y = f (x). |
(2.16) |
В уравнении (2.16) a1 (x), a 2 (x),K, a n (x) и f (x) заданы и непрерывны в некотором интервале (a, b).
Если f (x) ≠ 0, то уравнение (2.16) называется неоднородным или с правой частью. Если f (x) ≡ 0 , то уравнение называется однородным.
Зная одно частное решение y1 − линейного однородного уравнения, можно с помощью замены y = y1 ∫ z dx понизить порядок уравнения на единицу.
Если y1, y2 ,K, yn − линейно независимые частные решения однородного
уравнения |
|
y(n ) + a1 (x)y(n −1) + K+, a n (x)y = 0 , |
(2.17) |
Уравнение (2.22) является частным случаем уравнения вида (2.17) и его общее решение также имеет структуру (2.18)
y = C1 y1 + C2 y 2 K + C n y n ,
где y1 (x),K, yn (x) − частные решения уравнения (2.22), вронскиан которых не равен нулю; C1 , C2 ,K, Cn -произвольные константы. Такой набор функций называется фундаментальной системой решений. Для линейных однородных (с нулевой правой частью ) уравнений известно как построить фундаментальную систему решений, через которую можно записать общее решение уравнения
(2.22). Для этого выписываем так называемое |
характеристическое уравнение, |
|
соответствующее дифференциальному уравнению (2.22): |
|
|
k n + a1k n −1 + K + a n −1k + a n |
= 0 . |
(2.23) |
Характеристическое уравнение получается из уравнения |
(2.22) формальной |
|
заменой производных неизвестной функции и самой |
функции y на |
|
соответствующие степени параметра k : |
|
|
y(n ) → k n , y(n −1) → k n −1, K, y'→ k, y → 1.
Характеристическое уравнение (2.23) имеет ровно n корней (действительных или комплексных) среди которых могут быть совпадающие (кратные). Напомним, что кратностью корня k j называется степень двучлена
(k − k j )в разложении многочлена на множители.
|
Фундаментальная система решений строится в соответствии с корнями |
||||
характеристического уравнения (2.23) согласно следующей схеме: |
|
|
|||
1. |
каждому простому (кратности |
1) действительному корню k |
ставится в |
||
соответствие функция ekx ; |
|
|
|
|
|
2. |
каждому действительному корню k кратности m ставится в соответствие |
||||
m линейно независимых функций ekx , xekx ,K, xm −1ekx ; |
|
|
|||
3. |
каждой паре комплексно сопряженных корней |
k1 = α + iβ |
и |
k 2 = α − iβ |
|
ставится в соответствие пара функций |
eαx cos(βx), |
eαx sin(βx); |
|
|
|
4. |
каждой паре комплексно сопряженных корней |
k1 = α + iβ |
и |
k 2 = α − iβ |
|
кратности m ставится в соответствие |
2m функций eαx cos(βx), |
eαx sin(βx), |
xeαx cos(βx), xeαx sin(βx),K, xm −1eαx cos(βx), xm −1eαx sin(βx)
Таким образом, по n корням (с учетом кратности) мы построили n линейно независимых решений уравнения (2.22). Согласно (2.18), общее решение уравнения (2.22) выписывается как линейная комбинация этих всех функций.
|
|
Для выделения частного решения необходимо на функцию y и ее |
|||||||||||||
производные |
наложить |
n |
дополнительных |
условий, |
например |
||||||||||
y(x |
0 |
) = y |
00 |
, |
y'(x |
0 |
) = y |
, K, |
y(n −1) (x |
0 |
) = y |
n −1,0 |
(задача Коши). |
|
|
|
|
|
|
10 |
|
|
|
|
|
|
Примеры решения задач
ПРИМЕР 2.23. Найти общее решение уравнения y' '+2y'−3y = 0.
Решение. Заменяя в данном дифференциальном уравнении неизвестную функцию на единицу, ее производные на соответствующие степени параметра k ,
получаем |
|
характеристическое |
|
уравнение |
k 2 + 2k − 3 = 0 . |
|
Корни |
этого |
|||||||||||||||||||
уравнения k1 = −3, k2 = 1 действительны и различны. |
Фундаментальная система |
||||||||||||||||||||||||||
решений |
|
состоит |
из |
функций |
|
y1 |
= e−3x , y 2 |
= e x . |
|
Общее |
|
|
решение, |
||||||||||||||
следовательно, y(x) = C1e −3x |
+ C2 e x . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
ПРИМЕР 2.24. Найти общее решение уравнения y(4) + 2y ′′′+ 5y ′′= 0. |
|||||||||||||||||||||||||
|
|
Решение. Выписываем характеристическое уравнение |
|
k 4 + 2k3 + 5k2 = 0 , |
|||||||||||||||||||||||
его корни: k1 = 0 |
|
кратности 2, и пара простых комплексно сопряженных корней |
|||||||||||||||||||||||||
k3 = −1 + 2 i, |
k 4 = −1 − 2 i. В соответствии со схемой, |
фундаментальная система |
|||||||||||||||||||||||||
решений |
|
|
состоит |
из |
четырех |
функций |
y |
= e |
0 =1, y |
2 |
= xe0 = x, |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
y3 |
= e−x cos (2x), y4 = e−x sin (2x); |
общее |
|
решение |
|
есть |
их |
|
|
линейная |
|||||||||||||||||
комбинация y(x) = C1 + C 2 x + C3e −x cos(2x) + C 4 e −x sin(2x). |
|
|
|
|
|
|
|||||||||||||||||||||
|
|
ПРИМЕР 2.25. Найти частное решение уравнения |
y′′′ − y′′ + 4y′ − 4y = 0, |
||||||||||||||||||||||||
удовлетворяющее начальным условиям |
y(0) = 0, y′(0) = 2, y′′(0) = 10. . |
|
|
|
|||||||||||||||||||||||
|
|
Решение. Найдем сначала общее решение. Выписываем характеристическое |
|||||||||||||||||||||||||
уравнение |
k3 − k2 + 4k − 4 = 0 |
|
или |
|
k 2 (k −1)+ 4(k −1)= 0 , |
|
тогда |
||||||||||||||||||||
(k −1)(k 2 + 4)= 0 . |
Корни характеристического уравнения: |
k1 = 1, |
|
k 2 = 2 i, |
|||||||||||||||||||||||
k |
3 |
= −2 i. |
Общее |
решение |
|
y(x) = C ex |
+ C |
2 |
cos (2x) + C |
3 |
sin (2x). |
|
|
Найдем |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|||
производные |
|
|
y'(x) = C ex − 2C |
2 |
sin(2x) + 2C |
3 |
cos(2x), |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
y''(x) = C ex |
− 4C |
2 |
cos(2x) − 4C |
3 |
sin(2x) . |
Частное |
решение |
находится в |
|||||||||||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
x = 0 , |
y = 0 , |
y′ = 2 , |
|
y′′ = 10 : |
||||||||
соответствии |
с |
|
начальными |
|
условиями |
|
|
||||||||||||||||||||
0 = C1 + C 2 , |
|
|
|
2 = C1 + 2C3 , |
|
|
10 = C1 − 4C2 . |
|
Решая, |
|
|
получаем |
|||||||||||||||
C |
= 2,C |
2 |
= −2,C |
3 |
= 0 , тогда частное решение |
y(x) = 2ex − 2 cos(2x) . |
|
|
|||||||||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ПРИМЕР 2.26. Рассмотрим колебание груза массой m под действием
пружины жесткости k . Пусть |
x(t)− отклонение тела от положения равновесия |
x = 0, x′(t)− скорость тела, x′′(t)− ускорение. |
|
Решение. Запишем для |
этой системы закон Ньютона m x′′ = F . Если |
силами, действующими на тело, являются только возвращающая сила пружины
Fgh = −k x(t), k > 0, сила сопротивления |
Fсопр = −γ x′(t), γ > 0 , то получим |
уравнение линейного осциллятора при |
наличие сопротивления вида |
m x′′ + γ x′ + k x = 0 . Его характеристическое уравнение m k 2 + γ k + k = 0 имеет
корни k1,2 |
= |
− γ ± |
γ2 − 4m k |
|
|
|
. |
||
|
|
|||
|
|
|
2m |
|
Рассмотрим три случая: |
||||
1) |
Если |
γ2 − 4mk > 0 (сила сопротивления движению велика, |
возрастающая сила пружины мала), то корни действительны, различны и оба отрицательны.
Общее решение запишется в виде x(t)= C1ek1t + C2ek 2 t . Это случай так называемого апериодического решения. Точка асимптотически, без колебаний стремится к положению равновесия x = 0 при t → ∞ . Напомним k1, k2 < 0 .
2) |
Если γ2 − 4mk = 0, |
то |
корни |
характеристического уравнения |
||||
действительны и равны k1 = k 2 = − |
|
γ |
= −r < 0 . |
|||||
|
2m |
|||||||
|
|
|
|
|
|
|
|
|
Общее решение имеет вид |
x(t)= (C1 + C2 t)e−r t . |
|||||||
3) |
Наконец γ2 − 4mk < 0 . Корни в этом случае комплексны и сопряжены |
|||||||
|
|
|
|
|
|
|
|
|
k1 = −α + i β, k 2 = −α − i β , α = |
|
γ |
, β = |
|
γ2 − 4mk |
|||
|
|
|
|
. |
||||
|
2m |
|
|
|||||
|
|
|
|
|
|
2m |
||
Общее решение имеет вид |
x(t)= e−α t (C1 sin βt + C2 cosβt). Движение точки |
представляет собой колебания около положения равновесия с затухающей амплитудой. Отметим, что при отсутствии трения k = 0 движение будет
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
+ C2 |
|
k |
|
||||||
|
|
|
|
|||||||||
m t |
m t |
|||||||||||
периодическим x(t)= C1 sin |
|
|
. |
2.12 ЛИНЕЙНЫЕ НЕОДНОРОДНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ n −ГО ПОРЯДКА С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ. МЕТОД ВАРИАЦИИ ПРОИЗВОЛЬНЫХ ПОСТОЯННЫХ
В данном параграфе речь пойдет об уравнениях вида |
|
y(n ) + a1 y(n −1) + K + a n −1 y′ + a n y = f (x). |
(2.24) |
Если функция правой части f (x) не тождественный нуль, такие уравнения называют неоднородными или уравнениями с правой частью.
Коэффициенты a1 , a 2 ,K, a n − известные константы (в данном параграфе
это замечание не принципиально – все выкладки имеют место и |
в случае |
||
зависящих от x коэффициентов). |
|
||
Общее решение уравнения (2.24) представляет собой сумму |
|
||
y(x) = Y* (x) + |
|
|
(2.25) |
Y(x) , |
где Y* (x) - какое-либо частное решение неоднородного уравнения, Y(x) - общее решение однородного уравнения y(n ) + a1 y(n −1) + K + a n −1 y′ + a n y = 0 .
Пусть фундаментальная |
система решений однородного уравнения |
известна: y1 (x), y2 (x),K, yn (x) . |
Тогда общее решение однородного уравнения |
можно выписать Y(x) = C1y1 (x) + C2 y2 (x) + K + Cn yn (x). Таким образом, для построения общего линейное однородного уравнения решения необходимо найти какое-либо частное решение Y* (x).
Для нахождения частного решения неоднородного уравнения можно пользоваться методом вариации произвольных постоянных. Он заключается в
следующем. Решение уравнения (2.24) находится в виде линейной комбинации функций y1 (x), y2 (x),K, yn (x) с коэффициентами, уже зависящими от x :
Y* (x) = C (x)y (x) + C |
2 |
(x)y |
2 |
(x) + K + C |
n |
(x)y |
n |
(x). |
||||||||
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
||||
Производные |
от |
|
|
функций |
|
C1 (x), C2 (x),K, Cn (x) определяются из |
||||||||||
линейной алгебраической системы вида |
|
|
|
|
|
|
|
|||||||||
|
C′ |
(x)y |
+ C′ |
(x)y |
2 |
+ K + C′ |
(x)y |
n |
= 0, |
|||||||
|
1 |
|
|
1 |
|
2 |
|
|
|
n |
|
|
|
|
||
C′ |
(x)y′ |
+ C′ |
(x)y′ |
+ K + C′ |
(x)y′ |
= 0, |
||||||||||
|
1 |
|
|
1 |
|
2 |
|
2 |
|
|
n |
|
|
|
n |
(2.26) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M |
|
|
|
|
|
|
|
|
−1) + C′ (x)y(n −1) |
|
|
|
(x)y(n −1) = f (x), |
|||||||||||
C′ |
(x)y(n |
+ K + C′ |
||||||||||||||
1 |
1 |
|
|
|
2 |
|
2 |
|
|
|
n |
|
|
|
n |
|
y1 = cos 3x |
и |
|
y2 |
= sin 3x . Частное |
решение данного |
уравнения по методу |
|||||||||||||||||||||||||
вариации будет иметь вид Y* = C (x)cos 3x + C |
2 |
(x)sin 3x . |
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
C′ |
(x)cos 3x + C′ |
sin 3x = 0 |
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Составим систему |
− 3C |
(x)sin 3x + 3C′ cos 3x = |
9 . |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
sin 3x |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Вычислим главный определитель системы |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
cos 3x |
sin 3x |
|
|
= 3 ≠ 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
= W = |
|
|
|
|
система |
|
имеет |
|
единственное |
||||||||||||||||||||||
|
|
|
|
|
− 3sin x |
3cos 3x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
решение, по формулам Крамера имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
C′ = |
|
|
0 |
sin 3x |
|
|
= −9; |
C′ |
= |
|
|
cos 3x |
|
0 |
|
|
|
= −9ctg 3x . |
|||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
9 |
3cos 3x |
|
|
− 3sin 3x |
|
9 |
|
|
|
||||||||||||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
sin 3x |
|
|
|
|
|
|
|
|
|
sin 3x |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
C′ = |
|
C′ |
|
|
|
C′ = |
C′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
1 |
|
= −3, |
|
2 |
= 3ctg3x . |
Интегрируя последние равенства, |
||||||||||||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||||
1 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
найдем C1 (x) и |
|
C2 (x). |
C1 |
= −3x, |
C2 |
= ln |
|
sin 3x |
|
. |
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Запишем частное решение неоднородного уравнения
Y* = −3x cos 3x + ln sin 3x (sin 3x), тогда общее решение будет иметь вид y = C1 cos 3x + C2 sin 3x − 3x cos 3x + ln sin 3x (sin 3x).
2.13 НАХОЖДЕНИЕ ЧАСТНОГО РЕШЕНИЯ ЛИНЕЙНОГО НЕОДНОРОДНОГО УРАВНЕНИЯ С ПОСТОЯННЫМИ
КОЭФФИЦИЕНТАМИ СО СПИЦИАЛЬНОЙ ПРАВОЙ ЧАСТЬЮ. МЕТОД НЕОПРЕДЕЛЕННЫХ КОЭФИЦИЕНТОВ
Частное решение уравнения n −го порядка |
|
y(n ) + a1 y(n−1) + a 2 y(n −2) +K + a n y = f (x), |
(2.30) |
где |
|
f (x) = eα x (Pn (x)cos β x + Q m (x)sin βx), |
(2.31) |
а a1 , a 2 ,K, a n R следует искать в виде |
|
Y* (x) = x r eα x (PS (x)cos β x + QS (x)sin βx). |
(2.32) |