Vizualizare sursă pentru „Integrală multiplă”

[[Fișier:Volume under surface.png|right|thumb|O integrală dublă interpretată ca volumul de sub o suprafaţă. Regiunea dreptunghiulară de la baza corpului este domeniul de integrare, iar suprafaţa este graficul funcţiei de două variabile care este integrată.]]
Noțiunea de '''integrală multiplă''' este similară cu noțiunea de [[Integrală|integrală definită]], extinsă la [[funcție de mai multe variabile reale|funcții de mai multe variabile reale]], de exemplu, <math>f(x,y)</math> sau <math>f(x,y,z)</math>.

Așa cum integrala definită a unei funcții pozitive de o singură variabilă reprezintă aria suprafeței dintre graficul funcției și axa ''x'', '''integrala dublă''' a unei funcții pozitive de două variabile reprezintă [[volum]]ul regiunii de spațiu aflată între graficul funcției și planul care conține [[Domeniu (matematică)|domeniul de definiție]] al acesteia. (Același volum poate fi obținut prin calculul '''integralei triple''' — integrala unei funcții de trei variabile — a funcției constante ''f''(''x'', ''y'', ''z'') = 1 pe regiunea sus-menționată, dintre suprafață și plan.) Dacă există mai multe variabile, o integrală multiplă va da [[hipervolum]]ul unei funcții de mai multe variabile.

Integrala multiplă a unei funcții de <math>\;n</math> variabile: <math>\;f(x_1,x_2,\ldots,x_n)</math> pe un domeniu <math>\;D</math> este reprezentată cel mai adesea prin semne succesive de integrare în ordinea inversă a execuției (cel mai din stânga semn de integrare este calculat ultimul) urmate de funcție și [[argumentul unei funcții|argumentele]] integrand în ordinea corectă (cel mai din dreapta argument este ultimul calculat). Domeniul de integrare este fie reprezentat simbolic pentru fiecare integrand în dreptul fiecărui semn de integrare, fie este abreviat ca variabilă la cel mai din dreapta semn de integrare:

:<math>
\iint \ldots \int_\mathbf{D}\;f(x_1,x_2,\ldots,x_n) \;\mathbf{d}x_1 \mathbf{d}x_2\!\ldots\mathbf{d}x_n
</math>

Fiind imposibil de calculat [[Primitivă|primitivele]] unei [[funcție de mai multe variabile reale|funcții de mai multe variabile]], nu există integrale multiple nedefinite. Toate integralele multiple sunt integrale definite.

== Exemple ==
De exemplu, volumul [[paralelipiped]]ului de laturi 4×6×5 mai poate fi obținut în mai multe moduri:
* Prin integrala dublă
:: <math>\iint_D 5 \ dx\, dy</math> 
: a funcției ''f''(''x'', ''y'') = 5 calculată pe regiunea ''D'' din planul ''xy'' care este baza paralelipipedului.
* Prin integrala triplă
::<math>\iiint_\mathrm{paralelipiped} 1 \, dx\, dy\, dz</math> 
:a funcției constante 1 calculată pe paralelipiped.

== Definiție matematică ==

Fie ''n'' un număr întreg mai mare ca 1. Se consideră un așa numit ''dreptunghi n-dimensional'' ''semiînchis''. Pentru un [[Plan (matematică)|plan]], ''n'' = 2, iar ''integrala multiplă'' este o integrală dublă.

: <math>T=[a_1, b_1)\times [a_2, b_2)\times\cdots \times [a_n, b_n)\subset \mathbb R^n.</math>

Se împarte fiecare interval <math>[a_i, b_i)</math> într-un număr finit de subintervale disjuncte, fiecare subinterval închis la stânga, și deschis la dreapta. Se notează aceste subintervale cu <math>I_i.</math> Atunci, familia de dreptunghiuri de forma

: <math>C=I_1\times I_2\times \cdots\times I_n</math>

este o [[Partiție (matematică)|partiție]] a lui <math>T,</math> adică subdreptunghiurile <math>C</math> sunt disjuncte și reuniunea lor este <math>T.</math> ''Diametrul'' unui dreptunghi <math>C</math> este, prin definiție, cea mai mare din lungimile subintervalelor al căror produs este <math>C,</math> iar diametrul unei partiții date a lui <math>T</math> este definit ca cel mai mare diametru al subdreptunghiurilor din partiție.

Fie <math>f:T\to \mathbb R</math> o funcție definită pe un dreptunghi <math>T.</math> Se consideră o partiție

:<math>T=C_1\cup C_2\cup \cdots \cup C_m</math>

a lui <math>T</math> definită ca mai sus, unde <math>m</math> este un întreg pozitiv. O '''sumă Riemann''' este o sumă de forma

:<math>\sum_{k=1}^m f(P_k)\, \operatorname{m}(C_k)</math>

unde pentru fiecare <math>k</math>, punctul <math>P_k</math> este din <math>C_k,</math> și <math>\operatorname{m}(C_k)</math> este produsul lungimilor intervalelor al căror [[produs cartezian]] este <math>C_k.</math>

Despre funcția <math>f</math> se spune că este '''integrabilă Riemann''' dacă [[Limită (matematică)|limita]]

:<math>S = \lim_{\delta \to 0} \sum_{k=1}^m f(P_k)\, \operatorname{m}\, (C_k)</math>

există, unde limita este calculată peste toate partițiile posibile ale lui <math>T</math> de diametru cel mult <math>\delta.</math> Dacă <math>f</math> este integrabilă Riemann, <math>S</math> se numește '''integrala Riemann''' a funcției <math>f</math> pe mulțimea <math>T,</math> și se notează cu

: <math>\int_T\!f(x)\,dx.</math>

[[Integrală Riemann|Integrala Riemann]] a unei funcții definită peste o mulțime <math>n</math>-dimensională cu limite arbitrare, poate fi definită prin extinderea acelei funcții la o funcție definită pe un dreptunghi semiînchis ale cărei valori sunt zero în afara funcției originale. Atunci, integrala funcției originale pe domeniul original este definită ca integrala funcției extinse pe domeniul dreptunghiular, dacă aceasta există.

În cele ce urmează, integrala Riemann în ''n'' dimensiuni va fi numită '''integrală multiplă'''.

=== Proprietăți ===
Integralele multiple au multe din proprietățile integralelor funcțiilor de o variabilă. În plus, ca și în cazul cu o singură variabilă, se poate folosi integrala multiplă ca media unei funcții pe o mulțime dată. Anume, dată fiind o mulțime <math>D\subset \mathbb R^n</math> și o funcție integrabilă <math>f</math> pe <math>D</math>, valoarea medie a lui <math>f</math> pe domeniul de definiție este dată de

:<math>\bar{f} = \frac{1}{m(D)} \int_D f(x)\, dx,</math>

unde <math>\;m(D)</math> este [[Măsură (matematică)|măsura]] lui <math>\;D</math>.

=== Cazuri particulare ===

În cazul <math>T \subseteq \mathbb{R}^2,</math> integrala

:<math>\ell = \iint_T f(x,y)\, dx\, dy </math>

este '''integrala dublă''' a lui ''F'' pe ''T'', și dacă <math>T \subseteq \mathbb{R}^3</math> integrala

:<math>\ell = \iiint_T f(x,y,z)\, dx\, dy\, dz</math>

este '''integrala triplă''' a lui ''F'' pe ''T''.

Prin convenție, integrala dublă are două semne de integrare, iar cea triplă are trei; aceasta este doar o convenție de notare, și este utilă la calculul unei integrale multiple ca integrală iterativă.

== Metode de integrare ==
Rezolvarea problemelor cu integrale multiple constă în majoritatea cazurilor în găsirea modalității de reducere a integralei multiple la o serie de integrale de o singură variabilă, fiecare fiind rezolvabilă prin metode specifice.

=== Examinare directă ===
Uneori, este posibil să se obțină rezultatul integrării fără calcule directe.

==== Constante ====
În cazul unei [[Funcție constantă|funcții constante]], rezultatul este direct: se înmulțește măsura domeniului cu funcția constantă ''c''. Dacă ''c'' = 1, și integrarea se face pe o subregiune a lui '''R'''<sup>2</sup> rezultă aria acelei regiuni, iar în '''R'''<sup>3</sup> este volumul regiunii.
* De exemplu:

::<math>D = \{ (x,y) \in \mathbb{R}^2 \ : \ 2 \le x \le 4 \ ; \ 3 \le y \le 6 \}</math> and <math>f(x,y) = 2\,\!</math>

:Se integrează ''f'' peste ''D'':

::<math>\int_2^4 \int_3^6 \ 2 \ dx\, dy = \mbox{area}(D) \cdot 2 = (2 \cdot 3) \cdot 2 = 12.</math>

==== Folosirea unor posibile simetrii ====
În cazul unui domeniu în care există simetrii față de una dintre axe și unde funcția are cel puțin o [[Paritatea funcțiilor|paritate]] în raport cu o variabilă, integrala devine nulă (suma valorilor opuse și egale în modul este zero).

Este suficient ca - în funcții definite pe '''R'''<sup>''n''</sup> - valoarea dependentă este impară în raport cu axa de simetrie.
* Exemplu (1):
:Se dă <math>f(x,y) = 2 \ \sin(x) - 3 \ y^3 + 5</math> și <math>T = x^2 + y^2 \le 1</math> suprafața de integrare (discul de [[rază]] 1 cu centrul în origina axelor, cu marginea inclusă).

:Folosind proprietate de liniaritate, integrala poate fi descompusă în trei părți:

::<math>\iint_T (2 \ \sin(x) - 3 \ y^3 + 5) \ dx \, dy = \iint_T 2 \ \sin(x) \ dx \, dy - \iint_T 3 \ y^3 \ dx \ dy + \iint_T 5 \ dx \ dy</math>

:2 sin(''x'') și 3''y''<sup>3</sup> sunt ambele funcții impare și în plus este evident că discul ''T'' este simetric și în raport cu axa ''x'' și în raport cu axa ''y''; astfel, singura contribuție la rezultatul final al integralei o are funcția constantă ''5'' deoarece celelalte componente sunt nule.
* Exemplu (2):
:Se consideră funcția <math>f(x,y,z) = x \ e^{y^2 + z^2}</math> și - ca regiune de integrare - [[sfera]] de rază 2 centrată în originea axelor <math>T = x^2 + y^2 + z^2 \le 4</math>. "Bila" este simetrică în raport cu axele, dar este suficient să se integreze în raport cu axa ''x'' pentru a arăta că integrala este 0, deoarece funcția este impară în acea variabilă.

=== Formule de reducere ===
Formulele de reducere utilizează conceptul de [[domeniu simplu]] pentru a face posibilă descompunerea integralei multiple ca produs de alte integrale de o singură variabilă. Acestea trebuie să fie rezolvate de la dreapta la stânga, ținând cont că celelalte variabile sunt constante (aceeași procedură ca și la calculul [[Derivată parțială|derivatelor parțiale]]).

==== Domenii normale din '''R'''<sup>2</sup> ====
===== Axa ''x'' =====
Dacă ''D'' este un domeniu măsurabil perpendicular pe axa ''x'' și <math>f: D \longrightarrow \mathbb{R}</math> este o [[funcție continuă]]; atunci α(''x'') și β(''x'') (definite pe itnervalul [''a'',''b'']) sunt două funcții care determină ''D''. Atunci:

:<math>\iint_T f(x,y)\ dx\, dy = \int_a^b dx \int_{ \alpha (x)}^{ \beta (x)} f(x,y)\, dy.</math>

===== Axa ''y'' =====
Dacă ''D'' este un domeniu măsurabil perpendicular pe axa ''y'' și <math>f: D \longrightarrow \mathbb{R}</math> este o [[funcție continuă]]; atunci ''α(y)'' și β(''y'') (definite pe intervalul [''a'',''b'']) sunt două funcții care determină ''D''. Atunci:

:<math>\iint_T f(x,y)\ dx\, dy = \int_a^b dy \int_{ \alpha (y)}^{ \beta (y)} f(x,y)\, dx.</math>

[[Fișier:Esempio-formulediriduzione-r2.svg|thumb|160px|right|Exemplu: Regiunea D pentru integrale calculate prin metoda reducerii.]]

===== Exemplu =====
:Fie regiunea <math>D = \{ (x,y) \ : \ x=0, y=1, y=x^2 \}</math> (vezi graficul). Se calculează

::<math>\iint_D (x+y) \, dx \, dy.</math>

:Acest domeniu este perpendicular pe axele x și y. Pentru a aplica formulele, trebuie găsite funcțiile care determină ''D''.
:În acest caz, cele două funcții sunt:

::<math>\alpha (x) = x^2\,\!</math> și <math>\beta (x) = 1\,\!</math>

:intervalul este dat de intersecțiile funcțiilor cu dreapta <math>x = 0\,\!</math>, deci intervalul este <math>[a,b] = [0,1]\,\!</math> (normalitatea a fost aleasă în raport cu axa x).

:Acum se pot aplica formulele:

::<math>\iint_D (x+y) \, dx \, dy = \int_0^1 dx \int_{x^2}^1 (x+y) \, dy = \int_0^1 dx \ \left[xy \ + \ \frac{y^2}{2} \ \right]^1_{x^2}</math>

:(la început, se calculează a doua integrală, considerând ''x'' constant). Restul operațiunilor constau din aplicarea tehnicilor de bază de integrare:

::<math>\int_0^1 \left[xy \ + \ \frac{y^2}{2} \ \right]^1_{x^2} \, dx = \int_0^1 \left(x + \frac{1}{2} - x^3 - \frac{x^4}{2} \right) dx = \cdots = \frac{13}{20}.</math>

:Dacă se alege normalitatea în raport cu axa y atunci se calculează

::<math>\int_0^1 dy \int_0^{\sqrt{y}} (x+y) \, dx.</math>

:și se obține aceeași valoare.
[[Fișier:Dominio-normalità r3 esempio.svg|thumb|160px|right|Exemplu de domeniu normal în R<sup>3</sup> (planul ''xy'').]]

==== Domenii normale în R<sup>3</sup> ====
Extensia acestor formule la integralele triple este evidentă:

''T'' este un domeniu perpendicular pe planul ''xy'' în raport cu funcțiile α (''x'',''y'',''z'') și β(''x'',''y'',''z''). Atunci:

:<math>\iiint_T f(x,y,z) \ dx\, dy\, dz = \iint_D dx\, dy \int_{\alpha (x,y,z)}^{\beta (x,y,z)} f(x,y,z) \, dz</math>

(aceeași definiție există și pentru celelalte cinci domenii de normalitate din '''R'''<sup>3</sup>).

=== Schimbarea de variabilă ===

Limitele de integrare nu sunt de multe ori ușor de interschimbat (în absența normalității sau a unor formule complexe de integrare), și în acest caz se efectuează o ''schimbare de variabile'' pentru a rescrie integrala într-o regiune mai "comodă", descrisă de formule mai simple. Pentru a face aceasta, funcția trebuie să fie adaptată noilor coordonate.

:<u>Exemplu (1-a)</u>: 
::Funcția este <math>f(x, y) = (x-1)^2 +\sqrt y</math>;
::dacă se adoptă această substituție: <math>x' = x-1, \ y'= y \, \!</math> atunci <math>x = x' + 1, \ y=y' \,\!</math>
::se obține noua funcție <math>f_2(x,y) = (x')^2 +\sqrt y</math>. 
* Analog pentru domeniu, deoarece el este delimitat de variabilele originale care au fost transformate anterior (''x'' și ''y'' în exemplu).
* diferențialele ''dx'' și ''dy'' se transformă prin [[Jacobian|determinantul matricii jacobiene]] ce conține derivatele parțiale ale transformărilor privind noua variabilă (ca, de exemplu, transformarea în [[coordonate polare]]).

Există trei tipuri principale de schimbări de variabile (unul în '''R'''<sup>2</sup>, două în '''R'''<sup>3</sup>); totuși, se poate găsi o schimbare de variabilă potrivită pe același principiu la un mod mai general.

==== Coordonate polare ====
[[Fișier:Passaggio in coordinate polari.svg|thumb|270px|right|Transformarea din coordonate polare în coordonate carteziene.]]
În '''R'''<sup>2</sup>, dacă domeniul are [[simetrie circulară]] și funcția are anumite caracteristici deosebite, se poate aplica ''transformarea în coordonate polare'' ceea ce înseamnă că punctele generice ''P(x,y)'' în coordonate carteziene se transformă în puncte corespunzătoare lor în coordonate polare. Aceasta permite modificarea [[formă |formei]] domeniului și simplificarea calculelor.

Relația fundamentală pe care se bazează transformarea este următoarea:
:<math>f(x,y) \rightarrow f(\rho \ \cos \phi,\rho \ \sin \phi ).</math>

<u>Exemplu (2-a)</u>:<br />
:Funcția este <math>f(x,y) = x + y\,\!</math>

:șî la aplicarea transformării se obține
::<math>f(\rho, \phi) = \rho \cos \phi + \rho \sin \phi = \rho \ (\cos \phi + \sin \phi ).</math>

<u>Exemplu (2-b)</u>:
:Funcția este <math>f(x,y) = x^2 + y^2\,\!</math>

:În acest caz avem:
::<math>f(\rho, \phi) = \rho^2 (\cos^2 \phi + \sin^2 \phi) = \rho^2\,\!</math>

:folosind identitatea triunghiului (foarte utilă în simplificarea operației).

Transformarea domeniului domeniului se face definind circumgerința și amplitudinea unghiului descris pentru a defini intervalele ρ, φ pornind de la ''x'', ''y''.
[[Fișier:Esempio trasformazione dominio da cartesiano polare.svg|thumb|230px|right|Exemplu de transformare de domeniu de la sistemul cartezian la cel polar.]]

<u>Exemplu (2-c)</u>:
:Domeniul este <math>D = x^2 + y^2 \le 4\,\!</math>, adică un disc de rază 2; este evident că unghiul acoperit este întreg cercul, deci φ variază de la 0 la 2π, iar raza variază de la 0 la 2.

<u>Exemplu (2-d)</u>:
:Domeniul este <math>D = \{ x^2 + y^2 \le 9, \ x^2 + y^2 \ge 4, \ y \ge 0 \}</math>, adică coroana circulară din semiplanul din semiplanul pozitiv ''y''; se observă că φ descrie un unghi plan iar ρ merge de la 2 la 3. Deci domeniul transformat va fi următorul [[dreptunghi]]:
::<math>T = \{ 2 \le \rho \le 3, \ 0 \le \phi \le \pi \}</math>.

[[Jacobian|Determinantul jacobian]] al transformării este următorul:
:<math>\frac{\partial (x,y)}{\partial (\rho, \phi)} = 
\begin{vmatrix}
\cos \phi & - \rho \sin \phi \\
\sin \phi & \rho \cos \phi 
\end{vmatrix} = \rho
</math>

care a fost obținut introducând derivatele parțiale ale funcțiilor ''x' = ρ cos(φ), ''y'' = ρ sin(φ) în prima coloană în raport cu ρ șî în a doua în raport cu φ, deci diferențialele ''dx dy'' din tranfosrmare devin ρ ''d''ρ ''d''φ.

Odată ce funcția este transformată și domeniul este evaluat, se poate defini formula pentru schimbarea de variabile în [[coordonate polare]]:

:<math>\iint_D f(x,y) \ dx\, dy = \iint_T f(\rho \cos \phi, \rho \sin \phi) \rho \, d \rho\, d \phi.</math>

Se observă că φ este valid în intervalul [0, 2π] în timp ce ρ, deoarece este o măsură a lungimii, poate avea doar valori pozitive.

<u>Exemplu (2-e)</u>:
:Funcția este <math>f(x,y) = x\,\!</math> este același ca în exemplul 2-d.
:Din analiza anterioară a domeniului ''D'' știm că intervalele ρ (de la 2 la 3) și φ (de la 0 la π). Funcția se transformă:
::<math>f(x,y) = x \longrightarrow f(\rho,\phi) = \rho \ \cos \phi.</math>

:în cele din urmă se aplică formula de integrare:
::<math>\iint_D x \, dx\, dy = \iint_T \rho \cos \phi \ \rho \, d\rho\, d\phi.</math>

:Odată ce sunt cunoscute intervalele, se poate calcula
::<math>\int_0^{\pi} \int_2^3 \rho^2 \cos \phi \ d \rho \ d \phi = \int_0^{\pi} \cos \phi \ d \phi \left[ \frac{\rho^3}{3} \right]_2^3 = \left[ \sin \phi \right]_0^{\pi} \ \left(9 - \frac{8}{3} \right) = 0.</math>

[[Categorie:Calcul integral]]