№179.11.Ястребов
.pdf21
7.3. Уравнение, не содержащее явно независимую переменную x
Рассмотрим уравнение второго порядка вида |
F( y, y , y ) = 0 |
, не |
||||||||||
|
′ |
′′ |
|
|||||||||
содержащее явно независимую переменную x . |
|
|
|
|
|
|
||||||
Будем предполагать |
y(x) |
строго монотонной функцией. Тогда су- |
||||||||||
ществует обратная функция |
x = x( y) |
, и производные |
y , y |
′′ можно рас- |
||||||||
|
|
′ |
|
|||||||||
сматривать как сложные функции независимой переменной y : |
|
|
||||||||||
|
′ |
|
′ |
|
|
′′ |
′′ |
|
. |
|
|
|
|
y (x) = y (x( y)); |
|
y (x) = y (x( y)) |
|
|
|
|
|||||
Введем новую неизвестную функцию y′ = p( y) . По правилу диф- |
||||||||||||
ференцирования сложной функции |
|
|
|
|
|
|
|
|
||||
′′ |
′ ′ |
|
′ |
|
′ |
′ |
|
′ |
|
|
, |
|
yxx |
= ( yx ) x = ( p( y)) |
x = py ( y) yx (x) |
= p ( y) p( y) |
|
|
так что исходное уравнение второго порядка переходит в уравнение перво-
го порядка относительно новой неизвестной функции p( y) :
|
|
|
′ |
|
= 0 |
. |
|
|
(13) |
|
|
F1( y, p( y), p ( y)) |
|
|
|
|
|
||||
Если найден общий интеграл уравнения (13) |
|
|
|
|
|
|
||||
|
Φ( y, p,C1) = 0 , |
|
|
|
|
|
|
|
||
то, заменяя в нем |
p на y′ , приходим к уравнению первого порядка отно- |
|||||||||
сительно исходной неизвестной функции y(x) : |
|
|
|
|
|
|
|
|||
′ |
. |
|
|
|
|
|
|
|
|
|
Φ( y, y ,C1) = 0 |
|
|
|
|
|
|
|
|
|
|
Таким образом, решение уравнения второго порядка сводится к по- |
||||||||||
следовательному решению двух уравнений первого порядка. |
|
′ |
|
|
||||||
Пример. Рассмотрим уравнение |
y |
′′ |
′2 |
. Полагаем |
y |
= p( y) |
; |
|||
|
|
y = y |
|
|
|
|
|
|||
dp |
|
|
|
|
|
|
|
|
|
|
тогда y′′ = dy |
p. Исходное уравнение |
преобразуется |
|
|
к виду: |
|
|
|
|
|
|
22 |
|
|
|
|
|
|
|
|
|
dp py = p2 |
. Ограничиваясь случаем |
p ≠ 0 , получаем уравнение с раз- |
|||||||||||||
dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
деляющимися переменными: |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
dp |
= dy . |
|||||||||
|
|
|
|
|
p |
|
|
|
y |
||||||
Откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dp |
= |
∫ |
dy |
; |
ln |
|
|
p |
|
= ln |
|
y |
|
+ C . |
|
|
|
|
|
|||||||||||
|
∫ p |
|
y |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Произвольную константу интегрирования удобно записать в виде ln C1 ,
поскольку логарифмическая функция принимает все значения от −∞ до +∞ . Тогда
ln p = ln y + ln C1 = ln( C1 y ) . p = ±C1 y .
Поскольку здесь постоянный множитель при y , записанный в виде C1,
принимает, как и множитель ±C1, все вещественные значения, можно за-
писать: p = C1 y
Возвращаемся к исходной неизвестной функции y :
y′ = C1 y; |
dy |
= C1 y; |
dy |
= C1dx; |
∫ |
dy |
= |
∫ |
C1dx + C2 ; |
|
dx |
|
y |
|
y |
|
|
ln y = C1x + C2
— общий интеграл.
23
8. Линейное уравнение второго порядка
8.1. Основные понятия
Определение. Линейным уравнением второго порядка называется уравнение вида:
|
|
|
|
y′′ + p(x) y′ + q(x) y = |
f (x) |
(14) |
||||
с непрерывными на интервале (a, b) функциями p, q и f . |
|
|||||||||
|
|
Из теоремы 2, приведенной на с. 6, следует, что указанная непрерыв- |
||||||||
ность гарантирует при |
x0 (a, b) |
существование и единственность ре- |
||||||||
шения |
задачи |
Коши |
с |
любыми |
начальными |
данными |
||||
x0 |
|
(a, b); y0 |
′ |
R |
. |
|
|
|
|
|
|
, y0 |
|
|
|
|
|
Определение. Однородным линейным уравнением второго порядка
называется уравнение с нулевой правой частью:
y′′ + p(x) y′ + q(x) y = 0. |
(15) |
8.2. Свойства решений однородного линейного уравнения
Из свойств производной следует, что для любых функций y1, y2 и
любых вещественных чисел λ , μ :
(λ y1 + μ y2 )(n) = λ y1(n) + μ y2(n) .
Обозначим правую часть уравнений (14) и (15) через L( y) :
L( y) = y′′ + p(x) y′ + q(x) y .
Тогда эти уравнения принимают вид L( y) = f (x) и L( y) = 0 соот-
ветственно. При этом
L(λ y1 ± μ y2 ) = λ L( y1) ± μ L( y2 ) ,
|
24 |
|
Теорема 3: Если функции |
y1(x) и |
y2 (x) являются решениями од- |
нородного линейного уравнения |
(15), то функция λ y1(x) + μ y2 ( x) |
|
также является его решением. |
|
|
Доказательство. Пусть L( y1) = 0 |
и L( y2 ) = 0. Тогда |
L(λ y1 ± μ y2 ) = λ L( y1) ± μ L( y2 ) = 0 ± 0 = 0. ▄
8.3. Линейное уравнение с постоянными коэффициентами
Рассмотрим линейное уравнение с постоянными коэффициентами:
y′′ + py′ + qy = f (x) |
(16) |
( p и q – постоянные числа), и соответствующее ему однородное урав-
нение
y′′ + py′ + qy = 0 . |
(17) |
Определение. Характеристическим уравнением, соответствующим дифференциальному уравнению (17), называется алгебраическое квадратное уравнение
z2 + pz + q = 0. |
(18) |
Отметим, что второй производной y′′ дифференциального уравне-
ния соответствует в характеристическом уравнении z2 . Коэффициент при первой производной y′ переходит в коэффициент при первой степени z .
Наконец, коэффициент при y , то есть при производной нулевого порядка,
переходит в свободный член (коэффициент при нулевой степени z ). Примеры. 1. Для линейного однородного уравнения
y′′ − 3y′ − 10 y = 0 соответствующее характеристическое уравнение за-
писывается в виде z2 − 3z − 10 = 0 .
|
|
25 |
2. |
Уравнению |
y′′ + 2 y′ = 0 соответствует характеристическое |
уравнение z2 + 2z = 0. |
||
3. |
Уравнению |
y′′ + 7 y = 0 соответствует характеристическое урав- |
нение z2 + 7 = 0.
Общее решение однородного уравнения (17) можно получить, исходя из корней соответствующего характеристического уравнения (17). Здесь возможны три случая в соответствии с возможным значением его дискриминанта D = p2 − 4q .
А. Случай положительного дискриминанта. Пусть D > 0. В этом случае уравнение (18) имеет два различных вещественных корня z1 и z2 :
z |
= − |
p |
± |
D |
, |
|
|
||||
1,2 |
2 |
|
2 |
|
|
|
|
|
|
что соответствует разложению квадратного трехчлена на множители с вещественными коэффициентами:
z2 + pz + q = (z − z1)(z − z2 ).
Можно проверить, что в этом случае общее решение однородного уравнения (16) имеет вид:
|
|
|
|
y = C ez1x + C |
2 |
ez2 x |
. |
(19) |
|
|
|
|
|
|
1 |
|
|
|
|
Примеры. |
1. |
y′′ + 3y′ − 18 y = 0 ; начальные |
условия: |
||||||
x0 = 0; y(0) = 1; |
y (0) |
= −5 |
. Соответствующее характеристическое |
||||||
|
′ |
|
|
|
уравнение: z2 + 3z − 18 = 0 . Дискриминант D = 81. Корни квадратно-
го уравнения z1 = −6, z2 = 3. Общее решение имеет вид:
y = C1e−6 x + C2e3x .
26
Найдем частное решение для задачи Коши.
Дифференцируем общее решение: y′ = −6C1e−6 x + 3C2e3x . Под-
ставляем начальные условия в y и y′ (учитывая, что e0 = 1):
C1 + C2 = 1,
−6C1 + 3C2 = −5.
Решая эту систему, получаем: C1 = 89 ; C2 = 19 . Соответствующее решение задачи Коши: y = 89 e−6 x + 19 e3x .
2. |
y′′ + 2 y′ = 0 . Характеристическое уравнение z2 + 2z = 0 |
|
z(z + |
2) = 0. Корни квадратного уравнения z1 = − 2, z2 = 0 . Общее |
|
решение имеет вид: y = C e− |
2x + C . |
|
|
1 |
2 |
Б. Случай нулевого дискриминанта. Пусть D = 0. В этом случае характеристическое уравнение имеет один вещественный корень кратности 2:
z1 = z2 = − 2p ,
что соответствует разложению квадратного трехчлена на множители с вещественными коэффициентами:
z2 + pz + q = (z − z1)2 .
Можно проверить, что в этом случае общее решение однородного уравнения (17) имеет вид:
y = C ez1x + C |
2 |
xez2 x |
. |
(20) |
1 |
|
|
|
27
Пример. y′′ − 10 y′ + 25y = 0. Соответствующее характеристиче-
ское уравнение z2 − 10z + 25 = 0. Дискриминант D = 0. Кратный ко-
рень квадратного уравнения z1 = z2 = 5. Общее решение имеет вид:
y = C1e5x + C2 xe5x .
В. Случай отрицательного дискриминанта. Пусть D < 0 . В этом случае уравнение (17) имеет два различных комплексных корня z1 и z2 ,
которые задаются формулой:
z = α ± β i , где α = − p |
, β = |
|
D |
|
|
. |
|
|
|
||||||
|
|
|
|
||||
1,2 |
2 |
|
2 |
|
|
|
|
|
|
|
|
|
К этим значениям можно придти, используя формально выражение для корней, полученное в случае положительного дискриминанта, и помня,
что i обозначает «мнимую единицу» — комплексное число, для которого
i2 = −1:
z = − p ± D |
= − p ± |
|
|
D |
|
(−1) |
= − p ± |
|
|
D |
|
|
−1 = |
|||||||
|
|
|
|
|||||||||||||||||
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
1,2 |
2 |
2 |
|
|
2 |
|
|
2 |
2 |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|||||||||||||
|
|
= |
− p |
± |
|
|
|
D |
|
|
i . |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Можно проверить, что в этом случае общее решение однородного уравнения (17) имеет вид:
y = C eα x cos β x + C |
eα x sin β x |
. |
(21) |
|
1 |
2 |
|
|
|
Примеры. 1. y′′ − 12 y′ + 40 y = 0. Соответствующее характери-
стическое уравнение z2 − 12z + 40 = 0. Дискриминант D = −16 . Кор-
28
ни квадратного уравнения z1,2 |
= 6 ± 2i; α = 6, β = 2 . Общее решение |
||||||
имеет вид: |
y = C e6 x cos 2x + C e6 x sin 2x . |
|
|
|
|||
|
1 |
|
|
2 |
|
|
|
2. y′′ + 8 y = 0. Соответствующее характеристическое уравнение |
|||||||
z2 + 8 = 0. |
Корни |
квадратного уравнения |
z1,2 = ± −8 = 0 ± 8i; |
||||
α = 0, β = |
8 . Общее решение имеет вид: y = C cos 8x + C |
2 |
sin 8x . |
||||
|
|
|
|
1 |
|
|
|
Найдем частное решение для задачи Коши с начальными условиями |
|||||||
x0 = 0; y(0) = 0; |
y (0) = |
8 |
. Дифференцируем общее решение: |
||||
|
|
′ |
|
y′ = − 8C1 sin 8x + |
8C2 cos |
8x |
|
Подставляем начальные условия в |
y и |
y′ |
|
cos0 = 1, sin 0 = 0): |
|
|
|
C = 0 |
|
|
|
|
1 |
8. |
|
|
8C2 = |
|
.
(учитывая, что
Отсюда |
(C1 = 0; C2 = 1). Соответствующее частное решение |
y = sin |
8x . |
8.4. Структура общего решения неоднородного линейного уравнения
Определение. Для линейного уравнения второго порядка
y′′ + p(x) y′ + q(x) y = f ( x) L( y) = f (x)
соответствующим однородным уравнением называется уравнение
y′′ + p(x) y′ + q(x) y = 0 L(Y ) = 0 .
Можно доказать, что относительно структуры общего решения неоднородного линейного уравнения справедлива следующая
29
Теорема 4. Общее решение неоднородного уравнения yон предста-
вимо суммой общего решения соответствующего однородного уравнения yоо и частного решения неоднородного уравнения yчн: yон = yоо + yчн.
Доказательство. Поскольку
L( yоо + yчн) = L( yоо) + L( yчн) = 0 + f (x) = f (x),
то y он действительно является решением неоднородного уравнения. Да-
лее, y он зависит от произвольных постоянных C1 и C2 , поскольку от них зависит функция yоо.
Остается убедиться, что за счет выбора значений этих постоянных можно получить решение задачи Коши с любыми наперед заданными на-
чальными условиями |
y(x0 ) = y0 |
, y (x0 ) = y0 |
yоо |
— об- |
|
|
′ |
′ . Поскольку |
|
щее решение соответствующего однородного уравнения, то можно вы-
брать такие значения постоянных C1 и C2 , при которых
yоо(x0 ) = y0 − yчн(x0 ) ;
и
′ |
′ |
′ |
. |
yоо(x0 ) = y0 |
− yчн(x0 ) |
|
Тогда для функции yон при этих значениях C1 и C2 :
yон(x0 ) = yоо(x0 ) + yчн( x0 ) = ( y0 − yчн( x0 )) + yчн(x0 ) = y0
и |
|
|
|
|
|
|
|
′ |
′ |
′ |
′ |
′ |
′ |
′ |
. ▄ |
yон(x0 ) = yоо(x0 ) + yчн( x0 ) = ( y0 |
− yчн( x0 )) + yчн(x0 ) = y0 |
|
30
9. Метод вариации произвольных постоянных (метод Лагранжа)
Рассмотрим линейное неоднородное уравнение второго порядка с постоянными коэффициентами
y′′ + py′ + qy = f (x) .
Согласно изложенному в п. 8.4, его общее решение представимо в виде yон = yоо + yчн. Общее решение соответствующего однородного урав-
нения записывается в виде
|
|
|
|
|
yоо = C1 y1(x) + C2 y2 (x) . |
(22) |
|||||||
с произвольными постоянными C1, C2 . Функции |
y1(x), |
y2 (x) , в за- |
|||||||||||
висимости от корней характеристического уравнения, имеют вид: |
|||||||||||||
(а) y |
= ez1x , |
y |
2 |
= ez2 x |
— при D > 0, z |
≠ z |
2 |
; |
|
|
|||
1 |
|
|
|
|
1 |
|
|
|
|
|
|
||
(б) y |
= ez1x , |
y |
2 |
= xez1x — при D = 0, z |
|
= z |
2 |
; |
|
||||
1 |
|
|
|
|
1 |
|
|
|
|
|
|||
(в) y |
= eα x cos β x, y |
2 |
= eα x sin β x — при D < 0 , |
|
|||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
z1,2 = α ± β i .
Таким образом, остается найти какое-либо частное решение исход-
ного неоднородного уравнения yчн. Метод вариации произвольных по-
стоянных, предложенный Лагранжем, предполагает отыскание yчн в виде,
аналогичном (22), но уже с переменными множителями при y1 и y2 :
yчн = C1(x) y1(x) + C2 (x) y2 (x) .
Здесь C1(x), C2 (x) — подлежащие определению неизвестные функции.
Вычислим производные частного решения (аргумент x для кратко-
сти опускаем): |
|
|
|
|
|
|
|
′ |
′ |
y1 |
′ |
+ C2 |
′ |
y2 |
′ . |
yчн = C1 |
+ C1 y1 |
|
+ C2 y2 |