Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Специальные главы математического анализа сборник задач. Учебное пособие.pdf
Скачиваний:
0
Добавлен:
11.08.2026
Размер:
650 Кб
Скачать

Список литературы

1.Бицадзе А. В. Основы теории аналитических функций комплексного переменного / А. В. Бицадзе. Москва : Наука, 1969. 240 с.

2.Евграфов М. А. Аналитические функции / М. А. Евграфов. Москва : Наука, 1968. 472 с.

3.Сборник задач по теории аналитических функций / М. А. Евграфов, Ю. В. Сидоров, М. В. Федорук, М. И. Шабунин,

Ê.А. Бежанов. Москва : Наука, 1969. 388 с.

4.Пантелеев А. В. Теория функций комплексного переменного и операционное исчисление в примерах и задачах: учебное пособие / А. В. Пантелеев, А. С. Якимова. Москва : Высш. шк., 2001. 445 с.

5.Привалов И. И. Введение в теорию функций комплексного переменного / И. И. Привалов. Москва : Физматгиз, 1960. 444 с.

6.Лаврентьев М. А. Методы теории функций комплексного переменного / М. А. Лаврентьев, Б. В. Шабат. Москва : Наука, 1973. 736 с.

7.Лунгу К. Н. Сборник задач по высшей математике. 2 курс /

Ê.Н. Лунгу и др. Москва : Айрис-пресс, 2006. 592 с.

8.Письменный Д. Т. Конспект лекций по высшей математике / Д. Т. Письменный. Москва : Айрис-пресс, 2006.

9.Буров А. Н. Практикум по спецглавам математики: учебное пособие / А. Н. Буров, Н. Г. Вахрушева, С. В. Клишина Новосибирск: Èçä-âî ÍÃÒÓ, 2001. 102 ñ.

10.Воробь¼в Н. Н. Теория рядов / Н. Н. Воробь¼в. Москва : Наука, 1970. 204 с.

35

Приложение 1

Комплексные числа, операции над комплексными числами

Запись числа z в виде z = x + iy называют алгебраической

формой комплексного числа.

Запись числа z в виде z = |z|(cos φ + i sin φ) называется

тригонометрической формой.

Запись числа z в виде z = |z|eèëè z = |z|ei arg z называют

показательной (или экспоненциальной ) формой комплексного

числа.

 

 

 

 

 

 

p

 

 

φ = x.

Аргумент z можно

найти,

используя

 

 

Модуль определяется по формуле |z| =

x2 + y2.

 

y

Òàê êàê −π < arg z ≤ π, òî

 

 

 

 

 

формулу tg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ïðè x < 0, y > 0,

 

 

 

arctg y/x + π

 

 

 

φ = arg z =

arctg y/x ïðè x > 0,

 

 

 

 

 

y/x π

x < 0, y < 0,

 

 

 

 

arctg

 

 

ïðè

 

 

 

 

 

 

 

π/2

 

 

x = 0, y > 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−π/2

ïðè

x = 0, y < 0,

 

 

 

 

 

 

 

 

 

не определ¼н при x = 0, y = 0.

Сравнение: a + bi = c + di означает, что a = c и b = d

(два комплексных числа равны между собой тогда и только тогда, когда равны их действительные и мнимые части).

Сложение: (a + bi) + (c + di) = (a + c) + (b + d)i.

Вычитание: (a + bi) − (c + di) = (a − c) + (b − d)i.

Умножение:

(a + bi)(c + di) = ac + bci + adi + bdi2 = (ac − bd) + (bc + ad)i.

Деление:

a + bi

=

(a + bi)(c − di)

=

ac + bd

c + di

 

 

 

 

(c + di)(c − di) c2 + d2

 

 

1

a

b

В частности, a + bi = a2 + b2 a2 + b2 i.

bc − ad + c2 + d2 i.

36

Операции >, <, ≥, ≤ на множестве комплексных чисел не определены.

Формула Муавра для возведения комплексных чисел в натуральную степень:

zn = |z|n(cos nφ + i sin nφ).

Корень n-ой степени из комплексного числа

n z = |z|1/n

cos n

+ i sin

 

n

 

=

 

 

 

 

 

 

 

 

Arg z

 

Arg z

 

 

 

 

= |z|1/n cos

φ + 2πk

+ i sin

φ + 2πk

 

= |z|1/nei

φ+2πk

.

 

n

n

n

Элементарные функции комплексного переменного

 

 

 

 

 

z2

 

 

 

zn

 

 

 

zk

w = ez = 1 + z +

 

+ . . . +

 

+ . . . =

X

 

,

 

 

 

 

 

 

 

 

 

 

 

2!

 

 

 

 

 

n!

 

 

 

k=0

k!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k z2k

 

 

 

 

k z2k+1

X

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

cos z = (−1)

 

(2k)!

,

 

sin z =

(−1)

(2k + 1)!

.

 

k=0

 

 

 

 

 

 

 

 

 

k=0

 

 

 

 

 

 

 

 

 

Формулы Эйлера

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eiz = cos z + i sin z,

cos z =

eiz + e−iz

,

sin z =

eiz − e−iz

.

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2i

Гиперболические функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch z =

ez + e−z

, sh z =

ez − e−z

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

th z =

sh z

, cth z =

ch z

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch z

 

 

 

 

 

sh z

 

 

 

 

 

 

cos iz = ch z, ch iz = cos(−z) = cos z;

sh z = −i sin iz, sh iz = i sin z, sin iz = ish z.

37

Множество значений логарифма

Ln z = ln |z| + i(arg z + 2kπ), k Z.

При k = 0 главное значение

ln z = ln |z| + i arg z, −π < arg z < π.

Обратные тригонометрические и гиперболические функции

arccos z = i Ln (z +

z2

1);

 

 

 

 

π

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arcsin z =

 

 

− arccos z = −i Ln (iz + p1 − z2);

2

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

i

 

 

1 + iz

 

 

arctg z =

 

− arcctg z = −

 

Ln

 

 

;

 

2

2

 

1 − iz

 

 

 

 

arcctg z = −

 

i

 

 

z + i

 

 

 

 

 

 

 

 

 

 

 

Ln

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

2

z − i

 

 

 

 

 

 

arcsh z = Ln (z +

 

 

 

 

 

 

 

 

z2 + 1);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arcch z = Ln (z +

 

 

 

 

z2 − 1);

 

 

 

 

 

 

 

 

1

 

 

 

 

1 + z

 

 

 

 

 

 

 

 

 

 

arcth z =

 

Ln

p

 

;

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 − z

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

z + 1

 

 

 

 

 

 

 

 

 

 

arccth z =

 

 

Ln

 

.

 

 

 

 

 

 

 

 

 

 

2

z − 1

 

 

 

 

 

 

Общая степенная функция w = za,

 

a = α + i β

 

 

w = za = ea Ln z.

Учитывая z = r e, Ln = ln r + i(φ + 2kπ), получим

za = ea ln r−β(φ+2kπ)ei(α(φ+2kπ)+β ln r), k Z.

Дифференцирование функций комплексного переменного

Для того, чтобы функция f(z) = f(x, y) = u(x, y) + iv(x, y) была дифференцируемой в точке z = x + iy D, необходимо и достаточно, чтобы выполнялись следующие условия:

38

1)вещественнозначные функции u(x, y) и v(x, y) должны быть дифференцируемы в точке z;

2)в этой точке должны выполняться равенства:

∂u

=

∂v

,

∂u

= −

∂v

,

 

 

 

 

∂x

∂y

∂y

∂x

называемые условиями Коши Римана.

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

f(z) =

∂u

+ i

∂v

, f(z) =

∂v

− i

∂u

,

 

 

 

 

 

 

 

 

 

∂x

∂x

∂y

 

∂y

f(z) =

∂u

− i

∂u

, f(z) =

 

∂v

+ i

 

∂v

.

 

 

 

 

 

 

 

∂x

 

∂y

 

∂y

∂x

Нахождения аналитической функции по заданной е¼ действительной или мнимой части

Правило 1

1.Найти частные производные до второго порядка заданной функции двух переменных: u(x, y) или v(x, y). Проверить,

если требуется, что заданная функция гармоническая, т. е. в некоторой области D выполняется равенство ∆u = 0 или

∆v = 0.

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

(x,yZ )

v(x, y) = −∂u∂ydx + ∂u∂xdy + C.

(x0,y0)

Если выбрать кривую в виде ломаной со звеньями, параллельными осям, то

x

 

 

 

y

 

 

 

v(x, y) = xZ

 

x, y

)

dx + yZ

∂u(x, y)

 

∂u( 0

 

 

 

dy + C.

∂y

 

∂x

0

 

 

 

0

 

 

 

39

Èëè

(x,y)

 

 

Z

 

∂v

 

∂v

 

u(x, y) =

 

 

dx −

 

 

dy + C =

 

∂y

∂x

 

(x0,y0)

 

 

 

 

x

 

 

 

 

 

y

 

 

 

 

= xZ

x, y

)

dx − yZ

∂v(x, y)

 

∂v(

0

 

 

 

dy + C.

∂y

 

 

 

∂x

0

 

 

 

0

 

 

 

 

 

3.Записать функцию f(z) = u + iv.

4.Если задано дополнительное условие значение функции f(z) в некоторой точке, то следует использовать его для отыскания C.

Правило 2

Если заданы начальные условия, функция восстанавливается однозначно.

1. Äàíî: u = u(x, y), f(z0) = C0. Тогда

 

2

 

 

 

2i

 

0

 

 

z +

z

0

 

z −

z

0

 

 

 

 

 

f(z) = 2u

 

 

 

,

 

 

 

 

 

C

.

2. Äàíî: v = v(x, y), f(z0) = C0. Тогда

 

2

 

 

 

2i

 

0

 

 

z +

z

0

 

z −

z

0

 

 

 

 

 

f(z) = 2iv

 

 

 

,

 

 

 

 

+ C

.

Интегрирование функций комплексного переменного

Криволинейный интеграл от функции f(z) по кривой l

Z Z Z

f(z)dz = udx − vdy + i vdx + udy.

l l l

40

 

Свойства криволинейного интеграла

 

от функции комплексного переменного

1. Rl (C1f1(z) + C2f2(z))dz = C1 Rl

f1(z)dz + C2

Rl

f2(z)dz.

R

 

R

 

 

 

 

2.

f(z)dz = − f(z)dz.

 

 

 

AB

R

BA

R

 

 

 

R

 

 

 

 

3.

f(z)dz = f(z)dz +

f(z)dz.

 

 

AB

AC

CB

 

 

 

R

 

 

 

 

 

 

4.|dz| = lAB длина дуги AB.

 

AB

 

 

R

 

R

 

AB, ò. å.

 

R

|f(z)||dz|. В частности,

5.

l

f(z)dz

l

l

f(z)dz

≤ M lAB,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

если функция ограничена по модулю на кривой

|f(z)| ≤ M, z l.

R

6.dz = zB zA.

l

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

переменная, обозначим е¼ z. Тогда значение интеграла будет за-

R

висеть только от точки z, т. е. f(z)dz определяет некоторую

l

функцию в указанной области:

Zz

f(ξ)dξ = F (z).

z0

Функция F (z) интеграл с переменным верхним пределом .

F (z) имеет производную в любой точке е¼ определения и, следовательно, является аналитической. При этом

F (z) = f(z).

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

41

Функция F (z), для которой выполняется равенство

Zz

f(ξ)dξ = F (z),

z0

называется первообразной для функции f(z) в односвязной области, а совокупность (множество) первообразных Φ(z) = F (z) + c, где c = const, неопредел¼нным интегралом от функции f(z).

Теорема Ньютона Лейбница. Если f(z) аналитическая в односвязной области, F (z) е¼ первообразная, т. е. F (z) = f(z),

C ориентированная кривая с начальной точкой z1 и конечной точкой z2, то выполняется равенство

ZZz2

f(z)dz = f(z)dz = F (z2) − F (z1),

C z1

т. е. интеграл не зависит от пути, соединяющего точки z1 è z2, а зависит только от функции и конечных точек.

Интеграл

 

Z

f( )

1

 

 

 

ξ

dξ = F (z)

 

2πi

ξ

z

 

 

l

 

называется интегралом типа Коши.

Дифференцируя это равенство по переменной z получаем:

 

1

 

Z

 

 

 

)

 

 

F (z) =

 

 

 

 

f(ξ

dξ;

 

2πi

(ξ z)2

 

 

 

 

l

 

 

 

 

 

F (n)(z) =

 

n!

Z

 

 

f(ξ)

 

 

 

 

 

dξ.

2πi

z)n+1

 

 

 

 

l

 

 

 

 

 

42

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

Z

Первый способ. Вычисление интегралов f(z)dz от непрерывной

l

функции пут¼м сведения к криволинейным интегралам от функ-

ции действительных переменных с помощью формулы

Z Z Z

f(z)dz = udx − vdy + i vdx + udy.

l l l

1.Найти Re f(z) = u, Im f(z) = v.

2.Записать подынтегральное выражение

f(z)dz = (u + iv)(dx + idy) = udx − vdy + i(udy + vdx).

3. Вычислить криволинейные интегралы вида

Z P dx + Qdy ïî

 

l

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

Z

Второй способ. Вычисление интегралов f(z)dz от непрерыв-

l

ной функции пут¼м сведения к определ¼нному интегралу в случае параметрического задания пути интегрирования.

Если кривая l задана параметрически: z = z(t), α ≤ t ≤ β или

 

β

y = y(t) , α ≤ t ≤ β, то имеем Z f(z)dz = Z f(z(t))z(t)dt.

x = x(t)

 

l

α

1.Записать параметрическое уравнение кривой z = z(t) и из него определить пределы интегрирования: t = α соответствует начальной точке пути интегрирования, t = β конечной.

2. Найти дифференциал комплекснозначной функции z(t): dz = z(t)dt.

43

3. Подставить z(t) в подынтегральное выражение, преобразо-

вать интеграл к виду Zβ f(z(t))z(t)dt =

Zβ φ(t)dt.

α

α

4.Вычислить полученный определ¼нный интеграл от комплекснозначной функции действительной переменной.

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

1.Найти первообразную F (z), используя свойства интегралов,

табличные интегралы и методы, известные из действительного анализа.

2.Применить формулу Ньютона Лейбница:

ZZz2

f(z)dz = f(z)dz = F (z2) − F (z1).

C z1

Замечание

1.В случае многосвязной области проводятся разрезы так, чтобы можно было получить однозначную функцию F (z).

2.При интегрировании однозначных ветвей многозначных функций ветвь выделяется заданием значения функции в некоторой точке кривой интегрирования. Если кривая замкнутая, то начальной точкой пути интегрирования считается та точ- ка, в которой задано значение подынтегральной функции. Значение интеграла может зависеть от выбора точки.

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

 

1

I

f( )

f(z) =

 

 

ξ

dξ, z D.

2πi

ξ

z

 

 

C

 

44

ψ(z)

2.Аналитическая функция имеет производные любого порядка, для которых справедлива формула ( обобщ¼нная инте-

гральная формула Коши )

f(n)(z) =

n!

I

 

f(ξ)

 

 

 

 

dξ.

2πi

z)n+1

 

 

C

 

 

Вычисление интегралов по замкнутому контуру от функций

комплексного переменного

 

Будем рассматривать интегралы вида

I

ψ(z)dz, где функция

 

 

φ(z)

 

C

 

 

 

¯

 

 

 

 

φ(z) аналитическая в D, а ψ(z) многочлен, не имеющий нулей

на контуре C.

 

 

 

 

 

Правила

 

 

 

 

 

1. В области D нет нулей многочлена ψ(z). Тогда f(z) =

φ(z)

 

ψ(z)

 

 

 

 

 

функция аналитическая. Применяя основную теорему Коши,

имеем

I

φ(z)

dz = 0.

 

 

 

 

C

2.В области D расположен один простой нуль z = a много- члена ψ(z). Тогда записываем подынтегральную функцию

 

φ(z)

 

f(z)

 

 

 

 

¯

â âèäå:

ψ(z)

 

=

z − a

, где f(z) функция аналитическая в D.

Получим

 

I

φ(z)

 

I

f(z)

 

 

 

 

 

 

 

 

 

dz =

 

 

 

dz = 2πi f(a).

 

 

 

ψ(z)

(z

a)

 

 

 

C

 

 

 

 

C

 

 

 

3. В области D расположен один кратный нуль z = a много-

члена ψ(z) (кратности n). Тогда

φ(z)

=

 

 

f(z)

 

f(z)

ψ(z)

(z − a)n , ãäå

функция аналитическая в

¯

 

 

 

I

 

 

I

 

 

D:

 

 

 

 

 

 

 

 

φ(z)

f(z)

 

 

2 i

 

 

dz =

 

 

 

dz =

 

π

 

f(n−1)(a).

 

ψ(z)

(z

a)n

(n

1)!

 

C

C

 

 

 

 

 

 

 

 

 

 

45

4. В области D два нуля многочлена

ψ(z): z1 = a, z2 = b.

В этом случае интеграл запишем в виде:

I f(z)dz =

I f(z)dz +

I f(z)dz,

C

C1

C2

ãäå C1 è C2 границы непересекающихся окрестностей точек z1 = a, z2 = b. Каждый из полученных интегралов вычисляем в соответствии с п. 2 и 3. Очевидно, можно рассматривать случаи большего числа нулей ψ(z) в области D.

Функциональные ряды в комплексной области Основные табличные разложения (Маклорена):

 

z

 

z2

 

 

z3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zn

 

 

 

 

 

 

 

zk

ez = 1 + 1! +

2! +

 

3! + · · ·

 

 

+ n! + . . . =

 

P

 

k! , (R = ∞);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0

 

 

 

 

 

ln(1 + z) = z −

 

 

z2

 

 

 

 

z3

− · · · + (−1)

n+1 zn

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

+ . . . =

2

 

3

 

 

 

 

n

 

 

 

 

 

=

 

 

 

 

(−1)k−1zk , ( z

 

< 1);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=1

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin z = z

 

z3

 

 

+

z5

 

 

 

 

 

 

 

 

 

 

+

 

 

(−1)n

z2n+1

+ . . . =

3!

 

 

 

 

 

 

 

 

 

 

 

(2n + 1)!

 

 

 

 

 

 

 

 

 

 

 

 

5! − · · ·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

(

1)k

 

z2k+1, (R = );

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0 (2k + 1)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

cos z = 1 −

z

+

 

z

− · · · +

(−1)

z2n + . . . =

2!

4!

 

(2n)!

 

 

 

 

 

=

 

 

 

 

 

 

(

 

 

 

1)k

z2k, (R =

 

 

);

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0 (2k)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arctg z = z

 

z3

 

 

+

z5

 

 

 

 

 

 

 

+

 

(−1)n

z2n+1

+ . . . =

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 − · · ·

 

 

2n + 1

 

 

 

 

 

 

 

 

=

 

(−1)k z2k+1,

( z

 

 

 

1);

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

 

| ≤

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0 2k + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

46

 

z3

 

 

 

 

z5

 

 

 

 

 

1

 

 

2n+1 + . . . =

sh z = z +

 

+

 

 

+ · · · +

 

 

 

z

3!

5!

(2n + 1)!

=

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

P

 

(2k + 1)!z2k+1, (R = ∞);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2

z4

1

z2n + . . . =

ch z = 1 +

 

+

 

 

+ · · · +

 

2!

4!

(2n)!

 

=

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

(2k)!z2k, (R = ∞);

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

k=0

(1 + z)α = 1 + αz + α(α − 1)z2 + . . .

2!

+ α(α − 1) · · · (α − n + 1)zn + . . . , (|z| < 1) (биномиальный ряд); n!

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1 + z + z2

+ z3 + · · · + zn + . . . = zk, (|z| < 1);

1

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0

 

 

= 1 +

z

 

z2

 

+

z3

 

 

(−1)n(2n)!

zn + . . . =

1 + z

 

 

+

 

8

 

16 − · · ·

 

 

 

 

 

2

 

 

 

(1

2n)n!24n

 

 

 

 

 

 

 

 

 

(−1)

k

(2k)!

 

 

 

 

 

 

 

 

=

 

 

zk, ( z < 1).

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

(1

2k)k!24k

 

 

| |

 

 

 

 

 

 

 

 

 

k=0

 

 

 

 

 

 

 

 

 

 

 

 

Функция f(z) аналитическая в кольце, представляется в этом кольце сходящимся рядом по целым степеням, т. е. имеет место

равенство

 

 

 

 

 

X

cn(z − z0)n.

 

 

f(z) =

 

 

 

 

 

n=−∞

 

Коэффициенты ряда вычисляются по формуле

 

1

I

 

f(z)

 

cn =

 

 

 

 

dz, n = 0, ±1, ±2, . . . ,

2πi

(z

z0)n+1

 

 

γ

 

 

 

где γ произвольный контур, принадлежащий кольцу и охватывающий точку z0, в частности γ окружность |z − z0| = ρ, r < ρ < R.

47

Вычетом функции f(z) в изолированной особой точке

надлежащий окрестности точки zHγ

и охватывающий е¼.

¯

1

 

 

z0 C называется интеграл

2πi

f(z)dz, где γ контур, при-

 

0

 

Обход контура положительный, т. е. область, им ограниченная и принадлежащая окрестности z0, при обходе расположена слева: для z0 C обход против часовой стрелки, для z0 = ∞

по часовой. Обозначают вычет res f(z):

z0

I

res f(z) = 1 f(z)dz, γ O(z0) \ z0; 0 < |z − z0| < R;

z0 2πi

γ

I

res f(z) = 1 f(z)dz, γ O(∞) \ ∞; R < |z| < ∞.

2πi

γ

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

Теорема 1 (теорема о вычислении вычета). Вычет функции в изолированной особой точке равен коэффициенту c−1 ïðè

первой отрицательной степени в разложении функции в ряд

1

Лорана в окрестности этой точки, т. е. при z − z0 äëÿ z0 C, и этому коэффициенту, взятому с противоположным знаком,

äëÿ z0 = ∞:

res f(z) = c−1, z0 C;

z0

res f(z) = −c−1, z0 ∞.

При этом в зависимости от характера особой точки z0 верны следующие утверждения:

1)вычет в устранимой особой точке равен нулю;

2)для простого полюса (Π(1))

res f(z) = lim (f(z)(z − z0));

z0 z→z0

48

3) для полюса порядка n (Π(n), n > 1)

 

1

 

dn−1

n

res f(z) =

 

lim

 

(f(z)(z − z0) );

 

dzn−1

z0

(n − 1)! z→z0

4) если функция f(z) в окрестности точки z = z0 предста-

âèìà â âèäå f(z) = φ(z), прич¼м φ(z) и ψ(z) аналитиче-

ψ(z)

ские функции, такие что φ(z0 ̸= 0), ψ(z0) = 0, ψ(z0) ̸= 0 (z0 простой полюс), то

res f(z) =

φ(z0)

;

 

z0

ψ(z0)

5)для существенно особой точки вычет находится из разложения Лорана в окрестности точки z0:

res f(z) = c−1, z0 C;

z0

res f(z) = −c−1, z0 ∞.

С помощью вычетов можно записать в другой форме основную теорему Коши для сложного контура.

Теорема 2 (теорема Коши о вычетах). Если функция f(z)

аналитическая в

¯

 

 

 

 

D, за исключением конечного числа особых то-

÷åê zk D, то справедливо равенство

 

 

 

n

 

 

 

C

X

res f(z), z

 

D,

I f(z)dz = 2πi k=1

k

zk

 

где C граница области D.

Теорема 3 (обобщ¼нная теорема о вычетах). Сумма выче- тов функции f(z) во всех е¼ особых точках, включая бесконечно

удал¼нную точку, равна нулю:

n

X

res f(z) + res f(z) = 0,

zk

k=1

где C граница области D.

49

Вычисление несобственных интегралов с помощью вычетов

Рассмотрим интегралы вида

R(x)dx, ãäå R(x) =

Pn(x)

 

 

Q (x)

−∞

m

R

 

дробно-рациональная функция; Pn(x), Qm(x) многочлены степени n и m.

Pn(x)

Теорема 4. Пусть R(x) = Qm(x) удовлетворяет условиям:

1)m −n ≥ 2, т. е. степень знаменателя больше степени числителя, по крайней мере, на два;

2)Qm(x) ̸= 0 для x действительных, т. е. R(z) не имеет особых точек на действительной оси.

Тогда справедливы равенства:

p

 

 

 

−∞

X

 

 

 

Z

R(x)dx = 2 i

res R(z), Im z > 0,

π

k=1

zk

k

 

 

p

 

 

 

−∞

 

X

res R(z),

 

< 0,

Z

R(x)dx = −2πi k=1

Im zk

zk

 

ãäå zk, k = 1, 2, . . . , p все особые точки функции R(z), расположенные выше оси Ox (Im zk > 0) для первой формулы и ниже оси Ox (Im zk < 0) во втором случае.

Замечание

Если R(x) ч¼тная функция, то можно, используя эти фор-

мулы, вычислять интегралы вида

R0

f(x)dx, òàê êàê äëÿ ÷¼òíûõ

 

 

 

функций справедливо равенство RR f(x)dx = 2 R0

f(x)dx.

50

 

Алгоритм вычисления интегралов

R

R(x)dx с помощью

−∞

вычетов

1.Проверить условия применения формул: m − n ≥ 2, т. е.

степень знаменателя больше степени числителя, по крайней мере, на два; Qm(x) ̸= 0 для x действительных, т. е. R(z) не имеет особых точек на действительной оси.

2.Найти особые точки подынтегральной функции R(z).

3.Вычислить вычеты в особых точках функции R(z), расположенных

a) выше оси Ox, если применяется формула

 

 

 

 

 

p

 

 

 

 

 

 

 

−∞

 

 

 

X

 

 

 

 

> 0;

 

 

 

 

R(x)dx

 

 

res R(z), Im z

 

 

 

Z

 

= 2πi

k=1

zk

 

 

k

 

b) ниже оси Ox, если применяется формула

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

−∞

 

 

 

X

res R(z),

 

 

 

 

 

Z

R(x)dx = −2πi k=1

Im zk < 0.

 

 

zk

 

4. Записать результат по соответствующей формуле.

Рассмотрим интегралы вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

R(x) cos λxdx,

Z

 

R(x) sin λxdx,

Z

R(x)eiλxdx,

−∞

 

 

 

 

−∞

 

 

 

−∞

 

 

 

ãäå R(x) =

Pn(x)

, m − n ≥ 1 äëÿ x R.

 

 

 

 

Qm(x)

 

 

 

 

 

51

Теорема 5. Пусть функция R(x) рациональная функция, не имеющая особых точек на действительной оси (т. е. Q(x) ̸= 0 для x R), для которой z = ∞ нуль порядка не ниже первого (т. е. m − n ≥ 1). Тогда справедливы формулы:

1) ïðè λ > 0

p

R(x)eiλxdx = 2

 

i res R(z)eiλz, Im z

 

> 0;

Z

π

k=1

zk

k

 

−∞

 

X

 

 

 

2) ïðè λ < 0

p

R(x)eiλxdx =

2

 

i

 

res R(z)eiλz, Im z

 

< 0;

Z

π

 

k=1

zk

k

 

−∞

 

 

 

X

 

 

 

3) ïðè λ > 0

 

 

 

 

p

!

 

 

 

 

λ

 

− π

zk

 

 

k

 

−∞

xdx =

2

Im

X

 

, Im z

 

> 0;

R(x) cos

res R(z)eiλz

 

Z

 

 

 

 

k=1

 

 

 

 

 

 

 

 

 

p

!

 

 

 

 

−∞

λ

 

π

 

zk

 

k

 

 

 

xdx = 2

Re

X

, Im z

 

> 0.

R(x) sin

 

res R(z)eiλz

 

Z

 

 

 

 

k=1

 

 

 

 

 

Алгоритм вычисления интегралов

Z

R(x)eiλxdx

 

 

 

−∞

 

1.Проверить условия для функции R(x), сформулированные в теореме 5: R(x) рациональная функция, не имеет особых

точек на действительной оси (т. е. Q(x) ̸= 0 для x R), для которой z = ∞ нуль порядка не ниже первого (т. е. m − n ≥ 1).

52

2.Вычислить вычеты функции R(z)eiλz во всех е¼ особых точках:

a) лежащих выше оси Ox в случае λ > 0; b) лежащих ниже оси Ox в случае λ < 0.

3. Записать результат по формуле

 

p

 

 

 

 

 

 

 

 

−∞

iλx

X

 

 

 

iλz

 

 

 

 

Z

 

R(x)e dx = 2πi k=1

res

R(z)e , Im zk

> 0, ïðè

λ > 0

 

zk

 

 

 

и по формуле

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

−∞

 

iλx

X

 

 

iλz

 

 

 

 

Z

R(x)e dx = −2πi k=1

res

R(z)e , Im zk

< 0, ïðè

λ < 0.

zk

 

 

Z

 

4. Для вычисления интегралов

 

Z

R(x) cos λxdx è

R(x) sin λxdx

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

−∞

 

 

в результате, полученном в п. 3, отделить действительную и мнимую части и записать результаты соответственно по формулам:

 

 

 

 

p

!

 

 

 

 

λ

 

− π

zk

 

 

k

 

−∞

xdx =

2

Im

X

 

, Im z

 

> 0,

R(x) cos

res R(z)eiλz

 

Z

 

 

 

 

k=1

 

 

 

 

 

 

 

 

 

p

!

 

 

 

 

−∞

λ

 

π

 

zk

 

k

 

 

 

xdx = 2

Re

X

, Im z

 

> 0.

R(x) sin

 

res R(z)eiλz

 

Z

 

 

 

 

k=1

 

 

 

 

 

53

Операционное исчисление

Оригинал комплекснозначная функция f(t) действительного аргумента t, которая удовлетворяет следующим условиям:

1)f(t) = 0, t < 0;

2)f(t) локально интегрируема, т. е. на любом конечном про-

межутке [a, b], 0 ≤ a < b существует интеграл

Ra f(t)dt (èëè

 

b

на любом конечном отрезке [a, b] f(t) удовлетворяет условиям Дирихле: a) ограничена; b) либо непрерывна, либо имеет конечное число разрывов 1-го рода; c) имеет конечное число экстремумов);

3) f(t) имеет ограниченный рост, т. е. возрастает не быстрее по-

казательной функции: существуют такие постоянные M > 0 и σv≥ 0, что

|f(t)| < M eσvt, t > 0.

Изображение функции f(t) функция F (p) комплексного переменного p, определяемая равенством

 

 

 

 

F (p) = Z0

e−ptf(t)dt;

 

 

 

σv+i∞

 

f(t) = L−1 [F (p)] =

1

Z

eptF (p)dp.

 

2πi

σv−i∞

54

Свойства преобразования Лапласа

1.Линейность. Если f1(t), . . . , fn(t) оригиналы,

F1(p), . . . , Fn(p) соответствующие им изображения, то для лю-

ляется оригиналом и справедливо равенство

 

n

 

P

бых комплексных чисел ci, i = 1, . . . , n, функция

ckfk(t) ÿâ-

 

 

" n

ckfk(t)# = n ckL [fk(t)] =

 

 

k=1

 

L

n

ckFk(p),

 

 

X

X

X

 

 

 

 

k=1

k=1

k=1

 

 

 

 

 

Re p > max(σv1, . . . , σvn),

 

 

èëè

n

 

n

 

 

 

 

 

 

 

 

XX

ckfk(t)

ckFk(p), Re p > max(σv1, . . . , σvn).

k=1 k=1

Справедливо и обратное утверждение: если F1(p), . . . , Fn(p) изображения, то

L−1 [c1F1(p) + . . . + cnFn(p)] = c1f1(t) + . . . + cnfn(t).

2. Подобие (теорема подобия). Для любого a > 0 из

F (p) = L [f(t)] следует

L [f(at)] = 1F p , Re p > aσv,

a a

0

и обратно:

 

a .

L−1 [F (ap)] = af

1

 

 

t

3. Смещение (теорема смещения). При любом комплексном a из F (p) = L [f(t)] следует

L eatf(t) = F (p − a), Re(p − a) > σv,

0

т. е. умножение оригинала на eat соответствует смещению изображения на a.

55

4. Запаздывание (теорема запаздывания). Для любого τ > 0 из F (p) = L [f(t)] следует

L [f(t − τ)] = e−pτF (p), Re p > σv,

0

где f(t − τ) = f(t − τ) · η(t − τ), т. е. запаздывание оригинала на τ > 0 соответствует умножению изображения на e−pτ.

5. Дифференцирование оригинала. Если f(t), f(t), . . . , f(n)(t)

являются оригиналами и F (p) f(t), то

f(t) pF (p) − f(+0),

f′′(t) p2F (p) − pf(+0) − f(+0),

.

f(n)(t) pnF (p) − pn−1f(+0) − . . . − f(n−1)(+0),

ãäå f(k)(+0) = lim f(k)(t), k = 0, 1, 2, . . . , n − 1.

t→+0

6. Интегрирование оригинала. Если f(t) является оригиналом и F (p) f(t), то

Zt

f(τ)dτ F (p), Re p > σv,

p

0

0

т. е. интегрирование оригинала соответствует делению изображения на p.

7. Дифференцирование изображения. Если функция f(t)

является оригиналом и F (p) f(t), то

F (n)(p) (−1)ntnf(t), tnf(t) (−1)nF (n)(p).

Частные случаи: tf(t) −F (p), t2f(t) F ′′(p).

8. Интегрирование изображения. Если

f(t)

является ори-

t

гиналом и F (p) f(t), то

 

 

 

 

 

 

 

 

 

L f t

 

F (z)dz.

 

 

= Zp

 

 

 

(t)

 

 

 

 

 

56

9. Умножение изображений (теорема Бореля). Произведение двух изображений F1(p) è F2(p) также является изображением, прич¼м

F1(p) F2(p) f1(t) f2(t).

Функция f1(t) f2(t), определяемая формулой

t

t

f1(t) f2(t) = Z0 f1(τ)f2(t − τ)dτ =

Z0 f1(t − τ)f2(τ)dτ,

называется св¼рткой функций.

Теорема Бореля утверждает, что умножение изображений равносильно св¼ртке оригиналов: f(t) g(t) F (p) G(p).

10. Дифференцирование св¼ртки (интеграл Дюамеля): pF1(p)F2(p) f1(0)f2(t) + f1f2 =

= f2(0)f1(t) + f2f1 =

t

 

 

= f1(0)f2(t) + Z0

f1(τ)f2(t − τ)dτ =

= f2(0)f1(t) + Zt

f1(t − τ)f2(τ)dτ.

0

 

Начальные значения оригинала находятся по формуле

f(+0) = lim pF (p).

p→∞

Если существует конечный предел lim f(t) = f(+∞), то

t→+∞

f(+∞) = lim pF (p).

p→0

57

Нахождение изображений периодических функций

Пусть f(t) оригинал с периодом T , образованный повторением функции f0(t):

0, t ≤ 0,

f0(t) = φ(t), 0 < t ≤ T,0, t > T.

Алгоритм нахождения изображения F (p) периодической функции f(t)

1. Найти изображение функции f0(t):

f0(t) F0(p).

2. Найти изображение F (p) по формуле

F (p) = F0(p) .

1 − e−T p

Таблица преобразования Лапласа

 

f(t)

 

 

 

F (p)

 

f(t)

 

 

 

 

F (p)

ï/ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï/ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

η(t)

1

 

 

 

 

 

 

10

t sin ωt

 

 

 

 

2pω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

(p2 + ω2)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

C

 

 

 

 

 

C

 

11

t cos ωt

 

 

 

 

p2 − ω2

 

 

 

 

 

 

 

 

p

 

 

 

(p2 + ω2)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

t

1

 

 

 

 

 

 

 

eλt sin ωt

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

12

 

(p − λ)2 + ω2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

tn

 

 

 

 

 

n!

 

 

eλt cos ωt

 

 

 

 

p − λ

 

 

 

 

pn+1

13

 

(p − λ)2 + ω2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

δ(t)

1

 

 

 

 

 

 

14

t eλt sin ωt

 

 

 

2ω(p − λ)

 

 

 

 

 

 

 

((p − λ)2 + ω2)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

eλt

 

1

 

 

 

 

 

 

15

t eλt cos ωt

 

(p − λ)2 − ω2

 

 

 

p − λ

 

((p − λ)2 + ω2)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

tneλt

 

 

 

 

 

n!

 

sh ωt

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

(p − λ)n+1

16

 

 

 

 

p2 − ω2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

sin ωt

 

 

 

 

 

ω

 

ch ωt

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 + ω2

17

 

 

 

 

p2 − ω2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

cos ωt

 

 

 

 

 

 

p

 

t sh ωt

 

 

 

 

2pω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 + ω2

18

 

 

(p2 − ω2)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

59

 

 

 

 

 

 

 

f(t)

 

 

 

 

 

 

F (p)

 

 

 

 

 

 

 

 

 

 

 

 

 

f(t)

 

 

 

 

F (p)

ï/ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï/ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

19

 

 

 

 

 

t ch ωt

 

 

 

 

 

p2 + ω2

 

 

 

 

 

 

1

(eλt − 1 − λt)

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(p2 − ω2)2

 

 

 

28

 

 

 

λ2

 

 

p2(p − λ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Im(p + iω)n+1

 

 

 

 

 

 

 

 

 

 

λt2

 

 

 

 

 

p

 

 

 

 

 

20

 

 

 

 

tn sin ωt

 

 

 

 

 

 

 

 

 

 

n!

29

 

 

 

 

 

t +

 

eλt

 

 

 

 

 

 

 

 

 

 

 

 

(p2 + ω2)n+1

 

 

 

 

 

2

 

 

(p − λ)3

21

 

 

 

tn cos ωt

 

Re(p + iω)n+1

n!

30

 

 

 

 

 

 

(1 + λt) eλt

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

(p2 + ω2)n+1

 

 

 

 

 

 

 

 

(p − λ)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22

 

 

 

 

eλtsh ωt

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos2 ωt

 

 

p2 + 2ω2

 

 

 

 

 

 

 

 

(p − λ)2 − ω2

31

 

 

 

 

 

 

p(p2 + 4ω2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23

 

 

 

 

eλtch ωt

 

 

 

 

 

 

p − λ

 

 

 

 

 

 

 

 

 

 

 

 

 

sin2 ωt

 

 

 

 

2

 

 

 

 

 

 

 

(p − λ)2 − ω2

32

 

 

 

 

 

 

p(p2 + 4ω2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

t

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

ω3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

24

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

33

(sh ωt − sin ωt)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 +

λ

p

 

 

 

 

 

 

2

 

 

 

 

 

p4

4

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

ω2p

25

 

 

 

 

 

(eλt − 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

34

 

 

 

 

(ch ωt − cos ωt)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p(p

λ

)

 

 

 

 

 

2

 

 

 

 

 

p4

4

 

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

26

 

 

 

 

eλt − eμt

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

35

 

1

(sh ωt + sin ωt)

 

 

 

 

ωp2

 

 

 

 

 

 

 

(p − λ)(p − μ)

2

 

 

 

 

p4 − ω4

 

 

 

 

 

 

λ − μ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

27

1

 

(1 − cos ωt)

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

36

1

 

(ch ωt + cos ωt)

 

 

 

 

 

p3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω2

 

 

 

 

p(p2 + ω2)

 

 

 

2

 

 

 

 

 

p4 − ω4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60

полюсы

Нахождение оригинала по изображению

Если функция F (p) аналитична в некоторой окрестности бесконечно удал¼нной точки (в кольце R < |p| < ∞) и е¼ разложе-

1

ние в ряд Лорана по степеням p имеет вид

F (p) =

то функция

X

tn

f(t) = cn n!

n=0

X cn ,

n=0 pn+1

 

tn−1

X

 

=

cn

 

n=1

 

(n 1)!

 

 

является оригиналом, соответствующим изображению F (p). Если изображение F (p) является однозначной функцией и име-

ет лишь конечное число особых точек p1, p2, . . . , pn, лежащих в конечной части плоскости, в частности, если изображение F (p)

является дробно-рациональной функцией, прич¼м степень числителя меньше степени знаменателя и p1, p2, . . . , pn

этой функции, то его оригиналом служит функция

 

n

 

 

f(t) =

res(eptF (p)).

 

k=1

pk

 

X

 

 

Pm(p)

В случае F (p) = R(p) =

 

 

, ãäå Pm(p), Qn(p) многочлены

 

 

Qn(p)

степени m и n соответственно, не имеющие общих корней. Если

все полюсы p1, p2, . . . , pn функции F (p) простые (т. е. Π(1)), то справедлива формула

 

Pm(p)

Pm(pk)

res

 

ept =

 

 

epkt,

 

Qn

 

pk

Qn(p)

(pk)

n

f(t) = X Pm(pk)epkt.

k=1 Qn(pk)

61

Применение таблицы и свойств преобразования Лапласа для нахождения оригинала по изображению

1.Если изображение отличается от табличного на постоянный множитель, то следует умножить и одновременно поделить на этот множитель, а затем воспользоваться свойством линейности.

2.

Изображение, заданное в виде дроби

a ±

b

, разлагается на

c

 

 

сумму дробей.

 

 

 

3.

Если знаменатель дроби содержит квадратный тр¼хчлен, то

 

в н¼м выделяется полный квадрат: ap2 + bp + c =

= a(p ± α)2 ± ω2. При этом числитель дроби представляется

в виде многочлена от p ± α.

4.Если изображение представляет собой правильную рациональную дробь, следует разложить е¼ на простейшие дроби и для каждой из них найти оригинал. Также можно воспользоваться второй теоремой разложения и теоремой Бореля (теоремой о св¼ртке).

Алгоритм решения задачи Коши

1.Применить преобразование Лапласа: от неизвестных и известных функций перейти к их изображениям, записать уравнение (систему) в изображениях, соответствующее решаемой задаче.

2.Решить полученное уравнение (систему): найти изображение искомого решения.

3.Применить обратное преобразование Лапласа: найти оригинал для полученного в п. 2 изображения.

62

Алгоритм решения задачи Коши с помощью единичной переходной функции

1.Найти единичную переходную функцию h(t), соответствующую исходной задаче.

2.Найти производную h(t) единичной переходной функции.

3.Определить решение исходной задачи Коши по формуле

Zt

x(t) = h(τ)f(t − τ)dτ.

0

Алгоритм решения задачи Коши с помощью импульсной переходной функции

1.Найти импульсную переходную функцию k(t), решая вспомогательную задачу, соответствующую исходной.

2.Найти решение исходной задачи по формуле

Zt Zt

x(t) = k(τ)f(t − τ)dτ = k(t − τ)f(τ)dτ.

0 0

Связь между переходными функциями k(t) = h(t).

63

Ряды Фурье. Интеграл Фурье

Скалярное произведение функций в пространстве L2[a, b] определяется как

Zb

(f, g) = f(t) g(t)dt ( f, g L2[a, b])

a

или в более общем случае как

Zb

(f, g) = f(t) g(t) ρ(t)dt ( f, g L2[a, b]),

a

где ρ(t) > 0, t [a, b] некоторая определ¼нная функция, называемая весом.

Примеры норм

1)

Rn : x = s

k=1 ξk2

ïðè n = 1: x = |x|;

 

 

n

 

 

 

 

 

 

P

 

 

 

 

2)

C[0; 1] : x = maxt

 

|x(t)|;

 

 

Lp[0; 1] : x =

R

|x(t)|pdt

1

 

1

p

3)

0

.

Система функций {φk}k N некоторого функционального пространства F называется базисом, при выполнении условий:

1)любая конечная подсистема {φk}, k = 1, 2, . . . , n, образует линейно независимую систему;

2)система {φk}k N полна в F , т. е.

 

n

f F {ak}k N

X

: nlim ρ(f, akφk) = 0.

 

→∞

 

k=0

64

Система функций {φk}k N некоторого функционального пространства F называется ортогональной, в случае выполнения:

1)φk ̸= 0, k N;

2)i, φj) = 0 ïðè i ̸= j.

Если дополнительно φk = 1 для k N, то система называется ортонормированной.

1.Множество функций {sin nx}k N образует ортогональную на отрезке [0, π] систему функций.

2.Множество функций {1, cos nx, sin nx}k N образует ортого- нальную на отрезке [−π, π] систему функций.

Имеем:

π

 

 

 

 

 

cos nxdx =

2π, n̸= 0

 

Z

 

 

 

 

0, n = 0

 

−π

 

 

 

 

 

π

 

0, n = 0

 

Z

 

 

 

sin nxdx =

 

0, n ̸= 0

 

−π

 

 

 

 

Zπ

 

 

 

 

cos mx cos nxdx =

0, m ̸= n

 

 

 

 

 

π, m = n

−π

Zπ

sin mx sin nxdx =

0, m ̸= n π, m = n

−π

Zπ

sin mx cos nxdx = 0.

−π

3.Если совокупность степеней 1, t, t2, . . . , tn, . . . ортогонализовать в пространстве L2[a, b] с весом ρ(t), то прид¼м к системе многочленов: p0(t) = const, p1(t), p2(t), . . . , pn(t), . . .

65

При ρ(t) ≡ 1, a = −1, b = 1 получим с точностью до постоянных множителей систему многочленов Лежандра:

L0(x) = 1, Ln(x) =

1

 

n

2

1)

n

n N =

 

d

(x −

 

2nn!

 

dxn

 

 

 

 

 

 

 

 

 

 

=1; x; 32x2 12; 52x3 32x; 358 x4 154 x2 + 38; . . . ,

 

Lk+1

=

(2k + 1)xLk(x) − kLk−1(x)

.

 

 

 

 

 

 

 

 

 

 

 

 

k + 1

 

 

 

 

 

 

 

 

4. Полиномы Чебышева

 

 

 

 

 

 

dn

 

 

 

 

 

 

 

 

 

 

 

 

(−2)nn!

 

 

 

 

 

 

 

2n−1

 

 

 

 

 

 

 

 

1 − x2

T

(x) =

 

 

1

x2

(2n)!

 

 

 

 

 

 

 

 

n

 

 

p

 

 

 

dxn

 

 

 

 

1

 

образуют ортогональную систему с весом

ρ(x) =

 

1 − x2

в пространстве L2[−1; 1]:

T0(x) = 1,

T

(x) = x, T

(x) = 2x2

1,

T3(x) = 4x

3

− 3x,

4

1

 

2

2

 

 

 

T4(x) = 8x

− 8x

 

+ 1, . . . , Tk+1 = 2xTk(x) − Tk−1(x) .

Обобщ¼нным рядом Фурье функции f(x) L2[a, b] по ортого-

нальной системе функций {φk(x)}k N называется функциональный ряд

X

f(x) = c1φ1(x) + c2φ2(x) + . . . + cnφn(x) + . . . = ckφk(x),

k=1

коэффициенты которого находятся по формуле

ck = (f, φk2).

φk

66

Тригонометрическим рядом Фурье периодической функции f(x) с периодом 2π, определ¼нной на промежутке [−π, π], на-

зывается ряд

 

 

 

 

 

 

a0

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

+ (an cos nx + bn sin nx),

 

 

 

 

 

 

2

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

ãäå an =

1

 

Z

f(x) cos nxdx, (n = 0, 1, 2, . . .),

 

 

 

 

π

 

 

 

 

 

−π

 

 

 

 

 

 

 

 

π

 

 

 

bn =

1

 

 

Z

f(x) sin nxdx, (n = 1, 2, . . .),

 

 

π

−π

an, bn называются коэффициентами ряда Фурье.

Если этот ряд сходится, то его сумма S(x) есть периодическая функция с периодом 2π, т. е. S(x + 2π) = S(x).

Если функция f(x) задана на сегменте [−L, L], где L про-

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

 

a

 

 

 

 

 

an cos

 

nπx

 

nπx

 

20 + n=1

 

L

 

+ bn sin L ,

ãäå

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

nπx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an =

 

 

 

Z f(x) cos

 

 

 

dx, (n = 0, 1, 2, . . .),

L

 

 

L

 

 

 

 

 

 

−L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

bn =

1

Z

f(x) sin

nπx

dx, (n = 1, 2, . . .).

 

 

 

 

 

 

 

 

L

 

 

L

−L

67

и косинусы, т. е.

В случае когда f(x) четная функция, е¼ ряд Фурье содержит только свободный член a0

 

 

 

 

a0 +

 

 

 

 

 

 

 

 

an cos nπx,

 

 

 

 

 

 

X

 

 

 

 

 

 

 

2

n=1

 

L

 

 

 

 

 

 

 

 

 

ãäå

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

an =

2

Z0

f(x) cos

nπx

dx, (n = 0, 1, 2, . . .).

 

 

 

 

L

 

L

В случае когда f(x) нечетная функция, е¼ ряд Фурье содержит только синусы, т. е.

 

 

 

 

 

 

 

 

 

 

bn sin nπx,

 

 

 

X

 

 

 

 

 

 

 

n=1

 

 

L

 

 

 

 

 

 

 

ãäå

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

bn =

2

Z0

f(x) sin

nπx

dx, (n = 1, 2, . . .).

 

 

L

L

В случае когда f(x) задана на произвольном интервале (a, b),

b − a

обозначим L = 2 , зададим конкретное значение функции на одном из концов интервала (например, при x = b) и продолжим

данную функцию периодически, с периодом T = 2L на всю чис-

ловую ось. Полученная функция будет удовлетворять условиям теоремы Дирихле.

Коэффициенты ряда Фурье в этом случае можно определить следующим образом:

 

 

 

 

b

 

 

 

 

 

1

Za

nπx

 

 

an =

 

 

f(x) cos

 

 

dx, (n = 0, 1, 2, . . .),

L

 

L

 

 

 

 

 

 

b

 

 

 

 

 

1

 

Za f(x) sin

nπx

bn =

 

 

 

 

 

dx, (n = 1, 2, . . .).

 

L

L

 

68

Тогда ряд Фурье принимает вид

a

an cos

nπx

20 + n=1

L

 

X

 

+ bn sin

nπx

.

L

Если рассматривать f(t) как T -периодическую функцию с пе-

риодом T = 2L, ввести обозначение ω = T , учесть, что для T -периодической функции справедливо равенство

ZT c+ZT

f(t)dt = f(t)dt, c R,

0

c

то можно получить следующую форму ряда Фурье:

 

 

 

 

 

a0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

f(t)

2

+

(an cos nωt + bn sin nωt) =

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

a0

 

 

 

 

 

 

a0

 

 

 

 

 

 

X

 

 

 

 

 

 

 

X

 

 

 

=

2

 

+

 

An cos(nωt − φn) =

2

+

An sin(nωt + ψn),

 

 

 

 

 

n=1

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bn

 

ψ = arctgpn сдвиги по фазе соответствующих гармоник.

,

 

ãäå An =

 

an2

+ bn2 амплитуда,

nω частота, φn

= arctg

an

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n bn

Комплексная форма ряда Фурье

f(x) =

cneinx, cn =

1

Zπ

f(x)e−inxdx.

X

 

 

 

π

 

 

−∞

 

 

 

 

Если функция f(x)

задана на сегменте [−L, L], где L про-

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

 

1

L

 

 

 

 

 

X

 

 

 

 

f(x) =

cn e−inπx/L, cn =

2L

Z

f(x)e−inπx/Ldx.

−∞

 

 

 

L

 

69

Правая часть тождества Бесселя имеет вид

 

 

n

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

a02

 

 

T

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

T − 2

 

 

 

f(t) 2 k=1 ck φk(t) 2=

Z f2(t)dt − 4

k=1 Ak.

 

Относительная среднеквадратическая ошибка

 

n2 =

 

n2

может

 

 

T

быть записана через амплитуду отдельных гармоник:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

T

 

 

 

a02

 

 

 

 

1 n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T Z f2(t)dt

4 2 k=1 Ak.

 

 

 

 

 

 

Равенство Парсеваля имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

a02

 

 

 

 

T

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z f2(t)dt − 4 T = 2 k=1 Ak.

 

 

 

 

 

 

 

 

 

 

 

Энергия отдельной гармоники Ak cos(kωt − φk) за время T :

 

T

 

 

 

 

 

A2

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ek = Z0

 

 

 

 

 

 

Z0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

Ak2 cos2(kωt − φk)dt =

 

k

 

(1 + cos(2(kωt − φk))dt = Ak2

 

,

2

 

 

2

тогда энергия спектра может быть найдена по формуле

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E =

2

 

 

 

 

Ak2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя равенство Парсеваля, энергию

спектра

запишем

â âèäå:

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E = Z0 f2(t)dt −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

T.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

Оценить долю k-й гармоники в общей энергии спектра можно

по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

Ek

=

Ek

=

Ak2 T2

 

 

=

 

 

Ak2

 

 

.

 

T

a2

2

T

a2

 

 

E

R0 f2(t)dt −

0

T

 

 

 

R0 f2(t)dt −

0

 

 

 

 

4

 

 

T

2

70

Двойной интеграл Фурье

 

∞ ∞

f(t) cos ω(t − x)dt dω,

f(x) = π Z Z

1

 

 

0−∞

или интеграл Фурье для функции f(x):

Z

f(x) = (a(ω) cos ωx + b(ω) sin ωx) dω,

0

ãäå

 

a(ω) =

1

Z

f(t) cos ωtdt, b(ω) =

1

 

Z

f(t) sin ωtdt.

 

 

 

π

π

 

 

−∞

 

 

 

−∞

 

Комплексная форма интеграла Фурье

 

 

 

f(x) =

1

Z

Z

f(t)e−iω(t−x)dt.

 

 

−∞

−∞

 

Интеграл Фурье в комплексной форме можно записать в дру-

ãîì âèäå:

 

 

Z

 

 

 

 

 

 

 

 

f(t)e−iωtdt.

 

F (ω) =

 

 

−∞

 

 

 

 

 

 

f(x) =

1

Z

F (ω)eiωxdω.

 

−∞

Определение 1. Соответствие f(x) → F (ω) называется

преобразованием Фурье функции f(x).

Преобразование F (ω) → f(x) называется обратным преобразованием Фурье.

71

Основные теоремы преобразования Фурье

1.Теорема сложения

Åñëè f1(t) ↔ F1(ω) è f2(t) ↔ F2(ω),

òî α1f1(t) ± α2f2(t) ↔ α1F1(ω) ± α2F2(ω).

2.Теорема симметрии

Åñëè f(t) ↔ F (ω), òî F (±t) ↔ 2πf( ω).

3.Теорема подобия

Если f(t) ↔ F (ω) и α вещественная постоянная,

òî f(αt) ↔ 1F (ω/α).

α

4.Теорема о временном сдвиге

Åñëè f(t) ↔ F (ω), òî f(t ± α) ↔ e±iαωF (ω).

5.Теорема о частотном сдвиге

Åñëè f(t) ↔ F (ω), òî F (ω ± ω0) ↔ e 0tf(t).

6.Теоремы модуляции

1

f(t) cos(ω0t) ↔ 2(F (ω + ω0) + F (ω − ω0)),

i

f(t) sin(ω0t) ↔ 2(F (ω + ω0) + F (ω − ω0)).

7. Теорема о производной

Åñëè f(t)

F (ω), òî f(t)

iω F (ω) è F (ω)

it f(t),

f

(n)

 

n

F (ω) è F

(n)

 

 

n

f(t).

 

 

(t) ↔ (iω)

 

(ω) ↔ (−it)

 

 

 

8.

Теорема об интеграле

1

 

 

 

 

 

 

 

 

 

 

Åñëè f(t) ↔ F (ω), òî

f(t)dt ↔

 

F (ω).

 

 

 

 

−∞

 

 

 

 

 

 

 

 

R

 

 

 

 

 

9.

Теорема о св¼ртке

 

 

 

 

 

 

 

 

 

 

 

 

f1(t) f2(t) =

Z

f1(τ)f2(t − τ)dτ =

Z

f1(t − τ)f2(τ)dτ,

 

 

−∞

 

 

 

 

−∞

 

Åñëè f1(t) ↔ F1(ω) è f2(t) ↔ F2(ω),

òî f1(t) f2(t) ↔ F1(ω) F2(ω) è f1(t) f2(t) ↔ F1(ω) F2(ω), (f1(t) f2(t))= f1(t) f2(t) = f1(t) f2(t) ↔ iω F1(ω) F2(ω).

72

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