Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

14-es_1 / Высшая математика (РТФ) / умк_Цывис_Диф.уравнения_Ряды

.pdf
Скачиваний:
35
Добавлен:
18.05.2015
Размер:
5.68 Mб
Скачать

 

 

Пример 3. Вычислить 5 35

 

 

 

с точностью до 0,0001.

 

 

 

 

 

Решение. В силу равенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

4 × 32

 

 

 

 

 

 

 

 

 

 

35

 

 

3

 

 

1

 

 

 

3

 

 

 

 

 

3

 

9

 

5

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

35 = 25

 

 

 

= 2 1 +

 

 

 

= 2 1 +

 

 

×

 

 

 

-

 

 

 

 

 

 

= 2

+

 

-

 

= 2,0361,

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

32

 

 

32

 

 

5

 

 

 

32

21× 5

× 32

2

 

 

80

 

6400

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

следующий член будет

 

4 × 9 × 33

 

 

 

£

 

1

= 0,00004 .

 

 

 

 

 

3!× 53 × 323

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25000

 

 

 

 

 

 

 

 

 

 

 

0,25

2. Приближенное вычисление интегралов

Определенные интегралы можно вычислить с помощью рядов.

Пример 4. Вычислить приближенное значение интеграла

ex2 dx , взяв 3 члена разложения в ряд подынтегральной функции, най-

0

ти погрешность вычисления.

Решение. Разложив подынтегральную функцию в степенной ряд, получим

 

 

 

 

ex2

=1 - x2 +

x4

 

-

x6

+ ... , тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2!

 

3!

 

 

 

 

 

 

 

 

0,25

ex

2

0,25

- x2 +

x

4

 

 

x

6

 

 

 

x

3

 

x

5

 

14

dx =

1

 

-

 

 

... dx = x -

 

+

 

- ...

=

 

 

 

 

 

 

 

 

 

 

 

 

 

2! 3!

 

 

3

 

2!× 5

0

0

 

 

0

 

 

 

 

 

=

1

-

1

+

 

1

 

- ...

 

 

 

 

 

 

 

 

4

 

43 × 3 45 × 2!× 5

 

 

 

Полученный ряд знакопеременный.

 

 

 

 

Абсолютная погрешность будет меньше, чем первый отброшенный

член, т.е. меньше, чем

1

 

< 0,00001. Поэтому, производя вычисления

 

 

47 × 3!× 7

с точностью до 0,00001, будем иметь

 

 

 

 

 

0,250000 − 0,005208 + 0,000098 = 0,244890 .

0,25

 

 

 

 

 

 

 

 

 

 

Следовательно,

ex2 dx = 0,24489

(с точностью до 0,00001).

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

sin x

 

Пример 5. Вычислить интеграл

Six =

dx .

 

 

 

 

 

 

 

 

 

 

0

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

191

 

 

 

 

 

 

Решение. Так как sin x = x -

x3

+

x5

 

- ... + (-1)n

 

x2n +1

 

+ ... , откуда

 

 

 

 

 

 

(2n + 1)!

 

 

 

 

 

 

 

 

 

 

3!

5!

 

 

 

 

sin x

=1 -

x2

+

x4

- ... + (-1)n

x2n

 

 

, тогда

 

 

 

 

 

 

x

 

 

 

(2n +

 

 

 

 

 

 

 

 

3!

5!

 

 

1)!

 

 

 

 

x2n +1

 

 

 

 

 

 

 

 

 

x sin x

x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx = x -

 

+ ... +

(-1)

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

x

3!× 3

 

 

 

 

 

(2n + 1)!(2n

+1)

 

 

 

Подставляя в полученный ряд значения x, получим величину интеграла.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x sin x

 

 

 

 

 

 

 

Пример 6. Вычислить

I (x) =

 

 

 

dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

x

 

 

 

 

 

 

 

Решение.

 

 

 

 

 

 

 

 

 

 

 

 

= y ,

 

 

 

 

 

 

 

 

 

Полагая

 

 

 

x

приводим

 

данный интервал к виду

y

sin y2

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I (x) =

 

dy2 = 2sin y2dy , откуда

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

y

6

 

 

 

 

 

 

 

 

y

4n + 2

 

 

 

 

 

 

 

 

I (x) = 2

y2

-

 

+ ... + (-1)n

 

 

 

dy .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3!

 

 

 

 

 

 

 

(2n +1)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ряд, стоящий под знаком интеграла сходится при всех y, поэтому

 

 

 

 

 

 

 

 

y3

 

 

 

y7

 

 

(-1)n

 

 

 

y4n + 2

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I (x) = 2

 

 

 

-

 

 

 

+ ... +

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

3!× 7

 

 

 

 

 

 

 

 

 

 

(2n + 1)!(4n + 3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

x

 

 

 

x

8

 

 

 

 

 

 

 

 

 

x

2n +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

2 x

-

 

 

+ ... + (-1)n

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

3!× 7

 

 

 

 

 

 

 

 

 

 

 

 

+ 3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2n + 1)!(4n

 

3. Нахождение сумм числовых рядов

Используя известные разложения в степенной ряд, свойства сходя- щихся рядов, операции интегрирования и дифференцирования находим значение определенной функции в заданной точке.

Пример 7. Найти сумму ряда:

1)

1 + 2x + 3x2 + ... + (n + 1)xn + ... ;

 

 

 

 

 

 

 

2)

1× 2 + 2 × 3x + ... + (n + 1)(n + 2)xn + ...

 

 

 

 

 

 

 

Решение. 1) В силу равенств

 

 

 

 

 

 

 

 

1 + x + x2 + ... + xn + ... =

 

1

 

,

 

x

 

<1,

 

 

 

 

 

 

 

 

 

 

 

1 - x

 

 

 

 

 

 

 

 

 

 

x + x2 + x3 + ... + xn +1 =

 

x

,

 

x

 

<1,

 

 

 

 

 

 

 

- x

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

192

 

 

(x + x2 + x3 + ... + xn +1 )¢

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

x

<1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 - x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + 2x + 3x2 + ... + (n + 1)xn =

 

 

1

 

 

,

 

 

 

x

 

<1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 - x)2

 

 

 

2) В силу равенств

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + x + x2 + ... + xn + ... =

 

 

 

1

 

 

,

 

x

 

<1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 - x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 + x3 + x4 + ... + xn + 2 =

 

x2

 

= -(1 + x) +

 

 

1

,

 

 

 

 

- x

 

 

- x

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

и, дважды дифференцируя последнее равенство, получим

 

 

 

 

 

1× 2 + 3 × 2x + ... + (n + 1)(n + 2)xn =

 

 

 

 

 

 

 

2

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 - x)3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(-1)n

 

Пример 8. Найти сумму числового ряда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

(2n +

1)(2n + 3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n = 0

 

Решение. Данный числовой знакочередующийся ряд сходится абсо-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

лютно (сравнить с числовым сходящимся рядом

 

 

 

 

 

 

 

 

). Рассмотрим сте-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2n + 3

 

 

 

 

 

 

 

 

 

 

 

 

n =1 n2

 

 

 

 

 

 

 

 

(-1)

n

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пенной ряд

S (x) =

 

 

 

; S (0) = 0 , который равномерно схо-

 

 

 

 

 

 

 

 

 

 

n = 0

(2n + 1)(2n + 3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дится при x Î[-1;1].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ряд S ¢(x) = x

(-1)n x2n +1

 

= xarctgx равномерно сходится при x Î[0;1] .

 

 

 

 

 

 

 

 

 

2n + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S (0) = 0 ,

 

Решая

задачу

Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получим

 

 

 

S (x) = aarctgx,

 

 

 

S (x) =

1 + x2

arctgx -

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. Так как сумма исходного числового ряда равна

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(-1)

n

 

 

 

 

 

 

 

 

 

 

 

 

1

= p -

1

 

 

 

 

 

 

 

 

 

 

 

S(1), то имеем

 

 

 

 

 

 

 

 

 

= arctg1 -

 

.

 

 

 

 

 

 

 

 

 

+ 1)(2n +

 

 

 

 

 

 

 

 

 

 

n =0 (2n

3)

 

 

 

 

 

 

 

 

2

 

4

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

n +1

 

 

 

 

 

 

 

 

 

 

Пример 9. Найти сумму ряда (-1)n +1

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n =1

 

 

 

 

 

 

 

 

n(n + 1)

 

 

 

 

Решение. Радиус сходимости данного ряда

 

 

R = 1. Если

S (x) его

сумма, то эта функция удовлетворяет следующему дифференциальному

193

′′

 

1

 

 

 

 

 

 

 

 

 

 

 

 

= 1 + x с начальными условиями S (0)

= 0 .

 

уравнению S (x)

= S (0)

Решая

задачу Коши, получим, что S (x) = (x + 1)ln(1 − x) − (1 + x) + 1.

 

 

 

 

 

 

 

1

 

1

 

1

 

 

(−1)n +1

 

Пример 10.

Найти сумму ряда 1 −

 

+

 

 

 

+ ... +

 

 

 

+ ...

 

5

9

13

 

4n − 3

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Данный знакочередующийся ряд является сходящимся. Рассмотрим геометрическую прогрессию

1 − x4 + x8 x12 + x16 − ... + (−1)n x4n + ...,

 

ее сумма равна S (x) =

 

 

1

 

 

при

 

x

 

< 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + x4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрируя прогрессию, получим

 

 

 

 

 

 

 

 

 

x5

 

 

x

9

 

 

 

 

(−1)n +1 x4n + 3

 

 

 

x

 

 

dx

 

x

 

 

+

 

 

+ ...

 

 

 

 

 

 

 

 

 

 

+ ...+ =

 

 

 

.

5

 

9

 

 

4n + 3

1

+ x4

 

 

 

 

 

 

 

 

 

 

0

 

Интегрируя и полагая в последнем равенстве

 

 

x = 1, находим сумму

данного числового ряда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S =

 

2

ln(3 +

 

 

 

2

π .

 

 

 

 

 

 

 

 

 

 

2) +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

4. Решение дифференциальных уравнений с помощью рядов

В теоретических и прикладных исследованиях в качестве математи- ческих моделей реальных процессов часто используются такие дифферен- циальные уравнения, решения которых на основании общей теории диф- ференциальных существуют, но не являются элементарными функциями, или решение мы не можем найти. В этих случаях решения ищут в виде степенных рядов и некоторых элементарных функций. При таком подходе можно получить аналитические представления решений для широкого класса весьма сложных линейных и нелинейных дифференциальных урав- нений, в том числе уравнений в частных производных. Процесс нахожде- ния решения дифференциального уравнений с помощью степенных рядов обычно осуществляется в три этапа:

1)определяется вид степенного ряда;

2)определяем коэффициенты этого ряда (вычисляют их некоторым способом);

3)исследуют характеристики полученного ряда (прежде всего, ради- ус сходимости, область сходимости).

194

Отметим, что внутри области сходимости ряда решение y(x) ис-

ходного дифференциального уравнения может быть вычислено с любой наперед заданной точностью (в виде частичной суммы ряда с соответст- вующим числом членов).

Может оказаться, что решение исходного дифференциального урав- нения существует на промежутке более широком, чем область сходимости. Поэтому частичную сумму ряда можно использовать для нахождения при-

ближенных значений y(x) решения y(x) и за областью сходимости. Од-

нако в этом случае к полученным значениям y(x) следует относиться с

осторожностью, потому что их погрешность может оказаться неприемле- мой, т.к. вне области сходимости нет общих способов оценки отбрасывае- мый остаток ряда.

Рассмотрим два основных метода нахождения решений дифферен- циальных уравнений в виде степенных рядов на примере обыкновенных дифференциальных уравнений второго порядка.

4.1. Метод последовательного дифференцирования

Данный метод применяется для решения линейных и нелинейных дифференциальных уравнений и базируется на следующей теореме.

Теорема 1. Если в задаче Коши для дифференциального уравнения второго порядка

 

 

y′′ = f (x, y, y′)

 

(1)

 

 

 

 

 

 

y(x0 ) = y0

; y (x0 ) = y0 .

 

 

функция

в некоторой окрестности точки

имеет

f (x0 , y0 y )

(x0 , y0 , y0 )

ограниченные производные по x, y и y′ , то существует единственное реше-

ние задачи (1), которое можно представить в виде ряда Тейлора

 

 

 

 

 

y

(n) (x

0

)(x x

0

)n

 

 

 

 

 

 

 

y(x) = an (x x0 )n =

 

 

 

 

 

 

,

 

 

 

(2)

 

 

 

 

 

 

n!

 

 

 

 

 

 

 

 

 

n = 0

 

 

n = 0

 

 

 

 

 

 

 

 

 

 

 

сходящегося в некотором интервале

 

x x0

 

< R .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первые два коэффициента ряда (2) определяют по начальным условиям

 

y(x0 )

 

 

 

 

 

 

 

 

 

y

(n)

(x0 )

 

 

y (x0 )

 

 

 

 

 

 

 

 

 

a0 =

 

= y0 ; a1 =

 

= y0 , а остальные коэффициенты

an =

 

 

 

 

0!

1!

 

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вычисляют в два этапа: сначала находят

выражения

для

производных

y(n) (x) (последовательным дифференцированием по

x

обеих частей ис-

195

ходного уравнения), а затем в полученные выражения подставляют x = x0

и найденные значения y(x0 ), y¢(x0 ), y¢¢(x0 )...y(n) (x0 ) в ряд (2). Если удает-

ся получить y(n) (x) в общем виде тогда находим точное решение y(x) ,

и найти область сходимости полученного ряда где решение задачи Коши (1) в виде ряда существует. Если же коэффициенты ряда (2) в общем виде найти не удается, то выписывают несколько первых членов ряда Тейлора

частичную сумму, которую принимают за приближенное решение y(x)

исходной задачи. В этом случае остаются неизвестными как область схо- димости ряда, так и величина погрешности приближенного решения, до- пускаемая при отбрасывании остатка ряда.

y¢¢ = xy¢ + y

Пример 11. Найти решение задачи Коши

y(0) = y¢(0) =1.

Решение. Так как первые и вторые частные и смешанные производ- ные по x, y и yот правой части уравнения равномерно ограничены в

окрестности точки (0, 1, 1), а все частные и смешанные производные более высоких порядков равны нулю, то по теореме 1 решение данной задачи в виде ряда Тейлора (2) существует в некоторой окрестности точки x0 = 0 .

Из начальных условий имеем: a0 = y(0) = 1;

 

 

 

a1 = y (0) = 1.

 

 

Для нахождения остальных коэффициентов последовательно диффе-

ренцируем по x обе части исходного уравнения.

 

 

y′′ = xy′ + y

 

 

 

 

 

 

y′′′ = xy′′ + y′ + y′ = xy′′ + 2 y

 

 

y( IV ) = xy¢¢¢ + 3y¢¢

 

 

 

 

…………………………

 

 

 

 

 

 

y(n) = xy(n −1) + (n -1) y(n − 2)

 

 

Подставляя в полученные выражения

 

= 1,

на-

x = 0, y(0) = 1, y (0)

ходим y¢¢(0) =1, y¢¢¢(0) = 2 ×1, ...,

y(n) (0) = (n -1) × y(n − 2) (0) .

 

 

При вычислении y(n) (0) рассмотрим два случая:

 

 

n = 2k

 

y( 2k) = ( 2k -1)!,

 

 

n = 2k + 1 y(2k +1) = (2k)!!.

 

 

Тогда значения искомых коэффициентов ряда

 

 

a2k =

y

(2k ) (0)

=

(2k

− 1)!

 

 

 

 

 

 

 

,

 

 

 

(2k )!

 

 

 

 

 

 

 

(2k )!

 

 

196

y(x)

a2k +1 =

y(2k +1) (0)

=

 

(2k )!!

 

=

 

1

 

 

.

(2k + 1)!

(2k + 1)!

(2k

+

 

 

 

 

1)!!

Таким образом, решение задачи Коши имеет вид

 

x2k

 

 

x2k +1

 

 

 

 

y = 1 + x +

 

 

 

 

 

+

 

 

 

 

 

.

 

(2k )!!

(2k + 1)!!

 

 

 

 

 

 

 

k =1

 

 

 

 

k =1

 

 

 

 

 

 

Радиусы сходимости полученных рядов

R = ∞ , следовательно, полу-

ченное в виде ряда решение задачи Коши сходится на интервале (−∞;+∞) .

4.2. Метод неопределенных коэффициентов

Рассмотрим данный метод для нахождения общего решения линей- ного дифференциального уравнения второго порядка

y′′ + p(x) y′ + q(x) y = f (x) ,

(3)

а также для нахождения частных решений (задача Коши) для уравнения (3) при начальных условиях

(4)

y(x0 ) = y0 , y (x0 ) = y0 .

При этом степенной ряд с неопределенными коэффициентами под- ставляет в уравнение (3), а значит, исходя из соотношений, определяемых полученным равенством, определяют коэффициенты ряда. Вид и свойства полученного решения в большей мере зависит от свойств p(x), q(x) и f (x) .

Теорема 2. Если функции

p(x), q(x) и

 

f (x)

разложимы в сте-

пенные ряды по степеням

(x x0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0 )n , (5)

p(x) = pn (x x0 )n ;

q(x) =

qn (x x0 )n ;

f (x) =

fn (x

n = 0

 

 

 

 

 

n = 0

 

 

 

 

n = 0

 

сходящиеся в области

 

x x0

 

< R , то существует единственное решение

 

 

задачи Коши (3) – (4),

представимое в виде ряда Тейлора

 

y(x) =

x0 )n , an =

y

(n)

(x0 )

 

 

 

an (x

 

 

,

 

(6)

 

 

 

 

 

 

 

n = 0

 

 

 

n!

 

 

причем ряд (6) сходится, по крайней мере,

 

на том же

интервале

(x0 R, x0 + R) , что и ряд (5).

 

 

 

 

 

 

 

 

197

Пример 12. Решить задачу Коши y′′ = xy

y(0) = 1; y′(0) = 0 .

Решение. Будем искать решение данной задачи Коши относительно x

в виде степенного ряда y = a0 + a1x + ... + an xn + ... при начальных услови-

ях a0 = 1, a1 = 0 .

Дифференцируя дважды этот ряд и подставляя в данное уравнение, приравняв коэффициенты при одинаковых степенях, получим

 

2a2 = 0; 6a3 = 1; 12a4 = 0;

20a5 = a2 ; ...

Замечаем, что в этом случае имеем

 

 

 

 

 

 

 

 

 

 

 

a2 = a5 = ... = a3n + 2 = ... = 0 ,

 

 

 

 

 

 

 

 

 

a4 = a7 = ... = a3n +1 = ... = 0 ,

 

 

 

a3

=

1

 

,

a6 =

1

 

, ...,

 

a2n +3 =

1

 

, ...

 

 

5 ×

6

(3n + 2)(3n +

3)

 

2 ×

3

 

 

 

 

 

 

 

 

 

Таким образом, коэффициенты в исходном ряде будут

a

=1,

a

=

1

,

a

=

1× 4

,

... , a

=

1× 4 × 7...(3n - 2)

,

 

 

 

0

 

 

3

3!

6

6!

 

3n

 

(3n)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а остальные коэффициенты будут равны 0. Следовательно, получим ряд

 

x3

1× 4x6

1× 4 × 7...(3n - 2)x3n

y =1 +

 

+

 

 

+ ... +

 

,

3!

6!

 

 

 

 

(3n)!

который сходится везде.

198

§9. РЯДЫ ФУРЬЕ

В классе функциональных рядов важное значение имеют ряды Фу- рье, которые будут рассмотрены в этом параграфе.

9.1. Периодические процессы и периодические функции

Во всех случаях функцию, описывающую периодический процесс, представляют как сумму конечного или бесконечного числа простых пе- риодических функций A × sin(wx + j0 ) , где A, w, j0 постоянные величи-

ны. Такие функции называются гармониками, так как они описывают простейшее колебательное движение, называемое гармоническим. Посто-

янная A > 0 называется амплитудой колебания, w частотой колеба-

ния, wx + j0

фазой колебания, j0 начальной фазой.

 

 

 

Функция

y = f (x)

называется периодической, если для нее суще-

ствует такое число

T ¹ 0 , что выполняются следующие условия:

x - T и

1) при любом

x из области определения функции числа

x + T также принадлежат области определения;

 

 

 

 

2) f (x) = f (x + T ) = f (x - T ) .

 

 

 

 

 

Таким образом, функция

 

 

f(x) периодическая, тогда и только тогда,

когда $T ¹ 0; "x Î D( f );

(x ± T ) Î D( f ); f (x ± T ) = f (x) .

 

 

 

Число T ¹ 0

называют периодом данной функции. Если

T

явля-

ется периодом функции f(x), то и n ×T

период f(x), n Î N .

 

 

Если существует наименьший положительный период функции, то

его называют основным периодом. Например, функция

y = sin x

имеет

период 2p, y = cos 2x

π ;

sin

x

– 10p .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

Отметим, что если f(x)

имеет период T, то функция

f (w × x) –

тоже

периодическая, но с периодом

 

T

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

Заметим, что функция

 

f (x) = c,

c = const

с областью определения

D( f ) = R периодическая, с периодом

T = R,

где R любое число. Но

наименьшего (основного) периода функция не имеет. Напомним свойства периодических функций:

1. Сумма, разность, произведение и частное периодических функций периода T есть периодическая функция периода T;

199

2.

Если функция

 

 

y = f (x) имеет период T, то функция y = f (ax)

(a ¹ 0) имеет период

 

 

T

 

 

;

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x +T

T

3.

При любом x

 

справедливо равенство I =

f (t)dt = f (t)dt .

 

 

 

 

 

 

 

x

0

Действительно, так как (интеграл с переменным верхним пределом)

dI = f (x + T ) d (x + T ) f (x) dx = f (x + T ) − f (x) = 0 , dx dx dx

то I (x) = C , т.е. не зависит от x.

9.2. Тригонометрический ряд Фурье

Тригонометрическим рядом называется функциональный ряд вида

 

 

 

a0

 

 

 

 

 

 

 

+

(a cos nx + b sin nx) ,

(1)

 

 

 

 

 

 

2

 

n =1

n

n

 

 

 

 

 

 

 

где a0

, an , bn

(n = 1, 2 ...) –

действительные числа, называемые коэф-

фициентами ряда.

 

 

 

 

Тригонометрической системой функций называется множество функций

 

1,

cos x, sin x,

cos 2x, sin 2x, ..., cos nx, sin nx, ...

(2)

Теорема 1. Основная тригонометрическая система функций (2) об-

ладает следующими свойствами:

 

 

1)

интеграл по отрезку

[-p;p]

от произведения любых двух раз-

личных функций равен нулю, т.е.

 

 

 

 

π

 

 

 

 

 

 

1× sin mxdx = 0

 

 

 

 

−π

 

 

 

 

 

 

π

 

 

 

 

 

 

1× cos mxdx = 0

 

 

 

 

−π

 

 

 

 

 

 

π

 

 

 

 

 

 

cos mx cos nxdx = 0

(m ¹ n)

(3)

 

 

−π

 

 

 

 

 

 

π

 

 

 

 

 

 

sin mx × sin nxdx = 0

(m ¹ n)

 

 

 

−π

 

 

 

 

 

 

π

 

 

 

 

 

 

sin mx × cos nxdx = 0

(m, n =1, 2, 3 ...)

 

−π

200