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

Лекции по вышмату за 1 курс

.pdf
Скачиваний:
110
Добавлен:
15.03.2015
Размер:
1.53 Mб
Скачать

èëè ZZ

P (x; y; z) dydz + Q (x; y; z) dzdx + R (x; y; z) dxdy:

Равенство

ZZ

P (x; y; z) dydz + Q (x; y; z) dzdx + R (x; y; z) dxdy =

ZZ

=(P (x; y; z) cos + Q (x; y; z) cos + R (x; y; z) cos ) dS

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

94Свойства и правило вычисления поверхностного интеграла второго рода

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

2.Как уже было сказано, интеграл меняет знак, если изменить сторону поверхности.

3.Данный интеграл линеен относительно функции, т.е.

ZZ

~a +~b

~n0

dS = ZZ

(~a ~n0) dS + ZZ

~b ~n0

dS:

 

 

 

 

 

 

 

4.Интеграл аддитивен относительно области интегрирования, т.е., если поверхность представить в виде объединения двух частей 1 è 2,

RRRR RR

которые не имеют общих внутренних точек, то

=

+ :

 

1

2

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

z = ' (x; y), òî

ZZ

P (x; y; z) dydz + Q (x; y; z) dzdx + R (x; y; z) dxdy =

 

 

 

 

= ZZ

(P (x; y; z) cos + Q (x; y; z) cos + R (x; y; z) cos ) dS =

:

= ZZ

(P (x; y; ' (x; y)) cos + Q (x; y; ' (x; y)) cos + R (x; y; ' (x; y)) cos ) jcos j

 

 

 

dxdy

 

Однако, можно поступить и иначе. Разобьем интеграл

ZZ

(P (x; y; z) cos + Q (x; y; z) cos + R (x; y; z) cos ) dS

291

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

каждого из следующих уравнений z = ' (x; y), или y = (x; z), или

x = (y; z). Тогда

 

Q (' (y; z) ; y; z) cos cos = ZZ

 

ZZ Q (x; y; z) cos dS = ZZ

Q (' (y; z) ; y; z) dydz;

 

 

 

 

dydz

 

 

 

 

DOY Z

 

j

 

 

j

DOY Z

 

ZZ Q (x; y; z) cos dS = ZZ

Q (x; (x; z) ; z) cos cos

 

= ZZ

Q (x; (x; z) ; z) dzdx

 

 

 

 

dzdx

 

 

 

 

DOXZ

 

j

 

 

j

DOXZ

 

è

 

Q (x; y; ' (x; y)) cos cos

 

 

= ZZ

Q (x; y; ' (x; y)) dxdz;

ZZ Q (x; y; z) cos dS = ZZ

 

 

 

 

 

dxdy

 

 

 

 

 

DOXY

 

j

 

j

DOXY

 

где выбирается знак +, если соответствующий косинус положитель-

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

цателен, т.е. угол между нормалью и координатной осью тупой.

Пример 264 Вычислить RR y2 z2 dxdy, где - верхняя сторона ци-

линдра z =

 

a2 x2, отсеченная плоскостями y = 0 и y = b:

 

 

 

 

 

 

 

 

 

 

 

 

Решение Проекцией поверхности на плоскость ОХУ будет прямоуголь-

 

 

 

 

 

ник D. Так как нормаль к выбранной стороне поверхности образует острый

 

 

 

 

 

угол с осью OZ, то в формуле берется знак +. Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

ZZ

 

y2 z2

 

dxdy = ZZ

 

y2 a2 x2

 

 

 

 

a

 

 

 

b

 

 

 

 

 

 

 

 

 

 

a2 :

 

 

 

dxdy = Z a dx Z0

y2 + x2 a2 dy = 3ab 2b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

цилиндра

 

 

 

 

,

 

 

 

 

RR

 

 

 

 

 

 

è

 

 

 

, (

 

 

).

 

 

 

 

 

 

 

 

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

 

 

 

x2 + z2

+ 2y dxdz, где - внешняя часть

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отсеченная плоскостями

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = 2y

 

y = 2 z = x x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение Проекцией поверхности на плоскость OXZ будет треугольник

 

 

 

 

 

D. Нормаль к внешней стороне цилиндра составляет тупой угол с осью ОУ,

 

 

 

 

 

поэтому в формуле надо взять знак .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ZZ

x2 + z2 + 2y

dxdz = ZZ

 

x2 + z2 + x2

 

dxdz =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ZZ

 

2x2 + z2

 

dxdz = Z0

2

 

 

x

2x2 + z2

 

 

dz =

 

3

:

 

 

 

 

 

 

 

 

 

 

 

 

dx Z0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

28

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рона части сферы

 

 

 

RR

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

y2

+ z2 + 3x2

 

dydz, где - внутренняя сто-

 

 

 

 

 

 

 

 

 

 

 

 

 

x = p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

2

y

2

z

2 вырезанной поверхностью

 

y

2

+ z

2

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

292

Решение Проекцией поверхности на плоскость OYZ будет круг. Чтобы найти его радиус, рассмотрим систему уравнений наших поверхностей

(

x2 + y2 + z2 = R2

:

y2 + z2 = x2

Исключая из этой системы , получим уравнение окружности y2 + z2 =

R2

2 , которая является проекцией линии пересечения поверхностей на плоскость OYZ, т.е. границей области D. Нормаль к внутренней стороне состав-

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

ZZ

y2

+ z2

+ 3x2

 

dydz = ZZ

y2

+ z2

+ 3 R2

y2 z2

 

dydz = ZZ

3R2

2 y2

+ z2

dydz =

 

 

 

ZZ

 

2 ZZ

D

 

 

 

 

D

 

 

 

 

= 3R2

dydz +

y2 + z2 dydz

 

 

 

 

 

 

 

:

DD

RR

Первый из полученных интегралов равен площади области D : dydz =

D

R2

2 . Для вычисления второго интеграла перейдем к полярным координатам: x = r cos '; y = r sin '. Тогда

 

 

 

 

 

 

 

 

 

 

R/p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

R4

 

 

ZZ

 

y2

+ z2

 

dydz = Z0

d' Z0

 

 

 

r3dr =

 

:

 

 

 

 

 

8

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ZZ

 

 

 

 

 

 

 

R4

 

3 R4

 

5 R4

y2 + z2 + 3x2

dydz =

 

 

 

 

=

 

 

 

:

4

2

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

RR x2dydz + y2dzdx + z2dxdy, где - внешняя

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сторона сферы (x a)2 + (y b)2 + (z c)2 = R2:

 

 

 

 

 

Решение Разобьем интеграл на сумму трех слагаемых

ZZ

x2dydz + y2dzdx + z2dxdy = ZZ

x2dydz + ZZ

y2dzdx + ZZ

z2dxdy

и вычислим каждое из них.

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

части 1 è 2 плоскостью z = c. Тогда каждая из этих частей проектируется на плоскость ОХУ взаимно однозначно и дает в проекции круг (x a)2 +

(y b)2 = R2.

q

Кроме того, поверхность 1 задается уравнением z = c+ R2 (x a)2 (y b)2

и нормаль к ее внешней стороне образует острый угол с осью OZ, à ïî- q

верхность 2 задается уравнением z = c R2 (x a)2 (y b)2 è íîð- маль к ее внешней стороне образует тупой угол с той же осью. Отсюда

293

RR

RR

RR

 

 

 

 

 

 

 

 

 

 

 

=

+

è

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

ZZ z2dxdy = ZZ c +

 

 

 

 

 

 

 

2 dxdy =

 

 

R2 (x a)2

(y b)2

 

1

 

D

q

 

 

 

 

 

 

;

=

c2 + R2 (x a)2 + (y b)2

+ 2c R2

(x a)2 (y b)2

dxdy

ZZ

 

 

 

 

 

 

q

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

ZZ z2dxdy = ZZ c

 

 

 

2 dxdy =

 

R2 (x a)2 (y b)2

 

2

 

D

q

 

 

 

 

 

 

:

ZZ

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

2c R2 (x a)2 (y b)2

 

c2 + R2 (x a)2 + (y b)2

dxdy

D

 

 

 

 

 

 

 

 

 

 

 

 

 

Складывая эти интегралы, получим

 

 

 

 

 

 

ZZ

 

 

ZZ q

 

 

 

 

 

 

 

 

z2dxdy = 4c R2 (x a)2 (y b)2dxdy:

D

Последний интеграл вычислим с помощью перехода к смещенным полярным координатам: x a = r cos '; y b = r sin ' :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ZZ qR2 (x a)2 (y b)2dxdy = Z0

 

d' Z0

 

R2

r2

rdr =

2

 

 

:

 

 

 

 

 

 

3

 

 

 

 

 

 

 

D

 

 

 

 

RR

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

z2dxdy =

8 cR3

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

8 aR3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

y2dzdx = 8 bR3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

RR

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x dydz =

 

 

 

 

 

 

 

 

 

Другие интегралы вычисляются аналогично и равны

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

RRОкончательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ZZ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2dydz + y2dzdx + z2dxdy =

8

 

(a + b + c) :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

- внешняя сторона

 

RR

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, вырезанная

 

 

 

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

(y z) dydz + (z x) dzdx + (x y) dxdy , ãäå

 

 

 

 

 

 

 

 

 

 

 

 

 

части сферы

 

2

+ y

2

+ z

2

= 2Rx; z 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поверхностью x

+ y

 

 

= 2ax; R > a > 0:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение Воспользуемся формулой перехода от интеграла второго ро-

 

 

 

да к интегралу первого рода (90). Для этого сначала найдем орт внешней

 

 

 

нормали к поверхности: ~n = (2x 2R; 2y; 2z) и ~n0 =

xRR

;

y

;

z

. Тогда

 

 

 

R

R

 

 

 

ZZ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ZZ

 

 

 

 

R

 

 

 

 

 

 

 

 

R

 

R

 

 

(y

 

z) dydz + (z

 

 

x) dzdx + (x

 

y) dxdy =

 

 

 

 

(y

 

 

z) (x

 

R)

+

 

(z x) y

+

(x y) z

dS =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= R ZZ

R (z y) dS = ZZ (z y) dS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

294

Для вычисления последнего интеграла спроектируем поверхность на плоскость ОХУ. Получим

ZZ

 

 

ZZD

 

 

jcos j

ZZD

z

ZZD

 

 

ZZ

z

(z

 

y) dS =

(z

 

y)

dxdy

= R

 

z y

dxdy = R

dxdy

 

R

 

y

dxdy:

 

 

 

 

 

 

 

 

RR

Первый из полученных интегралов равен площади области D: dxdy =

D

a2. Во втором интеграле подынтегральная функция является нечетной

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

Окончательно,

ZZ

(y z) dydz + (z x) dzdx + (x y) dxdy = Ra2

95Скалярные и векторные поля

1.Пусть в пространстве R3 задана некоторая область (которая, в частности может совпадать с пространством R3). И пусть также каждой точке области сопоставляется некоторое число, т.е. на области за-

дана некоторая функция u = f (x; y; z). Тогда говорят, что в области

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

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

Уравнение f (x; y; z) = C, где - произвольная постоянная, задает по-

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

Если область ограничена цилиндрической поверхностью и в лю-

бой точке прямой, проходящей через область параллельно образующей этой поверхности значения функции f (x; y; z) равны между со-

бой, то поле называют плоским. Тогда обычно считают, что образующая цилиндра параллельна оси OZ, и задают поле функцией вида u = f (x; y) :

Если поле плоское, то уравнение f (x; y) = C задает множество линий, которые называются линиями уровня.

2.Если в пространстве R3 задана область и каждой точке этой области сопоставляется вектор, т.е. задана векторная функция

~a = ~a (x; y; z) = (P (x; y; z) ; Q (x; y; z) ; R (x; y; z)) ;

то говорят, что в области задано векторное поле.

Физики изучают силовые векторные поля, поля скоростей и т.д.

295

Линия в пространстве, касательная к которой в каждой точке совпадает с направлением поля, называется векторной линией этого поля. Если вспомнить, что координаты вектора касательной к кривой может быть заданы с помощью дифференциалов (dx; dy; dz), то ясно, что

векторная линия задается системой дифференциальных уравнений

dx

=

 

dy

=

dz

:

P (x; y; z)

 

Q (x; y; z)

R (x; y; z)

Можно доказать, что, если вектор ~a - ненулевой, то через каждую точку области проходит одна и только одна векторная линия.

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

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

Как и для скалярного поля определяется плоское векторное поле.

Пример 269 Найти линии уровня плоского поля f (x; y) = x2 + y2 2x:

Решение Это поле определено в любой точке плоскости. Линии уровня заданы уравнениями x2 + y2 2x = C. Очевидно, что при фиксированном

- это окружность. Для того, чтобы найти центр окружности и ее радиус, выделим полный квадрат: (x 1)2 + y2 = C + 1.

Следовательно, линиями óðîâня будет семейство окружностей с центром

 

в точке

(1; 0)

и радиусом p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C + 1; C 1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 270 Найти область , где определено скалярное поле u = arcsin

p

 

z

:

x2+y2

Решение Поле определено в области, где

p

 

1 èëè jzj

 

 

:

 

 

x2

+ y2

 

x2+y2

 

Эта область находится вне конуса

2

2

 

 

2

 

z

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

z = x

 

+ y

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поверхности уровня определяются уравнениями

 

 

 

 

 

z

 

= C;

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

arcsin px2+y2

2

 

C 2 èëè

p

 

= sin C. Åñëè

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

 

x2+y2

 

занном промежутке, то sin C тоже произвольная постоянная, лежащая в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

промежутке [ 1; 1]. Поэтому можно обозначить sin C = C и тогда поверх-

 

ности уровня будут задаваться уравнениями

 

~

 

 

2

 

 

 

2

; 1

~

 

 

 

 

 

что равносильно уравнениям

 

;

 

z = Cpx

 

+ y

 

C 1,

 

 

 

 

 

 

z2 = C~2 x2 + y2

C~

1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, поверхностями уровня данного

поля будет семейство кону-

 

ñîâ

 

 

 

 

z2 = C~2 x2 + y2

;

C~

1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ ~ ~

Пример 271 Найти векторные линии поля ~a = xi yj 2zk :

296

Решение Векторные линии определяются системой дифференциаль-

ных уравнений

dxx = dydy = dz2z :

Решить эту систему можно, решив по отдельности два уравнения

dxx = dyy è 2dyy = dzz :

Интегрируя первое из уравнений, получим xy = C1, а из второго - z = C2y2. Таким образом, искомыми векторными линиями будут линии пересечения поверхностей xy = C1 è z = C2y2:

96Дифференциальные характеристики скалярного поля

Пусть функция, определяющая скалярное поле дифференцируема. Тогда у этого поля существуют две дифференциальные характеристики.

1.Градиент.

Градиентом скалярного поля u = f (x; y; z) называют вектор

~g =

@x ;

@y ;

@z

:

 

 

@f

@f

@f

 

Обозначают этот вектор grad U:

2.Производная по направлению.

Пусть ~

l ненулевой вектор, задающий некоторое направление в про-

странстве, и ~

l0 = (cos ; cos ; cos ) - орт этого вектора. Возьмем в области , где задано скалярное поле, точку M0 (x0; y0; z0). Тогда про-

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

@U

= lim

f (x; y; z) f (x0; y0; z0)

;

@l

t

t!0

 

ãäå x = x0 + t cos ; y = y0 + t cos ; z = z0 + t cos , т.е. точка M (x; y; z)

берется на прямой, проходящей через точку M0 (x0; y0; z0) параллель-

но заданному вектору ~ l:

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

97Свойства градиента и производной по направлению

1.Рассмотрим поверхность уровня скалярного поля f (x; y; z) = C. Как известно, вектор

@f@x ; @f@y ; @f@z

297

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

2.Выведем формулу для вычисления производной по направлению. Для

этого рассмотрим вспомогательную функцию ' (t) = f (x0 + t cos ; y0 + t cos ; z0 + t cos ). Тогда по определению производной по направлению можно написать

 

 

 

 

@U

= lim

' (t) ' (0)

 

= '0 (0) :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@l

t!0

t

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

t=0 =

 

'0 (0) = (f (x0 + t cos ; y0 + t cos ; z0 + t cos ))0

t=0 = @x dt

+ @y dt +

@z dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@f dx

@f dy

@f dz

 

@f (x0; y0

; z0)

@f (x0

; y0; z0)

 

 

 

; z0)

 

 

 

 

 

:

 

@f

(x0; y0

 

 

 

 

 

 

=

 

 

cos +

 

 

 

 

cos +

 

 

 

 

cos = grad U

 

~l0

 

 

@x

 

 

 

 

 

 

@z

 

 

 

 

 

 

 

 

 

@y

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

@U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@l

= grad U l0:

 

 

 

 

(94)

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@U

= jgrad Uj jcos 'j ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где ' - угол между градиентом скалярного поля и вектором, задаю-

щим направление ~ l:

Отсюда получается очевидное неравенство

@U

jgrad Uj ;

@l

причем равенство достигается, когда jcos 'j = 1, т.е. когда направле-

íèå ~

l совпадает с градиентом скалярного поля или противоположно

åìó.

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

4.Производная по направлению не зависит от выбора системы координат. Это следует из ее определения, если переписать его следующим образом:

@U

= lim

f (M) f (M0)

:

@l

t

t!0

 

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

Пример 272 Найти градиент скалярного поля u = x2+y2 в точке M0 (1; 3) :

298

Решение Градиент данного скалярного поля в произвольной точке равен (2x; 2y). Значит, в точке M0 градиент будет равен grad U (1; 3) = (2; 6). Следует отметить, что этот вектор ортогонален окружности x2 + y2 = 10,

являющейся линией уровня скалярного поля, проходящей через точку M0.

Пример 273 Найти градиент потенциала электростатического поля U =

re ; ~r = (x; y; z) ; r = j~rj :

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

Далее, вычислим частные производные

 

 

 

 

 

 

 

 

 

 

U :

 

 

Решение Вычислим длину вектора ~r : j~rj =

x2

+ y2 + z2:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

функции

 

 

 

 

 

 

@U

=

dU @r

=

e

x

 

 

=

ex

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@x

dr

@x

r2

p

 

 

 

 

r3

 

Аналогично,

 

 

 

x2 + y2 + z2

 

 

 

 

@U

=

ey @U

=

ez

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@y

r3

@z

r3

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

e~r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

grad U (x; y; z) =

 

(x; y; z) =

 

 

:

 

 

 

 

r3

r3

 

 

Пример 274 Дано скалярное поле U = x3 3x2y + 3xy2 + 1 и две точки

A (3; 1)

è

B (6; 5)

. Вычислить производную по направлению

~

AB данного

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l =

!

 

 

поля в точке , пользуясь определением этой производной.

 

 

 

 

 

Решение Найдем сначала орт направления:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

~

 

 

 

~

 

 

3 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

AB = (3; 4) ; l = 5; l

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

5; 5 :

 

 

 

 

 

 

 

 

 

 

 

 

l = !

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Затем найдем значения функции в точке

 

: U (3; 1)

= 10 и в точке

M (x0 + t cos ; y0 + t cos ) = M 3 +

3 t; 1 +

 

4 t

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для упрощения вычислений функцию

запишем в виде

 

 

 

 

3

U = (x

 

y) +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y3 + 1. Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

3

4

 

 

3

 

 

 

 

 

 

 

 

 

 

54

 

63

 

 

 

 

 

U (M) = 2

 

 

 

 

+ 1 +

 

t

+ 1 = 10 +

 

 

 

t2

+

 

 

t3

 

 

 

 

5

 

5

25

125

 

 

По определению производной по направлению

 

 

 

 

 

 

 

 

 

 

 

 

 

@U (A)

= lim

f (M) U (A)

= lim

 

54t

+

 

63t2

 

= 0:

 

 

 

 

 

@l

 

 

 

 

 

125

 

 

 

 

 

 

 

 

t!0

 

 

t

 

 

 

 

t!0

25

 

 

 

 

 

 

 

 

Пример 275 Дано скалярное поле U = xy2 + z3 xyz. Найти @U@l

â òî÷-

~

ке A (1; 1; 2), если вектор l образует с координатными осями углы равные соответственно 60o; 45o; 60o:

299

Решение Найдем сначала орт направления:

 

 

 

 

 

 

 

1

 

p

 

 

1

!:

~l0 = (cos ; cos ; cos ) =

;

2

;

2

2

 

2

Затем найдем градиент данной функции в точке

:

grad U (1; 1; 2) =

y2 yz; 2xy xz; 3z2 xy (1;1;2) = ( 1; 0; 11) :

 

 

Для вычисления производной по направлению применим формулу ():

@U (1; 1; 2)

~

 

1

 

11

 

 

 

= grad U (1; 1; 2) l0

=

 

+

 

= 5:

@l

2

2

98Дифференциальные характеристики векторного поля

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

1. Дивергенция.

~

~

~

Дивергенцией векторного поля

 

~a = P (x; y; z) i+Q (x; y; z) j+R (x; y; z) k

в некоторой точке M0 (x0; y0; z0) области называется выражение

div~a = @P@x + @Q@y + @R@z ;

где производные вычисляются в этой точке.

2.Ротор.

Ротором векторного поля в некоторой точке M0 (x0; y0; z0) области

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

rot~a =

@R

 

@Q @P

 

@R @Q

 

@P

;

 

 

;

 

 

;

 

 

 

@y

@z

@z

@x

@x

@y

где производные вычисляются в этой точке.

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

 

~i

~j

~k

 

 

 

 

 

 

@ @ @

rot~a =

 

 

 

 

 

:

@x

@y

@z

 

 

 

 

 

 

 

 

P

 

Q

 

R

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

rot~a =~i

 

@

@

 

~j

 

@

@

 

+ ~k

 

@

 

@

 

 

 

Q

 

R

 

P

 

R

P

 

Q

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@y

@z

 

 

 

 

 

 

 

 

 

 

@x

@y

 

 

 

 

@x

 

@z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

300