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

Конспект лекций Высшая математика (Басканова)

.pdf
Скачиваний:
151
Добавлен:
17.01.2018
Размер:
5.08 Mб
Скачать

802

J a3 1 1 2a3 . 1 3 3 3

На линии L2: x2 z2 a2 , лежащей в плоскости xoy, z = 0, dz = 0, тогда

J2

(x2 y2 )dx (x2 y2 )dy .

 

L : x2

y2 a2

 

2

 

 

Так же как и в первом интеграле, уравнение окружности x2 z2

a2

возьмем в параметрической форме: х=acost; у=asint; dx=–asintdt; =acostdt. Перейдем к линейному интегралу

J2

a2 ( a sin tdt) a2 (cos2

t sin 2 t)a cos tdt ,

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J2

a3

sin tdt

(1 2sin2 t)d(sin t) ,

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

a

3

2 sin3 t

a

3

( 1 1) 0

2a

3

 

 

 

 

J2 a

 

cos t

0

sin t

 

3

 

 

 

.

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

Циркуляция поля по контуру L равна

 

 

 

 

 

 

 

 

 

 

 

 

 

ЦL

 

J1 J2

 

2a3

 

2a3

4a3

.

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

3

 

 

 

 

Теперь найдём циркуляцию данного поля по контуру L по теореме

Стокса

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

,

 

 

 

 

 

 

 

 

 

 

 

 

F

d ) (rotF

 

0 )d

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

L

 

q

 

 

 

 

 

 

 

 

 

 

 

Для этого нужно определить координаты ротора, т.е.

 

 

 

 

 

 

(2z

2

 

 

 

(x

2

y

2

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

)

 

 

ЦL F

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

z

 

cos(n, x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x2 y2 )

 

(2z2 )

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

x

 

cos(n, y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x

2

y

2

)

(x

2

y

2

)

 

 

 

 

 

d .

 

 

 

 

 

cos(n, z)

 

 

 

 

x

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

ЦL F 0 0 cos n, x 0 0 cos n, y 2x 2y cos n, x d ,

q

т.е. ротор данного поля во всех точках параллелен оси oz , поле является плоско-параллельным.

Поверхность q1 натянутая на контур L, однозначно проектируется на плоскость xoy, её проекцией является половина круга x2 z2 a2 . Циркуляция поля F равна двойному интегралу

803

ЦL F 2x 2 y dxdy .

qxoy

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

 

 

 

 

 

 

 

 

 

2 (x y)dxdy

 

x r cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЦL

 

 

y r sin

 

 

.

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

qxoy

 

 

 

 

 

 

 

dxdy r dr d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

ЦL

 

 

2 (r cos r sin )rdrd 2 (cos sin

)d r 2 dr ,

F

 

 

 

 

 

 

qxoy

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

a

2( 1 1) a

3

 

 

4a

3

 

 

 

 

 

 

 

 

 

 

 

 

 

ЦL

 

 

2 (sin cos )

0

 

 

 

 

 

 

.

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

0

3

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

804

ЛЕКЦИЯ 6.6. ВЕКТОРНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ ОПЕРАЦИИ ПЕРВОГО И ВТОРОГО ПОРЯДКА. СПЕЦИАЛЬНЫЕ ВИДЫ ВЕКТОРНЫХ ПОЛЕЙ И ИХ СВОЙСТВА

6.6.1. Векторные дифференциальные операции первого порядка

Напомним, что со скалярным полем u связывается векторное поле его градиента

gradu ux i uy j uz k ,

где x, y, z – декартовы координаты. Впрочем, градиент может быть определен и независимо от выбора системы координат, на основе понятия производной по направлению; поэтому градиент, как и дальнейшие рассматриваемые нами «производные поля », связан с заданным полем инвариантно. С векторным полем F связываются скалярное поле его дивергенции

 

 

 

div F

P

Q

R

 

 

 

и векторное поле его ротора

 

 

x

 

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

P

 

Q

P

 

 

R

 

Q

 

rot F

 

 

i

 

 

j

 

k ;

 

y

 

 

 

z

 

 

x

 

 

x

 

 

 

 

z

 

 

 

 

 

y

дивергенция также допускает инвариантное определение на основе по-

нятия потока, а ротор – на основе понятия циркуляции.

 

 

Гамильтон заметил, что операции grad u, div F, rot F можно более про-

сто записать, если ввести символ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

k

,

 

 

 

 

 

x

у

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

называемый набла (от греческого – арфа, форма которой напо-

минает значок ). Сам по себе этот символ представляет знак действия над полем, т.е. оператор. Этот оператор Гамильтона векторнодифференциальный и при своем действии обладает как свойствами вектора, так и свойствами оператора дифференцирования.

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

Рассмотрим эти операции. Пусть даны поля:

1.скалярное: u u (x, y, z),

2.векторное: F P (x, y, z),Q (x, y, z), R (x, y, z) .

1.Умножение вектора на скалярную функцию дает вектор, равный градиенту скалярного поля:

 

 

805

 

 

 

 

 

u (x, y, z)

u (x, y, z)

u (x, y, z)

,

u (x, y, z) i

x

j

y

k

z

 

 

 

 

 

grad u(P) u (x, y, z).

2. Применить символический вектор векторному полю можно двумя

способами: скалярно и векторно умножить на вектор F P,Q, R .

Скалярное произведение двух векторов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

k

 

z и F i P j Q k R

x

у

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

дает дивергенцию векторного поля

 

 

 

 

 

 

 

 

 

 

 

P x, y, z Q x, y, z

 

R x, y, z ,

 

F

 

 

 

 

 

 

x

 

 

y

 

 

 

z

div F P F .

3.Векторное произведение этих векторов дает ротор векторного поля

ij k

 

 

 

 

 

 

 

 

 

 

,

 

F

x

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

Q

 

R

 

rot(P) F .

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

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

избегать записи вида uv , которую более естественно понимать какu v (gradu)v=vgrad u, но иногда ее понимают как uv grad uv , а это,

конечно, не одно и то же.

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F iP

jQ kR i

 

 

j

 

k

 

 

P

 

Q

 

R

 

 

 

 

 

 

 

 

 

x

 

 

y

 

 

 

 

x

 

y

 

z

 

 

 

 

 

 

z

 

 

 

 

это скалярный дифференциальный оператор, который может действовать на скалярное или векторное поле. На основании формулы для производ-

ной по направлению этот оператор можно записать также в виде F , где

lF

lF – направление вектора F. В частности, для скорости изменения скалярного или векторного поля вдоль траектории получим выражение

 

 

 

806

 

 

 

 

 

d

v

 

.

 

 

 

dt

 

 

 

 

 

 

t

 

 

При действиях с оператором

 

надо пользоваться правилами вектор-

 

ной алгебры и правилами дифференцирования. Например,

 

rot(A B) (A B) A B rotA+rotB,

 

grad ( u) ( u) u gradu,

const ,

( )

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

 

 

 

( )

uv

u v v u

первое слагаемое получается, если в процессе дифференцирования считать v постоянным, а второе – если в этом процессе считать u постоянным. Поэтому дифференцирование ( ) можно выполнить так:

uv ucv uvc uc v u vc uv u v ,

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

grad (uv) (uv) uc v uvc u v v u ugradv+vgradu.

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

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

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

дексом с (const).

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

Пример 1. Используя символический метод, вычислить div a b . Решение. Воспользуемся свойствами векторного произведения:

div a b a b a bc ac b b a a b

b rot a a rotb .

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

807

6.6.2. Векторные дифференциальные операции второго порядка

Двукратное применение оператора (набла) к скалярным и векторным полям получило название векторных дифференциальных операций второго порядка.

Рассмотрим эти операции.

Нахождение градиента. Эту операцию можно применять только к ска-

лярным полям. После однократного применения оператора (набла), скалярной величиной является только дивергенция. Найдем градиент:

grad(div(F)) F .

Таким образом, двукратное применение оператора (набла) к векторному полю порождает новое векторное поле.

Следует отметить, что операция нахождения градиента от дивергенции на практике встречается крайне редко.

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

div(gradu) u u .

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

u x, y, z . Оператор 2 , равный скалярному произведению символи-

ческого вектора (набла) самого на себя, называют оператором Лапласа (лапласианом), его обозначают:

 

2

2

 

2

 

2

.

 

x2

y2

 

 

 

 

 

z2

Скалярный оператор Лапласа играет большую роль в математической физике. Найдем дивергенцию ротора:

div(rot(F))=( rot(F))

 

 

 

 

 

 

0 .

 

 

F

Это простое свойство имеет важное следствие. Именно, для любого поля F можно наряду с векторными линиями рассматривать вихревые линии, т.е. векторные линии поля rot F. Однако дивергенция любого векторного поля равна плотности источников векторных линий этого поля. Поэтому формула говорит, что вихревые линии не могут иметь ни источников, ни стоков, т.е. они не могут ни начинаться, ни кончаться.

Вычисление ротора. Ротор так же, как и дивергенция, характеризует векторное поле. Следовательно, его можно найти от градиента и ротора:

rot(grad u)=[ grad u] u 0.

Векторное произведение двух коллинеарных векторов равно нулю. Следовательно, поле градиентов не имеет вихрей:

rot(rot(F)) F .

808

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

rot(rot(F)) F F ; rot(rot(F))=grad(divF) F .

Таким образом, двукратное применение оператора (набла) к скалярным и векторным полям дает пять дифференциальных операций, две из них div(rot(F)) и rot(grad u) тождественно равны нулю. Остальные порождают новые векторные и скалярные поля, среди которых особое значение имеет скалярное поле лапласиана div(grad u).

Дифференциальные операции второго порядка удобно свести в табли-

цу:

 

 

 

 

 

 

 

 

Скалярное поле

 

 

 

 

 

Векторное поле

 

 

 

 

 

 

 

 

grad

 

 

 

div

 

rot

 

 

grad

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

grad(div a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div

 

div(gradu)= u

 

 

 

 

 

 

 

 

 

div(rot a)=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot

 

rot(grad u)=0

 

 

 

 

 

 

 

 

 

rot(rot a)=

 

 

 

 

 

 

 

 

 

 

 

 

grad(div a) – a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2.

Законы электромагнетизма

 

 

описываются

уравнениями Мак-

свелла в дифференциальной форме:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, H

 

 

 

,

 

 

 

 

 

 

 

0 ,

 

 

 

 

 

 

E

0.

 

 

H

 

E

 

E

 

H

 

с t

 

с t

 

 

 

 

 

 

 

Иначе:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

rotH,

(6.6.1)

 

 

 

 

 

 

 

 

 

 

 

с t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div

E

0 ,

 

 

 

 

 

(6.6.2)

 

 

 

 

 

 

 

 

 

 

 

H rotE,

(6.6.3)

 

 

 

 

 

 

 

 

 

 

 

 

с t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div

H

0

 

 

 

 

 

(6.6.4).

В данном случае нет зарядов и токов, а E , H – векторы напряженности электрического и магнитного полей; , – электрическая и магнитная

проницаемость; с – скорость света.

Если продифференцировать (6.6.1) по t и подставить Ht из (6.6.3), то получим

 

 

E

 

H

 

или

 

E E .

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

с t

2

 

t

 

с

2

t

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Преобразуем правую часть по формуле:

809

E , E E .

Итак, для векторного поля E имеем уравнение

2 E

 

c2

 

 

.

E

t 2

 

 

 

 

 

Мы получили одно из основных уравнений математической физики, называемое волновым уравнением.

Вывод: мы ввели дифференциальные операции первого и второго порядка и рассмотрели принципы работы с этими операторами.

6.6.3. Простейшие векторные поля

6.6.3.1. Потенциальное поле

К простейшим векторным полям относят потенциальное, соленоидальное и гармоническое. Рассмотрим свойства этих полей.

Определение 1. Векторное поле называют потенциальным, если оно образовано градиентами скалярной функции:

F grad u(x, y, z) .

Почти все силовые поля являются потенциальными (поле сил тяжести, электростатическое поле, магнитное и т.д.). Функцию u x, y, z , градиент ко-

торой образует векторное поле, называют его потенциалом. Координаты век-

тора потенциального поля

 

 

 

 

 

 

 

F

 

P x, y, z

 

Q x, y, z

 

R x, y, z

 

 

i

j

k

 

 

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

 

 

P x, y, z u ;

Q x, y, z u ;

 

R x, y, z

 

u .

 

 

x

 

 

y

 

 

 

z

Поэтому потенциальное поле можно задать одной скалярной функцией u x, y, z . Потенциальное поле обладает двумя очень важными свойствами.

Во-первых, оно не имеет вихрей, поскольку

rot(grad u) u 0,

т.е. ротор во всех точках потенциального поля равен нулю. Во-вторых, в потенциальном поле криволинейный интеграл по любому

замкнутому контуру равен нулю. Это следует из теоремы Стокса:

 

F

 

 

(rotF n0 )d 0 .

de

L

 

 

q

 

Из второго свойства вытекает важное следствие: в безвихревом потенциальном поле величина криволинейного интеграла по незамкнутой линии L не зависит от ее формы, а зависит только от координат начальной и конечной точек.

Пример 1. Проверить, будет ли потенциальным поле

 

 

 

 

 

 

 

 

 

 

 

810

 

 

 

 

 

 

 

 

 

 

 

 

10 xy3 12 x3 6

 

15 x2 y2 5 y

 

.

 

 

 

 

 

F

i

j

 

 

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

 

 

 

 

Q

 

P

 

 

 

15x2 y2 5y

 

10xy3

12x3

6

;

 

 

 

rot(F) k

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

x

 

 

 

y

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot(F)

 

30xy2 30xy2 0.

 

 

 

 

 

 

 

 

 

k

 

 

 

Данное поле является потенциальным, так как во всех точках ротор равен нулю.

Пример 2. Доказать, что поле

a 3x2 y2 z y2 z3 i 2x3 yz 2xyz3 j x3 y2 3xy2 z2 k

потенциально. Найти потенциал поля.

Решение. В потенциальном поле нет вихрей. Найдем ротор поля в произвольной точке:

 

 

 

 

 

 

i

 

 

 

 

 

j

 

 

 

 

 

k

 

 

 

 

 

 

rot a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3x2 y2 z y2 z3

2x3 yz 2xyz3 x3 y2

3xy2 z2

 

 

 

2x3 y 6xyz2

2x3 y 6xyz2

 

3x2 y2

3y2 z2

3x2 y2

3y2 z2

i

j

k 6x2 yz 2yz3 6x2 yz 2 yz3 i 0 j 0 k 0, rot a 0 .

А значит поле вектора a потенциально.

Для нахождения функции u x, y, z , которая является потенциалом

данного поля, возьмем криволинейный интеграл в этом поле по произвольной линии L от точки А до точки В:

B

 

 

 

 

B

 

 

 

a

 

dl

 

P xyz dx Q xyz dy R xyz dz .

A,поL

 

 

 

 

A

 

Так как поле потенциальное, его координаты равны частным производным от искомой функции u x, y, z :

P xyz

u

;

Q xyz

u

;

R xyz

u .

 

x

 

 

y

 

 

z

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

 

adl u dx u dy

u dz d u xyz .

B

 

 

 

 

B

 

 

B

A

 

 

 

 

A x

y

z

A

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

u B u A B adl .

A

811

Если верхнюю границу взять переменной B x, y, z , а нижнюю зафиксировать A x0 , y0 , z0 , то искомая функция

u xyz

B xyz

adl u x0 y0 z0

 

A( x0 y0z0 )

или

u xyz

B xyz

.

P xyz dx Q xyz dy R xyz dz u x0 y0 z0

A( x0 y0z0 )

От выбора точки A x0 y0 z0 зависит значение постоянной C u x0 y0 z0 . Возьмем за точку A начало координат: A 0,0,0 . В потенциальном поле кри-

волинейный интеграл не зависит от пути интегрирования, поэтому его можно вычислить как сумму трех интегралов по прямым L1 , L2 , L3 , параллельным

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

На прямой L1 : y 0; dy 0; z 0; dz 0.

На прямой L2 : x=const, dx 0; z 0; dz 0 .

 

 

 

 

 

 

На прямой L3 : x

const, dx 0; y=const, dy 0 .

 

 

 

 

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

 

равен сумме трех интегра-

a

лов:

B xyz

 

 

 

 

 

dy x3 y2 3xy2 z2 dz

u xyz

 

 

 

 

 

3x2 y2 z y2 z3 dx 2x3 yz 2xyz3

 

A 0,0,0

 

 

 

 

 

 

 

 

 

 

 

x

y

z

 

z

3

 

 

z

 

 

 

0 dx

0 dy x3 y2 3xy2 z2 dz x3 y2 z

3xy2

 

 

 

.

 

 

0

0

0

 

 

 

 

3

 

 

0

 

 

 

 

 

 

 

 

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

 

 

u x, y, z x3 y2 z xy2 z3 .

 

 

 

 

6.6.3.2. Соленоидальное поле

 

 

 

 

 

 

Определение 2.

Векторное поле

 

P,Q, R называют соленоидаль-

F

ным, если во всех его точках дивергенция равна нулю: div F p 0.

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

div(rotF) F 0 ,

а также поле линейных скоростей точек вращающегося тела. Соленоидальное поле обладает следующим свойством: поток вектора поля F P,Q, R

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