Теорема 3.13. Если ряды (3.59), (3.60) сходятся при |x−x0| < r, где r > 0, то для любых чисел y0, y00 , ..., y0(n−1) решение y(x) уравнения (3.14), удовлетворяющее начальным условиям
y(x0) = y0, y0(x0) = y00 , ..., y(n−1)(x0) = y0(n−1),
разлагается в степенной ряд
∞ |
|
X |
(3.61) |
y(x) = Ak(x − x0)k, |
k=0
который сходится при |x − x0| < r.
Эту теорему часто используют для нахождения решений уравнения (3.14). Для определения коэффициентов Ak в разложении (3.61) используют следующие 2 метода. Заметим, во-первых, что вследствие формулы Тейлора
A |
0 |
= y(x |
), A |
1 |
= y0(x |
), A |
2 |
= |
y00(x0) |
, ..., A |
n−1 |
= |
y(n−1)(x0) |
. |
|
|
|||||||||||||
|
0 |
|
0 |
|
2 |
|
|
(n − 1)! |
||||||
|
|
|
|
|
|
|
|
|
|
|
||||
Поэтому коэффициенты A0, ..., An−1 определяются сразу из начальных условий:
|
y00 |
|
y(n−1) |
|
A0 = y0, A1 = y00 , A2 = |
0 |
, ..., An−1 = |
0 |
. |
2 |
|
|||
|
|
(n − 1)! |
||
Первый метод состоит в следующем: ряд (3.61) подставляется в уравнение (3.14), приравниваются коэффициенты при одинаковых степенях (x− −x0) и из полученных уравнений последовательно находятся коэффициенты An, An+1, ... .
Пример. Найдем разложение в степенной ряд в окрестности x0 = 0 (с точностью до x4) общего решения уравнения
y00 − |
|
2 |
y0 |
+ |
|
1 |
y = 0. |
(3.62) |
1 |
x |
1 |
x |
|||||
|
|
− |
|
|
− |
|
||
Известно, что для нахождения общего решения достаточно отыскать ф.с.р. В данном случае будем искать ф.с.р. y1(x), y2(x) как решения, удовлетворяющие следующим начальным условиям:
(y10 |
(0) |
= 0, |
(y20 |
(0) |
= 1. |
|
|
y1 |
(0) |
= 1, |
y2 |
(0) |
= 0, |
1 |
0 |
|
|
|
|
|
Так как векторы 0 , |
1 линейно независимы, построенные решения дей- |
|||||
ствительно будут составлять ф.с.р.
65
Функции |
2 |
|
1 |
|
p1(x) = − |
, p2(x) = − |
|||
|
|
|||
1 − x |
1 − x |
разлагаются в ряды по степеням x, сходящиеся при |x| < 1. По теореме 3.13 решения y1(x), y2(x) будут разлагаться в ряды по степеням x, также сходящиеся при |x| < 1:
y1(x) = |
∞ |
(3.63) |
Ak(1)xk, |
||
|
X |
|
|
k=0 |
|
y2(x) = |
∞ |
(3.64) |
Ak(2)xk. |
||
|
X |
|
k=0
Из начальных условий находим:
A(1)0 = 1, A(1)1 = 0, A(2)0 = 0, A(2)1 = 1.
Продифференцируем ряды (3.63), (3.64) и подставим их в уравнение
(1 − x)y00 − 2y0 + y = 0,
равносильное уравнению (3.62) при x 6= 1:
∞ |
∞ |
∞ |
(1−x) |
k(k−1)Ak(i)xk−2 −2 |
kAk(i)xk−1 + Ak(i)xk = 0, i = 1, 2. (3.65) |
X |
X |
X |
k=2 |
k=1 |
k=0 |
Приравниваем коэффициент при x0 в левой части равенства (3.65) нулю:
2A(2i) − 2A(1i) + A(0i) = 0.
Отсюда
1 |
|
|
|
|
|
|
||||
A2(1) = − |
|
|
, A2(2) = |
1. |
||||||
2 |
||||||||||
Приравниваем нулю коэффициент при x1: |
|
|
|
|
||||||
6A3(i) − 2A2(i) − 4A2(i) + A1(i) = 6A3(i) − 6A2(i) + A1(i) = 0. |
||||||||||
Отсюда |
|
|
|
|
||||||
1 |
|
|
|
5 |
|
|||||
A3(1) = − |
|
, A3(2) = |
|
|
. |
|||||
2 |
6 |
|||||||||
Наконец, приравниваем нулю коэффициент при x2: |
||||||||||
12A4(i) − 6A3(i) − 6A3(i) + A2(i) = 12A4(i) − 12A3(i) + A2(i) = 0. |
||||||||||
Отсюда |
|
|
|
|
||||||
11 |
|
3 |
||||||||
A4(1) = − |
|
, A4(1) = |
|
|
. |
|||||
24 |
4 |
|||||||||
66
Таким образом, разложения для функций y1(x), y2(x) имеют вид:
|
x2 |
x3 |
|
11 |
4 |
|
|||||
y1(x) = 1 − |
|
− |
|
− |
|
|
x |
|
|
+ ..., |
|
2 |
2 |
24 |
|
|
|||||||
y2(x) = x + x2 + |
5x3 |
+ |
3x4 |
|
+ ..., |
||||||
6 |
|
||||||||||
|
|
|
|
4 |
|
|
|
||||
где многоточиями обозначены члены степенных рядов, содержащие x5, x6,
... .
Теперь получаем требуемое разложение общего решения уравнения
(3.62): |
+ c2x + c2 − 21 |
x2 + |
562 |
− 21 |
x3 + |
|
342 |
− |
241 |
x4 |
|
||||
y(x) = c1 |
+ ... |
||||||||||||||
|
|
c |
|
|
|
c |
|
c |
|
|
c |
|
11c |
|
|
Второй метод нахождения коэффициентов в разложении (3.61) состоит |
|||||||||||||||
в последовательном применении равенств |
|
|
|
|
|
|
|
||||||||
|
|
|
y(k)(x0) |
|
|
|
|
|
|
|
|
|
|||
|
Ak |
= |
|
|
, |
k = 0, 1, ..., |
|
|
|
|
|
|
|||
|
k! |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
вытекающих из формулы Тейлора. Как было показано, коэффициенты A0,
..., An−1 ряда (3.61) определяются из начальных условий. Так как искомая функция y(x) – решение уравнения (3.14), получаем равенство
y(n)(x0) = q(x0) − p1(x0)y(n−1)(x0) − ... − pn(x0)y(x0),
из которого находим y(n)(x0), а следовательно, и An. Дифференцируем уравнение (3.14) по x и подставляем x = x0:
y(n+1)(x0) = q0(x0) − p01(x0)y(n−1)(x0) − p1(x0)y(n)(x0) − ...
... − p0n(x0)y(x0) − pn(x0)y0(x0).
Все числа, входящие в правую часть этого равенства, уже известны. Из него находим y(n+1)(x0) и An+1. Продолжая дифференцировать уравнение (3.14), последовательно определим An+2, ... .
Пример. Найдем разложение в степенной ряд в окрестности x0 = 0 (с точностью до x5) общего решения уравнения
y00 − xy = 0. |
(3.66) |
Как и в предыдущем примере, будем искать ф.с.р. y1(x), y2(x) как решения, удовлетворяющие начальным условиям
(y10 |
(0) |
= 0, |
(y20 |
(0) |
= 1. |
y1 |
(0) |
= 1, |
y2 |
(0) |
= 0, |
67
Функции p1(x) = 0, p2(x) = x разлагаются в ряды по степеням x, сходящиеся при |x| < ∞, т. е. при |x| < R для любого R > 0. Поэтому по теореме 3.13 функции y1(x), y2(x) разлагаются в ряды, сходящиеся при |x| < R для любого R > 0, т. е. в ряды, сходящиеся при |x| < ∞.
Запишем представления (3.63), (3.64) для функций y1(x), y2(x). Из начальных условий следует, что
A(1)0 = 1, A(1)1 = 0, A(2)0 = 0, A(2)1 = 1.
Из (3.66) получаем при x = 0
y00(0) = 0.
Отсюда A(1)2 = 0, A(2)2 = 0. Дифференцируем равенство
|
|
y00 = xy |
|
|
|
|
|
|
||||
по x: |
|
y000 = xy0 |
|
|
|
|
(3.67) |
|||||
|
|
+ y. |
|
|
|
|||||||
Подставим в (3.67) x = 0: |
|
|
|
|
|
|
|
|
|
|
||
|
|
y000(0) = y(0). |
|
|
|
|
|
|
||||
Поэтому y000 |
(0) = 1, y000 |
(0) = 0, A(1) |
= |
1 |
|
, A(2) = 0. |
|
|
|
|
||
|
|
|
|
|
||||||||
1 |
2 |
3 |
6 |
3 |
|
|
|
|
|
|
||
Дифференцируем (3.67): |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||
|
|
y(0v) = xy00 |
+ 2y0. |
|
|
|
(3.68) |
|||||
При x = 0 |
|
y(0v)(0) = 2y0(0). |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|||||
Отсюда y1(0v)(0) = 0, y2(0v)(0) = 2, A4(1) = 0, A4(2) = |
|
2 |
= |
|
1 |
. |
||||||
4! |
|
|||||||||||
Дифференцируем (3.68): |
|
|
|
|
|
12 |
||||||
|
|
|
|
|
|
|
|
|
|
|||
|
|
y(v) = xy000 |
+ 3y00. |
|
|
|
|
|
|
|||
При x = 0
y(v)(0) = 3y00(0) = 0.
Отсюда A(1)5 = A(2)5 = 0 и функции y1(x), y2(x) представляются рядами
x3 y1(x) = 1 + 6 + ... ,
x4 y2(x) = x + 12 + ... ,
где многоточиями обозначены члены рядов, содержащие x6, x7, ... Общее решение уравнения (3.66), таким образом, имеет вид
y(x) = c1 + c2x + c1x3 + c2x4 + ... .
6 12
68
Список литературы
1.Бондарев А. С., Доценко М. Л., Фролова Е. В. Математический анализ (функции одной вещественной переменной): Учеб. пособие / Под ред.
А.Л. Белопольского; СПбГЭТУ. СПб., 1988. 100 с.
2.Колбина С. А., Пилюгин С. Ю. Линейная алгебра (дополнительные главы): Учеб. пособие. СПб.: Изд-во СПбГЭТУ “ЛЭТИ“, 2009. 60 с.
3.Белопольский А. Л., Бодунов Н. А., Червинская Н. М. Типовые расчеты по курсу “Алгебра и геометрия“: Учеб. пособие. СПб.: Изд-во СПбГЭТУ “ЛЭТИ“, 2008. 92 c.
69