Математический анализ
Двойные интегралы
Как Pixar вычисляет свет, попадающий в один пиксель? Двойным интегралом по видимому телесному углу, вычисляемым миллиарды раз за кадр.
- Pixar RenderMan: каждый пиксель — двойной интеграл от входящей радиации по 2D-телесному углу
- Аэродинамика Boeing: подъёмная сила считается как двойной интеграл от давления по поверхности крыла
- Реконструкция МРТ: интенсивность сигнала интегрируется по дискам и кольцам внутри пациента
- Страхование катастроф: ожидаемые убытки — двойной интеграл от плотности ущерба по зоне риска
Предварительные знания
- Одномерный определённый интеграл
- Частные производные и многомерные области
- Полярные координаты
Суммы Римана в 2D и повторные интегралы
RenderMan от Pixar усредняет входящий свет по 2D-телесному углу для каждого пикселя. Это двойной интеграл, вычисляемый миллиарды раз за кадр. Математика, стоящая за Toy Story 4, — это та же конструкция: разбиваем область, суммируем высоту на площадь, переходим к пределу.
При f = 1 получаем площадь D. При f = rho (поверхностная плотность) — полную массу. При f = h(x,y) (высота над плоскостью) — объём под графиком z = h(x,y) над D.
Перед интегрированием нарисуйте D. Картинка подсказывает, по какой переменной интегрировать первой и какими будут пределы внутреннего интеграла. Алгебра без эскиза — главный источник ошибок в пределах интегрирования.
Вычислите двойной интеграл от f(x,y) = x + y по прямоугольнику R = [0,1] x [0,2].
integral_0^1 integral_0^2 (x+y) dy dx = integral_0^1 (2x + 2) dx = 1 + 2 = 3.
Области I и II типа и перемена порядка интегрирования
Алгоритм штриховки AutoCAD должен вычислять площади произвольных замкнутых областей. Он раскладывает их на куски I и II типа и запускает двойной интеграл на каждом. Правильный выбор порядка превращает 12-секундное вычисление в 0.4-секундное. Тот же ответ, иной путь.
Перемена порядка работает только для абсолютно интегрируемых f (теорема Фубини). Для осциллирующих интегралов с компенсацией порядок может иметь значение. В вводных курсах все f удовлетворяют условию, но осознавайте предположение.
Эвристика: если внутренний интеграл тяжёлый — нарисуйте D и попробуйте другой порядок. Примерно один из трёх учебных интегралов сводится к элементарному после перемены порядка.
Область D ограничена кривыми y = x и y = x^2 при 0 <= x <= 1. Какой повторный интеграл правильно считает площадь D?
На [0,1] выполняется x >= x^2, поэтому y идёт от x^2 (парабола, снизу) до x (прямая, сверху). Площадь = integral_0^1 (x - x^2) dx = 1/2 - 1/3 = 1/6.
Двойные интегралы в полярных координатах
МРТ-сканеры реконструируют изображение, интегрируя интенсивность сигнала по дискам и кольцам внутри пациента. В декартовых это утраивает время счёта по сравнению с полярными. Полярные координаты — естественный язык для любой задачи с круговой симметрией, а якобиан r — плата за вход.
Знаменитый интеграл Гаусса integral по R^2 от e^{-(x^2+y^2)} = pi — это полярно-координатный трюк, доказывающий integral по R от e^{-x^2} dx = sqrt(pi). Каждый учебник теории вероятностей повторяет его.
Переходите к полярным всегда, когда область или подынтегральная функция содержат x^2 + y^2, круги, секторы или кольца. Якобиан r часто сокращает противный 1/r, дополнительно упрощая дело.
Чем заменяется dA при переходе от декартовых координат к полярным?
Определитель якобиана полярного отображения (r, theta) -> (r cos theta, r sin theta) равен r. Поэтому dA = dx dy = r dr d theta.
Объём, масса и среднее значение
Аэродинамики Boeing интегрируют давление по поверхности крыла, чтобы посчитать подъёмную силу. NOAA интегрирует осадки по водосборам для прогноза наводнений. Актуарии страховых компаний интегрируют плотность убытков по зонам риска для ценообразования катастроф-облигаций. Три индустрии — одна конструкция: двойной интеграл от плотности по области.
Момент инерции относительно оси — тоже двойной интеграл: I_x = двойной интеграл от y^2 * rho dA. Инженеры-механики используют его для вычисления изгибной жёсткости произвольных сечений.
Когда реальная величина описывается как 'суммарное что-то, распределённое по области' — почти всегда перед вами двойной интеграл. Опознайте плотность, область, выберите координаты.
Пластина занимает единичный диск x^2 + y^2 <= 1 с плотностью rho(x,y) = x^2 + y^2. Какова её полная масса?
Переходим в полярные: rho * dA = r^2 * r dr d theta = r^3 dr d theta. Масса = integral_0^{2pi} d theta * integral_0^1 r^3 dr = 2 pi * 1/4 = pi/2.
Куда ведут двойные интегралы
Двойные интегралы — точка входа во всё многомерное интегрирование: тройные интегралы, поверхностные интегралы и великие теоремы векторного анализа выстроены на них.
- Тройные интегралы — Следующий шаг: интегрирование по объёмным областям
- Замена переменных — Обобщает трюк с полярным якобианом до произвольных координат
- Теорема Грина — Превращает двойной интеграл по области в криволинейный по её границе
- Плотности вероятностей — Совместные распределения интегрируются двойными интегралами
Итоги
- Двойной интеграл = предел сумм Римана: двойной интеграл от f dA по R = lim sum f(x*, y*) Delta x Delta y
- Теорема Фубини: на прямоугольнике (или для абсолютно интегрируемой f) повторный интеграл равен двойному, порядок не важен
- Область I типа: y между кривыми g_1(x) и g_2(x); II типа: x между h_1(y) и h_2(y). Меняйте порядок, если внутренний интеграл не берётся
- Полярные координаты: dA = r dr d theta. Применяйте при круговой симметрии; превращают гауссиану по R^2 в однострочное доказательство
- Приложения: объём = двойной интеграл от h dA; масса = двойной интеграл от rho dA; среднее значение = (1/area) * двойной интеграл от f dA
Связанные уроки
- calc-20-extrema-multi — Частные производные и многомерные области закладывают почву для интегрирования в 2D
- calc-22-triple-integrals — Тройные интегралы расширяют ту же конструкцию на одну размерность выше
- calc-23-coordinate-changes — Замена переменных позволяет брать интегралы, которые в декартовых координатах выглядят страшно