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

Краткий курс математического анализа. Том 2

.pdf
Скачиваний:
329
Добавлен:
10.05.2015
Размер:
4.49 Mб
Скачать

270 Гл. 5. Интегральное исчисление функций многих переменных

и теорему 4 из п. 49.2 об интегрировании собственных интегралов, зависящих от параметра, получим

d

b

 

=

d

ηn

y) dx dy =

 

dy

,

 

 

 

,

 

c

a f (x

y) dx

(50.14) c n=1

ηn 1f (x

 

 

 

d

ηn

 

 

 

ηn

d

 

 

 

 

 

 

 

 

 

 

 

=

dy

f (x, y) dx =

 

dx

f (x, y) dy. (50.17)

 

 

n=1

 

 

(49.12)

n=1

 

 

 

 

ηn−1

 

 

 

 

c

 

 

c

 

 

 

ηn−1

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

ηn d

n

ηk d

 

 

 

dx

f (x, y) dy =

lim

dx

f (x, y) dy =

n=1

c

n→∞ k=1

 

c

 

 

ηn−1

ηk−1

 

 

 

 

 

 

 

ηn

d

 

 

 

 

=

lim

dx f (x, y) dy, (50.18)

 

 

 

 

 

n→∞

 

 

 

 

 

 

a

c

причем этот предел в силу (50.17) существует и равен интегралу

db

dy f (x, y) dx.

c a

Вспомнив теперь, что в качестве последовательности n} можно

взять любую последовательность такую, что lim ηn = b, a ηn < b,

 

 

 

n→∞

 

 

n = 1, 2, ..., согласно определению несобственного интеграла оконча-

тельно получим

 

 

 

 

 

ηn

d

b

d

 

 

lim

dx

f (x, y) dy =

dx f (x, y) dy.

(50.19)

n→∞

 

 

 

 

 

a

c

a

c

 

 

Из (50.17)–(50.19) следует (50.16).

 

 

 

Т е о р е м а 5 (правило Лейбница). Если функции f (x, y) и

∂f (x, y)

 

∂y

 

 

 

 

непрерывны на конечном или бесконечном прямоугольнике P , инте-

b

 

 

b

 

 

грал f (x, y) dx сходится,

а интеграл

∂f (x, y) dx равномерно схо-

 

 

 

∂y

 

 

a

 

 

a

 

 

 

 

 

b

 

 

дится на отрезке [c, d], то функция Φ(y) = f (x, y) dx дифференци-

a

§ 50. Несобственные интегралы, зависящие от параметра

271

руема на отрезке [c, d] и

bb

d

f (x, y) dx =

∂f (x, y)

dx.

dy

∂y

 

 

aa

Согласно формуле (50.14) имеем

b

ηn

 

f (x, y) dx =

f (x, y) dx,

an=1 ηn−1

b

 

ηn

 

 

∂f (x, y)

dx =

 

∂f (x, y)

dx,

 

 

 

 

∂y

 

 

∂y

 

 

 

 

 

a

 

n=1

 

 

 

ηn−1

 

(50.20)

(50.21)

(50.22)

причем рад (50.21) сходится, а ряд (50.22) равномерно сходится на отрезке [c, d]. В силу теоремы 5 п. 49.2 каждый член ряда (50.21) дифференцируем, и для его производной справедлива формула

 

ηn

 

ηn

 

 

d

 

f (x, y) dx =

 

∂f (x, y)

dx.

(50.23)

dy

 

∂y

 

 

 

 

 

ηn−1

 

ηn−1

 

 

Поэтому из теоремы о почленном дифференцировании рядов (теорема 9 из п. 31.4) следует, что сумма ряда (50.21), т. е. интеграл Φ(y) =

b

= f (x, y) dx, является дифференцируемой на отрезке [c, d] функцией,

a

и для его производной имеет место формула (50.20):

 

b

 

 

 

ηn

 

 

 

 

 

 

d

f (x, y) dx =

d

f (x, y) dx =

 

 

 

 

dy

 

 

 

 

 

 

 

 

(50.21) dy

n=1

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

ηn

 

 

ηn−1

ηn

 

b

 

 

 

 

 

 

 

 

d

 

 

 

 

∂f (x, y)

 

 

∂f (x, y)

dx.

 

=

f (x, y) dx =

 

dx =

 

 

dy

 

∂y

 

∂y

 

n=1

ηn−1

 

 

(50.23) n=1

 

 

 

 

 

 

 

 

 

 

ηn−1

 

a

 

 

50.3. Интегралы Эйлера. Интеграл

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

Γ(s) = xs−1e−x dx,

 

 

 

(50.24)

0

зависящий от параметра s, называется гамма-функцией, а интеграл

1

 

B(p, q) = xp−1(1 − x)q−1 dx,

(50.25)

0

272 Гл. 5. Интегральное исчисление функций многих переменных

зависящий от двух параметров p и q, называется бета-функцией. Гамма-функция и бета-функция называются эйлеровыми интеграла-

ми. Разобьем интеграл (50.24) на два интеграла:

+

1

+

 

 

 

 

Γ(s) =

xs−1e−x dx = xs−1e−x dx +

 

xs−1e−x dx,

(50.26)

0

0

1

 

 

 

 

и заметим, что

 

 

 

 

 

 

 

xs−1e−x xs−1 при

x → 0,

x

 

.

xs−1e−x = xs−1e−x/2e−x/2 = o(e−x/2)

при

→ ∞

 

 

 

 

 

1

Так как интеграл xs−1 dx сходится при s > 0 и расходится при

0

+

s 0, а интеграл e−x/2 dx сходится, то в силу равенства (50.26)

1

формула (50.24) для определения гамма-функции Γ(s) имеет смысл только при s > 0. Иначе говоря, если s > 0, то интеграл (50.24) сходится; если же s 0, то он расходится.

Аналогично,

1

B(p, q) = xp−1(1 − x)q−1 dx =

0

1/2

1

 

= xp−1(1 − x)q−1 dx +

xp−1(1 − x)q−1 dx, (50.27)

где

0

1/2

xp−1(1 − x)q−1 xp−1

при x → 0,

 

xp−1(1 − x)q−1 (1 − x)q−1 при x → 1.

1/2

Поскольку интеграл xp−1 dx сходится при p > 0 и расходится

0

1

при p 0, а интеграл (1 − x)q−1 dx сходится при q > 0 и расходится

1/2

при q 0, то в силу равенства (50.27) формула (50.25) для определения бета-функции B(p, q) имеет смысл только при p > 0 и q > 0, т. е. при этих значениях параметров интеграл (50.25) сходится, а если хотя бы одно из условий p > 0 и q > 0 не выполняется, то интеграл (50.25) расходится.

§ 50. Несобственные интегралы, зависящие от параметра

273

50.4. Интеграл Дирихле. Рассмотрим на примере интеграла

 

+

 

I(α) =

 

sin αx

dx

(50.28)

 

 

 

 

x

 

 

0

 

 

 

метод вычисления определенных интегралов с помощью интегралов, зависящих от параметра.

Можно показать, что первообразная для функции

sin αx

не выра-

x

 

 

жается через элементарные функции и тем самым интеграл (50.28) нельзя вычислить с помощью формулы Ньютона–Лейбница.

Интеграл (50.28) сходится при всех значениях α. В самом деле, если α = 0, то, очевидно, I(0) = 0. Если же α = 0, то, сделав замену переменной t = αx при α > 0 и t = −αx при α < 0, получим

 

+

 

sin t

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

dt = I(1),

если

α > 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I(α) =

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin t

 

 

 

 

 

 

 

 

 

 

 

0

 

t dt = −I(1),

если

α < 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интеграл же I(1) сходится (см. п. 29.6), поэтому сходится и инте-

грал I(α).

Интеграл (50.28) называют интегралом Дирихле или разрывным множителем Дирихле (поскольку этот интеграл как функция пара-

метра α, −∞ < α < +, имеет разрыв при α = 0).

Если формально продифференцировать по правилу Лейбница ин-

+

теграл (50.28), то получится расходящийся интеграл

cos αx dx.

0

Поэтому для вычисления интеграла I(α) рассмотрим более общий

интеграл

+

I(α, β) = e−βx sin αx dx. x

0

Продифференцировав по правилу Лейбница интеграл I(α, β) по параметру α, получим интеграл

+

e−βx cos αx dx,

(50.29)

0

который при любом фиксированном β > 0 равномерно сходится относительно параметра α, −∞ < α < +∞. Следовательно, к интегралу I(α, β) при β > 0 можно применить правило Лейбница дифференцирования по параметру α. Заметив, что интеграл (50.29)

274 Гл. 5. Интегральное исчисление функций многих переменных

вычисляется, например, двукратным интегрированием по частям (см. п. 22.4), причем

 

 

+

 

 

 

 

 

β

 

 

 

 

e−βx cos αx dx =

 

,

 

 

2

2

 

 

0

 

 

 

 

 

α + β

 

 

 

 

 

 

 

 

 

 

 

 

получим

∂I(α, β) =

β

 

. Отсюда при β > 0 можно найти I(α, β):

α2 + β2

 

∂α

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

I(α, β) =

β dt

+ C(β) = arctg

α

+ C(β).

 

t2 + β2

β

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Но I(0, β) = 0, следовательно, C(β) = 0. Итак,

 

 

I(α, β) = arctg

α

,

β > 0.

 

 

 

 

 

 

 

 

 

 

 

β

 

 

 

Нас, однако, интересует значение интеграла I(α, β) при β = 0. Покажем, что здесь возможен предельный переход под знаком интеграла при β → +0. Зафиксируем произвольно число b > 0 и покажем, что интеграл I(α, β) при любом фиксированном α = 0 равномерно сходится по параметру β на отрезке [0, b]. В самом деле, вычислив

по частям интеграл по полуинтервалу [η, +), η > 0, от функции

e−βx sin αx , получим x

+

e−βx

 

 

 

dx =

 

e−βx

α cos

 

 

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η

 

x

 

x

α2 + β2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin αx

 

 

 

 

1

 

 

 

 

 

αx +

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β sin αx

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α cos αx + β sin αx dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− e−βx

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α2 + β2

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Зададим произвольно ε > 0 и выберем ηε так, чтобы при η > ηε

выполнялись неравенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

αη + β sin αη

 

 

α + b 1

 

 

 

ε

 

 

 

 

 

 

 

 

 

 

 

 

e−ηβ

α cos

 

 

 

 

 

 

 

|

 

α| 2

 

 

<

 

 

,

 

 

 

 

 

 

 

 

 

 

η

 

 

 

α2 + β2

 

 

 

 

 

η

 

2

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

αx + β sin αx dx

 

 

 

α + b

 

 

 

ε

 

 

 

 

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η

 

 

 

 

 

 

 

 

 

 

 

e−βx α cos α2 + β2

 

 

 

 

x2

 

|

 

α|

2

 

x2

< 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда при

η > ηε для всех β

 

[0, b] получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

sin αx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e−βx

 

x

dx

< ε,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η

§ 50. Несобственные интегралы, зависящие от параметра

275

что и доказывает равномерную сходимость интеграла I(α, β) по параметру β на отрезке [0, b]. Теперь, в силу следствия теоремы 3 п. 50.2,

I(α) = I(α, 0) =

lim I(α, β) =

lim

 

arctg

 

α

=

π

sign α.

 

 

0

 

 

2

 

 

 

 

β→+

0

 

β→+

 

 

 

 

β

 

 

 

 

 

 

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

sin αx

π/2,

если

 

α > 0,

 

 

I(α) =

0,

если

 

α = 0,

 

 

x

 

 

 

 

0

 

 

 

π/2,

если

 

α < 0.

 

 

 

 

 

 

2

 

 

 

 

 

 

2 +sin xt

 

Отсюда между прочим следует, что

sign x =

 

 

I(x) =

 

 

 

 

dt,

π

π

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

т. е. функция sign x может быть задана не только словесным описанием, но и аналитическим выражением.

Г л а в а 6

ГАРМОНИЧЕСКИЙ АНАЛИЗ

§51. Тригонометрические ряды Фурье

51.1.Основные понятия. В этом параграфе будут изучаться

ряды вида

a0 + an cos nx + bn sin nx,

(51.1)

n=1

 

an R, n = 0, 1, 2, ..., bn R, n = 1, 2, ...,

 

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

1, cos x, sin x, ..., cos nx, sin nx, ...,

(51.2)

n = 1, 2, ...

(называемой тригонометрической системой или системой простых гармоник), умноженные на постоянные an, bn, называемые коэффициентами тригонометрического ряда (51.1).

Л е м м а 1. Функции тригонометрической системы (51.2) имеют следующие свойства:

π

 

m = n,

 

cos nx cos mx dx = 0,

m, n = 0, 1, 2, ...;

−π

 

 

 

π

 

m = n,

 

sin nx sin mx dx = 0,

m, n = 0, 1, 2, ...;

−π

 

 

(51.3)

π

 

 

 

 

 

cos nx sin mx dx = 0,

m, n = 0, 1, 2, ...;

−π

 

 

 

π

π

 

 

cos2 nx dx =

sin2 nx dx = π,

n = 1, 2, ...

−π −π

§ 51. Тригонометрические ряды Фурье

277

Действительно, например,

π

 

 

 

1

 

π

 

 

 

 

 

 

 

 

 

 

 

 

cos nx cos mx dx =

 

[cos (n + m)x + cos (n − m)x] dx =

2

−π

 

 

 

 

−π

 

π

 

 

 

 

 

 

 

π

 

 

 

 

 

=

sin (n + m)x

 

 

+ sin (n − m)x

 

= 0, n = m;

π

 

π

 

 

 

 

 

 

 

 

 

2(n + m)

 

 

π

 

2(n − m)

 

 

π

 

2

1

 

 

 

 

 

 

 

 

 

sin 2nx

 

 

 

 

 

 

 

cos nx dx =

2

 

(1 + cos 2nx) dx = π +

4n

 

 

π

= π,

n = 1, 2, ... .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−π

 

−π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично доказываются другие равенства (51.3). Прежде всего возникает вопрос, как найти коэффициенты триго-

нометрического ряда (51.1), если известна его сумма. Ответ на него легко дать, когда ряд (51.1) сходится равномерно на всей числовой оси.

Те о р е м а 1. Если тригонометрический ряд (51.1) равномерно сходится на всей числовой оси и f — его сумма, т. е.

 

 

 

f (x) = a0 +

an cos nx + bn sin nx,

(51.4)

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

то

a0

=

 

1

 

π

f (x) dx,

an =

1

 

π

f (x) cos nx dx,

 

 

 

2π

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−π

π

 

 

 

 

−π

(51.5)

 

 

bn =

1

f (x) sin nx dx,

 

n = 1, 2, ...

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−π

Поскольку ряд, стоящий в правой части равенства (51.4), равномерно сходится, то его можно почленно интегрировать (теорема 8 из п. 31.4) на отрезке [−π, π]. Поэтому, проинтегрировав обе части равенства (51.4), будем иметь

π

π

 

 

 

 

 

π a0

 

 

 

 

 

 

 

 

 

 

 

π f (x) dx =

 

+ n=1 an cos nx + bn sin nx dx =

 

 

 

 

π

 

π

π

 

 

 

= a0

dx +

an cos nx dx + bn

sin nx dx = 2πa0.

 

 

 

 

 

 

 

 

−π

n=1

−π

−π

 

 

 

 

 

 

Отсюда сразу следует первая формула (51.5).

Если обе части равенства (51.4) умножить на cos mx или на sin mx, m N, то, поскольку эти функции ограничены на числовой оси, в правой части будут снова стоять ряды, равномерно сходящиеся

278

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

на всей числовой оси (см. замечание 3 в п. 31.2). Для доказательства оставшихся формул (51.5) достаточно проинтегрировать обе части каждого из получившихся равенств и воспользоваться формулами (51.3). Например,

π

f (x) cos mx dx =

−π

π

 

 

 

 

 

=

 

π a0 cos mx + n=1 an cos nx cos mx + bn sin nx cos mx dx =

 

 

 

π

 

π

 

 

 

= a0

 

cos mx dx +

an cos nx cos mx dx +

 

 

 

 

 

 

 

 

−π

 

n=1

 

 

 

 

 

 

 

−π

π

 

 

 

 

 

 

 

+ bn

sin nx cos mx dx = πam.

 

 

 

 

 

 

 

 

(51.3)

 

 

 

 

 

 

 

 

−π

Отсюда am =

1

π

ϕ(x) cos mx dx.

 

π

 

 

 

 

 

 

 

 

 

 

−π

З а м е ч а н и е 1. Напомним (см. п. 29.5), что функция f называется абсолютно интегрируемой на некотором конечном или бесконечном промежутке с концами a и b, −∞ a < b +, если существует конечное число точек xj , j = 0, 1, ..., k, расширенной числовой прямой R таких, что:

1)a = x0 < x1 < ... < xj < ... < xk = b;

2)функция f интегрируема по Риману на любом отрезке [ξ, η]

(xj−1, xj ), j = 1, 2, ..., k;

 

xj

 

 

 

 

3) интегралы

f (x) dx абсолютно сходятся,

j = 1, 2, ..., k.

 

xj−1

 

 

 

 

При выполнении этих условий

 

 

 

 

b

k

xj

 

 

 

 

 

 

 

 

def

 

 

f (x) dx.

(51.6)

 

f (x) dx =

 

 

a

j=1 x

 

 

 

 

 

j−1

 

 

Всякое множество {x0, x1, ..., xk} указанных точек называется правильным разбиением на интервалы (xj−1, xj ) промежутка интегриро-

вания абсолютно интегрируемой функции f. Отметим, что интеграл (51.6) не зависит от выбора правильного разбиения промежутка интегрирования.

З а м е ч а н и е 2. Если функция f абсолютно интегрируема на отрезке [−π, π], то из неравенств |f (x) cos nx| |f (x)|, |f (x) sin nx|

§ 51. Тригонометрические ряды Фурье

279

|f (x)| следует, что интегралы

 

π

π

 

f (x) cos nx dx,

f (x) sin nx dx, n = 1, 2, ...,

 

−π

−π

 

также абсолютно (см. замечание из п. 29.5), а следовательно, и просто сходятся. Таким образом, формулы (51.5) имеют смысл для любой абсолютно интегрируемой на отрезке [−π, π] функции.

О п р е д е л е н и е 1. Если функция f абсолютно интегрируема на

отрезке [−π, π], то тригонометрический ряд (51.1), коэффициенты которого заданы формулами (51.5), называется тригонометрическим рядом Фурье (или, короче, рядом Фурье) функции f , его коэффициенты (51.5) — коэффициентами Фурье функции f (по тригонометри-

ческой системе функций (51.2)), а частичные суммы порядка n ряда Фурье — суммами Фурье порядка n функции f.

Если ряд (51.1) является рядом Фурье функции f , то пишут

f a0 + an cos nx + bn sin nx.

n=1

Используя введенные термины, теорему 1 можно перефразировать следующим образом.

Т е о р е м а 1. Равномерно сходящийся тригонометрический ряд является рядом Фурье своей суммы.

За м е ч а н и е 3. Если у интегрируемой (в собственном или не собственном смысле) функции изменить ее значения на конечном множестве точек, то она останется интегрируемой, и значение инте-

грала от нее по указанному промежутку не изменится. Поэтому две абсолютно интегрируемые на отрезке [−π, π] функции, отличающиеся друг от друга лишь на конечном множестве точек, имеют одинаковые ряды Фурье.

За м е ч а н и е 4. Напомним еще, что число T > 0 называется периодом функции f , если для любого числа x, принадлежащего области определения X R функции f , числа x + T и x − T также принад-

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

Если функция f определена, например, на отрезке [−π, π], то при f (−π) = f (π) ее нельзя продолжить на всю числовую ось так, чтобы получилась 2π-периодическая функция, а ее сужение на полуинтервале [−π, π) можно продолжить указанным образом: следует положить

def

[ π, π),

k = 0, 1, 2,

f (x + 2πk) = f (x), x

 

± ± ...

Функция f , очевидно, 2π-периодическая и на отрезке [−π, π] отличается от функции f , быть может, только в одной точке x = π.