Тройной интеграл

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

Тройной интеграл часто вычисляется путём сведения его к повторному интегралу, что позволяет свести задачу к последовательности более простых одномерных интегралов. Введение объясняет основные принципы разбиения области интегрирования, формулирует теорему о переходе к повторному интегралу и определяет цели статьи: изучение методов вычисления и их практическое применение.

Будем называть ограниченную замкнутую область  V  простой правильной \mathbf { \textit { V } } \textbf { простой правильной }, если выполняются два условия : проекция  V \mathbf { \textit { V } } на какую-либо координатную плоскость, например, на плоскость  Оху \mathbf { \textit { Оху } } - простая область  D \mathbf { \textit { D } }, и любая прямая, перпендикулярная этой плоскости и проходящая через внутреннюю точку  V \mathbf { \textit { V } }, пересекает границу  V \mathbf { \textit { V } } в двух точках. Такую область можно описать следующим образом:

V=((x,y,z)∣(x,y)∈D,;ψ_1(x,y)⩽z⩽ψ_2(x,y))V=\left( { (x,y,z)\vert (x,y)\in D,;\psi \_1 (x,y)\leqslant z\leqslant \psi \_2 (x,y) }\right) поверхность z=ψ_1(x,y)z=\psi \_1 (x,y) образована множеством нижних точек пересечения прямой, параллельной оси Oz \mathbf { \textit { Oz } }, с границей  V \mathbf { \textit { V } }; поверхность z=ψ_2(x,y)z=\psi \_2 (x,y) - множеством верхних точек пересечения).

Теорема

Если  V \mathbf { \textit { V } } - простая область с кусочно-гладкой границей,  f ( P )\mathbf { \textit { f } } (\mathbf { \textit { P } } ) - непрерывная функция, то ∭_Vf(P)dv=∬_Ddxdy∫_ψ_1(x,y)ψ_2(x,y)f(x,y,z)dz\iiint\limits\_V { f(P)dv } =\iint\limits\_D { dxdy\int\limits\_ { \psi \_1 (x,y) } ^ { \psi \_2 (x,y) } { f(x,y,z)dz } }.

 Доказать \textbf { Доказать }эту теорему можно так, как мы доказали теорему о переходе от двойного интеграла к повторному: установить, что для повторного интеграла в правой части формулы имеют место все свойства интеграла, разбить область  V \mathbf { \textit { V } } на подобласти  V _i( i =1,2,…, n )\mathbf { \textit { V } } \_ { i } (\mathbf { \textit { i } } = 1,2, { \ldots } ,\mathbf { \textit { n } } ), пользуюсь свойствами аддитивности и теоремой о среднем, представить повторный интеграл как интегральную сумму для тройного (∑_i=1nf(P_i)v(V_i))\left( { \sum\limits\_ { i=1 } ^n { f(P\_i )v(V\_i ) } }\right) и перейти к пределу при d=max⁡_i=1,2,…,ndiam(V_i)→0d=\mathop { \max } \limits\_ { i=1,2,\ldots ,n } diam(V\_i )\to 0.

Если расписать двойной интеграл по простой области  D (D=[(x,y)∣a⩽x⩽b,φ_1(x)⩽y⩽φ_2(x)])\mathbf { \textit { D } } \quad \left( { D=\left[{ (x,y)\left| { \begin{array} { l } a\leqslant x\leqslant b, \\ \varphi \_1 (x)\leqslant y\leqslant \varphi \_2 (x) \\ \end{array} }\right. }\right] }\right) в виде повторного, получим ещё более детализированную формулу для вычисления тройного интеграла: ∭_Vf(P)dv=∬_Ddxdy∫_ψ_1(x,y)ψ_2(x,y)f(x,y,z)dz=∫_abdx∫_φ_1(x)φ_2(x)dy∫_ψ_1(x,y)ψ_2(x,y)f(x,y,z)dz\iiint\limits\_V { f(P)dv } =\iint\limits\_D { dxdy\int\limits\_ { \psi \_1 (x,y) } ^ { \psi \_2 (x,y) } { f(x,y,z)dz } } =\int\limits\_a^b { dx\int\limits\_ { \varphi \_1 (x) } ^ { \varphi \_2 (x) } { dy\int\limits\_ { \psi \_1 (x,y) } ^ { \psi \_2 (x,y) } { f(x,y,z)dz } } }.

Можно также доказать, что тройной интеграл можно представить в виде повторного с другим порядком интегрирования. Обозначим z_min⁡=min⁡z_(x,y,z)∈V,;z_max⁡=max⁡z_(x,y,z)∈Vz\_ { \min } =\mathop { \min z } \limits\_ { (x,y,z)\in V } ,;z\_ { \max } =\mathop { \max z } \limits\_ { (x,y,z)\in V } т.е. минимальное и максимальное значения ординаты для точек области  V )\mathbf { \textit { V } } ), D_zD\_z - плоскую область, получающуюся при сечении  V \mathbf { \textit { V } } плоскостью  z \mathbf { \textit { z } } = const. Тогда ∭_Vf(P)dv=∫_z_min⁡z_max⁡dz∬_D_zf(x,y,z)dxdy\iiint\limits\_V { f(P)dv } =\int\limits\_ { z\_ { \min } } ^ { z\_ { \max } } { dz\iint\limits\_ { D\_z } { f(x,y,z) } } dxdy. Естественно, для конкретной задачи может оказаться предпочтительней проектировать  V \mathbf { \textit { V } } не на плоскость  Оху \mathbf { \textit { Оху } }, а на другую координатную плоскость.

Примеры решения задач

Пример 1

I=∭_Vdxdydz(1+x+y+z)3,;V:[x=0,;y=0,;z=0,x+y+z=1.I=\iiint\limits\_V { \frac { dxdydz } { \left( { 1+x+y+z }\right)^3 } } ,;V:\left[{ \begin{array} { l } x=0,;y=0,;z=0, \\ x+y+z=1. \\ \end{array} }\right.

Проекция области  V \mathbf { \textit { V } } на плоскость  Оху \mathbf { \textit { Оху } } - треугольник D:[x=0,;y=0,;x+y=1]D:\left[{ x=0,;y=0,;x+y=1 }\right], поэтому I=∭_Vdxdydz(1+x+y+z)3=∬_Ddxdy∫_01−x−ydz(1+x+y+z)3=∬_D[−12⋅1(1+x+y+z)2]∣_01−x−ydxdy=I=\iiint\limits\_V { \frac { dxdydz } { \left( { 1+x+y+z }\right)^3 } } =\iint\limits\_D { dxdy\int\limits\_0^ { 1-x-y } { \frac { dz } { (1+x+y+z)^3 } } } =\iint\limits\_D { \left. { \left[ { -\frac { 1 } { 2 } \cdot \frac { 1 } { (1+x+y+z)^2 } }\right] }\right|\_0^ { 1-x-y } dxdy } =

=−12∫_01dx∫_01−x[14−1(1+x+y)2]dy=−12∫_01[y4+11+x+y]∣_01−xdx=−12∫_01[1−x4+12−11+x]dx==−12∫_01[1−x4+12−11+x]dx=−12(3x4−x28−ln⁡(1+x))_01=−12(34−18−ln⁡2)=12(ln⁡2−58)=-\frac { 1 } { 2 } \int\limits\_0^1 { dx } \int\limits\_0^ { 1-x } { \left[ { \frac { 1 } { 4 } -\frac { 1 } { (1+x+y)^2 } }\right]dy } =-\frac { 1 } { 2 } \int\limits\_0^1 { \left. { \left[ { \frac { y } { 4 } +\frac { 1 } { 1+x+y } }\right] }\right|\_0^ { 1-x } dx } =-\frac { 1 } { 2 } \int\limits\_0^1 { \left[ { \frac { 1-x } { 4 } +\frac { 1 } { 2 } -\frac { 1 } { 1+x } }\right]dx } = \\ =-\frac { 1 } { 2 } \int\limits\_0^1 { \left[ { \frac { 1-x } { 4 } +\frac { 1 } { 2 } -\frac { 1 } { 1+x } }\right]dx } =-\frac { 1 } { 2 } \left( { \frac { 3x } { 4 } -\frac { x^2 } { 8 } -\ln (1+x) }\right)\_0^1 =-\frac { 1 } { 2 } \left( { \frac { 3 } { 4 } -\frac { 1 } { 8 } -\ln 2 }\right)=\frac { 1 } { 2 } \left( { \ln 2-\frac { 5 } { 8 } }\right)

Пример 2

I=∭_Vz⋅dxdydz,;V:[z2=h2R2(x2+y2),z=h.I=\iiint\limits\_V { z\cdot dxdydz } ,;V:\left[{ \begin{array} { l } z^2=\frac { h^2 } { R^2 } \left( { x^2+y^2 }\right), \\ z=h. \\ \end{array} }\right. Здесь  V \mathbf { \textit { V } } - внутренность конуса,  D \mathbf { \textit { D } } - проекция круга, получающегося при сечении этого конуса плоскостью  z = h \mathbf { \textit { z } } =\mathbf { \textit { h } } на  Оху \mathbf { \textit { Оху } }, т.е. круг, ограниченный кривой h2=h2R2(x2+y2)⇒x2+y2=R2h^2=\frac { h^2 } { R^2 } \left( { x^2+y^2 }\right)\Rightarrow x^2+y^2=R^2, поэтому I=∬_Ddxdy∫_hRx2+y2hzdz=∬_Dz22∣_hRx2+y2hdxdy=12∬_D[h2−h2R2(x2+y2)]dxdy=I=\iint\limits\_D { dxdy\int\limits\_ { \frac { h } { R } \sqrt { x^2+y^2 } } ^h { zdz } = } \iint\limits\_D { \left. { \frac { z^2 } { 2 } }\right|\_ { \frac { h } { R } \sqrt { x^2+y^2 } } ^h dxdy= } \frac { 1 } { 2 } \iint\limits\_D { \left[ { h^2-\frac { h^2 } { R^2 } \left( { x^2+y^2 }\right) }\right]dxdy } =

=h22R2∬_D[R2−(x2+y2)]dxdy==\frac { h^2 } { 2R^2 } \iint\limits\_D { \left[ { R^2-\left( { x^2+y^2 }\right) }\right]dxdy } = { переходим к полярным координатам } =h22R2∫_02πdφ∫_0R[R2−r2]rdr=h22R2∫_02π[R2r22−r44]_0Rdφ=h22R22πR44=πh2R24=\frac { h^2 } { 2R^2 } \int\limits\_0^ { 2\pi } { d\varphi } \int\limits\_0^R { \left[ { R^2-r^2 }\right]rdr } = \frac { h^2 } { 2R^2 } \int\limits\_0^ { 2\pi } { \left[ { R^2\frac { r^2 } { 2 } -\frac { r^4 } { 4 } }\right]\_0^R d\varphi } =\frac { h^2 } { 2R^2 } 2\pi \frac { R^4 } { 4 } =\frac { \pi h^2R^2 } { 4 }.

Вычислим тот же интеграл по другой формуле перехода к повторному интегралу: I=∫_0hdz∬_D_zzdxdy=∫_0hzdz∬_D_zzdxdy=∫_0hzdz∬_D_zdxdy=∫_0hzdz⋅s(D_z)=I=\int\limits\_0^h { dz\iint\limits\_ { D\_z } { zdxdy } } =\int\limits\_0^h { zdz\iint\limits\_ { D\_z } { zdxdy } } =\int\limits\_0^h { zdz\iint\limits\_ { D\_z } { dxdy } } =\int\limits\_0^h { zdz\cdot s(D\_z ) } = внутренний двойной интеграл - интеграл от функции, равной 1, поэтому он равен площади круга, получающегося при сечении конуса плоскостью z=constz=const, уравнение ограничивающей окружности, площадь s(D_z)=πz2R2h2)=∫_0hzπz2R2h2dz=πR2h2⋅z44∣_0h=πh2R24s(D\_z )=\frac { \pi z^2R^2 } { h^2 } )=\int\limits\_0^h { z\frac { \pi z^2R^2 } { h^2 } dz } =\frac { \pi R^2 } { h^2 } \cdot \left. { \frac { z^4 } { 4 } }\right|\_0^h =\frac { \pi h^2R^2 } { 4 }.

Это решение оказалось проще; мы сыграли на том, что подынтегральная функция не зависит от  х \mathbf { \textit { х } } и  у \mathbf { \textit { у } }.

Другие статьи из категории

Похожие статьи

Все статьи →

Тройной интеграл

Определение тройного интеграла. Теорема существования тройного интеграла

Тройной интеграл обобщает понятие двойного интеграла на случай трёхмерного пространства и является важным инструментом математического…

Тройной интеграл

Свойства тройного интеграла: анализ основных характеристик и их математическое обоснование

Тройной интеграл обладает рядом важных свойств, которые делают его удобным инструментом для решения задач анализа, геометрии и физики в…

Тройной интеграл

Замена переменных в тройном интеграле

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

Тройной интеграл

Примеры применения цилиндрических и сферических координат: практические задачи и их решения

Цилиндрические и сферические координаты являются мощным инструментом для работы с трёхмерными областями, обладающими определённой…

Тройной интеграл

Механические приложения тройного интеграла: масса, центр масс и моменты инерции тел

Тройной интеграл является мощным инструментом для решения задач механики, связанных с распределением массы в трёхмерных телах. Введение…

Нужен расчёт?

Пришлите исходные данные — уточним состав работы и ответим по срокам.

Связаться