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

pismennyy_d_t_konspekt_lekciy_po_vysshey_matematike_polnyy_k

.pdf
Скачиваний:
86
Добавлен:
12.05.2015
Размер:
11.99 Mб
Скачать

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

 

Можно показать, что полученное равенство справедливо при всех

х Е [-1; 1].

Ряды (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- -

42. 2

- ... - 4n+l

. n + 1 - ... '

если -1 < -~ ~ 1, т. е. -4 ~ х < 4.

 

470

Пример 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 + ... + апх?.

Точность этого равенства увеличивается с ростом п. Абсолютная по­

грешность этого приближенного равенства равна модулю остатка ряда,

т.е.

где

471

Таким образом, ошибку 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) = п!~п'

472

т. е. тп(х) < - 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)

 

 

 

473

Интегрируя обе части равенства (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 + ··

 

 

который сходится на всей числовой оси.

474

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) =

475

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уу',

у<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 = Уо·

476

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

ряд (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!

+ -

7! + ... ,

 

т.е. у= sinx.

 

Глава 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)

аа

Сдругой стороны,

 

J f(x).dx =

J f(x) dx +

J f(x) dx.

(66.2)

Ь

Ь

а+Т

 

478

ь+т

Но J

а+Т

ь

ь

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,... ,тридцатой

479