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

vm3

.pdf
Скачиваний:
24
Добавлен:
11.05.2015
Размер:
1.32 Mб
Скачать
= cos 8 + 165 sin 8 −

функции f (z), а потому все особые точки лежат внутри γ. Теперь

(4.13) следует из (4.12).

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

4.3. Приложение вычетов к вычислению интегралов

4.3.1. Вычисление интегралов по замкнутому контуру

Применяя основную теорему о вычетах, можно вычислять многие контурные интегралы.

 

 

1

I

sin 4zdz

Пример 4.12. Вычислить I =

 

 

.

2πi

(z − 2)2(z − 3)(z − 6)

 

 

 

 

|z|=5

 

 

Решение. По основной теореме о вычетах

I = 2πi

I

(z − 2)2(z − 3)(z − 6) = Res (z − 2)2(z − 3)(z − 6) ; z = 2 +

1

 

sin 4zdz

 

 

sin 4zdz

|z|=5

+Res (z − 2)2(z − 3)(z − 6) ; z = 3

 

sin 4zdz

+ lim

sin 4z

(z − 2)2(z − 6)

z→3

= lim

d

 

sin 4z

+

 

 

z→2 dz (z − 3)(z − 6)

 

sin 12 . 3

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

 

 

 

 

 

1

 

I

 

z15dz

 

 

 

 

 

 

Пример 4.13. Вычислить I =

 

 

 

 

 

.

 

 

 

 

 

 

 

2πi

z8 + 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|z|=3

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Видим, что I = −Res

 

z15dz

 

; z = ∞ . Так как

 

 

z8 + 2

 

 

 

 

z15

1

= z7 1 −

 

2

4

− · · ·

 

 

 

 

 

 

= z7

 

 

+

 

 

 

 

z8 + 2

1 + 2/z8

z8

z16

 

 

при |z| > 8

 

, то I = −2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

вычета в

 

, то

Если в этом примере не использовать понятие

 

8

 

 

 

пришлось бы вычислять восемь вычетов в точках zm =

 

.

 

 

2

 

 

101

 

 

 

 

 

 

 

 

 

2π

 

 

 

4.3.2. Вычисление интегралов типа I =

R(cos t, sin t)dt,

 

 

 

 

 

 

 

 

 

0

 

 

 

 

где R(u, v)

рациональная функция

своих аргументов

 

R

z + z ,

Пусть

z =

eit. Тогда sin t = 2i

z − z ,

cos t = 2

 

 

 

 

1

 

1

 

 

 

1

 

1

 

dz

. Когда t изменяется от 0 до , переменная z пробегает

dt = −i

 

z

окружность |z| = 1 в положительном направлении. Следовательно,

I =

2π

R(cos t, sin t)dt = −i

I

R z −2i

,

2

z =

0

 

R

 

 

 

1/z

z + 1/z

dz

I|z|=1

=R1(z)dz.

|z|=1

Через R1(z) обозначена некоторая рациональная функция от z. По

n

X

основной теореме о вычетах I = 2πi Res[R1(z); z = zm], где

m=1

zm полюсы функции R1(z), попавшие внутрь окружности |z| = 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

π

dt

 

 

 

 

 

 

Пример 4.14. Вычислить интеграл

I =

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

5 + 3 cos t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

Решение.

 

 

dz

 

 

 

 

 

 

 

R

 

 

dz

 

 

I = −i

I

 

 

 

 

 

= −i

I

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

[5 + (3/2)(z + 1/z)]z

5z + (3/2)z2 + 3/2)

 

|z|=1

 

dz

 

2i

 

|z|=1

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −2i

 

I

 

 

= −

 

I

 

 

 

 

 

=

 

 

 

 

3z2 + 10z + 3

3

(z + 3)(z + 1/3)

 

 

|z|=1

 

 

 

 

 

|z|=1

 

 

 

 

 

= 3(3 − 1/3) =

2 .

= 2πi

3

 

Res (z + 3)(z + 1/3) ; z = −3

 

 

2i

 

 

dz

 

 

 

1

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

+∞

 

 

 

 

 

 

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

f (x)dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

Теорема 4.9. Пусть функция f (z)

аналитичнаR

всюду в верхней

полуплоскости, кроме конечного числа особых точек, и аналитична

на вещественной оси. Если при больших

z

|

и на вещественной оси

 

 

 

 

c

|

 

выполняется неравенство |f (z)| <

 

 

, δ > 0, c = const, то инте-

 

z 1+δ

+∞

f (x)dx сходится и

+∞

| |

 

 

 

Res[f (z); z = zm], где

грал

f (x)dx = 2πi

 

−∞

 

−∞

 

 

 

 

 

m

 

 

 

 

 

X

R

 

R

 

 

 

 

 

 

zm все особые точки f (z), расположенные в верхней полуплоскости, т.е. Im zm > 0.

102

Доказательство. Сходимость интеграла

+∞

 

 

 

 

f (x)dx, и притом аб-

 

 

 

 

−∞

 

 

 

 

солютная,

следуетR

из неравенства

 

 

c

 

+∞

 

dx

 

 

 

 

 

R

 

 

 

|

|

 

 

|

|

 

|f (x)| < x 1+δ

и сходимости −∞

 

x 1+δ .

Возьмём окружность столь большого радиуса R, чтобы все особые точки f (z) попали внутрь неё. Составим контур L из верхней полуокружностии части вещественной оси (−R, R).

Тогда

 

 

n

R

 

 

 

 

 

 

 

I

f (z)dz = 2πi m=1Res[f (z); z = zm] = Z

f (x)dx+Z

f (z)dz. (4.14)

L

 

X

R

 

 

Докажем, что

→∞ R

f (z)dz = 0. Имеем

 

 

lim

 

 

R

 

f (z)dz <

 

|f (z)||dz| < Z

cds

=

 

cπR

=

πc

→ 0

 

 

 

 

 

 

 

z

1+δ

R1+δ

Rδ

R

→ ∞

 

 

R

 

| |

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

при R

 

так как δ > 0. Справедливость теоремы теперь следует

из соотношения (4.14) при R → ∞.

 

+∞

 

x2dx

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4.15. Вычислить интеграл

R

 

 

z2

 

.

−∞

(x2 + 1)(x2 + 9)

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Ясно, что для функции f (z) = (z2 + 1)(z2 + 9) все усло-

вия предыдущей теоремы выполнены. Поэтому

 

 

 

 

 

I = 2πi {Res[f (z); z = i] + Res[f (z); z = 3i]} =

16

 

4

 

 

2i(−1 + 9)

(−9 + 1)6i

16

 

 

= 2πi

−1

9

 

= 2π

1

+

3

 

=

π

.

 

 

 

 

 

 

 

4.3.4. Вычисление несобственных интегралов

 

 

 

+∞

 

 

 

 

 

 

 

 

 

 

типа

eiαxf (x)dx, α > 0

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

−∞

Часто используется следующая лемма.

Лемма (Жордана). Пусть дано некоторое семейство дуг окружностей CR : |z| = R, Im z > −a (R → ∞, a > 0 фиксировано).

Обозначим M (R) = sup|f (z)|, z CR. Если

(4.15)

 

 

lim M (R) = 0,

 

 

R→∞

 

то lim

eiαz f (z)dz = 0, где α действительно, α > 0.

 

R

→∞

R

 

 

 

CR

 

103

 

Доказательство.

Ограничим-

 

ся

случаем a = 0.

Пусть z =

 

= x + iy = re. Тогда на дуге CR

 

 

 

eiαz

= eαR sin ϕ.

(4.16)

 

 

 

 

 

 

 

 

 

 

Поскольку функция sin ϕ выпук-

 

ла, то её график

лежит выше хор-

 

ды, соединяющей точки

(0, 0) и

π

 

π

, 1 ,

 

2

 

 

 

2

 

 

 

 

 

 

 

 

поэтому

 

 

 

при 0 ≤ ϕ ≤ 2 справедливо неравенство sin ϕ ≥ π ϕ. Теперь из (4.16)

 

 

 

 

iαz

 

 

 

 

 

 

 

2

 

следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем

e

 

 

≤ exp −

π αRϕ ,

 

 

 

 

π/2

 

 

AB e

iαz

 

 

 

 

 

 

 

 

iαz

||f (z)|ds ≤ RM (R) Z

exp −

2

αRϕ dϕ =

 

f (z)dz

AB |e

 

 

π

R

 

 

 

 

 

R

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

π

 

 

αR

 

 

 

 

 

 

 

 

 

 

 

 

·

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= M (R)

 

 

1

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.15) следует, что

lim

eiαz f (z)dz = 0.

Отсюда и из условия

 

 

 

 

 

 

 

→∞

R

 

 

 

 

 

 

 

 

 

 

 

 

R AB

Если на дуге CB отсчитывать полярный угол от отрицательной оси по часовой стрелке, то для интеграла R eiαz f (z)dz получается

CB

такая же оценка, и утверждение леммы в случае a = 0 доказано.

В операционном исчислении применяется лемма Жордана в несколько видоизменённой форме.

Пусть дано некоторое семейство дуг окружностей CR : |p| = = R, Re p < a. Для любой функции F (p) такой, что величина

M (R) = sup |F (p)| → 0, p CR при R → ∞, и любого вещественного положительного t

R→∞ Z

 

lim

F (p)eptdp = 0.

(4.17)

CR

Справедливость этого утверждения следует из доказанной леммы, если положить в ней iz = p.

Теорема 4.10. Пусть функция f (z) аналитична в верхней полуплоскости, кроме конечного числа особых точек zm, аналитична на

вещественной оси и удовлетворяет условиям леммы Жордана. Тогда

+∞

интеграл R eiαxf (x)dx, α > 0 сходится и

−∞

X

I = 2πi Res eiαz f (z); z = zm , Im zm > 0.

m

104

Доказательство. Пусть L замкнутый контур, образованный верхней частью CR окружности |z| = R и интервалом (−R, R). Величину R возьмём столь большой, чтобы все особые точки zm по-

пали внутрь L. Тогда

L

eiαz f (z)dz = 2πi

m

Res

eiαz f (z); z = zm =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

R

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

eiαz f (z)dz +

 

eiαxf (x)dx, Im zm > 0. Устремим R → ∞. Пер-

 

CR

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

выйR

интеграл по леммеR

Жордана стремится к нулю, а второй к

искомому интегралу. Теорема доказана.

 

 

 

 

 

 

 

 

 

 

 

Замечание 1. Если α < 0, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Im zm < 0.

 

 

 

−∞ eiαxf (x)dx = −2πi X Res eiαz f (z); z = zm ,

 

 

 

R

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание 2. Если значения функции f (z) при вещественных z

вещественны, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+∞f (x) sin αxdx = Im "2πi

Res

eiαz f (z); z = zm # ,

(4.18)

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+∞f (x) cos αxdx = Re "2πi

Res

eiαz f (z); z = zm # ,

(4.19)

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

где Im zm > 0.

 

 

 

 

 

 

 

 

 

+∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin xdx

 

 

 

 

 

 

 

Пример 4.16. Вычислить I = Z

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 + 2x + 10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Сходимость данного интеграла следует из признака

+

 

 

 

 

 

 

 

 

 

 

 

 

 

sin x

 

 

 

 

 

 

1

 

 

 

сравнения, поскольку

x2

+ 2x + 10

x2

+ 2x + 10 , а интеграл

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 + 2x + 10

сходится. Условия леммы Жордана для функции

−∞

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (z) =

 

 

 

 

, очевидно, выполнены. По формуле (4.18) нахо-

 

 

z2 + 2z + 10

дим I = Im 2πiRes

 

 

 

 

eiz

 

 

 

 

 

 

 

; z = −1 + 3i =

 

(z + 1 + 3i)(z + 1 − 3i)

 

 

= Im 2πi lim

 

 

 

eiz

 

=

πe−3

Im ei =

πe−3

sin 1.

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

z

→−

1+3i

z + 1 + 3i

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

105

5. Интегралы, зависящие от параметра.- и B- функции. Функции Бесселя

m

X

Слагаемые в сумме f (n, y) представляют собой функции двух

n=1

аргументов: аргумента n, принимающего целые положительные значения, по которому производится суммирование, и аргумента y, выступающего как параметр. В результате суммирования по n получаем некоторую функцию ψ(y) параметра y. В интегралах, зависящих

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

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

также в данном разделе.

5.1. Свойства функций, заданных собственными интегралами, зависящими от параметра

Пусть в прямоугольнике П:{a ≤ x ≤ b, c ≤ y ≤ d} задана функция z = f (x, y), интегрируемая по x на [a, b] при любом значении y из [c, d]. Тогда на [c, d] определена функция

I(y) = Z

b

 

f (x, y)dx.

(5.1)

a

Говорят, что функция I(y) задана интегралом, зависящим от параметра. В тех случаях, когда интеграл (5.1) удаётся вычислить, получим явно заданную функцию I(y), и никаких новых задач не возникает. Мы будем изучать свойства функции I(y), исходя из свойств подынтегральной функции f (x, y), считая явный вид функции I(y)

неизвестным. Заметим, что интеграл (5.1) представляет собой новый способ задания функциональных зависимостей, широко применяемый в различных областях математики и её приложениях.

5.1.1. Непрерывность и переход к пределу под знаком интеграла

Теорема 5.1. Если функция f (x, y) непрерывна в прямоугольнике П, то функция I(y) непрерывна на отрезке [c, d].

106

Доказательство. Так как функция f (x, y) непрерывна в ограни-

ченном замкнутом прямоугольнике П, то она равномерно непрерывна в нём. Это означает, что для любого ε > 0 найдётся δ = δ(ε) > 0, не зависящее от выбора точек (x1, x2) и (y1, y2) из П, такое, что если

|x1 − x2| < δ, |y1 − y2| < δ, то

ε

 

|f (x1, y1) − f (x2, y2)| <

(5.2)

 

.

b − a

Пусть x любое из [a, b], а y любое из [c, d]. Дадим значению y приращение y такое, чтобы точка y + y не вышла из [c, d]. Поло-

жим в (5.2) x1 = x2 = x, y1 = y, y2 = y + y. Тогда

 

 

 

 

|I(y1) − I(y2)| = |I(y) − I(y + y)| = a [f (x1, y1) − f (x2, y2)]dx

b

 

 

b

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε

 

≤ |

|

|

 

|

 

 

 

 

 

b

a

R

f (x, y +

ε

f (x, y)

 

f (x, y +

y)

<

 

 

f (x, y)

y) dx. Так как

 

 

 

 

a

в силу (5.2), то |I(y)−I(y+ y)| < b − a ·(b−a) = ε, если | y| < δ. Это и означает непрерывность функции I(y) на [c, d]. Теорема доказана.

Следствие. Если функция f (x, y) непрерывна в прямоугольнике

П, то

b

 

b

 

 

yy0

Za

 

 

Za

yy0

(5.3)

lim

f (x, y)dx =

 

lim f (x, y)

dx,

где y0 любое из [c, d].

Действительно, в силу доказанной в теореме непрерывности I(y)

b

R

lim I(y) = I(y0) =

f (x, y0)dx. Так как функция f (x, y) непрерывна

y y0

a

 

в П, то f (x, y0) = lim f (x, y), и формула (5.3) получена.

yy0

Если выполняется соотношение (5.3), то говорят, что возможен предельный переход под знаком интеграла.

1

Пример 5.1. Найти lim R x2ex2y2 dx.

y→0 0

Решение. Функция f (x) = x2ex2y2 непрерывна в прямоугольнике П:{0 ≤ x ≤ 1, c ≤ y ≤ d}, где c < 0, d > 0. Поэтому справедливо соот-

 

1

 

2 x2y2

 

1

lim x2ex2y2

dx =

1

ношение (5.3), т.е.

lim

x

dx =

 

x2dx =

R

e

0

y→0

 

0

= 1/3.

 

 

 

R

 

 

R

 

y→0 0

 

 

 

 

 

 

 

Заметим, что условие непрерывности функции I(y) в теореме 5.1, является достаточным, но не необходимым, т.е. функция I(y) может

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

107

5.1.2. Дифференцируемость. Формула Лейбница

Теорема 5.2 (о дифференцировании по параметру под знаком интеграла). Если функции f (x, y) и fy(x, y) непрерывны в прямоуголь-

b

нике П:{a ≤ x ≤ b, c ≤ y ≤ d}, то функция I(y) = f (x, y)dx диф-

a

ференцируема на [c, d] и

при этом справедлива

формула Лейбница

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

dI

 

d

 

 

b

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

Za f (x, y)dx = Za

fy(x, y)dx.

 

 

(5.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

dy

 

 

 

Доказательство. Надо показать, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I(y + y) − I(y)

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

I

(y) =

lim

 

=

f

(x, y)dx.

 

 

 

 

 

 

 

 

 

 

y→0

 

 

 

 

y

 

Za

y

 

 

 

 

 

 

 

 

 

 

 

 

I(y + y) − I(y)

 

b

f (x, y + y) − f (x, y)

 

 

 

 

 

Находим

 

 

=

dx. По фор-

 

 

 

 

 

 

 

y

 

 

 

 

 

Za

 

 

 

y

 

 

f (x, y) = f (x, y+

муле Лагранжа можно записать f (x, y +

 

y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

y)

y, где

0

Θ

 

1. Следовательно,

I(y + y) − I(y)

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

b

f (x, y + Θ

 

 

 

 

 

 

 

I(y)

 

 

 

b

f

 

 

 

 

 

=

 

y)dx,

т.е.

= lim

 

(x, y + Θ

 

y)dx. Так

a

 

y

 

 

 

 

 

 

 

 

 

 

y→0 a

y

 

 

 

 

 

 

какR

функция f

(x, y + Θ

 

y) непрерывна

Rв П, то можно перейти

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к пределу под знаком интеграла согласно следствию из теоремы

5.1. Получаем I(y) =

R

b lim fy(x, y + Θ y) dx, следовательно,

 

a y→0

b

I(y) = R fy(x, y)dx. Теорема доказана.

a

Если имеет место формула (5.4), то говорят, что возможно дифференцирование по параметру под знаком интеграла. Формула (5.4) позволяет упростить вычисление некоторых интегралов.

Пример 5.2. Применяя дифференцирование по параметру, вычис-

 

1

1

1

 

 

R

 

лить интеграл I(m) = xm ln3 xdx, m > 0.

 

0

 

 

 

 

Решение. Имеем I(m) =

 

xmdx =

 

. Дифференцируя три-

0

m + 1

жды этот интеграл по

параметру m, на основании теоремы 5.2 по-

 

R

 

 

 

108

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

лучаем

R3

xm ln xdx = −

 

 

 

 

 

,

xm(ln x)2dx =

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(m + 1)2

(m + 1)3

 

 

 

 

 

 

1

m

 

 

6

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

(ln x) dx = −

(m + 1)4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y cos x)dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 5.3. Найти I(y) = Z0

 

ln(1 +

 

 

, |y| < 1.

 

 

 

 

 

 

 

 

 

cos x

 

 

 

 

 

 

 

 

 

 

 

Решение. Подынтегральная

функция

 

f (x, y) =

ln(1 + y cos x)

 

 

 

 

cos x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

устранимый

разрыв

 

 

в

точках

, y

,

так

 

 

 

как

 

 

 

 

 

 

 

 

имеет

 

 

 

ln(1 + y cos x)

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

lim

 

 

 

 

 

 

 

 

= y0. Этот разрыв можно устранить, по-

 

 

 

 

 

cos x

 

 

yy0,xπ/2

 

 

 

 

 

π

 

 

 

 

fy(x, y0) в

1

 

2

 

 

0

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

лагая f

 

π

, y

= y. Значение производной

 

 

 

 

 

 

точке

 

π

, y

 

 

 

 

 

 

 

 

 

 

 

 

 

определим

 

из

условия

 

fy

 

2 , y0

 

= yy0

,xπ/2 1 + y cos x

=

 

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

 

fy(x, y)

 

 

 

 

 

 

 

После

такого

доопределения

функции f (x, y)

и

 

будут

непрерывными в прямоугольнике {0 ≤ x ≤ π, −1 + ε ≤ y ≤ 1 − ε},

π

0 < ε < 1. Применяя формулу (5.4), получим I(y) = Z0

dx

 

.

1 + y cos x

Для вычисления интеграла I(y) применим теорию вычетов (см.

4.3.2). Если y = 0, то I(0) = π. Пусть y = 0. Положив z = eix,

 

 

 

 

 

1

π

 

 

dx

 

 

 

 

 

 

1

 

I

 

 

 

6

 

−idz

 

 

 

 

 

 

получим I

(y) =

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

2 Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + y cos x

 

 

 

2

 

 

 

z

 

1 + y

 

1

 

z +

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

·

2 ·

z

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

|z|=1

 

 

 

 

 

 

 

 

 

= 2πRes

" yz2 + 2z + y ; z = −y + p

 

y

y2

#

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

z=− y +

 

 

 

 

 

 

 

 

 

π

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2yz + 2

y y

2 =

 

 

 

1 y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1−

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I(y) = π Z

 

 

 

dy

 

 

 

Отсюда, интегрируя

 

по

y,

 

 

находим

 

 

 

 

+

 

 

 

 

 

 

 

1 − y2

+c

= π arcsin y + c. Так как I(0)

 

 

=

 

0, то

c

= 0,

 

p

 

 

 

 

 

 

 

 

 

 

 

 

следовательно,

I(y) = π arcsin y.

Рассмотрим случай, когда пределы интегрирования также зависят от параметра y. Пусть функция f (x, y) определена на замкну-

109

том множестве D : {a(y) ≤ x ≤ b(y), c ≤ y ≤ d}, причём кривые x = a(y) и x = b(y) не выходят за пределы прямоугольника П:{a ≤ x ≤ b, c ≤ y ≤ d}, и интегрируема по x на отрезке [a(y), b(y)] при любом фиксированном значении y из [c, d]. Тогда на сегменте

 

b(y)

 

[c, d] определена функция I(y) =

R

f (x, y)dx. Функция I(y) пред-

a(y)

ставлена интегралом, зависящим от параметра, пределы интегрирования которого также зависят от параметра y.

Теорема 5.3. Если функции f (x, y) и fy(x, y) непрерывны в прямоугольнике П, а функции a(y) и b(y) дифференцируемы на [c, d] и удовлетворяют условию a ≤ a(y) ≤ b, a ≤ b(y) ≤ b при c ≤ y ≤ d, то

 

b(y)

 

 

 

 

 

 

 

a R

f (x, y)dx дифференцируема на [c, d] и при этом

функция I(y) =

 

 

(y)

 

 

 

 

 

 

 

 

b(y)

 

 

 

 

 

 

I(y) =

Z

fy

(x, y)dx + f [b(y), y]

db

da

(5.5)

 

− f [a(y), y]

 

.

dy

dy

a(y)

Доказательство. Функцию I(y) можно считать сложной функцией от y : I(y) = Φ[y, a(y), b(y)]. Дифференцируя эту функцию по

правилу дифференцирования сложной функции от нескольких аргументов, получаем

 

∂I

 

 

∂Φ

 

∂Φ

 

da

 

∂Φ

 

∂b

(5.6)

 

 

=

 

 

 

+

 

 

·

 

+

 

·

 

.

∂y

 

∂y

 

∂a

dy

∂y

∂y

По теореме 5.2

 

 

 

 

 

 

b(y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Φ

=

 

Z

fy(x, y)dx.

 

 

 

(5.7)

 

 

 

 

 

 

 

 

 

 

 

 

∂y

 

 

 

 

a(y)

По правилу дифференцирования интеграла по верхнему пределу

∂Φ

= −f [a(y), y],

(5.8)

 

 

 

∂a

 

∂Φ

= f [b(y), y].

(5.9)

 

∂b

 

 

 

Внося (5.7), (5.8), (5.9) в (5.6), получаем (5.5). Теорема доказана.

cos y

Пример 5.4. Найти I(y), если I(y) = e 1−x2y dx.

 

 

 

 

 

 

 

 

 

 

sin y

 

 

 

Решение. По формуле (5.5)

 

R

 

 

 

 

cos y ey

 

 

 

 

 

 

ey

 

sin y ey

 

cos y.

 

 

 

 

 

 

 

1−cos2 y

1−sin2 y

I(y) =

1−x2

1

 

x2

dx

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin y

110

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