Формула замены переменной в кратном интеграле

Теорема (формула замены переменной в кратном интеграле)

Пусть отображение $F : \Omega \to \mathbb{R}^n$, где $\Omega \subset \mathbb{R}^n$ — открытое множество, заданное при помощи непрерывно дифференцируемых функций $x_i = \phi_i(u_1, \ldots, u_n), i = 1, \ldots, n$, является взаимно однозначным и удовлетворяет следующим условиям:

  1. производные $\frac{\partial \phi_i}{\partial u_i}$ ограничены в $\Omega$;
  2. производные $\frac{\partial \phi_i}{\partial u_i}$ равномерно непрерывны в $\Omega$;
  3. якобиан $J(u)$ отображения удовлетворяет при $u \in \Omega$ условию $\left|J(u)\right| \geq \alpha > 0$.

Тогда, если $G$ — измеримый компакт с кусочно-гладкой границей, лежащий во множестве $\Omega$ и $f(x)$ — непрерывна на множестве $G’ = F(G)$, то справедлива следующая формула замены переменных в кратном интеграле:
$$\int\limits_{G’} f(x)\,dx = \int\limits_G f(\phi_1(u), \ldots, \phi_n(u))\left|J(u)\right|\,du\quad(*),$$
где $x = (x_1, \ldots, x_n),\quad u = (u_1, \ldots, u_n)$.

Доказательство

Для начала рассмотрим еще 2 вспомогательных свойства:

  1. Если $L \subset \Omega$ есть непрерывно дифференцируемая кривая, то ее образ $L’ = F(L)$ есть непрерывно дифференцируемая кривая.
  2. Если $G$ — область и $\overline{G} \subset \Omega$ (где $\overline{G}$ — замыкание области $G$), тогда ее образ $G’ = F(G)$. Образ границы $\Omega$ есть граница $\Omega’$.

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

Рассмотрим доказательство для плоского случая (двойных интегралов). В силу свойств непрерывных функций образ $G’$ компакта $G$ при непрерывном и взаимно однозначном отображении $F$ является компактом, а по свойствам отображения $F$, указанным выше, граница компакта $G’$ является кусочно-гладкой кривой. Кусочно-гладкая кривая имеет жорданову меру нуль, а так как ограниченное множество измеримо по Жордану тогда и только тогда, когда его граница имеет жорданову меру нуль, то компакт $G’$ измерим, а оба интеграла в формуле $(*)$ существуют как интегралы от функций, непрерывных на компактах.

Поскольку компакт $G$ лежит в открытом множестве $\Omega$, то границы этих множеств не пересекаются. Так как граница любого множества замкнута и граница ограниченного множества ограничена, то расстояние между границами множеств $G$ и $\Omega$ есть положительное число $\delta$.

Примечание №1: под разбиением множества $A$ далее будем подразумевать совокупность измеримых множеств $\{A_1, \ldots, A_n\}$, таких что $A_1\cup \ldots \cup A_n = A$ и $A_i \cap A_j = \oslash, i \ne j$. Клеткой назовем множество вида $K = \{(x_1, \ldots, x_n)| a_i \leq x_i < b_i, 1 \leq i \leq n\}$, прямоугольником — клетку в пространстве $\mathbb{R}^2$.

Пусть $P$ есть замкнутый квадрат, содержащий компакт $G$. Если разбить стороны квадрата $P$ на равные части длины $h < \delta$ (чтобы отсутствовали квадраты, содержащие одновременно элементы границ $G$ и $\Omega$), то и сам квадрат $P$ окажется разбит на квадратные клетки с площадью $h^2$. Разбиение квадрата $P$ порождает разбиение $T$ компакта $G.$ Если малый квадрат со стороной $h$ целиком лежит внутри компакта $G$, то он является элементом разбиения $T$, а если он содержит граничные точки $G$, то соответствующим элементом разбиения является пересечение этого квадрата с компактом $G.$ Отображение $F$ порождает разбиение $T’$ компакта $G’ = F(G)$, причем элементами разбиения $T’$ являются образы элементов разбиения $T$. При отбрасывании в интегральной сумме слагаемых, которым отвечают квадраты, имеющие непустое пересечение с множеством жордановой меры нуль, характер соответствующего предела при мелкости разбиения, стремящемся к нулю, не изменится (о чем свидетельствует соответствующая лемма, см. примечание №2). А значит, при написании интегральных сумм можно учитывать только слагаемые, соответствующие целым квадратам и их образам при отображении $F$, остальные квадраты будут иметь непустое пересечение с границей $G$. Так как отображение $F$ равномерно непрерывно, то мелкость разбиения $T’$ стремится к нулю, когда стремится к нулю мелкость разбиения $T$.

Если малые квадраты $P_1, \ldots, P_n$ лежат внутри компакта $G$, то
их образы $P_1′, \ldots, P_n’$ лежат внутри $G’$. Пусть $(u_i, v_i)$ — координаты точки, лежащей в левом нижнем углу квадрата $P_i$, a $(\phi(u_i, v_i),\psi(u_i, v_i))$ — образ этой точки при отображении $F$.

Тогда можем записать интегралы, входящие в формулу $(*)$ как пределы интегральных сумм:
$\iint\limits_{G’}f(x, y)\,dxdy = \lim\limits_{h \to 0} \sum\limits_{i=1}\limits^n{f(x_i, y_i,)m(P_i’)},$
$\iint\limits_G f(\phi(u, v), \psi(u, v))\left|J(u, v)\right|\,dudv = $ $\lim\limits_{h \to 0} \sum\limits_{i=1}\limits^n{f(\phi(u_i, v_i), \psi(u_i, v_i))\left|J(u_i, v_i)\right|m(P_i)}.$

Для доказательства формулы $(*)$ покажем, что разность этих интегральных сумм стремится к нулю при $h \to 0$. В силу леммы о геометрическом смысле модуля якобиана отображения,
$\left|m(P_i’) — \left|J(u_i, v_i)\right|m(P_i)\right| \leq \alpha(h)m(P_i), \lim\limits_{h \to 0} \alpha(h) = 0$.
Принимая во внимание, что $\phi(u_i, v_i) = x_i, \psi(u_i, v_i) = y_i, \left|f(x, y)\right| < M$ (последнее в силу того, что функция [latex]f[/latex] непрерывна на компакте, а значит и ограниченна на нем), получаем оценку для разности интегральных сумм:
$\left|\sum\limits_{i=1}\limits^n{f(x_i, y_i)m(P_i’)} — \sum\limits_{i=1}\limits^n{f(\phi(u_i, v_i), \psi(u_i, v_i))\left|J(u_i, v_i)\right|m(P_i)}\right| \leq $ $\sum\limits_{i=1}\limits^n{\left|f(x_i, y_i)m(P_i’) — f(x_i, y_i)\left|J(u_i, v_i)\right|m(P_i)\right|} = $ $\sum\limits_{i=1}\limits^n{\left|f(x_i, y_i)\left|m(P_i’) — \left|J(u_i, v_i)\right|m(P_i)\right|\right|} \leq $ $M\sum\limits_{i=1}\limits^n{\alpha(h)m(P_i)} = $ $M\alpha(h)\sum\limits_{i=1}\limits^n{m(P_i)} \leq $ $M\alpha(h)m(G)$, из которой следует, что эта разность стремится к нулю при $h \to 0$ (т.к. $M$ и $m(G)$ — константы). Теорема доказана.

Примечание №2: о геометрическом смысле модуля якобиана отображения можно прочитать, например, в курсе лекций по мат. анализу В.И. Коляда, А.А. Кореновский (т.2, стр. 219) или в Тер-Крикоров А.М. и Шабунин М.И. «Курс математического анализа» (стр. 471). Лемма об отбрасывании слагаемых в интегральной сумме также присутствует и доказана, например, в учебнике Тер-Крикорова, стр. 458.

Замечание

Нарушение условия взаимной однозначности на множестве меры нуль и обращение якобиана отображения в
нуль на множестве меры нуль не влияют на справедливость формулы $(*)$ замены переменных в кратном интеграле. Такое множество $E$ меры нуль всегда можно накрыть клеточным множеством $A \subset G$ сколь угодно малой меры $\epsilon$, разбивающимся на квадраты. Из доказательства теоремы следует, что при отображении $F : G \to \mathbb{R}^n$ мера множества $A$ возрастет не более чем в [latex]c[/latex] раз, где $c$ — фиксированное. Поэтому найдутся такие константы $c_1$ и $c_2$, что $\left|\int\limits_{A’}f(x)\,dx\right| < c_1\epsilon, A’ = F(A),$ $\left|\int\limits_A f(\phi_1(u), \ldots, \phi_n(u))\left|J(u)\right|\,du\right| < c_2\epsilon, \forall \epsilon > 0.$ На множестве $G \setminus A$ условия теоремы выполнены, а так как интегралы $\int\limits_{G’}f(x)\,dx$ и $\int\limits_G f(\phi_1(u), \ldots, \phi_n(u))\left|J(u)\right|\,du$ отличаются на сколь угодно малое число, что следует из соответствующих неравенств, то они совпадают. Данное замечание является важным, в частности, при обосновании перехода к полярным, цилиндрическим и сферическим системам координат.

[свернуть]

Примеры

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

Тест: формула замены переменной в кратном интеграле

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


Таблица лучших: Замена переменной в кратных интегралах

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

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

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

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

Из курса аналитической геометрии известны следующие соотношения между декартовыми и полярными координатами: $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 баллов
Место Имя Записано Баллы Результат
Таблица загружается
Нет данных

Подстановки степени n

 Определение:

Любое взаимно однозначное отображение $latex A$ множества первых $latex n$ натуральных чисел на себя называется подстановкой $latex n$-й степени. Всякая подстановка $latex A$ может быть записана при помощи двух перестановок, подписанных одна под другой:

$latex
\begin{pmatrix}
i_{1} & i_{2} & \ldots & i_{n}\\
a_{i_{1}} & a_{i_{2}} & \ldots & a_{i_{n}}
\end{pmatrix} $
через $latex a_{i}$ здесь обозначается то число, в которое при подстановке $latex A$ переходит число $latex i$, $latex i = 1,2,$ $latex \ldots , n$.

Замечание:

От одной записи подстановки $latex A$ к другой можно перейти при помощи транспозиций столбиков. Любая подстановка $latex n$-й степени может быть записана в виде:

$latex \begin{pmatrix}
1 & 2 & \ldots & n\\
a_{1} & a_{2} & \ldots & a_{n}
\end{pmatrix} $
Т.е. с натуральным расположением чисел в верхней строке. При такой записи различные подстановки различаются друг от друга перестановками, стоящими в нижней строке, и поэтому число перестановок из $latex n$ чисел равно $latex n!$ .

Переход к любой другой записи подстановки $latex A$ можно осуществить, как мы знаем, путём последовательного выполнения нескольких транспозиций в верхней строке и соответствующих им транспозиций в нижней строке. Мы одновременно меняем чётность и поэтому сохраняем совпадение или противоположность этих чётностей.

Отсюда следует, что либо при всех записях подстановки $latex A$ чётности верхней и нижней строк совпадают, либо же при всех записях они противоположны. В первом случае подстановка $latex A$ называется чётной, а во втором — нечётной. В частности тождественная подстановка($latex E$) будет чётной:

$latex E =
\begin{pmatrix}
1 & 2 & \ldots & n\\
1 & 2 & \ldots & n
\end{pmatrix} $Число чётных подстановок равно числу нечётных, равно $latex {\frac12 n!}$

Пример

$latex \begin{pmatrix}
4 & 3 & 5 & 2 &1\\
3 & 5 & 2 & 1 &4
\end{pmatrix}$ всегда можно представить в виде $latex \begin{pmatrix}
1 & 2 & 3 & 4 &5\\
4 & 1 & 5 & 3 &2
\end{pmatrix}$

Подстановки степени n

Тест по теме «Подстановки степени $latex n$».

Таблица лучших: Подстановки степени n

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

Список литературы: