Тройные интегралы 🌐

Двойной интеграл научил тебя считать объём тела, накрытого поверхностью $z=f(x,y)$, — по сути, ты уже интегрировал по трёхмерному пространству, просто одна из трёх координат была не независимой переменной, а значением функции. Сегодня мы делаем следующий логичный шаг: убираем это ограничение полностью. Область интегрирования теперь — не плоская фигура, а настоящее тело в пространстве, произвольная трёхмерная область $V$, и мы интегрируем по ней функцию трёх независимых переменных $f(x,y,z)$. Формально это звучит как небольшое усложнение, но на практике тройной интеграл открывает доступ к вычислению объёма и массы объектов совершенно любой формы — от простого куба до неоднородной планеты с переменной плотностью породы.

Давай сразу обозначим, зачем тебе это в контексте машинного обучения. Тройные интегралы возникают всякий раз, когда нужно нормировать плотность вероятности, заданную в трёхмерном пространстве, — то есть подобрать множитель, при котором интеграл плотности по всему пространству (или по какой-то области) равен единице. Это не абстрактная задача: в физически обоснованных генеративных моделях (например, там, где данные — это трёхмерные облака точек, координаты атомов молекулы или объекты в трёхмерной сцене) плотность вероятности часто задаётся ненормированной «энергетической» формулой, и чтобы превратить её в настоящую плотность, нужно вычислить нормировочную константу — то есть взять тройной интеграл. А ещё сегодня мы явно свяжем переход к сферическим координатам с той же самой идеей, что используется при работе с многомерными гауссовскими распределениями в статистике: там тоже переходят к более удобной системе координат (обычно связанной с радиусом и углами), чтобы разложить сложный многомерный интеграл на произведение простых одномерных.

Представь себе, что тебе нужно взвесить не плоскую металлическую пластину (это была задача прошлого урока), а настоящее трёхмерное тело — скажем, каменный метеорит неправильной формы, плотность которого меняется от рыхлой оболочки к плотному металлическому ядру. Обычные формулы «масса равна плотности, умноженной на объём» здесь не работают, потому что плотность разная в разных точках. Единственный честный способ — разбить тело на бесконечно много бесконечно малых кубиков, в каждом из которых плотность можно считать практически постоянной, вычислить массу каждого кубика (плотность, умноженная на его крошечный объём) и сложить все эти вклады. Именно эта идея — предельный переход от суммы бесконечно малых вкладов к интегралу — и есть тройной интеграл.

В этом уроке мы пройдём полный путь: сначала дадим строгое определение тройного интеграла и разберём его геометрический смысл (объём тела) и физический смысл (масса тела через плотность). Затем научимся сводить тройной интеграл к трём последовательным однократным интегрированиям — обобщение теоремы Фубини на трёхмерный случай. После этого перейдём к двум важнейшим заменам переменных: цилиндрическим координатам (для тел с осевой симметрией — цилиндров, конусов) и сферическим координатам (для тел со сферической симметрией — шаров, слоёв), причём для обеих систем выведем якобиан не по памяти, а строгим вычислением определителя матрицы Якоби. Наконец, разберём, как именно сферические координаты помогают элегантно вычислить нормировочную константу трёхмерного гауссовского распределения — прямой мостик к тому, что тебе понадобится в статистике и машинном обучении.

История: откуда это взялось?

Идея разбиения тела на маленькие элементы для вычисления объёма и массы восходит ещё к методу Кавальери и работам Ньютона и Лейбница по интегральному исчислению XVII века, но полноценный, строго обоснованный аппарат кратных интегралов — двойных и тройных — оформился заметно позже, в XVIII–XIX веках, в первую очередь благодаря нуждам механики и физики. Задачи о массе неоднородных тел, о центре масс, о моменте инерции вращающегося твёрдого тела возникали у механиков постоянно: чтобы правильно рассчитать движение волчка, снаряда или планеты, нужно было учитывать не просто суммарную массу тела, а распределение этой массы по объёму. Леонард Эйлер и Жозеф Луи Лагранж, разрабатывая механику твёрдого тела в XVIII веке, регулярно сталкивались с интегралами по трёхмерным областям — и хотя формального названия «тройной интеграл» тогда ещё не было, вычисления по сути велись именно так: суммированием бесконечно малых вкладов массы по всему объёму тела.

Строгое обоснование кратных интегралов как предела интегральных сумм и, что важно для сегодняшнего урока, теорема о сведении кратного интеграла к повторному однократному интегрированию (обобщение того, что для двойных интегралов носит имя Гвидо Фубини, хотя сама техника использовалась математиками задолго до формальных доказательств Фубини в начале XX века) сформировались в рамках развития математического анализа XIX века, в частности в трудах Огюстена Коши (того самого автора градиентного спуска из прошлых уроков) и позже Бернхарда Римана, заложивших строгую теорию интеграла. Замена переменных в кратных интегралах и формула с якобианом — естественное развитие той же идеи, что мы уже видели для одной и двух переменных: при переходе к новым координатам объём (или площадь) элементарной ячейки растягивается или сжимается неравномерно, и якобиан — это в точности локальный коэффициент этого искажения.

Переход к сферическим координатам оказался особенно важен не только в геометрии, но и в физике: закон всемирного тяготения Ньютона, электростатическое поле точечного заряда, распределение температуры внутри сферически симметричного тела — всё это естественно описывается на языке радиуса $r$ и двух углов, а не на языке декартовых координат $x,y,z$. Собственно поэтому вывод якобиана $r^2\sin\varphi$ для сферических координат — не просто учебное упражнение, а инструмент, без которого невозможно было бы аналитически решить огромный класс задач классической физики XVIII–XIX веков — от расчёта гравитационного притяжения однородного шара (Ньютон решал эту задачу, по сути, эквивалентными геометрическими методами задолго до появления современной записи через тройной интеграл) до вычисления электрического поля заряженной сферы. Сегодня та же самая математика — переход к координатам, согласованным с симметрией задачи, — лежит в основе того, как статистики и специалисты по машинному обучению нормируют многомерные вероятностные распределения, обладающие сферической симметрией.

Определение и смысл: объём и масса тела в пространстве

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

Давай разберёмся с определением так же, как мы разбирались с двойным интегралом, только добавив ещё одно измерение. Возьми тело $V$ в пространстве — произвольную ограниченную трёхмерную область. Разбей его на множество маленьких прямоугольных параллелепипедов (кубиков) со сторонами $\Delta x, \Delta y, \Delta z$. В каждом таком кубике выбери произвольную точку $(x_i,y_j,z_k)$ и вычисли значение функции $f(x_i,y_j,z_k)$ в этой точке. Умножь это значение на объём кубика $\Delta V=\Delta x\,\Delta y\,\Delta z$ и сложи такие произведения по всем кубикам, покрывающим тело $V$. Устремляя размеры кубиков к нулю (а их количество — к бесконечности), ты получаешь тройной интеграл.

Если функция $f(x,y,z)\equiv1$ (то есть в каждой точке она равна единице), сумма вкладов $f\cdot\Delta V$ — это просто сумма объёмов всех кубиков, и в пределе она даёт в точности объём тела $V$. Если же $f(x,y,z)=\rho(x,y,z)$ — плотность вещества в точке $(x,y,z)$ (масса на единицу объёма), то вклад каждого кубика $\rho\cdot\Delta V$ — это приближённая масса этого маленького кусочка тела, а сумма всех вкладов в пределе даёт полную массу тела. Это ровно та же логика, что для двойного интеграла и массы плоской пластины, просто одним измерением больше.

Определение

Определение (тройной интеграл). Пусть функция $f(x,y,z)$ определена и ограничена на ограниченной области $V\subset\mathbb{R}^3$. Разобьём $V$ на $n$ маленьких элементов объёма $\Delta V_1,\dots,\Delta V_n$, в каждом выберем произвольную точку $(x_i,y_i,z_i)$. Тройным интегралом функции $f$ по области $V$ называется предел интегральной суммы при стремлении размеров всех элементов разбиения к нулю:

$$\iiint_V f(x,y,z)\,dV = \lim_{n\to\infty}\sum_{i=1}^n f(x_i,y_i,z_i)\,\Delta V_i$$

если этот предел существует и не зависит от способа разбиения и выбора точек. Здесь $dV$ — элемент объёма (в декартовых координатах $dV=dx\,dy\,dz$).

Геометрический смысл: если $f(x,y,z)\equiv1$, тройной интеграл равен объёму тела:

$$V_{\text{объём}} = \iiint_V dV$$

Физический смысл: если $f(x,y,z)=\rho(x,y,z)$ — плотность вещества в теле $V$, тройной интеграл равен массе тела:

$$m = \iiint_V \rho(x,y,z)\,dV$$

Обрати внимание на прямую аналогию с прошлым уроком: там двойной интеграл по плоской области давал площадь (при $f\equiv1$) или массу пластины (при $f=$ плотность); здесь тройной интеграл по трёхмерной области даёт объём (при $f\equiv1$) или массу тела (при $f=$ плотность). Эта закономерность продолжится и дальше в курсе: кратный интеграл любой размерности от единичной функции всегда даёт «объём» области интегрирования в соответствующей размерности, а от функции плотности — «массу», распределённую по этой области.

Разбор примеров

Пример 1 (лёгкий). Вычислить тройной интеграл $\iiint_V (x+y+z)\,dV$, где $V=[0,1]\times[0,1]\times[0,1]$ — единичный куб.

Разложим интеграл на сумму трёх слагаемых: $\iiint_V x\,dV+\iiint_V y\,dV+\iiint_V z\,dV$. По симметрии куба (все три переменные меняются в одинаковых пределах) все три интеграла равны между собой. Посчитаем первый:

Ответ: сумма трёх одинаковых слагаемых даёт $\iiint_V(x+y+z)\,dV=\dfrac12+\dfrac12+\dfrac12=\dfrac32$.

Пример 2 (средний). Найти объём тетраэдра, ограниченного координатными плоскостями $x=0$, $y=0$, $z=0$ и плоскостью $x+y+z=1$.

Область $V$: $x$ меняется от $0$ до $1$, при фиксированном $x$ переменная $y$ меняется от $0$ до $1-x$, а при фиксированных $x,y$ переменная $z$ меняется от $0$ до $1-x-y$ (это верхняя граница, заданная плоскостью). Объём — тройной интеграл единичной функции:

Внутренний интеграл по $y$:

Внешний интеграл по $x$:

Ответ: $V=\dfrac16$ — тот же результат, что даёт известная школьная формула объёма тетраэдра с тремя взаимно перпендикулярными рёбрами единичной длины ($V=\frac{abc}{6}$ при $a=b=c=1$), только полученный честным тройным интегрированием, без готовой формулы.

Пример 3 (сложный, машинное обучение). В энергетических моделях ненормированная плотность вероятности часто задаётся как $\tilde p(x,y,z)=e^{-E(x,y,z)}$, где $E$ — «энергия» состояния. Пусть на кубе $[0,1]^3$ энергия равна $E(x,y,z)=x+y+z$, то есть $\tilde p(x,y,z)=e^{-(x+y+z)}$. Найти нормировочную константу $Z=\iiint_{[0,1]^3}\tilde p\,dV$, чтобы $p_{\text{норм}}=\tilde p/Z$ была настоящей плотностью вероятности (интеграл по кубу равен единице).

Экспонента суммы раскладывается в произведение экспонент, а область интегрирования — это прямое произведение трёх одинаковых отрезков, поэтому тройной интеграл распадается на произведение трёх одинаковых однократных интегралов:

Численно $1-e^{-1}\approx1-0{,}3679=0{,}6321$, откуда $Z\approx0{,}6321^3\approx0{,}2526$.

Ответ: $Z=(1-e^{-1})^3\approx0{,}2526$. Именно на такую константу нужно поделить исходную ненормированную плотность $e^{-(x+y+z)}$, чтобы получить корректное распределение вероятностей на кубе. Это ровно та задача, с которой сталкивается специалист по машинному обучению, работая с энергетическими моделями или ненормированными плотностями: сама «форма» распределения (энергия $E$) может быть простой и хорошо интерпретируемой, а вот превращение её в настоящую вероятность требует вычисления интеграла нормировки — часто трёхмерного или ещё более многомерного.

Почему это важно

Тройной интеграл — это прямое, честное обобщение того, что ты уже умеешь делать с двойным интегралом, только теперь ты можешь работать с настоящими трёхмерными телами, а не только с поверхностями над плоскими областями. Геометрический смысл (объём) и физический смысл (масса через плотность) — это два языка описания одного и того же математического объекта, и очень часто задача, сформулированная как «найти объём», и задача «найти массу» решаются буквально одинаковыми выкладками с той лишь разницей, стоит ли под интегралом множитель плотности. А пример с нормировочной константой показывает: тот же самый инструмент, которым ты считаешь объём геометрического тела, работает и для нормировки вероятностных распределений — это одна из причин, почему тройные интегралы регулярно всплывают в теории вероятностей и машинном обучении, а не только в геометрии и механике.

Сведение к повторному интегралу: три интегрирования подряд

Интуиция: слой за слоем, как в двойном интеграле, только на уровень глубже

Вычислять тройной интеграл напрямую по определению (через предел интегральных сумм) практически никогда не приходится — точно так же, как для двойного интеграла, есть теорема, сводящая тройной интеграл к трём последовательным однократным интегрированиям. Идея прямое обобщение теоремы Фубини для двойных интегралов: если тело $V$ можно описать как множество точек $(x,y)$ из некоторой плоской области $D$ (это «тень» тела на плоскость $xy$), причём для каждой такой точки $(x,y)$ координата $z$ пробегает от нижней границы $z_1(x,y)$ до верхней границы $z_2(x,y)$, то сначала интегрируем по $z$ (получаем функцию от $x,y$ — по сути, «толщину» тела в этой точке, умноженную на значение функции), а затем берём получившийся двойной интеграл по плоской области $D$ уже знакомым способом из прошлого урока.

Определение

Теорема (сведение тройного интеграла к повторному). Пусть тело $V=\{(x,y,z):\ (x,y)\in D,\ z_1(x,y)\le z\le z_2(x,y)\}$, где $D$ — проекция тела $V$ на плоскость $xy$. Тогда

$$\iiint_V f(x,y,z)\,dV = \iint_D\left[\int_{z_1(x,y)}^{z_2(x,y)} f(x,y,z)\,dz\right]dA$$

Если далее область $D=\{(x,y):\ a\le x\le b,\ y_1(x)\le y\le y_2(x)\}$, двойной интеграл справа тоже сводится к повторному однократному интегрированию (теорема Фубини для двойных интегралов из прошлого урока), и в итоге получается три последовательных однократных интегрирования:

$$\iiint_V f(x,y,z)\,dV = \int_a^b\int_{y_1(x)}^{y_2(x)}\int_{z_1(x,y)}^{z_2(x,y)} f(x,y,z)\,dz\,dy\,dx$$

Важная деталь: порядок интегрирования по трём переменным можно менять (получится шесть возможных вариантов порядка — $dz\,dy\,dx$, $dz\,dx\,dy$, $dy\,dz\,dx$ и так далее), выбирая тот, для которого пределы интегрирования и сама подынтегральная функция описываются проще всего. Это ровно та же свобода выбора, что была для двойного интеграла, только теперь с тремя переменными вариантов заметно больше.

Разбор примеров

Пример 1 (лёгкий). Найти объём прямоугольного параллелепипеда $[0,2]\times[0,3]\times[0,4]$ через тройной интеграл единичной функции.

Все пределы — константы, поэтому интеграл сразу распадается на произведение трёх длин рёбер:

Ответ: $V=24$ — совпадает с элементарной формулой «длина умножить на ширину умножить на высоту», только теперь ты видишь, что эта формула — частный случай тройного интеграла по прямоугольной области.

Пример 2 (средний). Найти объём тела, ограниченного цилиндром $x^2+y^2\le1$ снизу плоскостью $z=0$ и сверху параболоидом $z=x^2+y^2$.

Здесь удобнее описать проекцию $D$ на плоскость $xy$ как единичный круг $x^2+y^2\le1$, а высота тела в каждой точке — от $z=0$ до $z=x^2+y^2$. В декартовых координатах:

Этот двойной интеграл проще всего взять в полярных координатах (приём из прошлого урока): $x^2+y^2=r^2$, $dA=r\,dr\,d\theta$, область $D$ становится прямоугольником $0\le r\le1$, $0\le\theta\le2\pi$:

Ответ: $V=\dfrac{\pi}2\approx1{,}571$. Заметь: даже когда мы формально ещё не вводили полноценные цилиндрические координаты для тройного интеграла, идея свести двойной интеграл по круглой проекции к полярным координатам уже работает — это прямой мостик к следующему разделу урока.

Пример 3 (сложный, машинное обучение). Известно, что стандартное трёхмерное нормальное распределение с независимыми координатами имеет ненормированную плотность $\tilde p(x,y,z)=e^{-(x^2+y^2+z^2)/2}$. Не переходя пока к специальным координатам (это мы сделаем в разделе про сферические координаты), покажи, что тройной интеграл этой функции по всему пространству $\mathbb{R}^3$ распадается на произведение трёх одинаковых одномерных гауссовых интегралов, и найди значение нормировочной константы, зная классический результат $\int_{-\infty}^{\infty} e^{-t^2/2}\,dt=\sqrt{2\pi}$.

Показатель экспоненты — это сумма трёх независимых слагаемых $-x^2/2-y^2/2-z^2/2$, поэтому сама экспонента раскладывается в произведение трёх множителей, каждый из которых зависит только от одной переменной:

Область интегрирования — всё пространство $\mathbb{R}^3=(-\infty,\infty)^3$, то есть тоже «прямое произведение» трёх одинаковых бесконечных отрезков. По той же логике, что и в примере 3 предыдущего раздела, тройной интеграл распадается в произведение трёх однократных интегралов:

Ответ: $Z=(2\pi)^{3/2}\approx15{,}749$. Это ровно та формула нормировочной константы стандартного трёхмерного гауссовского распределения, которую ты, возможно, уже видел в курсах статистики или машинного обучения в виде общей формулы $(2\pi)^{n/2}$ для размерности $n$ (при $n=1$ это классический $\sqrt{2\pi}$, при $n=2$ — это $2\pi$, знакомый из прошлого урока про двойные интегралы). Мы ещё вернёмся к этому же интегралу чуть позже, но уже через сферические координаты, — и убедимся, что оба способа дают одинаковый ответ.

Почему это важно

Сведение тройного интеграла к трём последовательным однократным интегрированиям — это единственный практический способ вычислять такие интегралы вручную (или программно, шаг за шагом). Ключевой навык здесь — правильно описать область интегрирования: понять, как устроена «тень» тела на плоскость и как меняются пределы внутренней переменной в зависимости от значений внешних. А третий пример показывает важный практический приём, который то и дело будет тебе полезен в статистике и машинном обучении: если подынтегральная функция и область интегрирования факторизуются (раскладываются на независимые множители по каждой переменной), тройной интеграл автоматически превращается в произведение трёх простых однократных интегралов — и тебе не придётся вычислять по-настоящему сложный трёхмерный интеграл вообще.

Цилиндрические координаты: вывод якобиана и применение

Интуиция: полярные координаты плюс обычная высота

Давай разберёмся, что делать, если тело обладает не прямоугольной, а осевой (вращательной) симметрией — например, это цилиндр, конус или тело вращения вокруг оси $z$. В таких случаях декартовы координаты $x,y,z$ неудобны: границы области получаются громоздкими (вспомни, как в прошлом уроке круглая область была неудобна в декартовых координатах, но становилась элементарной в полярных). Естественное решение — взять полярные координаты $(r,\theta)$ в горизонтальной плоскости $xy$ (то есть описывать положение точки не через $x,y$, а через расстояние от оси $z$ и угол поворота вокруг неё) и оставить вертикальную координату $z$ без изменений. Это и есть цилиндрические координаты — буквально полярные координаты, «поднятые» в третье измерение.

Определение и вывод якобиана

Определение (цилиндрические координаты). Положение точки задаётся тройкой $(r,\theta,z)$, где $r\ge0$ — расстояние от точки до оси $z$, $\theta\in[0,2\pi)$ — угол поворота в плоскости $xy$ (как в полярных координатах), а $z$ — обычная декартова координата высоты. Связь с декартовыми координатами:

$$x=r\cos\theta,\qquad y=r\sin\theta,\qquad z=z$$

Чтобы найти элемент объёма $dV$ в новых координатах, нужен якобиан — определитель матрицы частных производных старых координат по новым. Составим матрицу Якоби:

Третья строка почти пустая (только один ненулевой элемент), поэтому удобно разложить определитель именно по ней: определитель равен произведению единственного ненулевого элемента $(1)$ на его алгебраическое дополнение — определитель верхнего левого блока $2\times2$:

Значит $\det J=r$, и элемент объёма в цилиндрических координатах равен

Это ровно тот же множитель $r$, что и в полярных координатах для двойного интеграла из прошлого урока (недаром — цилиндрические координаты буквально «умножают» полярное разбиение плоскости $xy$ на дополнительный слой по оси $z$, который сам по себе якобиан не искажает).

Формула замены переменных. Если тело $V$ в цилиндрических координатах описывается неравенствами $r_1\le r\le r_2$, $\theta_1\le\theta\le\theta_2$, $z_1(r,\theta)\le z\le z_2(r,\theta)$, то

$$\iiint_V f(x,y,z)\,dV = \int_{\theta_1}^{\theta_2}\int_{r_1}^{r_2}\int_{z_1}^{z_2} f(r\cos\theta,\,r\sin\theta,\,z)\,r\,dz\,dr\,d\theta$$

Разбор примеров

Пример 1 (лёгкий). Найти объём цилиндра радиуса $R$ и высоты $h$ через тройной интеграл в цилиндрических координатах.

Область: $0\le r\le R$, $0\le\theta\le2\pi$, $0\le z\le h$ — все пределы постоянны.

Ответ: $V=\pi R^2 h$ — знакомая школьная формула объёма цилиндра, но теперь получена честным тройным интегрированием с явным учётом якобиана $r$.

Пример 2 (средний). Найти массу цилиндра радиуса $R=1$ и высоты $h=2$ (ось цилиндра совпадает с осью $z$, основание при $z=0$) с плотностью $\rho(r,z)=rz$ (плотность растёт и с удалением от оси, и с высотой).

Считаем множители по отдельности: $\displaystyle\int_0^{2\pi}d\theta=2\pi$, $\displaystyle\int_0^1 r^2\,dr=\frac13$, $\displaystyle\int_0^2 z\,dz=2$.

Ответ: $m=\dfrac{4\pi}3\approx4{,}189$.

Пример 3 (сложный, физика с прямым выходом на ML-интуицию). Найти момент инерции однородного сплошного цилиндра (плотность $\rho_0=\text{const}$, радиус $R$, высота $h$) относительно его оси симметрии (оси $z$), и выразить результат через массу цилиндра.

Момент инерции относительно оси $z$ — это интеграл от плотности, умноженной на квадрат расстояния точки до оси, а расстояние до оси в цилиндрических координатах — это в точности $r$:

Масса всего цилиндра $m=\rho_0\pi R^2h$ (пример 1 этого раздела, домноженный на постоянную плотность). Выразим момент инерции через массу:

Ответ: $I_z=\dfrac12 mR^2$ — это в точности классическая формула момента инерции сплошного однородного цилиндра относительно своей оси. Обрати внимание на приём, который здесь сработал: осевая симметрия задачи (и сама физическая величина, и форма тела зависят только от расстояния до оси, не от угла) сделала цилиндрические координаты идеальным выбором — интеграл по углу $\theta$ тривиально дал множитель $2\pi$, а вся содержательная работа свелась к одномерным интегралам по $r$ и по $z$. Ровно тот же принцип — использовать координаты, согласованные с симметрией задачи, — применяется и при работе с данными, которые естественно описываются в цилиндрической или сферической геометрии (например, облака точек с осевой или сферической симметрией в трёхмерных генеративных моделях).

Почему это важно

Цилиндрические координаты — это ровно тот инструмент, который превращает громоздкую задачу об интегрировании по цилиндру, конусу или телу вращения в набор простых, разделяющихся вычислений. Ключевое наблюдение, которое стоит унести с собой: если задача обладает симметрией относительно оси (то есть ничего не зависит от угла поворота вокруг этой оси), интеграл по углу $\theta$ почти всегда даёт тривиальный множитель $2\pi$ (или его часть, если рассматривается не полный оборот), и вся содержательная работа сосредотачивается в оставшихся двух переменных. Не забывай про якобиан $r$ — без него результат будет систематически неверным, причём именно в тех областях, что дальше от оси, ошибка будет больше.

Сферические координаты: вывод якобиана и применение

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

Для тел со сферической симметрией — шаров, сферических слоёв, конусов с вершиной в начале координат — идеальны сферические координаты. Идея похожа на то, как задаётся положение точки на глобусе: вместо трёх декартовых чисел $x,y,z$ используем расстояние от центра $r$ (насколько далеко точка от начала координат) и два угла — один отсчитывается от полярной оси $z$ (аналог широты, только отсчитываемый не от экватора, а от полюса), второй — это тот же азимутальный угол $\theta$, что и в цилиндрических координатах (аналог долготы).

Определение и вывод якобиана

Определение (сферические координаты). Положение точки задаётся тройкой $(r,\varphi,\theta)$, где $r\ge0$ — расстояние от точки до начала координат, $\varphi\in[0,\pi]$ — полярный угол, отсчитываемый от положительного направления оси $z$, $\theta\in[0,2\pi)$ — азимутальный угол в плоскости $xy$. Связь с декартовыми координатами:

$$x=r\sin\varphi\cos\theta,\qquad y=r\sin\varphi\sin\theta,\qquad z=r\cos\varphi$$

Составим матрицу Якоби, дифференцируя каждую из трёх формул по каждой из трёх новых переменных $r,\varphi,\theta$:

Третья строка удобна для разложения определителя: в ней есть ноль. Раскладываем по третьей строке, с учётом знаков алгебраических дополнений $(-1)^{3+j}$:

Минор $M_{31}$ (вычёркиваем третью строку и первый столбец, берём определитель из строк 1–2 и столбцов 2–3):

Минор $M_{32}$ (вычёркиваем третью строку и второй столбец, берём определитель из строк 1–2 и столбцов 1, 3):

Подставляем оба минора обратно:

Поскольку $\varphi\in[0,\pi]$, синус этого угла неотрицателен, поэтому модуль якобиана равен самому выражению без изменения знака:

Формула замены переменных. Если тело $V$ в сферических координатах описывается неравенствами $r_1\le r\le r_2$, $\varphi_1\le\varphi\le\varphi_2$, $\theta_1\le\theta\le\theta_2$, то

$$\iiint_V f(x,y,z)\,dV = \int_{\theta_1}^{\theta_2}\int_{\varphi_1}^{\varphi_2}\int_{r_1}^{r_2} f(r\sin\varphi\cos\theta,\,r\sin\varphi\sin\theta,\,r\cos\varphi)\;r^2\sin\varphi\;dr\,d\varphi\,d\theta$$

Разбор примеров

Пример 1 (лёгкий). Найти объём шара радиуса $R$ через тройной интеграл в сферических координатах.

Считаем множители: $\displaystyle\int_0^{2\pi}d\theta=2\pi$, $\displaystyle\int_0^{\pi}\sin\varphi\,d\varphi=\bigl[-\cos\varphi\bigr]_0^\pi=1-(-1)=2$, $\displaystyle\int_0^R r^2\,dr=\frac{R^3}3$.

Ответ: $V=\dfrac{4\pi R^3}3$ — знакомая школьная формула объёма шара, теперь строго выведенная через тройной интеграл с явным якобианом $r^2\sin\varphi$.

Пример 2 (средний). Найти массу шара радиуса $R=2$ с плотностью, растущей линейно с расстоянием от центра: $\rho(r)=r$.

Ответ: $m=16\pi\approx50{,}265$.

Пример 3 (сложный, машинное обучение — многомерный гауссиан). Вернёмся к трёхмерному стандартному гауссовскому распределению из примера 3 предыдущего раздела: ненормированная плотность $\tilde p(x,y,z)=e^{-(x^2+y^2+z^2)/2}$. На этот раз вычислим нормировочную константу $Z=\iiint_{\mathbb{R}^3}\tilde p\,dV$ не факторизацией на три одномерных интеграла, а напрямую через переход к сферическим координатам — точно так же, как в статистике многомерных гауссовских распределений часто переходят к радиусу и углам, чтобы отделить «радиальную» часть интеграла от «угловой».

В сферических координатах $x^2+y^2+z^2=r^2$, поэтому:

Угловые интегралы дают в точности $2\pi$ и $2$ — заметь, что вся угловая часть в сумме даёт множитель $4\pi$, который в физике и статистике часто называют «полным телесным углом» (площадью единичной сферы). Осталось вычислить радиальный интеграл $\int_0^\infty r^2 e^{-r^2/2}\,dr$. Известно (из теории гамма-функции, либо из свойства второго момента стандартного нормального распределения — половина полного интеграла по симметрии), что

Подставляем:

Ответ: $Z=(2\pi)^{3/2}\approx15{,}749$ — в точности то же самое значение, что мы получили в предыдущем разделе факторизацией на три независимых одномерных гауссовых интеграла. Это приятная проверка: два совершенно разных пути вычисления одного и того же тройного интеграла (факторизация по осям и переход к сферическим координатам) дают идентичный ответ. И это как раз тот самый концептуальный мостик к статистике: когда плотность зависит только от расстояния до центра (то есть обладает сферической симметрией, как стандартный многомерный гауссиан), переход к сферическим координатам разбивает интеграл на простую угловую часть (которая почти всегда даёт константу — площадь единичной сферы) и радиальную часть, которая часто сводится к табличному интегралу через гамма-функцию.

Почему это важно

Сферические координаты — это не просто ещё один способ вычислить объём шара; это универсальный инструмент для любой задачи с центральной (сферической) симметрией — будь то гравитационное поле, электростатический потенциал точечного заряда или, как мы только что увидели, нормировочная константа многомерного гауссовского распределения. Именно поэтому вывод якобиана $r^2\sin\varphi$ — не формальное упражнение, которое нужно вызубрить, а рабочий инструмент, к которому ты будешь возвращаться каждый раз, когда встретишь функцию, зависящую только от $\sqrt{x^2+y^2+z^2}$. А параллель со статистикой стоит запомнить особо: разложение интеграла на «угловую» часть (константа, связанная с площадью сферы) и «радиальную» часть (одномерный интеграл по $r$, часто через гамма-функцию) — это ровно тот же приём, которым в теории вероятностей выводят распределение хи-квадрат, вычисляют моменты многомерных нормальных распределений и работают с любыми сферически симметричными плотностями произвольной размерности.

Практика: 30 заданий

Базовые задания (1–10)

Задание 1: Вычислить $\iiint_V x\,dV$, где $V=[0,2]\times[0,1]\times[0,3]$.

Задание 2: Найти объём прямоугольного параллелепипеда $[0,2]\times[0,3]\times[0,4]$ через тройной интеграл.

Задание 3: Вычислить $\iiint_V (x+y)\,dV$, где $V=[0,1]^3$ — единичный куб.

Задание 4: Вычислить $\iiint_V z\,dV$, где $V$ — тетраэдр $x,y,z\ge0$, $x+y+z\le1$.

Задание 5: Найти массу куба $[0,1]^3$ с плотностью $\rho(x,y,z)=xyz$.

Задание 6: Найти объём цилиндра радиуса $R=3$, высоты $h=5$ через тройной интеграл в цилиндрических координатах.

Задание 7: Найти объём шара радиуса $R=2$ через тройной интеграл в сферических координатах.

Задание 8: Найти объём тела, ограниченного цилиндром $r\le1$ (то есть $x^2+y^2\le1$), плоскостью $z=0$ снизу и поверхностью $z=x^2+y^2$ сверху.

Задание 9: Записать формулу якобиана перехода к цилиндрическим координатам и объяснить, откуда берётся множитель $r$ (без полного вывода определителя — только итоговая идея).

Задание 10: Записать формулу якобиана перехода к сферическим координатам.

Средние задания (11–20)

Задание 11: Вычислить $\iiint_V xyz\,dV$, где $V=[0,1]\times[0,2]\times[0,3]$.

Задание 12: Найти объём тела, ограниченного цилиндром $x^2+y^2=4$, снизу плоскостью $z=0$, сверху плоскостью $z=x+5$.

Задание 13: Найти массу цилиндра ($R=2$, $h=3$, ось совпадает с осью $z$, основание при $z=0$) с плотностью $\rho(z)=z^2$.

Задание 14 (машинное обучение): Найти дисперсию по одной координате для стандартного трёхмерного гауссиана $f(x,y,z)=(2\pi)^{-3/2}e^{-(x^2+y^2+z^2)/2}$, зная, что по симметрии $E[x^2]=E[y^2]=E[z^2]=\frac13E[r^2]$, и что $E[r^2]=3$.

Задание 15: Найти объём тела между конусом $z=\sqrt{x^2+y^2}$ и плоскостью $z=1$.

Задание 16 (машинное обучение): Ненормированная сферически симметричная плотность $\tilde p(x,y,z)=e^{-r}$, где $r=\sqrt{x^2+y^2+z^2}$. Найти нормировочную константу $Z=\iiint_{\mathbb{R}^3}\tilde p\,dV$.

Задание 17: Вычислить $\iiint_V (x^2+y^2+z^2)\,dV$ по шару радиуса $R$ с центром в начале координат.

Задание 18: Вычислить $\iiint_V\sqrt{x^2+y^2}\,dV$, где $V$ — цилиндр $r\le2$, $0\le z\le3$.

Задание 19: Найти координату $z$ центра масс однородного полушария (верхняя половина шара радиуса $R$, $z\ge0$, постоянная плотность $\rho_0$).

Задание 20: В попытке найти объём шара радиуса $R$ кто-то вычислил $\displaystyle\int_0^{2\pi}\int_0^\pi\int_0^R dr\,d\varphi\,d\theta$ (забыв якобиан $r^2\sin\varphi$). Найти этот неверный результат и сравнить с правильным ответом.

Продвинутые задания (21–30)

Задание 21: Найти объём тела, ограниченного сферой $x^2+y^2+z^2=4$ и конусом $z=\sqrt{3(x^2+y^2)}$, лежащего внутри конуса (то есть ближе к оси $z$), при $z\ge0$.

Задание 22 (машинное обучение): Для стандартного трёхмерного гауссиана вычислить $E[r^2]=\iiint_{\mathbb{R}^3} r^2\cdot(2\pi)^{-3/2}e^{-r^2/2}\,dV$, используя табличный интеграл $\int_0^\infty r^4 e^{-r^2/2}\,dr = \dfrac32\sqrt{2\pi}$.

Задание 23: Найти объём шарового слоя между сферами радиусов $1$ и $2$ с общим центром.

Задание 24: Найти массу шарового слоя из задания 23 с плотностью $\rho(r)=\dfrac1{r^2}$ (типичная модель интенсивности излучения точечного источника, убывающей обратно пропорционально квадрату расстояния).

Задание 25: Найти объём тела, ограниченного снизу конусом $z=\sqrt{x^2+y^2}$, сверху сферой $x^2+y^2+z^2=2$, при $z\ge0$.

Задание 26 (машинное обучение): Для сферически симметричной плотности $p_{\text{норм}}(x,y,z)=\dfrac{e^{-r}}{8\pi}$ из задания 16 найти математическое ожидание расстояния до центра $E[r]$, зная $\int_0^\infty r^3 e^{-r}\,dr=\Gamma(4)=3!=6$.

Задание 27: Вычислить $\iiint_V z\,dV$ по верхней половине шара радиуса $R=3$ ($z\ge0$) двумя способами — через сферические и через цилиндрические координаты — и убедиться, что ответы совпадают.

Задание 28: Для «географической» параметризации сферы $x=r\cos\varphi\cos\theta$, $y=r\cos\varphi\sin\theta$, $z=r\sin\varphi$ (угол $\varphi\in[-\frac{\pi}2,\frac{\pi}2]$ отсчитывается от экваториальной плоскости, а не от полюса) найти якобиан.

Задание 29: Найти координату $z$ центра масс однородного конуса (плотность $\rho_0$, высота $h$, радиус основания $R$, вершина в начале координат, ось вдоль $z$, основание при $z=h$).

Задание 30 (машинное обучение): Используя тот же приём разложения на угловую и радиальную части (или прямую факторизацию), объяснить, почему нормировочная константа стандартного $n$-мерного гауссовского распределения равна $(2\pi)^{n/2}$, и проверить эту формулу при $n=1$ (обычный интеграл Гаусса), $n=2$ (двойной интеграл из прошлого урока) и $n=3$ (этот урок).

Частые ошибки

Разберём типичные ошибки, которые встречаются почти у каждого, кто только осваивает тройные интегралы.

• Забывают якобиан при замене переменных. В цилиндрических координатах $dV=r\,dr\,d\theta\,dz$, в сферических — $dV=r^2\sin\varphi\,dr\,d\varphi\,d\theta$. Пропуск множителя $r$ или $r^2\sin\varphi$ — самая частая и самая грубая ошибка на этой теме: как показано в задании 20, без якобиана результат может получиться не просто неточным, а даже неправильной размерности (не куб линейного размера, а первая степень).

• Путают полярный угол $\varphi$ и азимутальный угол $\theta$ в сферических координатах, а также путают пределы для $\varphi$. Полярный угол $\varphi$ отсчитывается от оси $z$ и меняется от $0$ до $\pi$ (а не от $0$ до $2\pi$ — это частая путаница с азимутальным углом), тогда как азимутальный угол $\theta$ в плоскости $xy$ меняется от $0$ до $2\pi$.

• Неверно описывают область интегрирования, не построив проекцию тела на плоскость. Тройной интеграл требует аккуратного описания того, как зависят пределы внутренней переменной ($z$, например) от значений внешних переменных ($x,y$) — без чёткого понимания геометрии тела легко перепутать верхнюю и нижнюю границы или упустить часть области.

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

• Забывают, что при замене переменных нужно не только подставить новые координаты в саму функцию, но и в пределы интегрирования, и добавить якобиан. Три отдельных шага, и пропуск любого из них даёт неверный результат: подстановка функции без пределов, пределы без якобиана или, наоборот, якобиан без корректной подстановки функции — распространённые частичные ошибки.

• Путают порядок множителей при разложении факторизуемого интеграла. Если подынтегральная функция и область распадаются на независимые части по каждой переменной (как в примерах с кубом и гауссианом), тройной интеграл превращается в произведение трёх однократных интегралов — но эта факторизация работает только тогда, когда и функция, и область действительно разделяются; для тетраэдра или шара, где пределы одной переменной зависят от других, так поступать нельзя.

Главное запомнить

• Тройной интеграл $\iiint_V f(x,y,z)\,dV$ — это предел суммы вкладов $f\cdot\Delta V$ по всем маленьким элементам объёма, на которые разбито тело $V$.

• Геометрический смысл: при $f\equiv1$ тройной интеграл равен объёму тела $V$.

• Физический смысл: при $f=\rho(x,y,z)$ (плотность) тройной интеграл равен массе тела.

• Тройной интеграл сводится к трём последовательным однократным интегрированиям — обобщение теоремы Фубини на трёхмерный случай; порядок интегрирования можно менять, выбирая самый удобный.

• В цилиндрических координатах $x=r\cos\theta$, $y=r\sin\theta$, $z=z$ якобиан равен $r$, элемент объёма $dV=r\,dr\,d\theta\,dz$ — используются для тел с осевой симметрией.

• В сферических координатах $x=r\sin\varphi\cos\theta$, $y=r\sin\varphi\sin\theta$, $z=r\cos\varphi$ якобиан равен $r^2\sin\varphi$, элемент объёма $dV=r^2\sin\varphi\,dr\,d\varphi\,d\theta$ — используются для тел со сферической симметрией.

• Оба якобиана получены прямым вычислением определителя матрицы Якоби, а не заучены по памяти — при другой параметризации углов (как в задании про «географические» координаты) якобиан меняется, и его нужно выводить заново.

• Если подынтегральная функция и область интегрирования факторизуются по осям (как прямоугольная область и сумма или произведение независимых слагаемых), тройной интеграл распадается на произведение трёх простых однократных интегралов.

• Нормировочная константа стандартного $n$-мерного гауссовского распределения равна $(2\pi)^{n/2}$ — при $n=3$ это ровно $(2\pi)^{3/2}$, что можно получить как факторизацией по осям, так и переходом к сферическим координатам, и оба пути дают одинаковый ответ.

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

Связь с другими темами курса

Этот урок напрямую опирается на двойные интегралы из урока 215: тройной интеграл — это то же самое понятие интегральной суммы, только на одно измерение больше, а переход к полярным координатам, освоенный в прошлом уроке, — это ровно тот же приём, что цилиндрические координаты применяют к горизонтальному сечению трёхмерного тела. Теорема Фубини, сводящая двойной интеграл к повторному, дословно обобщается на тройной интеграл — просто добавляется ещё одно интегрирование.

Вперёд этот урок ведёт к криволинейным интегралам — интегрированию вдоль кривых в пространстве, где вместо элемента объёма или площади фигурирует элемент длины дуги. Дальше в курсе последуют поверхностные интегралы (интегрирование по двумерным поверхностям, вложенным в трёхмерное пространство) и теоремы Гаусса-Остроградского и Стокса, которые связывают тройной интеграл по объёму с поверхностным интегралом по границе этого объёма — прямое развитие идей сегодняшнего урока. А для машинного обучения и статистики тройные интегралы (и их обобщение на произвольную размерность) — это тот самый инструмент, которым вычисляют нормировочные константы многомерных плотностей вероятности, моменты распределений и математические ожидания функций от случайных векторов; именно с этой техникой ты столкнёшься, изучая многомерное нормальное распределение, метод Монте-Карло или байесовский вывод.

Интересные факты

• Формула $E=mc^2$ здесь ни при чём, а вот вычисление гравитационного притяжения однородного шара — классическая задача, которую ещё Исаак Ньютон решил, по сути, эквивалентными геометрическими методами задолго до появления современной записи через тройной интеграл в сферических координатах; результат (что однородный шар притягивает внешнюю точку так, как будто вся его масса сосредоточена в центре) — прямое следствие сферической симметрии, той же самой, которую мы использовали сегодня для вычисления объёма и массы шара.

• Слово «якобиан» происходит от имени немецкого математика Карла Густава Якоба Якоби, который в XIX веке систематически изучал определители матриц частных производных при замене переменных в кратных интегралах — задолго до того, как эти определители стали стандартным инструментом в учебниках математического анализа по всему миру.

• Нормировочная константа $(2\pi)^{n/2}$ многомерного стандартного гауссовского распределения появляется не только в статистике, но и в физике — в частности, в статистической механике при вычислении статистических сумм для систем со многими степенями свободы, где каждая степень свободы вносит свой собственный гауссов множитель.

• В компьютерной графике и физически обоснованном рендеринге тройные интегралы по объёму используются для моделирования рассеяния света в туманах, дыме и других объёмных средах — интеграл считает суммарный вклад рассеянного света вдоль всего объёма, через который проходит луч.

Лайфхаки и полезные трюки

• Прежде чем садиться писать пределы интегрирования, нарисуй (хотя бы схематично) проекцию тела на плоскость $xy$ и подумай, как выглядит тело в разрезе на разной высоте $z$ — большая часть ошибок в тройных интегралах происходит именно от неверного представления геометрии тела, а не от вычислительных ошибок.

• Если тело или подынтегральная функция обладают осевой симметрией (зависят только от расстояния до оси, не от угла), сразу переходи к цилиндрическим координатам — интеграл по углу $\theta$ почти всегда даст тривиальный множитель $2\pi$ (или его долю), и работа сведётся к двум простым одномерным интегралам.

• Если тело или функция обладают сферической (центральной) симметрией — зависят только от расстояния до начала координат, — переходи к сферическим координатам: угловая часть (по $\varphi$ и $\theta$) почти всегда даёт константу $4\pi$ (полный телесный угол) или её долю, а вся содержательная работа сводится к одномерному интегралу по $r$.

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

• Хороший способ проверить себя после вычисления объёма или массы через замену координат — прикинуть порядок величины результата на глазок (например, для шара радиуса $2$ объём должен быть в районе $30$-$35$, поскольку $\frac{4}{3}\pi\cdot8\approx33{,}5$) — если получившееся число сильно отличается от прикидки, скорее всего где-то забыт якобиан или перепутаны пределы.

• Если ты сомневаешься в выведенном якобиане, проверь его на простом частном случае: для сферических координат подставь $\varphi=\frac{\pi}2$ (экваториальная плоскость) — якобиан $r^2\sin\varphi$ становится равным $r^2$, что совпадает с элементом площади $r\,dr\,d\theta$, умноженным на дополнительный радиальный множитель $r$ (интуиция: на экваторе сферические координаты почти совпадают с цилиндрическими).

Тройные интегралы завершают трилогию «интеграл по области возрастающей размерности», начатую обычным определённым интегралом и продолженную двойным интегралом, — и в этом смысле сегодняшний урок не столько добавляет что-то принципиально новое, сколько показывает, что вся эта идея прекрасно масштабируется на любое число измерений. Работа с объёмом и массой реальных трёхмерных тел, вывод якобианов цилиндрических и сферических координат, вычисление нормировочных констант многомерных вероятностных распределений — всё это варианты одной и той же математической идеи: разбить сложную область на бесконечно много бесконечно малых кусочков, посчитать вклад каждого и сложить всё вместе. Освоив эту идею по-настоящему на примере трёх измерений, ты без труда перенесёшь её и на четыре, и на сто измерений — а именно так выглядит математика, стоящая за многомерными вероятностными моделями в машинном обучении.