Q Положив в формуле (64.7) о:= -~ и замеrшв хна (-х2), получим
равенство
1 |
х2 |
|
1 · 3 4 |
|
1 · 3 · 5 |
6 |
+"., |
хЕ[-1;1]. |
~=l+2+2.4x +2·4·6х |
Тогда |
|
|
|
|
|
|
|
|
|
ж |
1 |
|
ж |
ж t2 |
ж 1 . 3 |
dt + " " |
! |
~dt |
=! dt + |
! 2 dt + |
! |
2 . 4 t |
|
|
|
|
|
|
|
|
4 |
|
о |
1-t |
о |
|
о |
|
о |
|
|
|
или |
arcsiп х = |
1 |
х3 |
1 · 3 |
х |
5 |
|
|
|
|
х + - ·- |
+ -- ·- |
+ ... |
|
|
|
2 |
3 |
2·4 |
5 |
|
Можно показать, что полученное равенство справедливо при всех
Ряды (64.4)-(64.14) в комбинации с правилами сложения, вычита
ния, умножения, дифференцирования, интегрирования степенных ря
дов (см. свойства степенных рядов) могут быть использованы при ра~
ложении (некоторых) других функций в ряд Маклорена (Тейлора).
|
Прu.мер 64.1. Разложить в ряд Маклорена функцию f(x) = зж. |
О Решение: Так как 3ж = e1n з~ |
= ежln 3 , то, заменяя х на х lп 3 в раз |
ложении (64.4), получим: |
|
|
|
|
ж |
= |
1 ln 3 |
ln2 3 2 |
ln3 3 3 |
lnn 3 |
п |
+"" |
|
3 |
+11х+Т!х |
+-згх |
+".+ 7х |
|
хЕ(-оо;оо). 8 |
|
Прu.м.ер 64.2. Выписать ряд Маклорена функции f(x) = ln(4-x), |
О Решение: Так как |
|
|
|
|
|
|
|
f(x) |
= ln(4- х) = ln4(1- ~) = ln4 + ln[l + (-~)], |
то, воспользовавшись формулой (64.9), в которой заменим хна (-i),
получим:
или |
|
|
|
1 |
1 х2 |
1 |
xn+l |
ln(4- х) = ln4- 4х - |
42. 2 |
- ... - 4n+l |
. n + 1 - ... ' |
если -1 < -~ ~ 1, т. е. -4 ~ х < 4. |
• |
|
Пример 64.з. Разложить в ряд Маклорена функцию
|
|
|
f(x) = - |
2 |
|
|
|
|
|
- . |
|
|
|
|
|
|
3 -х |
|
|
а РРшение: Воспользуемся формулой (64.8). Так как |
|
f(x) |
= _2_ = |
2 |
= ~. _1_, |
|
|
|
3-х |
3·(1-j) 3 |
1-j |
|
то, заменив хна ~ в формуле (64.8), получим: |
|
|
|
_2 = ~ . (1 + ::. + (=)2 + (::.) 3 + |
.) |
|
3-х |
3 |
3 |
3 |
3 |
·· |
' |
или |
|
|
|
|
|
|
|
2 |
2 2 т |
|
2 х2 |
2 х2 |
2 х3 |
|
2 xn |
3 - Х = 3 + 3 . 3 + 3 . 32 + 3 . 32 + 3 . 33 + ···+ 3 . 3n + ···' |
где -1 < ~ |
< 1, 1. е. |
-3 < х < 3. |
|
|
• |
|
|
|
§ 65. НЕКОТОРЫЕ ПРИЛОЖЕНИЯ СТЕПЕННЫХ РЯДОВ
65.1. Приближенное вычисление значений функции
Пусть требуется вычислить значение функции /(х) при х = Х1 с
заданной точностью е > О.
Если функцию /(х) в интервале (-R;R) можно разложить в сте
пенной ряд
f(x) = ао + aix + а2х2 + ... + anxn + ...
и х1 Е (-R; R), то точное значение /(х1) равно сумме этого ряда при
Х = Х1, Т. е.
f (х1) = ао + aix1 + a2xi + ... + anxf + ... ,
а приближенное - частичной сумме Sn(x1), т.е.
f(x1) ~ Sn(x1) = ао + aix1 + a2XI + ... + апх?.
Точность этого равенства увеличивается с ростом п. Абсолютная по
грешность этого приближенного равенства равна модулю остатка ряда,
т.е.
где
Таким образом, ошибку l/(x1 ) - Sn(x1)\ можно найти, оценив остаток тп(х1) ряда.
Для рядов лейбницевского типа
lrn(x1)I = lиnн(х1) + Un+2(x1) + Un+з(x1) +. · -1 < lиnн(x1)I
(см. п. 61.1).
В остальных случаях (ряд знакопеременный или знакоположи тельный) составляют ряд из модулей членов ряда и для него старают
ся найти (подобрать) положительный ряд с большими членами (обыч но это сходящийся ряд геометрической прогрессии), который легко бы
суммировался. И в качестве оценки \тn(х1)1 берут величину остатка
этого нового ряда.
Пример 65.1. Найти sin 1 с точностью до 0,001.
Q Решение: Согласно формуле (64.5),
sin 1 = 1 - _.!:.. 13 |
+ _.!:.. 15 - |
00 |
1 |
|
. · . = ~(-1) n+l |
. |
3! |
5! |
L-J |
(2п - |
1)! |
n=l
Стоящий справа ряд сходится абсолютно (проверить самостоятельно).
Так как J! ~ 0,008(3) > 0,001, а j, ~ 0,0002 < 0,001, то для нахождения
sin 1 с точностью до 0,001 достаточно первых трех слагаемых:
sin 1 ~ 1 - 11 |
+ 11 |
= 0,842. |
з. |
5. |
|
Допускаемая при этом ошибка меньше, чем первый отброшенный член
(т. е. меньше, чем 0,0002). Вычисленное микрокалькулятором значение
sin 1 примерно равно 0,84147. |
8 |
При.мер 65.2. Вычислить число е с точностью до 0,001.
Q Решение: Подставляя х = 1 в формулу (64.4), получим:
1 |
1 |
1 |
+ ... |
е = 1 + 1 |
+ 21 |
+ ... + 1 |
1. |
. |
п. |
|
Справа стоит знакоположительный ряд. Возьмем п слагаемых и оце
ним ошибку тп(х):
1 |
1 |
1 |
тп(х) = (п + 1)! |
+ (п + 2)! |
+ (п + 3)! + ... = |
=(п~1)! ( 1 +п:2+ (п+ 2)1(п+ 3) + ···)<
<(п~1)! ( 1+п:1+(п~1)2+···)= (п~1)! ( l - 1ntI) = п!~п'
т. е. тп(х) < - 1- . Остается подобрать наименьшее натуральное число |
1 |
|
n.·n |
< 0,001. |
п, чтобы выполнялось неравенство - - |
1 |
|
п.1 ·n |
|
Нетрудно вычислить, что это неравенство выполняется при п ~ 6. |
Поэтому имеем: |
|
1 |
1 |
1 |
1 |
1 |
517 |
~ 2,718. |
• |
е ~ 1+1 + 2! + |
З! + |
4! + 5! + 6! = 2720 |
|
Заме'Чание. Оценку остатка ряда можно производить с помощью
остаточного члена ряда Маклорена
с
где с находится между О и х1. В последнем примере R..(1) = (n ~ l)!'
О < с< 1. Так как ее < е1 < З, то R..(1) < (n ! l)!. При п = 6 имеем:
R 6 (1) < fi < 0,001, е ~ 1+1 + J"! + ... + ~ ~ 2,718.
65.2. Приближенное вычисление определенных
интеrралов
Бесконечные ряды применяются также для приближенного вычи
сления неопределенных и определенных интегралов в случаях, когда
первообразная не выражается в конечном виде через элементарные
функции (см. § 34) либо нахождение первообразной сложно.
ь
Пусть требуется вычислить Jf(x) dx с точностью до с:> О. Если
а
подынтегральную функцию /(х) можно разложить в ряд по степеням х и интервал сходимости (- R; R) включит в себя отрезок [а; Ь], то для
вычисления заданного интеграла можно воспользоваться свойством по
членного интегрирования этого ряда. Ошибку вычислений определяют так же, как и при вычислении значений функций.
|
1 |
|
4 |
Пример 65.3. |
Вычислить интегра.п Jе-х• dx с точностью до |
с:= 0,001. |
о |
Q Решение: Разложим подынтегральную функцию в ряд Маклорена,
заменяя хна (-х2 ) в формуле (64.4):
х2 |
х4 хв |
е-х• = 1 - 1f |
+ 2! - З! |
+ |
... ' |
х Е (-оо; оо). |
(65.1) |
|
|
|
Интегрируя обе части равенства (65.1) на отрезке [о; :i], лежащем вну
три интервала сходимости (-оо; оо), получим:
1 |
1 |
|
|
|
|
4 |
4 |
+ 2хТ4 |
|
хЗ!6 |
|
1е-х2 dx = / |
( 1 - 1хТ2 |
- |
+ ..'~ ) dx = |
оо
Получили ряд лейбницевского типа. Так как |
1 1 з = 0,0052 ... > |
1 |
1.. 3. 4 |
> 0,001, а 2! . 5 . 45 < 0,001, то с точностью до 0,001 имеем:
l |
|
|
|
|
|
|
4 |
|
|
1 |
|
1 |
|
f е-х |
2 |
dx ~ |
- |
= 0,245. |
|
4 |
192 |
о
Заме..,,ание. Первообразную F(x) для функции J(x) = е-х2
найти в виде степенного ряда, проинте1рировав равенство (65.1)
делах от О дох:
х 2 |
х ( |
t2 |
t4 |
) |
F(x) = Je-t dt = J 1 - |
l! + |
2! - . . . |
dt = |
оо
|
|
|
|
хз |
х5 |
х1 |
|
|
|
|
|
= х - 1! · 3 + 2! · 5 - 3! · 7 + · · ·' х Е (-оо; оо). |
Функции |
|
|
х2 |
и F(x) = |
:z: |
' |
|
|
|
важ- |
J(x) = Jn:e -т |
Jf(t) dt играют очень |
|
|
|
|
|
|
о |
|
|
|
|
|
ную роль в теории вероятностей. Первая - |
плотность стандартно- |
го распределенш~ вероятностей, вторая - |
фу·н:кч1J.Я Лапласа F(x) = |
х |
t2 |
|
|
|
|
|
|
|
|
|
|
== А;: Jе- 2 dt (или |
|
интеграл |
веро.ятносте11.). Мы |
получили, |
что |
о |
|
|
|
|
|
|
|
|
|
|
|
функция Лапласа представляется рядом |
|
|
|
|
|
|
|
1 |
( |
х - |
х3 |
х5 |
|
х1 |
|
) |
' |
|
|
F(x) = ../27Г |
|
2 · 3 + 22 2! · 5 - |
23 • 3! · 7 + ·· |
|
|
который сходится на всей числовой оси.
65.3.Приближенное решение дифференциальных уравнений
Если решение дифференциального уравнения не выражается че
рез элементарные функции в конечном виде или способ его решения слишком сложен, то для приближенного решения уравнения можно воспользоваться рядом Тейлора.
Познакомимся с двумя способами решения дифференциальных
уравнений с помощью степенных рядов.
Пусть, например, требуется решить уравнение |
|
у"= f(x;y;y'), |
(65.2) |
удовлетворяющРе начальным условиям |
|
Ylx=.1:0 = Уо, y'Jx=xo =У~· |
(65.3) |
Способ последовательного дифференцирования
Решение у= у(х) уравнения (65.2) ищем в виде ряда Тейлора:
у = у(хо) + -у'(-хо-)(х - |
Хо) + -у"(-хо-)( |
х - |
хо)2 |
+ ... |
|
|
11 |
21 |
y(n)(xo) |
n |
|
|
|
+ "., (65.4) |
|
"·+ |
1 |
(х - |
хо) |
|
|
|
п. |
|
|
|
при этом первые два коэффициента находим из начальных усло
вий (65.3). Подставив в уравнение (65.2) значения х = хо, у = Уо,
у' = уЬ, находим трРтий коэффициент: у"(хо) = /(хо; у0; уЬ). Значения
у"'(х0), у<4)(х0), ... находим путем последовательного дифференциро
вания уравнения (65.2) по х и вычисления производных при х = хо. Найденные значения производных (коэффициентов) подставляем в ра венство (65.4). Ряд (65.4) представляет искомое частное решение урав нения (65.2) для тех значений х, при которых он сходится. Частичная
сумма этого ряда будет приближенным решением дифференциального
уравнения (65.2).
Рассмотренный способ применим и для построения общего реше
ния уравнения (65.2), если у0 и у~ рассматривать как произвольные
ПОСТОЯННЫ{'.
Способ последовательного дифференцирования применим для ре
шения дифференциальных уравнений любого порядка.
Пример 65.4. Методом последовательного дифференцирования
найти пять первых членов (отличных от нуля) разложения в ряд ре-
шения уравнения у"= х2 + у2, у(-1) = 2, у'(-1) = ~·
Q Решение: Будем искать решение уравнения в виде
у= у(-1) + у'(~1) (х + 1) +у"(2~1) (х + 1)2 + у///(3~1) (х + 1)3 + ...
Здесь у(-1) = 2,у'(-1) = !·Находим у'~!-1), подставив х = -1 вис
ходное уравнение: у"(-1) = (-1)2 + 22 = 5. Для нахождения последу
ющих коэффициентов дифференцируем заданное дифференциальное
уравнt>ние:
ym = 2х + 2уу',
у<4> = 2 + 2(у')2 + 2уу", у(5 ) = 4у'у" + 2у'у" + 2уу/// = 6у'у" + 2уу111 ' •••
При х = -1 имеt>м:
у |
"'( |
-1 |
) |
1 |
|
|
=-2+2·2·2=0, |
1
у(4>(-1)= 2 + 2. 4 + 2. 2. 5 = 22,5,
у(5)(-1) = 6. ~. 5 + 2. 2. о= 15" ..
2
Подставляя найденные значения производных в искомый ряд, полу
чим:
Способ неопределенных коэффициентов
Этот способ приближенного решения наиболее удобен для инте грирования линейных дифференциальных уравнений с переменными коэффициентами.
Пусть, например, требуется решить уравнещ1е
у"+ Pi(x)y' + Р2(х)у = f(x) |
(65.5) |
с начальными условиями у(хо) = Уо, у'(хо) = у0.
Предполагая, что коэффициенты р1 (х), р2(х) и свободный член f(x) разлагаются в ряды по степеням х - х0, сходящиеся в некото
ром интервале (хо - R; х0 + R), искомое решение у= у(х) ищем в виде
степенного ряда
у =Со+ С1 (х - хо) + С2 (х - хо)2 + ... |
+ cn(x - xo)n + ... |
(65.6) |
с неопределенными коэффициентами.
Коэффициенты Со я с1 определяются при помощи начальных усло
вий Со= уо, С1 = Уо·
Для нахождения последующих коэффициентов дифференцируем
ряд (65.6) два раза (каков порядок уравнения) и подставляем выраже ния для функции у и ее производных в уравнение (65.5), заменив в нем р1 (х), Р2(х), /(х) их разложениями. В результате получаем тождество, из которого методом неопределенных коэффициентов находим недо стающие коэффициенты. Построенный ряд (65.6) сходится в том же интервале (хо - R; хо+ R) и служит решением уравнения (65.5).
При.мер 65.5. Найти решение уравнения
у"+ ху' +у= xcosx, у(О) =О, у'(О) = 1,
используя метод неопределенных коэффициентов.
О Решение: Разложим коэффициенты уравнения в степенные ряды:
Р1(х) = х, Р2 = 1,
f(x)=xcosx=x(1-~~ +~~ -".).
Ищем решение уравнения в виде ряда
У= Со+ С1Х + С2Х2 + С3Х3 + ...
Тогда
у1 = с1 + 2с2х + 3сзх2 + 4с4х3 + ... ,
У11 = 2с2 + 2 · 3 · С3Х + 3 · 4 · С4Х2 + .. .
Из начальных условий находим: Со= О, с1 = 1. Подставляем получен
ные ряды в дифференциальное уравнение:
(2с2 + 2 · 3 · сзх + 3 · 4 · с4х2 + ... ) + х(с1 + 2с2х + Зсзх2 + 4с4х3 + ... )+
+(ео+с1х+с2х2 +сзх3 +".)=х(1-~~ +~~ -:~ +".).
Приравниваем коэффициенты при одинаковых степенях х:
х0 : |
2с2 =О, |
|
|
|
х1 : |
2 · 3 · С3 + 2 = 1, |
|
|
х2: 3 · 4 · С4 + 2с2 + С2 = 0, |
1 |
, |
хэ: |
4 · 5 ·cs + Зсз + сз = - |
х4: |
5 · 6 · С5 + 4с4 + С4 = 0, |
2 |
|
|
Отсюда находим, что с2 |
= |
с4 |
= сб = |
. . . = |
О, |
сз = - j!, cs = J!, |
с7 = -;, ,... Таким образом, получаем решение уравнения в виде |
у = х - |
хз |
xs |
х1 |
|
|
3! |
+ 5Т - |
7! + ... , |
|
Глава XV. РЯДЫ ФУРЬЕ. ИНТЕГРАЛ ФУРЬЕ
1 Лекции 56-581
§бб. РЯДЫ ФУРЬЕ бб.1. Периодические функции.
Периодические процессы
При изучении разнообразных периодu'Ческих прочессов, т е про
цессов, которые через определf'нный промежуток времени повторяются
(встречаются в радиотехнике, электронике, теории упругости, теории и практике автоматического регулирования и т. д.), целесообразнее раз
лага1ь периодические функции, описывающие эти процессы, Hf'в сте
пенной ряд, а в так называемый тригонометрический ряд.
Напомним, что функция у = f(x), определенная на множестве
D, называется периодu'Ческоfi, (см. п. 14.3) с периодом Т > О, если
при каждом х Е D значение (х + Т) Е D и выполняе1ся равенство
f(x + Т) = f(x).
Для построения графика периодической функции периода Т досга
точно построить его на любом отрезке длины Т и периодически про должить его во всю область определения.
Отметим основные свойства периодической функции.
1.Алгебраическая сумма периодических функций, имеющих один
итот же период Т, есть периодическая функция с периодом Т.
2.Если функция f(x) имеет период Т, то функция f(ax) имеет
период f: действительно, f (а· ( х + f)) =!(ах+ Т) = f(ax)
3. |
Если функция f(x) имеет период Т и интегрируема на отрезке |
|
а+Т |
Ь+Т |
[хо;х1] |
Е IR, то J J(x) dx = |
J f(x) dx при любых а и Ь Е [хо; х1]. |
|
а |
Ь |
Q Пусть, например, О < а < Ь < Т, тогда
а+т |
ь |
а+т |
|
J f(x) dx = |
Jf(x) dx + |
J f(x) dx. |
(66.1) |
аа
Сдругой стороны,
4Т |
~Т |
4Т |
|
J f(x).dx = |
J f(x) dx + |
J f(x) dx. |
(66.2) |
Ь |
Ь |
а+Т |
|
ь |
ь |
f (х) dx = (подстановка х = и+Т) = Jf (и+Т) dи = |
Jf (х)dx. |
а |
а |
Подставляем полученный результат в (66.2) и, сравнивая с (66.1), име-
а+т |
ь+т |
|
|
ем J f(x) dx = |
j |
f(x) dx. |
|
8 |
а |
Ь |
|
|
|
|
т |
|
ь+т |
В частности, |
j f(x) dx = |
j |
f(x) dx. |
оь
Простейшими периодическими функциями являются тригономе-
трические функции sin х и cos х. Период этих функций равен 27r, т. е.
т = 27r.
Простейшим периодическим процессом (двюкением) является про стое гармони-ч.еское колебание (движение), описываемое функцией
у= А· sin(wt + <ро), |
(66.3) |
t ~ О. где А - амплитуда колебания, w |
-- -ч.астота, <р0 |
-- на-ч.альна.я |
фаза. |
|
|
Функцию такого вида (и ее график) называют простой гармони |
кой. Основным периодом функции (66 3) |
является Т = 27r, т. е. одно |
|
|
(;.} |
полное колебание совершается за промежуток времени 27r (w показы w
вает, сколько колебаний совершает точка в течение 211' единиц време
ни)
Проведем преобразование функции (66.3):
у = Аsin(wt + <ро) = А sin ""t cos <ро +А cos wt sin <ро = аcos""t + bsin""t,
(66.4)
где а = А sin <р0, Ь = А cos <ро. Отсюда видно, что простое гармоническое колебание описывается периодическими функциями sin""t и cos(l.lt.
Сложное гармони"tеское колебание, возникающее в результате на
ложения конечного (или бесконечного) числа простых гармоник, также
описывается функциями вида sin""t и cos(l.lt. Так, функция
<p(t) = Ао + А1 sin(t + <р1) + А2 sin(2t + <р2) + ... + А30 sin(30t + <р30) =
30
= Ао + L An siп(nt + <r?n)
n=l
30
или, что равносильно, функция <p(t) = Ао + Е (ancosnt + Ьnsinnt)
n=l
задает сложное гармоническое колебание. Так как период первой гар-
моники есть Т1 = 211', второй Т2 = |
2 |
2 |
{ , третьей Тз = |
:J,... ,тридцатой |