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

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

.pdf
Скачиваний:
11
Добавлен:
12.11.2023
Размер:
15.42 Mб
Скачать

Задача 1 1 , 10. Доказать,

что:

 

 

1)

§га<1 (<р +

Ф) =

§га<1 <р + $га<1 ф;

(11,27)

2 )

§га<1 («р • <|») =

? бгайф+.ф§гад<р)

(11,28)

3)

§га<1 Р (<р) =

Р' (<р) §га<1 <р,

(11,29)

где ® = <р (х, У. *)!

ф = ф (х,

«/,

г).

 

Р е ш е н и е . 1. Формула ( 11,27) следует прямо из формулы (11,4). Действительно,

ёга<1 (<р + Ф) = ^ (<Р + Ф) ‘ + 1; (<Р + ф) / + ^ (? + ф) к

поэтому

»*(*+«- (Й+йМй+йМй+й)*-

=(8?+Й/+Ж^+й?+Й*)-

=» §гас! <р + §гас1 ф

и формула (11,27) доказана. Несмотря на свою простоту она очень важна, так как с ее помощью решается вопрос о построении сум* мы векторных полей.

Частный случай. Если

ф =

С,

то,

поскольку в этом случае

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

 

 

 

 

дл

= а

=

а? -

п

дх

ду

 

Ъг

*

имеем

 

 

 

 

дгай (<р + С) = §га<1 ф.

2. Докажем теперь формулу (11,28)

§га<1 (<рф) = ^ (<рф) I + щ (<рф) / + ^ (<рф) \

Легко получаем:

ГуШ +

Умножая обе части каждого из этих равенств соответственно на (, и к и почленно складывая, имеем

0га<1 (<рф) = ф §га<1 <р + <р §га<1 ф.

Частные случаи: а) если 9 а ф, то

ёгаб <р* а 2<р§га<19}

б) если ф = С, то

§габ С9 = С§га<19.

3.

Формулу (11,29) докажите самостоятельно.

 

Задача

 

1 1 , 1 1 .

Доказать,

что,

если

г — радиус-вектор точки

А (*,

у, г), а <р =

<р(*, у, г) — функция трех независимых перемен-

-ных,

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а 3 ?

- §габ <р

 

 

 

(Нг • §га<1 9 — скалярное произведение векторов

3?

и §гай<р).

Р е ш е н и е .

Если функция

? =

<р(х, у,

г),

то

ее полный диф­

ференциал

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(А )

.Радиус-вектор

точки А(х, у, г)

 

 

 

 

 

 

 

 

 

 

~г = XI + у] + гк,

 

 

 

 

а его дифференциал

 

 

 

 

 

 

 

 

 

 

 

Зг == йх I + йу • / + & • Ъ.

 

 

Так как

вектор

 

 

 

 

 

 

 

 

 

 

 

 

б га ё ? = | Г + ! ? + ! * .

 

 

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

 

 

 

 

 

 

 

 

 

 

йг . 8га<1 9 =

 

+ р / у . +

р г 3 г .

 

-Сравнивая это выражение с выражением

(А) для

полного диффе­

ренциала функции

<р, получаем

требуемое равенство:

 

 

 

 

 

4<р =

йг ■§га<1 9.

 

 

 

 

Отсюда

можно сделать заключение, что, если

 

 

 

 

 

 

= йг>а,

 

 

 

(11,30)

то вектор

а

необходимо является

градиентом

некоторой функции

? = ?(*. У,

г):

 

 

 

 

 

 

 

 

а = ё,тад<?.

(11.31)

Задача 11,12 (для самостоятельного решения). Найти

0 , в ” '

+ „■ +

 

^ 1

 

~ ^

+ <«* + »*> *»•

Задача

11,13

(для самостоятельного

решения). Найти

 

 

 

 

йга<1 (ж* +

у* +

г2).

 

 

О т в е т .

2г.

 

 

 

 

 

 

 

 

Задача

11,14 (для

самостоятельного решения). Найти градиент

скалярного

произведения г» а, где Ъ— постоянный вектор, а г —

радиус-вектор точки А (х,_у,

г)

= х7+ у/ + &)■

 

Ответ,

бгас!(г а) = а.

 

 

 

 

 

дгаб (7 • а)

Полученный

результат указывает

на то,

что вектор

во всех точках

поля

сохраняет постоянное направление,

совпадаю­

щее с направлением

вектора а.

 

 

 

 

 

Кроме общего способа, основанного на

применении формулы

(11,4) к функции <р =

г а *

ах х +

ау у + аг • г, можно также

воспользоваться.результатом

задачи

(11,11). У нас функция <р ~

*- г а (о— постоянный вектор)

 

 

 

 

 

4<р= й (г . а) = а 4г.

Следовательно, на основании формулы (11,31)

а =* §га<1 =» §габ (7 ♦ о).

Читая это равенство справа налево, получаем требуемое

§га<1 (г .а ) — о.

 

Задача 11,15. Определить §га(1 <р(ы,

а), где и=*и(х, у,

г), о —

= ”(х. У, г).

 

 

 

 

опреде­

Р е ш е н и е . Чтобы воспользоваться формулой (11,4),

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

 

^

и

по формуле дифференци­

рования сложных функций

 

 

 

 

 

д у _д у д и . д у д о в

 

д х ~ д и д х ' до д х 9

 

д у __д у ‘ д и . д у д о #

 

д у ~ З а д у " * до д у 9

 

д у

д у

ди

. д у

д о

 

д г

д и

д г

до

д г *

 

Умножая обе части каждого из этих равенств соответственно на

I, / и к и почленно складывая их, получим

§гас! <?(и, о) = |? §габи + ^ §габV.

Задача 11,16. 1. Найти производную функции и = и(х, у, г) в направлении градиента функции V= V (х, у, г). 2. При каком усло­ вии эта производная равна нулю?

Р е ш е н и е . На основании формулы (11,4)

Проекция вектора §га(1о на оси прямоугольной системы координат

(8габп)ж = рх ; (%тайЬ)у = щ : (егас1ц)г = ^ .

Направляющие косинусы вектора а находят, как известно, по фор­ мулам

СОЗа ^ ; С058 = ^ ; созу = ^ .

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

 

 

 

 

(8га<1 о)х

 

до

 

соз (§га<1 V,

х) =

 

Т х

 

 

 

 

 

|ега<](*|

 

 

 

 

 

 

 

(вга<1о)у

 

до

 

соз (§га<1 V, у)

_

Т у

 

 

 

 

 

|вгай о |

 

| вга<1 о |

 

 

 

 

 

 

до

 

соз (дгаб о, г)

(«гай р)г

 

ъ

*

 

 

 

 

|ега<) о|

 

|йгас1«|

Производную функции

и = и(х, у, г)

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

7 = §га<1V, характеризуемому только что определенными косинусами,

найдем по формуле

(11,3)

 

 

 

 

 

 

 

д о

 

д о

 

, д а

д о

д и

 

д х

д и

д у

 

д г

Ж ~ ~ д х '

|егай о|

Гу' |вгай о

|

д г

| цгад о|

или

 

д и д о ,

д и д о

д и д о

 

 

 

 

 

д и

д х д х

д у д у ~ '

д г ^ г

 

 

В Т

 

|вгай о |

 

 

 

а так как

~ = Ыгайи)х;

=

(§га<1 и)„;

^

=

(дгас1 ы)^

3* = (бгас*°)*;

др =

(§га<1 у)у;

^

=

(§гас1 о)*,

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

ди_

&га<1и «%гас!о

 

ЗГ —

| {-гас! и |

*

Теперь переходим ко второму вопросу задачи.

Найденная производная ^ будет равна нулю, если числитель

последней дроби равен нулю, т. е. если

@габ и • §габ о = 0.

Из этого следует, что производная функции и = и (х, у, г) по направлению градиента функции о = о (*, у, г) равна нулю, если векторы §га(1 и и §габ V перпендикулярны.

С о д е р ж а н и е . Векторное поле. Потенциальные векторы.' Потенциал »бкторного поля. Циркуляция ректора. Линейный интеграл. Вихрь вектора.

ОСНОВНЫЕ СВЕДЕНИЯ ИЗ ТЕОРИИ

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

В случае векторного поля каждой точке Р пространства, в кото­ ром происходит образовавшее его физическое явление, ставится

в соответствие определенный вектор а{Р), являющийся функцией координат точки, т. е. а = а (х, у%г).

Вектор а(х, у, г) можно представить в виде

а(х, у\ г) = ах(х, у, г)1 + ау(х, у, г)] + аг(х, у, г ) \

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

х, у и г :

М * . у, г); а„(х, у, г); аг(х, у, г) —

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

1.Потенциальный вектор. Вектор в называется потенциаль­

ным, если он является градиентом некоторой скалярной функции

? (х, у, г) (а ~ §гас! <р). Поле потенциального вектора а = §габ <р на­ зывается потенциальным, а скалярная функция ? — потенциалом этого поля. В дальнейшем предполагается, что функция <р и ее частные производные до второго порядка включительно непрерывны.

Необходимым и достаточным условием потенциальности вектора а является выполнение равенств:

(12,1)

даи дах ~3х &

Пример. Потенциальным векторным полем является поле силы

Т — ^ г тяготения, вызываемого

материальной точкой

массы т,

помещенной в начале

координат.

Здесь г У хг + уг +

г* — рас­

стояние произвольной

точки пространства А (х, у, г) с массой /Пх=

= 1 до начала координат (г Ф 0). Потенциалом этого поля является

функция о = — (см. задачу 11,3).

2.Циркуляция вектора и линейный интеграл. Циркуляцией Гх.

вектора а =

а (х, у, г) по замкнутой линии I, называется интеграл

вида

ч .

 

 

Г,. = ^<а хйх + ауйу-\-агйг,

(1 2 ,2)

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

Если Ь — незамкнутая кривая, то интеграл в формуле (12,2) называется линейным интегралом вектораЪ и обозначается буквой и

и =

ахёх + ауйу -\-ахйг

(12,3)

 

(АВ)

 

или в векторной форме

и = $ а>0г,

(12,4)

(АВ)

где — дифференциал радиуса-вектора точки, движущейся по кри­ вой АВ

йг — йх ‘ 7 + Ау • / +

йг к.

Если вектор а — сила, то формула

(12,3) определяет работу

этой силы при перемещении точки по

дуге АВ.

3.Векторная линия. Векторной линией векторного поля вектора

аназывается кривая, в каждой точке которой касательная совпа­

дает с направлением вектора а. Через каждую точку Р векторного поля вектора а проходит по одной векторной линии.

Дифференциальные уравнения векторных линий записываются так:

 

й х

Л у _ й г

(12.5)

 

Ох

оу ~ Ог

4.

Вихрь вектора. Вихрем

вектора,

или ротором вектора а

= а (х,

у, г), обозначаемым

го! а,

называется вектор, определяе­

мый формулой

 

 

 

 

го1 а —

 

 

(12,6)

2В7

Проекции этого вектора на координатные оси равны

го1'0 =

;

 

(12.7)

. - _ д а

дах

Го1*а “ а

— -%•

Модуль вектора го! а определяется формулой

+ ( * - & + & - % ) '■ <вд

Формулу (12,6) можно записать в виде, удобном для запоми­ нания

 

1

]

к

 

 

 

 

 

го1 а =

д_

д

д

 

 

 

(12,9)

д х

5 у

Т

 

 

 

 

 

 

 

 

 

 

ах

ау

<

 

 

 

 

 

п

 

д

,

д

 

д

ау,

аг

При этом произведения символов ^

щ

 

и ^ на а,,

следует понимать так:

 

 

 

 

 

 

 

 

д . а __ дау.

д

а

_

дах

 

 

Т х а « ~ ' З Г ’

Т у

 

-

- 3 7 -

 

 

и т. д.

 

 

 

 

 

 

 

 

Задача 12 ,1 . Известно, что безвихревым

движением жидкости

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

точке жидкости

вы­

полняется условие

 

 

 

 

 

 

 

 

го1 о = 0,

 

 

 

 

 

 

где вектор о = о (дс, у, г) — скорость жидкости в рассматриваемой точке.

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

Р е ш е н и е . Обозначим проекции скорости о на оси прямо­

угольной системы координат через ох, уу, о2. Равенств*) го! о= 0 на основании формулы (12 ,6) можно переписать в виде

го1 о

Из этого следует, что

до.

до,

дох

д о .

д ° у

дох

Выполнение этих равенств является необходимым и достаточ­

ным условием для того,

чтобы

вектор о

был потенциальным —

см. определение на стр. 286 и формулы ( 12, 1).

 

Поскольку вектор о потенциальный, то

существует

такая ска­

лярная функция 9 = <р(х,

г)

потенциал скорости,

что

V §габ 9

(см. определение на стр. 263).

Отсюда

т. е. если

скорость о — потенциальный вектор,

а функция

9 =

= 9 (х, у,

г ) — потенциал скорости, то проекция

скорости на

коор­

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

производной

от потенциала скорости

по соответствующей

этой оси

координате.

 

 

 

 

 

Из (А) следует

 

 

 

 

V,

йх + Оуйу + огйг = д^йх + щду +

 

 

Линейный интеграл (12,3) по кривой АВ от

выражения охйх+

+ ОуЛу + огйг

 

 

 

 

5

ох4х + О у й у + огйг =

\

<*9 9 (Я) — 9 (Л))

(В)

(А В )

(А В )

 

 

 

это значит,

что он не зависит от вида

кривой,

а только

от коор­

динат ее концов А и В, так как 9 (А) н 9 (В) — значение потенциала 9 в точках А и В (предполагается однозначность потенциала ско­

рости —

функции 9).

 

по кривой АВ от

 

Таким образом, линейный интеграл

потен­

циального вектора равен разности значений потенциала

9 в кон­

цах кривой, вдоль которой ведется интегрирование.

 

 

Если

кривая АВ— замкнутая, то линейный интеграл

в

равен­

стве (В)

есть циркуляция

скорости, а

9 (В) — 9 (А) = 0,

так как

точки А и Я в случае замкнутости кривой АВ совпадут.

 

любой

Отсюда следует вывод,

что циркуляция скорости по

замкнутой кривой при наличии однозначного потенциала равна

нут, т. е.. если о * &гас! 9 и 9 — однозначная функция, то

ф'охйх + оус1у + огйг 0.

(Предположение об однозначности функции 9 является весьма суще­ ственным. Если 9 — неоднозначная функция, то сделанный вывод не всегда верен).

10 И. А. Капля*

289

З а м е ч а н и е . Следует иметь в виду, что теорема о равенстве нулю криволинейного интеграла от потенциального вектора по замк­ нутому контуру вёрна только тогда, когда пространство односвязно.

Задача 12 ,2 . Доказать, что если V— потенциальный вектор (о = = §гас!<р), то вихрь этого вектора равен нулю, т. е. что

 

 

го1(§гад 9) = 0.

 

Р е ш е н и е .

В предыдущей задаче было показано, что из равен­

ства го1 о = О

следует, что о есть

вектор

потенциальный, т. е.

V — §гас] 9 . Теперь решаем обратную задачу,

т. е. хотим доказать,

что' если

вектор V потенциальный,

то его вихрь равен нулю.

Итак,

дано,

что V= дгаб 9. Отсюда

 

Поскольку при

ограничениях, наложенных на функцию

9 , изме­

нение порядка дифференцирования

не изменяет

величины ее вто­

рой

смешанной

частной

производной, то

 

,

поэтому

° 1 у)х =

0. Аналогично найдем, что

(го! ь)у = 0

и

(го1 о), = 0 , а

отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

го1 о =

0.

 

 

 

 

Так

как

о = §гас!9 ,

то го1 (§гас1 9) =

0.

 

 

 

 

З а к л ю ч е н и е .

Вихрь

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

равен

нулю.

Задача 12,3.

Проекции

[ускорения

частицы

жидкости

на оси

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

 

Доказать, что при существовании потенциала скорости уско­

рение ТР также

является потенциальным

вектором.

Р е ш е н и е .

По условию о = дгас!9 ,

где 9 — потенциал ско­

рости. Отсюда

 

 

Соседние файлы в папке книги