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

Логинов Избранные разделы курса Векторный анализ (теория и примеры) 2009

.pdf
Скачиваний:
49
Добавлен:
17.08.2013
Размер:
1.11 Mб
Скачать

Рис. 4.2

 

4.1.2.

Дифференцирование интегралов,

 

 

 

зависящих от параметра

 

 

Теорема (Лейбница). Если f и f

непрерывны в [a,b]×[c,d] , то

 

 

y

 

 

 

b

 

 

 

b

 

F(y) = f ( x, y)dx дифференцируема на [c,d] и F

= f (x, y) dx.

a

 

 

y

a

y

 

 

 

 

Рассмотрим область типа В, указанную на рис. 4.3, и функцию f ,

определенную на прямоугольнике R, содержащем область D.

 

 

Рис. 4.3

 

 

 

Теорема.

Если f и ее производная

f непрерывны на R, x1(y),

 

 

 

y

 

 

x2(y) имеют

непрерывные на [c,d]

производные,

то

функция

 

 

81

 

 

 

x2 ( y)

F(y) = f ( x, y)dx также имеет производную

x1( y )

 

 

 

F

x2 ( y)

f ( x, y)

 

=

dx + f (x2 ( y), y )x2 ' ( y) - f (x1 ( y), y )x1 ' ( y) .

y

 

x1

( y )

y

 

 

 

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

4.2.1.Равномерная сходимость несобственного интеграла, зависящего от параметра

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

b

 

Φ( y) = f ( x, y)dx , − ∞ < a < b ≤ +∞ , y Y.

(4.1)

a

Предположим, что при некоторых y интеграл (4.1) является несобственным. Так, если − ∞ < a < b ≤ +∞ и при некотором y интеграл (4.1) имеет единственную особенность в b, то условием сходимости интеграла (4.1) будет существование конечного предела

η

lim f (x, y)dx .

η→b−0 a

Если при заданном y интеграл сходится, то для любого η [a,b)

b

интеграл f (x, y)dx (называемый остатком) будет существовать, и

η

b

условие сходимости можно записать в видеηlimb−0 f ( x, y)dx = 0 . В

η

случае расходимости этого интеграла естественно считать, что ус-

b

ловие ηlimb−0 f ( x, y)dx = 0 не выполнено. Таким образом, условие

η

сходимости будет в дальнейшем записываться в виде

82

b

lim f (x, y)dx = 0 .

η→b−0 η

b

Определение. Пусть интеграл с параметром f ( x, y)dx для

a

всех или для некоторых y Y имеет единственную особенность в точке b (если b конечное, интеграл 2-го рода) или в +∞ (интеграл 1-го рода). Сходящийся на Y интеграл называется равномерно сходящимся на Y, если выполнено условие

 

 

b

 

 

ε >0δ >0η(b-δ,b) y Y:

f (x, y)dx

< ε

(для интеграла 2-го

 

 

η

 

 

рода),

 

 

+∞

 

 

 

 

 

ε >0 Mη(M,+ ) y Y:

f ( x, y)dx

< ε

(для интеграла 1-го

 

η

 

 

рода).

 

Признак Вейерштрасса равномерной сходимости. Если суще-

ствует функция g(x), определенная на отрезке [a,b) (b – конечное или +), интегрируемая на любом отрезке [a, η), η (a,b) и такая,

b

что: 1) |f(x,y)| g(x), a ≤ x < b, y Y; 2) g ( x)dx сходится, то ин-

 

 

 

 

 

 

 

a

b

 

 

 

 

 

 

теграл f ( x, y)dx сходится равномерно на множестве Y.

a

 

 

 

 

 

 

Утверждение следует из неравенств

 

 

 

b

 

b

 

b

 

 

 

 

f (x, y)dx

 

f ( x, y)

 

dx g( x)dx .

 

 

 

 

 

 

 

η

 

η

 

η

Теорема (переход

к пределу под

 

знаком интеграла). Пусть

− ∞ < a < b ≤ +∞ и функция f(x,y) определена и непрерывна на от- резке [a,b) по x для всех y Y. Если для любых η(a,b) функция f(x,y)

83

равномерно сходится к функции g(x) на отрезке [a,b -η] при y→y0,

b

интеграл f ( x, y)dx равномерно сходится на Y и интеграл

a

b

g ( x)dx сходится, то имеет место равенство(рис. 4.4)

a

b b

lim f ( x, y)dx = g ( x)dx .

yy0

a

a

 

Рис. 4.4

Критерий Коши равномерной сходимости (интеграла 2-го

b

рода). Для равномерной сходимости интеграла f ( x, y)dx необ-

a

ходимо и достаточно, чтобы

η''

ε >0 δ>0 y Y η′,η′′ (b-δ,b): f (x, y)dx < ε .

η'

84

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

Теорема. Если f(x,y) определена и непрерывна на [a,b)×[c,d] и

b

интеграл Φ(y) = f ( x, y)dx сходится равномерно на [c,d] , то

a

этот интеграл является непрерывной функцией (рис. 4.5).

Рис. 4.5

4.2.3. Интегрирование интегралов, зависящих от параметра

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

b

[a,b)×[c,d], интеграл Φ(y) = f ( x, y)dx сходится равномерно на

a

bd

[c,d] , существует интеграл dxf ( x, y)dy , то

ac

d d b b d

Φ( y)dx = dyf ( x, y)dx = dxf ( x, y)dy .

c c a a c

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

b

[a,b)×[c,d), интеграл f ( x, y)dx сходится равномерно на любом

a

d

[c,η] , интеграл f ( x, y)dy сходится равномерно на любом [a,ξ] и

c

85

d b

существует один из повторных интегралов dyf ( x, y) dx ,

c a

bd

dxf ( x, y) dy , то существует и другой и выполняется равенст-

ac

во

d

b

b

d

dyf ( x, y)dx = dxf ( x, y)dy .

c

a

a

c

4.2.4.

Дифференцирование интегралов,

 

 

зависящих от параметра

Лемма. Если функция f(x,y) непрерывна на [a,b)×[c,d] , то схо-

 

b

 

димость интеграла

f ( x, y)dx эквивалентна условию: для любой

 

a

 

последовательности

ηk→b, η0=a , ηk [a,b)

сходится ряд

ηk +1

f ( x, y)dx . Аналогично для равномерной сходимости: равно-

k =0 ηk

b

мерная сходимость интеграла f ( x, y)dx на множестве Y эквива-

 

 

a

 

 

лентна условию: для любой

последовательности ηk→b, η0=a,

ηk [a,b)

равномерно на

Y сходится функциональный ряд

ηk +1

 

 

 

 

 

 

f ( x, y)dx .

 

 

 

 

 

k =0 η

 

 

 

 

 

 

k

 

 

 

 

 

 

Теорема. Пусть функции f(x,y) и

f

непрерывны на [a,b)×[c,d] .

 

 

 

 

y

 

 

b

 

 

 

b

f ( x, y)

 

Если f ( x, y)dx сходится для всех y, а

dx сходится рав-

 

a

 

 

 

a

y

 

 

 

 

 

86

b

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

 

 

 

 

a

на этом отрезке и

 

 

 

 

 

d

b

b

f ( x, y)

 

 

f ( x, y)dx =

dx .

 

 

 

 

dy

a

a

y

 

 

 

 

4.2.5. Функции Эйлера

Гамма-функция Эйлера Г(p) = x p−1ex dx , p > 0, непрерывна на

0

( 0 , ). Для доказательства этого утверждения разобьем интеграл

1 ∞

на два интеграла Г(p) = x p−1ex dx + x p−1ex dx . Докажем непре-

 

0

1

 

1

рывность функций x p−1ex dx ,

x p−1ex dx на ( 0 , ). Из нера-

 

0

1

1

1

 

венства x p−1ex dx xε−1ex dx , p , A], и сходимости интеграла

0

0

1

 

xε−1ex dx

по признаку Вейерштрасса следует, что интеграл

0

 

1

 

x p−1ex dx сходится равномерно на [ε , A] и, следовательно, явля-

0

ется непрерывной функцией на этом множестве [ε , A].

Для доказательства непрерывности второго интеграла выпишем

∞ ∞

неравенство x p−1ex dx x A−1ex dx , p , A]. Из сходимости

1

1

интеграла x A−1ex dx следует, что интеграл x p −1ex dx сходится

1 1

87

равномерно на [ε , A] и, следовательно, является непрерывной функцией на множестве [ε , A].

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

 

 

 

( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= x p−1exy dx.

 

 

 

 

 

 

 

 

(4.2)

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это равенство получается после замены x → xy . Действительно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(p) = x p−1ex dx = y p−1x p−1exy ydx = y p x p−1exy dx .

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Бета-функция Эйлера: В(p,q) = x p−1 (1 − x)q−1 dx , p > 0 , q > 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сделаем замену x =

 

y

 

,

dx =

 

 

 

dy

 

 

, тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + y

(1+ y)2

 

 

 

 

 

y

p−1

1

 

 

 

 

 

 

 

 

dy

 

 

 

 

y

p−1

В(p,q) = =

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

dy ,

(1

+ y)

p−1

 

(1 + y)

q−1

 

(1 + y)

2

(1 + y)

p+q

0

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда получается формула

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

p−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B( p, q) =

 

 

 

 

 

 

 

 

 

dy.

 

 

(4.3)

 

 

 

 

 

(1

+ y)

p+q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.2.6. Свойства функций Эйлера

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p + q)

 

 

 

 

Из формулы (4.2) следует, что

 

 

x p+q−1ex(1+ y) dx .

 

 

 

 

 

 

 

 

 

 

=

 

(1 + y)

p+q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Умножая обе части этого равенства на y p −1 , получим равенство

 

 

 

 

 

 

y p−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p + q)

 

 

 

 

 

 

 

 

= y p−1 x p+q−1exy ex dx .

 

(1 + y) p+q

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее интегрируем от 0 до :

88

 

y

p−1

( p + q)

 

 

dy = y p−1dy x p+q−1ex( y+1) dx =

(1

+ y) p+q

0

0

0

 

 

 

= xq−1ex dx ( xy) p−1 exy xdy .

0 0

Откуда, используя (4.3), получим: B( p, q) = ( p) (q) . В част-

 

 

 

 

( p + q)

y

p−1

 

π

 

ности, B( p,1 − p) = ( p) (1 − p) =

 

dy =

,0 < p < 1 .

1 + y

 

0

 

sin pπ

 

 

 

 

 

Непосредственно из определения следует, что Г(1) = 1, Г(p+1) =

=p Г(p).

Отметим, что из этой формулы следует, что гамма-функцию

достаточно знать на интервале (0, 1/2). Интеграл x p−1ex dx схо-

0

дится (p >0) , а интеграл x p−1 ln k x ex dx сходится равномерно на

0

любом отрезке [ε , A ] для 0 < ε < A. Поэтому интеграл можно дифференцировать по параметру. Докажем равномерную сходи-

мость интегралов x p−1 ln k x ex dx .

0

Вокрестности нуля |ln x| C1 (δ) для δ > 0 существует C1(δ). В xδ

окрестности бесконечности |ln x| C2 (δ)xδ для δ > 0 существует

C2(δ). Равномерная сходимость интеграла Г(k)(p)= x p−1 ln k x ex dx

 

 

 

 

0

на

любом отрезке [ε , A ]

следует из оценок

 

 

1

 

 

 

x p −1 ln k x ex dx

xε−1 | ln k x | ex dx + x A−1 | ln k x | ex dx

 

0

 

0

1

89

1 ∞

C1k (δ)xε−kδ−1ex dx + + C2k (δ) x A+kδ−1ex dx , для всех p , A].

0 1

Здесь для ε >0 следует выбрать δ так, чтобы ε - k δ оставалось больше нуля.

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

Формула Фруллани. Функция f(x) непрерывна и интеграл

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f ( x)

dx существует для любого A > 0, тогда справедлива формула

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Фруллани

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (ax) − f (bx)

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx = f(0) ln

, (a>0, b>0).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

x

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство. b f ' ( xy)dy =

1

b d f ( xy) =

f (bx) − f (ax)

. Ин-

 

x

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

b

 

 

 

 

 

 

f (bx) − f (ax)

dx = dxf '( xy)dy = dy f ' (xy)dx =

 

 

 

 

 

 

0

 

 

 

x

0

 

a

 

 

 

 

a

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

dy

b

 

 

 

 

 

 

dy

 

 

 

 

b

 

 

 

 

=

 

f ' (u)du = ( f (∞) − f (0))

= f (0) ln

,

 

 

 

 

 

 

 

 

 

a

y

0

 

a

 

 

 

 

 

 

y

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда и следует требуемое утверждение.

 

 

 

 

 

 

Интегрированием по частям вычисляются интегралы

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

β

 

 

 

 

e−αx cos βx dx =

 

, α ≥ 0,

e−αx sin βx dx =

, α ≥ 0 .

 

 

 

 

α2 + β2

α 2 + β2

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Другой способ вычисления этих интегралов: положим γ =–α+iβ, тогда

90

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