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

080100Экономика(МатАнализ и ЛинАлгебра) / Лекции_Математический анализ

.pdf
Скачиваний:
43
Добавлен:
26.02.2016
Размер:
1.35 Mб
Скачать

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

f (x) a0 + a1(x x0 ) +…+ an (x x0 )n .

Погрешность этого равенства равна значению суммы остатка ряда

α = rn (x) = an+1(x x0 )n+1 + an+2 (x x0 )n+2 +….

Погрешность легче оценить, если исходный степенной ряд удовлетворяет усло-

виям признака Лейбница. В этом случае α < an+1(x x0 )n+1 .

Примеры.

1) Найти приближенное значение sin18 и оценить погрешность. Решение. Для вычисления воспользуемся приближенным равенством

sin x x x3 , которое получается при оставлении в ряде Тейлора первых двух

3!

слагаемых.

Заметим, что эта формула для функции sin x будет справедлива,

только

если

аргумент

x

измеряется

в

радианах.

Так

как

18

 

= π 18

/ 180

= π / 10 0.31415 < 0.4 (рад.),

то

sin18

=sin(π/10)

 

 

π

 

π3

= 0.30899 .

 

 

 

 

 

 

10

6000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оценим погрешность, учитывая, что при x 1 ряд Тейлора для функции sin x удовлетворяет условиям признака Лейбница. Поэтому остаток ряда, можно оценить с помощью теоремы 3.2.2. По этой теореме сумма остатка не превосходит по модулю старшего из отброшенных членов. Следовательно, по-

грешность формулы

 

α

 

(0.4)5

= 0.000085 < 0.0001. Мы можем ручаться за

 

 

 

 

5!

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

четыре знака после запятой в нашем ответе. Следовательно, sin18

0.3090 .

2) Вычислить число e с точностью до 0.1.

 

 

 

Решение. Возьмем ряд Тейлора для функции ex и подставим в него x=1.

Тогда e =1 + 1 +

1

+

1

+

1

+

…=1 +1 + 1 +

1 +

1

+…. Возьмем

четыре пер-

 

 

 

24

1!

2!

3!

 

4!

 

2

6

 

 

вые слагаемые этого ряда. Тогда остаток ряда будет иметь требуемую оценку:

151

1

 

1

 

1

 

 

 

1

 

 

1

 

1

 

 

 

 

1

 

 

1

 

 

1

 

1

 

 

 

+

 

+

 

+…

=

 

 

 

+

 

 

+

 

 

+…<

 

1

+

 

 

+

 

 

+…+

 

+…

=

 

 

 

4!

4! 5

 

 

24

4

4

2

4n

4! 5! 6!

 

 

 

 

 

4! 5 6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

1

 

 

 

 

1

 

=

1

 

4

=

 

 

1

 

< 0.1.

 

 

 

 

 

 

 

 

24

1 1 / 4

24

3

18

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значит, достаточно взять сумму первых четырех слагаемых. Тогда e 1 +1 + 12 + 16 2.7 .

6.1.4.2. Вычисление интегралов

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

x

 

Si(x) = sin x dx ,

0

x

 

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

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

 

 

 

sin x

=1

x2

+

 

x4

... + (1)k

 

x2k

 

+....

 

 

 

 

x

 

 

 

(2k +1)!

 

 

 

3!

5!

 

 

 

 

 

Интегрируя обе части этого равенства, получим разложение интегрально-

го синуса в ряд:

 

 

 

 

 

 

 

 

 

 

 

 

 

x sin x

 

x3

x5

 

 

 

(1)k x2k+1

(1)k x2k+1

x

dx = x

 

+

 

... +

 

+... =

 

.

3! 3

5! 5

(2k +1)! (2k +1)

(2k +1)! (2k +1)

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0

 

 

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

6.1.4.3. Решение дифференциальных уравнений

Пусть требуется решить задачу Коши для дифференциального уравнения 1-го порядка

y′ = f (x, y); y(x0 ) = y0 .

152

Будем искать решение y(x)

в виде ряда по степеням (x x0 ), предпола-

гая, что такое разложение возможно в некоторой окрестности точки x0 :

y (x) = y (x ) +

y(x0 )

(x x ) +

 

y′′(x0 )

(x x )2 +

y′′′(x0 )

(x x )3 +….

 

 

 

 

0

1!

 

0

 

2!

0

3!

0

 

 

 

 

 

 

Искомые величины

y(x0 ) ,

y′′(x0 ) , y′′′(x0 ) , … находятся последователь-

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

Пример. Дано дифференциальное уравнение y′ = xy2 +1 с начальным условием y(1) =1. Найти четыре первых члена разложения решения в ряд.

Решение. Так как x0 =1, то по общей формуле имеем

y (x) = y (1) + y1!(1) (x 1) + y′′2!(1) (x 1)2 + y′′′3!(1) (x 1)3 +….

Значит для записи первых четырех членов ряда надо найти y(1) , y(1) , y′′(1) и

′′′

 

 

 

 

 

 

 

 

 

 

y (1) . Из начального условия знаем y(1) =1. Подставив значения x =1 и y =1 в

правую часть уравнения

y

= xy

2

+1 получим

2

 

 

y (1)

=1 1 +1 = 2. Далее, про-

дифференцируем обе части уравнения:

y′′ = y2 + 2xyyи подставим сюда значе-

ния x =1, y =1 и y

= 2 .

Тогда

 

′′

2

1 2

=5 . Аналогично, найдем

 

 

y (1)

=1 + 2 1

третью производную и ее значение при x=1:

y′′′(x) = 2yy′+ 2yy′ + 2x( y)2 + 2xyy′′, y′′′(1) = 2 1 2 + 2 1 2 + 2 1 22 + 2 1 1 5 = 26

и получим ответ:

y (x) 1 + 2(x 1) + 52 (x 1)2 +133 (x 1)3 .

6.1.5.Гармонический анализ

6.1.5.1.Периодические функции, их свойства, гармоники

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

153

ские ряды были введены Д. Бернулли в 1753 г. в связи с изучением колебаний струны. Однако вопрос о возможности такого разложения вызвал горячие споры среди математиков. Строго обоснован этот метод был в работах Фурье и Дирихле в двадцатых годах XIX века.

II.4.1.1.1.10. Определение периодической функции

Функция f (x) называется периодической, если существует постоянное число T > 0 , для которого f (x +T) = f (x) , каково бы ни было число x из обла-

сти определения этой функции. Число T c таким свойством называется периодом функции f (x) . Наиболее известными периодическими функциями являют-

ся тригонометрические функции: sin x , cos x , tg x , ….

II.4.1.1.1.11. Свойства периодической функции

Сумма, разность, произведение и частное функций периода T, очевидно, всегда дают функции того же периода.

Если построить график периодической функции y = f (x) для какого-

нибудь отрезка [a,a +T ] значений x, то полный график этой функции получит-

ся периодическим повторением построенного (рис. 6.1.1). y

a T

0

a

a + T

a + 2T x

РИС. 6.1.1

154

6.1.5.2. Гармоники, тригонометрический многочлен и тригонометрический ряд

II.4.1.1.1.12. Гармоники

Простейшей, и в то же время очень важной для приложений, является периодическая функция y = Asin(ωx + ϕ), где A, ω, ϕ – постоянные.

Эту функцию называют гармоникой с амплитудой A , частотой ω и

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

Гармоника имеет период T = 2πω. Пользуясь известной формулой синуса суммы, напишем:

Asin(ωx + ϕ) = A(cosωx sinϕ+ sinωx cosϕ) .

(85)

Положив

 

a = Asin ϕ, b = Acosϕ,

(86)

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

 

a cosωx + bsin ωx .

(87)

Обратно, всякая функция этого вида есть гармоника. Чтобы убедится в этом, достаточно найти А и ϕ из (86). При этом получим:

A = a2 +b2 , tg ϕ= ba .

В дальнейшем для гармоник мы будем пользоваться записью вида (87). Немного изменим запись гармоники, выразив частоту через период T. Положим T = 2 . Тогда вследствие равенства T = 2π / ω получим:

ω= 2Tπ = π ,

и, следовательно, гармоника с периодом T = 2 может быть записана в виде: a cos πx + bsin πx .

II.4.1.1.1.13. Тригонометрический многочлен и тригонометрический ряд

Рассмотрим последовательность гармоник

a

n

cos

πnx + b

sin πnx

, n =1,2,,

(88)

 

 

n

 

 

 

 

 

 

155

 

 

частоты которых кратны частоте первой из них. Их периоды T = 2π

ωn

= 2

n

.

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, T = 2

= nTn

 

и число T = 2 является периодом для всех гармо-

ник сразу.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Всякая сумма вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

πnx

+ bn sin

πnx

 

 

 

 

 

 

SN (x) = A + an cos

 

 

 

,

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

где A = const , будучи суммой функций периода

2

, есть функция того же периода. Функцию

SN (x)

будем называть тригонометрическим многочленом порядка N (периода 2 ). Число А принято обозначать

 

a0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Бесконечная сумма гармоник вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

πnx

 

 

πnx

 

 

 

 

 

 

0

+ an cos

 

 

+ bn sin

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

II.4.1.1.1.14. Свойства сходящегося тригонометрического ряда

Скалярным произведением двух функций f (x) и g(x) , определенных

b

на отрезке [a,b] , называется число f (x)g(x)dx . Две функции называются ор-

a

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

этой системы ортогональны.

Теорема. Система функций

1, cos πx , sin πx , , cos πnx , sin πnx ,

ортогональна на отрезке [, ]. Кроме того,

156

cos2 πnx dx = , sin2 πnx dx = .

6.1.5.3.Ряд Фурье периодической функции

II.4.1.1.1.15. Определение

Пусть функция f (x) разложена в тригонометрический ряд:

 

a

 

πnx

 

πnx

f (x) =

0

+ an cos

 

+ bn sin

.

2

 

 

n=1

 

 

 

 

Наша задача найти коэффициенты a0 , an , bn . Чтобы найти, например, an

умножим скалярно обе части последнего равенства на cos πnx . По свойству ор-

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

 

f (x) cos

πnx

 

πnx

cos

πnx

= an .

 

 

= an cos

 

 

 

 

 

 

 

 

 

 

Выразив отсюда коэффициент an , получим формулу

an = 1 f (x)cos πnx dx .

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

a0 = 1 f (x)dx,

an = 1

f (x)cos πnx dx, bn = 1 f (x)sin πnx dx .

 

 

 

Эти формулы называются формулами Эйлера-Фурье.

Тригонометрический ряд

 

 

 

a

 

πnx

 

πnx

f (x) =

0

+ an cos

 

+ bn sin

,

 

 

2

n=1

 

 

 

 

коэффициенты, которого определены формулами Эйлера-Фурье, называется рядом Фурье.

157

II.4.1.1.1.16. Ряды Фурье четной и нечетной функций

Рассмотрим теперь, как отыскиваются коэффициенты ряда Фурье, если разлагаемая функция обладает свойством четности или нечетности. Нетрудно показать следующие свойства:

– если F(x) – четная функция, то F(x)dx = 2F(x)dx ;

 

 

 

0

 

– если F(x) – нечетная функция, то

F(x)dx = 0.

 

 

 

 

 

 

а)

y

 

б)

y

 

 

 

 

 

0

x

 

0

x

 

 

 

 

 

РИС. 6.1.2

 

 

Геометрически эти утверждения очевидны, поскольку график четной функции симметричен относительно оси Oy (рис. 6.1.2, а), а нечетной функции

– относительно начала координат (рис. 6.1.2, б). Применим сказанное к формулам Эйлера-Фурье.

1)

 

Пусть периодическая функция

f (x)

четная. Тогда функция

f (x)cos

πnx

тоже четна (при смене знака

x оба множителя не меняются), а

 

функция

f (x)sin

πnx

 

нечетна (при смене знака

x второй множитель меняет

 

знак). Значит, bn = 0 ,

а в формулах для a0

и an

промежуток интегрирования

можно уменьшить вдвое.

Вывод. Если функция f (x) четна, то она разлагается в ряд косинусов,

(bn = 0) причем справедливы упрощенные формулы

158

a0 = 2 f (x)dx, an = 2 f (x)cos πnx dx .

(89)

2) Пусть теперь периодическая функция f (x) нечетная. Тогда функция

f (x)cos πnx также нечетна, а функция f (x)sin πnx четна (при смене знака x

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

a0 = 0 , an = 0 ), причем для вычисления bn справедлива упрощенная формула

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

(90)

 

Примеры.

1) Периодическая функция на интервале длиной в период задается равен-

ством f (x) =1 + 4x , x (2,2) . Разложить эту функцию в ряд Фурье.

Очевидно, T = 2 = 4 , = 2 . Функция не является четной или нечетной, поэтому используем общие формулы Эйлера-Фурье:

 

 

 

 

 

1

2

 

 

x

1

 

 

x2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

=

 

 

1

+

 

 

dx =

 

x +

 

 

 

 

 

= 2;

 

 

 

 

2

 

 

2

 

8

 

 

 

 

 

 

 

 

2

 

 

4

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

n

= 1

2 1

+

x

cos πnx dx , b

= 1

2 1 +

x

sin πnx dx .

 

2

 

 

4

 

 

2

 

n

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Для вычисления двух последних интегралов используем метод интегрирования по частям. В первом случае:

u =1

+

x

,

 

 

 

 

4

 

 

 

 

 

πnx

 

 

 

 

dx

dv = cos

2

 

 

 

 

Тогда

du

= dx

,

 

 

 

 

 

4

 

 

 

 

2

 

 

 

πnx

 

 

 

 

 

 

v =

 

 

 

sin

 

.

 

πn

2

 

 

 

 

 

 

 

1

 

 

 

x

2

 

πnx

 

2

1

 

2 2

 

πnx

 

1

 

 

1

 

πnx

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

a

=

 

 

1

+

 

 

 

sin

 

 

 

 

 

 

sin

 

dx =

 

0

+

 

cos

 

 

 

 

 

= 0 .

2

4 nπ

2

4

πn

2

2

π2n2

2

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При вычислении коэффициентов bn берем

159

u =1

+

x

,

 

 

 

 

4

 

 

 

 

 

πnx

 

 

 

 

dx

dv = sin

 

2

 

 

 

 

 

Получаем

du = dx

,

 

 

 

4

 

 

 

 

2

 

πnx

 

 

 

 

v = −

 

cos

 

.

πn

2

 

 

 

 

 

 

 

 

1

 

 

 

 

x

2

 

 

πnx

 

2

1

 

 

2 2

πnx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bn =

2

1

+

 

 

 

 

cos

 

2

 

 

 

+ 4

 

 

 

cos 2

 

dx

=

 

 

 

 

4

nπ

 

 

 

 

πn

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n

 

 

 

 

 

2

πnx

 

2

 

2

 

n

 

 

 

 

n+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

3

(1)

 

(1)

 

1

 

 

 

 

 

 

(1)

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

+

+

 

 

sin

 

 

 

= −

 

=

.

 

 

πn

 

 

 

πn

π2n2

2

 

 

πn

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

πn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Данная функция гладкая на отрезке [2,2], поэтому можно записать

(1)

n+1

 

πnx .

f (x) =1 +

 

sin

n=1

πn

 

2

Вынося за знак суммы множитель, независящий от n, получим

 

1

(1)

n+1

πnx .

f (x) =1 +

sin

 

π n=1

n

 

2

Заметим, что при x = 2 ряд Фурье имеет, очевидно, сумму 1, что соответствует теореме о достаточном условии разложимости функции в ряд Фурье. Из рис. 6.1.3 можно заметить, что

f (2 0) =

3

,

f (2 + 0)

=

1

S(2)

=

1

 

3

+

1

 

=1.

2

2

2

 

 

2

 

 

y

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

–2

0

 

 

2

 

 

6

 

 

 

 

 

8

 

x

 

 

 

 

 

РИС. 6.1.3

 

 

 

 

 

 

 

 

2) Периодическая функция на интервале длиной в период задается равенствами

x π,

при x (0,π),

 

 

2

 

f (x) =

3

 

 

π − x,

при x (π,2π).

 

2

 

 

 

160