Двойные интегралы 🗺️

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

Если обычный определённый интеграл $\int_a^b f(x)\,dx$ ты представляешь как площадь под кривой на плоскости, то двойной интеграл $\iint_D f(x,y)\,dA$ — это прямое обобщение той же идеи на одно измерение выше: объём тела, зажатого между плоскостью $xy$ и поверхностью $z=f(x,y)$, над некоторой областью $D$. Формально это снова предел интегральных сумм, только теперь суммируются не узкие вертикальные полоски высотой $f(x_i)$ и шириной $\Delta x$, а маленькие кусочки поверхности площадью $f(x_i,y_i)\cdot\Delta A$, где $\Delta A$ — площадь крошечного кусочка области $D$ на плоскости.

Для тебя как будущего специалиста по машинному обучению и анализу данных двойной интеграл — это далеко не абстрактная гимнастика. Совместная плотность вероятности двух случайных величин $f(x,y)$ описывает, как «размазана» вероятность по плоскости, а вероятность того, что пара случайных величин попадёт в конкретную область $D$ (скажем, в прямоугольник допустимых значений двух признаков), вычисляется ровно как двойной интеграл $P((X,Y)\in D)=\iint_D f(x,y)\,dx\,dy$. Условие нормировки плотности («суммарная вероятность равна единице») — это тоже двойной интеграл, только по всей плоскости. А знаменитый интеграл Гаусса $\int_{-\infty}^{\infty}e^{-x^2}dx=\sqrt\pi$, который лежит в основе нормального распределения — того самого распределения, которое встречается в статистике, инициализации весов нейросетей и байесовских моделях на каждом шагу, — вычисляется не напрямую, а хитрым трюком через переход к двойному интегралу и последующую замену на полярные координаты. Именно этот трюк мы подробно разберём сегодня.

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

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

Идея интегрирования по области на плоскости зрела в математике постепенно, вместе с общим развитием интегрального исчисления в XVIII веке. Ньютон и Лейбниц в конце XVII века построили аппарат для интегрирования функций одной переменной, но задачи, требующие суммирования по двумерной области — вычисление объёмов сложных тел, моментов инерции пластин, притяжения тел друг к другу по закону всемирного тяготения — оставались вне единой теории и решались случай за случаем, изобретательными, но разрозненными геометрическими приёмами.

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

Отдельного упоминания заслуживает история перехода к полярным координатам при вычислении интеграла Гаусса $\int_{-\infty}^{\infty}e^{-x^2}dx$. Этот интеграл принципиально не берётся в элементарных функциях: первообразную для $e^{-x^2}$ невозможно выразить через привычный набор функций — многочлены, экспоненты, логарифмы, тригонометрию — сколько ни старайся. Тем не менее сам определённый интеграл по всей числовой прямой имеет изящное точное значение $\sqrt\pi$, и первым, кто нашёл трюк для его вычисления — возвести искомый интеграл в квадрат, превратить квадрат в двойной интеграл по всей плоскости, а затем перейти к полярным координатам, где интеграл внезапно становится элементарно берущимся, — обычно называют Пьера-Симона Лапласа (около 1810-х годов), хотя схожие результаты в связи с нормальным распределением и теорией ошибок наблюдений примерно в то же время получал и Карл Фридрих Гаусс, чьим именем в итоге и назвали и сам интеграл, и построенное на нём распределение вероятностей. Этот трюк — один из самых красивых во всей математике: он показывает, что переход в более высокую размерность (от одной переменной к двум) иногда упрощает задачу, которая казалась неразрешимой в исходной постановке.

Двойной интеграл как предел сумм: определение и геометрический смысл

Интуиция: от узких полосок к маленьким кусочкам плоскости

Давай разберёмся, как обобщить обычный интеграл на случай двух переменных. Вспомни определение однократного интеграла: ты разбиваешь отрезок $[a,b]$ на узкие полоски шириной $\Delta x_i$, в каждой полоске выбираешь точку $x_i$, умножаешь значение функции $f(x_i)$ на ширину полоски $\Delta x_i$ и складываешь всё это по всем полоскам. При стремлении ширины полосок к нулю сумма стремится к площади под графиком — это и есть определённый интеграл.

Представь теперь, что вместо отрезка на прямой у тебя есть плоская область $D$ произвольной формы — скажем, круг, треугольник или что-то ещё более причудливое, — и над каждой точкой этой области подвешена высота $z=f(x,y)$, задающая некую поверхность. Разбей область $D$ мелкой сеткой на маленькие кусочки (пусть даже неправильной формы) площадью $\Delta A_i$ каждый. В каждом кусочке выбери произвольную точку $(x_i,y_i)$ и построй тонкий столбик с основанием $\Delta A_i$ и высотой $f(x_i,y_i)$ — его объём приближённо равен $f(x_i,y_i)\cdot\Delta A_i$. Сложи объёмы всех таких столбиков по всей области — получишь приближение к объёму тела, зажатого между плоскостью $xy$ и поверхностью $z=f(x,y)$. Чем мельче разбиение, тем точнее приближение, а предел этой суммы при стремлении диаметра всех кусочков к нулю и называется двойным интегралом.

Определение

Определение (двойной интеграл). Пусть функция $f(x,y)$ определена на ограниченной области $D$ плоскости. Разобьём область $D$ на $n$ малых частей площадью $\Delta A_1,\dots,\Delta A_n$, в каждой части выберем произвольную точку $(x_i,y_i)$ и составим интегральную сумму

$$S_n = \sum_{i=1}^n f(x_i,y_i)\,\Delta A_i$$

Если при стремлении наибольшего диаметра всех частей разбиения к нулю эта сумма стремится к одному и тому же конечному пределу независимо от способа разбиения и выбора точек $(x_i,y_i)$, этот предел называется двойным интегралом функции $f$ по области $D$:

$$\iint_D f(x,y)\,dA = \lim_{\max\Delta A_i\to0}\sum_{i=1}^n f(x_i,y_i)\,\Delta A_i$$

Величину $dA$ (или, что то же самое, $dx\,dy$) называют элементом площади. Функция $f$, для которой этот предел существует, называется интегрируемой на $D$. Достаточное условие интегрируемости — непрерывность $f$ на замкнутой ограниченной области $D$ с достаточно «хорошей» (кусочно-гладкой) границей.

Геометрический смысл двойного интеграла зависит от знака функции $f$. Если $f(x,y)\ge0$ на всей области $D$, интеграл $\iint_D f(x,y)\,dA$ равен объёму тела, ограниченного снизу плоскостью $z=0$, сверху — поверхностью $z=f(x,y)$, а с боков — вертикальной «стенкой», опирающейся на границу области $D$. Частный, но исключительно важный случай: если $f(x,y)\equiv1$ на всей области, интеграл $\iint_D 1\,dA$ в точности равен площади самой области $D$ — это прямой аналог того факта, что $\int_a^b1\,dx=b-a$ равен длине отрезка.

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

Пример 1 (лёгкий, вычисление площади через двойной интеграл). Показать, что двойной интеграл функции $f(x,y)=1$ по прямоугольнику $D=\{0\le x\le3,\ 0\le y\le2\}$ равен площади этого прямоугольника.

Разобьём прямоугольник на $n\times m$ мелких прямоугольничков со сторонами $\Delta x$ и $\Delta y$, тогда $\Delta A_i=\Delta x\,\Delta y$ для каждого кусочка. Интегральная сумма: $S=\sum_i 1\cdot\Delta A_i=\sum_i\Delta A_i$ — сумма площадей всех кусочков разбиения, которая по построению равна площади всего прямоугольника, независимо от мелкости разбиения. Значит, и предел (то есть сам двойной интеграл) равен этой же величине.

Ответ: $\iint_D 1\,dA = 3\cdot2=6$ — площадь прямоугольника. Этот пример — не просто формальность: он показывает, что двойной интеграл действительно обобщает понятие площади, а не просто напоминает его по аналогии.

Пример 2 (средний, интерпретация как объём). Функция $f(x,y)=5$ (постоянная высота) задана над кругом $D$ радиуса $2$ с центром в начале координат. Найти объём цилиндра, который представляет собой соответствующий двойной интеграл, не вычисляя его напрямую как предел сумм.

Поскольку $f(x,y)=5$ постоянна, тело под поверхностью $z=5$ над кругом радиуса $2$ — это просто прямой круговой цилиндр высотой $5$ с площадью основания $\pi\cdot2^2=4\pi$. Объём цилиндра — площадь основания, умноженная на высоту.

Ответ: $\iint_D 5\,dA=20\pi\approx62{,}83$. Обрати внимание: для постоянной функции двойной интеграл всегда сводится к произведению константы на площадь области — вычислять предел сумм заново для такого случая не нужно, это прямое следствие определения (постоянный множитель выносится из суммы, а сумма площадей кусочков разбиения по определению стремится к площади области).

Пример 3 (сложный, машинное обучение — нормировка совместной плотности). Пара случайных величин $(X,Y)$ имеет совместную плотность вероятности $f(x,y)=c$ (постоянную) на прямоугольнике $D=\{0\le x\le2,\ 0\le y\le5\}$ и $f(x,y)=0$ вне этого прямоугольника (равномерное распределение на прямоугольнике). Найти константу $c$, при которой $f$ действительно является плотностью вероятности.

Ключевое условие нормировки любой плотности вероятности: суммарная вероятность по всей плоскости должна быть равна единице, то есть $\iint_{\mathbb R^2} f(x,y)\,dA=1$. Поскольку $f=0$ вне прямоугольника $D$, это условие сводится к $\iint_D c\,dA=1$. Площадь прямоугольника $D$ равна $2\cdot5=10$, значит по формуле для постоянной функции:

Ответ: $c=0{,}1$. Это стандартная задача теории вероятностей: чтобы функция $f(x,y)$ была честной совместной плотностью, двойной интеграл от неё по всей области определения обязан равняться единице, — и в частном случае равномерного распределения это условие моментально сводится к простому геометрическому факту «площадь, умноженная на высоту постоянной плотности, равна единице».

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

Определение через предел интегральных сумм — это фундамент, но на практике вычислять двойные интегралы напрямую по определению (как предел сумм) почти никогда не приходится: это было бы так же неудобно, как вычислять обычный однократный интеграл через предел сумм Римана каждый раз заново. Нужен рабочий вычислительный инструмент — и таким инструментом служит сведение двойного интеграла к повторному, к которому мы переходим прямо сейчас.

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

Интуиция: нарезаем тело на тонкие пластинки

Давай разберёмся, как на практике вычисляют двойные интегралы. Представь снова тело под поверхностью $z=f(x,y)$ над прямоугольной областью $D=\{a\le x\le b,\ c\le y\le d\}$. Зафиксируй конкретное значение $x$ и разрежь тело вертикальной плоскостью, перпендикулярной оси $x$, в этой точке. Получившееся сечение — плоская фигура, ограниченная кривой $z=f(x,y)$ (где $x$ уже зафиксирован, а $y$ бежит от $c$ до $d$) — обычная одномерная функция от $y$, и площадь этого сечения вычисляется обычным однократным интегралом $\int_c^d f(x,y)\,dy$ (при фиксированном $x$).

Теперь представь, что ты «нарезаешь» всё тело на бесконечно много таких тонких пластинок толщиной $dx$, каждая площадью сечения $\int_c^d f(x,y)\,dy$, и складываешь их объёмы по всем значениям $x$ от $a$ до $b$ — то есть берёшь ещё один, внешний, интеграл по $x$. Это и есть повторный интеграл: сначала «внутри» интегрируешь по одной переменной при фиксированной другой (получаешь функцию, зависящую только от оставшейся переменной), а затем «снаружи» интегрируешь результат по этой оставшейся переменной.

Определение (теорема Фубини)

Теорема (сведение двойного интеграла к повторному, теорема Фубини). Пусть функция $f(x,y)$ непрерывна на прямоугольнике $D=\{a\le x\le b,\ c\le y\le d\}$. Тогда двойной интеграл равен повторному, причём порядок интегрирования можно выбирать свободно:

$$\iint_D f(x,y)\,dA = \int_a^b\left(\int_c^d f(x,y)\,dy\right)dx = \int_c^d\left(\int_a^b f(x,y)\,dx\right)dy$$

При вычислении внутреннего интеграла переменная внешнего интегрирования считается постоянной (в точности как при вычислении частной производной).

Обозначение повторного интеграла часто пишут без явных скобок: $\displaystyle\int_a^b\int_c^d f(x,y)\,dy\,dx$ — и здесь порядок записи дифференциалов важен: $dy\,dx$ означает, что сначала интегрируем по $y$ (внутренний интеграл, ближайший к функции), а уже потом по $x$ (внешний). Именно поэтому порядок дифференциалов в записи повторного интеграла нельзя переставлять небрежно — он однозначно говорит, какая переменная фиксируется первой.

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

Пример 1 (лёгкий). Вычислить $\displaystyle\iint_D (x+y)\,dA$ по прямоугольнику $D=\{0\le x\le2,\ 0\le y\le1\}$.

Сводим к повторному интегралу, интегрируя сначала по $y$ при фиксированном $x$:

Внутренний интеграл (переменная $x$ — константа):

Внешний интеграл:

Ответ: $\iint_D(x+y)\,dA=3$. Проверим порядок в обратную сторону — интегрируя сначала по $x$: внутренний интеграл $\int_0^2(x+y)\,dx=\left[\frac{x^2}2+xy\right]_0^2=2+2y$, внешний $\int_0^1(2+2y)\,dy=[2y+y^2]_0^1=2+1=3$ — тот же ответ, как и гарантирует теорема Фубини.

Пример 2 (средний). Вычислить $\displaystyle\iint_D xy^2\,dA$ по прямоугольнику $D=\{1\le x\le3,\ 0\le y\le2\}$.

Внутренний интеграл (константа $x$ выносится за знак интеграла по $y$):

Внешний интеграл:

Ответ: $\iint_D xy^2\,dA=\dfrac{32}3\approx10{,}667$.

Пример 3 (сложный, машинное обучение — вероятность попадания в прямоугольник). Пара случайных величин $(X,Y)$ имеет совместную плотность $f(x,y)=\dfrac{1}{4}xy$ на прямоугольнике $\{0\le x\le2,\ 0\le y\le2\}$ (и равна нулю вне него). Убедиться, что это действительно плотность вероятности (то есть двойной интеграл по всему прямоугольнику равен единице), а затем найти вероятность $P(0\le X\le1,\ 1\le Y\le2)$.

Проверка нормировки:

Отлично, нормировка подтверждена — это честная плотность. Теперь вычисляем вероятность попадания в прямоугольник $\{0\le x\le1,\ 1\le y\le2\}$ — это просто ещё один двойной интеграл той же самой функции, но по меньшей области:

Внутренний интеграл: $\displaystyle\int_1^2\frac{xy}4\,dy = \frac{x}4\left[\frac{y^2}2\right]_1^2 = \frac{x}4\cdot\frac{4-1}2=\frac{x}4\cdot\frac32=\frac{3x}8$.

Внешний интеграл: $\displaystyle\int_0^1\frac{3x}8\,dx = \frac38\left[\frac{x^2}2\right]_0^1 = \frac38\cdot\frac12=\frac3{16}=0{,}1875$.

Ответ: $P(0\le X\le1,\,1\le Y\le2)=\dfrac3{16}=0{,}1875$. Это в точности та операция, которую ты будешь выполнять при работе с многомерными случайными величинами: вероятность попадания в область — это двойной (а в общем случае — кратный) интеграл совместной плотности по этой области, и техника сведения к повторному интегралу здесь работает без каких-либо изменений.

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

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

Правильные области интегрирования и расстановка пределов

Интуиция: не всякая область — прямоугольник

Пока мы работали только с прямоугольными областями, где пределы интегрирования по обеим переменным — простые константы. Но на практике область $D$ почти никогда не бывает прямоугольником: это может быть треугольник, круг, область между двумя кривыми и так далее. Давай разберёмся, как расставлять пределы интегрирования в таком общем случае — и именно здесь чаще всего допускаются ошибки, поэтому отнесись к этому разделу особенно внимательно.

Представь область $D$ на плоскости и мысленно проведи через неё вертикальную линию (при фиксированном $x$). Если для каждого $x$ из некоторого интервала $[a,b]$ эта вертикальная линия входит в область $D$ и выходит из неё ровно один раз — то есть входит через нижнюю границу $y=g_1(x)$ и выходит через верхнюю границу $y=g_2(x)$ — область называется $x$-правильной (по-другому её называют элементарной относительно оси $x$, или областью первого типа). Аналогично, если каждая горизонтальная линия (при фиксированном $y$) пересекает границу области ровно один раз «снизу» и один раз «сверху» относительно $x$, область называется $y$-правильной.

Определение

Определение ($x$-правильная и $y$-правильная область). Область $D$ называется $x$-правильной, если она представима в виде

$$D = \{(x,y):\ a\le x\le b,\ \ g_1(x)\le y\le g_2(x)\}$$

где $g_1,g_2$ — непрерывные функции на $[a,b]$. Для такой области двойной интеграл сводится к повторному так:

$$\iint_D f(x,y)\,dA = \int_a^b\left(\int_{g_1(x)}^{g_2(x)} f(x,y)\,dy\right)dx$$

Область $D$ называется $y$-правильной, если она представима в виде

$$D = \{(x,y):\ c\le y\le d,\ \ h_1(y)\le x\le h_2(y)\}$$

и тогда

$$\iint_D f(x,y)\,dA = \int_c^d\left(\int_{h_1(y)}^{h_2(y)} f(x,y)\,dx\right)dy$$

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

Практический алгоритм расстановки пределов для $x$-правильной области такой: нарисуй область $D$ (это первый и самый важный шаг — без чертежа расставить пределы правильно почти невозможно); найди диапазон $[a,b]$, в котором вообще «живёт» переменная $x$ по всей области, — это пределы внешнего интеграла; затем для произвольного фиксированного $x$ из этого диапазона проведи вертикальную линию и определи, между какими кривыми $y=g_1(x)$ (снизу) и $y=g_2(x)$ (сверху) она входит и выходит из области, — это пределы внутреннего интеграла, зависящие от $x$.

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

Пример 1 (лёгкий, треугольник). Вычислить $\displaystyle\iint_D xy\,dA$, где $D$ — треугольник с вершинами $(0,0)$, $(2,0)$, $(0,2)$ (то есть область, ограниченная осями координат и прямой $x+y=2$).

Рисуем область: это треугольник, лежащий под прямой $y=2-x$, для $x$ от $0$ до $2$. Значит, область $x$-правильная: $0\le x\le2$, а для каждого фиксированного $x$ переменная $y$ бежит от $0$ (снизу, ось $x$) до $2-x$ (сверху, гипотенуза).

Внутренний интеграл: $\displaystyle\int_0^{2-x}xy\,dy = x\left[\frac{y^2}2\right]_0^{2-x} = \frac{x(2-x)^2}2$.

Внешний интеграл: раскроем $(2-x)^2=4-4x+x^2$, получим подынтегральное выражение $\dfrac{x(4-4x+x^2)}2=\dfrac{4x-4x^2+x^3}2=2x-2x^2+\dfrac{x^3}2$.

Ответ: $\iint_D xy\,dA=\dfrac23\approx0{,}667$.

Пример 2 (средний, область между параболой и прямой). Вычислить $\displaystyle\iint_D (x+2y)\,dA$, где $D$ — область, ограниченная параболой $y=x^2$ снизу и прямой $y=x+2$ сверху.

Сначала найдём точки пересечения границ: $x^2=x+2\Rightarrow x^2-x-2=0\Rightarrow(x-2)(x+1)=0\Rightarrow x=-1$ или $x=2$. При $x\in[-1,2]$ прямая $y=x+2$ проходит выше параболы $y=x^2$ (проверка в точке $x=0$: прямая даёт $2$, парабола даёт $0$ — да, прямая выше). Значит, область $x$-правильная: $-1\le x\le2$, $x^2\le y\le x+2$.

Внутренний интеграл: $\displaystyle\int_{x^2}^{x+2}(x+2y)\,dy = \Big[xy+y^2\Big]_{x^2}^{x+2}$.

Подставляем верхний предел: $x(x+2)+(x+2)^2 = x^2+2x+x^2+4x+4=2x^2+6x+4$.

Подставляем нижний предел: $x\cdot x^2+(x^2)^2=x^3+x^4$.

Разность: $2x^2+6x+4-x^3-x^4$.

Внешний интеграл:

При $x=2$: $-\dfrac{32}5-\dfrac{16}4+\dfrac{16}3+12+8 = -6{,}4-4+5{,}333+12+8=14{,}933$.

При $x=-1$: $-\dfrac{-1}5-\dfrac14+\dfrac{-2}3+3-4 = 0{,}2-0{,}25-0{,}667+3-4=-1{,}717$.

Разность: $14{,}933-(-1{,}717)=16{,}65$.

Ответ: $\iint_D(x+2y)\,dA\approx16{,}65$ (точнее, $\dfrac{999}{60}=\dfrac{333}{20}=16{,}65$ — можно проверить приведением всех дробей к общему знаменателю $60$).

Пример 3 (сложный, расстановка пределов двумя способами). Область $D$ — четверть круга радиуса $2$ в первой координатной четверти: $x^2+y^2\le4$, $x\ge0$, $y\ge0$. Расставить пределы интегрирования как для $x$-правильной, так и для $y$-правильной области (без вычисления самого интеграла — эту область мы позже пересчитаем через полярные координаты).

Как $x$-правильная область: $x$ пробегает от $0$ до $2$, а при фиксированном $x$ переменная $y$ бежит от $0$ (снизу, ось $x$) до $\sqrt{4-x^2}$ (сверху, дуга окружности, выраженная из уравнения $x^2+y^2=4$). Получаем:

Как $y$-правильная область: по симметрии ролей $x$ и $y$ в уравнении окружности, $y$ пробегает от $0$ до $2$, а при фиксированном $y$ переменная $x$ бежит от $0$ до $\sqrt{4-y^2}$:

Ответ: оба представления пределов корректны и дают одинаковый результат для любой непрерывной $f$ — это прямая иллюстрация того, что четверть круга является одновременно и $x$-правильной, и $y$-правильной областью, а выбор конкретного порядка интегрирования — вопрос удобства вычислений, а не корректности постановки задачи. При этом обрати внимание: подынтегральное выражение с корнями $\sqrt{4-x^2}$ часто оказывается неудобным для непосредственного вычисления в декартовых координатах — и именно для таких круговых областей ниже мы введём переход к полярным координатам, который полностью снимает эту проблему.

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

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

Замена порядка интегрирования

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

Мы уже видели в предыдущем разделе, что одну и ту же область можно описать и как $x$-правильную, и как $y$-правильную — и оба описания дают один и тот же двойной интеграл. Иногда, однако, порядок интегрирования, который кажется естественным на первый взгляд, приводит к внутреннему интегралу, который просто не берётся в элементарных функциях (например, $\int e^{-y^2}dy$ не выражается через элементарные функции ни при каких пределах, кроме особых случаев). В таких ситуациях спасает смена порядка интегрирования: заново нарисовать ту же область как правильную относительно другой переменной и переписать пределы.

Алгоритм смены порядка

Правило (смена порядка интегрирования). Чтобы поменять порядок интегрирования в повторном интеграле $\displaystyle\int_a^b\int_{g_1(x)}^{g_2(x)} f(x,y)\,dy\,dx$:

1. По исходным пределам восстанови и нарисуй саму область интегрирования $D$ на плоскости (это описание области как $x$-правильной).

2. Взгляни на ту же нарисованную область как на $y$-правильную: найди диапазон значений $y$ по всей области — это новые пределы внешнего интеграла — и для каждого $y$ выясни, между какими кривыми $x=h_1(y)$ и $x=h_2(y)$ заключена область — это новые пределы внутреннего интеграла.

3. Перепиши интеграл в новом порядке: $\displaystyle\int_c^d\int_{h_1(y)}^{h_2(y)} f(x,y)\,dx\,dy$. Значение интеграла при этом не меняется — меняется только порядок вычислений.

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

Пример 1 (лёгкий). Дан повторный интеграл $\displaystyle\int_0^1\int_0^x f(x,y)\,dy\,dx$. Изменить порядок интегрирования.

Восстанавливаем область: $0\le x\le1$, $0\le y\le x$ — это треугольник с вершинами $(0,0)$, $(1,0)$, $(1,1)$, ограниченный сверху прямой $y=x$, снизу прямой $y=0$, справа прямой $x=1$.

Взглянем как на $y$-правильную область: $y$ пробегает от $0$ до $1$. При фиксированном $y$ переменная $x$ бежит от $y$ (левая граница — прямая $y=x$, то есть $x=y$) до $1$ (правая граница).

Ответ: $\displaystyle\int_0^1\int_0^x f(x,y)\,dy\,dx = \int_0^1\int_y^1 f(x,y)\,dx\,dy$.

Пример 2 (средний, интеграл, который не берётся в исходном порядке). Вычислить $\displaystyle\int_0^1\int_x^1 e^{y^2}\,dy\,dx$, поменяв порядок интегрирования.

Внутренний интеграл $\int_x^1 e^{y^2}\,dy$ не выражается в элементарных функциях (первообразная $e^{y^2}$ не элементарна) — значит, вычислить в исходном порядке напрямую невозможно, и нужно менять порядок.

Восстанавливаем область: $0\le x\le1$, $x\le y\le1$ — треугольник с вершинами $(0,0)$, $(0,1)$, $(1,1)$ (ограничен слева осью $y$, снизу-справа прямой $y=x$, сверху прямой $y=1$).

Взглянем как на $y$-правильную область: $y$ пробегает от $0$ до $1$. При фиксированном $y$ переменная $x$ бежит от $0$ до $y$ (так как условие $x\le y$ переписывается как $0\le x\le y$).

Теперь внутренний интеграл берётся элементарно, поскольку $e^{y^2}$ здесь — константа относительно $x$:

Внешний интеграл — заменой $u=y^2$, $du=2y\,dy$:

Ответ: $\displaystyle\int_0^1\int_x^1 e^{y^2}\,dy\,dx = \dfrac{e-1}2\approx0{,}859$. Это классический пример того, зачем вообще нужна смена порядка интегрирования: в одном порядке интеграл принципиально не берётся в элементарных функциях, а в другом порядке — берётся буквально в два шага.

Пример 3 (сложный, область с двумя «кусками» после смены порядка). Дан интеграл $\displaystyle\int_0^1\int_{y}^{\sqrt y} f(x,y)\,dx\,dy$. Изменить порядок интегрирования.

Восстанавливаем область: $0\le y\le1$, $y\le x\le\sqrt y$. Заметим важную деталь — при $0\le y\le1$ всегда $y\le\sqrt y$ (например, при $y=0{,}25$: $y=0{,}25$, $\sqrt y=0{,}5$), так что пределы корректны и область не пустая: это область, зажатая между параболой $y=x^2$ (левая граница $x=\sqrt y$ переписывается как $y=x^2$) и прямой $y=x$ (правая граница $x=y$), для $x$ и $y$ от $0$ до $1$.

Точнее: условие $y\le x$ означает, что мы выше прямой $y=x$ снизу нет, наоборот — $x\ge y$ означает точка правее прямой (или прямая выше точки: $y\le x$). А условие $x\le\sqrt y$ равносильно $x^2\le y$ (при $x,y\ge0$), то есть мы ниже параболы $y=x^2$... Проверим на конкретной точке, чтобы не запутаться: возьмём $y=0{,}25$, тогда $x$ должен лежать между $0{,}25$ и $0{,}5$. Точка $x=0{,}4$, $y=0{,}25$: проверяем $y\le x$: $0{,}25\le0{,}4$ — да. Проверяем $x\le\sqrt y=0{,}5$: $0{,}4\le0{,}5$ — да. И относительно параболы $y=x^2$: при $x=0{,}4$ парабола даёт $y=0{,}16$, а у нас $y=0{,}25>0{,}16$ — то есть точка выше параболы. А относительно прямой $y=x$: при $x=0{,}4$ прямая даёт $y=0{,}4$, а у нас $y=0{,}25<0{,}4$ — то есть точка ниже прямой. Значит, область — это часть плоскости между параболой $y=x^2$ (снизу) и прямой $y=x$ (сверху), для $x$ от $0$ до $1$.

Взглянем как на $x$-правильную область: $x$ пробегает от $0$ до $1$, а при фиксированном $x$ переменная $y$ бежит от $x^2$ (снизу, парабола) до $x$ (сверху, прямая) — обе кривые пересекаются как раз в точках $x=0$ и $x=1$, что согласуется с диапазоном.

Ответ: $\displaystyle\int_0^1\int_{y}^{\sqrt y} f(x,y)\,dx\,dy = \int_0^1\int_{x^2}^{x} f(x,y)\,dy\,dx$. Этот пример показывает, насколько важно не полагаться на «автоматическое» жонглирование пределами, а каждый раз явно восстанавливать и проверять форму области хотя бы на одной пробной точке.

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

Смена порядка интегрирования — это не просто формальный трюк «на всякий случай», а рабочий инструмент, который регулярно спасает ситуацию, когда исходный порядок приводит к неберущемуся внутреннему интегралу (как в примере 2 выше, с $e^{y^2}$, — а такая ситуация в реальных вероятностных и физических задачах встречается сплошь и рядом, поскольку экспонента от квадрата переменной — это как раз ядро нормального распределения). Умение бегло перерисовывать одну и ту же область как $x$-правильную и как $y$-правильную — это практический навык, без которого целый класс задач попросту недоступен.

Замена переменных в двойном интеграле: якобиан и полярные координаты

Интуиция: растягиваем и искажаем координатную сетку

Иногда область интегрирования или подынтегральная функция настолько неудобны в декартовых координатах $x,y$, что имеет смысл перейти к совершенно другой системе координат — например, полярной. Но при замене переменных нельзя просто механически подставить новые выражения вместо старых: элемент площади $dA=dx\,dy$ тоже меняется, и меняется не произвольно, а по строгому правилу, которое учитывает, насколько новая система координат «растягивает» или «сжимает» пространство в окрестности каждой точки. Именно это растяжение и измеряет якобиан преобразования.

Представь декартову координатную сетку — квадратики размером $dx\times dy$. При переходе, скажем, к полярным координатам $(r,\theta)$ эти квадратики превращаются в маленькие искривлённые «кусочки кольца» — не квадраты, а фигуры, чья площадь зависит от того, насколько далеко от начала координат ты находишься (чем дальше от центра, тем больше «раздувается» кусочек при том же изменении угла $d\theta$). Якобиан — это в точности коэффициент, который переводит площадь маленького прямоугольника в новых координатах $dr\,d\theta$ в правильную площадь соответствующего кусочка на исходной плоскости $dx\,dy$.

Определение

Определение (якобиан преобразования). Пусть $x=x(u,v)$, $y=y(u,v)$ — гладкая замена переменных, переводящая область $D'$ в плоскости $(u,v)$ в область $D$ в плоскости $(x,y)$. Якобианом преобразования называется определитель матрицы частных производных:

$$J = \frac{\partial(x,y)}{\partial(u,v)} = \begin{vmatrix}\dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v}\\[2mm] \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v}\end{vmatrix} = \frac{\partial x}{\partial u}\frac{\partial y}{\partial v}-\frac{\partial x}{\partial v}\frac{\partial y}{\partial u}$$

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

$$\iint_D f(x,y)\,dx\,dy = \iint_{D'} f\big(x(u,v),y(u,v)\big)\,|J|\,du\,dv$$

где $|J|$ — модуль якобиана (площадь не может быть отрицательной, поэтому знак якобиана, отвечающий за ориентацию преобразования, отбрасывается).

Важнейший частный случай — переход к полярным координатам $x=r\cos\theta$, $y=r\sin\theta$, где $r\ge0$ — расстояние от начала координат, а $\theta$ — полярный угол. Вычислим якобиан этого преобразования явно:

Формула перехода к полярным координатам.

$$\iint_D f(x,y)\,dx\,dy = \iint_{D'} f(r\cos\theta,\,r\sin\theta)\,\,r\,dr\,d\theta$$

Множитель $r$ здесь — это в точности якобиан преобразования, и его ни в коем случае нельзя забывать: без него площадь кусочка сетки в полярных координатах вычислялась бы неверно.

Полярные координаты особенно удобны, когда область интегрирования — круг, кольцо или сектор (то есть когда естественные границы области описываются постоянными значениями $r$ и $\theta$, а не сложными кривыми в декартовых координатах), а также когда подынтегральная функция зависит от $x^2+y^2$ (поскольку в полярных координатах $x^2+y^2=r^2$ — простое и приятное упрощение).

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

Пример 1 (лёгкий, площадь круга через полярные координаты). Вычислить площадь круга радиуса $R$, используя переход к полярным координатам.

Площадь — это двойной интеграл функции $f\equiv1$ по кругу $D=\{x^2+y^2\le R^2\}$. В полярных координатах круг описывается предельно просто: $0\le r\le R$, $0\le\theta\le2\pi$ (никаких сложных зависимостей пределов друг от друга, в отличие от декартовых координат, где пришлось бы иметь дело с корнем $\sqrt{R^2-x^2}$).

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

Пример 2 (средний, объём тела вращения). Вычислить $\displaystyle\iint_D (x^2+y^2)\,dA$, где $D$ — круг $x^2+y^2\le4$ (это, например, момент инерции однородного диска относительно центральной оси при плотности $1$).

В полярных координатах $x^2+y^2=r^2$, область $D$: $0\le r\le2$, $0\le\theta\le2\pi$.

Внутренний интеграл: $\displaystyle\int_0^2 r^3\,dr=\left[\frac{r^4}4\right]_0^2=\frac{16}4=4$.

Внешний интеграл: $\displaystyle\int_0^{2\pi}4\,d\theta = 4\cdot2\pi=8\pi\approx25{,}13$.

Ответ: $\iint_D(x^2+y^2)\,dA=8\pi\approx25{,}13$. В декартовых координатах этот же интеграл потребовал бы возни с корнями $\sqrt{4-x^2}$ в пределах и тригонометрических подстановок при вычислении — переход к полярным координатам обошёл эту сложность полностью, заменив её на простое интегрирование степенной функции $r^3$.

Пример 3 (сложный, интеграл Гаусса через полярные координаты — прямая связь с нормальным распределением). Вычислить $\displaystyle I=\int_{-\infty}^{\infty} e^{-x^2}\,dx$ методом перехода к двойному интегралу и полярным координатам.

Этот интеграл не берётся напрямую (первообразная $e^{-x^2}$ не элементарна), но существует изящный обходной путь. Рассмотрим квадрат искомого интеграла $I^2$ и заметим, что его можно представить как произведение двух одинаковых интегралов с разными формальными именами переменных, а произведение интегралов, в свою очередь, переписать как один двойной интеграл по всей плоскости (по теореме Фубини, ведь подынтегральная функция распадается на произведение множителя от $x$ и множителя от $y$):

Ключевой шаг: переходим к полярным координатам по всей плоскости — область $\mathbb R^2$ описывается предельно просто как $0\le r<\infty$, $0\le\theta\le2\pi$, а подынтегральное выражение $x^2+y^2=r^2$ радикально упрощается:

Внутренний интеграл берём заменой $u=r^2$, $du=2r\,dr$:

Внешний интеграл: $\displaystyle\int_0^{2\pi}\frac12\,d\theta = \frac12\cdot2\pi=\pi$.

Значит, $I^2=\pi$, откуда (учитывая $I>0$, так как подынтегральная функция $e^{-x^2}$ всюду положительна):

Ответ: $I=\sqrt\pi\approx1{,}7725$. Это и есть тот самый знаменитый интеграл Гаусса — фундамент нормировки нормального распределения: стандартная нормальная плотность вероятности $\varphi(x)=\dfrac1{\sqrt{2\pi}}e^{-x^2/2}$ устроена именно так, чтобы после несложной замены переменной интегрировать её по всей числовой прямой давало ровно единицу, — и вычислить это можно исключительно благодаря трюку с переходом в двумерное пространство и полярными координатами, потому что напрямую, в одну переменную, этот интеграл не берётся никакими элементарными средствами.

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

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

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

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

Задание 1: Вычислить $\displaystyle\iint_D 1\,dA$, где $D=\{1\le x\le4,\ 2\le y\le5\}$ — прямоугольник.

Задание 2: Вычислить $\displaystyle\iint_D xy\,dA$ по прямоугольнику $D=\{0\le x\le2,\ 0\le y\le3\}$.

Задание 3: Вычислить $\displaystyle\iint_D (2x+3y)\,dA$ по прямоугольнику $D=\{0\le x\le1,\ 0\le y\le1\}$.

Задание 4: Изменить порядок интегрирования: $\displaystyle\int_0^2\int_0^{x} f(x,y)\,dy\,dx$.

Задание 5: Вычислить площадь области, ограниченной прямыми $y=0$, $x=2$ и $y=x$, представив её как двойной интеграл $\iint_D 1\,dA$.

Задание 6 (машинное обучение): Совместная плотность вероятности $f(x,y)=c$ (постоянная) задана на прямоугольнике $\{0\le x\le4,\ 0\le y\le5\}$. Найти константу $c$.

Задание 7: Вычислить $\displaystyle\iint_D e^{x+y}\,dA$ по прямоугольнику $D=\{0\le x\le1,\ 0\le y\le1\}$.

Задание 8: Найти площадь круга радиуса $3$ через переход к полярным координатам.

Задание 9: Вычислить якобиан преобразования $x=2u$, $y=3v$ (простое линейное растяжение осей).

Задание 10: Расставить пределы интегрирования для области $D$ — треугольника с вершинами $(0,0)$, $(1,0)$, $(0,1)$, как $x$-правильной области.

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

Задание 11: Вычислить $\displaystyle\iint_D xy\,dA$, где $D$ — треугольник с вершинами $(0,0)$, $(1,0)$, $(1,1)$.

Задание 12: Вычислить $\displaystyle\iint_D (x^2+y)\,dA$, где $D$ ограничена параболой $y=x^2$ и прямой $y=1$ (область под прямой и над параболой, $-1\le x\le1$).

Задание 13: Изменить порядок интегрирования: $\displaystyle\int_0^1\int_{\sqrt y}^1 f(x,y)\,dx\,dy$.

Задание 14 (машинное обучение): Совместная плотность $f(x,y)=xy$ на квадрате $\{0\le x\le\sqrt2,\ 0\le y\le\sqrt2\}$. Проверить, что это плотность (нормировка равна единице).

Задание 15: Вычислить $\displaystyle\iint_D (x^2+y^2)\,dA$, где $D$ — кольцо между окружностями радиусов $1$ и $2$ с центром в начале координат.

Задание 16: Вычислить якобиан преобразования к полярным координатам $x=r\cos\theta$, $y=r\sin\theta$ в общем виде и проверить его значение при $r=3$.

Задание 17: Вычислить $\displaystyle\iint_D e^{-(x^2+y^2)}\,dA$, где $D$ — круг радиуса $2$ с центром в начале координат.

Задание 18: Вычислить площадь области, ограниченной кривыми $y=x^2$ и $y=2-x^2$.

Задание 19 (машинное обучение): Совместная плотность двух признаков $f(x,y)=\dfrac{6}{5}(x+y^2)$ на прямоугольнике $\{0\le x\le1,\ 0\le y\le1\}$. Найти вероятность $P(X\le0{,}5)$ (то есть $\iint$ по подобласти $0\le x\le0{,}5$, $0\le y\le1$).

Задание 20: Изменить порядок интегрирования: $\displaystyle\int_0^1\int_0^{1-x^2} f(x,y)\,dy\,dx$.

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

Задание 21: Вычислить $\displaystyle\iint_D \frac{1}{(1+x^2+y^2)^2}\,dA$, где $D$ — вся плоскость $\mathbb R^2$ (несобственный интеграл).

Задание 22: Вычислить объём тела, ограниченного сверху параболоидом $z=4-x^2-y^2$ и снизу плоскостью $z=0$ (весь параболоид над плоскостью).

Задание 23: Вычислить $\displaystyle\iint_D xy\,dA$, где $D$ — область в первой четверти, ограниченная кривой $x^2+y^2=1$ (то есть четверть круга единичного радиуса).

Задание 24 (машинное обучение, двумерное нормальное распределение): Стандартная двумерная нормальная плотность (независимые компоненты) имеет вид $f(x,y)=\dfrac1{2\pi}e^{-(x^2+y^2)/2}$. Проверить нормировку — что $\iint_{\mathbb R^2} f(x,y)\,dA=1$.

Задание 25: Вычислить $\displaystyle\iint_D \sin(x^2+y^2)\,dA$, где $D$ — круг радиуса $\sqrt\pi$ с центром в начале координат.

Задание 26: Восстановить область по повторному интегралу $\displaystyle\int_{-1}^1\int_{-\sqrt{1-y^2}}^{\sqrt{1-y^2}} f(x,y)\,dx\,dy$ и переписать его как двойной интеграл по кругу в полярных координатах.

Задание 27 (продвинутое, замена переменных общего вида): Найти якобиан преобразования $x=u^2-v^2$, $y=2uv$ (комплексное отображение $z=w^2$, где $z=x+iy$, $w=u+iv$).

Задание 28: Вычислить объём тела, ограниченного цилиндром $x^2+y^2=9$ и плоскостями $z=0$, $z=x+5$ (при условии, что $x+5\ge0$ на всей области, то есть $x\ge-5$, что автоматически выполнено внутри круга радиуса $3$).

Задание 29 (машинное обучение, корреляция и эллиптическая область): Совместная плотность двух коррелированных признаков задана как $f(x,y)=\dfrac{2}{\pi}$ (постоянная) на эллиптической области $D=\{x^2+\dfrac{y^2}4\le1\}$. Проверить нормировку, используя замену переменных $x=u$, $y=2v$ (переводящую эллипс в круг единичного радиуса).

Задание 30 (продвинутое, обобщённый интеграл Гаусса): Используя результат интеграла Гаусса $\int_{-\infty}^{\infty}e^{-x^2}dx=\sqrt\pi$, вычислить $\displaystyle\int_{-\infty}^{\infty}e^{-ax^2}\,dx$ при $a>0$ через замену переменной (без повторного перехода к двойному интегралу).

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

Двойные интегралы — тема, где технически всё довольно просто (два обычных интеграла подряд), но концептуальные ошибки в постановке задачи встречаются очень часто. Вот самые типичные из них.

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

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

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

• Использование пределов внешнего интеграла, зависящих от переменной. Если во внешнем интеграле у тебя оказались пределы вида $\int_{f(y)}^{g(y)}$ вместо постоянных чисел — это верный признак того, что порядок интегрирования расставлен неверно: внешний интеграл суммирует «полосы» или «слои» всей области, и диапазон переменной внешнего интегрирования обязан охватывать весь диапазон значений этой переменной по всей области целиком, а значит быть просто числом, а не функцией.

• Игнорирование смены знака области при пересечении границ. Если границы области $g_1(x)$ и $g_2(x)$ где-то пересекаются внутри заявленного диапазона $x$ (то есть одна кривая «перескакивает» через другую), это признак того, что область на самом деле не является простой $x$-правильной областью на всём диапазоне и её нужно разбивать на несколько частей с отдельными пределами для каждой части.

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

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

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

• Двойной интеграл $\iint_D f(x,y)\,dA$ определяется как предел интегральных сумм по разбиению области $D$ — прямое обобщение однократного интеграла на плоский случай.

• Если $f(x,y)\ge0$, двойной интеграл равен объёму тела под поверхностью $z=f(x,y)$ над областью $D$; при $f\equiv1$ он превращается в площадь области $D$.

• Теорема Фубини сводит двойной интеграл к повторному: сначала интегрируем по одной переменной при фиксированной другой (внутренний интеграл), потом по оставшейся переменной (внешний интеграл).

• Область называется $x$-правильной, если для каждого $x$ из диапазона $[a,b]$ переменная $y$ ограничена снизу и сверху функциями $g_1(x)$ и $g_2(x)$; аналогично определяется $y$-правильная область.

• Пределы внешнего интеграла всегда постоянные числа; пределы внутреннего интеграла, как правило, зависят от переменной внешнего интегрирования — перед расстановкой пределов область обязательно нужно нарисовать.

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

• При замене переменных $x=x(u,v)$, $y=y(u,v)$ элемент площади меняется по формуле $dx\,dy=|J|\,du\,dv$, где $J$ — якобиан преобразования (определитель матрицы частных производных).

• Переход к полярным координатам $x=r\cos\theta$, $y=r\sin\theta$ даёт якобиан $J=r$, поэтому $dx\,dy=r\,dr\,d\theta$ — этот множитель $r$ никогда нельзя забывать.

• Двойной интеграл от совместной плотности вероятности $f(x,y)$ по области $D$ даёт вероятность попадания случайной пары $(X,Y)$ в эту область; интеграл по всей плоскости обязан равняться единице (условие нормировки).

• Интеграл Гаусса $\int_{-\infty}^{\infty}e^{-x^2}dx=\sqrt\pi$ вычисляется через возведение в квадрат, переход к двойному интегралу по всей плоскости и последующий переход к полярным координатам — этот результат лежит в основе нормировки нормального распределения.

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

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

Якобиан преобразования, который мы ввели сегодня как определитель матрицы частных производных, — это первое серьёзное практическое применение определителей матриц $2\times2$ в анализе; в следующем уроке про тройные интегралы этот же принцип обобщается на матрицу $3\times3$ частных производных при переходе к сферическим или цилиндрическим координатам. А связь с теорией вероятностей, которую мы разбирали через совместные плотности, станет ещё более явной в дальнейших курсах статистики и машинного обучения: нормировка многомерных плотностей, вычисление математических ожиданий функций от нескольких случайных величин и байесовский вывод в многомерных моделях — всё это в конечном счёте кратные интегралы ровно того типа, который мы разобрали сегодня.

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

• Интеграл Гаусса $\int_{-\infty}^{\infty}e^{-x^2}dx=\sqrt\pi$ невозможно вычислить, оставаясь в рамках одной переменной, — первообразная функции $e^{-x^2}$ доказуемо не выражается через элементарные функции. Тем не менее сам определённый интеграл по всей числовой прямой вычисляется элегантно и точно именно благодаря искусственному переходу в двумерное пространство и последующему переходу к полярным координатам — редкий и красивый пример того, как повышение размерности задачи упрощает, а не усложняет её решение.

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

• Плотность стандартного нормального распределения $\varphi(x)=\dfrac1{\sqrt{2\pi}}e^{-x^2/2}$ обязана своим знаменателем $\sqrt{2\pi}$ именно интегралу Гаусса: без строгого доказательства значения этого интеграла нельзя было бы утверждать, что площадь под колоколообразной кривой нормального распределения действительно равна единице, а значит, само распределение вообще было бы некорректно определено как плотность вероятности.

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

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

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

• Всегда начинай с чертежа области интегрирования, даже если задача кажется простой. Большинство ошибок с пределами интегрирования происходят именно от попытки расставить пределы «в уме», без визуализации формы области.

• Если подынтегральная функция явно зависит от $x^2+y^2$, а область — круг, кольцо или сектор, почти наверняка стоит сразу переходить к полярным координатам: это, как правило, кардинально упрощает и пределы интегрирования, и саму подынтегральную функцию.

• Если внутренний интеграл в исходном порядке интегрирования не берётся элементарно (особенно если видишь конструкции вроде $e^{y^2}$ или $\sin(y^2)/y$ без явного множителя $y$ или $y\,dy$ рядом), сразу пробуй сменить порядок интегрирования — очень часто это превращает неберущийся интеграл в элементарный.

• Если подынтегральная функция раскладывается в произведение множителя, зависящего только от $x$, и множителя, зависящего только от $y$ (например, $f(x,y)=g(x)h(y)$), а область — прямоугольник, двойной интеграл сразу распадается на произведение двух обычных однократных интегралов — не нужно даже формально сводить его к повторному в общем виде.

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

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

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

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