1. Уважаемые друзья!

    С 8 февраля 2018 года наш форум переходит в режим Элитарного Клуба.
    Теперь незарегистрированным посетителям запрещено подглядывать и подслушивать наши тайные переговоры, а чтобы зарегистрироваться, нужно... впрочем, если вы действительно достойны стать членом Клуба, то вы наверняка разберётесь, как это сделать.

    Возрадуйтесь, обладатели зарегистрированных аккаунтов! Обещаем вам чистки, репрессии и все остальные бонусы тоталитарного сообщества.

    Всегда ваша,
    Администрация Корума

МА. Вычисление двойного интеграла

Тема в разделе "Вопросы высшей математики", создана пользователем Schufter, 21 июл 2013.

Модераторы: onyx
  1. Schufter

    Schufter Мизантроп

    Теоретический минимум

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

    Определённый интеграл
    Интеграл функции [​IMG], вычисляемый по отрезку [​IMG]. Отрезок разбивается на малые необязательно равные частичные отрезки
    длиной [​IMG], внутри каждого такого отрезка выбирается точка [​IMG], и составляется интегральная сумма [​IMG].
    Затем осуществляется переход к пределу разбиения отрезка на бесконечно малые частичные отрезки, а именно требуют, чтобы длина самого
    большого частичного отрезка стремилась к нулю. Величина предела и является интегралом [​IMG].

    Двойной интеграл
    Интеграл функции [​IMG], вычисляемый по двумерной области [​IMG]. Область разбивается на малые части. Эти малые части
    характеризуются их площадями [​IMG] и длинами [​IMG] наибольшего отрезка, который можно провести внутри частичной области. Этот
    отрезок называется диаметром частичной области. В каждой частичной области выбирается точка [​IMG], и составляется интегральная сумма
    [​IMG]. Затем осуществляется переход к пределу разбиения области [​IMG] на бесконечно малые частичные области,
    а именно требуют, чтобы длина самого большого диаметра частичной области стремилась к нулю. Величина предела и является интегралом
    [​IMG].

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

    Примеры вычисления двойных интегралов сведением к повторному интегрированию

    Пример 1. Прямоугольная область интегрирования.
    Вычислить интеграл [​IMG] по области [​IMG].
    Это самый простой пример, однако он же и основополагающий с точки зрения теории, как станет ясно потом. Фактически здесь нужно сначала
    вычислить интеграл по одной из переменных, например, по переменной [​IMG], считая вторую переменную константой (есть некоторая
    аналогия с вычислением частных производных). Затем выполняется и второе интегрирование.
    [​IMG].

    После вычисления зададимся вопросом: имел ли значение порядок интегрирования? Ответ: если оба повторных интеграла существуют, то порядок
    интегрирования не важен, и его следует выбирать таким, чтобы вычисление проводилось по возможности удобнее.

    Пример 2. Область интегрирования не является прямоугольной.
    Вычислить интеграл [​IMG] по области, ограниченной параболой [​IMG] и прямой [​IMG].

    Сделав чертёж (см. рис. 1), оценим, в каком порядке легче интегрировать. Прежде всего, поясним, что мы планируем придерживаться прежнего плана
    действий: сначала проинтегрировать по одной переменной, считая вторую постоянной, а затем проинтегрировать по второй переменной. И тут возникает
    вопрос о расстановке пределов интегрирования. При интегрировании по прямоугольнику такого вопроса не было. Обратимся к построению двойного интеграла.
    Основным элементом этого построения является разбиение на малые (в пределе бесконечно малые) части области интегрирования. Мы умеем выполнять
    интегрирование по прямоугольнику. Разобьём область интегрирования на малые прямоугольники. Если они будут бесконечно малыми, то они покроют область
    с криволинейной границей полностью. Двойной интеграл обладает свойством аддитивности. Если область интегрирования разбить на две части, то интеграл
    будет равен сумме интегралов по полученным областям. Это верно для любого числа областей. Теперь посмотрим на такое разбиение с геометрической точки
    зрения. Если мы фиксируем, например, переменную [​IMG], считая её постоянной, то мы выделяем набор прямоугольников, выстроенных по прямой, параллельной
    оси абсцисс (см. рис. 2). Тогда понятно в каких пределах следует интегрировать по переменной [​IMG] при данном значении [​IMG]. Но при другом
    значении [​IMG] пределы интегрирования по переменной [​IMG] будут уже другими. Это значит, что эти пределы интегрирования являются функциями [​IMG].
    Т.е. интегрирование можно выполнять по формуле Ньютона-Лейбница, но пределы интегрирования будут переменными. Второе интегрирование выполняется
    уже стандартно, пределы определяются всеми возможными значениями [​IMG], исходя из области интегрирования.
    [​IMG]

    Применим эти рассуждения к нашему конкретному примеру. Допустим, мы хотим, чтобы внутреннее интегрирование было по переменной [​IMG]. Это значит,
    что во внутреннем интегрировании переменная [​IMG] играет роль константы. Постоянным значениям [​IMG] соответствуют горизонтальные прямые (теперь уже,
    конечно, нет рядов прямогольников). Проведём несколько таких прямых (см. рис. 3). Видно, что где бы ни была проведена прямая, она всегда «входит» в кривую
    [​IMG], а «выходит» из прямой [​IMG]. Вот так мы и расставим пределы интегрирования во внутреннем интеграле.
    Осталось лишь определить, в каких пределах изменяется [​IMG]. Для этого нужно решить уравнение
    [​IMG].
    Расставляем пределы интегрирования и вычисляем интеграл:
    [​IMG]
    [​IMG].

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

    Пример 3. Разбиение области интегрирования на части.
    Вычислить интеграл [​IMG] по области, ограниченной параболой [​IMG] и прямыми [​IMG], [​IMG].
    [​IMG]

    С целью определения пределов интегрирования снова будем проводить прямые, пересекающие область интегрирования. Только прямые
    на этот раз будут вертикальными, а следовательно, внутреннее интегрирование будет по переменной [​IMG]. Выбран именно этот порядок
    интегрирования, так как из уравнения параболы неудобно выражать [​IMG]. На чертеже (см. рис. 4) показаны четыре из возможных вертикальных
    прямых. Как видно, все они «выходят» из параболы, но две «входят» через прямую [​IMG], а две – через прямую [​IMG].
    Это наводит на мысль, что нужно разбивать область интегрирования на две части. Граница этого разбиения – часть оси ординат. К слову,
    другой порядок интегрирования в этом смысле был бы ничуть не лучше: область интегрирования всё равно пришлось бы разбивать.
    [​IMG]
    [​IMG]
    [​IMG]
  2. Schufter

    Schufter Мизантроп

    Теоретический минимум (продолжение)

    Обсудим тему расстановки пределов интегрирования в полярных координатах. И для начала напомним, как вводятся полярные координаты.
    Вводится точка, играющая роль начала координат - полюс - и выходящий из неё луч - полярная ось. Положение точки на плоскости характеризуется
    расстоянием до полюса [​IMG] и углом [​IMG], который составляет с полярной осью отрезок, соединяющий точку и полюс (см. рис. 5).
    [​IMG]
    Связь декартовых прямоугольных и полярных координат задаётся соотношениями
    [​IMG]
    Отсюда следует, что
    [​IMG],
    т.е. уравнение [​IMG] задаёт на плоскости окружность радиуса [​IMG], а уравнение [​IMG] задаёт на плоскости
    прямую, проходящую через начало координат. Полезно изобразить сеть координатных линий (т.е. линий, на которых та или иная координата постоянна)
    для случая прямоугольных и полярных координат (см. рис. 6).

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

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

    Пример 4. Кольцевая область интегрирования.
    Расставить пределы интегрирования в полярных координатах в интеграле [​IMG], где область интегрирования –
    кольцо [​IMG].
    [​IMG]

    Этот пример аналогичен примеру 1 в том смысле, что это простейший случай интегрирования в полярных координатах. Границу области
    интегрирования образуют координатные линии полярной системы координат (см. рис. 7). Поэтому пределы интегрирования расставляются элементарно:
    [​IMG]
    Так как подынтегральная функция не конкретизирована, то за порядком интегрирования не следим: если повторные интегралы существуют,
    то они равны. При конкретных вычислениях, разумеется, выбирается тот порядок интегрирования, при котором вычисления проще.

    Пример 5. Квадратная область интегрирования.
    Расставить пределы интегрирования в полярных координатах в интеграле [​IMG], где область интегрирования –
    квадрат [​IMG].

    При интегрировании в прямоугольных координатах этот пример был простейшим. Теперь он несколько усложнился. Мы расставим пределы для обоих
    порядков интегрирования. Но сначала запишем уравнения сторон квадрата [​IMG] в полярных координатах:
    [​IMG].
    [​IMG]

    Первым рассмотрим случай, когда внутреннее интегрирование выполняется по переменной [​IMG], тогда полярный угол фиксирован,
    т.е. для определения пределов интегрирования следует проводить лучи, исходящие из начала координат (см. рис. 8). Полярный угол в целом будет
    изменяться от нуля до [​IMG]. Но заметьте, что прямые могут «выходить» из вертикальной прямой, а могут – из горизонтальной. В связи с этим
    интеграл распадётся на два:
    [​IMG]
    Теперь проведём интегрирование в другом порядке: внутреннее интегрирование проведём по полярному углу, тогда переменная [​IMG] считается
    константой. Для определения пределов интегрирования проводим дуги окружностей (см. рис. 9). Как видно, до некоторого значения радиуса окружности
    она пересекает стороны квадрата, соответствующие значениям полярного угла 0 (горизонтальная сторона) и [​IMG] (вертикальная сторона).
    Когда радиус окружности становится равным единице, то она начинает пересекать две другие стороны. Наконец, после достижения радиусом окружности
    значения [​IMG] окружность вообще перестаёт пересекать стороны квадрата. Таким образом, расстановка пределов интегрирования следующая:
    [​IMG]
  3. Schufter

    Schufter Мизантроп

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

    Наконец, рассмотрим замену переменных при вычислении двойного интеграла. Цель замены переменных та же, что и в случае определённого интеграла -
    сделать вычисление технически проще. Кроме того, при замене может заметно упроститься расстановка пределов интегрирования. Это уже можно было
    видеть на примерах 1 и 5, когда при интегрировании по одной и той же области в одних координатах расстановка пределов была элементарной, а в
    других потребовалось разбивать интеграл на две части.

    Проследим аналогию в выполнении замены в определённом и двойном интеграле. Напомним процедуру замены переменной в определённом интеграле.
    Пусть требуется в интеграле [​IMG] перейти к переменной интегрирования [​IMG], причём [​IMG]:
    [​IMG].
    Выделим три момента: 1) выполнение замены в подынтегральной функции; 2) преобразование дифференциала; 3) преобразование пределов интегрирования.
    Первый и третий пункты особых комментариев не требуют, остановимся подробнее на втором. В определённом интеграле он достаточно очевиден,
    а потому именно здесь удобно отследить его геометрический смысл. Переход от дифференциала старой переменной к дифференциалу новой переменной
    происходит с помощью производной. Геометрически это связано с тем, что некоторому малому приращению старой переменной соответствует вполне определённое
    приращение новой переменной в силу наличия между ними функциональной зависимости. Это показано на рис. 10. Производная, появляющаяся при замене,
    представляет собой тангенс угла наклона касательной к графику функции [​IMG], а потому она позволяет выразить один катет в изображённом
    прямоугольном треугольнике (приращение [​IMG]) через другой катет (приращение [​IMG]).
    [​IMG]

    Обратимся к случаю двойного интеграла и начнём с частного случая перехода к полярным координатам. Замена в подынтегральной функции выполняется
    элементарно, расстановка пределов интегрирования обсуждалась ранее. Остановимся на преобразовании дифференциалов. В случае определённого
    интеграла преобразование дифференциала представляло собой переход к дифференциалу новой переменной с поправкой на изменение длины отрезка, который
    сначала представлял собой приращение [​IMG], а потом представляет собой приращение [​IMG]. В случае двойного интеграла
    нужно рассмотреть преобразование элемента площади [​IMG] при замене переменных. Рассмотрим подробно, как образуется
    элемент площади в координатах [​IMG]. Строятся координатные линии [​IMG]. Они и
    ограничивают элемент площади. Как это выглядит в случае полярных координат, показано на рис. 11.
    [​IMG]

    Легко вычислить площадь изображённого элемента в полярных координатах. Так как речь идёт о бесконечно малых приращениях переменных
    [​IMG] и [​IMG], то изображённую фигуру можно считать прямоугольником. Его стороны равны [​IMG] (такова разница
    в радиусах дуг окружностей) и [​IMG] (дуга окружности радиуса [​IMG], на которую опирается центральный угол [​IMG]).
    Итак, элемент площади в полярных координатах [​IMG].

    Понятно, что такое вычисление может быть не столь простым для других систем координат, и вообще хотелось бы иметь общую методику замены
    переменных в двойном интеграле. Можно доказать, что общая методика такова: вместо дифференциалов [​IMG] следует записать
    дифференциалы [​IMG] и умножить их на модуль якобиана перехода
    [​IMG].
    Этот якобиан как раз и выражает изменение элемента площади, связанное с его деформацией при замене переменных.
    Например, в случае полярных координат
    [​IMG],
    что согласуется с результатом нашего расчёта, сделанного с помощью чертежа.

    Итак, замена переменных в двойном интеграле проводится следующим образом:
    [​IMG].
    Изменение обозначения области интегрирования с [​IMG] на [​IMG] связано с необходимостью изменить пределы интегрирования.
    Запомнить формулу достаточно просто: во-первых, по сути она повторяет формулу замены в определённом интеграле, а во-вторых, порядок
    перечисления переменных в якобиане такой, что новые переменные [​IMG] должны "сокращаться", чтобы оставались старые переменные [​IMG].
    Подчеркнём, что это исключительно мнемоническое правило, позволяющее лучше запомнить формулу. Никакой доказательной силы оно, конечно,
    не имеет. Ещё одно замечание касается появления в формуле модуля. При замене в определённом интеграле мы не ставили производную под
    знак модуля. Однако знак автоматически учитывала расстановка пределов интегрирования, ведь в определённом интеграле существует общая
    методика изменения пределов интегрирования, а в двойном интеграле их приходится расставлять заново.

    Примеры

    Пример 1. Площадь эллипса.
    Вычислить интеграл [​IMG] по области, ограниченной эллипсом [​IMG].

    В соответствии с геометрическим смыслом двойного интеграла мы сейчас фактически вычислим площадь эллипса. Проводить вычисление
    в прямоугольных координатах неудобно, и принципиально этот метод не отличается от нахождения площади фигуры с помощью определённого
    интеграла. Здесь удобно перейти к так называемым обобщённым полярным координатам:
    [​IMG].
    В этих координатах уравнение эллипса принимает предельно простой вид [​IMG], полярный угол принимает все значения от нуля до [​IMG].
    Якобиан такого преобразования отличается от случая полярных координат только дополнительным множителем [​IMG]. Это достаточно очевидно,
    но можно и проверить это расчётом. Итак,
    [​IMG].

    Пример 2. Вид замены диктуется областью интегрирования.
    Вычислить интеграл [​IMG] по области, ограниченной кривыми [​IMG].
    [​IMG]

    Область интегрирования показана на рис. 12. Понятно, что вычисление в прямоугольных координатах очень неудобно. Можно ввести такую
    замену координат, чтобы область интегрирования существенно упростилась. В этом примере вид замены подсказывается уравнениями,
    определяющими границы области интегрирования:
    [​IMG].
    После применения такой замены нужно будет интегрировать по области [​IMG], т.е. по квадрату.
    Остаётся вычислить якобиан перехода:
    [​IMG].
    Обратите внимание: был вычислен якобиан обратного перехода (здесь такое вычисление удобнее), т.е. при подстановке в интеграл нужно взять
    обратную величину по отношению к этому якобиану, так как [​IMG].
    [​IMG]
Модераторы: onyx

Поделиться этой страницей