2.2.5. Нормальный вектор прямой
Или вектор нормали.
Что такое нормаль? Простыми словами, нормаль – это перпендикуляр. То есть, вектор нормали прямой перпендикулярен данной прямой. Очевидно, что у любой прямой их бесконечно много (так же, как и направляющих векторов), но нам хватит одного:
Если прямая задана общим уравнением в декартовой системе координат, то вектор
является вектором нормали данной прямой.
Обратите внимание, что это утверждение справедливо лишь для «школьной» системы координат; все предыдущие выкладки п. 2.2 работают и в общем аффинном случае.
Вектор нормали всегда ортогонален направляющему вектору прямой. Убедимся в ортогональности данных векторов с помощью скалярного произведения:
И тут всё ещё проще: если координаты направляющего вектора приходилось аккуратно «вытаскивать» из уравнения, то координаты вектора нормали
достаточно просто «снять».
Приведу примеры с теми же уравнениями, что и для направляющего вектора:
Можно ли составить уравнение прямой, зная одну точку и вектор нормали? Нутром чувствуется, можно. Ведь вектор нормали ортогонален направляющему вектору и образует с ним «жесткую конструкцию».
Нормальный вектор прямой, координаты нормального вектора прямой
Для изучения уравнений прямой линии необходимо хорошо разбираться в алгебре векторов. Важно нахождение направляющего вектора и нормального вектора прямой. В данной статье будут рассмотрены нормальный вектор прямой с примерами и рисунками, нахождение его координат, если известны уравнения прямых. Будет рассмотрено подробное решение.
Нормальный вектор прямой – определение, примеры, иллюстрации
Чтобы материал легче усваивался, нужно разбираться в понятиях линия, плоскость и определениями, которые связаны с векторами. Для начала ознакомимся с понятием вектора прямой.
Нормальным вектором прямой называют любой ненулевой вектор, который лежит на любой прямой, перпендикулярной данной.
Понятно, что имеется бесконечное множество нормальных векторов, расположенных на данной прямой. Рассмотрим на рисунке, приведенном ниже.
Получаем, что прямая является перпендикулярной одной из двух заданных параллельных прямых, тогда ее перпендикулярность распространяется и на вторую параллельную прямую. Отсюда получаем, что множества нормальных векторов этих параллельных прямых совпадают. Когда прямые a и а 1 параллельные, а n → считается нормальным вектором прямой a , также считается нормальным вектором для прямой a 1 . Когда прямая а имеет прямой вектор, тогда вектор t · n → является ненулевым при любом значении параметра t , причем также является нормальным для прямой a .
Используя определение нормального и направляющего векторов, можно прийти к выводу, что нормальный вектор перпендикулярен направляющему. Рассмотрим пример.
Если задана плоскость О х у , то множеством векторов для О х является координатный вектор j → . Он считается ненулевым и принадлежащим координатной оси О у , перпендикулярной О х . Все множество нормальных векторов относительно О х можно записать, как t · j → , t ∈ R , t ≠ 0 .
Прямоугольная система O x y z имеет нормальный вектор i → , относящийся к прямой О z . Вектор j → также считается нормальным. Отсюда видно, что любой ненулевой вектор, расположенный в любой плоскости и перпендикулярный О z , считается нормальным для O z .
Координаты нормального вектора прямой – нахождение координат нормального вектора прямой по известным уравнениям прямой
При рассмотрении прямоугольной системы координат О х у выявим, что уравнение прямой на плоскости соответствует ей, а определение нормальных векторов производится по координатам. Если известно уравнение прямой, а необходимо найти координаты нормального вектора, тогда необходимо из уравнения A x + B y + C = 0 выявить коэффициенты, которые и соответствуют координатам нормального вектора заданной прямой.
Задана прямая вида 2 x + 7 y — 4 = 0 _, найти координаты нормального вектора.
По условию имеем, что прямая была задана общим уравнением, значит необходимо выписать коэффициенты , которые и являются координатами нормального вектора. Значит, координаты вектора имеют значение 2 , 7 .
Бывают случаи, когда A или В из уравнения равняется нулю. Рассмотрим решение такого задания на примере.
Указать нормальный вектор для заданной прямой y — 3 = 0 .
По условию нам дано общее уравнение прямой, значит запишем его таким образом 0 · x + 1 · y — 3 = 0 . Теперь отчетливо видим коэффициенты, которые и являются координатами нормального вектора. Значит, получаем, что координаты нормального вектора равны 0 , 1 .
Если дано уравнение в отрезках вида x a + y b = 1 или уравнение с угловым коэффициентом y = k · x + b , тогда необходимо приводить к общему уравнению прямой, где можно найти координаты нормального вектора данной прямой.
Найти координаты нормального вектора, если дано уравнение прямой x 1 3 — y = 1 .
Для начала необходимо перейти от уравнения в отрезках x 1 3 — y = 1 к уравнению общего вида. Тогда получим, что x 1 3 — y = 1 ⇔ 3 · x — 1 · y — 1 = 0 .
Отсюда видно, что координаты нормального вектора имеют значение 3 , — 1 .
Ответ: 3 , — 1 .
Если прямая определена каноническим уравнением прямой на плоскости x — x 1 a x = y — y 1 a y или параметрическим x = x 1 + a x · λ y = y 1 + a y · λ , тогда получение координат усложняется. По данным уравнениям видно, что координаты направляющего вектора будут a → = ( a x , a y ) . Возможность нахождения координат нормального вектора n → возможно, благодаря условию перпендикулярности векторов n → и a → .
Имеется возможность получения координат нормального вектора при помощи приведения канонического или параметрического уравнений прямой к общему. Тогда получим:
x — x 1 a x = y — y 1 a y ⇔ a y · ( x — x 1 ) = a x · ( y — y 1 ) ⇔ a y · x — a x · y + a x · y 1 — a y · x 1 x = x 1 + a x · λ y = y 1 + a y · λ ⇔ x — x 1 a x = y — y 1 a y ⇔ a y · x — a x · y + a x · y 1 — a y · x 1 = 0
Для решения можно выбирать любой удобный способ.
Найти нормальный вектор заданной прямой x — 2 7 = y + 3 — 2 .
Из прямой x — 2 7 = y + 3 — 2 понятно, что направляющий вектор будет иметь координаты a → = ( 7 , — 2 ) . Нормальный вектор n → = ( n x , n y ) заданной прямой является перпендикулярным a → = ( 7 , — 2 ) .
Выясним, чему равно скалярное произведение. Для нахождения скалярного произведения векторов a → = ( 7 , — 2 ) и n → = ( n x , n y ) запишем a → , n → = 7 · n x — 2 · n y = 0 .
Значение n x – произвольное , следует найти n y . Если n x = 1 , отсюда получаем, что 7 · 1 — 2 · n y = 0 ⇔ n y = 7 2 .
Значит, нормальный вектор имеет координаты 1 , 7 2 .
Второй способ решения сводится к тому, что необходимо прийти к общему виду уравнения из канонического. Для этого преобразуем
x — 2 7 = y + 3 — 2 ⇔ 7 · ( y + 3 ) = — 2 · ( x — 2 ) ⇔ 2 x + 7 y — 4 + 7 3 = 0
Полученный результат координат нормального вектора равен 2 , 7 .
Ответ: 2 , 7 или 1 , 7 2 .
Указать координаты нормального вектора прямой x = 1 y = 2 — 3 · λ .
Для начала необходимо выполнить преобразование для перехода в общему виду прямой. Выполним:
x = 1 y = 2 — 3 · λ ⇔ x = 1 + 0 · λ y = 2 — 3 · λ ⇔ λ = x — 1 0 λ = y — 2 — 3 ⇔ x — 1 0 = y — 2 — 3 ⇔ ⇔ — 3 · ( x — 1 ) = 0 · ( y — 2 ) ⇔ — 3 · x + 0 · y + 3 = 0
Отсюда видно, что координаты нормального вектора равны — 3 , 0 .
Рассмотрим способы для нахождения координат нормального вектора при уравнении прямой в пространстве, заданной прямоугольной системой координат О х у z .
Когда прямая задается при помощи уравнений пересекающихся плоскостей A 1 x + B 1 y + C 1 z + D 1 = 0 и A 2 x + B 2 y + C 2 z + D 2 = 0 , тогда нормальный вектор плоскости относится к A 2 x + B 2 y + C 2 z + D 2 = 0 и A 2 x + B 2 y + C 2 z + D 2 = 0 , тогда получаем запись векторов в виде n 1 → = ( A 1 , B 1 , C 1 ) и n 2 → = ( A 2 , B 2 , C 2 ) .
Когда прямая определена при помощи канонического уравнения пространства, имеющего вид x — x 1 a x = y — y 1 a y = z — z 1 a z или параметрического, имеющего вид x = x 1 + a x · λ y = y 1 + a y · λ z = z 1 + a z · λ , отсюда a x , a y и a z считаются координатами направляющего вектора заданной прямой. Любой ненулевой вектор может быть нормальным для данной прямой, причем являться перпендикулярным вектору a → = ( a x , a y , a z ) . Отсюда следует, что нахождение координат нормального с параметрическими и каноническими уравнениями производится при помощи координат вектора, который перпендикулярен заданному вектору a → = ( a x , a y , a z ) .
Нормальный вектор прямой
Вы будете перенаправлены на Автор24
В аналитической геометрии часто требуется составить общее уравнение прямой по принадлежащей ей точке и вектору нормали к прямой.
Нормаль – синоним для слова перпендикуляр.
Общее уравнение прямой на плоскости выглядит как $Ax + By + C = 0$. Подставляя в него различные значениях $A$, $B$ и $C$, в том числе нулевые, можно определить любые прямые.
Можно выразить уравнение прямой и другим способом:
Это уравнение прямой с угловым коэффициентом. В нем геометрический смысл коэффициента $k$ заключается в угле наклона прямой по отношению к оси абсцисс, а независимого члена $b$ — в расстоянии, на которое прямая отстоит от центра координатной плоскости, т.е. точки $O(0; 0)$.
Рисунок 1. Варианты расположения прямых на координатной плоскости. Автор24 — интернет-биржа студенческих работ
Нормальное уравнение прямой можно выразить и в тригонометрическом виде:
$x cdot cos <alpha>+ y cdot sin <alpha>- p = 0$
где $alpha$ — угол между прямой и осью абсцисс, а $p$ — расстояние от начала координат до рассматриваемой прямой.
Возможны четыре варианта зависимости наклона прямой от величины углового коэффициента:
- когда угловой коэффициент положителен, направляющий вектор прямой идёт снизу вверх;
- когда угловой коэффициент отрицателен, направляющий вектор прямой идёт сверху вниз;
- когда угловой коэффициент равен нулю, описываемая им прямая параллельна оси абсцисс;
- для прямых, параллельных оси ординат, углового коэффициента не существует, поскольку тангенс 90 градусов является неопределенной (бесконечной) величиной.
Готовые работы на аналогичную тему
Чем больше абсолютное значение углового коэффициента, тем круче наклонен график прямой.
Зная угловой коэффициент, легко составить уравнение графика прямой, если дополнительно известна точка, принадлежащая искомой прямой:
$y — y_0 = k cdot (x — x_0)$
Таким образом, геометрически прямую на координатной всегда можно выразить с помощью угла и расстояния от начала координат. В этом и заключается смысл нормального вектора к прямой — самого компактного способа записи ее положения, если известны координаты хотя бы одной точки, принадлежащей этой прямой.
Вектором нормали к прямой, иначе говоря, нормальным вектором прямой, принято называть ненулевой вектор, перпендикулярный рассматриваемой прямой.
Для каждой прямой можно найти бесконечное множество нормальных векторов, равно как и направляющих векторов, т.е. таких, которые параллельны этой прямой. При этом все нормальные векторы к ней будут коллинеарными, хотя и не обязательно сонаправлены.
Обозначив нормальный вектор прямой как $vec(n_1; n_2)$, а координаты точки как $x_0$ и $y_0$, можно представить общее уравнение прямой на плоскости по точке и вектору нормали к прямой как
$n_1 cdot (x — x_n) + n_2 cdot (y — y_0) = 0$
Таким образом, координаты вектора нормали к прямой пропорциональны числам $A$ и $B$, присутствующим в общем уравнении прямой на плоскости. Следовательно, если известно общее уравнение прямой на плоскости, то можно легко вывести и вектор нормали к прямой. Если прямая, задана уравнением в прямоугольной системе координат
то нормальный вектор описывается формулой:
При этом говорят, что координаты нормального вектора «снимаются» с уравнения прямой.
Нормальный к прямой вектор и ее направляющий вектор всегда ортогональны по отношению друг к другу, т.е. их скалярные произведения равны нулю, в чем легко убедиться, вспомнив формулу направляющего вектора $bar
(-B; A)$, а также общее уравнение прямой по направляющему вектору $bar
(p_1; p_2)$ и точке $M_0(x_0; y_0)$:
В том, что вектор нормали к прямой всегда ортогонален направляющему вектору к ней можно убедиться с помощью скалярного произведения:
$bar
cdot bar = -B cdot A + A cdot B = 0 implies bar
perp bar$
Всегда можно составить уравнение прямой, зная координаты принадлежащей ей точки и нормального вектора, поскольку направление прямой следует из его направления. Описав точку как $M(x_0; y_0)$, а вектор как $bar(A; B)$, можно выразить уравнение прямой в следующем виде:
$A(x — x_0) + B(y — y_0) = 0$
Составить уравнение прямой по точке $M(-1; -3)$ и нормальному вектору $bar(3; -1)$. Вывести уравнение направляющего вектора.
Для решения задействуем формулу $A cdot (x — x_0) + B cdot (y — y_0) = 0$
Подставив значения, получаем:
$3 cdot (x — (-1)) — (-1) cdot (y — (-3)) = 0$ $3 cdot (x + 1) — (y + 3) = 0$ $3x + 3 — y — 3 = 0$ $3x — y = 0$
Проверить правильность общего уравнения прямой можно «сняв» из него координаты для нормального вектора:
$3x — y = 0 implies A = 3; B = -1 implies bar(A; B) = bar(3; -1),$
Что соответствует числам исходных данных.
Подставив реальные значения, проверим, удовлетворяет ли точка $M(-1; -3)$ уравнению $3x — y = 0$:
Равенство верно. Осталось лишь найти формулу направляющего вектора:
$bar
(-B; A) implies bar
(1; 3)$
Ответ: $3x — y = 0; bar
(1; 3).$
Получи деньги за свои студенческие работы
Курсовые, рефераты или другие работы
Автор этой статьи Дата последнего обновления статьи: 04 03 2022
http://zaochnik.com/spravochnik/matematika/prjamaja-ploskost/normalnyj-vektor-prjamoj-koordinaty-normalnogo-vek/
http://spravochnick.ru/matematika/normalnyy_vektor_pryamoy/
Векторы в пространстве и метод координат
Существует два способа решения задач по стереометрии
Первый — классический — требует отличного знания аксиом и теорем стереометрии, логики, умения построить чертеж и свести объемную задачу к планиметрической. Способ хорош тем, что развивает мозги и пространственное воображение.
Другой метод — применение векторов и координат. Это простые формулы, алгоритмы и правила. Он очень удобен, особенно когда времени до экзамена мало, а решить задачу хочется.
Если вы освоили векторы на плоскости и действия с ними — то и с векторами в пространстве разберетесь. Многие понятия окажутся знакомыми.
Система координат в пространстве
Выберем начало координат. Проведем три взаимно перпендикулярные оси X, Y и Z. Зададим удобный масштаб.
Получилась система координат в трехмерном пространстве. Теперь каждая его точка характеризуется тремя числами — координатами по X, Y и Z. Например, запись M(−1; 3; 2) означает, что координата точки M по X (абсцисса) равна −1, координата по Y (ордината) равна 3, а координата по Z (аппликата) равна 2.
Векторы в пространстве определяются так же, как и на плоскости. Это направленные отрезки, имеющие начало и конец. Только в пространстве вектор задается тремя координатами x, y и z:
Как найти координаты вектора? Как и на плоскости — из координаты конца вычитаем координату начала.
Длина вектора в пространстве – это расстояние между точками A и B. Находится как корень квадратный из суммы квадратов координат вектора:
Пусть точка M – середина отрезка AB. Ее координаты находятся по формуле:
Для сложения векторов применяем уже знакомые правило треугольника и правило параллелограмма
Сумма векторов, их разность, произведение вектора на число и скалярное произведение векторов определяются так же, как и на плоскости. Только координат не две, а три. Возьмем векторы и
.
Сумма векторов:
Разность векторов:
Произведение вектора на число:
Скалярное произведение векторов:
Косинус угла между векторами:
Последняя формула удобна для нахождения угла между прямыми в пространстве. Особенно если эти прямые – скрещиваются. Напомним, что так называются прямые, которые не параллельны и не пересекаются. Они лежат в параллельных плоскостях.
1. В кубе ABCDA1B1C1D1 точки E и K — середины ребер соответственно A1B1 и B1C1. Найдите косинус угла между прямыми AE и BK.
Если вам достался куб — значит, повезло. Он отлично вписывается в прямоугольную систему координат. Строим чертеж:
Длина ребра куба не дана. Какой бы она ни была, угол между AE и BK от нее не зависит. Поэтому возьмем единичный куб, все ребра которого равны 1.
Прямые AE и BK — скрещиваются. Найдем угол между векторами и
. Для этого нужны их координаты.
Запишем координаты векторов:
и найдем косинус угла между векторами и
:
2. В правильной четырехугольной пирамиде SABCD, все ребра которой равны 1, точки E, K — середины ребер SB и SC соответственно. Найдите косинус угла между прямыми AE и BK.
Лучше всего выбрать начало координат в центре основания пирамиды, а оси X и Y сделать параллельными сторонам основания.
Координаты точек A, B и C найти легко:
Из прямоугольного треугольника AOS найдем
Координаты вершины пирамиды:
Точка E — середина SB, а K — середина SC. Воспользуемся формулой для координат середины отрезка и найдем координаты точек E и K.
Найдем координаты векторов и
:
и угол между ними:
Покажем теперь, как вписать систему координат в треугольную призму.
3. В правильной треугольной призме ABCA1B1C1, все ребра которой равны 1, точка D — середина ребра A1B1. Найдите косинус угла между прямыми AD и BC1
Пусть точка A — начало координат. Возьмем ось X параллельно стороне BC, а ось Y перпендикулярно ей. Другими словами, на оси Y будет лежать отрезок AH, являющийся высотой треугольника ABC. Нарисуем отдельно нижнее основание призмы.
Запишем координаты точек:
Точка D — середина A1B1. Значит, пользуемся формулами для координат середины
отрезка.
Найдем координаты векторов и
, а затем угол между ними:
Смотрите, как легко с помощью векторов и координат найти угол между прямыми. А если требуется найти угол между плоскостями или между прямой и плоскостью? Для решения подобных задач нам понадобится уравнение плоскости в пространстве.
Плоскость в пространстве задается уравнением:
Здесь числа A, B и C — координаты вектора, перпендикулярного этой плоскости. Его называют нормалью к плоскости.
Вместо x, y и z можно подставить в уравнение координаты любой точки, принадлежащей данной плоскости. Получится верное равенство.
Плоскость в пространстве можно провести через любые три точки, не лежащие на одной прямой. Поэтому для того, чтобы написать уравнение плоскости, берем координаты трех принадлежащих ей точек. Подставляем их по очереди в уравнение плоскости. Решаем полученную систему.
Покажем, как это делается.
Напишем уравнение плоскости, проходящей через точки M (1; 0; 1), N (2; −2; 0) и K (4; 1; 2).
Уравнение плоскости выглядит так:
Подставим в него по очереди координаты точек M, N и K.
Для точки M:
То есть A + C + D = 0.
Для точки N:
Аналогично для точки K:
Получили систему из трех уравнений:
.
В ней четыре неизвестных: A, B, C и D. Поэтому одну из них мы выберем сами, а другие выразим через нее. Правило простое — вместо одной из переменных можно взять любое число, не равное нулю.
Пусть, например, D = −2. Тогда:
;
.
Выразим C и B через A и подставим в третье уравнение:
.
Решив систему, получим:
Уравнение плоскости MNK имеет вид:
Умножим обе части уравнения на −3. Тогда коэффициенты станут целыми:
Вектор — это нормаль к плоскости MNK.
Уравнение плоскости, проходящей через заданную точку имеет вид:
Угол между плоскостями равен углу между нормалями к этим плоскостям:
Не правда ли, знакомая формула? Скалярное произведение нормалей поделили на произведение их длин.
Заметим, что при пересечении двух плоскостей вообще-то образуется четыре угла.
Мы берем меньший из них. Поэтому в формуле стоит модуль скалярного произведения — чтобы косинус угла был неотрицателен.
4. В кубе ABCDA1B1C1D1 точки E и F — середины ребер соответственно A1B1 и A1D1. Найдите тангенс угла между плоскостями AEF и BDD1.
Строим чертеж. Видно, что плоскости AEF и BDD1 пересекаются где-то вне куба. В классическом решении пришлось бы строить линию их пересечения. Но векторно-координатный метод значительно всё упрощает. Не будем ломать голову над тем, по какой прямой пересекаются плоскости. Просто отметим координаты нужных нам точек и найдем угол между нормалями к плоскостям AEF и BDD1.
Сначала — нормаль к плоскости BDD1. Конечно, мы можем подставить координаты точек B, D и D1 в уравнение плоскости и найти коэффициенты, которые и будут координатами вектора нормали. А можем сделать хитрее — увидеть нужную нормаль прямо на чертеже. Ведь плоскость BDD1 — это диагональное сечение куба. Вектор перпендикулярен этой плоскости.
Итак, первый вектор нормали у нас уже есть:
Напишем уравнение плоскости AEF.
Берем уравнение плоскости и по очереди подставляем в него, вместо x, y и z, соответствующие координаты точек A, E и F.
Упростим систему:
.
Пусть С = -1. Тогда A = B = 2.
Уравнение плоскости AEF:
Нормаль к плоскости AEF:
Найдем угол между плоскостями:
5. Основание прямой четырехугольной призмы BCDA1B1C1D1 — прямоугольник ABCD, в котором AB = 5, AD = √33. Найдите тангенс угла между плоскостью грани AA1D1D и плоскостью, проходящей через середину ребра CD перпендикулярно прямой B1D, если расстояние между прямыми A1C1 и BD равно √3.
Эта задача наглядно показывает, насколько векторный метод проще классического. Попробуйте, для разнообразия, построить необходимые сечения и провести все доказательства — как это делается в «классике»
Строим чертеж. Прямую четырехугольную призму можно по-другому назвать «параллелепипед».
Замечаем, что длина и ширина параллелепипеда у нас есть, а вот высота — вроде не дана. Как же ее найти?
«Расстояние между прямыми A1C1 и BD равно √3». Прямые A1C1 и BD скрещиваются. Одна из них — диагональ верхнего основания, другая — диагональ нижнего. Вспомним, что расстояние между скрещивающимися прямыми равно длине их общего перпендикуляра. Общий перпендикуляр к A1C1 и BD — это, очевидно, OO1, где O — точка пересечения диагоналей нижнего основания, O1 — точка пересечения диагоналей верхнего. А отрезок OO1 и равен высоте параллелепипеда.
Итак, AA1 = √3
Плоскость AA1 D1 D — это задняя грань призмы на нашем чертеже. Нормаль к ней — это любой вектор, перпендикулярный задней грани, например, вектор или, еще проще, вектор
.
Осталась еще «плоскость, проходящая через середину ребра CD перпендикулярно прямой B1D». Но позвольте, если плоскость перпендикулярна прямой B1D — значит, B1D и есть нормаль к этой плоскости! Координаты точек B1 и D известны:
Координаты вектора — тоже:
Находим угол между плоскостями, равный углу между нормалями к ним:
Зная косинус угла, находим его тангенс по формуле
Получим:
Ответ:
Угол между прямой m и плоскостью α тоже вычисляется с помощью скалярного произведения векторов.
Пусть — вектор, лежащий на прямой m (или параллельный ей),
— нормаль к плоскости α.
Находим синус угла между прямой m и плоскостью α по формуле:
6. В кубе ABCDA1B1C1D1 точка E — середина ребра A1B1. Найдите синус угла между прямой AE и плоскостью BDD1.
Как всегда, рисуем чертеж и выбираем систему координат
Находим координаты вектора .
Нужно ли нам уравнение плоскости BDD1? В общем-то, без него можно обойтись. Ведь эта плоскость является диагональным сечением куба, а значит, нормалью к ней будет любой вектор, ей перпендикулярный. Например, вектор .
Найдем угол между прямой и плоскостью:
Ответ:
Расстояние от точки M с координатами x0, y0 и z0 до плоскости α, заданной уравнением Ax + By + Cz + D = 0, можно найти по формуле:
7. В основании прямоугольного параллелепипеда BCDA1B1C1D1 лежит прямоугольник ABCD со сторонами AB = , AD =
. Высота параллелепипеда AA1 =
. Найдите расстояние от точки A до плоскости A1DB.
Построим чертеж и выпишем координаты точек:
Запишем уравнение плоскости A1DB. Вы помните, как это делается — по очереди подставляем координаты точек A1, D и B в уравнение Ax + Be + Cz + D
Решим эту систему. Выберем
Тогда
Уравнение плоскости A1DB имеет вид:
Дальше все просто. Находим расстояние от точки A до плоскости A1DB:
В некоторых задачах по стереометрии требуется найти расстояние от прямой до параллельной ей плоскости. В этом случае можно выбрать любую точку, принадлежащую данной прямой.
Спасибо за то, что пользуйтесь нашими статьями.
Информация на странице «Векторы в пространстве и метод координат» подготовлена нашими авторами специально, чтобы помочь вам в освоении предмета и подготовке к ЕГЭ и ОГЭ.
Чтобы успешно сдать нужные и поступить в высшее учебное заведение или колледж нужно использовать все инструменты: учеба, контрольные, олимпиады, онлайн-лекции, видеоуроки, сборники заданий.
Также вы можете воспользоваться другими статьями из данного раздела.
Публикация обновлена:
08.05.2023
Вычисление координат вектора нормали. Рассматривая модели отражения света, вы,
наверное, обратили внимание
на то, что нормаль к поверхности — важный элемент. Определение вектора нормали
к поверхности в заданной точке может быть выполнено разными способами. В значительной степени это определяется типом модели
описания поверхности. Для поверхностей,
заданных в аналитической форме, известны методы дифференциальной геометрии,
которые основываются на вычислении частных производных функций описания. Например, если поверхность задана
параметрическими функциями
тогда
координаты вектора нормали можно вычислить так:
В
случае описания поверхности векторно-полигональной моделью для определения
нормалей можно использовать методы векторной алгебры.
Пусть в пространстве задана некоторая многогранная
поверхность. Рассмотрим одну ее плоскую грань, имеющую вид
треугольника (рис. 8.4).
Для
вычисления координат вектора нормали воспользуемся векторным произведением
любых двух векторов, которые лежат в плоскости грани. Такими векторами могут
служить и ребра грани, например, ребра 1-2 и 1-3. Однако формулы для векторного
произведения были определены нами только для
радиус-векторов. Чтобы перейти к радиус-векторам, введем новую систему
координат, центр которой совпадает с вершиной 1, а оси — параллельны осям бывшей системы. Координаты вершин в новой
системе:
Рис. 8.4. Одна грань поверхности
Радиус-векторы
Теперь назовем
ребро (1-2) вектором А, а ребро (1-3) — вектором В, как
показано на рис. 8.4. Таким образом, положение нормали к грани в пространстве будет
описываться радиус-вектором N. Его координаты в системе (х’, у’, г’) выразим
формулами для векторного произведения
Плоская грань
может быть изображена в разных ракурсах. В каждой конкретной ситуации необходимо выбирать направление нормали,
которое соответствует видимой стороне грани. Если плоская грань может быть видна с обратной
стороны, то тогда в расчетах отраженного света необходимо выбирать для нормали обратный
вектор, то есть (-N).
Если
полигональная поверхность имеет не треугольные грани, а, например, плоские четырехугольные, то расчет нормали можно выполнять
по любым трем вершинам грани.
В этом уроке я обьясню вам, как находить нормали полигонов.
Все математические функции мы вынесем в отдельный класс: 3dmath.h/.cpp
Отныне вы сможете простым подключением этих файлов использовать их функции в любых своих проектах.
Файлы 3dmath.h/.cpp будут содержать математические функции, не только для этого урока, но и на будущее. Математика — необходимая часть програмирования игр, и в особенности трёхмерных. Если вы не чувствуете себя уверенно в основах алгебры и геометрии, займитесь ими, пока не начнете ориентироватся. Вы не сможете стать эффективным 3д программистом без знания этих основ. Однако не гордитесь, если вы их знаете. Остается всё ещё очень много вещей помимо основ 😉
Итак, новый файл 3dmath.h
#ifndef _3DMATH_H
#define _3DMATH_H// Эта структура используется для хранения 3д точек и векторов. Она использовалась
// в наших уроках по камере, вернитесь назад, чтобы узнать подробности.
struct CVector3
{
public:
float x, y, z;
};
// Возвращает вектор, перпендикулярный двум переданным векторам (плоскости)
CVector3 Cross(CVector3 vVector1, CVector3 vVector2);
// Возвращает вектор между 2мя точками.
CVector3 Vector(CVector3 vPoint1, CVector3 vPoint2);
// Возвращает величину нормали или любого другого вектора
float Magnitude(CVector3 vNormal);
// Возвращает нормализованный вектор (его длинна становится равной 1)
CVector3 Normalize(CVector3 vNormal);
// Возвращает нормаль полигона (направление, куда повернут полигон)
CVector3 Normal(CVector3 vTriangle[]);
#endif
Файл 3dmath.cpp:
#include «3dmath.h»
#include <math.h>// * Нахождение нормали полигона *
// Чтобы найти нормаль полигона, нам нужно найти результат cross-a от двух
// векторов этого полигона. В общем, это всё, что нам нужно для получения направлений
// двух сторон треугольника. В конце концов, вектор — это только направление и длинна.
// Длинна вектора в нашем случае не важна. Нам нужно только узнать направление.
// Итак, имея 2 вектора треугольника, мы можем найти вектор, стоящий перпендикулярно
// к полигону.
// Теперь, в зависимости от порядка следования вершин, нормаль будет расположена с
// какой-то из сторон полигона. Вам остаётся только решить, в каком порядке отрисовывать
// вершины — ВСЕГДА запоминайте это.
// Обычно полигоны отрисовываются только с одной стороны. Никому не нужно второй раз
// отрисовывать сторону, которую не видно. Задумайтесь, если у вас есть какая-нибуть 3д
// модель, нужно ли вам отрисовывать внутренние стороны её полигонов? Конечно нет,
// это безсмысленно.
//
/////////////////////////////////////// CROSS \\*
/////
///// Возвращает вектор, перпендикулярный 2м переданным.
/////
/////////////////////////////////////// CROSS \\*
CVector3 Cross(CVector3 vVector1, CVector3 vVector2)
{
CVector3 vNormal; // результирующий вектор
// Еще раз, если у нас есть 2 вектора (2 стороны полигона), у нас есть плоскость.
// cross находит вектор, перпендикулярный плоскости, составляемой 2мя векторами.
// Формула в принципе проста, но сложна для запоминания:
// Значение X для векторы вычисляется так: (V1.y * V2.z) — (V1.z * V2.y)
vNormal.x = ((vVector1.y * vVector2.z) — (vVector1.z * vVector2.y));
// Значение Y для векторы вычисляется так: (V1.z * V2.x) — (V1.x * V2.z)
vNormal.y = ((vVector1.z * vVector2.x) — (vVector1.x * vVector2.z));
// Значение Z для векторы вычисляется так: (V1.x * V2.y) — (V1.y * V2.x)
vNormal.z = ((vVector1.x * vVector2.y) — (vVector1.y * vVector2.x));
return vNormal; // Возвращаем результат (направление, куда направлен полигон — нормаль)
}
/////////////////////////////////////// VECTOR \\*
/////
///// Возвращает вектор между 2мя точками.
/////
/////////////////////////////////////// VECTOR \\*
CVector3 Vector(CVector3 vPoint1, CVector3 vPoint2)
{
CVector3 vVector = {0};
// Чтобы получить вектор между 2 точками (направление), нужно вычесть вторую
// точку из первой.
vVector.x = vPoint1.x — vPoint2.x;
vVector.y = vPoint1.y — vPoint2.y;
vVector.z = vPoint1.z — vPoint2.z;
// Теперь возвращаем полученный результат
return vVector;
}
/////////////////////////////////////// MAGNITUDE \\*
/////
///// возвращает величину нормали
/////
/////////////////////////////////////// MAGNITUDE \\*
float Magnitude(CVector3 vNormal)
{
return (float)sqrt( (vNormal.x * vNormal.x) +
(vNormal.y * vNormal.y) +
(vNormal.z * vNormal.z) );
}
/////////////////////////////////////// NORMALIZE \\*
/////
///// возвращает нормализованный вектор (с длинной 1)
/////
/////////////////////////////////////// NORMALIZE \\*
CVector3 Normalize(CVector3 vNormal)
{
float magnitude = Magnitude(vNormal);
vNormal.x /= magnitude;
vNormal.y /= magnitude;
vNormal.z /= magnitude;
return vNormal;
}
/////////////////////////////////////// NORMAL \\*
/////
///// Возвращает нормаль полигона
/////
/////////////////////////////////////// NORMAL \\*
CVector3 Normal(CVector3 vTriangle[])
{
CVector3 vVector1 = Vector(vTriangle[2], vTriangle[0]);
CVector3 vVector2 = Vector(vTriangle[1], vTriangle[0]);
// В функцию передаются три вектора — треугольник. Мы получаем vVector1 и vVector2 — его
// стороны. Теперь, имея 2 стороны треугольника, мы можем получить из них cross().
// (*ЗАМЕЧАНИЕ*) Важно: первым вектором мы передаём низ треугольника, а вторым — левую
// сторону. Если мы поменяем их местами, нормаль будет повернута в противоположную
// сторону. В нашем случае мы приняли решение всегда работать против часовой.
CVector3 vNormal = Cross(vVector1, vVector2);
// Теперь, имея направление нормали, осталось сделать последнюю вещь. Сейчас её
// длинна неизвестна, она может быть очень длинной. Мы сделаем её равной 1, это
// называется нормализация. Чтобы сделать это, мы делим нормаль на её длинну.
// Ну а как найти длинну? Мы используем эту формулу: magnitude = sqrt(x^2 + y^2 + z^2)
vNormal = Normalize(vNormal);
// Теперь вернём «нормализованную нормаль» =)
// (*ПРИМЕЧАНИЕ*) если вы хотите увидеть, как работает нормализация, закомментируйте
// предидущую линию. Вы увидите, как длинна нормаль до нормалицации. Я стого рекомендую
// всегда использовать эту функцию. И запомните, неважно, какова длинна нормали
// (конечно, кроме (0,0,0)), если мы её нормализуем, она всегда будет равна 1.
return vNormal;
}
Теперь файл main.cpp; мы используем в нём описанные только что функции.
// Для начала после инклуда main.h подключим новый файл:
#include «3dmath.h»// Участок ниже — важен. Мы создаем массив трех векторов и сохраняем в него треугольник.
// Обратите внимание, мы инициализируем сначала левую точку треугольника, потом
// правую, потом верхнюю. Против часовой стрелки. Это важно, надо запомнить этот порядок.
// Иначе начнутся проблемы с направлением наших нормалей. Просто примите такой порядок
// как стандарт.
// Конечно, вы можете при желании изменить направление нормалей
// в OpenGL (см. glFrontFace() glCullFace().
// Разместите эту строку после включения хидеров:
CVector3 vTriangle[3] = { {—1, 0, 0}, {1, 0, 0}, {0, 1, 0} };
// Теперь перейдём к функции WinProc(). Чтобы сделать урок нагляднее,
// напишем обработку дополнительных клавиш:
// Если нажато «LEFT» или «RIGHT», двигаем треугольник
// и нормаль. Заметьте, что нормаль следует за центром треугольника, когда тот
// движется.
case WM_KEYDOWN:
switch(wParam)
{
case VK_ESCAPE:
PostQuitMessage(0);
break;
case VK_LEFT:
vTriangle[2].z -= 0.01f;
vTriangle[2].y -= 0.01f;
break;
case VK_RIGHT:
vTriangle[2].z += 0.01f;
vTriangle[2].y += 0.01f;
break;
}
break;
///////////////////////////////////
// И последний этап, переходим к отрисовке:
void RenderScene()
{
glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
glLoadIdentity();
// Изменяем позицию камеры, чтобы лучше видеть нормаль под углом. Смещаемся немного
// влево и вперёд от треугольника. И приподнимаем, чтобы смотреть немного сверху.
// Position View Up Vector
gluLookAt(—2.5f, 0, —0.5f, 0, 0.5f, 0, 0, 1, 0);
// Ниже рисуем треугольник, передавая данные из массива векторов.
glBegin (GL_TRIANGLES);
glColor3ub(255, 0, 0);
glVertex3f(vTriangle[0].x, vTriangle[0].y, vTriangle[0].z);
glColor3ub(0, 255, 0);
glVertex3f(vTriangle[1].x, vTriangle[1].y, vTriangle[1].z);
glColor3ub(0, 0, 255);
glVertex3f(vTriangle[2].x, vTriangle[2].y, vTriangle[2].z);
glEnd();
// Теперь, отобразив треугольник, давайте отобразим линию, идущую из него под
// углом 90 градусов, чтобы показать, на что похожа нормаль. Запомните, у нормали НЕТ
// позиции, только направление, но просто чтобы визуализировать, нарисуем её исходящей
// из центра треугольника.
// Получим нормаль треугольника, передав его массив.
CVector3 vNormal = Normal(vTriangle);
// Теперь просто для удобства найдем центр треугольника.
// Мы знаем, что середина по X = 0, так что её нам находить не надо.
CVector3 vCenter = {0};
// Ищем центр треугольника.
// Чтобы найти значение Y, просто делим верхнюю точку Y на 2, так как отсчет идёт от 0.
// Сделаем то же для Z, просто делим значение Z верхней точки на 2, т.к. отсчет идет от
// нуля по оси Z, и только верхняя точка треугольника будет двигатся.
// Если бы треугольник не был так просто расположен, пришлось бы сделать всё немного
// сложнее. Чтобы найти центр обьекта, сложите все значения X, Y, и Z, затем
// разделите каждое
// на количество вершин, и это даст вам центр. Например, после получения итога (total):
// total.x /= totalVertices; total.y /=totalVertices; total.z /= totalVertices;
// Я решил не делать отдельной функции для наших рассчетов, но вам — придется 😉
vCenter.y = vTriangle[2].y / 2;
vCenter.z = vTriangle[2].z / 2;
// Давайте нарисуем линию из центра треугольника.
// Первая точка будет центром, который мы только что вычислили.
// Следующая будет центром + нормаль треугольника.
glBegin (GL_LINES);
glColor3ub(255, 255, 0);
// Первая точка — центр треугольника:
glVertex3f(vCenter.x, vCenter.y, vCenter.z);
// Вторая — другой конец нормали:
glVertex3f(vCenter.x + vNormal.x, vCenter.y + vNormal.y, vCenter.z + vNormal.z);
glEnd();
// Это всё, теперь используйте клавиатуру, чтобы увидеть нормаль в движении.
SwapBuffers(g_hDC);
}
Вот и всё! Это основы нахождения нормалей. Ни о чем не приходится заботится, функции, написанные один раз, дальше всё делают сами. Нашу мы назвали «Normal()».
Не важно, сколько у полигона вершин, хоть 20 — всё, что нужно, это 3 его точки (первые 3). Этого достаточно, чтобы найти плоскость. Потом находим от неё нормаль.
Исходные коды к уроку
Векторный анализ — раздел математики, распространяющий методы математического анализа на векторы, как правило в двух- или трёхмерном пространстве. Объектами приложения векторного анализа являются: Векторные поля — отображения одного векторного пространства в другое.
Скалярное поле. Поверхности и линии уровня. Производная по направлению
Если в каждой точке пространства или части пространства определено значение некоторой величины, то говорят, что задано иоде данной величины. Поле называется скалярным, если рассматриваемая величина скалярна, т.е. вполне характеризуется своим числовым значением. Например, поле температур.
Скалярное поле задается скалярной функцией точки и = f(М). Если в пространстве введена декартова система координат, то и есть функция трех переменных х, у, z — координат точки М:
u = f(x,y,z). (1)
Определение:
Поверхностью уровня скалярного поля называется множество точек, в которых функция f(М) принимает одно и то же значение. Уравнение поверхности уровня
f(x, y, z) = с = const. (2)
Пример:
Найти поверхности уровня скалярного поля
Согласно определению уравнением поверхности уровня будет
Это уравнение сферы (с ≠ 0) с центром в начале координат.
Скалярное поле называется плоским, если во всех плоскостях, параллельных некоторой плоскости, поле одно и то же. Если указанную плоскость принять за плоскость хОу, го функция поля не будет зависеть от координаты г, т. е. будетфункцией только аргументов х и у,
u=f(x, y). (3)
Плоское поле можно характеризовать с помощью линий уровня — множества точек плоскости, в которых функция f(x, у) имеет одно и то же значение. Уравнение линии уровня —
f(х, у) = с = const. (4)
Пример:
Найти линии уровня скалярного поля
Линии уровня задаются уравнениями
При с = О получаем пару прямых у = х, у = -х. При с ≠ 0 получаем семейство гипербол (рис. 1).
Производная по направлению
Пусть имеется скалярное поле, определяемое скалярной функцией и = f(M). Возьмем точку М0 и выберем направление, определяемое вектором I. Возьмем другую точку М так, чтобы вектор М0М был параллелен вектору 1 (рис.2). Обозначим длину вектора МоМ через ∆l, а приращение функции f(М) — f(Mo), соответствующее перемещению ∆l, через ∆и. Отношение
определяет среднюю скорость изменения скалярного поля на единицу длины поданному направлению I.
Пусть теперь ∆l стремится к нулю так, чтобы вектор М0М все время оставался параллельным вектору I.
Определение:
Если при ∆l —> 0 существует конечный предел отношения (5), то его называют производной функции и = f(M) в данной точке М0 по данному направлению I и обозначают символом
Так что, по определению,
(6)
Это определение не связано с выбором системы координат, т. е. Hocит вариантный характер.
Найдем выражение для производной по направлению в декартовой системе координат. Пусть функция f(М) = f(х, у, z) дифференцируема в точке Мо(хо, yо, zо). Рассмотрим значение f(M) в точке М(х0 + ∆х,у0 + ∆y, zo + ∆z). Тогда полное приращение функции можно записать в следующем виде:
где ε = 0 при
а символы
означают, что частные производные вычислены в точке Мо. Отсюда (7)
Здесь величины суть направляющие косинусы вектора МоМ = ∆xi + ∆yj + ∆zk. Так как векторы МоМ и I сонаправлены (М0М ↑↑ I), то их направляющие косинусы одинаковы:
где
Так как M —» Mo, оставаясь все время на прямой, параллельной вектору I, то углы а, β, γ постоянны, а потому
Окончательно из равенств (7) и (8) получаем
Замечание:
Частные производные являются производными функции и по направлениям координатных осей Ox, Оу, Oz соответственно.
Пример:
Найти производную функции
в точке Mo(3,0,2) по направлению к точке M1(4,1, 3).
Имеем
Вектор МoМ = {1, 1, 1} имеет длину |МоМ| = /3. Его направляющие косинусы:
По формуле (9) будем иметь
Тот факт, что >0, означает, что скалярное поле в точке М0 в данном направлении возрастает.
Для плоского поля U = f(x, у) производная по направлению 1 в точке Мо(х0, у0) вычисляется по формуле (10)
где а — угол, образованный вектором I с осью Ох.
Замечание:
Формула (9) для вычисления производной по направлению I в данной точке М0 остается в силе и тогда, когда точка М стремится к точке Мо по кривой, для которой вектор I является касательным в точке Мо.
Пример:
Вычислить производную скалярного поля
и = arctg(xy)
в точке Mo(1, 1), принадлежащей параболе у = х2, по направлению этой кривой (в направлении возрастания абсциссы).
Направлением I параболы у = х2 в точке Mo{1, 1) считается направление касательной к параболе в этой точке (рис.3).
Пусть касательная к параболе в точке Мо образует с осью Ох угол a. Тогда tga == 2, откуда направляющие косинусы касательной
Вычислим значения в точке Mo(1, 1). Имеем
Теперь по формуле (10) получаем
Пример:
Найти производную скалярного поля и = In(xy + yz + zx) в точке Mo(0, 1, 1) по направлению окружности
Векторное уравнение окружности имеет вид
Находим единичный вектор т касательной к окружности
Точке Mo(0,1, 1) соответствует значение параметра t= π/2 Значение т в точке Мо будет равно
Отсюда получаем направляющие косинусы касательной к окружности в точке Mо: cos a = — 1, cos β = 0, cos γ = 0.
Вычислим значения частных производных данного скалярного поля в точке Mo(0, 1, 1)
Значит, искомая производная
Градиент скалярного поля
Пусть скалярное поле определяется скалярной функцией
u = f(x, y. z),
которая предполагается дифференцируемой.
Определение:
Градиентом скалярного поля u в данной точке М называется вектор, обозначаемый символом grad и и определяемый равенством
(1)
Ясно, что этот вектор зависит от функции f, так и от точки М, в которой вычисляется ее производная.
Пусть I° — единичный вектор в направлении I, т. е.
Тогда формулу для производной по направлению можно записать в следующем виде:
(3)
тем самым производная от функции и по направлению I равна скалярному произведению градиента функции u(M) на орт I° направления I.
Основные свойства градиента
Теорема:
Градиент скалярного поля перпендикулярен к поверхности уровня (или к линии уровня, если поле плоское).
Проведем через произвольную точку М поверхность уровня и = const и выберем на этой поверхности гладкую кривую L, проходящую через точку М (рис. 4). Пусть 1 — вектор, касательный к кривой L в точке М.
Так как на поверхности уровня и(М) = и(М1) для любой точки М1 ∈ L, то
С другой стороны, = (grad и, I°). Поэтому (grad и, I°) = 0. Это означает, что векторы grad и и I° ортогональны, grad u ⊥ I°.
Итак, вектор grad и ортогонален к любой касательной к поверхности уровня в точке М. Тем самым он ортогонален к самой поверхности уровня в точке М.
Теорема:
Градиент направлен в сторону возрастания функции поля.
Ранее мы доказали, что градиент скалярного поля направлен по нормали к поверхности уровня, которая может быть ориентирована либо в сторону возрастания функции и(М), либо в сторону ее убывания.
Обозначим через п нормаль к поверхности уровня, ориентированную в сторону возрастания функции и(М), и найдем производную функции и в направлении этой нормали (рис. 5). Имеем
Так как по условию и(М1) > и(М), то и(М1) — и(М) > 0, и поэтому
т.е. (grad и, п°) ≥ 0.
Отсюда следует, что grad и направлен в ту же сторону, что и выбранная нами нормаль п, т.е. в сторону возрастания функции и(М).
Теорема:
Длина градиента равна наибольшей производной по направлению в данной точке поля,
(здесь mах берется по всевозможным направлениям в данной точке М поля).
Имеем
где φ — угол между векторами I и grad n. Так как наибольшее значение cos φ равно 1, то наибольшим значением производной как раз и является |grad и|.
Пример:
Найти направление наибольшего изменения скалярного поля
и = ху + yz + zx
в точке Mо(1, 1, 1), а также величину этого наибольшего изменения в указанной точке.
Направление наибольшего изменения скалярного поля указывается вектором grad u(M). Имеем grad u(М) = (у + z)i + (х + г)j + (у + х)к, так что
grad и{Мo) = 2i + 2j + 2k.
Этот вектор определяет направление наибольшего возрастания поля в точке Мо(1,1,1). Величина наибольшего изменения поля в этой точке равна
Инвариантное определение градиента
Величины, характеризующие свойства изучаемого объекта и не зависящие от выбора системы координат, называются инвариантами данного объекта. Например, длина кривой — инвариант этой кривой, а угол касательной к кривой с осью Ох — не инвариант.
Основываясь на доказанных выше трех свойствах градиента скалярного поля, можно дать следующее инвариантное определение градиента.
Определение:
Градиент скалярного поля есть вектор, направленный по нормали к поверхности уровня в сторону возрастания функции поля и имеющий длину, равную наибольшей производной по направлению (в данной точке).
Пусть п° — единичный вектор нормали, направленный в сторону возрастания поля. Тогда
Пример:
Найти градиент расстояния
где Мo(хo,уo,zo) — некоторая фиксированная точка, а М(х,у,z) — текущая.
Имеем
где r° — единичный вектор направления MoM.
Правила вычисления градиента
- grad си(М) = с grad и{М), где с — постоянное число.
- grad(u + v) = grad и + grad v.
Приведенные формулы получаются непосредственно из определения градиента и свойств производных.
3. grad(u v) = v grad и+ и grad v.
По правилу дифференцирования произведения
Доказательство аналогично доказательству свойства 3.
Пусть F(u) — дифференцируемая скалярная функция. Тогда
grad F(u) = F'(u) grad и.
По определению градиента имеем
Применим ко всем слагаемым правой части правило дифференцирования сложной функции. Получим
В частности,
grad F(r) = F'(r) • p°. (6)
Формула (6) следует из формулы grad r = r°.
Пример:
Найти производную по направлению радиус-вектора r от функции u = sin r, где r = |r|. По формуле (3)
а по формуле (6) grad sin r = cos r • r° . В результате получим, что
Пример:
Пусть дано плоское скалярное поле
где r1, r2 — расстояния от некоторой точки Р(х,у) плоскости до двух фиксированных точек F1 и F2 этой плоскости, F1 ≠ F2.
Рассмотрим произвольный эллипс с фокусами F1 и F2 и докажем, что всякий луч света, вышедший из одного фокуса эллипса, после отражения от эллипса попадает в другой его фокус.
Линии уровня функции (7) суть
Уравнения (8) описывают семейство эллипсов с фокусами в точках F1 и F2.
Согласно результату примера 2 имеем
Тем самым градиент заданного поля равен вектору PQ диагонали ромба, построенного на ортах радиус-векторов, проведенных к точке Р(х,у) из фокусов F1 и F2, и значит, лежит на биссектрисе угла между этими радиус-векторами (рис. 6).
По теореме 1 градиент PQ перпендикулярен к эллипсу (8) в точке Р(х,у). Следовательно, нормаль к эллипсу (8) в любой его точке делит пополам угол между радиус-векторами, проведенными в эту точку. Отсюда и из того, что угол падения равен углу отражения, получаем: луч света, вышедший из одного фокуса эллипса, отразившись от него, непременно попадает в другой фокус этого эллипса.
Векторное поле. Векторные линии и их дифференциальные уравнения
Определение:
Если в каждой точке M(x,y,z) пространства или части пространства определена векторная величина
то говорят, что там задано векторное поле а.
Задание векторного поля равносильно заданию трех скалярных функций от трех переменных Р(х, у, z), Q(x, у, z), R(x, у, z).
Примерами векторных полей могут служить: силовое поле — поле некоторой силы F, поле скоростей v течения некоторой жидкости и др.
Для геометрической характеристики векторного поля служат векторные линии. Векторной линией векторного поля а называется кривая, касательная к которой в любой точке М имеет то же направление, что и вектор поля а в этой точке (рис. 7).
В силовом поле векторные линии называются силовыми линиями‘, в поле скоростей движения жидкости векторные линии называются линиями тока.
Дифференциальные уравнения векторных линий
Пусть векторное поле определяется вектор-функцией
где P(x,y,z), Q(x, у, z), R(x,y,z) — непрерывные функции переменных х, у, z, имеющие ограниченные частные производные первого порядка. Пусть
r(t) = x(t)i + y(t)j + z(t)k
— есть радиус-вектор текущей точки векторной линии векторного поля a (t — параметр). Из определения векторной линии следует, что вектор
а = Р(х, у, z)i + Q{x, у, z)j + R(x, у, z)k
и вектор касательной к этой кривой
должны быть коллинеарны в каждой точке векторной линии. Условием коллинеарности векторов является пропорциональность их координат:
Таким образом, мы получили для векторных линий систему дифференциальных уравнений в симметричной форме.
Допустим, что нам удалось найти два независимых интеграла системы (2): (3)
Система уравнений (3) определяет векторную линию как линию пересечения двух поверхностей. Произвольно меняя параметры c1 и c2 мы получаем семейство векторных линий как семейство с двумя степенями свободы.
Пример:
Hайти векторные линии векторного поля
а = хi + уj + 2zk.
Выписываем дифференциальные уравнения векторных линий, dx dy dz
или
Интегрируя эту систему, получим два уравнения
где c1, c2 — произвольные постоянные. Пересечение плоскостей у = c1х с параболическими цилиндрами z = c2x2 дает двухпараметрическое семейство векторных линий поля (рис. 8).
Определение:
Векторное поле называется плоским, если все векторы а параллельны одной и той же плоскости и в каждой плоскости, параллельной указанной, векторное поле одно и то же.
Посмотрим, как плоское векторное поле описывается в координатах. Если указанную в определении плоскость (или любую ей параллельную) принять за плоскость хОу, то векторы плоского поля не будут содержать компоненты по оси Oz и координаты векторов не будут зависеть от z:
Дифференциальные уравнения векторныхлиний плоского поля можно записать в следующем виде
или
(5)
Отсюда видно, что векторные линии плоского поля являются плоскими кривыми, лежащими в плоскостях, параллельных плоскости хОу.
Пример:
Найти векторные линии магнитного поля бесконечно длинного прямого провода.
Предположим, что проводник направлен вдоль оси Oz и по нему течет ток силы J, т. е. вектор тока
J = J • k.
Тогда вектор напряженности Н магнитного поля определяется по формуле
где
р = xi + yj + zk
— радиус-вектор точки М, р — расстояние от оси провода до точки М. Раскрывая векторное произведение (6), получим
Дифференциальные уравнения векторных линий:
Отсюда x = const, или xdx + ydy = 0. Окончательно имеем
т.е. векторные линии являются окружностями с центрами на оси Oz (рис.9).
Пример:
Найти векторные линии поля сил тяготения, образованного притягивающей материальной точкой массы т, расположенной в начале координат.
В данном случае сила F определяется так:
Дифференциальные уравнения векторных линий:
откуда, умножая каждую из дробей на получим
Чтобы получить уравнения векторных линий в параметрической форме, приравняем каждую из дробей величине . Имеем
Это — полупрямые, выходящие из начала координат.
Чтобы из семейства векторных линий выделить одну, надо задать точку М0(хо, yo, zо). через которую эта векторная линия должна проходить, и по координатам заданной точки определить величины С1, C2, C3.
Пусть, например, точка Мо имеет координаты хо = 3, yо = 5, zо = 7. Уравнение векторной линии, проходящей через точку Mo(3, 5, 7), можно записать так:
x = 3t, у — 5t, z = 7t.
Сама точка Мо получается при значении параметра t = 1.
Поток вектора через поверхность и его свойства
Рассмотрим сначала частный случай поля скоростей v течения жидкости. Выделим в поле некоторую поверхность Σ. Потоком жидкости через поверхность Σ называется количество жидкости, протекающее через поверхность Σ за единицу времени.
Этот поток легко вычислить, если скорость течения постоянна (v = const), а поверхность Σ —плоская. В этом случае поток жидкости равен объему цилиндрического тела с параллельными основаниями и образующими длины |v|, так как за единицу времени каждая частица перемещается на величину v (рис. 10),
П =Sh,
где S — площадь основания, h = npnv = (v, n°) — высота цилиндра и n — нормаль к его основанию, |п°| = 1.
Итак, при постоянной скорости v поток жидкости через плоскую поверхность Σ равен
(1)
Если скорость v изменяется непрерывно, а поверхность Σ — гладкая, то можно разбить поверхность Σ на столь малые части Σk (k = 1, 2,…, п), чтобы каждую часть Σk можно было приближенно считать плоской и вектор v на ней постоянным.
Так как поток жидкости через поверхность Σ равен сумме потоков жидкости через все ее части Σk, то мы получаем для вычисления потока приближенную формулу (2)
где п — общее число частей Σk, на которые разбита поверхность Σ, Рк — точка, лежащая на k-ой части, ∆σk — площадь части Σk поверхности, ( v, n°)pk означает скалярное произведение векторов v и п° в точке Pk ∈ Σk (рис. 11).
Назовем потоком жидкости через поверхность Σ предел суммы (2) при стремлении к нулю наибольшего из диаметров площадок Σk,
где d — наибольший из диаметров частей Σk (k= 1,2,…,п). Интеграл (3), определяющий поток жидкости, берется от скалярной функции (v, п°) по площади поверхности Σ.
Понятие потока произвольного вектора а через поверхность Σ вводится п о аналогии с введенным выше понятием потока жидкости через поверхность.
Определение:
Потоком вектора (векторного поля) а через поверхность Σ называется интеграл по поверхности Σ от проекции вектора а на нормаль к поверхности
Ясно, что интеграл (4) существует, если вектор а = Pi+Qj+Rk непрерывен, т. е. непрерывны его координаты Р(x, у, z), Q(x, у, z), R(x, y,z), и поверхность Σ — гладкая, т. е. имеет непрерывно меняющуюся касательную плоскость
Пример:
Поле создается точечным зарядом (электричесkое поле) или точечной маcсой (поле тяготения), помещенными в начале координат. Тогда вектор напряженности поля в любой точке Р будет равен
где q — величина заряда (массы), r = ОР — радиус-вектор точки Р. Требуется найти поток вектора напряженности Е через SR — сферу радиуса R с центром в начале координат.
Так как направление нормали к сфере совпадает с направлением радиус-вектора r, то п° = r° и поэтому
На сфере SR радиуса R имеем r = R, так что (Е, n°) = = const. Поэтому поток вектора через SR равен
Свойства потока вектора через поверхность
1. Линейность.
(5)
где λ и μ — постоянные числа.
2. Аддитивность. Если поверхность Σ разбита кусочно-гладкой кривой на две части Σ1 и Σ2, то поток через поверхность Σ равен сумме потоков через поверхности Σ1 и Σ2,
(6)
Это свойство позволяет распространить понятие потока на кусочно-гладкие поверхности Σ.
Понятие ориентации поверхности
Взяв, к примеру, цилиндрическую поверхность, замечаем, что если в некоторой ее точке М выбрать определенный (один из двух) единичный вектор нормали и непрерывно перемещаться затем по поверхности вместе с соответствующим вектором нормали по любому пути, не переходящему через край повержюсти, то при возвращении в точку М единичный вектор нормали совпадает с исходным (рис. 12).
Вместе с тем, существуют поверхности, для которых это не так. Примером такой поверхности может служить лист Мёбиуса (рис. 13). Существует путь (отмеченная на рисунке пунктиром средняя линия листа), перемещаясь по которому, мы возвратимся в начальную точку с единичным вектором нормали, противоположным исходному.
Описанное свойство разбивает все поверхности на два класса — двусторонние, или ориентируемые (плоскость, сфера, поверхность куба и т.п.), и односторонние, или неориентируемые (лист Мёбиуса).
3. Зависимость потока от ориентации поверхности (от ориентации вектора нормали к поверхности). Понятие потока вводится только для двусторонних поверхностей. Будем считать, что если в одной точке такой поверхности направление вектора нормали уже выбрано, то в любой другой ее точке берется тот вектор нормали, который получается из выбранного при непрерывном перемещении точки по поверхности (без перехода через границу). В частности, на замкнутой поверхности во всех точках берется либо внешняя нормаль, либо внутренняя (внутренняя нормаль направлена внутрь тела, ограниченного замкнутой поверхностью).
Обозначим через Σ+ ту сторону поверхности Σ, на которой выбран вектор нормали п+ = п, а через Σ- — сторону поверхности Σ, на которой берется вектор нормали (п_ = -п). Тогда получим
(7)
где п°_ = -п°+. Таким образом, при изменении ориентации поверхности (при изменении направления вектора нормали п° к поверхности Σ) поток вектора меняет знак на противоположный.
Пример:
Вычислить поток радиус-вектора r = хi + yj + zk через поверхность прямого кругового цилиндра высоты Н с радиусом основания R и осью Oz.
Поверхность Σ состоит из трех частей: боковой поверхности Σ1, верхнего основания Σ2 и нижнего основания Σ3 цилиндра. Искомый поток П в силу свойства аддитивности равен
П = П1 +П2 + П3,
где П1, П2, П3 — потоки данного поля через Σ1, Σ2 и Σ3 соответственно.
На боковой поверхности цилиндра вектор внешней нормали п°1 параллелен плоскости хОу, и поэтому
(см. рис. 14). Следовательно,
На верхнем основании Σ2 вектор нормали п°2 параллелен оси Оz, и поэтому можно положить п°2 = k. Тогда имеем
так что
На нижнем основании Σ3 вектор г перпендикулярен к вектору нормали п°3 = -k. Поэтому (r, п°3) = (r, -k) = 0 и
Значит, искомый поток
Здесь символ означает двойной интеграл по замкнутой поверхности.
Поток вектора через незамкнутую поверхность
Укажем некоторые способы вычисления потока вектора через незамкнутые поверхности.
Метод проектирования на одну из координатных плоскостей
Пусть поверхность S однозначно проектируется на область Dxy плоскости хОу. В этом случае поверхность S можно задать уравнением вида
z = f(x, у).
Орт п° нормали к поверхности S находится по формуле
Если в формуле (1) берется знак « -», то угол γ между осью Oz и нормалью п° —острый; если же знак «+», то угол γ — тупой.
Так как элемент площади dσ этой поверхности равен
то вычисление потока П через выбранную сторону поверхности S сводится к вычислению двойного интеграла по формуле
(3)
Символ
означает, что при вычислении в подынтегральной функции надо вместо z всюду поставить f(i, у).
Пример:
Найти поток вектора
через часть поверхности параболоида
отсеченной плоскостью z = 2. По отношению к области, ограниченной параболоидом, берется внешняя нормаль (рис. 15).
Данная поверхность проектируется на круг Dxy плоскости хОу с центром в начале координат радиуса R =. Находим орт п° нормали к параболоиду:
Согласно условию задачи вектор п° образует с осью Oz тупой угол γ, поэтому перед дробью следует взять знак минус. Таким образом,
Находим скалярное произведение
Если поверхность S проектируется однозначно на область Dyz плоскости yOz, то ее можно задать уравнением х = φ{у, z). В этом случае имеем
(6)
где
Знак «+» в последней формуле соответствует тому, что угол а между осью Ох и вектором нормали п° острый, и знак «-», если указанный угол тупой.
Наконец, если поверхность S проектируется однозначно на область Dxz плоскости xOz, то ее можно задать уравнением у = ψ(x, z) и тогда
Знак «+» перед дробью в формуле (10) означает, что угол β между осью Оу и вектором нормали п° — острый, а знак «—», что угол β — тупой.
Замечание:
Для нахождения потока вектора
а = Р(х, у, z)i + Q(z, у, z)j + R(х, у, k)
к через поверхность S, заданную уравнением z = f(x, у), методом проектирования на координатную плоскость хОу, не обязательно находить орт п° нормали, а можно брать вектор
Тогда формула (2) для вычисления потока П примет вид:
(11)
Аналогичные формулы получаются для потоков через поверхности, заданные уравнениями х = φ(у, z) или у = ψ(х, z).
Пример:
Вычислить поток вектора
а = хzi
через внешнюю сторону параболоида
ограниченного плоскостью z = 0 (рис. 16).
Имеем
n = ±(2ri + 2yj+k).
Так как угол γ — острый, следует выбрать знак «+». Отсюда
Искомый поток вычисляется так:
Переходя к полярным координатам х = р cos φ, y = p sin φ, 0 ≤ р ≤ 1. 0 ≤ φ < 2π, получим
Метод проектирования на все координатные плоскости
Пусть поверхность S однозначно проектируется на все три координатные плоскости. Обозначим через Dxy, Dxz, Dyz проекции S на плоскости хОу, xOz, yOz соответственно. В этом случае уравнение F(x, y, z) = 0 поверхности S однозначно разрешимо относительно каждого из аргументов, т. е.
x = x(y,z), y = y(x,z), z = z(x,y). (12)
Тогда поток вектора
а = Р(х, у, z)i + Q(x, у, z)j + R(x, у, z)k
через поверхность S, единичный вектор нормали к которой равен
можно записать так:
Известно, что
причем знак в каждой из формул (14) выбирается таким, каков знак cos a, cos β, cos γ на поверхности S. Подставляя соотношения (12) и (14) в формулу (13), получаем, что (15)
Пример:
Вычислить поток векторного поля
а = yi + zj + zk
через треугольник, ограниченный плоскостями z + y+ z = l (l>0), x=0, у — 0, z = 0 (угол γ — острый) (рис. 17).
Имеем
так что
Значит, перед всеми интегралами в формуле (15) следует взять знак « + ». Полагая Р = у, Q = z, R = х, получим
Вычислим первый интеграл в правой части формулы (16). Область Dyz — треугольник ВОС в плоскости yOz, уравнение стороны ВС: y+z = l, 0 ≤ у ≤ I. Имеем
Аналогично получим
Значит, искомый лоток равен
Метод введения криволинейных координат на поверхности
Если поверхность S является частью кругового цилиндра или сферы, при вычислении потока удобно, не применяя проектирования на координатные плоскости, ввести на поверхности криволинейные координаты.
А. Поверхность S является частью кругового цилиндра
ограниченного поверхностями z = f1(x,y) и z = f2(х. у), где f1(x. y) ≤ f2(x, y) (рис. 18). Полагая х = R cos φ, у = R sin φ, z = z, будем иметь
Элемент площади поверхности выражается так:
и поток вектора а через внешнюю сторону поверхности S вычисляется по формуле:
Пример:
Найти поток вектора
через внешнюю сторону поверхности цилиндра
ограниченной плоскостями
Имеем
Так как
то скалярное произведение (а, п°) на цилиндре (х = 2 cos φ, у = 2 sin φ, z = z) равно:
Тогда по формуле (18) получим
В. Поверхность S является частью сферы
ограниченной коническими поверхностями, уравнения которых в сферических координатах имеют вид и полуплоскостями
(рис. 19).Точки данной сферы описываются соотношениями
где Поэтому элемент площади
В этом случае поток векторного поля а через внешнюю часть поверхности S вычисляется по формуле
Пример:
Найти поток вектора
через внешнюю часть сферы
отсеченную плоcкостью z = 2 (рис. 20).
В данном случае имеем
Положим
Тогда скалярное произведение (а, п°) выразится так:
По формуле (21) получим
Замечание:
Здесь мы воспользовались формулой
Поток вектора через замкнутую поверхность. Теорема Гаусса—Остроградского
Теорема:
Если в некоторой области G пространства R3 координаты вектора
а = Р(х, у, z)i + Q(x, у, z)j + R(x, у, z)k
непрерывны и имеют непрерывные частные производные , то поток вектора а через любую замкнутую кусочно-гладкую поверхность S, лежащую в области G, равен тройному интегралу от
по области V, ограниченной поверхностью S:
Здесь п0 — орт внешней нормали к поверхности, а символ означает поток через замкнутую поверхность S. Эта формула называется формулой Гаусса—Остроградского.
Рассмотрим сначала вектор а, имеющий только одну компоненту а = R(x, у, z)k, и предположим, что гладкая поверхность S пересекается каждой прямой, параллельной оси Oz, не более чем в двух точках. Тогда поверхность S разбивается на две части S1 и S2, однозначно проектирующиеся на некоторую область D плоскости хОу (рис.21).
Внешняя нормаль к поверхности S2 образует острый угол γ с осью Oz, а внешняя нормаль к поверхности S1 образует тупой угол с осью Oz. Поэтому cos γ = (п°, к) > 0 на S2 и cos γ < 0 на S1, так что на S2 имеем cos γ = |cos γ|, а на S1 cos γ = —| cos γ|. В силу аддитивности потока имеем
Пусть dσ — элемент площади на поверхности S. Тогда
где dS — элемент площади области D. Сведем интегралы по поверхности к двойным интеграл ам по области D плоскости хОу, на которую проектируются поверхности S1 и S2. Пусть S2 описывается уравнением z = z2(x, у), а S, — уравнением z = z1(x, у). Тогда
Так как приращение непрерывно дифференцируемой функции можно представить как интеграл от ее производной
то для функции R(x, у, z) будем иметь
Пользуясь этим, получаем из формулы (3)
Если поверхность S содержит часть цилиндрической поверхности с образующими, параллельными оси Oz (рис. 22), то на этой части поверхности (Rk, п°) = 0 и интеграл ∫∫ (Rk, n°) dσ по ней равен нулю. Поэтому формула (4) остается справедливой и для поверхностей, содержащих указанные цилиндрические части.
Формула (4) переносится и на случай, когда поверхность 5 пересекается вертикальной прямой более, чем в двух точках (рис. 23).
Разрежем область V на части, поверхность каждой из которых пересекается вертикальной прямой не более чем в двух точках, и обозначим через Sp поверхность разреза. Пусть S1 и S2 — те части поверхности S, на которые она разбивается разрезом Sp,a V1 и V2 — соответствующие части области V, ограниченные поверхностями
Здесь Sp+ означает, что вектор нормали к разрезу Sp направлен вверх (образует с осью Oz острый угол), a Sp— — что этот вектор нормали направлен вниз (образует с осью Oz тупой угол). Имеем:
Складывая полученные равенства и пользуясь аддитивностью потока и тройного интеграла, получим
(интегралы по разрезу Sp взаимно уничтожаются). Рассмотрим, наконец, вектор
а = Р{х, у, z)i + Q(x, у, z)j + R(x, у, z)k.
Для каждой компоненты Pi, Qj, Rк мы можем написать формулу, аналогичную формуле (4) (все компоненты равноправны). Получим
Складывая эти равенства и пользуясь линейностью потока и тройного интеграла, получаем формулу Гаусса—Остроградского
Пример:
Вычислить поток вектора
а = 2xi — (z — 1)k
через замкнутую поверхность
1) по определению, 2) по формуле Остроградского.
1) Поток вектора а равен сумме
(рис. 24), где
Перейдем на цилиндре к криволинейным координатам
Тогда
Следовательно, П = -4π + 0 + 8π = 4π.
2) По формуле Гаусса—Остроградского имеем
Пример:
Вычислить поток радиус-вектора
r = xi + yj + zk
через сферу радиуса R с центром в начале координат:
1) по определению; 2) по формуле Остроградского.
1) Так как для сферы
2) Сначала находим
Отсюда
Пример:
Вычислить поток вектора
через замкнутую поверхность S, заданную условиями:
1) по определению; 2) по формуле Остроградстого (рис.25).
1} Имеем
где
(на S1 имеем z = 0),
Поэтому
Имеем
Поэтому
Переходя к цилиндрическим координатам
и замечая, что z = 9 — р на поверхности S, имеем
Замечание:
При вычислении потока через незамкнутую поверхность часто бывает удобно подходящим образом дополнить ее до замкнутой и воспользоваться формулой Гаусса—Оcтроградского.
Пример:
Вычислить поток вектора
через поверхность S:
Заданная поверхность S есть конус с осью Оу (рис. 26).
Замкнем этот конус куском Σ плоскости у = I. Тогда, обозначая через П1 искомый поток, а через П2 поток по поверхности Σ, будем иметь
где V — объем конуса, ограниченного поверхностями S и Σ.
Так как
т.к. на поверхности Σ выполняется равенство у = 1. Следовательно, П1 = π.
Дивергенция векторного поля. Соленоидальные (трубчатые) поля
Пусть S — замкнутая поверхность. Рассмотрим поле скоростей v течения жидкости и вычислим поток жидкости через поверхность 5. Если он положителен, то это означает, что из той части пространства, которая ограничена поверхностью S, вытекает больше жидкости, чем втекает в нее. В этом случае говорят, что внутри S имеются источники (выделяющие жидкость). Напротив, если поток отрицателен, то внутрь S втекает больше жидкости, чем вытекает из нее. В этом случае говорят, что внутри S имеются стоки (поглощающие жидкость).
Тем самым, величина
позволяет судить о природе части векторного поля, заключенного внутри поверхности S, а именно, о наличии источников или стоков внутри нее и их производительности (мощности).
Понятие о потоке вектора через замкнутую поверхность приводит к понятию дивергенции, или расходимости поля, которое дает некоторую количественную характеристику поля в каждой его точке.
Пусть М — изучаемая точка поля. Окружим ее поверхностью S произвольной формы, например, сферой достаточно малого радиуса. Область, ограниченную поверхностью 5, обозначим через (V), а ее объем через V.
Вычислим поток вектора а через поверхность S. Имеем
Составим отношение этого потока П к величине объема V,
Так как числитель представляет собой производительность источников (стоков) внутри области (V), то отношение (1) дает среднюю производительность единицы объема.
Определение:
Если отношение (1) имеет конечный предел, когда область (V) стягивается в точку М, то этот предел называют дивергенцией векторного поля (дивергенцией вектора а) в точке М и обозначают div а(М). То есть по определению
Дивергенция векторного поля есть скалярная величина (числитель и знаменатель дроби (2) суть скалярные величины).
Если diva(M) > 0, то в точке М расположен источник, если diva(M) < 0, то в точке М — сток.
Формула (2) позволяет сделать следующее заключение: дивергенция поля а в точке М есть объемная плотность потока вектора а в этой точке. Эта формула дает инвариантное определение дивергенции, не связанное с выбором систем координат — все величины, входящие в формулу (2), определяются непосредственно самим полем и от координатной системы не зависят.
Покажем, как вычисляется дивергенция в декартовых координатах при условии, что координаты вектора
а = Р(х, у, z)i + Q(x, у, z)j + R(x, y, z)k
непрерывны и имеют непрерывные частные производные в окрестности точки М. Тогда к потоку вектора а через любую замкнутую поверхность 5, расположенную в окрестности точки М, можно применить формулу Гаусса—Остроградского
Пользуясь теоремой о среднем для тройного интеграла, получим
Подставляя это выражение в формулу (2), определяющую дивергенцию, найдем
Когда область (V) стягивается в точку М, то и точка Мcp стремится к точке М и, в силу предположенной непрерывности частных производных, получаем
или, короче,
(3)
(все величины в формуле (3) вычисляются водной и той же точке).
Формула (3) дает выражение дивергенции в декартовых координатах. Попутно доказано само существование дивергенции вектора а при условии, что производные непрерывны.
Используя формулу (3) для дивергенции, запишем формулу Гаусса—Остроградского в векторной форме. Имеем
(4)
— поток вектора а через замкнутую поверхность S равен тройному интегралу от дивергенции вектора а по области (V), ограниченной поверхностью S.
Правила вычисления дивергенции
1, Дивергенция обладает свойством линейности
(5)
где С1,…, Сп — постоянные числа.
Пусть
а = Р(х, у, z)i + Q(x, у, z)j + R(x, у, z)k
и С — постоянное число. Тогда
Если
то
2. Дивергенция постоянного вектора с равна нулю
div e = 0. (6)
3. Дивергенция произведения скалярной функции и(М) на вектор а(М) вычисляется по формуле
div(ua) = u diva + (gprad u,a). (7)
В самом деле,
Пример:
Найти дивергенцию вектора
где r = |r| — расстояние от начала координат до переменной точки М(х,у,z),
По формуле (7) имеем
Так как r = xi + уj + zk. то
Трубчатое (соленоидальное) поле
Если во всех точках некоторой области G дивергенция векторного поля, заданного в этой области, равна нулю
div а ≡ 0, (8)
то говорят, что в этой области поле соленоидальное (или трубчатое).
Из формулы Гаусса—Остроградского вытекает, что в трубчатом поле поток вектора через любую замкнутую поверхность S, лежащую в этом поле, равен нулю
(9)
Свойства трубчатого поля
Рассмотрим в области, где задано поле вектора а, какую-нибудь площадку Σ (рис.27). Назовем векторной трубкой совокупность векторных линий, проходящих через границу γ = θΣ этой площадки. Пусть Σ1 — некоторое сечение векторной трубки. Выберем вектор нормали щ к сечению Σ1 так, чтобы он был направлен в ту же сторону, что и вектор а поля.
Теорема:
В трубчатом поле поток вектора а через любое сечение векторной трубки один и тот же.
Пусть Σ1 и Σ2 —непересекающиеся сечения одной и той же векторной трубки. Надо доказать, что
Обозначим через Σ3 часть поверхности векторной трубки, заключенную между сечениями Σ1 и Σ2. Поверхности Σ1, Σ2, Σ3 вместе образуют замкнутую поверхность Σ (рис.28).
Так как по условию поле вектора а — трубчатое, то
В силу аддитивности потока соотношение (10) можно переписать так:
В точках поверхности Σ3, составленной из векторных линий, имеем , так что (а, п°3) = 0 на Σз, и значит, последний интеграл в левой части (11) равен нулю. Таким образом, из (11) находим
Пусть поверхность Σ имеет ориентированный замкнутый контур L своей границей. Будем говорить, что поверхность Σ натянута на контур L. Вектор нормали п к поверхности Σ будем ориентировать так, чтобы из конца нормали обход контура L был виден против часовой стрелки (рис. 29).
Теорема:
В трубчатом поле поток вектора а через любую поверхность, натянутую на данный контур, один и тот же:
Замечание:
В трубчатом поле векторные линии могут быть либо замкнутыми кривыми, либо иметь концы на границе области, где поле задано.
Пример:
Рассмотрим силовое поле, создаваемое точечным зарядом q, помешенным в начале координат. Вычислим дивергенцию вектора Е напряженности
Имеем
где
Пользуясь формулой (7), получим
для r ≠ 0. Таким образом, поле вектора Σ, заданного формулой (13), будет трубчатым в любой области G, не содержащей точки O(0,0,0).
Вычислим поток вектора Σ через сферу Sr радиуса R с центром в начале координат O(0,0,0) (рис.30).
Имеем
Замечание:
Можно показать, что поток вектора (13) через любую замкнутую поверхность Σ, охватывающую точку O(0,0,0), всегда равен 4 πg.
Циркуляция векторного поля. Ротор вектора. Теорема Стокса
Пусть в некоторой области G задано непрерывное векторное поле
а(М) = Р(х, у, х)i + Q(x, у, z)j + R(х, у, z)k
и замкнутый ориентированный контур L.
Определение:
Циркуляцией вектора а по замкнутому контуру L называется криволинейный интеграл 2-го рода от вектора а по контуру L
Здесь dr — вектор, длина которого равна дифференциалу дуги L, а направление совпадаете направлением касательной к L, определяемым ориентацией контура (рис. 31) символ означает, что интеграл берется по замкнутому контуру L.
Пример:
Вычислить циркуляцию векторного поля
вдоль эллипса L:
По определению циркуляции имеем
Параметрические уравнения данного эллипса имеют вид:
и, значит, dx = -a sin tdt, dy = b cos tdt. Подставляя эти выражения в формулу (2), найдем
Ротор (вихрь) векторного поля
Рассмотрим поле вектора
а(М) = Р(х, у, z)i + Q(x, у, z)j + R(х, у, z)k,
Р, Q, R которого непрерывны и имеют непрерывные частные производные первого порядка по всем своим аргументам.
Определение:
Ротором вектора а(M) называется вектор, обозначаемый символом rot а и определяемый равенством
или, в символической, удобной для запоминания форме,
Этот определитель раскрывают по элементам первой строки, при этом операции умножения элементов второй строки на элементы третьей строки понимаются как операции дифференцирования, например,
Определение:
Если в некоторой области G имеем rot а = 0, то поле вектора а в области G называется безвихревым.
Пример:
Найти ротор вектора
Согласно формуле (3) имеем
Так как rot а — вектор, то мы можем рассматривать векторное поле — поле ротора вектора а. Предполагая, что координаты вектора а имеют непрерывные частные производные второго порядка, вычислим дивергенцию вектора rot а. Получим
т. е.
div rot a = 0. (3′)
Таким образом, поле вектора rot а соленоидально.
Теорема Стокса:
Циркуляция вектора а вдоль ориентированного замкнутого контура L равна потоку ротора этого вектора через любую поверхность Е, натянутую на контур L,
При этом предполагается, что координаты вектора а имеют непрерывные частные производные в некоторой области G пространства, содержащей поверхность Σ, и что ориентация орта нормали п° к поверхности Σ С G согласована с ориентацией контура L так, что из конца нормали обход контура в заданном направлении виден совершающимся против часовой стрелки.
Учитывая, что а = Pi + Qj + Rk, n° = cos ai + cos βj + cos γk, и пользуясь определением ротора (3), перепишем формулу (4) в следующем виде:
Рассмотрим сначала случай, когда гладкая поверхность Σ и ее контур L однозначно проектируются на область D плоскости хОу и ее границу — контур λ соответственно (рис. 32). Ориентация контура L порождает определенную ориентацию контура λ. Для определенности будем считать, что контур L ориентирован так, что поверхность Σ остается слева, так что веkтор нормали п к поверхности Σ составляет с осью Oz острый угол γ (cos γ > 0).
Пусть z = φ{х,у) — уравнение поверхности Σ и функция ф(х,у) непрерывна и имеет непрерывные частные производные в замкнутой области D. Рассмотрим интеграл
Линия L лежит на поверхности Σ. Поэтому, пользуясь уравнением этой поверхности z = φ(х, у),мы можем заменить z под знаком интеграла на φ(x, у). Координаты (х, у)
переменной точки кривой λ равны координатам соответствующей точки на кривой L, а потому интегрирование по L можно заменить интегрированием по λ,
Применим к интегралу, стоящему справа, формулу Грина. Имеем
Перейдем теперь от интеграла по области D к интегралу по поверхности Σ. Так как dS = cos γ • dσ,то из формулы (8) получим, что
Вектор нормали n° к поверхности Σ определяется выражением
или n° = cos a • i + cos β • j + cos γ • k. Отсюда видно, что
Поэтому равенство (9) можно переписать так:
Считая Σ гладкой поверхностью, однозначно проектирующейся на все три координатные плоскости, аналогично убеждаемся в справедливости формул
Складывая равенства (10), (11) и (12) почленно, получим формулу Стокса (5), или, короче,
Замечание:
Мы показали, что поле вектора rota — соленоидальное, и потому поток вектора rota не зависит от вида поверхности Σ, натянутой на контур L.
Замечание:
Формула (4) выведена в предположении, что поверхность Σ однозначно проектируется на все три координатные плоскости. Ecли это условие не выполнено, то разбиваем Σ на части так, чтобы каждая часть указан ному условию удовлетворяла, а затем пользуемся аддитивностью интегралов.
Пример:
Вычислить циркуляцию вектора
а = yi — xj + k
по линии L:
1) пользуясь определением; 2) по теореме Стокса.
1) Зададим линию L параметрически:
Тогда dx = -R sin t dt, dy = R cos t dt, H dz = 0, так что
2) Найдем rot a:
Натянем на контур L кусок плоскости z = H, так что п° = k. Тогда
Инвариантное определение ротора поля
Из теоремы Стокса можно получить инвариантное определение ротора поля, не связанное с выбором системы координат.
Теорема:
Проекция ротора а на любое направление не зависит от выбора системы координат и равна поверхностной плотности циркуляции вектора а по контуру площадки, перпендикулярной этому направлению,
Здесь ( Σ ) — плоская площадка, перпендикулярная вектору п; S — площадь этой площадки; L — контур площадки, ориентированный так, чтобы обход контура был виден из конца вектора п против хода часовой стрелки; ( Σ ) → М означает, что площадка ( Σ ) стягиваетcя к точке М, в которой рассматривается вектор rot а, причем вектор нормали п к этой площадке остается все время одним и тем же (рис. 33).
Применим сначала к циркуляции
вектора а теорему Стокса, а затем к полученному двойному интегралу — теорему о среднем значении:
(скалярное произведение (rot a, n°) берется в некоторой средней точке Mср площадки ( Σ )).
При стягивании площадки ( Σ ) к точке М средняя точка Мср тоже стремится к точке М и, в силу предполагаемой непрерывности частных производных от координат вектора а (а значит, и непрерывности rot а), мы получаем
Поскольку проекция вектора rot а на произвольное направление не зависит от выбора системы координат, то сам вектор rota инвариантен относительно этого выбора. Отсюда получаем следующее инвариантное определение ротора поля: ротор поля есть вектор, длина которого равна наибольшей поверхностной плотности циркуляции в данной точке, направленный перпендикулярно той площадке, на которой эта наибольшая плотность циркуляции достигается; при этом ориентация вектора rot a согласуется с ориентацией контура, при которой циркуляция положительна, по правилу правого винта.
Физический смысл ротора поля
Пусть твердое тело вращается вокруг неподвижной оси l с угловой скоростью w. Не нарушая общности, можно считать, что ось l совпадает с осью Oz (рис. 34). Пусть М(г) — изучаемая точка тела, где
r = xi + уj + zk.
Вектор угловой скорости в нашем случае равен w ≡ wk, вычислим вектор v линейной скорости точки М,
Отсюда
Итак, вихрь поля скоростей вращающегося твердого тела одинаков во всех точках поля, параллелен оси вращения и равен удвоенной угловой скорости вращения.
Правила вычисления ротора
1, Ротор постоянного вектора с равен нулевому вектору,
rot e = 0.
2. Ротор обладает свойством линейности
где c1, c2,…, cn — постоянные числа.
3. Ротор произведения скалярной функции и(М) на векторную а(М) вычисляется по формуле
rot(wa) = и rot а + [grad и, а].
В самом деле,
Независимость криволинейного интеграла от пути интегрирования
Определение:
Область G трехмерного пространства называется поверхностно односвязной, если на любой замкнутый контур, лежащий в этой области, можно натянуть поверхность, целиком лежащую в области G.
Например, внутренность сферы или все трехмерное пространство являются поверхностно односвязными областями; внутренность тора или трехмерное пространство, из которого исключена прямая, поверхностно односвязными областями не являются.
Пусть в поверхностно односвязной области G задано непрерывное векторное поле
а (М) = Р(М)i + Q(M)j + R(M) k.
Тогда имеет место следующая теорема.
Теорема:
Для того чтобы криволинейный интеграл
в поле вектора а не зависел от пути интегрирования, а зависел только от начальной и конечной точек пути (А и В), необходимо и достаточно, чтобы циркуляция вектора a вдаль любого замкнутого контура L, расположенного в области G, была равна нулю.
Необходимость. Пусть интеграл
не зависит от пути интегрирования. Покажем, что тогда
по любому замкнутому контуру L равен нулю.
Рассмотрим произвольный замкнутый контур L в поле вектора а и возьмем на нем произвольно точки A и В (рис.35).
По условию имеем
где L1 и L2 — различные пути, соединяющие точки А и В; откуда
Но L1 U L2 как раз и есть выбранный замкнутый контур L. Достаточность. Пусть
для любого замкнутого контура L. Покажем, что в этом случае интеграл
не зависит от пути интегрирования.
Возьмем в поле вектора а две точки А и В, соединим их произвольными линиями L1 и L2 к покажем, что
Для простоты ограничимся случаем, когда линии L1 и L2 не пересекаются. В этом случае объединение L1 ∪ L2 образует простой замкнутый контур L (рис. 36).
По условию
а по свойству аддитивности
Следовательно,
откуда справедливость равенства (2) и вытекает.
Теорема 9 выражает необходимое и достаточное условия независимости криволинейного интеграла от формы пути, однако эти условия трудно проверяемы. Приведем более эффективный критерий.
Теорема:
Для того, чтобы криволинейный интеграл
не зависел от пути интегрирования L, необходимо и достаточно, чтобы векторное поле а(М) = Р(X, у, z)i + Q(x, у, z)j + R(x, у, z)k было безвихревым,
rot a(M) = 0. (3)
Здесь предполагается, что координаты Р(х, у, z), Q(x, у, z), R(x, у, z) вектора а(М) имеют непрерывные частные производные первого порядка и область определения вектора а(М) поверхностно односвязна.
Замечание:
В силу теоремы 9 независимость криволинейного интеграла от пути интегрирования равносильна равенству нулю циркуляции вектора а вдоль любого замкнутого контура. Это обстоятельство мы используем при доказательстве теоремы.
Необходимость. Пусть криволинейный интеграл не зависит от формы пути, или, что то же, циркуляция вектора а по любому замкнутому контуру L равна нулю. Тогда
т. е. в каждой точке поля проекция вектора rot а на любое направление равна нулю. Это означает, что сам вектор rot а равен нулю во всех точках поля,
rot а ≡ 0.
Достаточность. Достаточность условия (3) вытекает из формулы Стокса, так как если rot а ≡ 0, то и циркуляция вектора по любому замкнутому контуру L равна нулю:
Ротор плоского поля a = P(x, y)i + Q(x, y)j равен
что позволяет сформулировать для плоского поля следующую теорему.
Теорема:
Для того, чтобы криволинейный интеграл
в односвязном плоском поле не зависел от формы линии L, необходимо и достаточно, чтобы соотношение
выполнялось тождественно во всей рассматриваемой области.
Если область неодносвязна, то выполнение условия
вообще говоря, не обеспечивает независимости криволинейного интеграла от формы линии.
Пример:
Пусть
Рассмотрим интеграл
Ясно, что подынтегральное выражение не имеет смысла в точке 0(0,0). Поэтому исключим эту точку. В остальной части плоскости (это будет уже не сщносвязная область!) координаты вектора а непрерывны, имеют непрерывные частные производные и
Рассмотрим интеграл (6) вдоль замкнутой кривой L — окружности радиуса R с центром в начале координат
Отличие циркуляции от нуля показывает, что интеграл (6) зависит от формы пути интегрирования.
Потенциальное поле
Определение:
Поле вектора а(М) называется потенциальным, если существует скалярная функция и{М) такая, что
grad и = a. (1)
При этом функция и{М) называется потенциалом поля ее поверхности уровня называются эквипотенциальными поверхностями.
Пусть
а = Р(х, у, z)i + Q(x, у, z)j + R(x, у, z)k.
Так как
то соотношение (1) равносильно следующим трем скалярным равенствам:
Заметим, что потенциал поля определяется с точностью до постоянного слагаемого: если grad и = а и grad v = а, то
и, следовательно, и = v + с, где с — постоянное число.
Пример:
Поле радиус-вектора г является потенциальным, так как
(напомним, что ). Потенциалом поля радиус-вектора является, следовательно,
Пример:
Поле вектора
является потенциальным.
Пусть функция φ(r) такая, что
найдена. Тогда
и
откуда
Значит,
— потенциал поля.
Теорема:
Для того чтобы поле вектора а было потенциальным, необходимо и достаточно, чтобы оно было безвихревым,
rot а ≡ 0, (2)
т. е. чтобы его ротор равнялся нулю во всех точках поля. При этом предполагается непрерывность всех частных производных от координат вектора а и поверхностная односвязность области, в которой задан вектор а.
Необходимость. Необходимость условия (2) устанавливается непосредственным подсчетом: если поле потенциально, т. е. а = grad и, то
в силу независимости смешанных производных от порядка дифференцирования.
Достаточность. Пусть поле вектора безвихревое (2). Для того чтобы доказать потенциальность этого поля, построим его потенциал и(М). Из условия (2) следует, что криволинейный интеграл
не зависит от формы линии L, а зависит только от ее начальной и конечной точек. Зафиксируем начальную точку Мо(xo, yо, zo), а конечную точку М(х, y, z) будем менять. Тогда интеграл (3) будет функцией точки М(х, у, z). Обозначим эту функцию через и(М) и докажем, что
grad u = а.
В дальнейшем будем записывать интеграл (3), указывая лишь начальную и конечную точку пути интегрирования,
Равенство grad и = а равносильно трем скалярны м равенства м
Докажем первое из них,
второе и третье равенства доказываются аналогично.
По определению частной производной имеем
Рассмотрим точку М1(х + ∆х, у, z), близкую к точке M(x,y,z). Так как функция и(М) определяется соотношением (4), в котором криволинейный интеграл не зависит от пути интегрирования, то выберем путь интегрирования так, как указано на рис.37.
Тогда
Отсюда
Последний интеграл берется вдоль отрезка прямой ММ1, параллельной оси Ох. На этом отрезке в качестве параметра можно принять координату х:
x = х, у = const, z = const.
Тогда dx = dx,dy = 0, dz = 0, так что
Применяя к интегралу в правой части (6) теорему о среднем, получаем
где величина ξ заключена между х и х + ∆х. Из формулы (7) вытекает, что
Так как ξ —► x при ∆x —» 0, то в силу непрерывности функции Р(х, у, z) получаем
Аналогично доказывается, что
Следствие:
Векторное поле является потенциальным тогда и только тогда, когда криволинейный интеграл в нем не зависит от пути.
Вычисление криволинейного интеграла в потенциальном поле
Теорема:
Интеграл в потенциальном поле а(М) равен разности значений, потенциала и(М) поля в конечной и начальной точках пути интегрирования,
Ранее былодоказано, что функция
является потенциалом поля.
В потенциальном поле криволинейный интеграл
не зависит от пути интегрирования. Поэтому, выбирая путь отточки М1 к точке М2 так, чтобы он прошел через точку Mo (рис. 38), получаем
или, меняя ориентацию пути в первом интеграле справа,
Так как потенциал поля определяется с точностью до постоянного слагаемого, то любой потенциал рассматриваемого поля можетбыть записан в виде
v(M) = u(M) + c, (10)
где с — постоянная.
Делая в формуле (10) замену u(M2) = v(M2) — с, и(М1) = v(M1) — с, получим для произвольного потенциала v(M) требуемую формулу
Пример:
В примере 1 было показано, что потенциалом поля радиус-вектора г является функция
Поэтому
где ri (i = 1,2) — расстояние от точки Mi(i = 1,2) до начала координат.
Вычисление потенциала в декартовых координатах
Пусть задано потенциальное поле
а(М) = Р(х, у, г)i + Q(x, у, z)j + R(x, у, z)k.
Ранее было показано, что потенциальная функция и(М) может быть найдена по формуле
Интеграл (11) удобнее всего вычислять так: зафиксируем начальную точку Мо(хо, yо, zо) и соединим ее с достаточно близкой текущей точкой M(x,y,z) ломаной М0М1М2М, звенья которой параллельны координатным осям, М0М1,||Ох, M1М2||Оу, М2М|| Oz (рис.39).
При этом на каждом звене ломаной изменяется только одна координата, что позволяет существенно упростить вычисления. В самом деле, на отрезке М0М1 имеем:
На отрезке М1М2:
х = const, dx = 0, у = у, dy = dy, z = z0 и dz = 0.
На отрезке M1M:
x = const, dx = 0, у = const, dy = 0, z = z и dz = dz.
Следовательно, потенциал u(M) равен
где x, у, z — координаты текущей точки на звеньях ломаной, вдоль которых ведется интегрирование.
Пример:
Доказать, что векторное поле
а = {у + z)i + (х + z)j + (х + y)k
является потенциальным, и найти его потенциал.
Проверим, будет ли поле вектора а(М) потенциально. С этой целью вычислим ротор поля. Имеем
Поле является потенциальным. Потенциал этого поля найдем с помощью формулы (12). Возьмем за начальную точку Mо начало координат О (так обычно поступают, если поле а(М) определено в начале координат). Тогда получим
Итак,
u(z, у, z) = ху + xz + yz + с,
где с — произвольная постоянная.
Потенциал этого поля можно найти и по-иному. По определению потенциал и(х, у, г) есть скалярная функция, для которой grad и = а. Это векторное равенство равносильно трем скалярным равенствам:
Интегрируя (13) по х, получим
где f(y,z) — произвольная дифференцируемая функция от у и z. Продифференцируем (16) по у:
откуда, учитывая (14), будем иметь
или
Проинтегрировав (17) по у, найдем
где F(z) — некоторая функция х. Подставив (18) в (16), получим
и(х, у, z) = ху + xz + у z + F(z).
Дифференцируя последнее равенство по z и учитывая соотношение (15), получим уравнение для F(z),
откуда = 0, так что F(z) = с = const. Итак,
u(x,y,z) = ху +yz + zx +с.
Оператор Гамильтона
Мы рассмотрели три основные операции векторного анализа: вычисление grad и для скалярного поля и = и(х, у, z) и div а и rot а для векторного поля а = а(x, у, z). Эти операции могут быть записаны в более простом виде с помощью символического оператора ∇ («набла»): (1)
Оператор ∇ (оператор Гамильтона) обладает как дифференциальными, так и векторными свойствами. Формальное умножение, например, умножение на функцию и(х, у), будем понимать как частное дифференцирование:
В рамках векторной алгебры формальные операции над оператором ∇ будем проводить так, как если бы он был вектором. Используя этот формализм, получим следующие основные формулы:
1, Если и = и(х, у, z) — скалярная дифференцируемая функция, то по правилу умножения вектора на скаляр получим
2. Если
a = P(x, y, z)i + Q(x, y, z)j + (x, y, z)k,
где P, Q, R — дифференцируемые функции, то по формуле для нахождения скалярного произведения получим
3. Вычисляя векторное произведение [ ∇, а], получим
Для постоянной функции и = с получим
∇c = 0,
а для постоянного вектора с будем иметь
( ∇, с) = 0 и [ ∇, с] = 0.
Из распределительного свойства для скалярного и векторного произведений получаем
( ∇, a + b) = ( ∇, а) + ( ∇, b),
Замечание:
Формулы (5) и (6) можно трактовать так же как проявление дифференциальных свойств оператора «набла»( ∇ — линейный дифференциальный оператор). Условились считать, что оператор ∇ действует на все величины, написанные за ним. В этом смысле, например,
( ∇, а) ≠ (а, ∇ ),
ибо ( ∇, а) = div а есть функция в то время как
— скалярный дифференциальный оператор.
Применяя оператор ∇ к произведению каких-либо величин, надо иметь в виду обычное правило дифференцирования произведения.
Пример:
Доказать, что
grad(u v) = v Brad и + u grad v. (7)
По формуле (2) с учетом замечания 1 получаем
∇(uv) = v∇u + u ∇v,
или
grad(uv) = v grad u + u grad v.
Чтобы отметить тот факт, что «набла» не действует на какую-либо величину, входящую в состав сложной формулы, эту величину отмечают индексом с («const»), который в окончательном результате опускается.
Пример:
Пусть и(x,y,z) — скалярная дифференцируемая функция, a(x,y,z) — векторная дифференцируемая функция. Доказать, что
div(ua) =u diva + (a, grad u). (8)
Перепишем левую часть (8) в символическом виде
div(ua) = ( ∇, ua).
Учитывая дифференциальный характер оператора ∇, получаем
Так как ис — постоянный скаляр, то его можно вынести за знак скалярного произведения, так что
(на последнем шаге мы опустили индекс с).
В выражении ( ∇, иас) оператор ∇ действует только на скалярную функцию и, поэтому
В итоге получаем
div(ua) = u div а + (a, grad и).
Замечание:
Используя формализм действий с оператором ∇ как с вектором, надо помнить, что ∇ не является обычным вектором — он не имеет ни длины, ни направления, так что, например, вектор ( ∇, а} не будет, вообще говоря, перпендикулярным вектору а (впрочем, для плоского поля а = Р(х, y)i + Q(x, y)j вектор
перпендикулярен плоскости хОу, а значит, и вектору а).
Не имеет смысла и понятие коллинеарности по отношению к символическому вектору ∇. Например, выражение [∇ φ, ∇ ψ] где φ и ψ — скалярные функции, формально напоминает векторное произведение двух кoллинеарных векторов, которое всегда равно нулю. Однако в общем случае это не имеет места. В самом деле, вектор ∇ φ = grad φ направлен по нормали к поверхности уровня φ = const, а вектор ∇ ψ = grad ψ определяет нормаль к поверхности уровня ψ = const. В общем случае эти нормали не обязаны быть коллинеарными (рис. 40). С другой стороны, в любом дифференцируемом скалярном поле φ (х, у, z) имеем [∇ φ, ∇ ψ] = 0.
Эта примеры показывают, что с оператором «набла» нужно обращаться с большой осторожностью и при отсутствии уверенности в полученном результате его следует проверить аналитическими методами.
Дифференциальные операции второго порядка. Оператор Лапласа
Дифференциальные операции второго порядка получаются в результате двукратного применения оператора ∇.
1, Пусть имеем скалярное поле и = и(x,y,z). В этом поле оператор ∇ порождает векторное поле
∇u = grad и.
В векторном поле grad и можно определить две операции:
( ∇, ∇u) = div grad u, (1)
что приводит к скалярному полю, и
[ ∇, ∇m] = rot grad u, (2)
что приводит к векторному полю.
2. Пусть задано векторное поле а = Pi + Qj + Rk. Тогда оператор (2) порождает в нем скалярное поле
(∇, а) = div a.
В скалярном поле div а оператор ∇ порождает векторное поле
∇ (∇,a) = grad div а. (3)
3. В векторном поле а = Pi + Qj + Rк оператор ∇ порождает также векторное поле
[∇, а] = rot a.
Применяя к этому полю снова оператор ∇, получим:
а) скалярное поле
(∇, [∇, а]) = div rot а, (4).
б) векторное поле
(∇, [∇, а]) = rot rot а. (5)
Формулы (1)-(5) определяют так называемые дифференциальные операции второго порядка.
Выберем в пространстве прямоугольную декартову систему координат Oxyz и рассмотрим каждую из формул (1)-(5) более подробно.
1, Предполагая, что функция и(х, у, z) имеет непрерывные вторые частные производные по х, у и z, получим
Символ
называется оператором Лапласа, или лапласианом. Его можно представить как скалярное произведение оператора Гамильтона ∇ на самого себя, т.е.
Оператор ∆ (дельта) играет важную роль в математической физике. Уравнение (6)
называется уравнением Лапласа. С его помощью описывается, например, стационарное распределение тепла.
Скалярное поле и(х, у, z), удовлетворяющее условию ∆и = 0, называется лапла-совым или гармоническим полем.
Например, скалярное поле и = 2х2 + Зу — 2x2 является гармоническим во всем трехмерном пространстве: из того, что
2. Пусть функция u(z, у, z) имеет непрерывные частные производные второго порядка включительно. Тогда
rot grad u ≡ 0. (7)
В самом деле, действуя формально, получим
rot grad и = [∇, ∇u] = [ ∇, ∇ ] u = 0,
ибо [∇, ∇] = 0 как векторное произведение двух одинаковых «векторов».
Тот же результат можно получить, используя выражения градиента и ротора в декартовых координатах
3. Пусть задано векторное поле
а = Р(х, у, z)i + Q(x, у, z)j + R(x, у, z)k,
координаты которого P, Q, R имеют непрерывные частные производные второго порядка. Тогда получим
4. При тех же условиях, что и в пункте 3, имеем (9)
Это соотношение уже было доказано ранее путем непосредственных вычислений. Здесь мы приведем его формальное доказательство, используя известную формулу из векторной алгебры
(А, [В, С]) = (С,[А,В])= (В, [С, А]).
Имеем
div rot а = (∇, [∇, а]) = (а, [∇, ∇]) = О,
так как [∇, ∇] = 0 как векторное произведение двух одинаковых «векторов».
5. Покажем, наконец, что при тех же условиях, что и ранее,
rot rot а = grad div а — ∆а. (10)
Так как
rot rot а = [ ∇, [∇,a]),
то, полагая в формуле для двойного векторного произведения [А, [В, С]] = В(А, С) — (А, В)С,
А = ∇, B = ∇, С = а,
получим
Но ( ∇, а) = div а и ( ∇, ∇) = ∆. Поэтому окончательно будем иметь
rot rot а = grad div а — ∆а,
где grad div а выражается по формуле (8), а ∆а для вектора а = Pi + Qj + Rk надо понимать так:
∆а = ∆Р • i + ∆Q • j + ∆R • k.
В заключение приведем таблицу дифференциальных операций второго порядка.
Заштрихованные прямоугольники означают, что соответствующая операция не имеет смысла (например, градиент от rot а).
Понятие о криволинейных координатах
Во многих задачах бывает удобно определять положение точки простр анства не декартовыми координатами (х, у, z), а тремя другими числами (q1, q2, q3), более естественно связанными с рассматриваемой частной задачей.
Если задано правило, согласно которому каждой точке М пространства отвечает определенная тройка чисел (q1, q2, q3) и, обратно, каждой такой тройке чисел отвечает единственная точка М, то говорят, что в пространстве задана криволинейная координатная система. В этом случае величины q1, q2, q3 называют криволинейными координатами точки М.
Координатными поверхностями в системе криволинейных координат q1, q2, q3 называются поверхности
На координатных поверхностях одна из координат сохраняет постоянное значение. Линии пересечения двух координатных поверхностей называются координатными линиями.
В качестве примеров криволинейных координат рассмотрим цилиндрические и сферические координаты.
Цилиндрические координаты
В цилиндрических координатах положение точки М в пространстве определяется тремя координатами:
Координатные поверхности:
р = const — круговые цилиндры с осью Оz;
φ = const — полуплоскости, примыкающие к оси Oz;
z = const — плоскости, перпендикулярные оси Oz (рис. 41).
Координатные линии:
1) линии (р) — лучи, перпендикулярные оси Oz и имеющие начало на этой оси, т. е. линии пересечения координатных поверхностей φ = const, z = const;
2) линии (φ) — окружности с центрами на оси Oz, лежащие в плоскостях, перпендикулярных оси Oz;
3) линии (z) — прямые, параллельные оси Oz.
Связь декартовых координат точки (х, у, z) с цилиндрическими координатами (р, φ, z) задается формулами
x = p cos φ, y = p sin φ, z = z. (2)
Сферические координаты
В сферических координатах положение точки М в пространстве определяется следующими координатами:
Координатные поверхности (рис. 42):
r = const — сферы с центром в точке О;
θ = const — круговые полуконусы с осью Oz;
φ = const — полуплоскости, примыкающие к оси Oz.
Координатные линии:
1) линии (г) — лучи, выходящие из точки О;
2) линии (θ) — меридианы на сфере;
3) линии (φ) — параллели на сфере.
Связь декартовых координат (х, у, z) точки М с ее сферическими координатами (r, θ, φ) задается формулами
х = r cos φ sin θ,
у = r sin φ sin θ, (4)
z = r cos θ.
Введем единичные векторы e1, е2, е3 (орты), направленные по касательным к координатным линиям(q1),(q2),(q3)в тoчке М в сторону возрастания переменных q1,q2,q3 соответственно.
Определение:
Система криволинейных координат называется ортогональной, если в каждой точке М орты e1, е2, е3 попарно ортогональны.
В такой системе ортогональны и координатные линии, и координатные поверхности.
Примерами ортогональных криволинейных координат служат системы цилиндрических и сферических координат. Мы ограничимся рассмотрением только ортогональных систем координат.
Пусть r = r(q1, q2, q3) — радиус-вектор точки М. Тогда можно показать, что
(5)
где
— коэффициенты Ламэ данной криволинейной системы координат. Вычислим коэффициенты Ламэ для цилиндрических координат
Так как х = р cos φ, у = р sin φ, z = z, то
(6)
Аналогично для сферических координат имеем
(7)
Величины
являются дифференциалами длин дуг соответствующих координатных линий.
Основные операции векторного анализа в криволинейных координатах
Дифференциальные уравнения векторных линий
Рассмотрим поле вектора
Уравнения векторных линий в криволинейных координатах q1,q2, q3 имеют вид
В цилиндрических координатах (q1 = р, q2= φ, q3 = z)
(1)
в сферических координатах (q1 = r, q2 = θ, q3 = φ)
(2)
Градиент в ортогональных координатах
Пусть и = u(q1, q2, q3) — скалярное пoле. Тогда
В цилиндрических координатах (q1 = р, q2 = φ, q3 = z)
(3)
в сферических координатах (q1 = r, q2 = θ, q3 = φ) (4)
Ротор в ортогональных координатах
Рассмотрим векторное поле
и вычислим rot а. Имеем
В цилиндрических координатах
(5)
в сферических координатах
(6)
Дивергенция в ортогональных координатах
Дивергенция div а векторного поля
вычисляется по формуле
(7)
В цилиндрических координатах
в цилиндрических координатах
в сферических координатах
Применяя формулу (7) к единичным векторам е1, е2, е3, получим
Вычисление потока в криволинейных координатах
Пусть S — часть координатной поверхности q1 = с = const, ограниченная координатными линиями
Тогда поток вектора
через поверхность S в направлении вектора e1 вычисляется по формуле
(8)
Аналогично вычисляется поток через часть поверхности q2 = с, а также через часть поверхности д3 = с, где с = const.
Пример:
Найти поток П векторного поля
через внешнюю сторону верхней полусферы S радиуса R с центром в начале координат.
Полусфера S есть часть координатной поверхности r = const, а именно r = R. На полусфере S имеем
Учитывая, что в сферических координатах
по формуле (8) найдем
Вычисление потенциала в криволинейных координатах
Пусть в некоторой области Ω задано потенциальное векторное поле
т. e. rot а = 0 в области Ω.
Для нахождения потенциала и = и(q1, q2, q3) этого векторного поля запишем равенство а(М) = grad u(M) в следующем виде:
Отсюда следует, что
(9)
Интегрируя систему дифференциальных уравнений с частными производными (9), найдем искомый потенциал
где с — произвольная постоянная.
В цилиндрических координатах
система (9) принимает вид
В сферических координатах
система (9) имеет вид
Пример:
Найти потенциал векторного поля, заданного в цилиндрических координатаx
Убедимся, что rot a = 0. По формуле (S) получим
т.е. данное поле потенциально.
Искомый потенциал u = и(р, φ, z) является решением следующей системы дифференциальных уравнений с частными производными (см. формулу (10)):
Интегрированием по р из первого уравнения находим
Дифференцируя соотношение (11) no φ и используя второе уравнение, получим
или = 0, откуда с = c1(z). Таким образом,
Дифференцируя это соотношение по z и используя третье уравнение, получим
или c1`(z) =0, откуда c1(z) = с. Итак, потенциал данного поля
Линейный интеграл и циркуляция в ортогональных криволинейных координатах
Пусть векторное поле
определено и непрерывно в области Ω изменения ортогональных криволинейных координат q1, q2, q3. Так как дифференциал радиус-вектора r любой точки M(q1, q2, q3) ∈ Ω выражается формулой
то криволинейный интеграл вектора а(М) по ориентированной гладкой или кусочно-гладкой кривой L ⊂ Ω будет равен
(13)
В частности, для цилиндрических координат (q1 = р, q2 = φ, q3 = z, Н1 = 1, Н2 = р, Н3=1) будем иметь
Отсюда по формуле (13) получим
(14)
Аналогично для сферических координат (q1 = r, q2 = θ, q3 = φ, Н1 = 1, Н2 = r, H3 = r sin θ будем иметь
Отсюда по формуле (13) получим
(15)
Если кривая L замкнута (начальная и конечная точки кривой L совпадают), то циркуляция Ц векторного поля а (М) в криволинейных координатах q1, q2, q3 вычисляется по формуле (13), а в случае цилиндрических или сферических координат — по формулам (14) или (15) соответственно.
Пример:
Вычислить циркуляцию векторного поля, заданного в цилиндрических координатах
по замкнутой кривой L,
Координаты данного вектора равны соответственно
Контур L представляет собой замкнутую кривую, расположенную в плоскости z = 0 (рис. 43).
Подставляя координаты данного вектора в формулу .(14), получим
На кривой L имеем
Искомая циркуляция будет равна
Оператор Лапласа в ортогональных координатах
Используя формулы (16) и (17), для оператора Лапласа ∆ получим следующее выражение:
В цилиндрических координатах
1) получим
В сферических координатах
будем иметь
Пример:
Найти все решения уравнения Лапласа ∆и = 0, зависящие только от расстояния r.
Так как искомое решение и должно зависеть только от расстояния точки М от начала координат г, т. е. и = и (r), то уравнение Лапласа ∆и = 0 в сферических координатах будет иметь вид
Отсюда так что
где С1 и С2 — постоянные.
Решение заданий и задач по предметам:
- Математика
- Высшая математика
- Математический анализ
- Линейная алгебра
Дополнительные лекции по высшей математике:
- Тождественные преобразования алгебраических выражений
- Функции и графики
- Преобразования графиков функций
- Квадратная функция и её графики
- Алгебраические неравенства
- Неравенства
- Неравенства с переменными
- Прогрессии в математике
- Арифметическая прогрессия
- Геометрическая прогрессия
- Показатели в математике
- Логарифмы в математике
- Исследование уравнений
- Уравнения высших степеней
- Уравнения высших степеней с одним неизвестным
- Комплексные числа
- Непрерывная дробь (цепная дробь)
- Алгебраические уравнения
- Неопределенные уравнения
- Соединения
- Бином Ньютона
- Число е
- Непрерывные дроби
- Функция
- Исследование функций
- Предел
- Интеграл
- Двойной интеграл
- Тройной интеграл
- Интегрирование
- Неопределённый интеграл
- Определенный интеграл
- Криволинейные интегралы
- Поверхностные интегралы
- Несобственные интегралы
- Кратные интегралы
- Интегралы, зависящие от параметра
- Квадратный трехчлен
- Производная
- Применение производной к исследованию функций
- Приложения производной
- Дифференциал функции
- Дифференцирование в математике
- Формулы и правила дифференцирования
- Дифференциальное исчисление
- Дифференциальные уравнения
- Дифференциальные уравнения первого порядка
- Дифференциальные уравнения высших порядков
- Дифференциальные уравнения в частных производных
- Тригонометрические функции
- Тригонометрические уравнения и неравенства
- Показательная функция
- Показательные уравнения
- Обобщенная степень
- Взаимно обратные функции
- Логарифмическая функция
- Уравнения и неравенства
- Положительные и отрицательные числа
- Алгебраические выражения
- Иррациональные алгебраические выражения
- Преобразование алгебраических выражений
- Преобразование дробных алгебраических выражений
- Разложение многочленов на множители
- Многочлены от одного переменного
- Алгебраические дроби
- Пропорции
- Уравнения
- Системы уравнений
- Системы уравнений высших степеней
- Системы алгебраических уравнений
- Системы линейных уравнений
- Системы дифференциальных уравнений
- Арифметический квадратный корень
- Квадратные и кубические корни
- Извлечение квадратного корня
- Рациональные числа
- Иррациональные числа
- Арифметический корень
- Квадратные уравнения
- Иррациональные уравнения
- Последовательность
- Ряды сходящиеся и расходящиеся
- Тригонометрические функции произвольного угла
- Тригонометрические формулы
- Обратные тригонометрические функции
- Теорема Безу
- Математическая индукция
- Показатель степени
- Показательные функции и логарифмы
- Множество
- Множество действительных чисел
- Числовые множества
- Преобразование рациональных выражений
- Преобразование иррациональных выражений
- Геометрия
- Действительные числа
- Степени и корни
- Степень с рациональным показателем
- Тригонометрические функции угла
- Тригонометрические функции числового аргумента
- Тригонометрические выражения и их преобразования
- Преобразование тригонометрических выражений
- Комбинаторика
- Вычислительная математика
- Прямая линия на плоскости и ее уравнения
- Прямая и плоскость
- Линии и уравнения
- Прямая линия
- Уравнения прямой и плоскости в пространстве
- Кривые второго порядка
- Кривые и поверхности второго порядка
- Числовые ряды
- Степенные ряды
- Ряды Фурье
- Преобразование Фурье
- Функциональные ряды
- Функции многих переменных
- Метод координат
- Гармонический анализ
- Вещественные числа
- Предел последовательности
- Аналитическая геометрия
- Аналитическая геометрия на плоскости
- Аналитическая геометрия в пространстве
- Функции одной переменной
- Высшая алгебра
- Векторная алгебра
- Векторы
- Скалярное произведение векторов
- Векторное произведение векторов
- Смешанное произведение векторов
- Операции над векторами
- Непрерывность функций
- Предел и непрерывность функций нескольких переменных
- Предел и непрерывность функции одной переменной
- Производные и дифференциалы функции одной переменной
- Частные производные и дифференцируемость функций нескольких переменных
- Дифференциальное исчисление функции одной переменной
- Матрицы
- Линейные и евклидовы пространства
- Линейные отображения
- Дифференциальные теоремы о среднем
- Теория устойчивости дифференциальных уравнений
- Функции комплексного переменного
- Преобразование Лапласа
- Теории поля
- Операционное исчисление
- Системы координат
- Рациональная функция
- Интегральное исчисление
- Интегральное исчисление функций одной переменной
- Дифференциальное исчисление функций нескольких переменных
- Отношение в математике
- Математическая логика
- Графы в математике
- Линейные пространства
- Первообразная и неопределенный интеграл
- Линейная функция
- Выпуклые множества точек
- Система координат