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

Конспект лекций по математике для студентов инженерно-технических специальностей. В 4 ч. Ч. 3

.pdf
Скачиваний:
1
Добавлен:
24.11.2025
Размер:
2.02 Mб
Скачать

В данном случае оценка производных f ( n) (x) затруднительна, поэтому воспользуемся другим, более общим приемом. Функция f (x) = (1 + x)α

удовлетворяет равенству

(1 + x) f ′(x) = αf (x)

(16.39)

и условию f (0) =1.

Так как решение задачи Коши для дифференциального уравнения (16.39)

единственно, то если найдется степенной ряд 1 + an xn , удовлетворяющий

n=1

уравнению (16.39), то он будет искомым разложением заданной функции. Итак

(1 + x) (a1 + 2a2 x + 3a3 x2 + ... + nan xn1 + ...) = α(1 + a1 x + a2 x2 + ... + an xn + ...)

или

 

 

 

 

 

 

 

 

 

 

 

 

a + (a + 2a

)x + (2a

2

+ 3a

)x2

+ ... + (na

n

+ (n +1)a

n+1

)xn =

1

1

2

 

3

 

 

 

 

 

 

 

= α + αa x + αa

2

x2

+ ... + αa

xn + ...

 

 

 

 

 

1

 

 

 

n

 

 

 

 

Приравнивая коэффициенты при одинаковых степенях x , имеем

a1

= α

=

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1 1)

 

 

α(α 1)

 

 

 

 

 

 

 

 

 

 

 

 

a1

+ 2a2 = αa1 a2 =

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a2 +

3a3 = αa2 a3

 

 

a

2)

 

 

α(α 1)(α 2)

 

 

=

2

 

 

 

=

 

 

 

 

 

3

 

 

3!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an

=

α(α 1)(α 2)...(α n + 1)

, n N .

n!

 

 

 

193

Таким образом,

 

 

 

 

 

 

 

 

 

(1 + x)α =1 + αx +

 

α(α 1)

x2

+ ... +

α(α 1)(α 2)...(α n +1)

xn +

 

 

 

 

 

 

2!

 

 

 

 

n!

+ ...+ =1 + α(α

 

 

 

(16.40)

1)(α 2)...(α n +1) xn , x (1;1)

 

 

 

 

 

 

 

 

n=1

 

 

n!

 

 

 

поскольку радиус сходимости полученного ряда равен 1, что несложно проверить, использовав признак Д’Аламбера. Ряд, стоящий в правой части формулы (16.40), называется биномиальным.

При α = n N все коэффициенты ряда (16.40), начиная с номера n +1, обращаются в нуль, и степенной ряд преобразуется в бином Ньютона.

(1 + x)n =1 + nx +

n(n 1)

x2

+ ... + x4 = n

Cnk xk .

 

 

n!

k =0

 

16.15. Приложение степенных рядов

Приближенное вычисление значений функций

Для нахождения приближенного значения функций f (x) в точке x0 с

заданной точностью поступим следующим образом. Разложим функцию f (x) в

ряд по степеням x x1 с интервалом сходимости, содержащим точку x0 , где x1

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

 

 

 

 

an (x0

x1 )n

оставим только члены,

гарантирующие

заданную точность

n=0

 

 

 

 

вычислений.

Минимальное число n0

таких членов

ряда определим из

194

соответствующей оценки либо остатка Rn (x0 ) формулы Тейлора, либо остатка rn (x0 ) ряда Тейлора, так как в случае сходимости степенного ряда функции f (x)

они равны между собой.

Пример 1 6 . 1 0 . Вычислить с точностью ε = 0,01 число e .

Решение. Так как

 

 

ex = n

xk

 

+ Rn (x) =

n

xk

+

ec

 

 

xn+1 , 0 < c < x, x R ,

 

 

k!

 

 

(n +

1)!

 

 

 

 

k =0

 

 

 

 

 

 

 

 

 

 

k =0

k!

 

 

 

то из оценки

 

R (1)

 

=

 

 

ec

 

 

<

 

 

 

3

 

0,01 следует, что n 5, т. е. n

0

=5 . Итак,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

(n +1)!

 

 

 

(n

+1)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e 1 +1 +

 

1

 

+

 

1

+

1

 

 

+

1

 

2 + 0,500 + 0,167 + 0,042 + 0,008 = 2,717

 

 

 

 

 

 

 

4!

 

 

 

 

2!

 

3!

 

5!

 

 

 

 

 

 

 

 

 

 

или

e 2,72 (ε = 0,01) .

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

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

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

sin x

с точностью ε = 0,001.

x

0

 

Решение. Используя разложение (16.37) для sin x , получаем

sin x

x

2n

 

= (1)n

 

 

, x R

x

(2n +1)!

n=0

 

195

1

 

sin x

 

 

(1)

n

 

 

1

 

 

(1)

n

 

 

 

dx =

 

 

 

 

 

x2n dx =

 

 

 

 

 

 

 

 

(2n +1)!

(2n +1)!(2n +1)

0

 

 

x

 

 

 

 

n=0

0

 

 

n=0

 

 

 

 

rn

 

 

1

 

 

 

 

 

 

 

0,01 n > 2 n0

= 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2n + 3)!(2n + 3)

 

1 sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

dx 1

 

 

 

+

 

 

1

0,0555 + 0,0016 = 0,9461

 

 

x

 

 

 

 

3 3!

5 5!

или

 

 

 

 

 

 

 

 

0

 

sinx x

0,946

(ε = 0,001) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисление сумм числовых рядов

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

Пример 1 6 . 1 2 . Найти сумму сходящегося числового ряда

(1)

n

 

 

 

.

(16.41)

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

n=0

 

 

Решение. Рассмотрим степенной ряд

(1)

n

x

2 n+3

 

S(x) =

 

 

, S(0) = 0,

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

n=0

 

сходящийся (равномерно) при x [− 1;1]. Ряд

196

(1)

n

x

2n+1

 

S(x) = x

 

 

= xarctg x

2n +

1

n=0

 

сходится равномерно на [0;1], поскольку в этом случае он является рядом Лейбница и

r (x)

 

 

 

x

 

2n+2

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

2n

+1

 

2n +1

 

 

 

 

Тогда ряд S(x) можно дифференцировать почленно на [0;1]. Решая дифференциальное уравнение S(x) = xarctgx с начальным условием S(0) = 0,

находим

S(x) = x tarctgtdt =

1

((1 + x2 )arctg x x), x [0;1].

2

0

 

Так как сумма числового ряда (16.41) равна S(1) , то

(1)

n

 

1

 

π

 

1

 

 

=

(2arctg1 1) =

.

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

2

4

2

n=0

 

 

 

 

Интегрирование дифференциальных уравнений

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

Пример 1 6 . 1 3 . Найти решение

yy= sin y,

(16.42)

197

удовлетворяющее начальному условию y(0) = π .

 

2

 

Решение. Уравнение (16.42) допускает разделение переменных:

 

 

ydy

= dx,

(16.43)

 

 

 

sin y

 

однако интеграл от левой части уравнения (16.43) не выражается в элементарных функциях. Будем искать решение в виде ряда Маклорена

( n)

(0)x

n

y(x) =

y

 

 

.

 

 

n!

 

n=0

 

 

 

Так как y(0) = π2 , а

y

=

sin y

,

(16.44)

 

y

то y(0) = π2 . Дифференцируя по x обе части равенства (16.44), находим

 

 

y′′ =

( ycos y) y ysin y

=

y( y cos y sin y)

,

(16.45)

 

 

y2

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

y(0)sin

 

 

 

2

 

3

 

 

 

 

откуда y′′(0)

 

 

 

 

 

2

 

 

 

 

 

 

=

 

 

 

 

 

=

 

.

 

 

 

 

π

2

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

Дифференцируя обе части равенства (16.45), находим

y′′′(0) . Продолжая

этот процесс, можно получить любое число членов разложения в ряд Маклорена искомого решения y = y(x) :

198

y = π + 2 x 22 x2 + ...

2 π π3

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

РЯДЫ ФУРЬЕ

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

Ряды Фурье используются для описания периодических процессов, решения дифференциальных уравнений, приближения периодических и непериодических функций. В этих случаях, функцию, описывающую периодический процесс, представляют как сумму конечного и бесконечного числа периодических функций вида Asin(ωx + ϕ0 ) , где A, ω, ϕ0 постоянные. Эти функции описывают простейшее колебательное движение и называются гармониками. A > 0

называется амплитудой, ω − частотой колебания, ωx + ϕ0 фазой колебания,

ϕ0 начальной фазой.

Используя формулы тригонометрии, можно записать

Asin(ωx + ϕ0 ) = a cos ωx + bsin ωx,

где a = Asin ϕ0 , b = Acosϕ0 .

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

199

(an cosωn x + bn sin ωu x)

n

Во всех этих случаях имеем дело с периодическими функциями. Напомним свойства этих функций.

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

2.Если функция f(x) имеет период Т, то функция f(ax) (a ≠ 0) имеет период T / a .

3.Определенный интеграл от периодической функции f(x) с периодом Т по любому отрезку длиной Т имеет одно и то же значение, т. е. для любых х справедлива формула

x+T

T

 

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

(16.46)

x

0

 

Действительно, т. к. производная от интеграла

x+T

f (x + T )

f (x) = 0

 

f (t)dt

=

 

x

 

 

 

 

т. е. не зависит от х.

16.17. Ортогональные системы функций

Определение 16.8. Функция f(x) называется кусочно-непрерывной на отрезке [a,b], если она непрерывна на этом отрезке, за исключением, быть может,

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

200

Пусть f(x) – кусочно-непрерывная на [a,b] функция и x0 [a,b] точка разрыва первого рода. Тогда в этой точке существуют односторонние пределы f (x0 ± 0) , поэтому на каждом участке непрерывности существуют определенные

интегралы Римана b

f (x)dx и b

f 2 (x)dx . В этом случае функция f(x) называется

a

a

 

функцией с интегрируемым квадратом.

Так как на множестве кусочно-непрерывных функций определены линейные операции, удовлетворяющие аксиомам линейного пространства, то это множество образует линейное пространство. Введем в нем операцию скалярного произведения функций ϕ(x) и ψ(x) .

Определение 16.9. Скалярным произведением функций ϕ(x) и ψ(x) на отрезке [a,b] называется число

(ϕ,ψ) = b ϕ(x)ψ(x)dx

(16.47)

a

 

Скалярное произведение функций обладает следующими свойствами:

1)(ϕ,ψ) = (ψ,ϕ);

2)(ϕ1 + ϕ2 ,ψ) = (ϕ1 ,ψ) + (ϕ2 ,ψ);

3)(αϕ,ψ) = α(ϕ,ψ) α R;

4)(ϕ,ϕ) 0, (ϕ,ϕ) = 0 ϕ = 0, т. е. удовлетворяет аксиомам евклидова пространства.

Множество всех кусочно-непрерывных функций на [a,b] со скалярным произведением, введенным по формуле (16.47), обозначают L2 [a,b] и называют пространством L2 , которое бесконечномерное.

201

Определение 16.10. Неотрицательное число ϕ = b ϕ2 (x)dx называется

a

нормой функции ϕ(x) в L2 [a,b]. Учитывая, что b ϕ2 (x)dx = (ϕ,ϕ) , можно записать

a

 

ϕ

 

 

 

= (ϕ,ϕ) .

(16.48)

 

 

 

Функция ϕ(x) называется нормированной, если ее норма равна единице.

Определение 16.11. Две функции ϕ(x) L2 [a,b] и ψ(x) L2 [a,b] называются ортогональными на [a,b], если их скалярное произведение на [a,b] равно нулю,

b

т. е. (ϕ,ψ) = ϕ(x)ψ(x)dx = 0.

a

Определение 16.12. Система функций (ϕn (x)) = (ϕ1 (x),ϕ2 (x),...ϕn (x)...)

(конечная или бесконечная) называется ортогональной на отрезке [a,b], если все функции этой системы попарно ортогональны на [a,b], т. е.

(ϕm ,ϕn ) = 0 m n, m, n N .

Определение 16.13. Ортогональная система функций (ϕn (x)) на отрезке

b

[a,b] называется ортонормированной, если ϕ 2 = (ϕn ,ϕn ) = ϕ2n (x)dx =1n N .

a

Любую ортогональную на [a,b] систему функций (ϕn (x)) с ϕn 0 можно нормировать. Для этого достаточно разделить каждую функцию системы на ее норму. Получим ортогональную систему функций (ϕn (x) /ϕn ).

202

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]