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

Теорема (формула замены переменной в кратном интеграле)

Пусть отображение F:ΩRn, где ΩRn — открытое множество, заданное при помощи непрерывно дифференцируемых функций xi=ϕi(u1,,un),i=1,,n, является взаимно однозначным и удовлетворяет следующим условиям:

  1. производные ϕiui ограничены в Ω;
  2. производные ϕiui равномерно непрерывны в Ω;
  3. якобиан J(u) отображения удовлетворяет при uΩ условию |J(u)|α>0.

Тогда, если G — измеримый компакт с кусочно-гладкой границей, лежащий во множестве Ω и f(x) — непрерывна на множестве G=F(G), то справедлива следующая формула замены переменных в кратном интеграле:
Gf(x)dx=Gf(ϕ1(u),,ϕn(u))|J(u)|du(),


где x=(x1,,xn),u=(u1,,un).

Доказательство

Для начала рассмотрим еще 2 вспомогательных свойства:

  1. Если LΩ есть непрерывно дифференцируемая кривая, то ее образ L=F(L) есть непрерывно дифференцируемая кривая.
  2. Если G — область и ¯GΩ (где ¯G — замыкание области G), тогда ее образ G=F(G). Образ границы Ω есть граница Ω.

Первое свойство является простым следствием правила нахождения производной сложной функции, а второе — теоремы о неявных функциях.

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

Поскольку компакт G лежит в открытом множестве Ω, то границы этих множеств не пересекаются. Так как граница любого множества замкнута и граница ограниченного множества ограничена, то расстояние между границами множеств G и Ω есть положительное число δ.

Примечание №1: под разбиением множества A далее будем подразумевать совокупность измеримых множеств {A1,,An}, таких что A1An=A и AiAj=,ij. Клеткой назовем множество вида K={(x1,,xn)|aixi<bi,1in}, прямоугольником — клетку в пространстве R2.

Пусть P есть замкнутый квадрат, содержащий компакт G. Если разбить стороны квадрата P на равные части длины h<δ (чтобы отсутствовали квадраты, содержащие одновременно элементы границ G и Ω), то и сам квадрат P окажется разбит на квадратные клетки с площадью h2. Разбиение квадрата P порождает разбиение T компакта G. Если малый квадрат со стороной h целиком лежит внутри компакта G, то он является элементом разбиения T, а если он содержит граничные точки G, то соответствующим элементом разбиения является пересечение этого квадрата с компактом G. Отображение F порождает разбиение T компакта G=F(G), причем элементами разбиения T являются образы элементов разбиения T. При отбрасывании в интегральной сумме слагаемых, которым отвечают квадраты, имеющие непустое пересечение с множеством жордановой меры нуль, характер соответствующего предела при мелкости разбиения, стремящемся к нулю, не изменится (о чем свидетельствует соответствующая лемма, см. примечание №2). А значит, при написании интегральных сумм можно учитывать только слагаемые, соответствующие целым квадратам и их образам при отображении F, остальные квадраты будут иметь непустое пересечение с границей G. Так как отображение F равномерно непрерывно, то мелкость разбиения T стремится к нулю, когда стремится к нулю мелкость разбиения T.

Если малые квадраты P1,,Pn лежат внутри компакта G, то
их образы P1,,Pn лежат внутри G. Пусть (ui,vi) — координаты точки, лежащей в левом нижнем углу квадрата Pi, a (ϕ(ui,vi),ψ(ui,vi)) — образ этой точки при отображении F.

Тогда можем записать интегралы, входящие в формулу () как пределы интегральных сумм:
Gf(x,y)dxdy=limh0ni=1f(xi,yi,)m(Pi),
Gf(ϕ(u,v),ψ(u,v))|J(u,v)|dudv= limh0ni=1f(ϕ(ui,vi),ψ(ui,vi))|J(ui,vi)|m(Pi).

Для доказательства формулы () покажем, что разность этих интегральных сумм стремится к нулю при h0. В силу леммы о геометрическом смысле модуля якобиана отображения,
|m(Pi)|J(ui,vi)|m(Pi)|α(h)m(Pi),limh0α(h)=0.
Принимая во внимание, что ϕ(ui,vi)=xi,ψ(ui,vi)=yi,|f(x,y)|<M (последнее в силу того, что функция f непрерывна на компакте, а значит и ограниченна на нем), получаем оценку для разности интегральных сумм:
|ni=1f(xi,yi)m(Pi)ni=1f(ϕ(ui,vi),ψ(ui,vi))|J(ui,vi)|m(Pi)| ni=1|f(xi,yi)m(Pi)f(xi,yi)|J(ui,vi)|m(Pi)|= ni=1|f(xi,yi)|m(Pi)|J(ui,vi)|m(Pi)|| Mni=1α(h)m(Pi)= Mα(h)ni=1m(Pi) Mα(h)m(G), из которой следует, что эта разность стремится к нулю при h0 (т.к. M и m(G) — константы). Теорема доказана.

Примечание №2: о геометрическом смысле модуля якобиана отображения можно прочитать, например, в курсе лекций по мат. анализу В.И. Коляда, А.А. Кореновский (т.2, стр. 219) или в Тер-Крикоров А.М. и Шабунин М.И. «Курс математического анализа» (стр. 471). Лемма об отбрасывании слагаемых в интегральной сумме также присутствует и доказана, например, в учебнике Тер-Крикорова, стр. 458.

Замечание

Примеры

Основными примерами использования данной формулы являются переход к полярным, цилиндрическим и сферическим координатам для вычисления двойных и тройных интегралов.

Тест: формула замены переменной в кратном интеграле

Для закрепления материала, рекомендуется пройти тест по данной теме.


Таблица лучших: Замена переменной в кратных интегралах

максимум из 3 баллов
Место Имя Записано Баллы Результат
Таблица загружается
Нет данных

Использование полярных, цилиндрических и сферических координат для вычисления кратных интегралов

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

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

Из курса аналитической геометрии известны следующие соотношения между декартовыми и полярными координатами: x=rcosϕ,y=rsinϕ().
При этом, r0,0ϕ<2π. Рассмотрим вспомогательную плоскость ROΦ, где r и ϕ являются декартовыми координатами, и определим на ней множество точек G, такое, что: G={(r,ϕ)|r>0,0ϕ<2π}.

Тогда формулы () определяют непрерывно дифференцируемое отображение F:G~XOY, где ~XOY=XOY{(0,0)}.

По определению полярных координат, в декартовой системе координат XOY r задает радиус окружности с центром в начале координат, а ϕ определяет луч, исходящий из центра координат, такой что угол между лучом и положительным направлением оси OX равен ϕ. С геометрической точки зрения очевидно, что они пересекаются в единственной точке.

Таким образом, любую точку P=(x0,y0) из ~XOY можно однозначно определить пересечением луча, направленного под углом ϕ0 и окружности радиусом r0, и тогда точка P=(r0,ϕ0) будет единственным прообразом P в G. Очевидно, что любой элемент из G служит прообразом, и что двум различным точкам из G будут соответствовать 2 различные точки из ~XOY. Таким образом, отображение F между точками плоскостей G и ~XOY взаимно однозначное:

kursach

Якобиан полученного отображения будет равен:
JF=xrxϕyryϕ=cosϕrsinϕsinϕrcosϕ=r

Теперь рассмотрим множество точек G, полученное добавлением к множеству G отрезка r=0, т.е. G={(r,ϕ)|r0,0ϕ<2π}. G уже является прообразом всей плоскости XOY, но на отрезке r=0,0ϕ<2π не достигается взаимная однозначность, а |JF|=0. Обратим внимание, что его Жорданова мера равна нулю.

Наконец, пусть дана область ΩXOY и функция f, непрерывная на измеримом множестве ¯Ω. Ее прообразом при отображении F, заданного формулами (), будет некоторая область ΩG. Если область Ω не содержит точки O — начала координат, то выполнены все условия теоремы о замене переменной в кратных интегралах, и справедлива формула:
Ωf(x,y)dxdy=Ωf(rcosϕ,rsinϕ)rdrdϕ


Если же точка OΩ, то взаимная однозначность и не обращение якобиана в нуль не выполняются на множестве r=0, что не влияет на справедливость данной формулы (следует из замечания к указанной теореме).

Пример №1

Использование цилиндрических и сферических координат

Рассмотрим теперь пространство R3, в котором задана декартова система координат OXYZ. Цилиндрические координаты связанны с декартовыми следующим образом:
x=rcosϕ,y=rsinϕ,z=t(),
где r0,0ϕ<2π,tR (величины r и ϕ для любой точки A=(x,y,z) определяются таким же образом, как и в полярных координатах для ее проекции P=(x,y,0) на XOY). Теперь, аналогично случаю с полярными координатами, рассмотрим вспомогательное пространство ORΦT, где r,ϕ,t — декартовы координаты, а в нем — множество точек G={(r,ϕ,t)|r0,0ϕ<2π,tR}.

Отображение F:GOXYZ, определяемое формулами (), является непрерывно дифференцируемым.
JF=xrxϕxtyryϕytzrzϕzt=cosϕrsinϕ0sinϕrcosϕ0001=r

Очевидно, что как и в случае с полярными координатами, отображение F — взаимно однозначное, и его якобиан не равен нулю. Данные условия не выполняются только при r=0, т.е. на множестве L={(r,ϕ,t)|r=0,0ϕ<2π,tR}. Пересечение такого множества с любым другим ограниченным множеством есть ограниченное линейное множество, и жорданова мера этого пересечения равна нулю.

Тогда, если дана область ΩOXYZ, и функция f непрерывна на измеримом множестве ¯Ω, а ΩG — прообраз данной области при отображении F, то выполнены все условия теоремы о замене, и справедлива следующая формула:
Ωf(x,y,z)dxdydz=Ωf(rcosϕ,rsinϕ,t)rdrdϕdt

Наконец, рассмотрим сферические координаты, связанные с декартовыми следующими соотношениями: x=rcosϕcosψ,y=rsinϕcosψ,z=rsinψ(),
где r0,0ϕ<2π,π2ψπ2. Введем вспомогательное пространство ORΦΨ, где r,ϕ,ψ — декартовы координаты, а в нем рассмотрим множество точек G={(r,ϕ,ψ)|r0,0ϕ<2ππ2ψπ2}.

Отображение F:GOXYZ, определяемое формулами (), непрерывно дифференцируемо.
JF=xrxϕxψyryϕyψzrzϕzψ=cosϕcosψrsinϕcosψrcosϕsinψsinϕcosψrcosϕcosψrsinϕsinψsinψ0rcosψ= r2cosψ.

Взаимная однозначность данного отображения устанавливается по тем же рассуждениям, что и в предыдущих двух случаях, и не выполняется только при r=0,ψ=π2,ψ=π2, когда и якобиан равен нулю. Однако любое подмножество множества, задаваемого такими равенствами, будет представлять собой ограниченную часть плоскости с жордановой мерой нуль в пространстве OXYZ, что не помешает совершить замену.

Тогда, при соответствующих условиях, справедлива формула замены переменной (ΩOXYZ,ΩG):
Ωf(x,y,z)dxdydz=

Ωf(rcosϕcosψ,rsinϕcosψ,rsinψ)r2cosψdrdϕdψ

Пример №2

Тест: Использование полярных, цилиндрических и сферических координат для вычисления кратных интегралов

Для закрепления материала, рекомендуется пройти тест по данной теме.


Таблица лучших: Переход к полярным, цилиндрическим и сферическим координатам при вычислении кратных интегралов

максимум из 7 баллов
Место Имя Записано Баллы Результат
Таблица загружается
Нет данных

M1630. О количестве способов представления натурального числа в виде суммы степеней двойки

Условие

nN обозначим через f(n) число способов представления числа n в виде суммы целых неотрицательных степеней числа 2. Преставления, отличающиеся лишь порядком слагаемых, считаются одинаковыми. Например, f(4)=4, так как число 4 может быть представлено следующими четырьмя способами:
4=4;
4=2+2;
4=2+1+1;
4=1+1+1+1.
Докажите, что nN,n3:2n2/4<f(2n)<2n2/2.

Решение

Если n=2k+1 — любое нечетное число, большее 1, то каждое его представление в требуемом по условию задачи виде содержит 1 в качестве слагаемого. Убрав эту единицу, мы получим представление числа 2k. Очевидно и обратное. Следовательно, f(2k+1)=f(2k). ()

Если n=2k — любое четное число, то каждое его представление в требуемом виде принадлежит к одному из двух типов: либо оно содержит слагаемое 1, либо не
содержит таких слагаемых. В первом случае, убрав одно слагаемое 1, мы получим представление числа 2k1. Как и выше, легко заметить, что есть взаимно однозначное соответствие между всеми представлениями числа 2k1 и представлениями числа 2k первого типа. Во втором случае мы можем разделить все слагаемые
на 2 и получить представление числа k. Это соответствие также взаимно однозначно. Итак, f(2k)=f(2k1)+f(k). ()

Обе полученные формулы выполнены для всех натуральных k>1. Очевидно, что f(1)=1. Пусть по определению f(0)=1, тогда формула () выполнена и при k=0. Заметим еще, что из () и () следует, что f(n) не убывает.

Согласно (), число f(2k1) в () можно заменить на
f(2k2), откуда f(2k)f(2k2)=f(k) k1. Суммируя эти равенства от 1 до n, получаем, что f(2n)=ni=0f(i). ()
В правой части () каждое слагаемое не больше последнего, а так как 2=f(0)+f(1)f(n) n2, то f(2n)=2+ni=2f(i)2+(n1)f(n)f(n)+(n1)f(n)nf(n) n2. Следовательно, f(2n)2n1f(2n1)2n12n2f(2n2)2(n1)+(n2)++2+1f(2)=2n(n1)2+1.
И так как n(n1)2+1<n22, т.е. 2n(n1)2+1<2n22 n3, верхняя оценка для f(2n) получена.

Чтобы получить нижнюю оценку, докажем сначала, что f(b+1)f(b)f(a+1)f(a) (4), если 0ab — целые числа одинаковой четности. Действительно, если a и b четны, то из () следует,что каждая часть неравенства (4) обращается в нуль, если они оба нечетны, то из () следует, что f(b+1)f(b)=f(b+12), f(a+1)f(a)=f(a+12), при чем, как было показано ранее, f не убывает.

Возьмем произвольные r и k, rk, r — четное. Подставим в неравенство (4) a=r+j,b=rj,j=0,1,,k1. Получаем следующий набор неравенств:
j=0:f(r+1)_f(r)f(r+1)f(r)_
j=1:f(r+2)_f(r+1)_f(r)_f(r1)_
j=2:f(r+3)_f(r+2)_f(r1)_f(r2)_
(…)
j=k2:f(r+k1)_f(r+k2)_f(r+k3)_f(r+k2)_
j=k1:f(r+k)f(r+k1)_f(r+k2)_f(r+k1)
Просуммировав левые и правые части неравенств соответственно, легко увидеть, что все промежуточные слагаемые сокращаются, и остается следующее: f(r+k)f(r)f(r+1)f(rk+1), а в силу того, что r — четное и f(r+1)=f(r), получим f(r+k)+f(rk+1)2f(r). Так как rk, просуммируем уже такие неравенства по k=1,2,,r, и получим 2ri=1f(i)2rf(r), или, что тоже самое 2ri=0f(i)f(0)2rf(r). В силу () 2ri=0f(i)=f(4r), а f(0)=1. Получаем следующее неравенство: f(4r)2rf(r)+1>2rf(r). Возьмем r=2m2, тогда f(2m)>2m1f(2m2). (5)
(Очевидно, что mN,m>2:2m2 — четно. Также, непосредственно подстановкой легко убедиться, что неравенство (5) справедливо и при m=2:f(4)=4>2=f(2)f(1).)

Наконец. пусть nN,n>1. Если l — положительное целое такое, что 2ln, то применяя неравенство (5) к m=n,n1,,n2l+2, получаем:
f(2n)>2n1f(2n2)>2n12n3f(2n4)> >2(n1)+(n3)++(n2l+1)f(2n2l)=2l(n2l)f(2n2l) (слагаемых в показателе степени — l штук, а 1+3++(2l1)=l2, тогда (n1)+(n3)++(n2l+1)=nll2=l(nl)).

Теперь, пусть n — четно, l=n2. Получаем, f(2n)>2n2(nn2)f(20)=2n24f(1)=2n24.
Если же n — нечетно, то положим l=n12. Тогда, f(2n)>2n12(nn12)f(21)=2(n1)(n+1)4f(2)=2n214+1=2n2+34>2n24.
Таким образом, нужный результат доказан n>1. Непосредственно проверяется, что и для n=1 соответствующее неравенство справедливо.

Тогда получаем, что nN,n3:2n2/4<f(2n)<2n2/2, что и требовалось доказать.