Intégrale paramétrique

En mathématiques, et plus précisément en analyse, une intégrale paramétrique (également appelée intégrale à paramètre) est une fonction d'une variable, définie à partir d'une fonction de deux variables – la variable d'intégration et le paramètre – par intégration sur un ensemble fixe par rapport à la variable d'intégration. Un exemple typique est la fonction gamma d'Euler définie par $\Gamma (x)=\int _{0}^{+\infty }t^{x-1}\exp(-t)~\mathrm {d} t$, pour laquelle on intègre sur $[0,+\infty [$ par rapport à la variable d'intégration $t$, la variable $x$ étant ici le paramètre.

Les deux variables, ainsi que les valeurs de la fonction, sont souvent choisies dans un espace euclidien. Une classe importante d'exemples est l'ensemble des transformées, dont la transformée de Fourier.

Définition formelle

Soient T un ensemble, $(\Omega ,{\mathcal {A}},\mu )$ un espace mesuré et

$f:T\times \Omega \to \mathbb {R} ^{n}$

une application telle que pour tout élément t de T, l'application

$\Omega \to \mathbb {R} ^{n},\quad \omega \mapsto f(t,\omega )$ soit intégrable.

Alors l'application F définie par :

$F:T\to \mathbb {R} ^{n},\quad F(t)=\int f(t,\omega )~\mathrm {d} \mu (\omega )$

est appelée une intégrale paramétrique.

Le plus souvent, dans les applications :

• l'entier naturel n est égal à 1 ;
• T est un ouvert de ℝ ;
• $(\Omega ,{\mathcal {A}},\mu )$ est une partie d'un espace euclidien, implicitement munie des tribu et mesure de Lebesgue ou de Borel.
• les fonctions sont continues et les intégrales sont considérées au sens de Riemann, mais la théorie générale de Lebesgue s'applique à ce cas particulier : sur un segment, une fonction bornée est Riemann-intégrable si et seulement si elle est continue presque partout, et toute fonction Riemann-intégrable est Lebesgue-intégrable.

Exemples

Transformée de Fourier

Soit g une fonction intégrable de ℝⁿ dans ℂ, la transformée de Fourier de g est la fonction de ℝⁿ dans ℂ définie par :

${\widehat {g}}(t)={\frac {1}{\sqrt {2\pi }}}\int _{\mathbb {R} ^{n}}\exp[-2i\pi (t|\omega )]g(\omega )\,\mathrm {d} \omega ,$

où $(\cdot |\cdot )$ désigne le produit scalaire usuel.

Fonction gamma d'Euler

Le graphe de la fonction gamma. Sur l'axe des abscisses, on fait varier le paramètre $x$.

La fonction gamma d'Euler est définie entre autres pour tout réel x strictement positif, par :

$\Gamma (x)=\int _{0}^{+\infty }t^{x-1}\exp(-t)~\mathrm {d} t.$

Potentiel du champ de gravitation

Le potentiel du champ de gravitation V(x) créé par un corps matériel M de densité variable ρ en un point x de ℝ³ extérieur à M est donné par :

$V(x)=-G\int _{M}{\rho (y) \over \|x-y\|_{2}}~\mathrm {d} y,$

où G désigne la constante de gravitation et $\|\cdot \|_{2}$ la norme euclidienne.

Limite

Reprenons la définition formelle ci-dessus en supposant de plus que T est une partie de ℝ, que x est un réel adhérent à T, et que :

• $\forall \omega \in \Omega \quad \lim _{t\to x}f(t,\omega )=\varphi (\omega )$ pour une fonction $\varphi$ ;
• il existe une application intégrable $g:\Omega \to \mathbb {R}$ telle que

$\forall t\in T\quad \forall \omega \in \Omega \quad \|f(t,\omega )\|\leq g(\omega )$.

Alors, le théorème de convergence dominée permet de prouver que φ est intégrable et que

$\lim _{t\to x}F(t)=\int \varphi (\omega )~\mathrm {d} \mu (\omega ),$

soit encore :

$\lim _{t\to x}\left(\int _{\Omega }f(t,\omega )~\mathrm {d} \mu (\omega )\right)=\int \left(\lim _{t\to x}f(t,\omega )\right)~\mathrm {d} \mu (\omega ).$

Remarques.

• La première hypothèse peut être affaiblie en supposant que la limite existe seulement pour presque tout ω ∈ Ω, sous réserve que l'espace mesuré $(\Omega ,{\mathcal {A}},\mu )$ soit complet (ce qui est le cas pour les tribu et mesure de Lebesgue).
• La seconde hypothèse peut être doublement affaiblie en supposant seulement qu'il existe une fonction intégrable g telle que pour chaque élément t de T appartenant à un certain voisinage de x on ait : $\|f(t,\cdot )\|\leq g$ presque partout.
• Les énoncés des sections suivantes possèdent des variantes analogues.
• L'énoncé ci-dessus, même ainsi renforcé, reste vrai quand T et x sont une partie et un élément d'un espace métrique autre que ℝ (par exemple ℝ ou ℝ²).

Démonstration

Soit $(t_{n})_{n\in \mathbb {N} }$ une suite dans T qui converge vers x. La suite $(f(t_{n},\cdot ))_{n\in \mathbb {N} }$ de fonctions intégrables converge simplement vers φ et l'on a, par la seconde hypothèse :

$\forall n\in \mathbb {N} \quad \forall \omega \in \Omega \quad \|f(t_{n},\omega )\|\leq g(\omega )$.

Le théorème de convergence dominée entraîne alors l'intégrabilité de φ et les relations :

$\lim _{n\to \infty }F(t_{n})=\lim _{n\to \infty }\int f(t_{n},\omega )\mathrm {d} \mu (\omega )=\int \varphi (\omega )~\mathrm {d} \mu (\omega )$.

Continuité

Continuité locale : si l'on reprend la section précédente en supposant de plus que x appartient à T (donc pour tout ω ∈ Ω, $f(\cdot ,\omega )$ est continue au point x et $\varphi (\omega )=f(x,\omega )$), on en déduit que F est continue en x.

Continuité globale : par conséquent, si f est continue sur T × Ω avec T partie ouverte (ou plus généralement : localement compacte) de ℝ et Ω fermé borné d'un espace euclidien, alors F est définie et continue sur T.

Démonstration

Pour tout élément t de T, $f(t,\cdot )$ est continue sur le compact Ω, donc intégrable sur Ω pour la mesure de Lebesgue, si bien que F est définie sur T.

Soit x ∈ T. Pour tout ω ∈ Ω, $f(\cdot ,\omega )$ est continue sur T. De plus, si K est un voisinage compact de x dans T alors, par continuité de f, il existe une constante M telle que :

$\forall t\in K\quad \forall \omega \in \Omega \quad \|f(t,\omega )\|\leq M.$

En prenant g = M dans la proposition précédente, cela prouve que F est continue en x.

Dérivabilité

La règle de dérivation sous le signe d'intégration est connue sous le nom de règle de Leibniz (pour d'autres règles portant ce nom, voir Règle de Leibniz ).

Étude locale

Reprenons la définition formelle ci-dessus en supposant de plus que T est un intervalle de ℝ et que :

• pour tout ω ∈ Ω, $f(\cdot ,\omega )$ est dérivable sur T ;
• il existe une application intégrable g : Ω → ℝ telle que

$\forall t\in T,\quad \left\|{\partial f \over \partial t}(t,\cdot )\right\|\leq g$.

Alors, pour tout $t\in T$ l'intégrale paramétrique F est dérivable au point t, l'application ${\partial f \over \partial t}(t,\cdot )$ est intégrable, et :

$F'(t)=\int {\partial f \over \partial t}(t,\omega )~\mathrm {d} \mu (\omega ).$

Démonstration

Fixons t ∈ T et posons, pour tout ω ∈ Ω et tout réel h non nul tel que t+h ∈ T :

$k(h,\omega )={\frac {f(t+h,\omega )-f(t,\omega )}{h}}.$

On a alors :

• $\lim _{h\to 0}k(h,\omega )={\partial f \over \partial t}(t,\omega )$ ;
• $\|k(h,\omega )\|\leq g(\omega )$ (d'après l'inégalité des accroissements finis).

L'énoncé de la section « Limite » permet de conclure.

Étude globale

Avec les mêmes hypothèses que dans l'énoncé « Continuité globale » (f est continue sur T × Ω avec T partie localement compacte de ℝ et $\Omega$ fermé borné d'un espace euclidien), si l'on suppose de plus que ${\partial f \over \partial t}$ est définie et continue sur T × Ω, alors F est de classe C¹ sur T et pour tout x ∈ T, on a :

$F'(x)=\int {\partial f \over \partial t}(x,\omega )~\mathrm {d} \mu (\omega ).$

Démonstration

Soit K un compact de T. Par continuité de ${\partial f \over \partial t}$ sur le compact T × Ω, il existe une constante M telle que :

$\forall x\in K\quad \forall \omega \in \Omega \quad \left\|{\partial f \over \partial t}(x,\omega )\right\|\leq M.$

En prenant g = M dans la proposition précédente, cela prouve que F est dérivable (avec la formule annoncée) sur tout compact K de T, donc sur T.

La continuité de F' résulte alors de l'énoncé « Continuité globale ».

Forme générale unidimensionnelle

Le résultat suivant peut être vu comme une généralisation du premier théorème fondamental de l'analyse et peut s'avérer utile dans le calcul de certaines intégrales réelles.

Soit f : ℝ² → ℝⁿ telle que f et ${\partial f \over \partial x}$ soient continues sur ℝ², et soient a et b deux fonctions dérivables de ℝ dans ℝ. Alors, l'« intégrale paramétrique » (généralisée) F définie sur ℝ par :

$F(x)=\int _{a(x)}^{b(x)}f(x,y)~\mathrm {d} y$

est dérivable et

$F'(x)=f(x,b(x))b'(x)-f(x,a(x))a'(x)+\int _{a(x)}^{b(x)}{\partial f \over \partial x}(x,y)~\mathrm {d} y.$

Remarque : pour une fonction f qui ne dépend que de la seconde variable, on retrouve bien le théorème fondamental de l'analyse en posant a(x) = a et b(x) = x.

Démonstration

Soit H l'application de ℝ³ dans ℝ définie par :

$H(u,v,x)=\int _{u}^{v}f(x,y)~\mathrm {d} y.$

Du premier théorème fondamental de l'analyse, on déduit :

${\partial H \over \partial v}(u,v,x)=f(x,v)\quad {\text{et}}\quad {\partial H \over \partial u}(u,v,x)=-f(x,u),$

et en appliquant la règle de Leibniz, on a :

${\partial H \over \partial x}(u,v,x)=\int _{u}^{v}{\partial f \over \partial x}(x,y)~\mathrm {d} y.$

Comme $F(x)=H(a(x),b(x),x)$, le théorème de dérivation des fonctions composées donne :

$F'(x)={\partial H \over \partial v}(a(x),b(x),x)b'(x)+{\partial H \over \partial u}(a(x),b(x),x)a'(x)+{\partial H \over \partial x}(a(x),b(x),x).$

En remplaçant les trois dérivées partielles de $H$ par leurs valeurs respectives, on trouve l'identité annoncée.

Théorème de Fubini

Soient par exemple X une partie de ℝ^p, Y une partie de ℝ^q, et

$f:X\times Y\to \mathbb {R} ^{n}$

une application intégrable. Alors, d'après le théorème de Fubini, la fonction $f(x,\cdot )$ est intégrable pour presque tout x de X, l'intégrale paramétrique F définie par

$F(x)=\int _{Y}f(x,y)~\mathrm {d} y$

est intégrable sur X, et l'on a :

$\int _{X\times Y}f=\int _{X}F$

(et même chose en intervertissant les rôles de x et y).

Exemples de calcul

Calculs élémentaires

Exemple :

On peut vérifier en utilisant la règle de Leibniz que pour tous réels a et b strictement positifs :

$\int _{0}^{+\infty }{\exp(-ax)-\exp(-bx) \over x}\mathrm {d} x=\ln {b \over a}$.

Démonstration

Fixons a > 0, et soient F et g définies sur ]0,+∞[ par :

$F(b)=\int _{0}^{+\infty }{\exp(-ax)-\exp(-bx) \over x}\mathrm {d} x,\quad g(b)=\ln {b \over a}$.

On a clairement F(a) = g(a) = 0. Il suffit donc de montrer que leurs dérivées sont égales pour tout b > 0 pour vérifier l'identité.

En appliquant la règle de Leibniz pour F, on a :

$F'(b)=\int _{0}^{+\infty }{\partial \over \partial b}{\exp(-ax)-\exp(-bx) \over x}~\mathrm {d} x=\int _{0}^{+\infty }\exp(-bx)~\mathrm {d} x={1 \over b}=g'(b)$.

Exemple :

Soient X = [0 ; 2], Y = [1 ; 3] et f définie sur X × Y par f(x,y) = x² + y. Elle est intégrable sur X × Y puisqu'elle est continue. Par le théorème de Fubini, son intégrale se calcule donc de deux façons :

$\int _{X\times Y}f=\int _{0}^{2}\left(\int _{1}^{3}(x^{2}+y)~\mathrm {d} y\right)~\mathrm {d} x=\int _{0}^{2}(2x^{2}+4)~\mathrm {d} x={40 \over 3}$

et

$\int _{X\times Y}f=\int _{1}^{3}\left(\int _{0}^{2}(x^{2}+y)~\mathrm {d} x\right)~\mathrm {d} y=\int _{1}^{3}\left({8 \over 3}+2y\right)~\mathrm {d} y={40 \over 3}$.

Technique de Feynman

La technique de Feynman (Feynman's technique ou Feynman's integral trick) est une méthode de calcul intégral décrite par Richard Feynman, qui consiste à définir une intégrale comme une valeur particulière d'une fonction sous forme d'une intégrale paramétrique pour donner une expression plus simple de cette fonction et donner la valeur de l'intégrale comme valeur de la fonction évaluée en un point .

Exemples

Intégrale de Serret

On veut calculer l'intégrale :

$$I=\int _{0}^{1}{\frac {\ln(1+x)}{1+x^{2}}}\mathrm {d} x$$

Par la technique de Feynman, on considère la fonction :

$$I(a)=\int _{0}^{1}{\frac {\ln(1+ax)}{1+x^{2}}}\mathrm {d} x$$

On a alors I(0) = 0 et I(1) = I. On calcule la dérivée de la fonction, en dérivant sous le signe intégral puis en décomposant la fraction en éléments simples :

${\begin{aligned}I'(a)&=\int _{0}^{1}{\frac {\partial }{\partial a}}\left({\frac {\ln(1+ax)}{1+x^{2}}}\right)\mathrm {d} x\\&=\int _{0}^{1}{\frac {x}{(1+ax)(1+x^{2})}}\mathrm {d} x\\&={\frac {1}{1+a^{2}}}\int _{0}^{1}\left[{\frac {a}{1+x^{2}}}+{\frac {x}{1+x^{2}}}-{\frac {a}{1+ax}}\right]\mathrm {d} x\\&={\frac {1}{1+a^{2}}}\left[a\arctan(x)+{\frac {1}{2}}\ln(1+x^{2})-\ln(1+ax)\right]_{x=0}^{x=1}\\&={\frac {1}{1+a^{2}}}\left[a{\frac {\pi }{4}}+{\frac {1}{2}}\ln(2)-\ln(1+a)\right].\end{aligned}}$

On revient donc au calcul de l'intégrale par intégration :

${\begin{aligned}I=I(1)-I(0)=\int _{0}^{1}I'(a)\mathrm {d} a&=\int _{0}^{1}\left[{\frac {\pi }{4}}{\frac {a}{1+a^{2}}}+{\frac {1}{2}}\ln(2){\frac {1}{1+a^{2}}}-{\frac {\ln(1+a)}{1+a^{2}}}\right]\mathrm {d} a\\&=\left[{\frac {\pi }{4}}{\frac {1}{2}}\ln(1+a^{2})+{\frac {1}{2}}\ln(2)\arctan(a)\right]_{a=0}^{a=1}-I\\&={\frac {\pi }{4}}{\frac {1}{2}}\ln(2)+{\frac {1}{2}}\ln(2){\frac {\pi }{4}}-I.\end{aligned}}$

Ce qui permet de conclure $I={\frac {\pi \ln(2)}{8}}$.

Intégrale de Gauss

L'intégrale de Gauss joue un rôle important en analyse et en calcul des probabilités, elle est définie par :

$\int _{-\infty }^{+\infty }\mathrm {e} ^{-x^{2}}~\mathrm {d} x={\sqrt {\pi }}.$

Cette égalité peut s'obtenir de plusieurs façons, dont une faisant intervenir les intégrales paramétriques

$F(t)=\int _{0}^{t}\mathrm {e} ^{-x^{2}}\,\mathrm {d} x\quad {\text{et}}\quad G(t)=\int _{0}^{1}{\frac {\mathrm {e} ^{-t^{2}(1+x^{2})}}{1+x^{2}}}\,\mathrm {d} x$.

Voir aussi

Article connexe

Produit de convolution

Bibliographie

• Jean Mawhin, Analyse, fondements, techniques, évolution, De Boeck Université, 1997, 2^e éd., 808 p. (ISBN 978-2-8041-2489-2)
• (en) « Differentiation under the integral sign », sur PlanetMath