%\input{/home/moi/Lycee/makalice}\setlength{\parindent}{0pt}
\documentclass[a4paper,10pt,latin1]{polymaths}
%\documentclass[xcolor=pst,dvips]{beamer}
%\usepackage{BeamerPolymaths}
\entete{\textbf{\NomChap\ } - \NomSection }{ \NoLecon\   Aventure\  - \NomLecon}{IUT de Nantes - SGM1}
\renewcommand\thesubsection{}
\renewcommand{\e}{{\rm e}}
%\setcounter{page}{66}

\begin{document}

%\tableofcontents

\setlength{\parindent}{0pt}
\LECON{VII}{Septième}{INT\'EGRALES DOUBLES }{  }{Intégrales doubles}{}






\begin{Abstract}
Nous allons découvrir dans ce chapitre une notion aux innombrables applications physiques. Compte tenu du peu d'heures dont nous disposons, nous nous contenterons d'un survol  rapide en sacrifiant la cohérence et la rigueur du raisonnement mathématique.
\end{Abstract}



\SECTION{Définition}{Définition}

Nous avions introduit rigoureusement le calcul intégral lors du chapitre V. Grosso modo, nous avions interprété l'intégrale en termes de limite d'une somme et par suite en termes d'aire. Nous allons appliquer le même  principe aux fonctions de deux variables.

Nous travaillerons sur des domaines du plan limités par une courbe fermée, qui sont  en quelque sorte pour $\bbr^2$   ce que sont les intervalles fermés bornés (les segments ) pour $\bbr$. 

Considérons donc un domaine $\DR$ dont les abscisses extrêmes sont respectivement $a$ et $b$ ($a<b)$ d'une part et $c$ et $d$ ($c<d$) d'autre part, et une fonction $f$ de deux variables \textbf{définie sur } $\Mathbold{\DR}$ \textbf{ et nulle ailleurs}.

Comme nous l'avions fait pour les fonctions d'une variable, nous allons subdiviser les segments $[a\pv b]$ et $[c\pv d]$ 
\[ a=x_1<x_2<\cdots<x_n=b\] \[ c=y_1<y_2<\cdots<y_p=d\]

On découpe ainsi notre domaine $\DR$ en petits rectangles $r_{ij}$ définis par 
\[ x_i<x<x_{i+1}\ \ \ y_j<y<y_{j+1}\]
et d'aire $\Delta x_i\times \Delta y_j$ avec  \[ \Delta x_i=x_{i+1}-x_i\ \ \ \ \Delta y_j=y_{j+1}-y_j\] 

On va prendre un point $M_{ij}$ de coordonnées $(\al_i,\beta_j)$  appartenant à chacun de ces rectangles et qui les représentera. 

\ENC{\Def}{On dira alors que la fonction de deux variables $f$ est intégrable si et seulement si le nombre $\Sum_{i=1}^{n}\Sum_{j=1}^pf(\al_i,\beta_j)\times\Delta x_i\times \Delta y_j$, c'est à dire la somme de tous les produits $f(\al_i,\beta_j)\times\Delta x_i\times \Delta y_j$ sur tous les rectangles $r_{ij}$ tend vers un nombre fini constant $\IR(f,\DR)$ quand les $\Delta x_i$ et les $\Delta y_j$ tendent vers 0.\\
\sm
Dans ce cas, on note \[ \ds\IR(f,\DR)=\intt_{\DR}f(x,y)\d x\d y \]
}


On admettra le résultat fondamental suivant 

\ENC{\Thm}{Si $f$ est continue sur $\DR$, alors elle est intégrable sur $\DR$}

...sauf qu'on ne sait pas trop ce qu'est une fonction de 2 variables continue. Donc on supposera que tous les exercices proposées et toutes les situations physiques rencontrées utiliseront des fonctions intégrables là où il faut : fermez les yeux et faites-nous confiance...mouais. 
 

\SECTION{Interprétations graphiques}{Interprétations graphiques}


\SUBSECTION{Aire du domaine $\Mathbold{\DR}$}

D'après notre approche de la définition d'une fonction intégrable, $\intt_{\DR}\d x\d y$ est la limite de la somme des aires des <<~petits rectangles~>>, donc c'est l'aire du domaine $\DR$.

\SUBSECTION{Volume d'un cylindre} 
Chaque produit $f(\al_i,\beta_j)\times\Delta x_i\times \Delta y_j$ est égal au volume du pavé droit de base le rectangle $r_{ij}$ et de hauteur $f(\al_i,\beta_j)$\footnote{Voir la figure faite au tableau pendant votre partie de bataille navale}. Or dans un RON $\oijk$, $z=f(x,y)$ est l'équation d'une surface $\SR$. Alors l'intégrale double $\intt_{\DR}f(x,y)\d x\d y $ est en fait la mesure du volume limité par la surface $\SR$, le plan $(xOy)$ et le cylindre de section droite le domaine $\DR$ : c'est une généralisation de l'interprétation de l'intégrale (simple) mais avec une dimension de plus.

\begin{description}
\item[IMPORTANT] comme la double intégration est une généralisation de la simple, on admettra que les propriétés vues à la leçon V s'appliquent aux intégrales doubles, à savoir la linéarité, la positivité, la conservation de l'ordre, etc.
\end{description}

\Exot{}
\'Enocez ces propriétés dans le cas de l'intégrale double en les transcrivant sous forme de ... formules.

\SECTION{Calcul des intégrales doubles en coordonnées cartésiennes}{Calcul des intégrales doubles en coordonnées cartésiennes}

\SUBSECTION{Le domaine $\Mathbold{\DR}$ est rectangulaire}

C'est le cas le plus simple, car le plus proche de l'intégration simple\footnote{Sachant toujours que nous ne considérons que des fonctions continues, et que toutes les fonctions rencontrées seront continues là où il faut, donc nous ne nous poserons plus de question à ce sujet pourtant fondamental.}.  

Soit donc $\DR$ le domaine $[a\pv b]\times [c\pv d]$. 

L'intégrale $\Int{c}{d}f(x,y)\d y$ est une fonction ne dépendant que de $x$ : notons $F(x)=\ds\int_c^df(x,y)\d y$. 
 
On peut alors calculer l'intégrale $I_1=\ds\int_a^bF(x)\d x=\int_a^b\pa{\int_c^df(x,y)\d y}\d x.$

On calcule de même l'intégrale $I_2=\ds\int_c^dG(y)\d y=\int_c^d\pa{\int_a^bf(x,y)\d x}\d y.$

Des théorèmes qui sortent du cadre de notre cours permettent de démontrer qu'en fait $I_1=I_2=\ds\intt_{\DR}f(x,y)\d x\d y$ :

\ENC{\Thm\ \ Théorème de Fubini dans le cas d'un domaine rectangulaire}{
\[ \ds\intt_{\DR}f(x,y)\d x\d y=\int_a^b\pa{\int_c^df(x,y)\d y}\d x=\int_c^d\pa{\int_a^bf(x,y)\d x}\d y \]}

\begin{description}
\item[Exemple :] Soit $\DR$ le domaine défini par $0\ie x\ie 2$ et $0\ie y\ie 1$ et $f~:~(x,y)\mapsto f(x,y)=x\e^{xy}+3xy^2$.

\begin{center}
\begin{arrayl}
\ds\intt_{\DR}f(x,y)\d x\d y&=&\ds\int_0^2\pa{\ds\int_0^1x\e^{xy}+3xy^2\d y}\d x\\
                                         &=&\ds\int_0^2\pa{\cro{\fr{x}{x}\e^{xy}+\fr{3xy^3}{3}}_0^1}\d x\\
                                        &=&\ds\int_0^2\pa{\e^x-1+x}\d x\\
                                        &=&\cro{\e^x-x+\fr{x^2}{2}}_0^2\\
                                        &=&\e^2-2+2-1=\e^2-1
\end{arrayl}
\end{center}
\end{description}

\Exot{}
Reprendre l'exemple précédent, mais en commençant par intégrer par rapport à $x$ (Reconnaître une forme $u/v$ : difficile...).

\Exot{}
On considère le   domaine $\DR=\big\{(x,y)\in\bbr^2\,|\, 1\ie x\ie 2,\ 1\ie y\ie 2\big\}$. Calculer $\ds\intt_{\DR}f(x,y)\d x\d y$ dans les cas suivants
\begin{nliste}{2}
\item $f(x,y)=\fr{x^2}{1+y^2}$
\item $f(x,y)=\fr{1}{\pa{x+y}^2}$
\end{nliste}


\SUBSECTION{Domaine quelconque}

Considérons une courbe fermée un peu spéciale : des parallèles aux axes de coordonnées ne la coupent qu'en deux points au maximum\footnote{Réveillez-vous ! Le dessin est au tableau...} (on dit qu'elle est convexe). Elle est comme d'habitude comprise dans le rectangle $[a\pv b]\times [c\pv d]$.

Alors pour une valeur de $x\in[a\pv b]$ donnée, on peut définir les deux ordonnées des points d'intersection de la courbe avec la droite <<~verticale~>> passant par le point de coordonnées $(x,0)$. Notons ces coordonnées $y_1(x)$ et $y_2(x)$.

On peut faire de même avec  toute valeur de $y$ dans $[c\pv d]$.

On montre alors que le théorème de Fubini s'applique encore

\ENC{\Thm\ \ Théorème de Fubini dans le cas d'un domaine <<~convexe~>>}{
\[ \ds\intt_{\DR}f(x,y)\d x\d y=\int_a^b\pa{\int_{y_1(x)}^{y_2(x)}f(x,y)\d y}\d x=\int_c^d\pa{\int_{x_1(y)}^{x_2(y)}f(x,y)\d x}\d y \]}

Lorsque $\DR$ n'est plus convexe, on découpe si possible  le domaine en parties convexes disjointes.


\begin{description}
\item[Exemple : ]  Calculons $I=\ds\intt_{\DR}xy\d x\d y$ avec $\DR=\ac{(x,y)\in\bbr^2\ |\ x\se 0,\ 0\ie y\ie 1,\ x+y\ie 3}$.

Il faut commencer par dessiner le domaine. Il s'agit d'un trapèze rectangle.  Le plus simple est de <<~balayer~>> le domaine à $y$ constant.
Alors \[I=\ds\int_0^1y\pa{\int_0^{3-y}x\d x}\d y=\fr{1}{2}\int_0^1y(3-y)^2\d y=\fr{1}{2}\cro{\fr{y^4}{4}-2y^3+\fr{9y^2}{2}}_0^1=\fr{11}{8}\]

Essayons maintenant de balayer à $x$ constant. Nous allons découper le trapèze  en deux parties par une droite parallèle à $(Ox)$ :  cela donne un rectangle $\RR$ et un triangle rectangle $\TR$ (voir la figure au tableau). 

D'une part $\ds\intt_{\RR}xy\d x\d y=\int_0^2 x\pa{\int_0^1 y\d y}\d x=1$

D'autre part $\ds\intt_{\TR}xy\d x\d y=\int_2^3 x\pa{\int_0^{3-x} y\d y}\d x=\fr{1}{2}\int_2^3x(3-x)^2\d x=\ud \cro{\fr{x^4}{4}-2x^3+\fr{9x^2}{2}}_2^3=\fr{3}{8}$

Finalement, on obtient $I=1+\ofr{3}{8}=\ofr{11}{8}$.




\end{description}



\Exot{}
Calculer $\ds\intt_{\DR}f(x,y)\d x\d y$ dans les cas suivants
\been
\item $f(x,y)=x^2y$ et $\DR=\ac{(x,y)\in\bbr^2\,|\,0\ie x\ie1, \, y\se 0,\, x-y\se 0}$
\item $f(x,y)=(x^2+y)$ et $\DR$ est le triangle de sommets $O(0,0)$, $A(2,1)$ et $B(2,-2)$.
\item $f(x,y)=x\e^{x+2xy}$, $\DR=\ac{(x,y)\in\bbr^2\,|\,0\ie x\ie 2\, ,0\ie y\ie 2\, , xy\se2}$
\een  


\SECTION{Passage en coordonnées polaires}{Passage en coordonnées polaires}

Lorsque le domaine étudié est une portion de disque ou d'une surface engendrée par une courbe en coordonnées polaires, il est utile d'effectuer un changement de variables. Le domaine $\DR$ en fonction de $x$ et $y$ devient un domaine $\Delta$ en fonction de $r$ et $\theta$ qu'on espère plus simple à <<~man\oe uvrer~>>. 

Nous admettrons la formule suivante qui utilise un clacul de déterminant que vous découvrirez en algèbre linéaire.

\ENC{\Thm\ \ Passage en coordonnées polaires}{\[ \ds\intt_{\DR}f(x,y)\d x\d y=\int_{\Delta}f(r\cos \theta, r\sin\theta) r\, \d r\d \theta\]} 




\begin{description}
\item[Exemple : ] Calculons $I=\ds\intt_{\DR} (x^2+y^2)\d x\d y$, avec $\DR$ le quart de cercle de centre O(0,0), de rayon 2 avec $x\se 0$ et $y\se 0$.

D'abord, un petit dessin. 
 
Sachant que $x^2+y^2=r^2$, $I=\ds\intt_{\Delta}r^2\times r\,\d r\d \theta$, avec $\Delta=\ac{(r,\theta)\,|\, 0\ie r\ie2,\, 0\ie\theta\ie \ofr{\pi}{2}}$

\[I=\ds\int_0^{\efr{\pi}{2}}\d \theta\int_0^2r^3\d r=\fr{\pi}{2}\times 4=2\pi\]

\end{description}


\Exot{}
\been
\item Calculer $\ds\intt_{\DR}\fr{1}{1+x^2+y^2}\d x\d y$ avec $\DR=\ac{(x,y)\,|\, x^2+y^2\ie 1}$
\item Calculer $\ds\intt_{\DR}x^2+y^2\d x\d y$ avec $\DR=\ac{(x,y)\,|\, x^2+y^2\ie x}$
\item Calculer $\ds\intt_{\DR}(x^2+y^2)\d x\d y$ avec $\DR=\ac{(x,y)\,|\, x^2+y^2\ie x,\, x^2+y^2\se y}$
\een

\Exot{}

Vous allez bientôt utiliser en probabilité l'intégrale  $I=\ds\int_{-\infty}^{+\infty}\e^{-t^2}\d t$. Nous n'avons pas étudié ce genre d'intégrale avec une borne infinie, mais étant donné que le résultat que nous allons établir va vous être très rapidement utile et qu'une méthode de l'obtenir utilise les intégrales doubles, nous allons malgré tout nous lancer.

Nous admettrons que I existe bel et bien et que de plus $I=\Lim[a]{+\infty}I(a)$ avec $I(a)=\ds\int_{-a}^a\e^{-t^2}\d t$
\been
\item Montrez que $J(a)=\ds\intt_{\DR(a)}\e^{-(x^2+y^2)}\d x\d y=I^2(a)$ avec $\DR(a)$ le carré $[-a,a]\times[-a,a]$.
\item Nous ne connaissons pas de primitive de $t\mapsto  \e^{-t^2}$, alors le <<~$x^2+y^2$~>> nous incite à vouloir passer en coordonnées polaires. Malheureusement, l'expression polaire du  domaine $\DR(a)$ semble problématique. Alors nous allons ruser. 
Appelons $\CR(a)$ le disque   de centre O et de rayon $a$ et $G(a)=\ds\intt_{\CR(a)}\e^{-(x^2+y^2)}\d x\d y$. Montrez que $G(a)\ie J(a)\ie G(a\sqrt{2})$.  
\item Faites un dessin, et calculez $G(a)$ à l'aide d'un passage en coordonnées polaires.
\item Déduisez-en $\Lim[a]{+\infty}G(a)$, puis $\Lim[a]{+\infty}J(a)$. Déduisez-en I.
\een

\SECTION{Applications physiques}{Applications physiques}

\SUBSECTION{Masse d'une <<~plaque~>>}

Si on connait la masse surfacique $f(x,y)$ d'une <<~plaque~>> homogène $\DR$, alors sa masse vaut   $\ds\intt_{\DR}f(x,y)\d x\d y$

\SUBSECTION{Moments d'inertie}

Avec les mêmes hypothèses que précédemment, les moments d'inertie de la plaque par rapport à $(Ox)$ et $(Oy)$ valent respectivement $\ds\intt_{\DR}y^2\,f(x,y)\d x\d y$ et $\ds\intt_{\DR}x^2\,f(x,y)\d x\d y$

\SUBSECTION{Centre de gravité}

Toujours dans les mêmes conditions, les coordonnées du centre de gravité de la plaque sont $\fr{1}{\MR}\ds\intt_{\DR}x\,f(x,y)\d x\d y$ et $\fr{1}{\MR}\ds\intt_{\DR}y\,f(x,y)\d x\d y$, avec $\MR$ la masse de la plaque dont le calcul a été explicité plus haut.


\Exot{}
Déterminez le centre de gravité du domaine limité par le cercle de centre O, de rayon 1 et par les demi-droites d'équations $\theta=\ofr{\pi}{6}$ et $\theta=\ofr{\pi}{2}$.

\Exot{} 
Déterminez le moment d'inertie d'un disque homogène de centre O et de rayon 1 inclus dans le plan ($xOy)$ par rapport à l'axe $(Oy)$


\end{document} 