Использование полярных, цилиндрических и сферических координат для вычисления кратных интегралов

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

Использование полярных координат

Из курса аналитической геометрии известны следующие соотношения между декартовыми и полярными координатами: $x = r\cos\phi,\quad y = r\sin\phi\quad(*)$.
При этом, $r \geq 0, 0 \leq \phi <2\pi$. Рассмотрим вспомогательную плоскость $RO\Phi$, где $r$ и $\phi$ являются декартовыми координатами, и определим на ней множество точек $G$, такое, что: $G = \{(r, \phi)| r > 0, 0 \leq \phi < 2\pi\}$.

Тогда формулы $(*)$ определяют непрерывно дифференцируемое отображение $F : G \to \widetilde{XOY}$, где $\widetilde{XOY} = XOY \setminus\{(0, 0)\}$.

По определению полярных координат, в декартовой системе координат $XOY$ $r$ задает радиус окружности с центром в начале координат, а $\phi$ определяет луч, исходящий из центра координат, такой что угол между лучом и положительным направлением оси $OX$ равен $\phi$. С геометрической точки зрения очевидно, что они пересекаются в единственной точке.

Таким образом, любую точку $P = (x_0, y_0)$ из $\widetilde{XOY}$ можно однозначно определить пересечением луча, направленного под углом $\phi_0$ и окружности радиусом $r_0$, и тогда точка $P’ = (r_0, \phi_0)$ будет единственным прообразом $P$ в $G$. Очевидно, что любой элемент из $G$ служит прообразом, и что двум различным точкам из $G$ будут соответствовать 2 различные точки из $\widetilde{XOY}$. Таким образом, отображение $F$ между точками плоскостей $G$ и $\widetilde{XOY}$ взаимно однозначное:

kursach

Якобиан полученного отображения будет равен:
$J_F = \begin{array}{|cc|} \frac{\partial x}{\partial r} & \frac{\partial x}{\partial \phi} \\ \frac{\partial y}{\partial r} & \frac{\partial y}{\partial \phi} \end{array} = \begin{array}{|cc|} \cos{\phi}& -r\sin{\phi} \\ \sin{\phi} & r\cos{\phi} \end{array} = r$

Теперь рассмотрим множество точек $G’$, полученное добавлением к множеству $G$ отрезка $r = 0$, т.е. $G’ = \{(r, \phi)| r \geq 0, 0 \leq \phi < 2\pi\}.$ $G’$ уже является прообразом всей плоскости $XOY$, но на отрезке $r = 0, 0 \leq \phi < 2\pi$ не достигается взаимная однозначность, а $\left|J_F\right| = 0$. Обратим внимание, что его Жорданова мера равна нулю.

Наконец, пусть дана область $\Omega \subset XOY$ и функция $f$, непрерывная на измеримом множестве $\overline{\Omega}$. Ее прообразом при отображении $F$, заданного формулами $(*)$, будет некоторая область $\Omega’ \subset G’$. Если область $\Omega$ не содержит точки O — начала координат, то выполнены все условия теоремы о замене переменной в кратных интегралах, и справедлива формула:
$$\iint\limits_{\Omega} f(x, y)\,dxdy = \iint\limits_{\Omega’} f(r\cos{\phi}, r\sin{\phi})r\,drd\phi$$
Если же точка $O \in \Omega$, то взаимная однозначность и не обращение якобиана в нуль не выполняются на множестве $r = 0$, что не влияет на справедливость данной формулы (следует из замечания к указанной теореме).

Пример №1

Вычислить интеграл:
$\iint\limits_{\Omega}(x^2 + y^2)\,dxdy, \Omega = \{(x, y)| y \geq 0, x^2 + y^2 \leq a^2\}.$
Заметим, что в полярных координатах полукруг [latex]\Omega[/latex] будет представлять из себя более простую область интегрирования:

example

Поэтому, воспользуемся формулой замены переменной и перейдем к полярным координатам:
$\iint\limits_{\Omega}(x^2 + y^2)\,dxdy = \iint\limits_{\Omega’}r^2r\,drd\phi = \int\limits_0\limits^{\pi}\,d\phi\int\limits_0\limits^ar^3\,dr = $ $ \int\limits_0\limits^{\pi}\frac{a^4}{4}\,d\phi = \frac{\phi a^4}{4}|_0^{\pi} = \frac{\pi a^4}{4}$.

[свернуть]

Использование цилиндрических и сферических координат

Рассмотрим теперь пространство $\mathbb{R}^3$, в котором задана декартова система координат $OXYZ$. Цилиндрические координаты связанны с декартовыми следующим образом:
$x = r\cos\phi,\quad y = r\sin\phi,\quad z = t\quad(**),$
где $r \geq 0, 0 \leq \phi <2\pi, t \in \mathbb{R}$ (величины $r$ и $\phi$ для любой точки $A = (x, y, z)$ определяются таким же образом, как и в полярных координатах для ее проекции $P’ = (x, y, 0)$ на $XOY$). Теперь, аналогично случаю с полярными координатами, рассмотрим вспомогательное пространство $OR\Phi T$, где $r, \phi, t$ — декартовы координаты, а в нем — множество точек $G = \{(r, \phi, t)| r \geq 0, 0 \leq \phi <2\pi, t \in \mathbb{R}\}$.

Отображение $F : G \to OXYZ$, определяемое формулами $(**)$, является непрерывно дифференцируемым.
$J_F = \begin{array}{|ccc|} \frac{\partial x}{\partial r} & \frac{\partial x}{\partial \phi} & \frac{\partial x}{\partial t} \\ \frac{\partial y}{\partial r} & \frac{\partial y}{\partial \phi} & \frac{\partial y}{\partial t} \\ \frac{\partial z}{\partial r} & \frac{\partial z}{\partial \phi} & \frac{\partial z}{\partial t} \end{array} = \begin{array}{|ccc|} \cos{\phi}& -r\sin{\phi} & 0 \,\\ \sin{\phi} & r\cos{\phi} & 0 \,\\ 0 & 0 & 1\,\end{array} = r$

Очевидно, что как и в случае с полярными координатами, отображение $F$ — взаимно однозначное, и его якобиан не равен нулю. Данные условия не выполняются только при $r = 0$, т.е. на множестве $L = \{(r, \phi, t)| r = 0, 0 \leq \phi <2\pi, t \in \mathbb{R}\}$. Пересечение такого множества с любым другим ограниченным множеством есть ограниченное линейное множество, и жорданова мера этого пересечения равна нулю.

Тогда, если дана область $\Omega \subset OXYZ$, и функция $f$ непрерывна на измеримом множестве $\overline{\Omega}$, а $\Omega’ \subset G$ — прообраз данной области при отображении $F$, то выполнены все условия теоремы о замене, и справедлива следующая формула:
$$\iiint\limits_{\Omega} f(x, y, z)\,dxdydz = \iiint\limits_{\Omega’} f(r\cos{\phi}, r\sin{\phi}, t)r\,drd\phi dt$$

Наконец, рассмотрим сферические координаты, связанные с декартовыми следующими соотношениями: $x = r\cos{\phi} \cos{\psi},\quad y = r\sin{\phi} \cos{\psi},\quad z = r\sin{\psi}\quad (***),$
где $r \geq 0, 0 \leq \phi < 2\pi, -\frac{\pi}{2} \leq \psi \leq \frac{\pi}{2}$. Введем вспомогательное пространство $OR\Phi\Psi$, где $r, \phi, \psi$ — декартовы координаты, а в нем рассмотрим множество точек $G = \{(r, \phi, \psi)| r \geq 0, 0 \leq \phi < 2\pi -\frac{\pi}{2} \leq \psi \leq \frac{\pi}{2}\}$.

Отображение $F : G \to OXYZ$, определяемое формулами $(***)$, непрерывно дифференцируемо.
$J_F = \begin{array}{|ccc|} \frac{\partial x}{\partial r} & \frac{\partial x}{\partial \phi} & \frac{\partial x}{\partial \psi} \\ \frac{\partial y}{\partial r} & \frac{\partial y}{\partial \phi} & \frac{\partial y}{\partial \psi} \\ \frac{\partial z}{\partial r} & \frac{\partial z}{\partial \phi} & \frac{\partial z}{\partial \psi} \end{array} = \begin{array}{|ccc|} \cos{\phi}\cos{\psi}& -r\sin{\phi}\cos{\psi} & -r\cos{\phi}\sin{\psi} \,\\ \sin{\phi}\cos{\psi} & r\cos{\phi}\cos{\psi} & -r\sin{\phi}\sin{\psi} \,\\ \sin{\psi} & 0 & r\cos{\psi}\,\end{array} = $ $ r^2\cos{\psi}$.

Взаимная однозначность данного отображения устанавливается по тем же рассуждениям, что и в предыдущих двух случаях, и не выполняется только при $r = 0, \psi = -\frac{\pi}{2}, \psi = \frac{\pi}{2}$, когда и якобиан равен нулю. Однако любое подмножество множества, задаваемого такими равенствами, будет представлять собой ограниченную часть плоскости с жордановой мерой нуль в пространстве $OXYZ$, что не помешает совершить замену.

Тогда, при соответствующих условиях, справедлива формула замены переменной ($\Omega \subset OXYZ, \Omega’ \subset G$):
$$\iiint\limits_{\Omega} f(x, y, z)\,dxdydz = $$ $$\iiint\limits_{\Omega’} f(r\cos{\phi}\cos{\psi}, r\sin{\phi}\cos{\psi}, r\sin{\psi})r^2\cos{\psi}\,drd\phi d\psi$$

Пример №2

Вычислить интеграл $\iiint\limits_{\Omega} e^{{(x^2 + y^2 + z^2)}^{\frac{3}{2}}}\,dxdydz$, где граница области $\Omega$ задается уравнением $x^2 + y^2 + z^2 = 1$.

Область интегрирования представляет собой шар радиуса [latex]1[/latex] с центром в начале координат. Следовательно, будет удобно воспользоваться переходом к цилиндрической системе координат. В ней новая область интегрирования [latex]\Omega'[/latex] будет определятся следующими неравенствами: $0 \leq \phi \leq 2\pi,\quad -\frac{\pi}{2} \leq \psi \leq \frac{\pi}{2},\quad 0 \leq r \leq 1$. Воспользуемся формулой замены переменной для сферических координат:
$\iiint\limits_\Omega e^{(x^2 + y^2 + z^2)^{\frac{3}{2}}}\,dxdydz = $ $\iiint\limits_{\Omega’} e^{r^{2^{\frac{3}{2}}}}r^2\cos{\psi} \,drd\phi d\psi = \int\limits_0\limits^{2\pi}\,d\phi\int\limits_0\limits^1e^{r^3}r^2\,dr\int\limits_{-\frac{\pi}{2}}\limits^{\frac{\pi}{2}}\cos{\psi}\,d\psi = $ $\int\limits_0\limits^{2\pi}\,d\phi\int\limits_0\limits^1e^{r^3}r^2\,dr \cdot (\sin{\psi})|_{-\frac{\pi}{2}}^{\frac{\pi}{2}} = 2\int\limits_0\limits^{2\pi}\,d\phi\int\limits_0\limits^1\frac{1}{3}e^{r^3}\,d(r^3) = $ $\frac{2}{3}\int\limits_0\limits^{2\pi}\,d\phi \cdot e^{r^3}|_0^1 = \frac{2}{3}(e — 1)\int\limits_0\limits^{2\pi}\,d\phi = \frac{2}{3}(e — 1) \cdot \phi|_0^{2\pi} = \frac{4\pi}{3}(e — 1)$

[свернуть]

Тест: Использование полярных, цилиндрических и сферических координат для вычисления кратных интегралов

Для закрепления материала, рекомендуется пройти тест по данной теме.


Таблица лучших: Переход к полярным, цилиндрическим и сферическим координатам при вычислении кратных интегралов

максимум из 7 баллов
Место Имя Записано Баллы Результат
Таблица загружается
Нет данных

Использование полярных, цилиндрических и сферических координат для вычисления кратных интегралов: 2 комментария

  1. — Почему не разбить на три статьи? Аргументируйте.
    — Дифференциал в подинтегральном выражении лучше отделять маленьким пробелом (\,).

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *