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

лекции вышмат

.pdf
Скачиваний:
93
Добавлен:
01.06.2024
Размер:
5.41 Mб
Скачать

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

xds ( L )

, где (L) – часть окружности

x

2

+ y

2

 

 

=

1

,

находящаяся в 1-ой четверти системы

координат.

Решение. Перепишем уравнение окружности в виде

x = cos t, y = sin t, z

=0, t 0,

2

. Тогда по формуле (6)

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

xds = cos t

(sin t)

2

+ (cos t)

2

dt = cos tdt = sin t

/ 2

 

 

0

( L )

0

 

 

 

 

0

 

= 1.

130

ЛЕКЦИЯ 5

Криволинейный интеграл второго рода

Определение 1. Пусть в пространстве заданы кривая

(AB) и на ней некоторая функция

P(M ) = P(x, y, z) , где

M (x, y, z) (AB) . Разобьем

кривую

(AB)

на части точками

A

x

, y

, z

i ). На каждой части (

A A

возьмем произвольную точку

M

i (

x

, y

, z

i )

(рис. 1).

i (

i

i

 

i i+1 )

 

i

i

 

Рисунок1

Пусть

xi

=

xi+1

xi

, т.е.

xi

это проекция дуги

(A

, A

)

i

i+1

 

на ось 0x (длина этой проекции со

 

 

 

 

n1

 

 

знаком). Составим

интегральную

сумму

σ = lim P(Mi ) xi . Обозначим через

λ

 

 

 

λ0

i=0

 

 

 

 

 

 

 

 

максимальную по i

длину дуги (

A A

λ = max si . Если существует предел наших

i i+1 ):

 

 

 

i

 

 

 

интегральных сумм при λ 0 , который не зависит от разбиения (AB) на части (

A A

и

i i+1 )

от выбора точек

M

i

 

(

A A

 

 

 

i i+1 ), то этот предел

второго рода от

функции

P(M ) = P(x, y, z)

обозначается

 

P(x, y, z)dx . Таким образом,

 

 

( AB)

 

 

 

 

 

 

 

 

 

n1

i

i

 

 

P(x, y, z)dx =lim

 

 

 

 

 

P(M

) x

 

 

( AB)

λ0

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

называется криволинейным интегралом

по кривой

(AB)

по

переменной х и

 

n1

 

 

 

 

 

 

= lim

 

P(x , y

, z

) x

(1)

λ0

i

i

i

 

i .

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Q(x, y, z)dy

( AB)

 

n1

 

 

= lim Q(Mi ) yi

λ0

i=0

 

 

 

yi+1

yi

 

 

,

 

 

n

 

 

 

R(x, y, z)dz = lim R(Mi ) zi

( AB)

λ0

i=0

 

 

 

z

z

 

 

 

 

 

i+1

i

и

 

Pdx + Qdy + Rdz

( AB)

 

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

Определение 1 аналогично определению обычного определенного интеграла, и свойства криволинейного интеграла второго рода аналогичны свойствам определенного интеграла и доказываются таким же образом. Из определения 1 также следует, что, в отличие от

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

 

P(x, y, z)dx = −

P(x, y, z)dz

( AB)

( BA)

 

(при составлении интегральных

сумм для интеграла справа точки Ai и

учетом этого свойства,

P(x, y, z)dx

 

( AB)

Ai+1

меняются местами, и xi xi+1 = −(xi+1 xi ) ). С

=

 

P(x, y, z)dx + P(x, y, z)dx , независимо от

( AC )

(CB)

взаимного расположения точек A, B и C (см. рис. 2). Аналогично для интегралов по y и по z и для суммы таких интегралов.

131

Рисунок 2

Если контур L замкнут, то определяется так: пусть A, В,

криволинейный интеграл вдоль С, D L – произвольные точки,

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

любой

=

( L)

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

 

+

 

(рис. 3).

 

 

( ACB)

 

( BDA)

 

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

не изменится.

Рисунок 3

Если (L) – замкнутый контур, то под интегралом Pdx + Qdy + Rdz , при отсутствии

( L)

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

Pdx + Qdy +

( L)

Rdz

.

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

Теорема 1. Пусть кривая

(AB)

задана параметрически:

x = x(t), y = y(t),

z =

z(t)

, t α,β

и при изменении параметра

t от α до β кривая описывается в направлении от A

t = α ,

t = β

получаем

соответственно

 

точки

A и B .

Пусть

P(x, y, z),Q(x, y, z), R(x, y, z)

непрерывны вдоль нашей кривой (AB) , а функции

z(t)

непрерывно дифференцируемы на

 

α,β

. Тогда

криволинейный

 

 

P(x, y, z)dx + Q(x, y, z)dy + R(x, y, z)dz существует и

 

 

 

 

 

к B ; при функции

х (t), у (t),

интеграл

( AB)

P(x, y, z)dx + Q(x, y, z)dy + R(x, y, z)dz =

( AB)

β

= P(x(t), y(t), z(t))dx(t) + Q(x(t), y(t), z(t))dy(t) + R(x(t), y(t), z(t))dz(t) =

α

β

 

 

 

 

 

 

(2)

= P(x(t), y(t), z(t))x (t) +Q(x(t), y(t), z(t)) y (t) + R(x(t), y(t), z(t))z (t) dt

α

132

(т.е. для вычисления криволинейного интеграла мы должны всюду – и в подынтегральных функциях, и под знаками дифференциалов – заменить x, y, z на их выражения через

параметр t).

Доказательство существования интеграла мы здесь опустим и проверим лишь справедливость формулы (2). Докажем одну часть равенства (2) (остальные части доказываются аналогично):

 

P(x, y, z)dx =

( AB)

 

β P(x(t),

α

y(t), z(t))x (t)dt

.

(3)

В соответствии с определением 1, левая часть этого равенства равна

 

 

n1

i

i

i

i

 

 

 

P(x, y, z)dx =lim

 

P(x , y

, z

) x

( AB)

λ0

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1

i

 

 

= lim

 

P(x ,

λ0

i=0

 

 

 

y

, z

) (x(t

i+1

)

i

i

 

 

x(t

))

i

 

.

Согласно теореме Лагранжа, последнее равенство можно переписать в виде

 

P(x, y, z)dx

( AB)

 

 

n1

i

 

 

=lim

 

P(x ,

λ0

i=0

 

 

 

y

, z )x (t

) t

i

i

i

i

 

, где

t

i

 

– некоторое промежуточное значение

между

точек

t

i

 

M

 

и

ti

i

(x

,

 

i

 

+1 .

y

,

i

 

Так как интеграл существует, то последний предел не зависит от выбора

z

i ), поэтому для вычисления интеграла мы можем выбрать эти точки

 

наиболее удобным для нас способом: так, чтобы они соответствовали значениям

t

=

t

i

.

 

 

Тогда

 

 

n1

 

P(x, y, z)dx =lim P(x(ti ), y(ti ), z(ti ))x (ti ) ti

( AB)

λ0

i=0

 

 

 

.

Сумма в правой части этого равенства является интегральной суммой для интеграла в правой части формулы (3) и в пределе как раз дает этот интеграл. ■

Пример 1. Вычислить

 

xdx + ydy + ( x + y 1)dz , где L – отрезок прямой от А (1, 1, 1) до В (2, 3, 4).

 

 

( L )

 

x Решение. Уравнения прямой имеют вид: 2

уравнения этой прямой. При t = 0 получаем точку В. Следовательно,

1

=

 

1

 

x =

y 1

=

z

 

 

3 1

 

4

1, y = 1, z

1

=

 

1

 

 

= 1

,

t x = t + 1, y

т.е. точку А, при

= 2t t =

+ 1, z = 3t +

1 получаем x

1

=

– параметрические

2, y = 3, z = 4

, т.е.

 

1

 

 

 

 

 

 

 

xdx + ydy + ( x + y 1)dz =

(t + 1)dt + (2t + 1)d (2t + 1) + (t + 1 + 2t

+ 1 1)d (3t + 1) =

 

 

( L )

0

 

 

 

 

 

 

 

 

 

 

 

 

2

1

 

 

1

1

1

 

t

 

 

 

 

 

 

 

 

= (t + 1)dt + (2t + 1)2dt + (3t + 1)3dt = (t + 1 + 4t + 2 + 9t + 3)dt = (14t + 6)dt = 14

 

 

 

+ 6 t

1

2

 

0

0

0

0

 

0

 

 

 

 

 

 

 

= 7 + 6

=

13

.

Пример 2. Вычислить

Решение.

(L )

ydx +

xdy

, где L – путь от A(0, 0) до B(1, 1) вдоль:

y = x

и

y =

x

2

 

(рис.4).

Рисунок 4

133

y

y

=

=

 

 

 

 

 

 

 

1

 

 

 

x

: берем за параметр

x 0,1

 

ydx + xdy =

 

xdx + xdx =

 

 

 

 

 

 

 

( L )

 

0

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

: берем за параметр x 0,1

 

 

 

2

 

2

x

 

ydx + xdy = x

dx + xdx

 

 

 

 

 

 

( L )

 

0

 

 

 

1

 

 

 

2

1

 

x

 

 

 

 

 

2 xdx = 2

 

 

=

0

 

2

0

1

 

+ 2x

 

 

)dx

= (x

2

2

 

 

 

0

 

 

 

 

 

1

=

;

1

 

 

3

1

 

x

 

2

 

 

 

3 x

dx = 3

3

=

0

 

0

 

 

 

1

.

В обоих случаях мы получили один и тот же ответ. Встает вопрос: правилом независимости криволинейного интеграла (второго интегрирования? Приведенный ниже пример показывает, что, правила нет.

не является ли этот факт рода) от формы пути вообще говоря, такого

Пример 3. Для тех же путей, что в примере 2, вычислить

2 ydx +

xdy

.

 

 

 

 

 

1

 

 

 

1

 

2

1

 

 

 

 

 

 

 

 

 

 

x

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

y = x :

2 ydx + xdy =

2xdx + xdx = 3

xdx = 3

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( L )

 

0

 

 

 

0

2

0

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

1

 

 

y = x2

:

2 ydx + xdy = 2x2 dx + xdx2

= (2x2

+ 2x2 )dx = 4 x2 dx

 

( L )

 

 

0

 

 

0

 

 

 

 

0

 

 

( L )

;

= 4 x3 1

3 0

=4 .

3

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

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

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

Далее будем рассматривать криволинейный интеграл

 

Pdx + Qdy +

( AB)

 

Rdz

, в котором

функции P(x, y, z), Q(x, y, z), R(x, y, z) непрерывны в некоторой области (Т) трехмерного пространства.

Теорема 2. Для того чтобы

 

Pdx + Qdy + Rdz

при любых кривых (АВ) (Т) не зависел от

( AB)

пути интегрирования (АВ), а зависел только от положения

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

 

Pdx + Qdy + Rdz = 0

 

 

( L)

 

(L) (T ) .

 

 

начальной и конечной точек А и для любого замкнутого контура

Необходимость. Пусть наш интеграл не зависит от формы пути интегрирования. Тогда для произвольного замкнутого контура (L) (T ) , изображенного на рис. 5,

Рисунок 5

134

имеем:

Pdx + Qdy + Rdz =

 

Pdx + Qdy + Rdz +

 

Pdx + Qdy + Rdz =

 

 

 

( L)

 

( ACB)

 

( BDA)

 

=

 

Pdx + Qdy + Rdz

 

Pdx + Qdy + Rdz = 0,

 

 

 

( ACB)

 

 

( ADB)

 

 

 

 

так как интеграл не зависит от пути интегрирования.

Достаточность.

теоремы

 

Pdx

 

 

( ACB)

 

В этой части теоремы нам нужно доказать, что при выполнении условий

+ Qdy + Rdz =

 

Pdx + Qdy + Rdz .

 

 

( ADB)

 

Докажем, что разность левой и правой частей этого равенства равна 0:

 

Pdx + Qdy + Rdz

 

Pdx + Qdy + Rdz =

( ACB)

 

( ADB)

 

=

 

Pdx + Qdy + Rdz +

 

Pdx + Qdy + Rdz =

 

Pdx + Qdy + Rdz

 

( ACB)

 

( BDA)

 

( ACBDA)

 

=

0

,

как интеграл по замкнутому контуру. ■

 

 

 

 

Теорема 3. Для того чтобы

Pdx + Qdy + Rdz не зависел от формы пути интегрирования

 

 

( AB)

 

 

 

 

 

для любой кривой (АВ) (Т), необходимо и достаточно, чтобы выражение

Pdx +Qdy + Rdz

являлось

в области

(Т) дифференциалом от

некоторой

функции трех переменных

u(x, y, z)

(т.е.

чтобы

существовала

такая

функция

u(x, y, z) ,

что

du(x, y, z) = P(x, y, z)dx +Q(x, y, z)dy + R(x, y, z)dz

).

Необходимость. Пусть

 

Pdx + Qdy + Rdz

 

 

 

 

( AB)

 

 

 

 

 

 

( x, y,z )

 

Тогда введем функцию

u(x, y, z) =

 

 

 

 

 

 

 

 

 

( x

, y

,z

)

 

 

0

0

0

 

не зависит от формы пути интегрирования.

Pdx + Qdy + Rdz , где

(x0 , y0 , z0 ) (T )

фиксированная

соединяющему

через точку

(x,

( y

точка и

криволинейный

x0 , y0 , z0 )

и (х, у, z). Тогда,

, z) , имеем:

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

любому

(x + x,

y

пути в (Т), , z) проходит

 

 

( x+ x, y,z )

( x, y,z )

 

xu = u(x + x, y, z) u(x, y, z) =

 

 

Pdx + Qdy + Rdz

 

 

Pdx + Qdy

 

 

( x

, y

,z )

( x

, y

, z

)

 

 

0

0

0

0

0

0

 

 

( x+ x, y,z )

 

 

 

 

 

 

 

=

 

Pdx + Qdy + Rdz

 

 

 

 

 

 

 

 

; последний интеграл берется по любому пути;

 

( x, y,z )

 

 

 

 

 

 

 

+ Rdz =

возьмем за такой

путь y =

отрезок, соединяющий точки (х, у, z) и (x + x, y, z = z – постоянные, следовательно,

y,

z)

; его

уравнения: х = t, dy = dz =

t

[x, x + x];

0. Отсюда:

( x+ x, y,z )

 

 

 

 

x+ x

 

 

теорема о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

среднем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

), где с между х и

 

Pdx + Qdy + Rdz

=

P(t, y, z)dt

 

 

=

P(c, y, z)

x + x x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( x, y,z )

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

u(x, y, z)

 

u

 

 

 

 

 

Р(x,y,z) -

 

 

 

 

 

 

 

 

 

 

 

 

непрерывна

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

= P(c, y, z)

 

 

 

 

= lim

x

 

 

= lim P(c, y, z)

=

P(x, y, z)

 

 

 

 

 

 

 

 

 

x

 

 

x

 

x

 

 

 

 

 

 

 

 

x0

 

x0

 

 

 

 

 

 

 

Аналогично

u(x, y, z)

= Q(x, y, z) и

u(x, y, z)

= R(x, y, z) . Следовательно,

 

 

 

y

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

du = ux dx + uy dy + uz dz = Pdx + Qdy + Rdz .

x + x

135

Достаточность. Пусть

Pdx +Qdy + Rdz = du

уравнениями:

x = x(t), y = y(t), z = z(t),t

 

α,β

 

 

 

 

 

 

P =

u

, Q =

u

, R =

u

и

Pdx + Qdy + Rdz =

 

x

y

z

 

 

 

 

 

( AB)

 

 

 

 

 

 

 

 

 

 

 

 

 

β

 

 

 

 

 

 

 

 

 

 

= [P(x(t), y(t), z(t))x (t) + Q(x(t), y(t), z(t)) y (t)

и кривая (АВ) задана параметрическими

. Тогда

+ R(x(t), y(t), z(t))z (t)]dt =

α

 

 

 

 

 

 

β

u

 

u

 

u

 

= [

(x(t), y(t), z(t))x (t) +

(x(t), y(t), z(t)) y (t) +

(x(t), y(t), z(t))z (t)]dt

x

y

z

α

 

 

 

 

 

 

 

 

 

По правилу дифференцирования сложной функции, это выражение равно

.

β

d

 

 

 

 

 

 

 

 

 

 

 

u(x(t), y(t), z(t)) dt =

α

dt

 

 

 

 

 

 

 

 

 

 

 

Итак, если

u

dx +

u

dy +

x

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( AB)

u (x(t), y(t), z(t))

u

dz = du

 

z

, то

 

Pdx + Qdy + Rdz

t =β

u (x(t), y(t), z

 

= u(B) u( A) = u

(t))

B

.

A

 

t =

= u(B) u( A)

.

(4)

Правая часть формулы (4) не зависит от формы пути из А в В. Этим мы и доказали независимость интеграла от формы пути интегрирования. ■

Заодно мы получили формулу для вычислений, и для примера 2 выше:

(1,1)

 

(1,1)

 

 

 

 

ydx + xdy =

 

d (xy) = (xy)

(1,1)

= 1

(0,0)

( 0 , 0 )

 

( 0 , 0 )

 

 

 

для любого пути.

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

y x

области (D) на плоскости 0xy . Для того чтобы

 

Pdx + Qdy

 

 

( AB)

 

непрерывны в некоторой

не зависел от формы пути

интегрирования для любой кривой (АВ) (D), необходимо, а в предположении односвязности области, и достаточно, чтобы при (x, y) D

P(x, y)

=

Q(x, y)

.

(5)

y

x

 

 

 

 

Напомним, что область

(D)

на плоскости называется односвязной, если вместе с каждым

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

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

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

Теорема 5. Пусть функции

P(x, y, z),Q(x, y, z), R(x, y, z)

и их частные производные

P P Q Q R R

y , z , x , z , x , y непрерывны в некоторой области (Т) трехмерного пространства.

136

Для того чтобы

 

Pdx + Qdy + Rdz

( AB)

 

не зависел от формы пути интегрирования для любой

кривой (АВ) (Т), необходимо, а в предположении односвязности области, и достаточно, чтобы при (x, y, z) (T )

P

=

Q

;

Q

=

R

;

R

=

P

.

(6)

y

x

z

y

x

z

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 4 является частным случаем теоремы 5.

Необходимость. Пусть криволинейный интеграл

 

Pdx + Qdy + Rdz

( AB)

 

не зависит от

формы пути интегрирования. Тогда по теореме 3 выражение дифференциалом от некоторой функции трех

Pdx +Qdy + Rdz = du

Pdx +Qdy + Rdz

переменных

является

u(x, y, z) :

P =

u

,Q =

u

, R =

u

 

x

y

z

 

 

 

 

производных. Отсюда

P

=

y

 

P

y

 

u

2

 

y x

=Qx

;

.

Q

 

 

u

 

 

u

 

 

u

 

 

 

2

 

 

2

 

 

2

 

 

x

=

x y

. Но

y x

=

x y

по теореме о смешанных

 

 

 

 

Аналогично

проверяются другие равенства:

Q

 

2

u

 

R

 

2

u

 

Q

 

R

 

R

 

2

u

 

P

 

2

u

 

R

=

 

;

=

 

 

=

;

=

 

;

=

 

 

z

z y

y

y z

z

y

x

x z

z

z x

x

 

 

 

 

 

 

 

 

 

 

Тем самым, мы доказали необходимость и для теоремы 4.

Достаточность условий (5) и (6) будет проверена позднее. ■

=Pz

.

137

ЛЕКЦИЯ 6

ТЕОРИЯ ПОЛЯ

Скалярное поле. Производная по направлению. Градиент

Определение 1. Пусть каждой точке М некоторой области (Т) трехмерного пространства соответствует скалярная (числовая) величина U =U (M ) . Тогда говорят, что в области (Т)

задано скалярное поле. Если положение точки М определять ее координатами в некоторой прямоугольной декартовой системе координат M (x, y, z) , то U =U (M ) =U (x, y, z) . Мы

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

Определение 2. Рассмотрим произвольную точку

M0 (x0 , y0 , z0 ) (T ) . Проведем через эту

точку направленную прямую (ось)

l .

Пусть

M (x, y, z) l – произвольная

точка

оси.

Обозначим M 0 M – расстояние от M0

до

M со знаком: «+» при смещении

M от

M0 в

сторону направления оси, « » при смещении

M от M0

в сторону,

противоположную

направлению оси. Производной скалярного поля

u = u(M )

в точке M0

по направлению l

(обозначение:

u

lu

l

(M 0 ) ) называется

 

 

 

 

 

 

 

 

 

 

 

 

 

u(M ) u(M

)

 

 

u(x, y, z) u(x

, y

, z

)

 

(M

) =

lim

 

0

 

=

lim

 

0

0

0

 

,

0

 

M0 M 0

M

M

 

 

M0M 0

M

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

(1)

если этот предел существует.

 

Пусть ось l составляет с положительными направлениями осей координат 0х, 0у, 0z

углы

α, β, γ соответственно. Обозначим M0M = t . Так как координаты любого вектора –

это длины его проекций на оси координат со знаком, то

M

0

M =

t cos α,t cos β,t cos γ

 

 

 

:

Рисунок 1

(горизонтальная прямая на рис. 1 – это прямая параллельная одной из осей координат, а отмеченный угол равен α или β , или γ соответственно).

Но координаты вектора также есть разность координат его конца и начала:

M 0 M = x x0 , y y0 , z z0

x x0

= t cos α, y y0

= t cosβ, z z0

= t cos γ

,

x = x0 +t cos α, y

Теперь u = u(x,

= y0 +t cosβ, y, z) = u(x(t),

z = z0 +t cos γ , y(t), z(t)) = φ(t

т.е. x, y, и z есть некоторые функции t.

) – некоторая функция от t , φ(0) = u(x0 , y0 , z0 )

 

u

 

 

u(x, y, z) u(x

, y

, z

)

и

 

(M0 ) =

lim

 

 

0

0

0

 

l

M

 

M

 

 

 

 

 

M0 M 0

0

 

 

 

 

 

 

 

 

 

 

 

 

= lim

φ(t) φ(0)

t

t 0

=

φ (0)

.

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

(t) = ux dxdt + uy dydt + uz dzdt = ux cos α + uy cosβ + uz cos γ .

138

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

 

u

,

u

,

u

 

будут браться в точке M0

и

 

x

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

(M

 

)

=

u

(M

 

) cos α

+

u

(M

 

) cosβ

+

u

(M

 

) cos γ

.

l

0

x

0

y

0

z

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 1. Если на множестве

T

функция

 

u(M ) = u(x, y, z)

имеет непрерывные частные

производные, то в любой точке

 

M 0 T

 

эта функция имеет производную по любому

направлению l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

(M

 

) =

u

(M

 

) cos α +

u

(M

 

) cosβ +

u

(M

 

) cos γ

,

 

(2)

l

0

x

0

y

0

z

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

α, β, γ

- углы

l

с осями координат

0x

,

0y

,

0z

соответственно.

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

Для ответа на этот вопрос сначала заметим, что вектор

e =

cos α, cosβ, cos γ

– это

 

 

 

единичный вектор (т.е. вектор длины 1), направленный вдоль оси

l

вектора как раз и равны косинусам его углов с положительными координат). Введем еще один вектор:

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

Определение

 

u

,

grad u =

x

 

 

3.

uy

Градиентом функции

u(M ) = u(x, y, z)

,

u

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

 

 

z

 

 

в точке

M0 ).

M

0

 

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

Правая часть формулы (2) равна скалярному произведению этих двух векторов; так как

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

ab = a

b cos φ

,

где

φ

угол между

a

 

 

u

= grad u e = grad u

e cos φ = grad u cos φ ,

где

φ

угол между grad u

l

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

последнее выражение будет наибольшим при cos φ =1, т.е. φ = 0 .

иb ,

иосью

то

l ;

Таким образом, производная функции в точке

M

0

 

будет наибольшей по направлению

вектора

grad u

в этой точке, и эта наибольшая производная по направлению равна

grad

u(M

0

 

)

.

Замечание. Из формулы (1) следует, что значение производной по направлению не зависит

от выбора системы координат

xyz . Тогда вектор

grad u(M

0

)

тоже не зависит от выбора

 

 

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

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

u =

наибольшую производную по направлению в точке Решение. M 0 M = 2 1, 4 2, 6 3 = 1, 2, 3 ; M

xyz

в точке

M0 .

 

 

M

=

1

+

 

 

2

 

0

 

 

 

M 0 (1, 2, 3) по направлению

2

2

+ 3

2

=

14

; из рис. 1

 

 

M 0 M к

точке M (2, 4, 6)

 

и

видно,

что cos =

 

1

 

,

 

 

 

 

 

14

 

 

cos =

2

 

, cos =

3

 

;

u

= yz,

u

= xz,

u

= xy;

u

(M

0 ) = 6,

u

(M 0 ) = 3,

u

(M

0 ) = 2;

 

 

 

 

 

 

x

y

 

 

y

 

 

 

 

 

 

14

 

 

14

 

 

 

 

z

x

 

 

z

 

139