Кратные интегралы и их вычисление

Зачем нужны кратные интегралы

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

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

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

Пусть f(x,y) определена и ограничена на замкнутой области D ⊂ ℝ². Разбиваем D на малые кусочки с площадями ΔAₖ, выбираем в каждом точку (ξₖ, ηₖ) и формируем интегральную сумму Σ f(ξₖ, ηₖ) ΔAₖ. Если этот предел существует при измельчении разбиения и не зависит от способа разбиения и выбора точек, он называется двойным интегралом:

∬_D f(x,y) dA = lim Σ f(ξₖ, ηₖ) ΔAₖ.

Геометрический смысл: если f(x,y) ≥ 0, то двойной интеграл равен объёму под поверхностью z = f(x,y) над областью D. При f = 1 получаем площадь области: S(D) = ∬_D dA.

Физический смысл: если ρ(x,y) — поверхностная плотность (кг/м²), то масса тела M = ∬_D ρ(x,y) dA. Координаты центра масс: x̄ = (1/M) ∬_D x ρ dA, ȳ = (1/M) ∬_D y ρ dA. Момент инерции относительно начала координат: I₀ = ∬_D (x² + y²) ρ dA.

Теорема Фубини: сведение к повторному интегралу

Если область D является «правильной» (x-простой или y-простой) и f непрерывна на D, то двойной интеграл сводится к последовательным обычным интегралам.

Для x-простой области D = {(x,y): a ≤ x ≤ b, φ₁(x) ≤ y ≤ φ₂(x)}:

∬D f dA = ∫ₐᵇ [∫{φ₁(x)}^{φ₂(x)} f(x,y) dy] dx.

Смысл: фиксируем x (вертикальный «слой»), интегрируем по y внутри слоя, затем суммируем все слои по x.

Развёрнутый пример. Вычислим ∬_D (x² + y²) dA, где D — треугольник с вершинами (0,0), (1,0), (1,1).

Описание области: 0 ≤ x ≤ 1, 0 ≤ y ≤ x.

∫₀¹ [∫₀ˣ (x² + y²) dy] dx = ∫₀¹ [x²y + y³/3]₀ˣ dx = ∫₀¹ (x³ + x³/3) dx = ∫₀¹ (4x³/3) dx = [x⁴/3]₀¹ = 1/3.

Смена порядка интегрирования. Иногда интеграл в одном порядке не берётся аналитически, но легко вычисляется в другом. Тот же треугольник как y-простая область: 0 ≤ y ≤ 1, y ≤ x ≤ 1.

∫₀¹ [∫_y¹ (x² + y²) dx] dy = ∫₀¹ [x³/3 + y²x]_y¹ dy = ∫₀¹ (1/3 + y² - y³/3 - y³) dy = 1/3 + 1/3 - 1/12 - 1/4 = 1/3. ✓

Замена переменных и якобиан

При замене (x,y) = Φ(u,v) формула имеет вид:

∬D f(x,y) dx dy = ∬{D'} f(x(u,v), y(u,v)) |J| du dv,

где якобиан J = det(∂(x,y)/∂(u,v)) = ∂x/∂u · ∂y/∂v − ∂x/∂v · ∂y/∂u — коэффициент «растяжения» площади при отображении. Абсолютная величина |J| обеспечивает положительность площади.

Полярные координаты: x = r cos θ, y = r sin θ. Якобиан: J = r.

∬D f(x,y) dx dy = ∬{D'} f(r cos θ, r sin θ) · r dr dθ.

Множитель r перед dr dθ — «цена» полярного элемента площади: маленький сектор на расстоянии r от начала имеет площадь r·Δr·Δθ.

Классический пример: Гауссов интеграл ∫_{-∞}^{∞} e^{-x²} dx = √π.

Вычислим I² = (∫{-∞}^{∞} e^{-x²} dx)(∫{-∞}^{∞} e^{-y²} dy) = ∬_{ℝ²} e^{-(x²+y²)} dA.

В полярных координатах:

I² = ∫₀^{2π} ∫₀^{∞} e^{-r²} r dr dθ = 2π · [-e^{-r²}/2]₀^{∞} = 2π · 1/2 = π.

Значит, I = √π — результат, фундаментальный для теории вероятностей и статистики.

Тройные интегралы и криволинейные координаты

Тройной интеграл ∭_V f(x,y,z) dV определяется аналогично. Для вычисления используют:

Цилиндрические координаты: x = r cos θ, y = r sin θ, z = z; dV = r dr dθ dz. Удобны для тел с осью симметрии (цилиндры, конусы).

Сферические координаты: x = ρ sin φ cos θ, y = ρ sin φ sin θ, z = ρ cos φ; dV = ρ² sin φ dρ dφ dθ. Удобны для тел с центральной симметрией (шары, части шаров).

Объём шара радиуса R:

∭_{ρ≤R} dV = ∫₀^{2π} ∫₀^π ∫₀^R ρ² sin φ dρ dφ dθ = 2π · 2 · R³/3 = 4πR³/3. ✓

Момент инерции шара относительно оси z при равномерной плотности ρ₀:

I_z = ρ₀ ∭_{ρ≤R} (x² + y²) dV = ρ₀ ∫₀^{2π} ∫₀^π ∫₀^R ρ² sin²φ · ρ² sin φ dρ dφ dθ = 2MR²/5,

где M — масса шара. Это формула, используемая в механике для расчёта кинетической энергии вращения.

Вероятностные приложения

Двойные и тройные интегралы — язык многомерных вероятностных распределений. Если (X, Y) — непрерывный случайный вектор с плотностью p(x,y), то:

P((X,Y) ∈ D) = ∬_D p(x,y) dA.

Маргинальная плотность X: p_X(x) = ∫_{-∞}^{∞} p(x,y) dy — интегрирование «по ненужной» переменной. Совместное нормальное распределение, распределения Дирихле, вычисление нормировочных констант — всё это задачи на кратные интегралы.

Вопрос для размышления: задача вычисления кратного интеграла по сложной области часто сводится к смене порядка интегрирования или замене координат. Попробуйте вычислить ∬_D y e^{x²} dA на треугольнике D: 0 ≤ y ≤ x ≤ 1, сначала в порядке dy dx, а затем меняя порядок.