Высшая математика. В 7 ч. Ч. 7. Элементы теории поля
.pdf
Видно, что в данном направлении функция убывает (производная отрицательна).
Направление наибыстрейшего возрастания функции определяется градиентом gradu(M ) 4i 2 j 2k , поэтому косинус искомого угла равен:
|
|
|
|
e |
|
|
|
u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
cos |
|
grad u(М) |
|
|
|
|
r |
|
М |
|
|
|
|
|
4 |
|
|
|
2 |
, |
||
|
|
|
|
|
|
|||||||||||||||||
|
grad u(М) |
|
|
|
|
grad u(М) |
|
|
( 4) |
2 |
2 |
2 |
( 2) |
2 |
3 |
|||||||
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где учтено, что скалярное произведение градиента на орт направления равно
производной по этому направлению: |
u |
|
gradu(M ) e . |
|
|||
|
r |
|
М |
|
|
5. Найти скорость и направление наибыстрейшего возрастания функции u xsin z y cos x в начале координат.
Решение. Найдем градиент функции в начале координат:
u |
|
|
sin z ysin x |
|
M |
0, |
u |
|
cos x |
|
M |
1, |
u |
|
xcos z |
|
M |
0, |
gradu 1 j . |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|||||||||||||||
x |
|
M |
|
|
|
y |
|
M |
|
|
z |
|
M |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, в начале координат наибыстрейшее возрастание функции происходит со скоростью gradu 1 в положительном направлении оси y.
6. Найти единичный вектор, который перпендикулярен поверхности уровня скалярного поля u xy xz yz в точке M 1; 1; 1 .
Решение. Так как в каждой точке поля градиент перпендикулярен поверхности уровня, проходящей через эту точку, для решения задачи достаточно найти единичный вектор в направлении градиента поля в точке M 1; 1; 1 :
u |
|
|
y z |
|
M |
2, |
u |
|
|
|
x z |
|
M |
2, |
|
u |
|
x y |
|
M |
2, |
|
gradu 2i 2 j 2k, |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
x |
|
M |
|
|
|
|
|
|
y |
|
M |
|
|
|
|
|
z |
|
M |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
j 2k |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
grad u(M ) |
|
|
|
|
|
2i 2 j 2k |
|
|
2i |
2 |
|
1 |
|
|
1 |
|
|
1 |
|
|||||||||||||
|
|
|
n0 |
|
|
grad u(M ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
j |
|
|
k . |
||||||
|
|
|
|
|
|
22 22 22 |
|
|
|
2 |
3 |
|
|
3 |
3 |
3 |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
10
7. Вычислить координаты единичного вектора, направленного по нормали к поверхности z2 x2 xyz y5 5 в точке M 1; 1; 2 .
Решение. По нормали к поверхности направлен вектор градиента, который при неявном задании поверхности уравнением вида F (x, y, z) 0 вычисляется по формуле grad F (Fx , Fy , Fz ) .
Найдем градиент и его орт:
F' |
|
|
yz3 3x2 yz |
|
2, |
F' |
|
|
xz(z2 x2) 5y4 |
|
|
|
1, F |
' |
|
|
|
|
3z2xy x3 y |
|
11, |
|||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
x |
|
M |
|
|
|
M |
|
y |
|
M |
|
|
|
|
|
|
|
|
M |
|
z |
|
M |
|
|
|
M |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
grad F 2i 1 j |
11k , |
|
grad F |
|
|
22 12 112 |
|
|
|
|
9 14 3 14 , |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
e |
grad F |
|
|
2i 1 j 11k |
|
|
|
2 |
i |
|
1 |
j |
|
|
|
11 |
k . |
|
|
||||||||||
|
|
|
grad F |
|
|
|
3 14 |
3 |
14 |
3 |
14 |
|
|
|||||||||||||||||||
|
|
|
|
|
3 14 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
Таким образом, существуют два взаимно противоположных единичных вектора, направленных по нормали к заданной поверхности: n0 е.
8. Найти векторные линии векторного поля a(M ) 4zj 9 yk . Найти поле дивергенции и поле ротора.
Решение. Для заданного векторного поля P 0, Q 4z, R 9 y , и его векторные линии определяются системой дифференциальных уравнений
d0x d4yz dz9 y .
Из уравнения d4yz d9zy следует уравнение 4zdz 9 ydy 0, интегрируя ко-
торое получим векторные линии в виде эллипсов, центры которых лежат на оси х:
z2 y2 C1. 32 22
11
Из уравнения d0x d4yz следует dx 0 , что дает x C2 . Таким образом,
векторные линии данного поля – это эллипсы, лежащие в плоскостях, перпендикулярных оси х.
Найдем дивергенцию и ротор в произвольной точке поля:
|
i |
j |
k |
9 4 i 0 0 j 0 0 k 13i . |
div a 0, |
rot a x |
y |
z |
0 4z 9 y
Таким образом, дивергенция и ротор поля постоянны, причем дивергенция равна нулю.
2. ЗАДАЧИ ДЛЯ САМОСТОЯТЕЛЬНОГО РЕШЕНИЯ
Найти линию уровня плоского скалярного поля, проходящую через заданную точку:
1. z |
x2 y2 , M 2; 2 3 . |
2. |
u(x, y) |
|
x |
, M 1; 1 . |
|
x2 |
3y2 |
||||||
|
|
|
|
|
Найти поверхность уровня скалярного поля, проходящую через заданную точку:
3.u (x 2)2 ( y 5)2 (z 3)2 , M 2; 2; 1 .
4.u 2x2 y2 z2 , A 1; 0; 1 .
5. Найти градиент скалярного поля u |
x2 y2 z2 |
в точке M 1; 1; 1 . |
|||||
6. Показать, что в точках M1 1; 0; 1 |
и M2 1; 1; 0 |
|
градиенты скалярной |
||||
функции u x2 yz xy2z xyz2 ортогональны. |
|
|
|
|
|
||
|
|
|
|
1 |
|
|
|
7. Определить точки, в которых градиент поля u ln x |
равен |
i j . |
|||||
|
|||||||
|
|
|
|
y |
|
||
12
8. |
Рассчитать производную скалярного поля |
u |
xy |
в точке M1 1; 1; 1 |
|
z |
|||||
|
|
|
|
||
по направлению к точке M2 3; 2; 3 . |
|
|
|
||
9. |
Найти производную поля u arctg x2 y2 в точке M 1; 1 в направлении |
||||
четвертого координатного угла.
10. Найти угол между направлением скорости наибыстрейшего возраста-
ния функции u xy в точке M 0; 1 и отрицательным направлением оси х.
11.Показать, что скорости наибольшего возрастания функции u arctg(xy)
вточках M1 1; 1 и M2 1; 1 равны по модулю и противоположно на-
правлены.
12. |
Найти единичные векторы нормали к поверхности уровня скалярного |
|
поля u ln(xy xz yz) в точке M 1; 0; 1 . |
||
13. |
Определить направления векторов нормали к поверхности уровня |
|
плоского скалярного поля u 2x2 3y2 3x 9 y 3 в точке M 0; 1 . |
||
14. |
Найти |
единичный вектор нормали к поверхности ln(1 x2 ) yz 0 |
в точке |
M 1; |
0; 1 , составляющий острый угол с положительным направле- |
нием оси х.
15.Найти векторные линии векторного поля a(M ) xi 4zk .
16.Найти векторные линии поля gradu , если u x y2 .
17.Найти дивергенцию векторного поля a x2 yx i y2 yx j .
18.Найти ротор векторного поля a xyi y2 j zxk .
19.Найти grad(div a) , где a x3i y3 j z3k .
20.Доказать, что для любого векторного поля справедливо равенство div(rot a) 0 .
13
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ответы |
|
|
|
|
|
|
|
|
1) окружность |
|
x2 y2 |
16 . |
2) эллипс |
x 2 2 |
|
y2 |
|
1. |
3) сфера |
||||||||||||||||
|
|
4 |
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 3 |
|
|
|
|||
(x 2)2 ( y 5)2 (z 3)2 25 . 4) однополостный гиперболоид |
2x2 y2 |
z2 |
1. |
||||||||||||||||||||||||
|
|
|
|
|
|
1 |
|
|
|
1 |
|
|
|
1 |
|
M1 2; 1 , M2 0; 1 . |
|
|
|
|
|
|
|||||
5) |
gradu |
|
|
|
|
|
|
; |
|
|
|
|
; |
|
|
. |
7) |
8) 1. |
|
9) 0. |
10) |
. |
|||||
|
|
3 |
|
3 |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
1 |
|
|
|
|
2 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
||||
12) |
е |
|
|
; |
|
|
|
|
|
|
; |
|
|
|
|
. |
13) биссектрисы 1-го и 3-го координатных углов. |
||||||||||
6 |
|
|
|
|
6 |
|
|
|
6 |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
14) |
п |
|
1 |
; 2 |
; |
2 |
. |
15) z C x4 |
, y C . |
16) z C |
, x 1 ln y C |
2 |
. 17) x 3y . |
|
|
|
3 |
3 |
|
3 |
|
1 |
2 |
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
18)zj . 19) 6xi 6 yj 6zk .
3.ПОТОК, РАБОТА И ЦИРКУЛЯЦИЯ ВЕКТОРНОГО ПОЛЯ. ТЕОРЕМА ОСТРОГРАДСКОГО-ГАУССА. ТЕОРЕМЫ СТОКСА
ИГРИНА. ВИДЫ ВЕКТОРНЫХ ПОЛЕЙ
Потоком векторного поля a P x, y, z i Q x, y, z j R x, y, z k через выбранную сторону поверхности S называется интеграл по поверхности S от скалярного произведения вектора поля a(М) на единичный вектор нормали n0 (М) к поверхности, т. е. поверхностный интеграл 2-го рода:
|
|
|
|
an0dS . |
|
|
S |
|
При изменении стороны поверхности на противоположную, т. е. при за- |
||
мене n0 на n0 , |
поток изменяет знак. С учетом n0 |
cos i cos j cos k |
и cos dS dydz, |
cos dS dxdz, cos dS dxdy, |
поверхностный интеграл |
2-го рода может быть записан в скалярной форме и вычислен с помощью двойных интегралов по проекциям Dyz, Dxz, Dxy заданной поверхности на координатные плоскости yOz, xOz, xOy соответственно:
14
|
|
|
|
|
an0dS (P cos Q cos R cos )dS |
||
|
S |
S |
|
|
P(x( y, z), y, z)dydz |
|
Q(x, y(x, z), z)dxdz R(x, y, z(x, y))dxdy, |
Dyz |
|
Dxz |
Dxy |
где функции x( y, z), y(x, z), z(x, y) получены разрешением уравнения поверхности относительно соответствующей координаты. Если S задана в явном виде функцией z z(x, y) , то справедлива формула
|
|
a n0dS |
a n dxdy |
|
|
S |
Dxy |
|
|
zxP(x, y, z(x, y) zyQ x, y, z(x, y) R x, y, z(x, y) dxdy, |
|
|
Dxy |
|
|
где n zx , |
zy , 1 – вектор нормали к поверхности. Аналогичным образом |
||
поток может вычисляться по проекции на плоскости yOz (если поверхность явно задана функцией x x( y, z) и имеет нормаль n (1, xy , xz ) ) или xOz (при
задании функцией y y(x, z) и нормали n ( yx , 1, yz ) ).
Если S замкнута и ограничивает объем V, а функции координат P, Q, R непрерывны вместе со своими частными производными 1-го порядка, то по теореме Остроградского-Гаусса поток вектора a через внешнюю сторону S равен тройному интегралу от дивергенции этого поля по области V:
|
|
|
|
P |
|
Q |
|
R |
an0dS divadV |
x |
y |
dxdydz. |
|||||
S |
V |
|
V |
|
|
z |
||
Заметим, что вычислить тройной интеграл часто оказывается значительно
проще. |
|
|
|
|
|
|
|
|
|
|
|
|
Напомним, |
что |
|
– скалярная величина. Если угол |
(а, n0 ) |
|
, то |
||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
0; если |
|
|
, |
то |
0; если |
|
|
, то |
0 . В гидродинамике поток |
|||
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
вектора скорости несжимаемой жидкости через заданную поверхность дает
15
объем жидкости, протекающей через нее в единицу времени, 0 означает, что из области V вытекает больше жидкости, чем втекает, т. е. внутри V имеются источники; 0 означает, что из V вытекает меньше жидкости, чем втекает, т. е. внутри V имеются стоки; при 0 в область V втекает столько же жидкости, сколько вытекает. В общем случае поток векторного поля пропорционален количеству векторных линий, пронизывающих поверхность в направлении вектора нормали в единицу времени.
Работой векторного поля a(М) по пути L называется криволинейный интеграл
A a 0dl Pdx Qdy Rdz,
L L
где L – кусочно-гладкая кривая;
0 (М) – единичный вектор в направлении касательной к ней.
С физической точки зрения данный интеграл равен работе переменной силы F Pi Qj Rk , действующей на материальную точку (физическое тело), при ее перемещении вдоль кривой L. Если кривая L замкнута, данный интеграл называется циркуляцией поля a по контуру L. Направление обхода контура либо указывается заранее, либо по умолчанию берется положительным, т. е. против часовой стрелки. Отметим, что вдоль замкнутых векторных линий циркуляция всегда отлична от нуля, так как на векторной линии скалярное произведение a 0 сохраняет знак: положительный, если направления вектора поля и обхода векторной линии совпадают, и отрицательный, если эти направления противоположны.
Согласно теореме Стокса, циркуляция вектора a по замкнутому контуру L равна потоку вектора rot a через незамкнутую поверхность S, ограниченную этим контуром (или натянутую на этот контур):
|
|
|
a 0dl rot a |
n0dS, |
|
L |
S |
|
16
где 0, n0 – единичные векторы, направленные соответственно по касательной
кконтуру и по нормали к поверхности.
Вскалярной форме:
|
|
R |
|
Q |
P |
|
R |
|
Q |
|
P |
|
P dx Q dy Rdz |
y |
dydz |
z |
dzdx |
x |
dxdy. |
||||||
L |
S |
|
z |
|
|
x |
|
|
y |
|||
В приведенных выше формулах при обходе контура выбранная сторона поверхности должна оставаться слева или, что то же самое, вектор нормали к поверхности должен быть направлен так, чтобы видимый из его конца обход контура совершался против часовой стрелки. Если контур плоский, то в качестве поверхности часто выбирают ограниченную им область на плоскости.
Как видим, интегральные и дифференциальные характеристики векторного поля взаимосвязаны: циркуляция по замкнутому контуру выражается через ротор этого же поля, а поток через замкнутую поверхность – через дивергенцию. Не вдаваясь в подробности, отметим, что строгие определения ротора и дивергенции векторного поля даются через его циркуляцию и поток соответственно.
Очевидно, что формула Грина
|
|
Q |
|
P |
|
P dx Q dy |
dxdy |
||||
L |
D |
x |
|
y |
|
является частным случаем формулы Стокса, когда контур и ограниченная им область лежат в координатной плоскости xOy. Формула Грина позволяет вычислить площадь S плоской области D, ограниченной контуром L:
SD dxdy xdy |
y dx |
1 |
xdy ydx. |
||
D |
L |
L |
|
2 |
L |
|
|
||||
Рассмотрим некоторые специальные виды векторных полей. Поле a назы-
вается потенциальным (или градиентным, или безвихревым), если в каждой
17
точке M x, y, z оно является |
градиентом некоторой скалярной функции |
||||||||
u x, y, z : |
|
|
|
|
|
|
|
|
|
|
u(М) |
|
u(М) |
|
u(М) |
|
|||
а(М) gradu(М) |
x |
i |
|
y |
j |
|
z |
k |
, |
которая называется потенциалом поля. Чтобы векторное поле a было потенциальным, необходимо и достаточно, чтобы во всех точках поля его ротор был равен нулю:
|
|
|
R Q |
|
P R |
|
Q P |
|
rot a |
0 |
|
y z |
0, |
z x |
0, |
x y |
0. |
Потенциальное векторное поле a обладает следующими свойствами:
1) работа поля a не зависит от формы пути или, что равносильно, циркуляция поля по любому замкнутому контуру равна нулю:
a 0dl P dx Q dy Rdz 0;
LL
2)работа поля при перемещении из точки A в точку B не зависит от пути
иопределяется приращением потенциала поля:
A a 0dl u(B) u( A);
AB
3) потенциал поля вычисляется с помощью криволинейного интеграла:
M |
|
(x, y,z) |
|
u(M ) |
a 0dl |
|
P dx Q dy Rdz С, |
M0 |
|
(x0 , y0 ,z0 ) |
|
где точки M0 x0, y0, z0 и M x, y, z |
принадлежат полю; |
||
С – произвольная постоянная. |
|
||
18
Аналогичная ситуация имеет место, если векторное поле является плоским. В частности, если в поле a(М) P(x, y)i Q(x, y) j выполняется условие
|
|
|
|
|
|
|
Q P |
, |
|
|
|
|
||
|
|
|
|
|
|
|
x |
|
y |
|
|
|
|
|
|
|
Q |
|
P |
|
, то выражение Pdx Qdy является пол- |
||||||||
означающее, что rot a |
|
x |
|
|
k 0 |
|||||||||
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
ным дифференциалом потенциала u(x, y) , вычисляемого по формулам |
||||||||||||||
x |
|
|
y |
|
|
|
|
|
x |
|
|
|
y |
|
u(M ) P(x, y0 )dx Q(x, y)dy С P(x, y)dx Q(x0, y)dy С, |
||||||||||||||
x0 |
|
y0 |
|
|
|
|
|
x0 |
|
|
y0 |
|||
где точки M0 x0, y0 |
и M x, |
y принадлежат полю; |
|
|
||||||||||
С – произвольная постоянная. |
|
|
|
|
|
|
|
|||||||
Векторное поле |
a |
называется соленоидальным (или трубчатым), если |
||||||||||||
в каждой точке поля его дивергенция равна нулю: |
|
|
|
|||||||||||
|
|
|
|
|
P(М) |
|
Q(М) |
|
R(М) |
|
||||
div a |
0 |
|
x |
|
|
y |
|
z |
|
0. |
||||
Свойства соленоидального векторного поля:
1) поток вектора поля через любую замкнутую поверхность равен нулю:
П an0dS 0;
S
2)поле не содержит ни источников, ни стоков, векторные линии поля непрерывны и замкнуты;
3)поток вектора поля через любое сечение векторной трубки есть величина постоянная (отсюда название – «трубчатое»);
19
