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

А. Р. Лакерник. Высшая математика.Краткий курс

.pdf
Скачиваний:
5
Добавлен:
10.10.2023
Размер:
3.65 Mб
Скачать

IV. Функции нескольких переменных. Интегралы, зависящие от п араметра

x - x0 = t cosa, y - y0 = t cosb, z - z0 = t cos g,

откуда x = x0 + t cosα, y = y0 + t cosβ, z = z0 + t cosγ, т.е. x, y, и z есть некоторые функции t.

Теперь u = u(x,y,z) = u(x(t),y(t),z(t)) = j(t) – некоторая функция от t, ϕ(0) = u(x0,y0,z0) è

¶u

(M

0

) =

lim

u(x,y,z) - u(x0,y0,z0 )

 

= lim

j(t) - j(0)

 

= j¢(0).

¶l

M M

t

 

 

M0M 0

t0

 

 

 

 

 

 

0

 

 

 

 

 

По правилу дифференцирования сложной функции (см. (12.9))

j¢(t) =

¶u dx +

¶u dy

+

¶u dz =

¶u cos a +

¶u cosb +

¶u cos g.

 

¶x dt

¶y dt

 

¶z dt

¶x

 

 

¶y

¶z

При t = 0 частные производные ¶u

, ¶u

,

¶u

берутся в точке M0 è

 

 

 

 

¶x

¶y

 

¶z

 

 

ul (M0 ) = ux (M0 )cosa + uy (M0)cosb + uz (M0 )cos g.

Таким образом, нами доказана следующая теорема.

Теорема 12.14. Если на множестве T функция u(M ) = u(x,y,z) имеет непрерывные частные производные, то в любой точке M0 T эта функция имеет производную по любому направлению l :

¶u

(M0) =

¶u

(M0 )cosa +

¶u

(M0 )cosb +

¶u

(M0 )cos g, (12.41)

¶l

 

¶x

 

¶y

 

¶z

 

где a, b, g – углы l с осями координат 0x, 0y, 0z соответственно.

Как известно, производная функции – это скорость изменени я этой функции. Производная функции по направлениюl – это скорость изменения функции по этому направлению. По какому направл ению функция изменяется быстрее всего, т.е. по какому направлен ию производная будет наибольшей?

Для ответа на этот вопрос сначала заметим, что вектор e = {cos a, cosb, cos g} – это единичный вектор (т.е. вектор длины 1), направленный вдоль оси l (координаты такого вектора как раз и равны косинусам его углов с положительными направлениями осей координат). Введем еще один вектор.

242

12. Функции нескольких переменных

Определение 12.32. Градиентом функции u(M ) = u(x,y,z) в точке

M0

 

 

ì¶u

,

¶u

,

¶uü

(все производные берутся

 

 

называется вектор grad u = н

¶x

¶y

ý

в точке M0).

î

 

 

¶z þ

 

 

 

 

 

 

 

 

 

 

Правая часть формулы (12.41) равна скалярномó ïðîèзведению

этих двух векторов; таê как скалярное произведение ab = a b cos j, где j – угол между a и b, то

¶u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= grad u ×

e

=

grad u

 

 

e

 

cosj =

grad u

cosj,

¶l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{

 

 

 

 

 

1

 

 

 

 

 

 

где j – угол между grad u и осью l. Последнее выражение будет наибольшим при cosj = 1, т.е. j = 0.

Таким образом, производная функции в точке M0 будет наибольшей по направлению вектора grad u в данной точке и эта наибольшая

производная по направлению равна grad u(M0) .

Замечание. Из формулы (12.40) следует, что значение производной по направлению не зависит от выбора системы координат xyz. Тогда вектор

grad u(M0) тоже не зависит от выбора системы координат: он направлен по тому направлению, производная по которому наибольшая, и длина его равна этой наибольшей производной по направлению.

Пример. Найти производную функции u = xyz

в точке M0 (1, 2, 3) ïî

направлению M0M к точке M (2, 4, 6) и наибольшую производную по на-

правлению в точке M0 .

 

 

 

 

Ð å ø å í è å

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= {2 -1,4 - 2,6 -3} = {1,2,3};

 

 

 

 

=

12 +22 +32

=

14.

 

 

M0M

M0M

 

Из рис. 75 видно, что cosa =

 

1

, cosb =

 

2

, cos g =

3

 

;

 

 

 

14

14

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¶u = yz,

¶u = xz,

¶u = xy;

¶u

(M0) = 6,

¶u

(M0) = 3,

¶u (M0) = 2;

 

¶x

¶y

¶z

¶x

 

 

 

 

 

 

 

¶y

 

 

 

 

 

 

¶z

 

 

 

 

 

 

 

 

(12.41)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= {6,3,2}

Þ ¶u(M0) = 6

1

+3

 

2

 

 

+2

3

=

 

18

 

»4,81.

 

grad u(M0)

 

 

 

 

14

 

 

 

14

 

 

 

 

 

 

 

 

 

¶l

 

 

14

 

 

 

 

 

 

 

14

 

 

 

 

Наибольшая производная по направлению в точке M0

 

 

 

 

 

 

 

 

 

=

 

{6,3,2,}

 

= 62 + 32 + 22

 

=

49 = 7.

 

 

 

 

 

 

 

 

 

grad u(M0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

243

IV. Функции нескольких переменных. Интегралы, зависящие от п араметра

13.ИНТЕГРАЛЫ, ЗАВИСЯЩИЕ ОТ ПАРАМЕТРА

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

Определение 13.1. Рассмотрим функцию, которая задается в виде интеграла

b

 

I(a) = ò f (x,a)dx,

(13.1)

a

где f (x,a) – некоторая функция двух переменных. Такой интеграл называется интегралом, зависящим от параметра, а переменная α называется параметром.

При различных значениях α правая часть формулы (13.1) дает различные значения интеграла (если он вообще существует), поэ тому этот интеграл действительно является функцией параметра a .

Свойства интегралов, зависящих от параметра

Обозначим через P замкнутый прямоугольник на плоскости xα :

{-¥ < a £ x £ b < +¥, - ¥ < b £ a £ g < +¥}.

Теорема 13.1. Пусть функция f (x,a) непрерывна на прямоугольнике P. Тогда функция I(a) непрерывна на отрезке [b,g] .

¡ Пусть a0 О[b, g]. Надо доказать, что

lim I(a) = I(a0) , ò. å.

 

 

 

 

 

 

 

 

 

 

 

 

α→α0

"e > 0

$d = d(e) > 0: "a,

 

a - a0

 

< d

 

I(a) - I(a0)

 

< e , ò.å. â ñèëó (13.1):

 

 

 

 

 

 

b

 

b

 

 

0

 

 

 

b

0 û

 

 

 

 

 

 

 

 

 

 

 

 

ò

 

ò

 

 

 

 

 

ò ë

 

 

 

 

f (x,a)dx -

 

 

f (x,a

 

 

)dx

=

é f (x,a)

- f (x,a ùdx

< e .

 

 

a

 

a

 

 

 

 

 

 

 

a

 

 

 

 

Возьмем произвольное e > 0 . Так как функция f (x,a) непрерывна на замкнутом прямоугольнике P, то согласно теореме 12.4 она равномерно непрерывна на этом прямоугольнике. Значит, $d > 0:

äëÿ "M1,M2 ÎP , r(M1,M2 ) < d будет выполняться неравенство

 

f (M2 ) - f (M1)

 

<

 

e

.

 

 

 

 

b

- a

 

 

 

 

 

 

 

 

Íî åñëè

 

a - a0

 

< d , òî

 

 

 

(r(M(x,a),M0(x,a0)) = (x - x)2 + (a - a0)2 = a - a0 < d Þ

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

 

 

 

 

 

 

b

 

 

 

 

 

b

 

 

 

 

 

 

 

I(a) - I(a

0

)

 

£

é f (x,a) - f (x,a

0

)ùdx

£

ò

 

f (x,a) - f (x,a

0

)

 

dx £

 

 

 

 

 

 

 

 

 

ò ë

 

 

û

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

£

e

(b - a) = e . x

 

 

 

 

 

 

 

 

 

 

 

 

b - a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 13.2. Пусть функция f (x,a) непрерывна на прямоугольнике P. Тогда

γ

γ

éb

ù

b é

γ

ù

òI(a)da = ò êò f (x,a)dxú da = ò ê

ò f (x,a)daúdx.

β

β

ê

ú

ê

β

ú

 

 

a ë

 

û

 

 

ëa

û

 

¡ Равенство следует из того, что оба повторных интеграла в левой и правой его частях равны двойному интегралу òò f (x,a)dxda

(ñì. ãë. 19). x

P

Теорема 13.3. Пусть функции f (x,a) и ¶f (x,a) непрерывны на

¶a

прямоугольнике P. Тогда при a О [b,g] существует

 

d b

b

¶f (x,a)

I ¢(a) =

 

ò f (x,a)dx = ò

 

dx.

da

¶a

 

 

a

a

 

 

¡ Согласно определению производной, применяя теорему Лагранжа 5.4 и обозначая промежуточную точку между α и a + Da через c = a + qDa, 0 < q <1, имеем

I ¢( )

 

lim I (a + Da) - I(a)

=

lim

1

 

 

éb

f (x,

 

)dx

-

b f (x,

 

)dx

ù

=

 

 

 

 

 

 

ê

 

 

ú

 

 

a =

 

 

Da

 

Δα→0 Da

 

a + Da

 

 

ò

 

a

 

 

 

 

Δα→0

 

 

êò

 

 

 

 

 

 

 

 

 

 

 

ú

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ëa

 

 

 

 

 

 

a

 

 

 

 

 

û

 

 

 

 

 

b f (x,a + Da) - f (x,a)

 

 

 

 

 

 

 

b

¶f (x,a + qDa)

 

 

 

 

 

 

 

 

 

 

 

= Δα→lim0 ò

 

 

 

 

 

 

 

 

 

dx

= Δα→lim0

ò

 

¶a

 

 

dx.

 

 

 

 

 

 

 

 

 

 

 

Da

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

¶f (x,a)

 

 

 

 

Нам надо доказать, что последний предел равен ò

 

 

 

dx , ò.å.

 

 

¶a

 

 

"e > 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

$d = d(e) > 0:"Da,| Da |< d будет выполняться равенство

 

 

 

 

 

 

 

b

¶f (x,a + qDa)

b

¶f (x,a)

 

=

 

b

é

¶f (x,a + qDa)

-

¶f

(x,a)ù

 

 

< e.

 

 

 

 

 

 

ò

 

 

¶a

dx - ò

¶a

 

dx

 

ò

ê

 

¶a

 

 

 

¶a

údx

 

a

 

 

a

 

 

 

 

 

 

a

ë

 

 

 

 

 

û

 

 

 

 

 

244

245

IV. Функции нескольких переменных. Интегралы, зависящие от п араметра

Возьмем произвольное e > 0. Так как функция

¶f (x,a)

 

непрерыв-

 

¶a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

на на замкнутом прямоугольнике P, то она равномерно непрерывна на

этом прямоугольнике, значит, $d > 0:"M1,M2 ÎP,r(M1,M2 ) < d

 

 

 

 

 

 

 

 

¶f (M

2

) - ¶f (M

1

)

 

<

 

 

e

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¶a

 

 

 

¶a

 

 

 

 

 

 

 

 

 

 

 

b - a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Íî åñëè

 

Da

 

< d , то | qDa |< d , тогда, как при доказательстве теоре-

 

 

ìû 13.1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

é

¶f (x,a + qDa)

 

¶f (x,a)ù

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(r(M(x,a + qDa),M0(x,a)) < d Þ

ò

ê

 

 

 

 

 

 

 

 

 

 

 

 

-

 

 

 

údx

£

 

 

 

 

 

 

 

¶a

 

 

 

 

¶a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

ë

 

 

 

 

 

 

 

 

 

 

 

 

û

 

 

 

b

 

¶f (x,a + qDa)

 

 

¶f

(x,a)

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£ ò

 

 

 

¶a

 

 

-

 

 

¶a

 

 

 

 

 

 

dx <

 

 

(b - a) = e. x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b - a

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 (x

 

 

+ a

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ð å ø å í è å

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

dx

 

 

1

 

 

 

 

x

 

b

 

 

1

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ò

 

 

=

 

 

arctg

 

 

 

 

=

 

arctg

 

.

 

 

 

 

 

 

 

 

 

 

 

 

x2

+ a2

 

 

a

a

0

a

a

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дифференцируем обе части по a, используя теорему 13.3:

b

 

 

2a

 

 

 

 

 

 

1

 

 

b

 

1

 

 

 

1

 

æ

 

b ö

 

-ò

 

 

 

dx

= -

 

arctg

+

 

 

 

 

-

Þ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

 

÷

(x

2

+ a

2 2

a

2

a

a

 

 

 

b

2

a

2

0

 

)

 

 

 

 

 

 

 

 

 

1

 

 

è

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ a2

 

 

 

 

 

 

 

 

b

 

 

dx

 

 

 

 

1

 

 

b

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

ò

 

 

 

 

=

 

 

arctg

+

 

 

 

 

 

 

 

 

.

 

 

 

 

(x

2

2

2

 

 

 

3

a

 

 

2

2

 

2

 

 

 

 

0

 

+ a

)

 

 

2a

 

 

 

 

2a

 

(a + b )

 

 

 

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

Теоремы 13.1–13.3 применять к несобственным интегралам возможно только при определенных условиях. Для формулировки этих условий нужно ввести понятие равномерной сходимости несобственных интегралов, зависящих от параметра.

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

b

Рассмотримнесобственныйинтеграл I(a) = ò f (x,a)dx соднойосо-

a

бенностью в точке b. В этом и следующем разделах через P обозначим прямоугольник {-¥ < a £ x £ b £ +¥, - ¥ < b £ a £ g < +¥}.

Как известно, I(a) сходится при a О[b,g], если существует конеч-

 

 

c

 

 

 

 

èëè "e > 0 $h = h(e,a)Î[a,b): "cО(h,b) выполня-

íûé limcb

ò f (x,a)dx

 

c<b

a

 

 

c

 

 

b

 

 

 

 

 

 

 

 

 

ется условие

ò f (x,a)dx - ò f (x,a)dx

< e . По свойству аддитивности не-

 

 

 

 

 

a

 

 

a

 

собственных

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

неравенство можно пер еписать в

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

âèäå

ò f (x,a)dx

< e . Существенным здесь является то, что h выбира-

 

c

 

 

 

 

a , т.е. зависит не только от ε , íî è îò a .

åòñÿ

при каждом

Определение 13.2. Интеграл I(a) называется равномерно сходя-

щимся на

отрезке

[b,g] , åñëè "e > 0 $h = h(e)Î[a,b): "cÎ(h,b) è

 

 

æ

 

b

 

 

 

 

ö

 

 

 

 

 

 

 

 

 

 

ç

 

ò

f (x,a)dx

 

÷

 

"a Î[b,g] ç

 

 

 

< e÷ .

 

 

è

 

c

 

 

 

 

ø

 

В отличие от предыдущего в этом определении h зависит от e , но

не зависит от a , т.е. годится сразу для всех a О[b,g].

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

Теорема 13.4 (признак Вейерштрасса). Пусть существует функция g(x) ³ 0, x О[a,b), такая, что для каждого x О[a,b) и каждого a О[b,g]

b

выполняется неравенство f (x,a) £ g(x) , причем интеграл ò g(x)dx ñõî-

b

 

a

 

 

дится. Тогда I(a) = ò f (x,a)dx равномерно сходится на отрезке [b,g] .

a

 

b

 

 

¡ Возьмем произвольное e > 0 . Так как ò g(x)dx сходится, то

 

b

a

 

 

$hÎ[a,b) :"cÎ(h,b) (ò g(x)dx < e) Þ "cÎ(h,b), "a Î [b,g]

 

c

 

246

247

IV. Функции нескольких переменных. Интегралы, зависящие от п араметра

æ

 

b

 

 

b

 

b

ö

 

 

 

 

ç

 

ò

f (x,a)dx

£

ò

|f (x,a)|dx £

ò

÷

ç

 

 

 

 

g(x)dx < e÷

è

 

c

 

 

c

 

c

ø

(при доказательстве были использованы замечание к теорем е 11.3 и свойство 5. из разд. 11.2). x

13.3. Свойства несобственных интегралов, зависящих от параметра

Теорема 13.5. Пусть функция f (x,a) непрерывна на прямоуголь-

b

нике P, и интеграл I(a) = ò f (x,a)dx равномерно сходится на отрезке

a

[b,g] . Тогда функция I(a) непрерывна на отрезке [b,g] .

¡ Пусть a0 О[b,g]. Надо доказать, что lim I(a) = I(a0) , ò.å. ÷òî

α→α0

"e > 0 $d = d(e) > 0:"a, a - a0 < d ( I(a) - I(a0) < e). Возьмем произ-

вольное e > 0. Так как интеграл I(a) равномерно сходится на отрезке

[b,g] , òî $hÎ[a,b): "cÎ(h,b)

è "a Î[b,g]

æ

 

b

 

 

 

 

 

 

 

 

e ö

 

 

 

 

 

 

 

 

 

 

 

çç

 

ò f (x,a)dx

<

 

 

÷ч. Фиксиру-

 

3

ем одно такое c.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

c

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как по теореме 13.1 функция ò f (x,a)dx непрерывна на отрез-

[

 

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, а значит, и в точке

 

 

 

 

 

 

$d >

0:

 

,

a - a0

< d

 

 

 

 

 

 

 

 

êå

b,g

 

 

 

0 , òî

 

"a

 

 

 

 

 

 

 

 

 

 

 

æ

 

c

f (x,a)dx

-

c

f (x,a

 

 

)dx

 

=

 

c

 

 

 

 

 

 

 

 

 

 

 

 

)ùdx

 

< e

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

ò

ò

 

 

 

 

 

é f (x,a) - f (x,a

 

÷ Þ

 

 

 

 

ç

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

ò ë

 

 

 

 

 

 

 

 

 

 

0

û

 

 

 

 

3

÷

 

 

 

 

è

 

a

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

"a a - a < d I(a) - I(a ) =

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

é f (x,a) - f (x,a )ùdx

 

 

 

 

 

 

 

 

,

 

0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

ò ë

 

 

 

 

 

 

 

0 û

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

0

 

û

 

 

 

 

b

 

 

 

 

 

 

 

 

b

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

ò ë

 

 

 

 

 

 

 

 

 

 

 

+

ò

f (x,a)dx -

ò

 

 

 

 

 

 

 

£

 

 

 

 

 

 

 

 

é f (x,a) - f (x,a

 

 

)ùdx

 

 

f (x,a

)dx

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

< e

+ e

+ e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£

 

é f (x,a) - f (x,a

 

)ùdx

+

 

ò

f (x,a)dx

+

ò

 

f (x,a

 

)dx

= e. x

 

 

ò ë

 

 

 

 

 

 

 

 

 

0

û

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

3

3

3

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Теорема 13.6. Пусть функция f (x,a) непрерывна на прямоуголь-

 

 

b

 

 

 

нике P и интеграл I(a) = ò f (x,a)dx

равномерно сходится на отрезке

[b,g] . Тогда

 

a

 

 

 

 

 

 

 

 

γ

γ

éb

ù

b éγ

ù

òI(a)da = ò êò

f (x,a)dxú da = ò êò f (x,a)daúdx.

β

β

ê

ú

êβ

ú

 

 

a ë

û

 

 

ëa

û

γ éb

 

ù

 

 

¡ Интеграл ò ê

ò f (x,a)dxъda в левой части последней формулы

β ê

 

 

ú

 

 

ëa

 

û

 

 

существует, так как по теореме 13.5 функция I(a) непрерывна на [b,g] .

Нам надо доказать, что существует и равен интегралу слева следующий интеграл:

 

 

b é

γ

ù

 

 

c é

γ

ù

 

 

ò ê

ò f (x,a)daúdx c=b

ò ê

ò f (x,a)daúdx,

 

ê

β

ú

c<b

ê

β

ú

 

 

 

a ë

 

û

a ë

 

û

 

ò.å. "e > 0

$h = h(e)Î[a,b): "c Î(h,b)

 

 

 

 

 

 

 

c é

γ

ù

γ

éb

 

ù

 

< e .

 

 

 

 

 

 

ò ê

ò

f (x,a)daúdx

- ò

ê

ò f

(x,a)dxú da

 

 

 

ê

β

ú

β

ê

 

 

ú

 

 

 

 

a ë

 

 

û

 

 

 

 

 

 

 

 

 

 

ëa

 

û

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

Возьмем произвольное e > 0 . Так как ò f (x,a)dx равномерно схо-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

дится на [b,g] , òî $hÎ[a,b) :"cÎ(h,b),"a Î[b,g]

 

 

 

 

æ

 

b

 

 

 

 

 

e

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

ò

f (x,a)dx

<

 

 

 

 

÷ Þ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

g - b ÷

 

 

 

 

è

 

c

 

 

 

 

 

 

ø

 

 

 

"cÎ(h,b) по теореме 13.2 имеем:

 

 

 

 

 

 

 

 

 

 

 

c é

γ

 

 

ù

γ éb

 

ù

 

=

 

 

 

 

 

 

ò ê

ò f (x,a)daú dx - ò ê

ò f (x,a)dxú da

 

 

ê

β

 

 

ú

β

 

 

 

ú

 

 

 

 

a ë

 

 

 

û

ê

 

 

 

 

 

 

 

 

 

 

 

ëa

 

û

 

 

 

=

 

γ

éc

 

 

ù

 

γ

éb

 

ù

 

 

=

 

 

 

 

 

 

 

 

ò

êò f (x,a)dxúda - ò

êò

f

(x,a)dxúda

 

 

 

β

ê

 

 

ú

 

β

ê

 

ú

 

 

 

 

 

 

 

ëa

 

 

û

 

 

ëa

 

û

 

 

 

248

249

IV. Функции нескольких переменных. Интегралы, зависящие от п араметра

 

γ é

b

 

 

 

ù

 

 

 

 

 

γ

 

b

 

 

 

 

 

 

 

 

 

e

 

 

(g - b) = e. x

 

 

 

 

 

=

ò ê

-ò

f (x,a)dxúda

£

ò

ò f (x,a)dx

da <

 

 

 

 

 

 

 

 

g - b

 

 

 

 

 

 

β ê

c

 

 

 

ú

 

 

 

 

 

β

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ë

 

 

 

û

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

è

¶f (x,a)

непрерывны на

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

¶a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[

 

]

 

прямоугольнике P, интеграл I(α) =

ò

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x,α)dx сходится при α β,γ

 

, à

 

 

 

= b f

(x,α)

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

интеграл Φ α

 

равномерно сходится на отрезке

β,

γ

]

. Òîã-

(

)

 

ò

∂α

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[

 

 

 

]

 

 

[

 

 

 

 

 

 

I(

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

да функция

 

 

α a дифференцируема на отрезке β,

γ

 

è

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

b

 

 

 

 

 

 

b f (x,α)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I (α) =

 

 

 

ò f (x,α)dx = ò

 

 

 

 

dx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dα

 

∂α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¡ Надо доказать, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I ¢( )

lim

1

 

éb

f

(x,

 

 

 

)dx

 

b f (x,

 

 

)dx

ù

 

 

 

 

 

 

 

 

 

 

 

 

 

ê

 

 

 

 

-

 

 

ú =

 

 

 

 

 

 

 

 

 

 

 

a =

 

 

 

 

 

 

 

 

 

a + Da

ò

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

Δα→0 Da

êò

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ú

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ëa

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

û

 

 

 

 

 

 

 

 

 

b

 

f (x,α + Δα) f (x,α)

 

 

òò.54..45

 

 

 

b

 

f (x,α + θΔα)

 

 

 

 

 

 

= Δα→lim0 ò

 

 

 

 

 

 

 

 

 

 

 

dx

= Δα→lim0 ò

 

 

 

 

 

 

 

 

 

dx,

 

 

 

 

 

 

 

Δα

 

 

 

 

 

 

 

 

 

 

 

 

 

∂α

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ãäå 0 < θ < 1,

Δα, Δα < δ

действительно равен òb f (x,α)dx, ò. å."e > 0 δ = δ(ε) > 0: a ∂α

æ

b

¶f (x,a + qDa)

b

¶f (x,a)

 

 

ö

ç

ò

dx –ò

 

< e

÷

ç

¶a

¶a

dx

÷ .

è

a

 

a

 

 

 

ø

 

 

 

 

b

f (x,α)

 

 

 

 

 

 

 

 

Возьмем произвольное e > 0. Так как ò

 

 

 

dx равномерно схо-

∂α

 

 

 

 

 

 

a

 

 

 

æ

 

b

 

 

 

 

ö

 

 

]

 

η a,b) : c (η,b),α β,γ

]

 

 

¶f (x,a)

 

 

e

дится на

β,γ

, òî

 

ç

 

ò

dx

<

÷ .

 

 

¶a

 

 

[

 

[

[

 

ç

 

 

3

÷

 

 

 

 

 

 

 

 

è

 

c

 

 

 

 

ø

Фиксируем такое c. Так же как при доказательстве теоремы 13.3,

δ >

0:

Δα

Δα < δ è x a,c

]

будет выполняться неравенство

 

,|

|

 

 

[

 

 

 

 

 

 

 

 

 

 

æ

 

¶f (x,a + qDa)

-

¶f (x,a)

 

<

e

ö

Þ

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

 

 

÷

 

 

 

¶a

 

¶a

 

 

 

 

è

 

 

 

 

 

 

3(c - a) ø

 

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

Δα, Δα < δ будет справедливо, что

 

 

 

 

 

 

 

 

 

 

 

b

 

¶f (x,a + qDa)

 

 

 

 

 

 

 

 

b

¶f

(x,a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ò

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

- ò

 

 

 

 

 

 

 

dx

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¶a

 

 

 

 

 

 

 

 

 

 

 

 

¶a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c ¶f (x,a + qDa)

 

 

 

 

 

 

 

c

 

¶f (x,a)

 

 

 

 

 

 

 

 

b ¶f

(x,a + qDa)

 

 

 

 

 

 

b

 

¶f (x,a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

ò

 

 

 

 

 

 

dx - ò

 

 

 

 

 

 

 

 

 

dx +

ò

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

-

ò

 

 

 

 

 

dx

£

 

 

 

 

¶a

 

 

 

 

¶a

 

 

 

 

 

 

 

 

 

 

 

 

 

¶a

 

 

 

 

 

 

 

 

¶a

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c é¶f (x,a + qDa)

 

 

 

¶f

(x,a)ù

 

 

 

 

 

 

 

 

 

 

 

b

 

¶f (x,a + qDa)

 

 

 

 

 

 

 

 

b

¶f (x,a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£

ò ê

 

 

 

 

 

 

 

 

-

 

 

 

 

 

 

 

údx

+

 

 

ò

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

+

ò

 

 

 

 

 

 

dx

 

£

 

¶a

 

 

 

 

 

 

 

 

 

¶a

 

 

 

 

 

 

 

 

¶a

 

 

 

 

 

 

¶a

 

 

 

 

 

 

a ë

 

 

 

 

 

 

 

 

 

 

 

 

û

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

¶f (x,a + qDa)

 

 

 

¶f

(x,a)

 

 

 

 

 

 

 

 

 

 

b ¶f (x,a + qDa)

 

 

 

 

 

 

 

 

b

¶f (x,a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£ ò

 

 

¶a

 

 

 

 

 

 

-

 

 

 

 

¶a

 

 

dx +

ò

 

 

 

 

 

 

¶a

 

 

 

 

dx

+

ò

 

 

 

¶a

 

 

dx

<

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

ε

 

 

ε

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<

 

 

 

 

 

(c a)

+

 

 

+

 

 

= ε. x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3(c a)

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+∞ 1e−αx

 

 

 

 

 

 

 

 

 

 

 

Пример. Вычислить несобственный интеграл I(α) = ò

 

 

 

 

 

dx.

 

 

 

 

 

 

 

 

 

 

xex

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ð å ø å í è å

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Этот интеграл имеет одну особенность +∞ , так как в точке 0 функция

имеет конечный предел по правилу Лопиталя:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1e−αx

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

αeαx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

 

 

 

 

 

 

= lim

 

 

 

 

lim

 

 

 

= α.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xex

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

x0 ex

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1e−αx

 

 

 

 

 

 

 

 

+∞ 1 e−αx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее имеем I(α) = ò

 

 

xex

 

 

dx + ò

 

 

xex

 

dx.

Первый из этих интегралов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

не является несобственным, т.е. равен некоторому числу, а вт орой есть раз-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ö

ность сходящегося интеграла

ò

 

 

 

dx

ç 1

 

 

£

1

 

 

= e

–x

è

ò

e

–x

dx =– e

–x

 

= e

1

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xex

 

ç xex

 

 

ex

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e–ax

 

 

1

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и интеграла ò

 

dx = ò

 

 

, который сходится при α + 1 > 0, α > −1

xex

xe(a + 1)x

 

11

æ

1

 

1

 

 

 

 

ö

 

ç

è

ò

e

(α+1)x

dx

÷

α

ç

 

 

÷

 

 

ç

xe(α+1)x

 

e(α+1)x

 

сходится – см. ниже÷

, значит, I( ) сходится

è

 

 

 

 

 

ø

 

 

 

 

 

 

1

 

 

 

 

 

ïðè α > −1:

250

251

IV. Функции нескольких переменных. Интегралы, зависящие от п араметра

òò..1313..77+∞

æ

1- e

−αx

¢

+∞

xe

−αx

 

 

 

+∞

1

 

 

+∞

 

1

 

 

I ¢(a) = ò

 

ö

dx = ò

 

 

dx

=

ò e(α+1)xdx = -

 

e(α+1)x

 

=

 

Þ

ç

 

x

÷

 

 

x

 

 

 

 

 

xe

xe

a +1

0

a +1

0

è

 

ø α

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

I(a) = ò

da

 

 

= ln(a +1) + c,

 

 

 

 

 

 

 

 

 

 

 

 

 

a +1

 

 

 

 

 

 

 

 

 

à òàê êàê I(0) = 0 = ln1+ c, òî c = 0 è I(a) = ln(a +1).

Для правомерности применения теоремы 13.7 достаточно прове рить

+∞

 

равномерную сходимость интеграла ò e(α+1)xdx ïðè a ³ -1+ d, ãäå d ³ 0 –

0

 

произвольное фиксированное число. Тогда равенство I(a) = ln(a +1)

будет

верно при a ³ -1+ d, а значит, и при a > -1 . Эта равномерная сходимость следует из признака Вейерштрасса 13.4: e(α+1)x £ e−δx , и сходится интеграл

+∞

1

 

+∞

1

ò

 

 

e−δxdx = - d e−δx

 

0

= d.

0

 

 

 

 

13.4. Гамма-функция

Существует ряд функций, которые задаются (сходящимися) не собственными интегралами. В этом разделе мы рассмотрим одну из таких функций.

Определение 13.3. Гамма-функцией называется функция Γ(α), которая задается следующим несобственным интегралом:

+∞

 

Γ(α) = ò xα−1exdx.

(13.2)

0

 

Данный интеграл является несобственным. У него две «особе нности»: бесконечный предел ( +∞ ) и разрыв подынтегральной функции при х = 0 (если a -1< 0 Ю a < 1). Докажем, что этот интеграл сходится при всех a > 0 (т.е. формула (13.2) определяет гамма-функцию при a > 0 ). По определению интеграла с несколькими «особенностями» (см. разд. 11.1), его надо представить в виде суммы интегралов с одной «о собенностью» в каждом:

1+

Ã(a ) = òxα – 1e – x dx+ òxα – 1e– x dx.

0

1

1442443 1442443

I 1 (α )

I 2 (α )

252

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

1.Согласно теореме сравнения в предельной форме 11.2 сравним

I1(a) ñ

1 xα−1dx =

1

dx

:

lim

xα−1ex

= lim ex = 1 è

1 xα−1dx =

1

dx

ò x1– a

xα−1

ò x1−α

 

ò

 

x→+0

x→+0

ò

 

0

0

 

 

 

 

 

0

0

 

сходится при 1 – a < 1, a > 0. Следовательно I1(a) сходится при a > 0 (и

расходится при a £ 0).

2. Докажем, что I2(a) сходится при всех a (фактически это следует из того, что при x ® +¥ величина e –x стремится к нулю быстрее, чем

x в любой степени). Сравним, например, I2(a) со сходящимся интег-

+∞ dx

 

lim

xα−1ex

= lim

xα+1

= 0

 

 

ралом ò x2

:

1

 

ex

(åñëè a + 1

£ 0, то последнее

x→+∞

 

x→+∞

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

равенство очевидно; если же a + 1 > 0, то нужно достаточное количество раз применить правило Лопиталя (см. теорему 5.5), при этом зна менатель не меняется, а из числителя x в конце концов «уйдет», что и даст нужный нам результат). Cледовательно, по теореме сравнения в предельной форме 11.2 (случай K = 0), I2(α) сходится при всех a.

3. То есть интеграл (13.2) сходится при a > 0 ирасходитсяпри a £ 0.

Нахождение значений гамма-функции ( a > 0 )

1. Формула приведения имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Γ(α + 1) = αΓ(α).

 

 

 

(13.3)

¡ Применяя формулу интегрирования по частям, имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ã(a ) = òxα – 1e – x dx Þ Ã(a+1) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

=

ò

xα e– x dx = –xα e– x

 

++

ò

e– x axα – 1dx = aÃ(a ).

 

 

 

 

 

 

 

 

 

 

 

{123

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

0 u

dv

 

 

14243 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0

 

 

1442443

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

αà (α )

 

 

 

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u = xα

 

 

 

 

dv = exdx

 

 

 

 

 

 

 

 

du = axα−1dx v = -e x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

α

 

 

 

 

 

 

 

α −

 

 

xα

axα−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

e

 

x

 

 

 

= 0;x

 

e

 

x

 

 

= lim

x e

 

x =

lim

 

= lim

 

=

 

 

 

 

 

 

x=0

 

 

 

 

 

 

x=+∞

x→+∞

 

 

 

x→+∞ ex

x→+∞

ex

253

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

IV. Функции нескольких переменных. Интегралы, зависящие от п араметра

= lim

a(a -1)xα−2

= ... = 0. x

 

ex

 

x→∞

 

 

 

 

2. Гамма-функция от 1:

 

 

 

 

 

G(1) = 1.

(13.4)

 

+∞

 

 

0+∞ = - (xlim→+∞ ex - 1) = -(0 - 1) = 1.

 

¡ G(1) = ò exdx = -ex

 

x

 

 

 

 

0

 

 

 

 

3. Гамма-функция натурального аргумента:

 

 

 

 

 

G(n) = (n - 1)!

(13.5)

¡Пусть a = n – натуральное число, тогда по формуле приведения

Ã(n ) = (n – 1)Ã(n – 1) = (n – 1)(n – 2 )Ã(n – 2 ) =

=(n – 1)(n – 2 )(n – 3)Ã(n – 3) = ... =

=(n – 1)(n – 2 )(n – 3)× ... ×1×Ã(1) Þ Ã(n ) = Ã(n – 1)! x

{

1

В частности, G(2) = 1! = 1;G(3) = 2! = 2;G(4) = 3! = 6;G(5) = 4! = 24.

1

4. Гамма-функция от 2

 

 

 

 

 

 

 

 

 

 

æ 1

ö

= p.

 

(13.6)

 

 

 

 

 

 

 

 

 

 

G ç

2

÷

 

 

 

 

 

 

 

 

 

 

 

è

ø

 

 

 

æ

1 ö

=

+∞

x

1

x

dx =

+∞ ex

 

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

2

, dx = 2tdt,

¡ G ç

2

÷

ò

2e

 

ò

dx.

 

è

ø

 

0

 

 

 

 

0

x

 

 

 

 

 

x = t. Тогда

æ 1 ö

 

+e–t 2

 

+

–t 2

è

2

ø

=

ò

t

×2tdt = 2

ò

e

dt = p.

Ãç

 

÷

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

14243

π

2

Замечание. При доказательстве было использовано равенство

+∞

 

2

 

p

 

ò e

 

 

t

 

dt =

 

.

(13.7)

 

2

0

 

 

 

 

 

 

Это так называемый интеграл Эйлера–Пуассона; вычислить е го в данный момент достаточно сложно. Формула (13.7) будет доказана в гл. 19 , где она

254

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

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

5. Гамма-функция полуцелого аргумента. Введем следующие обозначения:

 

(2n)!! = 2× 4 ×6×...×(2n);

 

 

(2n -1)!! = 1×3×5×...×(2n - 1).

 

 

Тогда справедлива следующая формула:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

1 ö

 

 

 

(2n - 1)!!

 

 

 

 

 

 

 

 

 

 

 

 

 

(13.8)

 

 

 

 

 

 

 

G ç n

+

 

 

÷

 

=

 

 

 

 

 

 

 

 

 

 

p, nÎ Z .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

ø

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

1 ö

(13.3)

æ

 

 

1 ö

 

 

æ

 

 

 

1

ö

 

 

æ

 

 

1

öæ

 

 

æ

 

 

¡

Gçn +

 

÷

 

=

 

çn

-

 

÷ G

çn -

 

 

 

÷

 

=

çn -

 

 

֍n -

 

 

÷ G

çn -

 

÷

= ... =

 

 

 

2

2

 

 

2

2

 

 

è

2 ø

 

 

 

è

 

 

ø è

 

 

 

 

ø

è

 

 

øè

 

ø è

2 ø

 

 

æ

1 öæ

 

3

öæ

 

 

5

ö

 

 

 

1

×

 

æ

 

1 ö

 

(2n – 1)(2n – 3)(2n – 5 )× ... ×1

p =

= ç n –

÷ç n –

 

÷ç n –

 

÷× ...×

 

 

Ãç

 

÷

=

 

 

 

 

 

 

 

 

2 n

 

 

 

 

 

è

2 øè

 

2

øè

 

 

2

ø

 

 

 

2

 

 

è

2 ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1444442444443 123

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n сомножителей

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

(2n - 1)!!

 

 

 

p. x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Примеры. Вычислить значения G çæ 7

 

÷ö

 

è Gçæ 9

÷ö .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è 2

 

ø

 

 

è 2

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ð å ø å í è å

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

7

ö

 

 

æ

 

 

1

ö

n=3

 

5!!

 

 

 

 

 

1×3×

5

 

15

 

 

 

 

 

 

 

 

 

à ç

÷

= G ç

3+

2

÷

=

 

3

 

 

p =

 

8

 

p =

 

8

p;

 

 

 

 

 

 

 

 

è

2

ø

 

 

è

 

 

ø

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

9

ö

 

 

æ

 

 

1

ö n=4

7!!

 

 

 

 

 

1×3×5× 7

 

 

105

 

 

 

 

 

 

 

à ç ÷

= G ç 4

+

2

÷

 

=

2

4

 

p =

 

 

16

 

p =

 

16

 

p.

 

 

 

 

 

 

 

è

2

ø

 

 

è

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6. Гамма-функция отрицательного аргумента.

До сих пор мы определяли гамма-функцию для положительных a

+∞

(лишь для них сходится ò xα−1exdx ). Формула приведения

0

G(a + 1) = aG(a) может служить определением гамма-функции отрицательного аргумента. Из этой формулы

G(a) =

G(a + 1)

.

(13.9)

 

 

a

 

 

 

 

255

IV. Функции нескольких переменных. Интегралы, зависящие от п араметра

Теперь, зная G(a + 1), находим G(a). Если a О(-1,0) , то a + 1О(0,1), а на этом интервале гамма-функция уже известна.

Таким же методом определим гамма-функцию для a О(-1,0). Точ- но так же, зная гамма-функцию при a О(-1,0), по формуле (13.9) найдем гамма-функцию для a О(-2,-1) и т.д. В итоге мы определим гам- ма-функцию во всех отрицательных нецелых точках.

Получим также формулу для гамма-функции в отрицательных

полуцелых точках:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

æ

 

 

 

 

 

 

3 ö

 

 

 

 

 

æ

-n +

5

ö

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

1 ö

 

 

 

ç -n +

 

÷

 

 

 

Gç

2

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

2 ø

 

 

 

 

 

è

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

à ç -n +

 

 

÷

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ...=

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

1

 

 

æ

 

 

 

1 öæ

 

 

 

3

ö

 

 

 

 

 

 

 

è

 

 

 

ø

 

 

 

 

 

 

-n +

 

 

 

 

-n +

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

ç

 

 

֍ -n

2

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

2 øè

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

æ 1

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Gç

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p × 2

n

 

 

 

 

 

=

 

 

 

 

 

 

è 2

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

æ

 

 

1 öæ

-n +

3

ö

 

 

 

æ

-

1 ö

 

(-2n +1)(-2n + 3)...(-2n + (2n -1))

 

ç -n +

2

֍

2

÷

...ç

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

øè

 

 

 

 

ø

 

 

 

è

 

 

2 ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

2n p

 

 

 

 

 

 

 

 

=

 

(-1)n 2n p

Þ

 

 

 

 

 

 

 

 

(-1)n(2n -1)(2n - 3)×...×1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2n -1)!!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

 

1

ö

 

 

(-1)n2n p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Gç

-n +

 

 

÷

=

 

 

 

 

 

 

 

,

nÎZ .

(13.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

(2n -1)!!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

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

точках.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ð å ø å í è å

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

1

ön=1

(-1)2

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.

Gç

-

 

÷

=

 

 

 

 

 

 

 

 

= -2

 

p;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

(2 -1)!!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

-

3

ön=2

22 p

=

4

 

 

p

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Gç

2

÷

=

 

3!!

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

5

ön=3

-23 p

 

 

 

 

8 p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.

Gç

-

 

÷

=

 

 

 

 

 

 

 

= -

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

5!!

 

 

 

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

7

ön=4

24 p

 

 

 

16

 

p

 

 

 

 

 

16

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.

Gç

-

 

÷

=

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

7!!

 

 

7

×5

×3×1

 

105

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

График гамма-функции

Можно показать, что гамма-функция G(a) непрерывна при a > 0 и a < 0 , a ¹ -n , где nО Z.

Исходя из формулы приведения G(a) = G(a + 1) , имеем a

G(+0) = lim G(a) = lim

G(a +1)

=

G(1)

 

=

1

 

 

= +¥;

a

 

lim a

 

lim a

α→+0

α→+0

 

 

 

 

 

 

 

 

 

 

α→+0

 

 

 

α→+0

 

 

G(-0) = lim G(a) = lim

G(a +1)

=

G(1)

 

=

1

 

 

= -¥;

a

 

lim a

 

lim a

α→−0

α→−0

 

 

 

 

 

 

 

 

 

 

α→−0

 

 

 

α→−0

 

 

G(-1+ 0) = lim

G(a) = lim

G(a + 1)

=

G(+0)

= -¥;

 

a

-1

 

α→−1+0

 

α→−1+0

 

 

 

 

 

 

G(-1- 0) = lim G(a) =

lim

G(a +1)

=

G(-0)

= +¥ è ò.ä.

 

 

α→−10

α→−10

a

 

 

 

 

-1

 

 

 

 

Теперь, учитывая эти пределы и посчитанные значения гамма -фун- кции в некоторых точках, можем изобразить график гамма-фу нкции (рис. 76). При этом отметим, что мы не искали точки экстремума ф ункции (что совсем не так легко).

y

4 3 2 1 0

1

2

3 4 x

Ðèñ. 76

256

257

c = y0ekx0 è y = y0ekx0ekx = y0
âèäà y = cx
xc2 , - xc2 = - cxx .

5. Дифференциальные уравнения

V

ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ

«

14. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА И УРАВНЕНИЯ, ДОПУСКАЮЩИЕ ПОНИЖЕНИЕ ПОРЯДКА

14.1. Примеры задач, приводящих к дифференциальным уравнениям

Задача 1. Найти все кривые, обладающие

y

 

следующим свойством: если в любой точке

 

кривой провести касательную, то отрезок ка-

A

 

сательной, заключенный между осями коор-

 

 

 

динат, будет делиться в точке касания попо-

 

 

ëàì (ðèñ. 77).

M

 

Решение

B

 

 

x

Пусть y = f (x) искомая кривая. Обозна-

 

чая через X и Y координаты касательной к

 

Ðèñ. 77

этой кривой, запишем уравнение этой каса-

 

 

тельной в точке (x, y): Y - y = y¢(X - x).

Найдем точки пересечения касательной с осями координат:

A : X = 0 ÞY = y - y¢x; B :Y = 0 Þ X - x = - yy¢ , X = x - yy¢.

Так как точка касания M должна являться серединой отрезка [A, B], то ее координаты должны быть средними арифметическими со ответ-

 

 

 

 

 

 

 

 

 

0 + x

y

 

 

 

 

ствующих координат точек A и B, т.е. x =

y

, èëè

-

y

 

= x,

 

 

2

 

 

 

 

y

 

y yx + 0

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

= − x

è y =

2

, èëè -y¢x = y, y¢ = - x .

 

 

 

 

 

 

14. Дифференциальные уравнения первого порядка...

Таким образом, наша кривая должна удовлетворять соотноше нию y′ = − xy . Такого типа соотношения, куда входят производные искомой

функции, называются дифференциальными уравнениями (или, более строго, обыкновенными дифференциальными уравнениями).

Убедимся, что нашему соотношению удовлетворяют все функц ии , где c – произвольная постоянная: y′ =

Ниже будет показано, что других функций, удовлетворяющих данно-

му уравнению, нет, т.е. функция y = xc дает все решения уравнения y¢ = - xy .

Задача 2. Известно, что скорость распада радия прямо пропорциональна имеющемуся количеству радия. Пусть в момент времен иx = x0 имелся y0 г радия. Найти y (x) – количество радия в любой момент времени x > x0 .

Решение

По условию задачи y¢(x) = ky(x) ( k < 0 – коэффициент пропорциональности) и y(x0) = y0 . Это есть дифференциальное уравнение и так называемое начальное условие, т.е. задается значение функ ции при некотором значении аргумента.

Убедимся, что нашему уравнению удовлетворяют все функции вида y = cekx : y¢ = ckekx , ckekx = kcekx (с – произвольная постоянная). Ниже будет показано, что других решений у данного уравнен ия нет.

Теперь используем начальное условие y(x0) = y0 : cekx0 = y0 , ò.å. k(xx0) – единственное решение нашей

задачи.

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

(в примере 1, если y(x0) = y0

, то из формулы y =

c

ïðè x = x0 имеем

x

 

 

=

c

, c = x

 

 

 

 

 

y

 

y

 

).

 

 

 

0

x

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

258

259

5.Дифференциальные уравнения

14.2.Дифференциальные уравнения произвольного

èпервого порядков

Определение 14.1. Дифференциальным уравнением называется соотношение, связывающее независимую переменную x, искомую функцию y = f (x) и производные этой функции y¢,y¢¢,...,y(n) . В общем слу-

чае это соотношение можно записать в виде

 

F (x,y,y¢,y¢¢,...,y(n) ) = 0,

(14.1)

где F – некоторая функция (n + 2)x переменных.

Порядком дифференциального уравнения называется наивысший порядок производной, входящей в это уравнение (т.е. формула (14.1) задает дифференциальное уравнение n-го порядка).

Решением (или частным решением) дифференциального уравн ения(на некотором множестве) называется всякая функция y = j(x) , которая при подстановке в уравнение обращает его в тождество (на э том множестве).

Если решение уравнения задано в неявной форме Ô (x,y)= 0, то такое равенство называют интегралом дифференциального уравнения.

График решения дифференциального уравнения называется интегральной кривой.

Решить дифференциальное уравнение – значит найти все его решения.

Определение 14.2. Дифференциальным уравнением первого порядка называется уравнение

F (x,y,y) = 0,

(14.2)

где F – некоторая функция трех переменных. Если из этого уравнения можно выразить y, то оно примет вид

y′ = f (x,y),

(14.3)

где f (x,y) – некоторая функция двух переменных. Уравнение (14.3) называется уравнением, разрешенным относительно произв одной. Мы будем в основном рассматривать именно такие уравнения .

Определение 14.3. Задачей Коши для уравнения (14.3) называется задача

14. Дифференциальные уравнения первого порядка...

ìy¢ = f (x,y);

 

íy(x

) = y

0

.

(14.4)

î

0

 

 

 

Здесь x0 è y0 – некоторые числа.

Требуется найти решение дифференциального уравнения (14.3) , удовлетворяющее начальному условию y(x0) = y0 .

Теорема 14.1 (существования и единственности решения задачи

Коши). Пусть функция f (x,y) и ее частная производная f (x,y) íå-

y

прерывны в некоторой области D на плоскости 0xy и точка (x0,y0)ОD . Тогда задача Коши (14.4) имеет решение, и притом единственное . Это решение определено в некоторой окрестности точки x0.

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

Геометрический смысл теоремы заключается в том, что через точку (x0, y0) проходит единственная интегральная кривая.

Определение 14.4. Пусть в области D выполняются условия те-

оремы 14.1. Функция

 

f = j(x,c),

(14.5)

где c – постоянная, называется общим решением уравнения первого порядка (14.3) в некоторой окрестности точки x0 U(x0) , åñëè:

1)ïðè "x ÎU (x0) и "cОC , где C – некоторое множество (в простых случаях c вообще любое) функция (14.5) является решением уравнения (14.3);

2)любое решение уравнения (14.3), график которого лежит в обла -

ñòè D ( x U (x0)), получается из формулы (14.5) при некотором значе- нии c ОC .

Покажем, что условие 2 в этом определении можно заменить на ус-

ловие 2¢. Для любого начального условия y(x0) = y0 , ãäå (x0,y0) D , существует значение постоянной c0 ОC , при котором функция (14.5) удовлетворяет этому начальному условию: j(x0,c0) = y0 .

¡ Пусть функция (14.5) удовлетворяет условию 2. Зададим произ-

вольное начальное условие y(x0) = y0 , ãäå (x0,y0) D . Тогда по теореме 14.1 существует единственное решение уравнения (14.3), удовле т- воряющее данному начальному условию. Согласно условию 2, э то ре-

260

261