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

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

350

ЛЕКЦИЯ 3.9. ВЫЧИСЛЕНИЕ КРИВОЛИНЕЙНОГО, ДВОЙНОГО И ТРОЙНОГО ИНТЕГРАЛОВ ПУТЕМ СВЕДЕНИЯ К ЛИНЕЙНОМУ

3.9.1 Уравнения линий в полярной системе координат

Полярная система координат состоит: из начала отсчета или полюса (точка О) (Рис.3.9.1а) и полярной оси Ох. Возьмем произвольную точку на плоскости – М. Расстояние от этой точки до полюса О называют полярным радиусом, обозначают – r. Угол, на который нужно повернуть полярную ось до совпадения с полярным радиусом, называют полярным углом, обозначают – φ. Таким образом, положение любой точки на плоскости в полярной системе, так же как и в декартовой, однозначно задается двумя числами: расстоянием до полюса – r и полярным углом – φ.

Совместим начало отсчета декартовой и полярной систем и найдем координаты одной и той же точки М в этих системах (Рис.3.9.2а). Тогда, декартовы координаты

точки М будут связаны с ее полярными координатами соотношениями: x r cos ,

y r sin .

Из последних равенств легко получить выражения через х и y для полярных координат r и φ точки М:

r x2 y2

arctg xy

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

Аналогично, уравнение линии в полярной системе отсчета должно связывать полярные координаты точек (r,φ), принадлежащих данной линии, т.е.: r = f(φ) или r = r(φ), и также является функцией одного переменного φ. Чтобы построить график линии, заданной уравнением r = r(φ), нужно для каждого значения полярного угла φ(выбранного произвольно) найти соответствующее значение r из уравнения r = r(φ), и отложить его на полярной оси, повернутой на соответствующий угол φ. Соединяя построенные таким образом точки, получают уравнение линии.

Пример 1. Построить график линии, заданной уравнением:

351

r (1 cos )

Решение. Будем изменять полярный угол φ в пределах от 0 до 2π с шагом 4 . Соответствующие значения полярного радиуса r занесем в таблицу.

кардиоидой.

Взаключении скажем о координатных линиях, т.е. о линиях, на которых одна или другая координата сохраняют постоянное значение.

Вдекартовой системе координатные

линии

x= const и y = const образуют два семейства прямых, параллельных одной или другой из осей координат. В полярной системе линии r

=const представляют собой семейство концентрических окружностей, а линии φ = const – семейство лучей, выходящих из полюса.

352

3.9.2 Вычисление криволинейного интеграла

Криволинейный интеграл по кривой L от заданной на ней функции f(p) записывают следующим образом:

f ( p)dl , (3.9.1а)

L

где L – конечная часть линии, определяющая область интегрирования; dℓ - длина (мера) бесконечно малых элементов линии;

f(p) – подынтегральная функция, значения которой берутся в точках L. Криволинейный интеграл вычисляют путем сведения его к линейному

интегралу, при этом должно быть задано:

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

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

а. Пусть линия L – является плоской кривой и задана в декартовой системе координат уравнением y = y(x), которое представляет собой непрерывную и дифференцируемую функцию на интервале [α,b]. Тогда часть ли-

нии между точками А и В (обозначим ее LAB) не будет иметь изломов.

В точках LAB задана некоторая функция

z f ( p) f (x, y)

(например, плотность массы или заряда), график которой схематически изображен на рисунке 3.9.4а. Эту функцию мы интегрируем, и ее не следует путать с уравнением линии y = y(x).

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

353

Как следует из формулы (3.9.2а), чтобы свести криволинейный интеграл по плоской кривой к линейному, нужно проделать две операции:

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

во-вторых, дифференциал длины dℓ следует также выразить через уравнение линии по формуле

dl 1 yx2 dx

При этом нужно помнить, что элемент длины dℓ выражается через уравнение линии L только в том случае, когда это уравнение является диф-

ференцируемой функцией на интервале [α,b], т.е. линия LAB гладкая и не имеет изломов.

б. А теперь рассмотрим криволинейный интеграл:

f (x, y)dl

LAB

по линии LAB, которая задана параметрически уравнениями: x = x(t)

y= y(t).

Вее начальной точке А(x1,y1) – параметр принимает значения – t1, в конечной В(x1,y1) – t2. Найдем выражение для dℓ через параметрическое уравне-

Заменим у подынтегральной функции f(x,y) переменные x и y их зависимостью от параметра t из уравнения линии x = x(t), y = y(t), а вместо дифференциала dℓ подставим найденное для него выражение через это же уравнение. В результате получим линейный интеграл от функции одного переменного t:

в. И наконец, пусть линия LAB задана в полярной системе координат

уравнением:

354

r = r(φ)

Полярные координаты начальной и конечной точек линии соответственно равны

A(r1, 1), B(r2 , 2 ) .

Выразим дифференциал длины линии dℓ через ее уравнение r = r(φ).

Полярные координаты точек, принадлежащих линии LAB, будут связаны с их декартовыми координатами следующими отношениями:

x r( cos (3.9.4а) y r( sin

где r(φ) – полярное уравнение LAB. Полярный угол φ можно рассматривать как параметр, а соотношение (3.9.4а) – параметрическими уравнениями

линии LAB, тогда:

Найдем производные x и y :

x r cos r sin y r sin r cos

Возьмем в квадрат обе части последних равенств и сложим: x 2 y 2 r 2 r2 (3.9.6а)

Подставим в выражение для dℓ (3.9.5а) вместо x 2 y 2 правую часть равенства (3.9.6а), получим:

dl r 2 r2 d (3.9.7а)

Чтобы перейти к линейному интегралу, нужно заменить дифференциал длины dℓ через уравнение линии r = r(φ) по формуле (3.9.7а), а переменные y и x у функции f(x,y) – их зависимостью от φ из соотношений (3.9.4а):

Сведение криволинейного интеграла к линейному во всех трех случаях по идее очень близко к замене переменной в линейном интеграле. Однако, нужно иметь в виду одно отличие. После замены переменной в линейном интеграле может оказаться, что нижний предел больше верхнего. При вычислении же криволинейного интеграла всегда нижний предел должен быть меньше верхнего. Это обусловлено тем, что элемент длины dℓ не может быть отрицательным, поэтому в формулах

355

dx, dt, dφ тоже всегда больше нуля, а это значит, что при переходе от криволинейного интеграла к линейному переменная интегрирования пробегает свои значения в сторону возрастания.

Пример 2. Найти массу одной четверти окружности x2 + y2 = α2, расположенной в первом квадранте, если плотность массы задана функцией ρ = xy.

Решение. Обозначим массу кривой – М. Так как кривая неоднородна, ее масса равна криволинейному интегралу:

М (x, y)dl xydl

найденное для него значение, получим:

Масса четверти данной окружности зависит от ее радиуса α.

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

L

арка циклоиды:

x (t sin t)

y (1 cost)

Решение. Уравнение кривой L задано параметрически. Найдем значение параметра t в начальной и конечной точках первой арки циклоиды. Для

этого построим ее график. Будем изменять параметр t через интервал 2 . Соответствующие значения координат x и y занесем в таблицу:

График первой арки циклоиды (Рис.3.9.6а) получается при изменении параметра t от 0 до 2π.

Циклоида представляет собой траекторию точки окружности, катящейся без скольжения по оси Ох.

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

кривой:

xt (1 cost)

yt sin t

Подставляя производные в выражение для dℓ, получим

Заменим у подынтегральной функции переменную y ее зависимостью от t из уравнения циклоиды и подставим найденное значение dℓ:

2 2 ydl 2 32 [t sin t]02 4 32 .

L

Окончательный результат получен по формуле Ньютона-Лейбница. Пример 4. Вычислить криволинейный интеграл (x y)dl ,

L

где L: r = α сosφ.

Решение. Найдем пределы изменения полярного угла φ. Построим график линии L.

357

r cos

Уравнение r = αcosφ задает окружность с радиусом 2 и центром, смещенным от полюса по полярной оси вправо на 2 (Рис.3.9.7а). При этом по-

лярный угол φ изменяется от 0до π.

Перейдем от криволинейного интеграла к линейному по формуле:

Подставляя вместо r(φ) уравнение линии и производную от него, полу-

чим:

2 cos2 d 2 cos sin d

Криволинейный интеграл будет равен:

358

В разобранных примерах мы пытались показать, как осуществляют переход от криволинейного интеграла по длине дуги к линейному в случаях различного задания уравнения линии L, по которой ведется интегрирование. Помимо криволинейных интегралов по длине дуги рассматривают криволинейные интегралы по координатам (второго рода). Эти интегралы подробно разобраны в разделе «Теория поля». Поэтому здесь на них мы останавливаться не будем.

3.9.3. Объем тел с известным поперечным сечением

Двойной интеграл по плоской области D, от заданной на ней функции z f (x, y) записывают так:

f (x, y)ds,

D

где ds – мера бесконечно малых элементов области D. Вычисление двойного интеграла сводят к последовательному вычислению двух линейных интегралов по переменным x и y. Чтобы получить формулу перехода к линейным интегралам, решим задачу о вычислении объема тела с известным поперечным сечением.

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

S = S(x).

Нужно найти объем этого тела. Спроектируем крайние точки тела на ось Ox. Пусть их координаты равны х = α, х = b.

Далее, разобьем тело плоскостями, перпендикулярными оси Ox, на n – частей.

Объем i-ой части (Рис. 3.9.1) можно приближенно принять за объем прямого

цилиндра с высотой xi и основанием S(xi):

Vi s(xi ) xi

Тогда объем всего тела приближенно будет равен сумме объемов элементарных частей: