
080100Экономика(МатАнализ и ЛинАлгебра) / Лекции_Математический анализ
.pdf
Приближенно значения функции можно вычислить по формуле
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 = 2∫F(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

