6. Двойные интегралы

Интеграл от функции двух переменных по плоской области: объём под поверхностью, площадь, масса пластинки. Главный навык главы — по рисунку области расставить пределы в повторном интеграле, поменять порядок интегрирования и вовремя перейти к полярным координатам.

КР2, задание 1 (пример в РПД)ДЗ2, задачи 1-3практика 1 (пример в РПД)РГР1, задание 1вопросы экзамена №14-20, 23-24

Откуда известно, что это нужно

В РПД курса двойной интеграл стоит в разделе 2 (темы 2.1-2.5 и 2.7: мера Жордана, двойной интеграл, криволинейные координаты, замена переменных, несобственные интегралы, приложения). В примерах из РПД: задание 1 КР2 — площадь фигуры двойным интегралом с проверкой сменой порядка; ДЗ2 — вычисление двойного интеграла, смена порядка, обобщённые полярные координаты; пример работы на практике 1 — площади и интегралы в полярных координатах; задание 1 типового расчёта — площадь области двойным интегралом с рисунком. В примере КР1 из РПД есть задание на переход к полярным координатам с рисунком. Все эти задачи разобраны ниже.

Что нужно помнить из прошлых семестров: определённый интеграл и его геометрический смысл (площадь под графиком), формула Ньютона-Лейбница, табличные интегралы, замена переменной и интегрирование по частям, $\int\cos^2 t\,dt = \frac t2 + \frac{\sin 2t}{4}$, $\int_0^{\pi/2}\cos^4 t\,dt = \frac{3\pi}{16}$. Полярные координаты: $x = r\cos\varphi$, $y = r\sin\varphi$, $x^2+y^2 = r^2$.

1. Идея: объём, масса, интегральные суммы

вопрос экзамена №14вопрос экзамена №15

Две задачи, которые приводят к двойному интегралу

Объём. Над областью $D$ на плоскости $xOy$ лежит поверхность $z = f(x,y) \ge 0$. Нужен объём тела между $D$ и поверхностью (такое тело называют цилиндрическим брусом). Делаем как с площадью под графиком в первом семестре: режем $D$ на маленькие кусочки $D_i$ площади $\Delta S_i$, в каждом берём точку $P_i$, ставим на кусочек столбик высоты $f(P_i)$. Объём столбика $f(P_i)\,\Delta S_i$, объём всего тела примерно равен сумме столбиков. Чем мельче кусочки, тем точнее.

Масса пластинки. Тонкая пластинка имеет форму $D$, плотность $\rho(x,y)$ (масса на единицу площади) меняется от точки к точке. На маленьком кусочке плотность почти постоянна, поэтому масса кусочка $\approx \rho(P_i)\,\Delta S_i$, а вся масса $\approx \sum \rho(P_i)\,\Delta S_i$. Та же самая сумма.

Ниже поверхность $z = 3 - \frac{x^2+y^2}{4}$ над квадратом $[0,2]\times[0,2]$ (полупрозрачная) и столбики интегральной суммы с точкой $P_i$ в центре клетки. Ползунок $k$ — сколько клеток по каждой стороне. При росте $k$ ступеньки прижимаются к поверхности.

Определение (двойной интеграл)

Пусть $D$ — замкнутая ограниченная область, у которой есть площадь (см. ниже про меру Жордана), $f(x,y)$ задана на $D$.

  1. Разбиение $T$: делим $D$ на части $D_1,\dots,D_n$, у которых есть площади $\Delta S_i$ и которые не перекрываются (общими могут быть только границы).
  2. Диаметр части $d(D_i)$ — наибольшее расстояние между её точками; мелкость разбиения $\lambda(T) = \max_i d(D_i)$.
  3. В каждой части выбираем точку $P_i(\xi_i,\eta_i) \in D_i$ и составляем интегральную сумму $\sigma = \sum_{i=1}^n f(\xi_i,\eta_i)\,\Delta S_i$.

Если существует конечный предел $\sigma$ при $\lambda(T)\to 0$, не зависящий ни от способа разбиения, ни от выбора точек, он называется двойным интегралом:

$$\iint_D f(x,y)\,dx\,dy = \lim_{\lambda\to 0}\sum_{i=1}^n f(\xi_i,\eta_i)\,\Delta S_i,$$

а $f$ — интегрируемой на $D$. Обозначения $dS$ и $dx\,dy$ — одно и то же (элемент площади).

Как понимать

Мера Жордана: что значит «у области есть площадь»

Для многоугольника площадь понятна. Для области, ограниченной кривыми, её надо определить. Идея Жордана: накрыть плоскость сеткой квадратов со стороной $h$ и посчитать два числа:

Всегда $s_h \le S_h$. Разница $S_h - s_h$ — это площадь квадратов, задевающих границу. На картинке они закрашены; уменьши $h$ и посмотри, как полоска сужается.

Определение (мера Жордана)

Ограниченное множество $D$ называется измеримым по Жордану (квадрируемым), если $\lim\limits_{h\to 0} s_h = \lim\limits_{h\to 0} S_h$. Это общее значение и есть площадь (мера) $D$. Числа $\sup s_h$ и $\inf S_h$ называют внутренней и внешней мерой.

Критерий: $D$ измеримо тогда и только тогда, когда его граница имеет меру нуль (её можно накрыть клетками сколь угодно малой суммарной площади).

Что из этого нужно на практике

Суммы Дарбу и существование интеграла

Определение через предел сумм неудобно проверять: надо перебрать все разбиения и все точки. Суммы Дарбу убирают выбор точек: берём в каждом куске самое маленькое и самое большое значение.

Определение (суммы Дарбу)

Пусть $f$ ограничена на $D$, $T$ — разбиение, $m_i = \inf\limits_{D_i} f$, $M_i = \sup\limits_{D_i} f$. Нижняя и верхняя суммы Дарбу:

$$s_T = \sum_{i=1}^n m_i\,\Delta S_i, \qquad S_T = \sum_{i=1}^n M_i\,\Delta S_i.$$

Для любого выбора точек $s_T \le \sigma \le S_T$. При измельчении разбиения $s_T$ не убывает, $S_T$ не возрастает, и любая нижняя сумма не больше любой верхней.

Теорема (критерий интегрируемости)

Ограниченная на измеримой области $D$ функция $f$ интегрируема тогда и только тогда, когда $\lim\limits_{\lambda(T)\to 0}(S_T - s_T) = 0$, то есть $\lim\limits_{\lambda\to0} \sum (M_i - m_i)\,\Delta S_i = 0$.

Достаточные условия. Интегрируемы:

Как понимать

$S_T - s_T = \sum (M_i - m_i)\Delta S_i$ — суммарный объём «зазоров» между столбиками, построенными по максимуму и по минимуму. Для непрерывной функции на мелком куске $M_i - m_i$ мал сразу для всех кусков (равномерная непрерывность на замкнутой ограниченной области), поэтому зазор $\le \varepsilon\cdot S(D) \to 0$. Это и есть доказательство первого достаточного условия в одну строку.

Чему равен $\iint_D 1\,dx\,dy$, если $D$ — круг радиуса 3?
Интеграл от единицы равен площади области: $\pi\cdot 3^2 = 9\pi$. $6\pi$ — длина окружности.

2. Свойства двойного интеграла, теорема о среднем

вопрос экзамена №15вопрос экзамена №16

Свойства те же, что у определённого интеграла, и доказываются одинаково: переходом к пределу в интегральных суммах. Все функции ниже интегрируемы на $D$.

Свойства
  1. Площадь: $\iint_D dx\,dy = S(D)$.
  2. Линейность: $\iint_D (\alpha f + \beta g)\,dS = \alpha\iint_D f\,dS + \beta\iint_D g\,dS$.
  3. Аддитивность по области: если $D = D_1\cup D_2$, причём $D_1$ и $D_2$ пересекаются только по границе, то $\iint_D f\,dS = \iint_{D_1} f\,dS + \iint_{D_2} f\,dS$.
  4. Монотонность: если $f \le g$ на $D$, то $\iint_D f\,dS \le \iint_D g\,dS$. В частности, $f\ge 0 \Rightarrow \iint_D f\,dS \ge 0$.
  5. Модуль: $|f|$ тоже интегрируема и $\left|\iint_D f\,dS\right| \le \iint_D |f|\,dS$. Произведение интегрируемых функций интегрируемо.
  6. Оценка: если $m \le f \le M$ на $D$, то $m\,S(D) \le \iint_D f\,dS \le M\,S(D)$.
  7. Значения $f$ на множестве меры нуль (например, на кривой) не влияют на интеграл.
Теорема о среднем

Если $f$ непрерывна на замкнутой ограниченной связной измеримой области $D$, то найдётся точка $P_0(x_0,y_0)\in D$, для которой

$$\iint_D f(x,y)\,dx\,dy = f(x_0,y_0)\cdot S(D).$$

Число $\bar f = \frac{1}{S(D)}\iint_D f\,dS$ называют средним значением $f$ на $D$.

Доказательство. По теореме Вейерштрасса $f$ достигает на $D$ минимума $m$ и максимума $M$. Из оценки $m \le \bar f \le M$. Непрерывная функция на связной области принимает все промежуточные значения, значит, и значение $\bar f$ в некоторой точке $P_0$.

Как понимать

Тело под поверхностью можно «разровнять» в цилиндр с тем же основанием $D$ и той же высотой $\bar f$; объём не изменится. Теорема говорит, что такая высота реально достигается в какой-то точке области.

Пример 1 (оценка интеграла)

Оцените $I = \iint_D (x^2 + 4y^2 + 9)\,dx\,dy$, где $D$: $x^2+y^2 \le 4$, не вычисляя его. Затем вычислите точно (после раздела 5) и найдите среднее значение функции.

Решение

Площадь $S(D) = 4\pi$. Ищем наименьшее и наибольшее значения $f = x^2+4y^2+9$ в круге. Наименьшее: в начале координат, $f(0,0) = 9$. Наибольшее: $f = (x^2+y^2) + 3y^2 + 9 \le 4 + 3\cdot 4 + 9 = 25$, достигается в точках $(0,\pm 2)$. По свойству оценки

$$9\cdot 4\pi \le I \le 25\cdot 4\pi, \qquad 36\pi \le I \le 100\pi.$$

Точно (полярные координаты, $dx\,dy = r\,dr\,d\varphi$):

$$I = \int_0^{2\pi}\!\!d\varphi\int_0^2 \bigl(r^2\cos^2\varphi + 4r^2\sin^2\varphi + 9\bigr)r\,dr = \int_0^{2\pi}\bigl(4\cos^2\varphi + 16\sin^2\varphi + 18\bigr)d\varphi = 4\pi + 16\pi + 36\pi = 56\pi.$$

Здесь $\int_0^2 r^3 dr = 4$, $\int_0^2 9r\,dr = 18$, $\int_0^{2\pi}\cos^2\varphi\,d\varphi = \int_0^{2\pi}\sin^2\varphi\,d\varphi = \pi$. Оценка выполнена. Среднее значение $\bar f = 56\pi/4\pi = 14$.

Ответ: $36\pi \le I \le 100\pi$; $I = 56\pi$, $\bar f = 14$.

Пример 2 (среднее значение)

Найдите среднее значение $f(x,y) = x^2 + y$ на прямоугольнике $D = [0,2]\times[0,1]$ и укажите точку, где оно достигается.

Решение

Прямоугольник — простейшая область: пределы постоянные (раздел 3).

$$\iint_D (x^2+y)\,dx\,dy = \int_0^2 dx\int_0^1 (x^2+y)\,dy = \int_0^2\Bigl(x^2 + \frac12\Bigr)dx = \frac83 + 1 = \frac{11}{3}.$$

$S(D) = 2$, поэтому $\bar f = \frac{11}{6}$. Точка: нужно $x^2 + y = \frac{11}{6}$, например $x = 1$, $y = \frac56$ — она лежит в $D$.

Ответ: $\bar f = \frac{11}{6}$, достигается, например, в точке $\left(1, \frac56\right)$ (такая точка не единственна).

3. Вычисление: сведение к повторному интегралу

вопрос экзамена №16главный навык главы

Прямоугольник

Объём тела можно считать «ломтиками». Зафиксируем $x$ и разрежем тело плоскостью, параллельной $yOz$. В сечении получится криволинейная трапеция площади $A(x) = \int_c^d f(x,y)\,dy$. Объём — сумма ломтиков толщины $dx$, то есть $\int_a^b A(x)\,dx$.

Теорема (о сведении к повторному, прямоугольник)

Пусть $f$ интегрируема на $D = [a,b]\times[c,d]$ и при каждом $x\in[a,b]$ существует $\int_c^d f(x,y)\,dy$. Тогда

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

Для непрерывной $f$ можно и в другом порядке: $\int_c^d dy\int_a^b f\,dx$ — ответ тот же.

Запись $\int_a^b dx\int_c^d f\,dy$ читается справа налево: сначала внутренний интеграл по $y$ (при этом $x$ — константа), получаем функцию от $x$, потом её интегрируем по $x$.

Пример 3

Вычислите $\iint_D (x^2 + xy)\,dx\,dy$, $D = [0,1]\times[0,2]$.

Решение

Внутренний интеграл по $y$ от 0 до 2, $x$ — константа:

$$\int_0^2 (x^2 + xy)\,dy = \Bigl(x^2 y + \frac{x y^2}{2}\Bigr)\Big|_{y=0}^{y=2} = 2x^2 + 2x.$$

Внешний: $\int_0^1 (2x^2 + 2x)\,dx = \frac23 + 1 = \frac53$.

Проверка в другом порядке: $\int_0^1 (x^2 + xy)\,dx = \frac13 + \frac y2$, $\int_0^2\bigl(\frac13 + \frac y2\bigr)dy = \frac23 + 1 = \frac53$.

Ответ: $\frac53$.

Области I и II типа

Определение

Область I типа (элементарная относительно оси $Oy$, «простая в направлении $Oy$»):

$$D = \{(x,y):\ a \le x \le b,\ \ \varphi_1(x) \le y \le \varphi_2(x)\},$$

где $\varphi_1, \varphi_2$ непрерывны. Любая вертикальная прямая $x = x_0$, $a\lt x_0\lt b$, пересекает $D$ по одному отрезку: входит через нижнюю границу $y = \varphi_1(x)$, выходит через верхнюю $y = \varphi_2(x)$.

Область II типа (элементарная относительно оси $Ox$):

$$D = \{(x,y):\ c \le y \le d,\ \ \psi_1(y) \le x \le \psi_2(y)\}.$$

Любая горизонтальная прямая пересекает $D$ по одному отрезку: вход через левую границу $x = \psi_1(y)$, выход через правую $x = \psi_2(y)$.

Теорема (формулы сведения)

Если $f$ непрерывна на $D$, то

$$\text{I тип:}\quad \iint_D f\,dx\,dy = \int_a^b dx\int_{\varphi_1(x)}^{\varphi_2(x)} f(x,y)\,dy, \qquad \text{II тип:}\quad \iint_D f\,dx\,dy = \int_c^d dy\int_{\psi_1(y)}^{\psi_2(y)} f(x,y)\,dx.$$

Двигай ползунок $p$: вертикальный отрезок при $x = p$ всегда входит в область через нижнюю кривую и выходит через верхнюю. Пределы внутреннего интеграла — уравнения этих кривых, решённые относительно $y$.

Алгоритм расстановки пределов
  1. Нарисуй область. Каждую линию из условия построй по отдельности, найди точки пересечения (решая системы уравнений). Без рисунка пределы почти всегда ставят неправильно.
  2. Выбери порядок. Смотри, какой тип у области: если вертикальные прямые входят и выходят всегда через одни и те же кривые — I тип, внешний интеграл по $x$. Если так с горизонтальными — II тип, внешний по $y$. Если годятся оба — бери тот, где проще интегралы.
  3. Внешние пределы — числа: крайние значения $x$ (для I типа) или $y$ (для II типа) на области. Обычно это абсциссы или ординаты точек пересечения.
  4. Внутренние пределы — функции: проведи прямую $x = \mathrm{const}$ (или $y = \mathrm{const}$) через область. Нижний предел — кривая, через которую входишь, верхний — через которую выходишь. Уравнения кривых реши относительно внутренней переменной.
  5. Если входная или выходная кривая меняется (у границы есть «излом» по ходу прямой), разбей область вертикальной (горизонтальной) прямой через точку излома и сложи интегралы.
  6. Проверка: внешние пределы — константы, внутренние зависят только от внешней переменной, нижний предел $\le$ верхнего на всём промежутке. Подставь одно промежуточное значение и проверь по рисунку.
Ловушка

Внешние пределы никогда не зависят от переменных: $\int_0^{x} dx\ldots$ — бессмыслица. После вычисления двойного интеграла должно получиться число (или выражение от параметров $a$, $R$ и т.п.), а не функция от $x$ или $y$.

Пример 4

Вычислите $\iint_D xy\,dx\,dy$, где $D$ — треугольник $x\ge 0$, $y \ge 0$, $x + y \le 1$.

Решение

$x$ меняется от 0 до 1. При фиксированном $x$ вертикальная прямая входит в треугольник на оси $Ox$ ($y = 0$) и выходит на гипотенузе $y = 1 - x$. Это область I типа:

$$\iint_D xy\,dx\,dy = \int_0^1 dx\int_0^{1-x} xy\,dy = \int_0^1 x\cdot\frac{(1-x)^2}{2}\,dx = \frac12\int_0^1 (x - 2x^2 + x^3)\,dx = \frac12\Bigl(\frac12 - \frac23 + \frac14\Bigr) = \frac{1}{24}.$$

Ответ: $\frac{1}{24}$.

Пример 5

Вычислите $\iint_D x^2 y\,dx\,dy$, где $D$ ограничена линиями $y = x^2$ и $y = \sqrt x$. Сделайте это в обоих порядках.

Решение

Точки пересечения: $x^2 = \sqrt x \Rightarrow x^4 = x \Rightarrow x = 0$ или $x = 1$. На $(0,1)$ выше лежит $\sqrt x$ (проверка при $x = 0{,}25$: $\sqrt x = 0{,}5 \gt x^2 = 0{,}0625$).

I тип (внешний по $x$):

$$\int_0^1 dx\int_{x^2}^{\sqrt x} x^2 y\,dy = \int_0^1 x^2\cdot\frac{x - x^4}{2}\,dx = \frac12\Bigl(\frac14 - \frac17\Bigr) = \frac{3}{56}.$$

II тип (внешний по $y$): решаем уравнения относительно $x$: $y = \sqrt x \Leftrightarrow x = y^2$ (это левая граница), $y = x^2 \Leftrightarrow x = \sqrt y$ (правая). $y$ от 0 до 1:

$$\int_0^1 dy\int_{y^2}^{\sqrt y} x^2 y\,dx = \int_0^1 y\cdot\frac{y^{3/2} - y^6}{3}\,dy = \frac13\Bigl(\frac{2}{7} - \frac18\Bigr) = \frac13\cdot\frac{9}{56} = \frac{3}{56}.$$

Здесь $y\cdot y^{3/2} = y^{5/2}$ и $\int_0^1 y^{5/2}dy = \frac27$. Ответы совпали.

Ответ: $\frac{3}{56}$.

Пример 6 из примера ДЗ2 в РПД

Вычислите $\iint_G (4 - y)\,dx\,dy$, где область $G$ ограничена кривыми $x^2 = 4y$, $y = 1$, $x = 0$ ($x \gt 0$).

Решение

Парабола $y = \frac{x^2}{4}$, прямая $y = 1$ пересекает её при $x = 2$ (берём $x \gt 0$). Область — криволинейный треугольник с вершинами $(0,0)$, $(0,1)$, $(2,1)$.

Удобнее II тип: при фиксированном $y\in[0,1]$ горизонтальная прямая входит на оси $Oy$ ($x = 0$) и выходит на параболе $x = 2\sqrt y$. Подынтегральная функция зависит только от $y$, поэтому внутренний интеграл тривиален:

$$\int_0^1 dy\int_0^{2\sqrt y}(4 - y)\,dx = \int_0^1 (4-y)\,2\sqrt y\,dy = \int_0^1 \bigl(8y^{1/2} - 2y^{3/2}\bigr)dy = 8\cdot\frac23 - 2\cdot\frac25 = \frac{16}{3} - \frac45 = \frac{68}{15}.$$

Проверка I типом: $x\in[0,2]$, $y$ от $\frac{x^2}{4}$ до 1:

$$\int_0^2 dx\int_{x^2/4}^{1}(4-y)\,dy = \int_0^2\Bigl(\frac72 - x^2 + \frac{x^4}{32}\Bigr)dx = 7 - \frac83 + \frac{1}{5} = \frac{68}{15}.$$

Ответ: $\frac{68}{15}$.

Пример 7 (один порядок удобнее другого)

Вычислите $\iint_D y\,dx\,dy$, где $D$ — треугольник, ограниченный прямыми $y = x$, $y = 2 - x$, $y = 0$.

Решение

Вершины: $(0,0)$, $(2,0)$ и пересечение $x = 2 - x$, то есть $(1,1)$.

II тип — один интеграл. $y\in[0,1]$; горизонтальная прямая входит на прямой $x = y$ и выходит на прямой $x = 2 - y$:

$$\int_0^1 dy\int_y^{2-y} y\,dx = \int_0^1 y(2 - 2y)\,dy = 1 - \frac23 = \frac13.$$

I тип — два интеграла. Вертикальная прямая всегда входит на $y = 0$, а выходит при $x\lt 1$ через $y = x$, при $x \gt 1$ через $y = 2 - x$. Верхняя граница меняется в точке $x = 1$, поэтому разбиваем:

$$\int_0^1 dx\int_0^x y\,dy + \int_1^2 dx\int_0^{2-x} y\,dy = \int_0^1\frac{x^2}{2}dx + \int_1^2\frac{(2-x)^2}{2}dx = \frac16 + \frac16 = \frac13.$$

Ответ: $\frac13$. Вывод: выбирай порядок, при котором область не приходится резать.

Пример 8 (парабола «на боку»)

Вычислите $\iint_D x\,dx\,dy$, где $D$ ограничена параболой $x = y^2$ и прямой $x = y + 2$.

Решение

Пересечения: $y^2 = y + 2 \Rightarrow y^2 - y - 2 = 0 \Rightarrow y = -1$ или $y = 2$. Точки $(1,-1)$ и $(4,2)$.

II тип (удобно). $y\in[-1,2]$, слева парабола $x = y^2$, справа прямая $x = y + 2$ (при $y = 0$: $0 \le x \le 2$ — верно):

$$\int_{-1}^{2}dy\int_{y^2}^{y+2}x\,dx = \frac12\int_{-1}^{2}\bigl((y+2)^2 - y^4\bigr)dy = \frac12\Bigl(\frac{(y+2)^3}{3} - \frac{y^5}{5}\Bigr)\Big|_{-1}^{2} = \frac12\Bigl(\frac{64 - 1}{3} - \frac{32 + 1}{5}\Bigr) = \frac12\Bigl(21 - \frac{33}{5}\Bigr) = \frac{36}{5}.$$

I тип (для сравнения). Решаем относительно $y$: парабола даёт две ветви $y = \pm\sqrt x$, прямая $y = x - 2$. При $x\in[0,1]$ вертикаль входит через $y = -\sqrt x$, при $x\in[1,4]$ — через прямую $y = x - 2$; выходит всегда через $y = \sqrt x$:

$$\int_0^1 dx\int_{-\sqrt x}^{\sqrt x} x\,dy + \int_1^4 dx\int_{x-2}^{\sqrt x} x\,dy = \int_0^1 2x^{3/2}dx + \int_1^4\bigl(x^{3/2} - x^2 + 2x\bigr)dx = \frac45 + \Bigl(\frac{62}{5} - 21 + 15\Bigr) = \frac{36}{5}.$$

Ответ: $\frac{36}{5}$.

Пример 9 (пределы в обоих порядках)

Область $D$ ограничена параболой $y = 2 - x^2$ и прямой $y = x$. Расставьте пределы в $\iint_D f\,dx\,dy$ в обоих порядках и вычислите $\iint_D x\,dx\,dy$.

Решение

Пересечения: $2 - x^2 = x \Rightarrow x^2 + x - 2 = 0 \Rightarrow x = -2,\ x = 1$. Точки $(-2,-2)$ и $(1,1)$. Вершина параболы $(0,2)$.

I тип: $x\in[-2,1]$, снизу прямая, сверху парабола — один интеграл:

$$\iint_D f\,dx\,dy = \int_{-2}^{1}dx\int_{x}^{2-x^2}f(x,y)\,dy.$$

II тип: парабола относительно $x$: $x = \pm\sqrt{2 - y}$. При $y\in[-2,1]$ горизонталь входит через левую ветвь $x = -\sqrt{2-y}$ и выходит через прямую $x = y$. При $y\in[1,2]$ прямую она уже не пересекает: выходит через правую ветвь $x = \sqrt{2-y}$. Значит,

$$\iint_D f\,dx\,dy = \int_{-2}^{1}dy\int_{-\sqrt{2-y}}^{y}f\,dx + \int_{1}^{2}dy\int_{-\sqrt{2-y}}^{\sqrt{2-y}}f\,dx.$$

Вычисляем I типом: $\int_{-2}^1 x(2 - x^2 - x)\,dx = \Bigl(x^2 - \frac{x^4}{4} - \frac{x^3}{3}\Bigr)\Big|_{-2}^{1} = \Bigl(1 - \frac14 - \frac13\Bigr) - \Bigl(4 - 4 + \frac83\Bigr) = \frac{5}{12} - \frac{32}{12} = -\frac94$.

Отрицательный ответ нормален: большая часть области лежит при $x \lt 0$. Площадь (при $f = 1$) получилась бы $\frac92$ в обоих порядках.

Ответ: $-\frac94$.

Как выглядит тело из примера 4: поверхность $z = xy$ над треугольником. Двойной интеграл $\frac1{24}$ — его объём.

Область между $y = x^2$ и $y = 4$. Какая запись верна?
$x$ меняется от $-2$ до $2$ (пересечение $x^2 = 4$), снизу парабола, сверху прямая. Во втором варианте пределы внутреннего интеграла по $x$ зависят от $x$ — так нельзя. В четвёртом потеряна левая половина: правильно $\int_0^4 dy\int_{-\sqrt y}^{\sqrt y} f\,dx$.

4. Смена порядка интегрирования

ДЗ2, задача 2КР2, задание 1 (проверка)

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

Алгоритм смены порядка
  1. По пределам выпиши неравенства, задающие область. Для $\int_a^b dx\int_{\varphi_1(x)}^{\varphi_2(x)} f\,dy$ это $a\le x\le b$, $\varphi_1(x)\le y\le\varphi_2(x)$. Если интегралов несколько — область каждого.
  2. Нарисуй область: линии $x = a$, $x = b$, $y = \varphi_1(x)$, $y = \varphi_2(x)$. Отметь угловые точки.
  3. Найди наименьшее и наибольшее значение новой внешней переменной на области — это внешние пределы.
  4. Реши уравнения границ относительно новой внутренней переменной. Проведи прямые нового направления: через какую кривую входят, через какую выходят. Если кривые входа/выхода меняются, разбей интеграл на части в угловых точках.
  5. Проверка: посчитай площадь (интеграл с $f = 1$) в обоих порядках — числа должны совпасть.
Пример 10

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

Решение

Область: $0 \le x\le 1$, $x^2\le y\le x$ — между параболой $y = x^2$ (снизу) и прямой $y = x$ (сверху), они пересекаются в $(0,0)$ и $(1,1)$.

Новый внешний — $y\in[0,1]$. Относительно $x$: прямая $x = y$, парабола $x = \sqrt y$. Горизонталь входит через прямую (она левее: при $y = 0{,}25$ прямая даёт $x = 0{,}25$, парабола $x = 0{,}5$), выходит через параболу:

$$\int_0^1 dx\int_{x^2}^{x} f\,dy = \int_0^1 dy\int_{y}^{\sqrt y} f\,dx.$$

Проверка на $f = x + y$: слева $\int_0^1\bigl(x(x - x^2) + \frac{x^2 - x^4}{2}\bigr)dx = \frac{3}{20}$, справа $\int_0^1\bigl(\frac{y - y^2}{2} + y(\sqrt y - y)\bigr)dy = \frac{3}{20}$.

Ответ: $\int_0^1 dy\int_{y}^{\sqrt y} f(x,y)\,dx$.

Пример 11 из примера ДЗ2 в РПД

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

$$\int_3^7 dx\int_{9/x}^{3} f(x,y)\,dy + \int_7^9 dx\int_{9/x}^{10-x} f(x,y)\,dy.$$
Решение

Шаг 1. Области. Первый интеграл: $3\le x\le 7$, $\frac9x\le y\le 3$. Второй: $7\le x\le 9$, $\frac9x\le y\le 10 - x$.

Шаг 2. Рисунок. Нижняя граница везде гипербола $y = \frac9x$. Верхняя: прямая $y = 3$ до $x = 7$, дальше прямая $y = 10 - x$ (при $x = 7$ она тоже даёт $y = 3$ — части стыкуются). Угловые точки: $(3,3)$ — гипербола и $y = 3$; $(7,3)$; $(9,1)$ — гипербола и $x+y = 10$ ($9/x = 10 - x \Rightarrow x^2 - 10x + 9 = 0 \Rightarrow x = 9$). Получился криволинейный треугольник. Верхняя граница из двух кусков — поэтому в исходном порядке и было два интеграла.

Шаг 3. По $y$ область занимает $[1, 3]$.

Шаг 4. Относительно $x$: гипербола $x = \frac9y$, прямая $x = 10 - y$. Горизонталь при любом $y\in[1,3]$ входит через гиперболу, выходит через наклонную прямую (прямая $y=3$ — это верхний край, а не правый). Один интеграл:

$$\int_1^3 dy\int_{9/y}^{10-y} f(x,y)\,dx.$$

Шаг 5. Проверка площадью. $\int_1^3\bigl(10 - y - \frac9y\bigr)dy = 20 - 4 - 9\ln 3 = 16 - 9\ln 3$. Исходные интегралы при $f=1$: $\int_3^7\bigl(3 - \frac9x\bigr)dx + \int_7^9\bigl(10 - x - \frac9x\bigr)dx = (12 - 9\ln\frac73) + (4 - 9\ln\frac97) = 16 - 9\ln 3$. Совпало.

Ответ: $\displaystyle\int_1^3 dy\int_{9/y}^{10-y} f(x,y)\,dx$.

Пример 12 (после смены порядка — два интеграла)

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

Решение

Область: $0\le y\le 1$, слева дуга $x = -\sqrt{1 - y^2}$ (левая половина окружности $x^2 + y^2 = 1$), справа прямая $x = 1 - y$. Это четверть круга во второй четверти плюс треугольник с вершинами $(0,0)$, $(1,0)$, $(0,1)$.

Новый внешний $x\in[-1,1]$. Вертикаль всегда входит через $y = 0$, а выходит при $x\in[-1,0]$ через окружность $y = \sqrt{1 - x^2}$, при $x\in[0,1]$ — через прямую $y = 1 - x$. Верхняя граница меняется при $x = 0$:

$$\int_{-1}^{0}dx\int_0^{\sqrt{1-x^2}}f\,dy + \int_0^1 dx\int_0^{1-x}f\,dy.$$

Проверка площадью: $\frac\pi4 + \frac12$ в обоих порядках.

Ответ: $\displaystyle\int_{-1}^{0}dx\int_0^{\sqrt{1-x^2}}f\,dy + \int_0^1 dx\int_0^{1-x}f\,dy$.

Пример 13 (угол между прямыми)

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

Решение

Область: треугольник с вершинами $(0,0)$, $(2,2)$, $(2,4)$. По $y$ от 0 до 4. Относительно $x$: $y = 2x\Rightarrow x = \frac y2$ (левая граница), $y = x \Rightarrow x = y$ и отрезок $x = 2$ (правые).

Горизонталь входит всегда через $x = \frac y2$. Выходит: при $y\le 2$ через прямую $x = y$, при $y\ge 2$ через $x = 2$. Разбиваем по $y = 2$:

$$\int_0^2 dy\int_{y/2}^{y} f\,dx + \int_2^4 dy\int_{y/2}^{2} f\,dx.$$

Проверка на $f = xy$: исходный $\int_0^2 x\cdot\frac{4x^2 - x^2}{2}dx = \frac32\cdot 4 = 6$; новый $\int_0^2 y\cdot\frac{y^2 - y^2/4}{2}dy + \int_2^4 y\cdot\frac{4 - y^2/4}{2}dy = \frac32 + \frac92 = 6$.

Ответ: $\displaystyle\int_0^2 dy\int_{y/2}^{y} f\,dx + \int_2^4 dy\int_{y/2}^{2} f\,dx$.

Пример 14 (обратные тригонометрические)

Измените порядок: $\displaystyle\int_0^\pi dx\int_0^{\sin x} f(x,y)\,dy$.

Решение

Область под аркой синусоиды, $y$ от 0 до 1. Уравнение $y = \sin x$ на $[0,\pi]$ даёт два решения: $x = \arcsin y$ (левая половина арки, $x\in[0,\frac\pi2]$) и $x = \pi - \arcsin y$ (правая). Ловушка: ответ $x = \arcsin y$ для правой половины неверен, $\arcsin$ возвращает значения только из $[-\frac\pi2,\frac\pi2]$.

$$\int_0^\pi dx\int_0^{\sin x} f\,dy = \int_0^1 dy\int_{\arcsin y}^{\pi - \arcsin y} f\,dx.$$

Проверка на $f = y$: слева $\int_0^\pi\frac{\sin^2 x}{2}dx = \frac\pi4$; справа $\int_0^1 y(\pi - 2\arcsin y)\,dy = \frac\pi2 - 2\cdot\frac\pi8 = \frac\pi4$ (здесь $\int_0^1 y\arcsin y\,dy = \frac\pi8$, по частям).

Ответ: $\displaystyle\int_0^1 dy\int_{\arcsin y}^{\pi - \arcsin y} f(x,y)\,dx$.

Пример 15 (без смены порядка не взять)

Вычислите $\displaystyle\int_0^1 dx\int_x^1 e^{y^2}\,dy$.

Решение

Внутренний интеграл $\int e^{y^2}dy$ не выражается через элементарные функции. Меняем порядок. Область: $0\le x\le 1$, $x\le y\le 1$ — треугольник с вершинами $(0,0)$, $(0,1)$, $(1,1)$. В другом порядке: $y\in[0,1]$, $0\le x\le y$:

$$\int_0^1 dy\int_0^y e^{y^2}dx = \int_0^1 y\,e^{y^2}dy = \frac12 e^{y^2}\Big|_0^1 = \frac{e - 1}{2}.$$

Внутренний интеграл по $x$ от функции, не зависящей от $x$, дал множитель $y$ — ровно то, что нужно для замены $t = y^2$.

Ответ: $\frac{e-1}{2}$.

Ловушка

При смене порядка нельзя просто переставить пределы местами: $\int_0^1 dx\int_{x^2}^{x}f\,dy \ne \int_{x^2}^{x}dy\int_0^1 f\,dx$. Новые пределы всегда получаются только из рисунка.

5. Криволинейные координаты и замена переменных

вопрос экзамена №17вопрос экзамена №18практика 1, ДЗ2 задача 3, КР1 задание 5в (примеры в РПД)

Криволинейные координаты и якобиан

В определённом интеграле замена $x = x(t)$ даёт $dx = x'(t)\,dt$: множитель $x'(t)$ показывает, во сколько раз растягиваются маленькие отрезки. В двойном интеграле то же самое, только растягиваются маленькие площадки, и коэффициент растяжения — определитель.

Определение

Пусть формулы $x = x(u,v)$, $y = y(u,v)$ взаимно однозначно отображают область $\Delta$ плоскости $(u,v)$ на область $D$ плоскости $(x,y)$, функции непрерывно дифференцируемы. Тогда пару $(u,v)$ называют криволинейными координатами точки $(x,y)$. Линии $u = \mathrm{const}$ и $v = \mathrm{const}$ — координатные линии (для полярных: окружности $r = \mathrm{const}$ и лучи $\varphi = \mathrm{const}$).

Якобиан (функциональный определитель) отображения:

$$J(u,v) = \frac{\partial(x,y)}{\partial(u,v)} = \begin{vmatrix} x'_u & x'_v \\ y'_u & y'_v \end{vmatrix} = x'_u y'_v - x'_v y'_u.$$
Площадь в криволинейных координатах

Маленький прямоугольник $[u, u+du]\times[v, v+dv]$ переходит в криволинейный четырёхугольник, который почти совпадает с параллелограммом, натянутым на векторы $(x'_u, y'_u)\,du$ и $(x'_v, y'_v)\,dv$. Площадь параллелограмма — модуль определителя из координат векторов, то есть $|J|\,du\,dv$. Отсюда

$$dS = |J(u,v)|\,du\,dv, \qquad S(D) = \iint_\Delta |J(u,v)|\,du\,dv.$$
Теорема (замена переменных в двойном интеграле)

Пусть отображение $x = x(u,v)$, $y = y(u,v)$ взаимно однозначно переводит замкнутую измеримую область $\Delta$ в $D$, непрерывно дифференцируемо и $J(u,v)\ne 0$ в $\Delta$ (допускается обращение в нуль или нарушение однозначности на множестве меры нуль, например в одной точке или на отрезке границы). Тогда для непрерывной на $D$ функции $f$

$$\iint_D f(x,y)\,dx\,dy = \iint_\Delta f\bigl(x(u,v),\,y(u,v)\bigr)\,|J(u,v)|\,du\,dv.$$
Полезный факт

Якобиан обратного отображения — обратное число: $\frac{\partial(x,y)}{\partial(u,v)} = 1\Big/\frac{\partial(u,v)}{\partial(x,y)}$. Если новые переменные заданы как $u = u(x,y)$, $v = v(x,y)$, проще посчитать $\frac{\partial(u,v)}{\partial(x,y)}$ и перевернуть, чем выражать $x$ и $y$ (примеры 24, 25).

Полярные координаты

$x = r\cos\varphi$, $y = r\sin\varphi$, $r\ge 0$, $\varphi\in[0, 2\pi)$ (или $(-\pi,\pi]$, как удобнее).

$$J = \begin{vmatrix} \cos\varphi & -r\sin\varphi \\ \sin\varphi & r\cos\varphi \end{vmatrix} = r\cos^2\varphi + r\sin^2\varphi = r, \qquad \boxed{dx\,dy = r\,dr\,d\varphi}.$$

Геометрически: «полярный прямоугольник» $[r, r+dr]\times[\varphi,\varphi+d\varphi]$ — кусок кольца со сторонами $dr$ (вдоль луча) и $r\,d\varphi$ (дуга). Его площадь $\approx r\,dr\,d\varphi$. Чем дальше от центра, тем длиннее дуга при том же $d\varphi$ — отсюда множитель $r$.

Когда и как переходить к полярным
  1. Когда выгодно: область ограничена окружностями (с центром в начале координат или проходящими через него) и лучами из начала; подынтегральная функция содержит $x^2 + y^2$.
  2. Подставь $x = r\cos\varphi$, $y = r\sin\varphi$ в уравнения границ и в $f$. Не забудь множитель $r$.
  3. Внешний интеграл по $\varphi$: найди лучи из начала координат, между которыми лежит область (крайние углы).
  4. Внутренний по $r$: проведи луч под углом $\varphi$; $r$ на входе в область — нижний предел, на выходе — верхний. Если начало координат внутри области или на её границе, вход при $r = 0$.
  5. Если начало координат внутри области, $\varphi$ пробегает полный оборот $[0,2\pi]$.
Линия в декартовыхВ полярныхКомментарий
$x^2 + y^2 = R^2$$r = R$окружность с центром в $O$
$x^2 + y^2 = 2ax$$r = 2a\cos\varphi$, $\varphi\in[-\frac\pi2,\frac\pi2]$окружность с центром $(a,0)$, радиус $a$, проходит через $O$
$x^2 + y^2 = 2by$$r = 2b\sin\varphi$, $\varphi\in[0,\pi]$окружность с центром $(0,b)$, радиус $b$
$y = kx$ (луч)$\varphi = \arctg k$ (или $+\pi$)при $k = 1$: $\varphi = \frac\pi4$
$x = a$$r = \frac{a}{\cos\varphi}$вертикальная прямая
$y = b$$r = \frac{b}{\sin\varphi}$горизонтальная прямая
$\frac{x^2}{a^2}+\frac{y^2}{b^2} = 1$$r = 1$ в обобщённых полярныхсм. ниже

Вывод второй строки: $x^2+y^2 = 2ax \Rightarrow r^2 = 2ar\cos\varphi \Rightarrow r = 2a\cos\varphi$ (сократили на $r$, точку $r = 0$ не теряем: она получается при $\varphi = \pm\frac\pi2$). Условие $r\ge 0$ даёт $\cos\varphi \ge 0$, отсюда промежуток для $\varphi$. Выделив полный квадрат, $(x - a)^2 + y^2 = a^2$ — центр и радиус.

Обобщённые полярные координаты

Для эллипсов: $x = a\,r\cos\varphi$, $y = b\,r\sin\varphi$. Тогда $\frac{x^2}{a^2} + \frac{y^2}{b^2} = r^2$, эллипс превращается в $r = 1$, а

$$J = \begin{vmatrix} a\cos\varphi & -ar\sin\varphi \\ b\sin\varphi & br\cos\varphi \end{vmatrix} = abr, \qquad dx\,dy = ab\,r\,dr\,d\varphi.$$

Важно: здесь $\varphi$ — не настоящий угол точки (луч $y = x$ при $a\ne b$ не соответствует $\varphi = \frac\pi4$), это просто параметр.

Пример 16 из примера КР1 в РПД, задание 5в

Плоская область $D$ ограничена линиями $y = \sqrt3\,x$, $y = -\frac{x}{\sqrt3}$, $y = 1$. Получите систему неравенств, описывающую её в полярных координатах. Изобразите область. (Дополнительно: найдите площадь двумя способами.)

Решение

Лучи. $y = \sqrt3\,x$: $\tg\varphi = \sqrt3$. Область лежит ниже $y = 1$ и выше обеих прямых, значит в верхней полуплоскости, поэтому $\varphi = \frac\pi3$. $y = -\frac{x}{\sqrt3}$: $\tg\varphi = -\frac1{\sqrt3}$, в верхней полуплоскости $\varphi = \frac{5\pi}{6}$. Прямые перпендикулярны: разность углов $\frac\pi2$.

Прямая $y = 1$: $r\sin\varphi = 1 \Rightarrow r = \frac{1}{\sin\varphi}$.

Луч из начала координат под углом $\varphi\in[\frac\pi3, \frac{5\pi}{6}]$ стартует в $O$ (вершина треугольника) и выходит на прямой $y = 1$:

$$D:\quad \frac\pi3 \le \varphi \le \frac{5\pi}{6}, \qquad 0 \le r \le \frac{1}{\sin\varphi}.$$

Проверка площадью. Вершины: $(0,0)$, $(\frac{1}{\sqrt3}, 1)$, $(-\sqrt3, 1)$, основание на $y = 1$ длины $\frac1{\sqrt3} + \sqrt3 = \frac{4}{\sqrt3}$, высота 1, $S = \frac{2}{\sqrt3}$. В полярных:

$$S = \int_{\pi/3}^{5\pi/6}d\varphi\int_0^{1/\sin\varphi} r\,dr = \int_{\pi/3}^{5\pi/6}\frac{d\varphi}{2\sin^2\varphi} = -\frac12\ctg\varphi\Big|_{\pi/3}^{5\pi/6} = \frac12\Bigl(\sqrt3 + \frac{1}{\sqrt3}\Bigr) = \frac{2}{\sqrt3}.$$

Ответ: $\frac\pi3\le\varphi\le\frac{5\pi}6$, $0\le r\le\frac1{\sin\varphi}$; $S = \frac{2}{\sqrt3}$.

Пример 17 (кольцо)

Вычислите $\iint_D \ln(x^2+y^2)\,dx\,dy$, $D$: $1\le x^2+y^2\le 4$.

Решение

Кольцо: $\varphi\in[0,2\pi]$, $1\le r\le 2$, $\ln(x^2+y^2) = \ln r^2 = 2\ln r$. Пределы постоянные, интеграл распадается в произведение:

$$\int_0^{2\pi}d\varphi\int_1^2 2\ln r\cdot r\,dr = 2\pi\Bigl(r^2\ln r - \frac{r^2}{2}\Bigr)\Big|_1^2 = 2\pi\Bigl(4\ln 2 - 2 + \frac12\Bigr) = \pi(8\ln 2 - 3).$$

$\int 2r\ln r\,dr$ — по частям: $u = \ln r$, $dv = 2r\,dr$, $v = r^2$, получаем $r^2\ln r - \int r\,dr$.

Ответ: $\pi(8\ln 2 - 3)\approx 8{,}00$.

Пример 18 (окружность через начало координат)

Вычислите $\iint_D y\,dx\,dy$, где $D$: $x^2+y^2\le 2y$.

Решение

$x^2 + (y-1)^2\le 1$ — круг радиуса 1 с центром $(0,1)$, касается оси $Ox$ в начале координат. В полярных: $r^2\le 2r\sin\varphi \Rightarrow r\le 2\sin\varphi$, нужно $\sin\varphi\ge 0$, то есть $\varphi\in[0,\pi]$. Луч выходит из $O$ (на границе) и выходит из круга при $r = 2\sin\varphi$:

$$\int_0^\pi d\varphi\int_0^{2\sin\varphi} r\sin\varphi\cdot r\,dr = \int_0^\pi \sin\varphi\cdot\frac{8\sin^3\varphi}{3}\,d\varphi = \frac83\int_0^\pi\sin^4\varphi\,d\varphi = \frac83\cdot\frac{3\pi}{8} = \pi.$$

Здесь $\sin^4\varphi = \frac{3}{8} - \frac12\cos2\varphi + \frac18\cos 4\varphi$, на $[0,\pi]$ косинусы дают ноль.

Проверка по смыслу: $\iint_D y\,dS = y_c\cdot S$, центр круга $y_c = 1$, площадь $\pi$ (раздел 7).

Ответ: $\pi$.

Пример 19 из примера практики 1 в РПД

Перейдя к полярным координатам, вычислите $\iint_G (x^2 + y^2)\,dx\,dy$, где $G$ ограничена кривыми $x^2 + y^2 = ax$, $x^2+y^2 = 2ax$, $y = 0$ ($y \gt 0$).

Решение

Обе окружности проходят через начало координат: $x^2+y^2 = ax$ — центр $(\frac a2, 0)$, радиус $\frac a2$; $x^2+y^2 = 2ax$ — центр $(a,0)$, радиус $a$. Малая лежит внутри большой, касаясь её в точке $O$. Область — верхняя половина «полумесяца» между ними.

В полярных: $r = a\cos\varphi$ и $r = 2a\cos\varphi$. Условие $y \gt 0$ и $\cos\varphi\ge0$ дают $\varphi\in[0,\frac\pi2]$. Луч входит в область на малой окружности, выходит на большой:

$$\iint_G(x^2+y^2)\,dx\,dy = \int_0^{\pi/2}d\varphi\int_{a\cos\varphi}^{2a\cos\varphi} r^2\cdot r\,dr = \int_0^{\pi/2}\frac{16 - 1}{4}a^4\cos^4\varphi\,d\varphi = \frac{15a^4}{4}\cdot\frac{3\pi}{16} = \frac{45\pi a^4}{64}.$$

$\int_0^{\pi/2}\cos^4\varphi\,d\varphi = \frac{3\pi}{16}$: $\cos^4\varphi = \frac38 + \frac12\cos2\varphi + \frac18\cos4\varphi$.

Ответ: $\frac{45\pi a^4}{64}$.

Пример 20 из примера практики 1 в РПД

Найдите площадь фигуры, ограниченной кривыми $x^2+y^2 = 2ax$, $x^2+y^2 = 2bx$, $y = x$, $y = 0$ ($0 \lt a \lt b$).

Решение

Окружности $r = 2a\cos\varphi$ и $r = 2b\cos\varphi$ касаются в начале координат, меньшая внутри. Прямые $y = 0$ и $y = x$ — лучи $\varphi = 0$ и $\varphi = \frac\pi4$. Площадь:

$$S = \int_0^{\pi/4}d\varphi\int_{2a\cos\varphi}^{2b\cos\varphi}r\,dr = \int_0^{\pi/4}2(b^2 - a^2)\cos^2\varphi\,d\varphi = 2(b^2-a^2)\Bigl(\frac\varphi2 + \frac{\sin2\varphi}{4}\Bigr)\Big|_0^{\pi/4} = 2(b^2-a^2)\Bigl(\frac\pi8 + \frac14\Bigr).$$

Ответ: $S = \frac{(\pi + 2)(b^2 - a^2)}{4}$.

Пример 21 из примера практики 1 в РПД

Найдите площадь фигуры, ограниченной кривыми $r = a(1 + \cos\varphi)$ и $r = a\cos\varphi$ ($a \gt 0$).

Решение

$r = a(1+\cos\varphi)$ — кардиоида, $\varphi\in[0,2\pi]$. $r = a\cos\varphi$ — окружность $x^2+y^2 = ax$ диаметра $a$, $\varphi\in[-\frac\pi2,\frac\pi2]$. При каждом $\varphi$ $a\cos\varphi \lt a(1+\cos\varphi)$, так что круг лежит внутри кардиоиды (касание в $O$). Фигура — кардиоида без круга.

Способ 1 (разность площадей). Площадь кардиоиды: $\int_0^{2\pi}d\varphi\int_0^{a(1+\cos\varphi)}r\,dr = \frac{a^2}{2}\int_0^{2\pi}(1 + 2\cos\varphi + \cos^2\varphi)\,d\varphi = \frac{a^2}{2}(2\pi + \pi) = \frac{3\pi a^2}{2}$. Площадь круга $\frac{\pi a^2}{4}$. Разность $\frac{5\pi a^2}{4}$.

Способ 2 (сразу двойным интегралом). Луч под углом $\varphi\in[-\frac\pi2,\frac\pi2]$ входит в фигуру на окружности, выходит на кардиоиде; при $\varphi\in[\frac\pi2,\frac{3\pi}{2}]$ окружности нет, вход при $r = 0$:

$$S = \int_{-\pi/2}^{\pi/2}d\varphi\int_{a\cos\varphi}^{a(1+\cos\varphi)}r\,dr + \int_{\pi/2}^{3\pi/2}d\varphi\int_0^{a(1+\cos\varphi)}r\,dr = \frac{a^2}{2}(\pi + 4) + \frac{a^2}{2}\Bigl(\frac{3\pi}{2} - 4\Bigr) = \frac{5\pi a^2}{4}.$$

В первом слагаемом $(1+\cos\varphi)^2 - \cos^2\varphi = 1 + 2\cos\varphi$; во втором $\int_{\pi/2}^{3\pi/2}\cos\varphi\,d\varphi = -2$, $\int_{\pi/2}^{3\pi/2}\cos^2\varphi\,d\varphi = \frac\pi2$.

На рисунке к задаче в РПД заштрихована только верхняя половина фигуры ($y\ge 0$). Если спрашивают её, ответ вдвое меньше: $\frac{5\pi a^2}{8}$ — фигура симметрична относительно $Ox$.

Ответ: $S = \frac{5\pi a^2}{4}$ (верхняя половина $\frac{5\pi a^2}{8}$).

Пример 22 из примера ДЗ2 в РПД

Вычислите двойной интеграл, перейдя к обобщённым полярным координатам $x = ar\cos\varphi$, $y = br\sin\varphi$:

$$\iint_G\frac{dx\,dy}{\sqrt{c^2 - \frac{x^2}{a^2} - \frac{y^2}{b^2}}}\quad (c \gt 1),$$

где область $G$ ограничена эллипсом $\frac{x^2}{a^2}+\frac{y^2}{b^2} = 1$.

Решение

После замены $\frac{x^2}{a^2} + \frac{y^2}{b^2} = r^2$, эллипс: $0\le r\le 1$, $0\le\varphi\le 2\pi$, $dx\,dy = ab\,r\,dr\,d\varphi$. Подкоренное выражение $c^2 - r^2 \gt 0$ при $c \gt 1$, интеграл собственный.

$$\iint_G = \int_0^{2\pi}d\varphi\int_0^1\frac{ab\,r\,dr}{\sqrt{c^2 - r^2}} = 2\pi ab\Bigl(-\sqrt{c^2 - r^2}\Bigr)\Big|_0^1 = 2\pi ab\bigl(c - \sqrt{c^2 - 1}\bigr).$$

Внутренний интеграл: $\int\frac{r\,dr}{\sqrt{c^2 - r^2}} = -\sqrt{c^2 - r^2}$ (замена $t = c^2 - r^2$).

Ответ: $2\pi ab\bigl(c - \sqrt{c^2-1}\bigr)$.

Пример 23 (площадь эллипса)

Найдите площадь эллипса $\frac{x^2}{a^2}+\frac{y^2}{b^2}\le 1$.

Решение

Обобщённые полярные: $S = \int_0^{2\pi}d\varphi\int_0^1 ab\,r\,dr = 2\pi\cdot\frac{ab}{2} = \pi ab$. При $a = b = R$ получаем площадь круга $\pi R^2$.

Ответ: $\pi ab$.

Пример 24 (линейная замена)

Вычислите $\iint_D (x+y)^2\,dx\,dy$, где $D$ — параллелограмм, ограниченный прямыми $x + y = 1$, $x + y = 3$, $x - y = 0$, $x - y = 2$.

Решение

В декартовых координатах пришлось бы резать на три куска. Вводим $u = x + y$, $v = x - y$: область становится прямоугольником $1\le u\le 3$, $0\le v\le 2$.

Якобиан проще считать обратный: $\frac{\partial(u,v)}{\partial(x,y)} = \begin{vmatrix}1 & 1\\ 1 & -1\end{vmatrix} = -2$, значит $\left|\frac{\partial(x,y)}{\partial(u,v)}\right| = \frac12$. (Можно и прямо: $x = \frac{u+v}{2}$, $y = \frac{u - v}{2}$, определитель $-\frac12$.)

$$\iint_D (x+y)^2\,dx\,dy = \int_1^3 du\int_0^2 u^2\cdot\frac12\,dv = \int_1^3 u^2\,du = \frac{27 - 1}{3} = \frac{26}{3}.$$

Ответ: $\frac{26}{3}$.

Пример 25 (площадь в криволинейных координатах)

Найдите площадь области в первой четверти, ограниченной гиперболами $xy = 1$, $xy = 2$ и прямыми $y = x$, $y = 4x$.

Решение

Границы сами подсказывают координаты: $u = xy\in[1,2]$, $v = \frac yx\in[1,4]$. В плоскости $(u,v)$ область — прямоугольник.

$$\frac{\partial(u,v)}{\partial(x,y)} = \begin{vmatrix} y & x\\ -\frac{y}{x^2} & \frac1x\end{vmatrix} = \frac yx + \frac yx = \frac{2y}{x} = 2v, \qquad |J| = \frac{1}{2v}.$$ $$S = \int_1^2 du\int_1^4\frac{dv}{2v} = 1\cdot\frac{\ln 4}{2} = \ln 2.$$

Проверено численно в декартовых координатах (там три куска по $x$).

Ответ: $S = \ln 2$.

Ловушки замены

6. Несобственные двойные интегралы

вопрос экзамена №19вопрос экзамена №20

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

Определение (неограниченная область)

Пусть $D$ неограниченна, $f$ интегрируема на каждой ограниченной измеримой части $D$. Возьмём последовательность ограниченных измеримых областей $D_1\subset D_2\subset\ldots\subset D$, «исчерпывающую» $D$: любая ограниченная часть $D$ рано или поздно попадает в $D_n$. Если предел

$$\iint_D f\,dx\,dy = \lim_{n\to\infty}\iint_{D_n} f\,dx\,dy$$

существует, конечен и не зависит от выбора последовательности $D_n$, интеграл называется сходящимся.

Определение (неограниченная функция)

Пусть $f$ неограниченна в окрестности точки $P_0\in D$ (или линии) и интегрируема вне любой её окрестности. Вырезаем из $D$ малую окрестность $U_\varepsilon$ особой точки и полагаем

$$\iint_D f\,dx\,dy = \lim_{\varepsilon\to0}\iint_{D\setminus U_\varepsilon} f\,dx\,dy,$$

если предел конечен и не зависит от формы вырезаемых окрестностей.

Что полезно знать
Пример 26 (интеграл Пуассона)

Вычислите $\iint_{\R^2} e^{-x^2-y^2}\,dx\,dy$ и выведите отсюда $\int_{-\infty}^{+\infty}e^{-x^2}dx = \sqrt\pi$.

Решение

Шаг 1. $f \gt 0$, исчерпываем плоскость кругами $D_R$: $r\le R$:

$$\iint_{D_R}e^{-x^2-y^2}dx\,dy = \int_0^{2\pi}d\varphi\int_0^R e^{-r^2}r\,dr = 2\pi\cdot\Bigl(-\frac12 e^{-r^2}\Bigr)\Big|_0^R = \pi\bigl(1 - e^{-R^2}\bigr)\xrightarrow[R\to\infty]{}\pi.$$

Шаг 2. Тот же интеграл по квадратам $K_a = [-a,a]^2$ (для $f\ge0$ можно, ответ тот же). На квадрате переменные разделяются:

$$\iint_{K_a}e^{-x^2}e^{-y^2}dx\,dy = \int_{-a}^{a}e^{-x^2}dx\cdot\int_{-a}^{a}e^{-y^2}dy = \Bigl(\int_{-a}^{a}e^{-x^2}dx\Bigr)^2.$$

Если не хочется ссылаться на независимость от последовательности, есть прямое сжатие: $D_a\subset K_a\subset D_{a\sqrt2}$, поэтому $\pi(1 - e^{-a^2})\le\bigl(\int_{-a}^a e^{-x^2}dx\bigr)^2\le\pi(1 - e^{-2a^2})$, и обе границы стремятся к $\pi$.

Шаг 3. При $a\to\infty$: $\bigl(\int_{-\infty}^{+\infty}e^{-x^2}dx\bigr)^2 = \pi$, интеграл положителен, значит $\int_{-\infty}^{+\infty}e^{-x^2}dx = \sqrt\pi$.

Через $e^{-x^2}$ первообразную не выразить, а двойной интеграл с множителем $r$ из якобиана берётся сразу. Этот интеграл нужен в теории вероятностей (нормальное распределение).

Ответ: $\pi$; $\int_{-\infty}^{+\infty}e^{-x^2}dx = \sqrt\pi$.

Пример 27 (неограниченная функция)

Исследуйте на сходимость и вычислите $\iint_D\frac{dx\,dy}{\sqrt{x^2+y^2}}$, $D$: $x^2+y^2\le 1$.

Решение

В точке $O$ функция неограниченна. Вырезаем круг $r\le\varepsilon$:

$$\iint_{\varepsilon\le r\le 1}\frac{r\,dr\,d\varphi}{r} = 2\pi(1 - \varepsilon)\xrightarrow[\varepsilon\to0]{}2\pi.$$

Интеграл сходится и равен $2\pi$: множитель $r$ из якобиана «гасит» особенность. Для $\frac{1}{x^2+y^2}$ уже было бы $2\pi\int_\varepsilon^1\frac{dr}{r} = -2\pi\ln\varepsilon\to\infty$ — расходится ($\alpha = 2$ — граница).

Ответ: сходится, $2\pi$.

Пример 28 (внешность круга)

Вычислите $\iint_D\frac{dx\,dy}{(x^2+y^2)^2}$, $D$: $x^2+y^2\ge 1$.

Решение

Исчерпываем кольцами $1\le r\le R$:

$$\int_0^{2\pi}d\varphi\int_1^R\frac{r\,dr}{r^4} = 2\pi\Bigl(-\frac{1}{2r^2}\Bigr)\Big|_1^R = \pi\Bigl(1 - \frac{1}{R^2}\Bigr)\to\pi.$$

Здесь $\alpha = 4 \gt 2$ — сходится, как и предсказывает эталон.

Ответ: $\pi$.

7. Приложения: площадь, объём, масса, центр масс, моменты

вопрос экзамена №23вопрос экзамена №24КР2, задание 1РГР1, задание 1

Формулы
ВеличинаФормула
Площадь плоской фигуры$S = \iint_D dx\,dy$
Объём тела между $z = z_1(x,y)$ (снизу) и $z = z_2(x,y)$ (сверху) над $D$$V = \iint_D\bigl(z_2(x,y) - z_1(x,y)\bigr)\,dx\,dy$
Масса пластинки с плотностью $\rho(x,y)$$m = \iint_D\rho\,dx\,dy$
Статические моменты относительно осей$M_x = \iint_D y\,\rho\,dx\,dy,\quad M_y = \iint_D x\,\rho\,dx\,dy$
Центр масс$x_c = \frac{M_y}{m},\quad y_c = \frac{M_x}{m}$
Моменты инерции относительно осей $Ox$, $Oy$ и начала координат$I_x = \iint_D y^2\rho\,dS,\quad I_y = \iint_D x^2\rho\,dS,\quad I_O = \iint_D(x^2+y^2)\rho\,dS = I_x + I_y$
Как запомнить
Теорема Штейнера

Пусть $l_c$ — ось через центр масс пластинки, $l$ — параллельная ей ось на расстоянии $d$. Тогда

$$I_l = I_{l_c} + m\,d^2.$$

Доказательство для осей, параллельных $Ox$: пусть центр масс на оси $l_c$: $y = y_c$, ось $l$: $y = y_c + d$. Тогда $I_l = \iint (y - y_c - d)^2\rho\,dS = \iint (y-y_c)^2\rho\,dS - 2d\iint (y - y_c)\rho\,dS + d^2\iint\rho\,dS$. Средний интеграл равен $M_x - y_c m = 0$ по определению центра масс. Остаётся $I_{l_c} + md^2$.

Следствие: из всех параллельных осей момент инерции минимален относительно оси через центр масс.

Про потенциал (вопрос №23)

В вопросе 23 упомянут потенциал гравитационного поля. Для массы, распределённой по пластинке $D$ с плотностью $\rho$, потенциал в точке $M_0(x_0,y_0,z_0)$ вне пластинки (в нормировке, где потенциал точечной массы $m$ равен $\frac{\gamma m}{r}$, $\gamma$ — гравитационная постоянная) — это сумма вкладов маленьких кусочков:

$$u(M_0) = \gamma\iint_D\frac{\rho(x,y)\,dx\,dy}{\sqrt{(x - x_0)^2 + (y - y_0)^2 + z_0^2}}.$$

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

Пример 29 из примера КР2 в РПД, задание 1

Фигура ограничена линиями $x + 2y - 2 = 0$, $x = (y-1)^2$. При помощи двойного интеграла найдите её площадь. Проверьте полученный результат, изменив порядок интегрирования в повторном интеграле.

Решение

Рисунок и точки пересечения. $x = (y-1)^2$ — парабола с вершиной $(0,1)$, ветви вправо. Прямая: $x = 2 - 2y$, или $y = 1 - \frac x2$. Пересечение: $(y-1)^2 = 2 - 2y \Rightarrow y^2 - 2y + 1 = 2 - 2y \Rightarrow y^2 = 1$, $y = \pm1$. Точки $(0,1)$ и $(4,-1)$.

Способ 1 (внешний по $y$). При $y\in[-1,1]$ горизонталь входит через параболу $x = (y-1)^2$, выходит через прямую $x = 2 - 2y$ (при $y = 0$: $1\le x\le 2$, верно):

$$S = \int_{-1}^{1}dy\int_{(y-1)^2}^{2-2y}dx = \int_{-1}^{1}\bigl(2 - 2y - y^2 + 2y - 1\bigr)dy = \int_{-1}^{1}(1 - y^2)\,dy = 2 - \frac23 = \frac43.$$

Способ 2 (внешний по $x$, проверка). Решаем относительно $y$: парабола $y = 1\pm\sqrt x$, прямая $y = 1 - \frac x2$. При $x\in[0,4]$ вертикаль входит через нижнюю ветвь $y = 1 - \sqrt x$ и выходит через прямую (при $x = 1$: от $y = 0$ до $y = \frac12$; верхняя ветвь $y = 1+\sqrt x$ лежит выше прямой и область не ограничивает):

$$S = \int_0^4 dx\int_{1-\sqrt x}^{1 - x/2}dy = \int_0^4\Bigl(\sqrt x - \frac x2\Bigr)dx = \frac23\cdot 8 - 4 = \frac43.$$

Результаты совпали. Если бы фигура лежала по разные стороны от вершины параболы по $x$, во втором способе пришлось бы разбивать на части — здесь повезло.

Ответ: $S = \frac43$.

Пример 30 (объём под параболоидом)

Найдите объём тела, ограниченного параболоидом $z = 4 - x^2 - y^2$ и плоскостью $z = 0$.

Решение

Параболоид пересекает плоскость $z = 0$ по окружности $x^2+y^2 = 4$. Проекция тела на $xOy$ — круг $r\le 2$, верх $z_2 = 4 - r^2$, низ $z_1 = 0$:

$$V = \iint_{r\le2}(4 - x^2 - y^2)\,dx\,dy = \int_0^{2\pi}d\varphi\int_0^2(4 - r^2)\,r\,dr = 2\pi\Bigl(2r^2 - \frac{r^4}{4}\Bigr)\Big|_0^2 = 2\pi\cdot 4 = 8\pi.$$

Ответ: $V = 8\pi$.

Пример 31 (объём между двумя поверхностями)

Найдите объём тела, ограниченного параболоидами $z = x^2 + y^2$ и $z = 2 - x^2 - y^2$.

Решение

Проекция. Поверхности пересекаются, где $x^2+y^2 = 2 - x^2 - y^2$, то есть $x^2+y^2 = 1$ (на высоте $z = 1$). Проекция тела — круг $r\le 1$. Внутри круга сверху $2 - r^2$, снизу $r^2$ (проверка в центре: $2 \gt 0$).

$$V = \iint_{r\le 1}\bigl((2 - r^2) - r^2\bigr)\,dx\,dy = \int_0^{2\pi}d\varphi\int_0^1(2 - 2r^2)\,r\,dr = 2\pi\Bigl(r^2 - \frac{r^4}{2}\Bigr)\Big|_0^1 = \pi.$$

Ответ: $V = \pi$.

Пример 32 (объём внутри цилиндра)

Найдите объём тела, ограниченного поверхностями $z = x^2+y^2$, $z = 0$ и цилиндром $x^2 + y^2 = 2x$.

Решение

Проекция — круг $(x-1)^2 + y^2\le 1$, в полярных $r\le 2\cos\varphi$, $\varphi\in[-\frac\pi2,\frac\pi2]$. Сверху $z = r^2$, снизу $z = 0$:

$$V = \int_{-\pi/2}^{\pi/2}d\varphi\int_0^{2\cos\varphi}r^2\cdot r\,dr = \int_{-\pi/2}^{\pi/2}4\cos^4\varphi\,d\varphi = 4\cdot\frac{3\pi}{8} = \frac{3\pi}{2}.$$

Ответ: $V = \frac{3\pi}{2}$.

Пример 33 (центр тяжести полукруга)

Найдите центр тяжести однородного полукруга $x^2+y^2\le R^2$, $y\ge 0$.

Решение

Плотность $\rho = 1$ (для однородной пластинки константа сократится). По симметрии относительно $Oy$: $x_c = 0$. Масса $m = \frac{\pi R^2}{2}$.

$$M_x = \iint_D y\,dS = \int_0^\pi d\varphi\int_0^R r\sin\varphi\cdot r\,dr = \frac{R^3}{3}\int_0^\pi\sin\varphi\,d\varphi = \frac{2R^3}{3}, \qquad y_c = \frac{M_x}{m} = \frac{2R^3/3}{\pi R^2/2} = \frac{4R}{3\pi}.$$

Ответ: $\left(0, \frac{4R}{3\pi}\right)$, $y_c\approx 0{,}42R$.

Пример 34 (масса и центр масс неоднородной пластинки)

Пластинка — треугольник $x\ge0$, $y\ge0$, $x + y\le 1$ с плотностью $\rho = x$. Найдите массу и центр масс.

Решение

Пределы как в примере 4: $\int_0^1 dx\int_0^{1-x}\ldots\,dy$.

$$m = \int_0^1 dx\int_0^{1-x}x\,dy = \int_0^1 x(1-x)\,dx = \frac16,$$ $$M_y = \int_0^1 dx\int_0^{1-x}x\cdot x\,dy = \int_0^1 x^2(1-x)\,dx = \frac1{12}, \qquad M_x = \int_0^1 dx\int_0^{1-x}y\cdot x\,dy = \int_0^1\frac{x(1-x)^2}{2}dx = \frac{1}{24}.$$

$x_c = \frac{1/12}{1/6} = \frac12$, $y_c = \frac{1/24}{1/6} = \frac14$. Для однородного треугольника центр был бы в $\left(\frac13,\frac13\right)$; плотность растёт с $x$, и центр масс сдвинулся вправо.

Ответ: $m = \frac16$, центр масс $\left(\frac12, \frac14\right)$.

Пример 35 (моменты инерции, Штейнер)

Однородный диск радиуса $R$, плотность $\rho = 1$, центр в начале координат. Найдите $I_O$, $I_x$ и момент инерции относительно касательной $y = -R$.

Решение

$m = \pi R^2$.

$$I_O = \int_0^{2\pi}d\varphi\int_0^R r^2\cdot r\,dr = \frac{\pi R^4}{2} = \frac{mR^2}{2}.$$

По симметрии $I_x = I_y$, а $I_x + I_y = I_O$, значит $I_x = \frac{\pi R^4}{4} = \frac{mR^2}{4}$ (можно и напрямую: $\int_0^{2\pi}\sin^2\varphi\,d\varphi\int_0^R r^3dr = \pi\cdot\frac{R^4}{4}$).

Касательная $y = -R$ параллельна $Ox$ (ось через центр масс) на расстоянии $d = R$. По Штейнеру:

$$I_l = \frac{mR^2}{4} + mR^2 = \frac{5mR^2}{4} = \frac{5\pi R^4}{4}.$$

Проверка напрямую: $\iint_D (y + R)^2\,dS = I_x + 2R\iint_D y\,dS + R^2 m = \frac{\pi R^4}{4} + 0 + \pi R^4 = \frac{5\pi R^4}{4}$.

Ответ: $I_O = \frac{\pi R^4}{2}$, $I_x = \frac{\pi R^4}{4}$, $I_l = \frac{5\pi R^4}{4}$.

8. Как это выглядит в РГР1

РГР1, задание 1

В типовом расчёте (Попов, «Типовой расчёт», 3 семестр) двойным интегралам посвящено задание 1: «Плоская область $D$ ограничена заданными линиями. 1) Сделайте схематический рисунок области $D$. 2) С помощью двойного интеграла найдите площадь области $D$». В вариантах встречаются параболы и полупараболы ($x = \sqrt y$, $x = -\sqrt{2y}$), кубические корни, показательные функции $2^x$, $\left(\frac12\right)^x$, логарифмы, «локон Аньези» $y = \frac{2}{x^2+1}$, прямые. Объём в РГР считается тройным интегралом (задание 2, глава 7).

Оформление (по РПД)

Ниже полностью разобран вариант 1.

Пример 36 РГР1, задание 1, вариант 1

Плоская область $D$ ограничена линиями $y = \frac{2}{x^2+1}$, $y = \sqrt x$, $x = -\sqrt y$.

1) Сделайте схематический рисунок области $D$. 2) С помощью двойного интеграла найдите площадь области $D$.

Решение

1) Линии.

  • $y = \frac{2}{x^2+1}$ — «колокол» (локон Аньези): чётная функция, максимум $y(0) = 2$, $y\to 0$ при $x\to\pm\infty$, $y(\pm1) = 1$.
  • $y = \sqrt x$ — правая верхняя ветвь параболы $x = y^2$; определена при $x\ge 0$, $y\ge 0$.
  • $x = -\sqrt y$ — это $y = x^2$ при $x\le 0$: левая половина параболы $y = x^2$.

Точки пересечения.

  • Колокол и $y = \sqrt x$: $\frac{2}{x^2+1} = \sqrt x$. Подбором $x = 1$: $1 = 1$. Других корней нет: при $x\ge0$ левая часть убывает, правая возрастает. Точка $B(1,1)$.
  • Колокол и $y = x^2$ ($x\le 0$): $\frac{2}{x^2+1} = x^2 \Rightarrow x^4 + x^2 - 2 = 0 \Rightarrow (x^2 + 2)(x^2 - 1) = 0 \Rightarrow x^2 = 1$, $x = -1$. Точка $A(-1,1)$.
  • $y = \sqrt x$ и $x = -\sqrt y$: обе проходят через $O(0,0)$ и больше не пересекаются (одна при $x\ge0$, другая при $x\le 0$).

Область — криволинейный треугольник $OAB$ с вершинами $O(0,0)$, $A(-1,1)$, $B(1,1)$: сверху колокол, снизу слева $y = x^2$, снизу справа $y = \sqrt x$.

2) Площадь. $S = \iint_D dx\,dy$.

I способ (внешний по $x$). $x\in[-1,1]$. Вертикаль всегда выходит через колокол $y = \frac{2}{x^2+1}$, а входит при $x\in[-1,0]$ через $y = x^2$, при $x\in[0,1]$ через $y = \sqrt x$. Нижняя граница меняется при $x = 0$, делим $D$ на $D_1$ (слева) и $D_2$ (справа):

$$S = \int_{-1}^{0}dx\int_{x^2}^{\frac{2}{x^2+1}}dy + \int_0^1 dx\int_{\sqrt x}^{\frac{2}{x^2+1}}dy = \int_{-1}^{0}\Bigl(\frac{2}{x^2+1} - x^2\Bigr)dx + \int_0^1\Bigl(\frac{2}{x^2+1} - \sqrt x\Bigr)dx.$$ $$\int_{-1}^{0}\Bigl(\frac{2}{x^2+1} - x^2\Bigr)dx = \Bigl(2\arctg x - \frac{x^3}{3}\Bigr)\Big|_{-1}^{0} = \frac\pi2 - \frac13,\qquad \int_0^1\Bigl(\frac{2}{x^2+1} - \sqrt x\Bigr)dx = \Bigl(2\arctg x - \frac23x^{3/2}\Bigr)\Big|_0^1 = \frac\pi2 - \frac23.$$ $$S = \frac\pi2 - \frac13 + \frac\pi2 - \frac23 = \pi - 1.$$

II способ (внешний по $y$, проверка). Выражаем $x$: $y = x^2\ (x\le0) \Rightarrow x = -\sqrt y$; $y = \sqrt x \Rightarrow x = y^2$; колокол $\Rightarrow x^2 = \frac2y - 1$, $x = \pm\sqrt{\frac2y - 1}$ (определено при $0 \lt y\le 2$). По $y$ область занимает $[0,2]$, вершина колокола $(0,2)$.

  • $y\in[0,1]$ ($D_3$): горизонталь входит через $x = -\sqrt y$, выходит через $x = y^2$;
  • $y\in[1,2]$ ($D_4$): входит через левую половину колокола $x = -\sqrt{\frac2y - 1}$, выходит через правую $x = \sqrt{\frac2y-1}$.
$$S = \int_0^1 dy\int_{-\sqrt y}^{y^2}dx + \int_1^2 dy\int_{-\sqrt{2/y - 1}}^{\sqrt{2/y - 1}}dx = \int_0^1\bigl(y^2 + \sqrt y\bigr)dy + 2\int_1^2\sqrt{\frac{2 - y}{y}}\,dy.$$

Первый интеграл: $\frac13 + \frac23 = 1$. Во втором замена $y = 2\sin^2 t$, $dy = 4\sin t\cos t\,dt$, $\sqrt{\frac{2 - 2\sin^2 t}{2\sin^2 t}} = \frac{\cos t}{\sin t}$; при $y = 1$: $t = \frac\pi4$, при $y = 2$: $t = \frac\pi2$:

$$2\int_{\pi/4}^{\pi/2}\frac{\cos t}{\sin t}\cdot 4\sin t\cos t\,dt = 8\int_{\pi/4}^{\pi/2}\cos^2 t\,dt = 4\Bigl(t + \frac{\sin 2t}{2}\Bigr)\Big|_{\pi/4}^{\pi/2} = 4\Bigl(\frac\pi4 - \frac12\Bigr) = \pi - 2.$$

$S = 1 + \pi - 2 = \pi - 1$. Совпало с I способом. Видно, что II способ здесь заметно тяжелее: пределы по $y$ для колокола неудобные. В РГР основным берём I способ, II — как проверку.

Ответ: $S = \pi - 1\approx 2{,}14$ (кв. ед.).

Чек-лист задания 1 РГР
  1. Полное условие в тетради.
  2. Каждую линию разобрать отдельно: что это за кривая, какая её часть (например, $x = -\sqrt y$ — только левая половина параболы).
  3. Найти все точки пересечения аналитически; если корень угадан подбором, обосновать единственность (монотонность).
  4. Рисунок на миллиметровке: подписанные кривые, точки, штриховка, линии разбиения.
  5. Повторный интеграл с пояснением, через какую кривую входим и выходим; при смене границы — разбиение.
  6. Вычисление, проверка вторым порядком (или хотя бы грубой оценкой по рисунку: здесь треугольник $OAB$ имеет площадь 1, «шапка» колокола над $y = 1$ — чуть больше 1, так что $\pi - 1\approx 2{,}14$ правдоподобно).
  7. Ответ с единицами.

9. Тренировка

Сначала реши сам, потом открывай решение. Задачи идут от простых к уровню КР2 и ДЗ2.

Задача 1

Вычислите $\iint_D (x + y)\,dx\,dy$, $D$ ограничена линиями $y = x^2$, $y = 1$.

Решение

Пересечение $x^2 = 1$: $x = \pm1$. I тип: $x\in[-1,1]$, $x^2\le y\le 1$.

$$\int_{-1}^{1}dx\int_{x^2}^{1}(x+y)\,dy = \int_{-1}^{1}\Bigl(x(1 - x^2) + \frac{1 - x^4}{2}\Bigr)dx = 0 + \frac12\Bigl(2 - \frac25\Bigr) = \frac45.$$

Первое слагаемое — нечётная функция на симметричном отрезке, даёт 0.

Ответ: $\frac45$.

Задача 2

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

Решение

Область: $0\le y\le 1$, $y\le x\le 2-y$ — треугольник с вершинами $(0,0)$, $(2,0)$, $(1,1)$ (тот же, что в примере 7). По $x$ от 0 до 2, снизу $y = 0$, сверху при $x\le1$ прямая $y = x$, при $x\ge 1$ прямая $y = 2 - x$:

$$\int_0^1 dx\int_0^x f\,dy + \int_1^2 dx\int_0^{2-x}f\,dy.$$

Проверка площадью: $\frac12 + \frac12 = 1 = \int_0^1(2 - 2y)\,dy$.

Задача 3

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

Решение

Верхний полукруг радиуса 1. По $y$ от 0 до 1, горизонталь входит через левую дугу $x = -\sqrt{1 - y^2}$, выходит через правую:

$$\int_0^1 dy\int_{-\sqrt{1-y^2}}^{\sqrt{1-y^2}}f(x,y)\,dx.$$

Проверка площадью: $\frac\pi2$ в обоих порядках.

Задача 4

Вычислите $\iint_D\sqrt{x^2+y^2}\,dx\,dy$, $D$: $x^2+y^2\le 2x$.

Решение

Круг $(x-1)^2+y^2\le 1$, в полярных $0\le r\le 2\cos\varphi$, $\varphi\in[-\frac\pi2,\frac\pi2]$. Подынтегральная функция $r$, плюс якобиан $r$:

$$\int_{-\pi/2}^{\pi/2}d\varphi\int_0^{2\cos\varphi}r^2\,dr = \frac83\int_{-\pi/2}^{\pi/2}\cos^3\varphi\,d\varphi = \frac83\cdot\frac43 = \frac{32}{9}.$$

$\int_{-\pi/2}^{\pi/2}\cos^3\varphi\,d\varphi = \int(1 - \sin^2\varphi)\,d(\sin\varphi) = \Bigl(s - \frac{s^3}{3}\Bigr)\Big|_{-1}^{1} = \frac43$.

Ответ: $\frac{32}{9}$.

Задача 5

Найдите площадь фигуры, ограниченной параболами $y = x^2$ и $y = 2 - x^2$.

Решение

Пересечение: $x^2 = 2 - x^2 \Rightarrow x = \pm1$.

$$S = \int_{-1}^{1}dx\int_{x^2}^{2 - x^2}dy = \int_{-1}^{1}(2 - 2x^2)\,dx = 4 - \frac43 = \frac83.$$

Ответ: $\frac83$.

Задача 6

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

Решение

Проекция — треугольник $x, y\ge 0$, $x+y\le1$; сверху $z = x^2+y^2$, снизу $z = 0$.

$$V = \int_0^1 dx\int_0^{1-x}(x^2 + y^2)\,dy = \int_0^1\Bigl(x^2(1-x) + \frac{(1-x)^3}{3}\Bigr)dx = \frac1{12} + \frac1{12} = \frac16.$$

Ответ: $V = \frac16$.

Задача 7

Найдите массу и момент инерции $I_x$ круглой пластинки $x^2+y^2\le 1$ с плотностью $\rho = \sqrt{x^2+y^2}$ (растёт от центра к краю).

Решение
$$m = \int_0^{2\pi}d\varphi\int_0^1 r\cdot r\,dr = \frac{2\pi}{3},\qquad I_x = \int_0^{2\pi}\sin^2\varphi\,d\varphi\int_0^1 r^2\cdot r\cdot r\,dr = \pi\cdot\frac15 = \frac\pi5.$$

В $I_x$ три множителя $r$ после $r^2\sin^2\varphi$: $y^2 = r^2\sin^2\varphi$, плотность $r$, якобиан $r$ — итого $r^4\sin^2\varphi$.

Ответ: $m = \frac{2\pi}{3}$, $I_x = \frac\pi5$.

Задача 8

Вычислите $\iint_D e^{-(x^2+y^2)}\,dx\,dy$, $D$: $x^2+y^2\le R^2$. Что будет при $R\to\infty$?

Решение
$$\int_0^{2\pi}d\varphi\int_0^R e^{-r^2}r\,dr = \pi\bigl(1 - e^{-R^2}\bigr).$$

При $R\to\infty$ стремится к $\pi$ — это интеграл Пуассона из примера 26.

Ответ: $\pi\bigl(1 - e^{-R^2}\bigr)$.

Задача 9

Вычислите $\iint_D\Bigl(\frac{x^2}{4} + y^2\Bigr)dx\,dy$, $D$: $\frac{x^2}{4}+y^2\le 1$.

Решение

Обобщённые полярные с $a = 2$, $b = 1$: $x = 2r\cos\varphi$, $y = r\sin\varphi$, $|J| = 2r$, функция равна $r^2$:

$$\int_0^{2\pi}d\varphi\int_0^1 r^2\cdot 2r\,dr = 2\pi\cdot\frac24 = \pi.$$

Ответ: $\pi$.

Задача 10

Вычислите $\displaystyle\int_0^{\sqrt\pi}dx\int_x^{\sqrt\pi}\sin(y^2)\,dy$.

Решение

$\int\sin(y^2)\,dy$ не берётся. Область: треугольник $0\le x\le y\le\sqrt\pi$. Меняем порядок: $y\in[0,\sqrt\pi]$, $0\le x\le y$:

$$\int_0^{\sqrt\pi}dy\int_0^y\sin(y^2)\,dx = \int_0^{\sqrt\pi}y\sin(y^2)\,dy = -\frac12\cos(y^2)\Big|_0^{\sqrt\pi} = -\frac12(\cos\pi - 1) = 1.$$

Ответ: 1.

Область $1\le x^2+y^2\le 9$, $y\ge 0$. Какие пределы в полярных координатах верны?
Верхняя половина кольца: $y\ge0 \Leftrightarrow \varphi\in[0,\pi]$. Радиусы — корни из 1 и 9, то есть 1 и 3, а не сами числа.

10. Типичные ошибки

  1. Забыли якобиан $r$ в полярных координатах (или $abr$ в обобщённых). Проверка: площадь круга должна выйти $\pi R^2$.
  2. Неверные пределы по $\varphi$. Для $x^2+y^2 = 2ax$ угол от $-\frac\pi2$ до $\frac\pi2$, а не от 0 до $2\pi$. Для области в одной четверти — проверь знаки $x$ и $y$. Луч $y = -\frac{x}{\sqrt3}$ в верхней полуплоскости — это $\varphi = \frac{5\pi}{6}$, а не $-\frac\pi6$.
  3. Путают порядок: пишут $\int dx\int\ldots dy$, а внутренние пределы выражены через $y$. Внутренние пределы зависят только от внешней переменной.
  4. При смене порядка переставляют пределы механически, без рисунка.
  5. Не разбивают область, когда кривая входа или выхода меняется (примеры 8, 9, 12, 13 и РГР).
  6. Путают нижнюю и верхнюю границу: получается отрицательная площадь. Проверяй подстановкой одной промежуточной точки.
  7. Теряют половину области при решении относительно $x$: из $y = x^2$ получается $x = \pm\sqrt y$, а не только $\sqrt y$. И наоборот, $x = -\sqrt y$ из условия — это только левая половина параболы.
  8. $\arcsin$ для правой половины арки: $y = \sin x$ на $[\frac\pi2,\pi]$ обращается как $x = \pi - \arcsin y$.
  9. Путают $M_x$ и $M_y$: $M_x = \iint y\rho\,dS$ (плечо до оси $Ox$ — это $y$).
  10. Объём как $\iint z_2\,dS$, когда тело ограничено снизу не плоскостью $z = 0$: нужно $\iint (z_2 - z_1)\,dS$, и область $D$ — проекция линии пересечения поверхностей.

Мини-шпаргалка главы

Смысл: $\iint_D f\,dS = \lim\limits_{\lambda\to0}\sum f(P_i)\Delta S_i$; $f\ge0$ — объём под поверхностью, $f = 1$ — площадь, $f = \rho$ — масса.

Мера Жордана: площадь есть, если граница имеет площадь нуль. Интегрируемы непрерывные функции на замкнутых измеримых областях.

Дарбу: $s_T = \sum m_i\Delta S_i\le\sigma\le S_T = \sum M_i\Delta S_i$; интегрируема $\Leftrightarrow S_T - s_T\to0$.

О среднем: $\iint_D f\,dS = f(P_0)\,S(D)$; оценка $mS\le I\le MS$.

I тип: $\displaystyle\int_a^b dx\int_{\varphi_1(x)}^{\varphi_2(x)}f\,dy$ — вертикаль: вход снизу, выход сверху.

II тип: $\displaystyle\int_c^d dy\int_{\psi_1(y)}^{\psi_2(y)}f\,dx$ — горизонталь: вход слева, выход справа.

Смена порядка: пределы → неравенства → рисунок → новые пределы; граница меняется — разбей. Проверка площадью.

Замена: $dx\,dy = |J|\,du\,dv$, $J = \begin{vmatrix}x'_u & x'_v\\ y'_u & y'_v\end{vmatrix}$, $\frac{\partial(x,y)}{\partial(u,v)} = 1\big/\frac{\partial(u,v)}{\partial(x,y)}$.

Полярные: $dx\,dy = r\,dr\,d\varphi$. $x^2+y^2=2ax \to r = 2a\cos\varphi$, $\varphi\in[-\frac\pi2,\frac\pi2]$; $x^2+y^2 = 2by\to r = 2b\sin\varphi$, $\varphi\in[0,\pi]$; $x = a\to r = \frac{a}{\cos\varphi}$.

Обобщённые полярные: $x = ar\cos\varphi$, $y = br\sin\varphi$, $dx\,dy = ab\,r\,dr\,d\varphi$, эллипс $\to r\le1$.

Несобственные: предел по исчерпывающим областям; $\iint_{\R^2}e^{-x^2-y^2} = \pi$, $\int_{-\infty}^{\infty}e^{-x^2}dx = \sqrt\pi$; $\iint\frac{dS}{r^\alpha}$: у нуля сходится при $\alpha\lt2$, на бесконечности при $\alpha\gt2$.

Приложения: $V = \iint_D(z_2 - z_1)\,dS$; $m = \iint\rho$; $M_x = \iint y\rho$, $M_y = \iint x\rho$; $x_c = \frac{M_y}{m}$, $y_c = \frac{M_x}{m}$; $I_x = \iint y^2\rho$, $I_y = \iint x^2\rho$, $I_O = I_x + I_y$; Штейнер $I_l = I_c + md^2$.

Полезные интегралы: $\int_0^{\pi/2}\cos^2 = \frac\pi4$, $\int_0^{\pi/2}\cos^4 = \frac{3\pi}{16}$, $\int_{-\pi/2}^{\pi/2}\cos^3 = \frac43$, $\int_0^{2\pi}\cos^2 = \pi$.