Материал: Diff_lektsii__1

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

Теорема 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