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

книги / Числовые и степенные ряды. Ряды Фурье

.pdf
Скачиваний:
2
Добавлен:
20.11.2023
Размер:
972.67 Кб
Скачать

Вычислим коэффициенты:

 

 

c0 = −2,

c1 = 1.

 

 

 

 

 

 

 

 

 

3 4с + с = 0, c = −

c1

 

 

= −

1

.

 

 

 

 

 

 

 

 

 

 

 

2

1

 

 

2

 

 

 

 

 

 

3 4

 

 

3 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 6с + с = 0, c = −

 

c2

=

 

1

.

 

 

 

 

 

 

 

 

 

3

2

 

3

 

 

 

 

5

6 3 4 5 6

 

 

 

 

 

 

 

 

 

 

 

 

 

7 8c + c = 0, c = −

c3

 

 

 

= −

 

 

1

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

3

 

4

 

7

8

 

 

 

3 4 5 6 7 8

 

 

 

 

 

 

 

 

 

 

 

 

9 10c5 + c4 = 0, c5

= −

 

c4

 

 

=

 

 

 

 

 

1

 

 

 

.

9

10

 

3 4

5 6 7 8 9 10

 

 

 

 

 

 

 

 

 

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

(2n + 2)(2n +1)cn+1 + cn = 0,

cn+1

= −

c

=

(1)n+2

.

n

 

 

 

 

 

 

 

 

 

(2n + 2)

(2n + 1)

 

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

 

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

Преобразуем выражения для найденных коэффициентов. Умножимчислительи знаменателькаждойдробина2, получим

 

 

 

c

 

= −

 

2

= −

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

3 4

4!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

=

 

 

2

 

=

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

2 3

4 5

6 6!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c = −

 

 

 

 

2

 

 

 

= −

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

2 3 4

5 6

7 8

 

 

 

 

8!

 

 

 

 

 

 

 

 

 

 

c =

 

 

 

 

 

 

2

 

 

 

 

 

 

 

=

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

2

3

4 5 6 7 8 9 10

10!

 

 

. . . . . . . . . . . . . . . . . . .

61

cn

cn+1 =

=

(1)n+1 2

=

(1)n+1 2

.

2 3 4 5 (2n 1)2n

(2n)!

(1)n+2 2

=

(1)n+2 2

 

 

.

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

(2n + 2)!

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

Подставим найденные коэффициенты в степенной ряд:

y = c0 + c1x + c2 x2 + … + cn xn + cn+1xn+1 + cn+ 2 xn+ 2 + … .

Получим

y = −2 +

2

x

2

 

x

2

2

x

3

2

x

4

 

2

x

5

 

 

(1)n+1 2

x

n

 

 

 

 

 

 

 

+

 

 

 

 

 

 

+

 

 

… +

 

 

+ … =

 

2!

4!

 

6!

 

8!

 

10!

 

(2n)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

= −21

1

x +

1

x2

1

x3 +

1

x4

1

x5 … + (1)

xn + … = −2 (1)

xn .

 

 

 

 

 

 

2!

4!

6!

8!

 

10!

 

 

(2n)!

 

 

n=0 (2n)!

 

 

 

 

 

 

 

Решение данного уравнения есть степенной ряд

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = −2 n=1 (2n)! xn .

 

 

 

 

 

 

 

 

Область сходимости степенного ряда x (−∞;).

 

 

 

 

 

Замечание. Решение уравнения 4xy′′ + 2 y′ + y = 0

можно най-

ти иначе, используя знак суммы ряда:

 

 

 

 

 

 

 

 

 

 

y = c0 + c1x + c2 x2 + … + cn xn + cn+1 xn+1 + cn+ 2 xn+ 2 + … = cn xn .

n=0

Дважды продифференцируем ряд:

y′ = c1 + 2 c2 x + ... + n cn xn1 + (n + 1)cn+1 xn + ... = (n + 1)cn+1 xn .

n=0

62

y′′ = 2 c2 + 3 2 c3 x + ... + n(n 1) cn xn2 +

+ (n + 1) ncn+1 xn1 + (n + 2)(n + 1) cn+ 2 xn

 

 

 

 

 

+ ... = (n + 2)(n + 1) cn+ 2 xn .

 

 

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = cn xn ,

 

 

 

 

 

 

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 y′ = 2(n + 1)cn+1 xn ,

 

 

 

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4xy′′ = 4x (n + 2)(n + 1)cn+ 2 xn = 4(n + 2)(n + 1)cn+ 2 xn+1 =

 

n=0

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 4(n + 1)n cn+1 xn .

 

 

 

 

 

 

 

n=0

 

 

 

 

 

 

 

Подставим в уравнение 4xy′′ + 2 y′ + y = 0

найденные суммы ря-

да:

 

 

 

 

 

 

 

 

 

 

(n + 1)n cn+1xn +

 

 

 

 

 

4

2(n + 1)cn+1xn + cn xn = 0

n=0

 

 

 

n=0

 

 

n=0

 

или

 

 

 

 

 

 

 

 

 

 

4(n + 1)n cn+1 xn + 2(n + 1)cn+1 xn + cn xn = 0.

 

n=0

 

 

 

 

 

 

 

 

 

 

4(n + 1)n c

+ 2(n + 1)c

+ c xn = 0.

 

 

n+1

 

n+1

n

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

= 0.

 

 

(4(n + 1)n +

2(n + 1))cn+1 + cn x

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

= 0.

 

((2n + 2) 2n + (2n + 2))cn+1 + cn x

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

(2n + 2)(2n + 1)c

+ c xn

= 0.

 

 

 

 

 

n+1

n

 

 

 

 

 

n=0

63

Получили коэффициент при xn :

(2n + 2)(2n +1)cn+1 + cn = 0.

Выразим коэффициент cn +1 :

cn+1

= −

 

cn

 

(2n + 2)(2n

+ 1)

 

 

или

 

 

 

.

cn = −

cn1

 

 

 

 

 

 

 

2n(2n 1)

 

Получили рекуррентную формулу для вычисления коэффициентов разложения. Из начальных условий следует, что

c = y = −2, c = y′ = 1.

0

0

1

0

При n =2,3,… найдём соответствующие коэффициенты.

64

ГЛАВА 3. РЯДЫ ФУРЬЕ

3.1. Основные понятия

Определение 8. Функциональный ряд вида

a0 + a cosx + b sinx + a cos2x + b sin2x + a cos3x + b sin3x + …

2

1

1

2

2

3

3

 

…+ an cos(nx) + bn sin(nx) +…

 

 

 

(35)

называется тригонометрическим рядом, где постоянные числа a0 ,

an и bn (n = 1,2,3) – коэффициенты ряда.

 

Функции

cos(nx)

и

sin(nx)

(n =1,2,3)

– периодические,

с общим периодом 2π .

Если ряд (35) сходится, то его

сумма является периодической

функцией f (x)

с периодом 2π ,

т.е.

f (x) = f (x + 2π). Пусть на ин-

тервале (−π;π)

периодическую

функцию f (x) с периодом 2π

можно представить тригонометрическим рядом, сходящемся к данной функции:

 

a0

 

 

 

 

f (x) =

+ (an cos(nx) + bn sin (nx)).

(36)

 

2

n=1

 

 

 

 

Функции cos(nx)

и sin(nx) – ограниченные

 

cos(nx)

 

≤ 1 и

 

 

sin (nx) ≤ 1 . Рассмотрим числовой ряд, составленный их коэффициентов тригонометрического ряда:

a0 + a + b + a

+ b + a + b + …a

+ b + … .

(37)

2

1 1 2

2 3 3

n

n

 

 

 

 

 

 

Числовой ряд сходится абсолютно, если сходится ряд из абсо-

лютных величин:

 

 

 

 

 

a0

 

+

 

a

 

+

 

b

 

+

 

a

 

+

 

b

 

+

 

a

 

+

 

b

 

+ …

 

a

 

+

 

b

 

n

+ … .

(38)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

1

 

 

 

1

 

 

 

2

 

 

 

2

 

 

 

3

 

 

 

3

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

65

Следовательно, ряд (35) мажорируемый, его можно почленно интегрировать на отрезке [π;π] .

Если 2π – периодическую функцию f (x) можно представить в видетригонометрическогоряда(36), товозникаютследующиевопросы:

1)Каким условиям должна удовлетворять функция, чтобы её можно было разложить в тригонометрический ряд?

2)Как найти коэффициенты разложения в тригонометрический ряд?

Вначале ответим на второй вопрос.

3.2. Нахождение коэффициентов тригонометрического ряда

Пусть 2π – периодическая функция f (x) задана на отрезке [π;π] , тогда коэффициенты разложения в тригонометрический ряд можно вычислить по формулам

 

 

 

 

1

π

 

 

 

a0

=

f ( x)dx,

(39)

 

 

 

 

 

 

 

π − π

 

 

 

1

π

 

( x) cos(nx)dx,

 

an

=

f

(40)

 

 

π − π

 

 

 

 

 

 

1

π

 

(x) sin (nx)dx,

 

bn

=

f

(41)

 

 

π

− π

 

 

 

 

Коэффициенты a0 ,

an

и bn (n = 1,2,3) ,

найденные по фор-

мулам (39)–(41), называются коэффициентами Фурье для функции f (x), тригонометрический ряд (35) – рядом Фурье.

3.3. Достаточное условие разложимости функции в ряд Фурье

Пусть 2π – периодическая функция f (x), задана на отрезке

[π;π] .

Теорема 14 (Теорема Дирихле). Пусть 2π – периодическая функция f (x), на отрезке [π;π] удовлетворяет двум условиям:

66

1)функция f (x) кусочно-непрерывна, т.е. непрерывна или имеет конечное число точек разрыва 1-го рода;

2)функция f (x) кусочно-монотонна, т.е. монотонна на всём

отрезке, либо этот отрезок можно разбить на конечное число интервалов так, что на каждом из них функция монотонна.

Тогда соответствующий функции f (x) ряд Фурье сходится на этом отрезке и при этом

1) в точках непрерывности функции сумма ряда S (x) совпадает с самой функцией: S (x) = f (x) ;

2) в каждой точке x0 разрыва функции сумма ряда равна среднему арифметическому односторонних пределов:

S (x0 ) =

 

f (x0 0) + f ( x0

+ 0)

;

 

2

 

 

 

 

 

 

3) на концах отрезка в точках x = −π и x = π сумма ряда

S (−π) = S (

π) =

f (−π + 0) + f (π − 0)

.

 

 

2

 

 

 

Если 2π – периодическая функция

f (x) на отрезке [π;π]

удовлетворяет условиям 1 и 2 теоремы Дирихле, то её можно разложить в ряд Фурье

f ( x) = a0

+ (an cos(nx) + bn sin (nx))

 

 

2

n=1

с коэффициентами, определяемыми по формулам (39)–(41). Поскольку исходная функция f (x) и сумма ряда Фурье являются пе-

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

Замечание

1. Теорема Дирихле даёт достаточное условие разложимости функции, но не необходимое. Существуют функции, не удовлетво-

67

ряющие условиям теоремы Дирихле, но разложимы в тригонометрический ряд Фурье.

2. Функция f (x) имеет период 2π : f (x) = f (x + 2π) = = f (x + 2πk).

Коэффициенты ряда Фурье можно вычислить по любому промежутку длины периода функции, если на нём выполнены условия теоремы Дирихле:

a0 = π1 a+ 2π f (x)dx,

a

an = π1 a+ 2π f ( x) cos(nx)dx,

a

bn = π1 a+ 2π f ( x) sin (nx)dx.

a

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

sin(nπ) = 0 , cos(nπ) = (1)n .

Пример 41. Разложить в ряд Фурье функцию f (x) периода 2π , заданную на отрезке [π;π]:

x 1,

−π ≤ x < 0,

f ( x) =

0 x ≤ π.

x + 1,

На рис. 1 изображен график функции.

Функция f (x) удовлетворяет условиям теоремы Дирихле. Найдём коэффициенты разложения:

 

 

1

π

1

0

 

1

π

a0

=

f (x)dx =

 

(x 1)dx +

( x + 1)dx =

 

 

π

− π

π

− π

 

π

0

68

 

 

 

 

 

 

 

1

 

x2

 

 

 

0

 

x2

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

+ x

 

 

−π

+

 

 

+ x

 

0

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

1

 

 

 

(−π)2

 

+ (−π)

 

+

 

π2

 

 

 

 

=

 

 

 

 

 

 

 

0

 

 

 

 

 

 

+ π

 

 

 

 

 

 

 

π

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

1

 

π2

 

 

 

 

π2

 

 

 

=

1

2π = 2.

 

 

 

 

 

2

− π +

2

 

+ π

π

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

π

(x) cos(nx)dx =

 

 

 

 

 

 

 

 

 

 

 

an =

f

 

 

 

 

 

 

 

 

 

 

 

 

 

π

− π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0

(x 1) cos

(nx)dx +

1

π

 

 

 

 

 

 

 

 

 

 

 

 

=

 

(x + 1) cos(nx)dx =

 

 

 

π

− π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0

 

( x + 1) cos(nx)dx + 1

π

 

 

 

 

 

 

 

 

 

 

 

= −

 

( x + 1) cos(nx)dx.

 

 

 

 

π

− π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2y

 

 

 

 

 

 

 

 

 

 

 

 

y = f(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

-–π

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

–3π

–2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

2π

3π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

–1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 1

Вычислим неопределённый интеграл по частям:

(x + 1) cos(nx)dx =

 

u = x + 1,

dv = cos(nx)dx,

 

=

 

 

 

du = dx,

v =

1 sin (nx).

 

 

 

 

 

n

 

 

69

=(x + 1) 1n sin (nx) 1n sin (nx)dx =

=(x + 1) 1n sin (nx) 1n sin (nx)dx =

=(x + 1) 1n sin (nx) + 1n 1n cos(nx) =

=1n (x + 1)sin (nx) + n12 cos(nx) + C.

Нашли неопределённый интеграл:

( x + 1) cos (nx)dx = 1n ( x + 1)sin (nx) + n12 cos (nx) + C.

Тогда

 

 

 

= − 1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an

 

(x + 1)

cos(nx)dx +

 

 

(x + 1) cos(nx)dx =

 

 

 

 

 

 

π − π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

 

 

 

 

(x + 1)sin (nx) +

 

 

 

 

 

 

cos(nx)

 

−π

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

2

 

 

 

 

 

 

 

 

 

 

 

π n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

( x

+ 1)sin (nx) +

 

 

 

 

 

 

 

cos(nx)

 

 

 

 

=

 

 

 

 

 

 

 

 

 

π

 

 

 

 

n

2

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

= −

 

 

sin 0

+

 

 

 

 

cos 0

 

 

 

 

(−π + 1)sin (nπ) +

 

 

 

 

 

 

 

cos(nπ)

+

 

n

2

 

 

 

 

 

n

2

 

 

π n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

1

1

(π + 1)sin (nπ) +

1

cos(nπ)

1

sin 0 +

 

 

1

 

cos0

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

2

 

n

2

 

 

 

 

π n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

1

1

 

 

1

 

(

1)

n

+

1

1

 

(1)

n

 

1

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

2

 

n

2

 

 

 

 

 

 

 

2

 

 

n

2

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

1

 

 

1

+

 

1

(1)

n

+

 

1

 

(1)

n

1

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

n

2

n

2

 

n

2

 

 

n

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

70