Eclats de vers : Matemat 12 : Analyse complexe

Index des Grimoires

Retour à l’accueil

Table des matières

\( \newcommand{\parentheses}[1]{\left(#1\right)} \newcommand{\crochets}[1]{\left[#1\right]} \newcommand{\accolades}[1]{\left\{#1\right\}} \newcommand{\ensemble}[1]{\left\{#1\right\}} \newcommand{\identite}{\mathrm{Id}} \newcommand{\indicatrice}{\boldsymbol{\delta}} \newcommand{\dirac}{\delta} \newcommand{\moinsun}{{-1}} \newcommand{\inverse}{\ddagger} \newcommand{\pinverse}{\dagger} \newcommand{\topologie}{\mathfrak{T}} \newcommand{\ferme}{\mathfrak{F}} \newcommand{\img}{\mathbf{i}} \newcommand{\binome}[2]{ \left\{ \begin{array}{c} #1 \\ #2 \\ \end{array} \right\} } \newcommand{\canonique}{\mathfrak{c}} \newcommand{\tenseuridentite}{\boldsymbol{\mathcal{I}}} \newcommand{\permutation}{\boldsymbol{\epsilon}} \newcommand{\matriceZero}{\mathfrak{0}} \newcommand{\matriceUn}{\mathfrak{1}} \newcommand{\christoffel}[2]{ \left\{ \begin{array}{c} #1 \\ #2 \\ \end{array} \right\} } \newcommand{\lagrangien}{\mathfrak{L}} \newcommand{\sousens}{\mathfrak{P}} \newcommand{\partition}{\mathrm{Partition}} \newcommand{\tribu}{\mathrm{Tribu}} \newcommand{\topologies}{\mathrm{Topo}} \newcommand{\setB}{\mathbb{B}} \newcommand{\setN}{\mathbb{N}} \newcommand{\setZ}{\mathbb{Z}} \newcommand{\setQ}{\mathbb{Q}} \newcommand{\setR}{\mathbb{R}} \newcommand{\setC}{\mathbb{C}} \newcommand{\corps}{\mathbb{K}} \newcommand{\boule}{\mathfrak{B}} \newcommand{\intervalleouvert}[2]{\relax \ ] #1 , #2 [ \ \relax} \newcommand{\intervallesemiouvertgauche}[2]{\relax \ ] #1 , #2 ]} \newcommand{\intervallesemiouvertdroite}[2]{[ #1 , #2 [ \ \relax} \newcommand{\fonction}{\mathbb{F}} \newcommand{\bijection}{\mathrm{Bij}} \newcommand{\polynome}{\mathrm{Poly}} \newcommand{\lineaire}{\mathrm{Lin}} \newcommand{\continue}{\mathrm{Cont}} \newcommand{\homeomorphisme}{\mathrm{Hom}} \newcommand{\etagee}{\mathrm{Etagee}} \newcommand{\lebesgue}{\mathrm{Leb}} \newcommand{\lipschitz}{\mathrm{Lip}} \newcommand{\suitek}{\mathrm{Suite}} \newcommand{\matrice}{\mathbb{M}} \newcommand{\krylov}{\mathrm{Krylov}} \newcommand{\tenseur}{\mathbb{T}} \newcommand{\essentiel}{\mathfrak{E}} \newcommand{\relation}{\mathrm{Rel}} \newcommand{\strictinferieur}{\ < \ } \newcommand{\strictsuperieur}{\ > \ } \newcommand{\ensinferieur}{\eqslantless} \newcommand{\enssuperieur}{\eqslantgtr} \newcommand{\esssuperieur}{\gtrsim} \newcommand{\essinferieur}{\lesssim} \newcommand{\essegal}{\eqsim} \newcommand{\union}{\ \cup \ } \newcommand{\intersection}{\ \cap \ } \newcommand{\opera}{\divideontimes} \newcommand{\autreaddition}{\boxplus} \newcommand{\autremultiplication}{\circledast} \newcommand{\commutateur}[2]{\left[ #1 , #2 \right]} \newcommand{\convolution}{\circledcirc} \newcommand{\correlation}{\ \natural \ } \newcommand{\diventiere}{\div} \newcommand{\modulo}{\bmod} \newcommand{\pgcd}{pgcd} \newcommand{\ppcm}{ppcm} \newcommand{\produitscalaire}[2]{\left\langle #1 \left|\right\relax #2 \right\rangle} \newcommand{\scalaire}[2]{\left\langle #1 \| #2 \right\rangle} \newcommand{\braket}[3]{\left\langle #1 \right| #2 \left| #3 \right\rangle} \newcommand{\orthogonal}{\bot} \newcommand{\forme}[2]{\left\langle #1 , #2 \right\rangle} \newcommand{\biforme}[3]{\left\langle #1 , #2 , #3 \right\rangle} \newcommand{\contraction}[3]{\left\langle #1 \odot #3 \right\rangle_{#2}} \newcommand{\dblecont}[5]{\left\langle #1 \right| #3 \left| #5 \right\rangle_{#2,#4}} \newcommand{\major}{major} \newcommand{\minor}{minor} \newcommand{\maxim}{maxim} \newcommand{\minim}{minim} \newcommand{\argument}{arg} \newcommand{\argmin}{arg\ min} \newcommand{\argmax}{arg\ max} \newcommand{\supessentiel}{ess\ sup} \newcommand{\infessentiel}{ess\ inf} \newcommand{\dual}{\star} \newcommand{\distance}{\mathfrak{dist}} \newcommand{\norme}[1]{\left\| #1 \right\|} \newcommand{\normetrois}[1]{\left|\left\| #1 \right\|\right|} \newcommand{\adh}{adh} \newcommand{\interieur}{int} \newcommand{\frontiere}{\partial} \newcommand{\image}{im} \newcommand{\domaine}{dom} \newcommand{\noyau}{ker} \newcommand{\support}{supp} \newcommand{\signe}{sign} \newcommand{\abs}[1]{\left| #1 \right|} \newcommand{\unsur}[1]{\frac{1}{#1}} \newcommand{\arrondisup}[1]{\lceil #1 \rceil} \newcommand{\arrondiinf}[1]{\lfloor #1 \rfloor} \newcommand{\conjugue}{conj} \newcommand{\conjaccent}[1]{\overline{#1}} \newcommand{\division}{division} \newcommand{\difference}{\boldsymbol{\Delta}} \newcommand{\differentielle}[2]{\mathfrak{D}^{#1}_{#2}} \newcommand{\OD}[2]{\frac{d #1}{d #2}} \newcommand{\OOD}[2]{\frac{d^2 #1}{d #2^2}} \newcommand{\NOD}[3]{\frac{d^{#3} #1}{d #2^{#3}}} \newcommand{\deriveepartielle}[2]{\frac{\partial #1}{\partial #2}} \newcommand{\dblederiveepartielle}[2]{\frac{\partial^2 #1}{\partial #2 \partial #2}} \newcommand{\dfdxdy}[3]{\frac{\partial^2 #1}{\partial #2 \partial #3}} \newcommand{\dfdxdx}[2]{\frac{\partial^2 #1}{\partial #2^2}} \newcommand{\gradient}{\mathbf{\nabla}} \newcommand{\combilin}[1]{\mathrm{span}\{ #1 \}} \newcommand{\trace}{tr} \newcommand{\proba}{\mathbb{P}} \newcommand{\probaof}[1]{\mathbb{P}\left[#1\right]} \newcommand{\esperof}[1]{\mathbb{E}\left[#1\right]} \newcommand{\cov}[2]{\mathrm{cov} \left( #1 , #2 \right) } \newcommand{\var}[1]{\mathrm{var} \left( #1 \right) } \newcommand{\rand}{\mathrm{rand}} \newcommand{\variation}[1]{\left\langle #1 \right\rangle} \newcommand{\composante}{comp} \newcommand{\bloc}{bloc} \newcommand{\ligne}{ligne} \newcommand{\colonne}{colonne} \newcommand{\diagonale}{diag} \newcommand{\matelementaire}{\mathrm{Elem}} \newcommand{\matpermutation}{permut} \newcommand{\matunitaire}{\mathrm{Unitaire}} \newcommand{\gaussjordan}{\mathrm{GaussJordan}} \newcommand{\householder}{\mathrm{Householder}} \newcommand{\rang}{rang} \newcommand{\schur}{\mathrm{Schur}} \newcommand{\singuliere}{\mathrm{DVS}} \newcommand{\convexe}{\mathrm{Convexe}} \newcommand{\petito}[1]{o\left(#1\right)} \newcommand{\grando}[1]{O\left(#1\right)} \)

1. Exponentielle complexe

1.1. Dépendances

  • Chapitre \ref{chap:complexe} : Les complexes

1.2. Introduction

AFAIRE : ARRANGER LE CHAPITRE

Considérons l'unique solution \(x\) de l'équation différentielle :

\( \OD{x}{t}(t) = z \cdot x(t) \\ x(0) = 1 \)

où \(z \in \setC\). On définit alors l'exponentielle d'un nombre complexe par :

\[\exp( z \cdot t ) = x(t)\]

On a donc simplement :

\[\exp(z) = x(1)\]

1.3. Additivité

Comme les fonctions :

\( s(t) = \exp( (z_1 + z_2) \cdot t ) \\ p(t) = \exp(z_1 t) \cdot \exp(z_2 t) \)

vérifient la même équation différentielle :

\( \OD{s}{t}(t) = (z_1 + z_2) \ s(t) \\ s(0) = 1 \)

\( \OD{p}{t}(t) = z_1 \ p(t) + z_2 \ p(t) = (z_1 + z_2) \ p(t) \\ p(0) = 1 \)

on a \(s(t) = p(t)\) pour tout \(t\). En considérant le cas \(t=1\), on arrive à la propriété d'additivité des exponentielles :

\[\exp(z_1 + z_2) = \exp(z_1) \cdot \exp(z_2)\]

1.4. Lien avec les fonctions trigonométriques

Considérons maintenant le cas particulier où \(u(t) = \exp(\img t)\) :

\( \OD{u}{t}(t) = \img u(t) \\ u(0) = 1 \)

On a alors :

\[\frac{d^2 u}{dt^2}(t) = - u(t) \qquad u(0) = 1 \qquad \OD{u}{t}(0) = \img\]

Soit la fonction \(u : \setR \mapsto \setR\) définie pour tout réel \(t\) par :

\[u(t) = \cos(t) + \img \sin(t)\]

On voit que :

\[\OD{u}{t} = -\sin(t) + \img \cos(t)\]

vérifie la même équation. Par unicité de la solution, on a :

\[\exp(\img t) = \cos(t) + \img \sin(t)\]

On en déduit directement que :

\( \exp(\img \pi/2) = \img \\ \exp(\img \pi) = -1 \\ \exp(3\pi\img/2) = \img \\ \exp(2\pi\img) = 1 \)

ainsi que la périodicité :

\[\exp(t+2\pi) = \exp(t)\]

pour tout \(t\in\setR\).

En utilisant l'additivité, on note que :

\[\exp(z) = \exp(a) \cdot \exp(\img b) = \exp(a) \cdot (\cos(b) + \img \sin(b))\]

De la même façon, si nous prenons \(u(t) = \exp(-i.t)\), nous en déduisons :

\( \frac{d^2 u}{dt^2}(t) = - u(t) \\ u(0) = 1 \qquad \OD{u}{t}(0) = -i \)

Donc :

\[\exp(-\img t) = \cos(t) - \img \sin(t)\]

En additionnant, puis en soustrayant les relations :

\( \exp( \img u ) = \cos(u) + \img \sin(u) \\ \exp( - \img u ) = \cos(u) - \img \sin(u) \)

on obtient :

\( \cos(u) = \frac{ \exp(\img u) + \exp(-\img u) }{2} \\ \sin(u) = \frac{ \exp(\img u) - \exp(-\img u) }{2 \img} \)

Choisissant un angle \(\theta\in\setR\) tel que :

\( \cos(\theta) = \frac{\Re(z)}{\abs{z}} \\ \sin(\theta) = \frac{\Im(z)}{\abs{z}} \)

on peut réexprimer \(z\) comme :

\[z = \abs{z}(\cos(\theta)+\img\sin(\theta)) = \abs{z}\exp(\img\theta)\]

Comme les fonctions \(\cos\) sont \(2\pi\) périodiques, il existe une infinité d'angles \(\theta\) vérifiant cette propriété. On définit l'argument de \(z\) comme l'unique \(\theta\) vérifiant cette propriété et se trouvant dans l'intervalle \([0,2\pi)\).

\( θ = arg(z) \quad⇔\quad

\begin{cases} z = \abs{z}\exp(\img\theta) \\ \theta\in [0,2\pi) \end{cases}

\)

Inspiré par relation :

\[z = \abs{z}\exp(\img\arg(z))\]

et cherchant à étendre la propriété :

\[\ln(x \cdot y) = \ln(x) + \ln(y)\]

du logarithme sur \(\setR\), on définit le logarithme d'un complexe par :

\( \ln(z) = \ln(\abs{z}) + \ln(\exp(\img\arg(z))) \\ \ln(z) = \ln(\abs{z}) + \img\arg(z) \)

Pour un \(z \in \setC\) donné, l'ensemble des \(y \in \setC\) tels que :

\[\exp(y) = z\]

peut s'écrire :

\[\mathcal{Y} = \{ y_k = \ln(\abs{z}) + \img \arg(z) + 2 \pi \img k : k \in \setZ \}\]

Les complexes permettent de retrouver aisément les propriétés fondamentales des fonctions trigonométriques \(\cos\), \(\sin\). En se rappelant que :

\[\exp(-\img t) \cdot \exp(i t) = \exp(i t - i t) = \exp(0) = 1\]

on retombe sur :

\( ( \cos(t) + \img \sin(t) ) ( \cos(t) - \img \sin(t) ) = 1 \\ \cos(t)^2 + \sin(t)^2 = 1 \)

En dérivant la relation reliant l'exponentielle aux fonctions trigonométriques, on obtient :

\( \OD{}{t}( \cos(t) + \img \sin(t) ) = \OD{}{t}\exp(\img t) \\ \OD{}{t} \cos(t) + \img \OD{}{t}\sin(t) = \img \exp(\img t) \\ \OD{}{t} \cos(t) + \img \OD{}{t}\sin(t) = \img \cos(t) - \sin(t) \)

on a bien :

\( \OD{}{t} \cos(t) = -\sin(t) \\ \OD{}{t}\sin(t) = \cos(t) \)

On a aussi :

\( \cos(u+v) + \img \sin(u+v) = \exp\left[ \img ( u + v ) \right] \\ \cos(u+v) + \img \sin(u+v) = \exp\left[ \img u \right] \cdot \exp\left[ \img v \right] \\ \cos(u+v) + \img \sin(u+v) = (\cos(u) + \img \sin(u))(\cos(v) + \img \sin(v)) \)

ce qui donne, tous calculs faits :

\( \cos(u+v) = \cos(u)\cos(v) - \sin(u)\sin(v) \\ \sin(u+v) = \sin(u)\cos(v) + \sin(v)\cos(u) \)

2. Dérivation et intégration dans le plan complexe

2.1. Dépendances

  • Chapitre \ref{chap:differ} : Les différentielles
  • Chapitre \ref{chap:integ} : Les intégrales

2.2. Dérivées

On définit la dérivée d'une fonction \(f : \setC \mapsto \setC\) par analogie avec la dérivée des fonctions réelles :

\[\OD{f}{z}(z) = \lim_{ \substack{ h \to 0 \\ h \in\setC } } \frac{f(z+h) - f(z)}{h}\]

Si \(f\) est dérivable sur \(A\subseteq\setC\), on dit qu'elle est analytique sur \(A\).

On définit les fonctions \(u,v : \setR^2 \mapsto \setR\) associées à \(f\) par :

\( u(x,y) = \Re\left(f(x + \img y) \right) \\ v(x,y) = \Im\left(f(x + \img y) \right) \)

On a alors \(f(z) = u(x,y) + \img v(x,y)\). Si \(f\) est dérivable, la limite ne peut pas dépendre du chemin suivi pour arriver en \((x,y)\). On peut donc choisir \(h\) de la forme \(\Delta x\) ou \(\img \Delta y\) avec \(\Delta x,\Delta y\in\setR\). Cela donne :

\( \OD{f}{z}(z) = \lim_{\Delta x \to 0} \frac{\PD{u}{x}\Delta x + \img \PD{v}{x}\Delta x}{\Delta x} \\ \OD{f}{z}(z) = \lim_{\Delta y \to 0} \frac{\PD{u}{y}\Delta y + \img \PD{v}{y}\Delta y}{\img \Delta y} \)

et par comparaison :

\( \PD{u}{x} = \PD{v}{y} \\ \PD{u}{y} = -\PD{v}{x} \)

En dérivant de nouveau ces équations, on obtient :

\( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0 \\ \frac{\partial^2 v}{\partial x^2} + \frac{\partial^2 v}{\partial y^2} = 0 \)

Les puissances de \(z\) sont dérivables. On vérifie que :

\[\OD{z^n}{z}=n z^{n-1}\]

2.3. Extension de l'intégrale au plan complexe

Soit une fonction \(f : A \mapsto \setC\). Comme \(u = \Re(f)\) et \(v = \Im(f)\) sont des fonctions à valeurs réelles :

\[u,v : A \mapsto \setR\]

leur intégrale est bien définie. On définit alors l'intégrale de \(f\) par :

\[\int_A f d\mu = \int_A u d\mu + \img \int_A v d\mu\]

2.4. Intégrale de ligne

Soit une fonction \(\gamma : [a,b] \mapsto \setC\) et la courbe :

\[l = \{ \gamma(t) : t \in [a,b] \} \subset \setC\]

On définit l'intégrale de ligne de \(f\) sur \(l\) par :

\[\int_l f(z) dz = \int_a^b (f \circ \gamma)(t) \OD{\gamma}{t}(t) dt\]

On peut également écrire :

\( \int_l f(z) dz = \int_l (u + \img v)(dx + \img dy) \\ \int_l f(z) dz = \int_l (u dx - v dy) + \img \int_l (u dy + v dx) \)

où les intégrales des membres de droites sont des intégrales de ligne classiques de fonctions \(\setR^2 \mapsto \setR\).

2.4.1. Théorème fondamental

Soit une fonction analytique \(F : \setC \mapsto \setC\). La dérivée d'une composée de fonction nous donne :

\[\OD{}{t}(F \circ \gamma)(t) = (\OD{F}{z} \circ \gamma)(t) \OD{\gamma}{t}(t)\]

On en déduit, en appliquant le théorème fondamental que :

\[\int_l \OD{F}{z}(z) dz = F(\gamma(b))-F(\gamma(a))\]

2.4.2. Contour fermé

Dans la suite, lorsque \(l\) est une courbe fermée et est donc la frontière d'un certain ensemble \(l = \partial S\) inclus dans \(\setC\), nous notons :

\[\oint_l f(z) dz = \int_{\partial S} f(z) dz\]

Soit \(f\) une fonction analytique et \(l = \partial S\) une courbe fermée. On a :

\[\oint_l f(z) dz = \oint_l (u dx - v dy) + \img \oint_l (u dy + v dx)\]

En appliquant le théorème de Stokes, il vient :

\( \oint_l f(z) dz = \int_S (-\PD{v}{x}-\PD{u}{y}) dS + \img \int_S (\PD{u}{x}-\PD{v}{y}) dS \\ \)

Mais comme \(f\) est analytique :

\( \PD{u}{x} = \PD{v}{y} \\ \PD{u}{y} = -\PD{v}{x} \)

Les termes du membre de droite s'annulent donc et :

\[\oint_l f(z) dz = 0\]

2.4.3. Cercle

Soit \(C\) le cercle de centre \(a\) et de rayon \(R\) :

\[C = \{ \gamma(\theta) = a + R \exp(\img\theta) : \theta \in [ 0 , 2 \pi ) \}\]

On a alors :

\[\OD{\gamma}{\theta}(\theta) = R \img \exp(\img\theta)\]

et :

\[(z - a)^k = \left( R \exp(\img\theta) \right)^k\]

pour tout \(z\in C\) et \(k \in \setZ\). Donc :

\( \oint_C \frac{dz}{(z-a)^k} = \int_0^{2\pi} \frac{R \img \exp(\img\theta)}{R^k \exp(\img k\theta)} \\ \oint_C \frac{dz}{(z-a)^k} = \img R^{1-k} \int_0^{2\pi} \exp[\img(1-k)\theta] d\theta \)

Mais les propriétés de périodicité des fonctions trigonométriques et de l'exponentielle associée nous permettent de vérifier que :

\[\int_0^{2\pi} \exp[\img l \theta] d\theta = 2 \pi \delta_{l,0}\]

pour tout \(l\in\setZ\). On a donc :

\[\oint_C \frac{dz}{(z-a)^k} = 2 \pi \img \delta_{k,1}\]

pour tout \(k \in \setZ\). L'intégrale s'annule pour tout les \(k \ne 1\).

Vu que la fonction \(f(z) = (z-a)^{-k}\) est analytique partout sauf en \(a\), on obtient le même résultat pour toute courbe fermée entourant \(a\) :

\[\oint_{\partial S} \frac{dz}{(z-a)^k} = 2 \pi \img \delta_{k,1} \delta_S(a)\]

2.4.4. Théorème de Cauchy

Si \(f\) est analytique, la fonction

\( g(z) =

\begin{cases} \frac{f(z)-f(a)}{z-a} & \mbox{si } z \ne a \\ \OD{f}{z}(a) & \mbox{si } z = a \end{cases}

\)

est analytique aussi. On a donc :

\[\oint_l g(z) dz = 0\]

Soit \(a\in\setC\) et une courbe fermée entourant \(a\) telle que \(a\notin l\). On a alors :

\[\oint_l \frac{f(z)-f(a)}{z-a} dz = \oint_l g(z) dz = 0\]

On en déduit que :

\[f(a) \oint_{\partial S} \frac{1}{z-a} dz = \oint_{\partial S} \frac{f(z)}{z-a} dz\]

ce qui nous donne la valeur de \(f(a)\) en fonction d'une intégrale :

\[f(a) = \frac{1}{2 \pi \img} \oint_{\partial S} \frac{f(z)}{z-a} dz\]

En dérivant \(k\) fois cette dernière relation par rapport à \(a\), on obtient :

\[\OD{f}{z}(a) = \frac{k !}{2 \pi \img} \oint_{\partial S} \frac{f(z)}{(z-a)^{k+1}} dz\]

2.4.5. Séries de Laurent

Supposons que \(f\) puisse s'écrire comme une combinaison linéaire de puissances entière de \(z-a\). On a :

\[f(z) = \sum_{n=-\infty}^{+\infty} \alpha_n (z-a)^n\]

Multipliant le tout par \((z-a)^{-k-1}\) et intégrant sur un contour entourant \(a\), on obtient :

\[\oint_l f(z)(z-a)^{-k-1} dz = \sum_{n=-\infty}^{+\infty} \alpha_n \oint_l (z-a)^{n-k-1} dz\]

Mais nous avons vu que seule l'intégrale de \((z-a)^{-1}\) ne s'annule pas. On en déduit que :

\[\alpha_k = \frac{1}{2 \pi \img} \oint_l \frac{f(z)}{(z-a)^{k+1}} dz\]

2.4.6. Théorème des résidus

Considérons le cas où :

\[f(z) = \sum_{n=-p}^{+\infty} \alpha_n (z-a)^n\]

pour un certain \(p \in \setZ\). On dit alors que \(f\) à un pôle d'ordre \(p\) en \(a\). En multpliant le tout par \((z-a)^p\) et en dérivant \(p-1\) fois, on obtient :

\[\frac{d^{p-1}}{dz^{p-1}}[(z-a)^p f(z)] = (p-1) ! \left[a_{-1} + a_0 (z-a) + a_1 (z-a)^2 + ...\right]\]

On a donc :

\[\lim_{z \to a} \frac{d^{p-1}}{dz^{p-1}}[(z-a)^p f(z)] = (p-1) ! a_{-1}\]

Mais comme \(a_{-1}\) n'est rien d'autre, à un facteur \(2\pi\img\) près, que l'intégrale de \(f\) sur un contour entourant \(a\), on a : \[\] \ointl f(z) dz = 2 π \img limz → a \frac{1}{(p-1)!} \frac{dp-1}{dzp-1}[(z-a)p f(z)]

3. Polynômes et exponentielles

\label{chap:polyexpo}

3.1. Chebyshev

Les polynômes de Chebyshev \(T_n\) sont défini par :

\( T_n : x\mapsto\cos(n \arccos(x) ) \\ T_n(\cos(\theta)) = \cos(n \theta) \)

Il est clair que :

\( T_0(x) = 1 \\ T_1(x) = x \)

Les propriétés des fonctions trigonométriques :

\( \cos[(n+m)\theta]= \cos[n\theta]\cos[m\theta]+\sin[n\theta]\sin[m\theta] \\ \cos[(n-m)\theta]= \cos[n\theta]\cos[m\theta]-\sin[n\theta]\sin[m\theta] \)

nous montrent que :

\[T_{n+m}(x) + T_{n-m}(x) = 2 T_n(x) T_m(x)\]

On en déduit entre autre que :

\[T_{n+1}(x) + T_{n-1}(x) = 2 x T_n(x)\]

Les fonctions \(T_n\) sont donc des polynômes. Comme :

\[T_n \left[ \cos\left(\frac{2k+1}{2n}\right) \right] = 0\]

les racines sont données par :

\[x_k = \cos\left(\frac{2k+1}{2n}\right)\]

pour \(k=0,1,...,n-1\). On déduit aussi la propriété suivante de la définition :

\[(T_n \circ T_m)(x) = T_{m \cdot n}(x) = (T_m \circ T_n)(x)\]

Nous allons voir que les \(T_n\) sont orthogonaux moyennant un certain poids. Si \(m = n = 0\), il est clair que :

\[\int_0^\pi T_0(\cos(\theta)) T_0(\cos(\theta)) d\theta = \int_0^\pi d\theta = \pi\]

Si \(m, n \ne 0\), on a :

\( \int_0^\pi T_m(\cos(\theta)) T_n(\cos(\theta)) d\theta = \int_0^\pi \cos(m\theta) \cos(n\theta) d\theta \\ \int_0^\pi T_m(\cos(\theta)) T_n(\cos(\theta)) d\theta = \unsur{2} \int_0^\pi \cos[(m+n)\theta] d\theta + \unsur{2} \int_0^\pi \cos[(m-n)\theta] d\theta \)

Si \(m = n\), on a alors :

\[\int_0^\pi T_m(\cos(\theta)) T_n(\cos(\theta)) d\theta = \unsur{2n} [\sin(2 n \pi) - \sin(0)] + \frac{\pi}{2}\]

La première intégrale s'annule donc par périodicité de la fonction \(\sin\). A présent, si \(m \ne n\), on a :

\( \int_0^\pi T_m(\cos(\theta)) T_n(\cos(\theta)) d\theta = \\ \unsur{2(m+n)} [\sin[(m+n)\pi] - \sin(0)] + \\ \unsur{2(m-n)} [\sin[(m-n)\pi] - \sin(0)] = 0 \)

Rassemblant ces résultats, on obtient :

\( ∫0^π Tm(cosθ) Tn(cosθ) dθ =

\begin{cases} \pi & \mbox{si } m = 0 \\ \pi/2 & \mbox{si } m = n \ne 0 \\ 0 & \mbox{si } m \ne n \end{cases}

\)

Le changement de variable \(x = \cos(\theta)\),

\( dx = -\sin(\theta) d\theta \\ d\theta = \frac{dx}{\sqrt{1-x^2}} \)

nous donne :

\[\int_0^\pi T_m(\cos\theta) T_n(\cos\theta) d\theta = \int_{-1}^1 \frac{T_m(x) T_n(x)}{\sqrt{1-x^2}}dx\]

On en déduit que :

\[\int_{-1}^1 \frac{T_m(x) T_n(x)}{\sqrt{1-x^2}} dx = \delta_{mn} \frac{\pi}{2} (1+\delta_{m,0})\]

pour tout \(m,n\in\setN\).

3.2. Hermite

Les polynômes de Hermite sont orthogonaux pour le produit scalaire :

\[\int_{-\infty}^{+\infty} H_n(x) H_m(x) \exp(-x^2) dx = 2^n n ! \sqrt{\pi} \delta_{mn}\]

Ils obéissent à la récurrence :

\( H_0(x) = 1 \\ H_1(x) = 2 x \\ H_{n+1}(x) = 2 x H_n(x) - 2 n H_{n-1}(x) \\ \)

3.3. Laguerre

Les polynômes de Laguerre sont orthogonaux pour le produit scalaire :

\[\int_{0}^{+\infty} L_n(x) L_m(x) \exp(-x) dx = \delta_{mn}\]

Ils obéissent à la récurrence :

\( L_0(x) = 1 \\ L_1(x) = 1-x \\ (n+1) L_{n+1}(x) = (2 n + 1 - x) L_n(x) - n L_{n-1}(x) \)

4. Analyse de Fourier

4.1. Séries de Fourier

ARRANGER LE CHAPITRE

Constatons tout d'abord que :

\( ∫0 exp(\img k x) dx =

\begin{cases} 2 \pi & \mbox{si } k = 0 \\ 0 & \mbox{si } k \ne 0 \end{cases}

\)

pour tout \(k\in\setZ\). En effet, si \(k=0\), l'intégrale devient :

\[\int_0^{2\pi} \exp(0) dx = \int_0^{2\pi} 1 dx = 2 \pi\]

Par contre, si \(k\ne 0\), le changement de variable :

\( s = \img k x \\ ds = \img k dx \)

nous mène à :

\[\int_0^{2\pi} \exp(\img k x) dx = \unsur{\img k} \int_0^{2\pi\img k} \exp(s) ds\]

le théorème fondamental nous donne alors :

\[\int_0^{2\pi} \exp(\img k x) dx = \unsur{\img k} [\exp(2\pi\img k)-\exp(0)] = 0\]

pour tout \(k\in\setZ\), par périodicité des exponentielles imaginaires. On généralise ce résultat en utilisant le changement de variable :

\[t = a + \frac{T}{2\pi} x\]

ce qui nous donne :

\( ∫aa+T exp\left(\frac{2\pi\img k x}{T}\right) dx =

\begin{cases} T & \mbox{si } k = 0 \\ 0 & \mbox{si } k \ne 0 \end{cases}

\)

Définissons à présent le produit scalaire :

\[\scalaire{u}{v} = \int\limits_{-T/2}^{T/2} \bar{u}(x) v(x) dx\]

Soient les fonctions \(e_k\) :

\[e_k : x \mapsto \exp\left(\frac{2\pi\img k x}{T}\right)\]

où \(k\in\setZ\). Notons que :

\( \bar{e}_m e_n = \exp\left(-\frac{2\pi\img m x}{T}\right)\exp\left(\frac{2\pi\img n x}{T}\right) \\ \bar{e}_m e_n = \exp\left(\frac{2\pi\img (n-m) x}{T}\right) \)

La suite des \(e_k\) est donc orthonormée :

\[\scalaire{e_m}{e_n} = \int\limits_{-T/2}^{T/2} \exp\left(\frac{2\pi\img (n-m) x}{T}\right)dx = 2 \pi \delta_{mn}\]

Soit à présent \(u\in\ev{e_k : k\in\setZ}\) :

\[u(x) = \sum_{k=-\infty}^{+\infty} \hat{u}_k e_k\]

Nous avons vu au chapitre \ref{chap:vector} que les composantes d'un vecteur (ici la fonction \(u\)) par rapport à une base orthonormée s'écrivent :

\[\hat{u}_k = \scalaire{e_k}{u} = \int\limits_{-T/2}^{T/2} u(x) \exp\left(-\frac{2\pi\img k x}{T}\right)dx\]

Soit à présent \(v\) :

\[v(x) = \sum_{k=-\infty}^{+\infty} \hat{v}_k e_k\]

Utilisant à nouveau les propriétés des bases orthonormées, on a :

\[\scalaire{u}{v} = \sum_{k=-\infty}^{+\infty} \conj(\hat{u}_k) \hat{v}_k\]

4.2. Transformeée de Fourier discrète

Dans le cas où :

\( a = \exp(\img x) \\ a^k = \exp(\img k x) \)

l'équation de la progression géométrique du chapitre \ref{chap:real} :

\[\sum_{i=0}^{n} a^i = \frac{1-a^{n+1}}{1-a}\]

devient :

\[\sum_{k=0}^{n} \exp(\img k x) = \frac{1-\exp(\img (n+1) x)}{1-\exp(\img x)}\]

On voit alors que si \((n+1) x\) est un multiple entier de \(2\pi\), la somme s'annule par périodicité sauf dans le cas où \(x\) est un multiple entier de \(2\pi\), où elle vaut \(n+1\). Posant \(N = n - 1\) et \(x = 2\pi l/N\), on obtient :

\( ∑k=0N-1 exp \left( \frac{2 \pi \img k l}{N} \right) =

\begin{cases} N & \mbox{si } m = l/N \in\setZ \\ 0 & \mbox{si } m = l/N \notin\setZ \end{cases}

\)

pour tout \(k,l,m,N \in \setZ\). Dans la suite nous utilisons la notation :

\[e(k,l) = \exp \left( \frac{2 \pi \img k l}{N} \right)\]

Nous allons en déduire la forme de la transformeée de Fourier discrète. Soit la série \(u_k : k = 0,...,N-1\) :

\[u_k = \sum_{n=-N/2}^{N/2-1} U_n e(k,n)\]

où nous supposons que \(N/2 \in\setN\). Multipliant par \(e(k,-m)\) où \(m \in \{-N/2,...,N/2-1\}\), et utilisant l'additivité des exponentielles, il vient :

\( u_k e(k,-m) = \sum_{n=-N/2}^{N/2-1} U_n e(k,n) e(k,-m) \\ u_k e(k,-m) = \sum_{n=-N/2}^{N/2-1} U_n e(k,n-m) \)

Sommons à présent sur \(k\) :

\[\sum_{k=0}^{N-1} u_k e(k,-m) = \sum_{n=-N/2}^{N/2-1} U_n \sum_{k=0}^{N-1} e(k,n-m)\]

Mais la somme sur \(k\) du membre de droite s'annule sauf lorsque \(n-m\) est multiple entier de \(N\). Comme \(m,n\) sont dans \(\{-N/2,...,N/2-1\}\), le seul cas possible est ici \(m = n\). On a donc :

\[U_m = \unsur{N} \sum_{k=0}^{N-1} u_k e(k,-m)\]

Nous avons donc obtenu une bijection entre les \(u_k\) et les \(U_n\) :

\( u_k = \sum_{n=-N/2}^{N/2-1} U_n \exp\left(\frac{2 \pi \img k n}{N}\right) \\ U_n = \unsur{N} \sum_{k=0}^{N-1} u_k \exp\left(-\frac{2 \pi \img k n}{N}\right) \)

4.3. Transformeée de Fourier

On peut présenter le résultat précédent sous une autre forme. On peut poser :

\( x_k = k\Delta x \\ y_n = n\Delta y \)

ainsi que :

\( u_k = u(x_k) \\ U_n = U(y_n) \)

pour des fonctions \(u,v \in X \subset \fonction(\setR,\setR)\). Soit :

\[u(x_k) = \sum_{n=-N/2}^{N/2-1} U(y_n) \exp(2 \pi \img x_k y_n) \Delta y\]

où \(k = -N/2 , ... ,N/2 - 1\). Suivant le même procédé que ci-dessus, on aboutit à la relation inverse :

\[U(y_n) = \sum_{k = -N/2}^{N/2-1} u(x_k) \exp(-2 \pi \img x_k y_n) \Delta x\]

pour autant que :

\[\Delta x \Delta y = \frac{1}{N}\]

Choisissons :

\[\Delta x = \Delta y = \frac{1}{\sqrt{N}}\]

On voit que \(\Delta x\) et \(\Delta y\) tendent alors vers \(0\) lorsque \(N\) tend vers l'infini. Les sommes ci-dessus se rapprochent donc de plus en plus d'intégrales sur l'intervalle \([-\sqrt{N}/2,\sqrt{N}/2]\), qui tend lui-même vers \((-\infty,+\infty) = \setR\). On arrive donc aux relations :

\( u(x) = \int\limits_{-\infty}^{+\infty} U(y) \exp(2\pi\img x y) dy \\ U(y) = \int\limits_{-\infty}^{+\infty} u(x) \exp(-2\pi\img x y) dx \)

On définit alors la transformeée de Fourier \(\mathcal{F} : X \mapsto X\) et son inverse par :

\( U(y) = \mathcal{F}(u)(y) = \int\limits_{-\infty}^{+\infty} u(x) \exp(-2\pi\img x y) dx \\ u(x) = \mathcal{F}^{-1}(U)(x) = \int\limits_{-\infty}^{+\infty} U(y) \exp(2\pi\img x y) dy \)

4.3.1. Delta de Dirac

On déduit des relations ci dessus que :

\( u(x) = \int\limits_{-\infty}^{+\infty} \exp(2\pi\img x y) dy \int\limits_{-\infty}^{+\infty} u(z) \exp(-2\pi\img z y) dz \\ u(x) = \int\limits_{-\infty}^{+\infty} u(z) dz \int\limits_{-\infty}^{+\infty} \exp(2\pi\img (x-z) y) dy \)

On en déduit la relation fondamentale :

\[\int\limits_{-\infty}^{+\infty} \exp(2\pi\img (x-z) y) dy = \delta(x-z)\]

qui est l'extension de l'orthonormalité des bases discrètes. Ce n'est donc pas une intégrale au sens classique du terme, mais une distribution.

Cette relation nous montre aussi, lorsque \(z=0\), que :

\[\mathcal{F}(1)(x) = \delta(x)\]

Inversément, on a :

\[\mathcal{F}(\delta)(x) = \int\limits_{-\infty}^{+\infty} \delta(x) \exp(-2\pi\img x y) dx = \exp(0) = 1\]

4.3.2. Produit scalaire

Considérons à présent le produit scalaire :

\[\scalaire{u}{v} = \int\limits_{-\infty}^{+\infty} \bar{u}(x) v(x) dx\]

et examinons \(\scalaire{\mathcal{F}(u)}{\mathcal{F}(v)}\). En utilisant la propriété d'orthonormalité, on arrive à :

\( \int\limits_{-\infty}^{+\infty} \conj\left[\mathcal{F}(u)(y)\right] \mathcal{F}(v)(y) dy = \int_{\setR^2} \bar{u}(x) v(z) \exp\left(2\pi\img (x-z) y\right) dx dz \\ \int\limits_{-\infty}^{+\infty} \conj\left[\mathcal{F}(u)(y)\right] \mathcal{F}(v)(y) dy = \int_{\setR^2} \bar{u}(x) v(z) \delta(x-z) dx dz \)

et donc :

\( \int\limits_{-\infty}^{+\infty} \conj\left[\mathcal{F}(u)(y)\right] \mathcal{F}(v)(y) dy = \int\limits_{-\infty}^{+\infty} \bar{u}(x) v(x) dx \)

La transformeée de Fourier possède donc la propriété de conserver le produit scalaire :

\[\scalaire{\mathcal{F}(u)}{\mathcal{F}(v)} = \scalaire{u}{v}\]

4.3.3. Convolution

La transformeèe d'un produit de convolution s'écrit :

\( \mathcal{F}(u \star v)(y) = \int\limits_{-\infty}^{+\infty} \exp(-2\pi\img x y) dx \int\limits_{-\infty}^{+\infty} u(x-z) v(z) dz \)

Considérons le changement de variable :

\( \xi = x - z \\ \eta = z \)

On a alors \(x = \xi + \eta\) et :

\( \mathcal{F}(u \star v)(y) = \int\limits_{-\infty}^{+\infty} u(\xi) \exp(-2\pi\img \xi y) d\xi \int\limits_{-\infty}^{+\infty} v(\eta) \exp(-2\pi\img \eta y) d\eta \)

c'est-à-dire :

\[\mathcal{F}(u \star v) = \mathcal{F}(u)\mathcal{F}(v)\]

4.3.4. Dérivées

Soit \(u\) une fonction qui s'annule à l'infini. La transformeée de Fourier de sa dérivée s'écrit :

\( \mathcal{F}\left(\OD{u}{x}\right)(y) = \int\limits_{-\infty}^{+\infty} \OD{u}{x}(x)\exp(-2\pi\img x y) dx \)

En intégrant par partie, et en tenant compte du fait que les limites à l'infini sont nulles, il vient :

\( \mathcal{F}\left(\OD{u}{x}\right)(y) =

  • (-2π\img y) ∫\limits-∞+∞ u(x)exp(-2π\img x y) dx

\)

et finalement :

\[\mathcal{F}\left(\OD{u}{x}\right)(y) = 2 \pi \img y \mathcal{F}(y)\]

4.3.5. Déphasage

Considérons l'opérateur de translation :

\[t_a(u)(x) = u(x-a)\]

La transformeée de Fourier d'une translatée de la fonction \(u\) s'écrit :

\( \mathcal{F}\left[t_a(u)\right](y) = \int\limits_{-\infty}^{+\infty} u(x-a) \exp(-2\pi\img x y) dx \\ \mathcal{F}\left[t_a(u)\right](y) = \exp(-2\pi\img a y) \int\limits_{-\infty}^{+\infty} u(x-a) \exp(-2\pi\img (x-a) y) dx \)

Mais comme \(dx = d(x-a)\), on en déduit que :

\[\mathcal{F}\left[t_a(u)\right](y) = \exp(-2\pi\img a y)\mathcal{F}(u)(y)\]

4.3.6. Dilatation

Soit \(d_a\) l'opérateur de dilatation :

\[d_a(u)(x) = u(a x)\]

où \(a > 0\) est un réel strictement positif.

La transformeée de Fourier de la fonction dilaté \(d_a(u)\) s'écrit :

\( \mathcal{F}\left[d_a(u)\right](y) = \int\limits_{-\infty}^{+\infty} u(a x) \exp(-2\pi\img x y) dx \)

Considérons le changement de variable \(z = a x\). On a alors :

\( \mathcal{F}\left[d_a(u)\right](y) = \unsur{a}\int\limits_{-\infty}^{+\infty} u(z) \exp\left(-2\pi\img z \frac{y}{a}\right) dz \\ \mathcal{F}\left[d_a(u)\right](y) = \unsur{a}\mathcal{F}(u)(y/a) \)

c'est-à-dire :

\[\mathcal{F}\left[d_a(u)\right] = \unsur{a}\mathcal{F}\left[d_{1/a}(u)\right]\]

4.3.7. Gaussienne

Définissons la gaussienne \(G\) :

\[G(x) = \exp(-\pi x^2)\]

Sa transformeée s'écrit :

\[\mathcal{F}(G)(y) = \int\limits_{-\infty}^{+\infty} \exp(-\pi x^2 - 2\pi\img x y) dx\]

Il est clair que :

\( \exp\left[ -\pi (x+\img y)^2 \right] = \exp(-\pi x^2 - 2 \pi\img x y)\exp(\pi y^2) \)

Il vient donc, en effectuant le changement de variable \(\xi = x + \img y\) :

\[\mathcal{F}(G)(y) = \exp(-\pi y^2) \int\limits_{-\infty}^{+\infty} \exp(-\pi \xi^2) d\xi\]

Mais nous avons vu au chapitre \ref{chap:integ} que cette dernière intégrale vaut \(1\). Donc :

\[\mathcal{F}(G) = G\]

La transformeée de Fourier laisse \(G\) inchangée.

4.4. Transformeée de Fourier

La transformeée de Fourier (voir chapitre \ref{chap:fourier}) est définie par :

\[\mathcal{F}(u)(y) = \int_{\setR} u(x) \exp(2\pi \img x y) dx\]

où \(\img = \sqrt{-1}\). On a :

\( \int_{\setR} \mathcal{F}(u)(y) v(y) dy = \int_{\setR} u(x) v(y) \exp(2\pi \img x y) dx dy \\ \int_{\setR} \mathcal{F}(u)(y) v(y) dy = \int_{\setR} u(x) \mathcal{F}(v)(x) dx \)

et donc :

\[\forme{\mathcal{F}(u)}{v} = \forme{u}{\mathcal{F}(v)}\]

5. Ondelettes

\label{chap:ondelet}

Les ondelettes constitue un outil d'analyse de signal puissant. L'idée est de décomposer une fonction en différentes échelles, des échelles grossières jusqu'à des résolution très fines. Contrairement aux séries de Fourier, elles ont une influence locale et des relations simples entres les différentes échelles représentées, ce qui permet de passer rapidement (en \(\mathcal{O}(N)\)) d'une représentation à une autre.

Elles sont fort utilisées en compression de donnée (son, image, …) où les coefficients les plus négligeables sont ignorés.

5.1. Introduction

AFAIRE : ARRANGER LE CHAPITRE

Le plus simple est de commencer avec un exemple célèbre : les ondelettes de Haar :

\( \varphi(x) = \indicatrice_{[0,1[}(x) \\ \psi(x) = \indicatrice_{[0,1/2[}(x) - \indicatrice_{[1/2,1[}(x) \)

La fonction \(\varphi\) est appelée fonction d'échelle, tandis que \(\psi\) est l'ondelette proprement dite. La propriété principale d'un système d'ondelette est qu'il permet facilement de changer d'échelle. Ainsi, on a :

\( \varphi(x) = \varphi(2x) + \varphi(2x-1) = a_0^0 \varphi(2x) + a_1^0 \varphi(2x-1) \\ \psi(x) = \varphi(2x) - \varphi(2x-1) = a_0^1 \varphi(2x) + a_1^1 \varphi(2x-1) \)

où l'on a posé \(a_0^0 = 1\), \(a_1^0 = 1\), \(a_0^1 = 1\), \(a_1^1 = -1\). En inversant les deux relations ci-dessus, on arrive à :

\( \varphi(2x) = \frac{1}{2} ( \varphi(x) + \psi(x) ) \\ \varphi(2x-1) = \frac{1}{2} ( \varphi(x) - \psi(x) ) \)

Les deux représentations ci-dessous sont donc équivalentes :

\( f(x) = f_1 \varphi(2x) + f_2 \varphi(2x-1) \\ f(x) = \frac{f_1+f_2}{2} \varphi(x) + \frac{f_1 - f_2}{2} \psi(x) \)

La différence est que l'on représente \(f\) tantôt comme une somme de fonctions de base de même échelle, tantôt comme une superposition d'échelles différentes. En effet, $ϕ(2x), ϕ(2x-1),ψ(x) $ représentent l'échelle fine \(1/2\), tandis que \(\varphi(x)\) représente l'échelle grossière \(1\) (un petit dessin peut aider).

Une autre propriété importante est que le support des fonctions \(\varphi\) et \(\psi\) est borné, et donc la fonction \(\varphi_{jk} : x \mapsto \varphi( 2^j x - k )\) à un support \(S_j\) de plus en plus fin lorsque \(j\) augmente. On dit que les ondelettes ont une influence locale.

5.2. Définition

Nous allons construire un système d'ondelettes générique à partir des coefficients :

\[a_k^i \quad : \quad i = 0,1,...,m-1 \qquad k, m \in \setZ, m > 1\]

Adoptant la convention \(\sum_k = \sum_{k \in \setZ}\), nous demandons que ces coefficients vérifient les propriétés de normalisation :

\( \sum_k a^0_{p+m k} = 1 \\ \sum_k a_k^i = m \indicatrice_{i,0} \)

et d'orthogonalisation suivantes :

\( \sum_k a_{k + m r}^i a_{k + m s}^j = m \indicatrice_{ij} \indicatrice_{rs} \\ \sum_k\sum_{i=0}^{m-1} a^i_{r+m k} a^i_{s+m k} = m \indicatrice_{rs} \)

Par ailleurs, afin d'assurer le caractère compact des fonctions associées, on impose la nullité des coefficients en dehors d'un certain domaine :

\( a^i_k = 0 \qquad \forall k \notin \{0,1,2,.... m g - 1 \} \)

où \(g > 0\) est un entier fixé.

Nous allons voir dans la suite que ces propriétés permettent de montrer d'importants résultats sur le comportement des ondelettes associées. Remarquons déjà les conséquences directes suivantes :

\( \sum_k a^0_{p-m k} = 1 \\ \sum_k\sum_{i=0}^{m-1} a^i_{r-m.k} a^i_{s-m.k} = m \indicatrice_{rs} \\ \)

Dans la suite, nous utilisons le produit scalaire :

\[\scalaire{u}{v} = \int_{\setR} u(x) v(x) dx\]

5.3. Construction

La fonction d'échelle s'obtient en partant de :

\[\varphi^{(0)}(x) = \indicatrice_{[0,1[}(x)\]

et en itérant :

\[\varphi^{(n)}(x) = \sum_k a_k^0 \varphi^{(n-1)}(m x - k)\]

Supposant que l'algorithme converge vers un point fixe (voir chapitre \ref{chap:ode}), on définit la fonction d'échelle :

\[\varphi = \lim_{n\rightarrow\infty} \varphi^{(n)}\]

Qui possède la propriété :

\[\varphi(x) = \sum_k a_k^0 \varphi(m x - k)\]

Le système d'ondelettes est définit finalement par :

\( \varphi(x) = \sum_k a_k^0 \varphi(m x - k) \\ \psi^i(x) = \sum_k a_k^i \varphi(m x - k) \)

5.4. Propriétés

Les propriétés des \(a^i_k\) permettent de montrer, par récurrence sur \(n\), que :

\( \int_{\setR} \varphi(x) dx = 1 \\ \int_{\setR} \psi^i(x) dx = 0 \)

En effet, on a \(\int_{\setR} \varphi^{(0)} = 1\) et si on suppose que \(\int_{\setR} \varphi^{(n-1)} = 1\) alors :

\( \int_{\setR} \varphi^{(n)}(x) dx = \sum_k a_k^0 \int_{\setR} \varphi^{(n-1)}(m x - k) dx \\ \int_{\setR} \varphi^{(n)}(x) dx = \left( \sum_k a_k^0 \right) \frac{1}{m} \int_{\setR} \varphi^{(n-1)}(\xi) d\xi = 1 \)

Par ailleurs on a :

\( \int_{\setR} \psi^i(x) dx = \sum_k a_k^i \int_{\setR} \varphi(m x - k) dx \\ \int_{\setR} \psi^i(x) dx = \left( \sum_k a_k^i \right) \frac{1}{m} \int_{\setR} \varphi(\xi) d\xi = 0 \)

On peut aussi montrer par récurrence :

\[\support(\varphi) \subseteq \left[0, 1 + (g-1) \sum_{k=0}^{+\infty} \left( \frac{1}{m} \right)^k \right] = \left[0, 1 + (g-1) \frac{m}{m-1} \right]\]

Le support de \(\varphi\) est donc bien borné.

Une autre propriété importante qui va nous permettre de prouver la convergence des développements en ondelettes est la partition de l'unité :

\[\sum_k \varphi(x-k) = 1 \qquad \forall x \in \setR\]

En effet, il est clair que \(\sum_k \varphi^{(0)}(x-k) = 1\). Partant de l'hypothèse de récurrence \(\sum_k \varphi^{(n-1)}(x-k) = 1\), on obtient :

\( \sum_k \varphi^{(n)}(x-k) = \sum_k\sum_l a^0_l \varphi^{(n-1)}\left( m (x - k) - l \right) \\ \sum_k \varphi^{(n)}(x-k) = \sum_l\sum_k a^0_l \varphi^{(n-1)}(m x - m k - l) \)

Posant \(p = m k + l\), cette dernière équation devient :

\[\sum_k \varphi^{(n)}(x-k) = \sum_p \varphi^{(n-1)}(m x - p) \sum_k a^0_{p - m k}\]

Mais comme $ ∑k a0p - m k = 1$ quel que soit la valeur de \(p\), on en déduit que :

\[\sum_k \varphi^{(n)}(x-k) = \sum_p \varphi^{(n-1)}(m x - p) = 1\]

et la propriété de partition de l'unité est démontréé en passant à la limite :

\[\sum_k \varphi(x-k) = \lim_{n \rightarrow +\infty} \sum_k \varphi^{(n)}(x-k) = 1\]

5.5. Niveaux

On définit la fonction d'échelle et les ondelettes du niveau \(j \in \setN\) comme étant les dilatations translatées des fonctions \(\varphi,\psi^i\) :

\( \varphi_{jk}(x) = m^{j/2} \varphi(m^j x - k) \\ \psi^i_{jk}(x) = m^{j/2} \psi^i(m^j x - k) \)

Le coefficient \(m^{j/2}\) permet d'obtenir les normalisations : \(\scalaire{\varphi_{jk} }{\varphi_{jk} } = 1\) et \(\scalaire{\psi^i_{jk} }{\psi^i_{jk} } = 1\)

La définition du système d'ondelettes permet d'obtenir les relations suivantes entres échelles voisines :

\( \varphi_{jk} = \frac{1}{\sqrt{m}} \sum_l a^0_l \varphi_{j+1,m k + l} = \frac{1}{\sqrt{m}} \sum_p a^0_{p - m k} \varphi_{j+1,p}\\ \psi^i_{jk} = \frac{1}{\sqrt{m}} \sum_l a^i_l \varphi_{j+1,m k + l} = \frac{1}{\sqrt{m}} \sum_p a^i_{p - m k} \varphi_{j+1,p}\\ \)

5.6. Développement en série et algorithme de Mallat

5.6.1. Convergence

Commençons par remarquer que :

\[\sup_{\xi\in\setR} \varphi_{jk}(\xi) = m^{j/2} \sup_{\xi\in\setR} \varphi(\xi) = m^{j/2} M\]

où \(M\) ne dépend pas de \(j\).

Soit la mesure de Lebesgue \(\mu\) définie par :

\[\mu([a,b]) = \abs{b-a}\]

On a alors :

\[\mu(\support(\varphi_{jk})) = m^{-j} \mu(\support(\varphi)) = m^{-j} \indicatrice_0\]

où \(\indicatrice_0\) de dépend pas de \(j\). Le support des fonctions \(\varphi_{jk}\) est donc de plus en plus fin lorsque \(j\) augmente.

Nous définissons la famille d'ensembles :

\[K_j(x) = \{ k \in \setZ : \varphi_{jk}(x) \ne 0 \}\]

où \(x \in \setR\). On obtient, à partir des propriétés du support de \(\varphi\), que :

\[m^j x - k \in \suppport(\varphi) \Rightarrow m^j x - 2 g \le k \le m^j x\]

pour tout \(k\in K_j(x)\). On a donc :

\[\# K_j(x) \le 2 g\]

Pour un réel \(x\) fixé, le nombre de \(\varphi_{jk}(x)\) non nuls est donc borné, et ce quelle que soit la valeur de \(j \in \setN\).

Choisissons \(x,y \in \suppport(\varphi_{jk})\). On a alors :

\[\abs{x-y} \le \indicatrice_0 m^{-j}\]

et par continuité de \(f\) :

\[\abs{f(x)-f(y)} \le \epsilon_j\]

avec :

\[\lim_{j \rightarrow +\infty} \epsilon_j = 0\]

Utilisant un simple changement de variable \(\xi = m^j y - k\), on obtient :

\( \int_{\setR} \varphi_{jk}(y) dy = \int_{\setR} m^{j/2} \varphi\left( m^j y - k \right) dy \\ \int_{\setR} \varphi_{jk}(y) dy = \int_{\setR} m^{j/2} \varphi(\xi) m^{-j} d\xi \\ \int_{\setR} \varphi_{jk}(y) dy = m^{-j/2} \int_{\setR} \varphi(\xi) d\xi = m^{-j/2} \)

Comme \(m^j x \in \setR\), on a aussi :

\[\sum_k \varphi_{jk}(x) = m^{j/2} \sum_k \varphi( m^j x - k ) = m^{j/2}\]

Donc :

\( \sum_k \varphi_{jk}(x) \int_{\setR} \varphi_{jk}(y) dy = m^{-j/2} \sum_k \varphi_{jk}(x) \\ \sum_k \varphi_{jk}(x) \int_{\setR} \varphi_{jk}(y) dy = m^{-j/2} m^{j/2} = 1 \)

On en déduit en particulier que :

\[f(x) = \sum_k \varphi_{jk}(x) \int_{\setR} f(x) \varphi_{jk}(y) dy\]

On aimerait bien montrer que les fonctions échelles de niveau \(j\) quelconque forment un cadre (voir chapitre \ref{chap:tensor}), c'est-à-dire que :

\[f = \sum_k\scalaire{f}{\varphi_{jk}}\varphi_{jk}\]

Pour \(x\) fixé, seuls les \(k\) appartenant à \(K_j(x)\) vont donner une contribution non nulle à la somme :

\[\sum_k\scalaire{f}{\varphi_{jk}}\varphi_{jk}(x)} = \sum_{k \in K_j(x) } \scalaire{f}{\varphi_{jk}}\varphi_{jk}(x)}\]

Rassemblant tous les résultats ci-dessus, on en déduit que :

\( \abs{f(x) - \sum_k\scalaire{f}{\varphi_{jk}}\varphi_{jk}(x)} = \abs{\sum_{k \in K_j}\varphi_{jk}(x) \int_{\setR} \left( f(x) - f(y) \right) \varphi_{jk}(y) dy } \\ \abs{f(x) - \sum_k\scalaire{f}{\varphi_{jk}}\varphi_{jk}(x)} \le 2 g M m^{j/2} \epsilon_j \indicatrice_0 m^{-j} M m^{j/2} = 2 g M^2 \indicatrice_0 \epsilon_j \)

qui tend bien vers 0 lorsque \(j\) tend vers l'infini.

5.6.2. Algorithme de Mallat

Considérons à présent l'approximation suivante :

\[f(x) \approx \tilde{f}(x) = \sum_k c_{Rk} \varphi_{Rk}(x)\]

où \(c_{Rk} = \scalaire{f}{\varphi_{Rk}}\). Nous allons maintenant montrer par récurrence que pour tout \(J \le R\) :

\[\tilde{f}(x) = \sum_k c_{Jk} \varphi_{Jk}(x) + \sum_k \sum_{i=1}^{m-1} \sum_{j=J}^{R} d^i_{jk} \psi^i_{jk}(x)\]

où :

\( c_{jk} = \scalaire{ f }{ \varphi_{jk} } \\ d^i_{jk} = \scalaire{ f }{ \psi^i_{jk} } \)

Il suffit pour cela de constater que :

\( \sum_k c_{jk} \varphi_{jk} + \sum_k\sum_{i=1}^{m-1} d^i_{jk} \psi^i_{jk} = \sum_k \scalaire{ f }{ \varphi_{jk} } \varphi_{jk} + \sum_k\sum_{i=1}^{m-1} \scalaire{ f }{ \psi^i_{jk} } \psi^i_{jk} \\ = \frac{1}{m}\sum_{p,q}\sum_k \sum_{i=0}^{m-1} \scalaire{ f }{ \varphi_{j+1,p} } \varphi_{j+1,q} a^i_{p - m k} a^i_{q - m k} \\ = \sum_{p,q} \scalaire{ f }{ \varphi_{j+1,p} } \varphi_{j+1,q} \frac{1}{m} \sum_k \sum_{i=0}^{m-1} a^i_{p - m k} a^i_{q - m k} \\ = \sum_{p,q} \scalaire{ f }{ \varphi_{j+1,p} } \varphi_{j+1,q} \indicatrice_{pq} \\ = \sum_p \scalaire{ f }{ \varphi_{j+1,p} } \varphi_{j+1,p} \)

en utilisant les propriétés des \(a^i_k\). On peut donc construire le développement du niveau \(j+1\) à partir des fonctions d'échelles et des ondelettes du niveau \(j\). La définition des \(c_{jk}\), \(d^i_{jk}\) combinée aux équations reliant les échelles voisines nous donne :

\( c_{jk} = \scalaire{ f }{ \varphi_{jk} } = \frac{1}{\sqrt{m}}\sum_p\scalaire{ f }{ \varphi_{j+1,p} } a^0_{p - m k} \\ c_{jk} = \frac{1}{\sqrt{m}}\sum_p c_{j+1,p} a^0_{p-m.k} = \frac{1}{\sqrt{m}}\sum_l c_{j+1,l + m k} a^0_{l} \)

ainsi que :

\( d^i_{jk} = \scalaire{ f }{ \psi^i_{jk} } = \frac{1}{\sqrt{m}}\sum_p\scalaire{ f }{ \varphi_{j+1,p} } a^i_{p - m k} \\ d^i_{jk} = \frac{1}{\sqrt{m}}\sum_p c_{j+1,p} a^i_{p - m k} = \frac{1}{\sqrt{m}}\sum_l c_{j+1,l + m k} a^i_{l} \)

La relation inverse s'obtient à partir de :

\( ∑p \scalaire{ f }{ ϕj+1,p } ϕj+1,p = ∑k cjk ϕjk + ∑k∑i=1m-1 dijk ψijk
p cj+1,p ϕj+1,p = \frac{1}{\sqrt{m}}∑p ϕj+1,p \left[ ∑k cjk a0p - m k

  • k∑i=1m-1 dijk aip - m k \right]

\)

En comparant les coefficients des \(\varphi_{j+1,p}\), et comme cette équation doit être valable pour tout \(x\), on a :

\( c_{j+1,p} = \frac{1}{\sqrt{m}}\sum_k c_{jk} a^0_{p-m.k} + \sum_k\sum_{i=1}^{m-1} d^i_{jk} a^i_{p - m k} \\ c_{j+1,p} = \frac{1}{\sqrt{m}}\sum_k c_{jk} a^0_{p-m.k} + \sum_k\sum_{i=1}^{m-1} d^i_{jk} a^i_{p - m k} \\ \)

On voit donc que pour passer d'une représentation mono-échelle

\[f(x) = \sum_k c_{Rk} \varphi_{Rk}(x)\]

à une représentation multi-échelle

\[f(x) = \sum_k c_{Jk} \varphi_{Jk}(x) + \sum_k \sum_{i=1}^{m-1} \sum_{j=J}^{R} d^i_{jk} \psi^i_{jk}(x)\]

(ou l'inverse) il suffit d'appliquer successivement les équations ci-dessus reliant entre-eux les coefficients \(c_{j+1,p}\), \(c_{jk}\), \(d^i_{jk}\), pour \(j = R-1,R-2,...,J\). En pratique, on approxime les coefficients du niveau de résolution maximal \(R\) par :

\[c_{Rk} \approx m^{R/2} f \left(\frac{k}{m^R}\right)\]

Lorsque \(f\) est assez régulière pour avoir la convergence :

\[f(x) = \lim_{R \to +\infty} \sum_k c_{Rk} \varphi_{Rk}(x)\]

on écrit :

\[f(x) = \sum_k c_{Jk} \varphi_{Jk}(x) + \sum_k \sum_{i=1}^{m-1} \sum_{j=J}^{+\infty} d^i_{jk} \psi^i_{jk}(x)\]

Auteur: chimay

Created: 2023-05-10 mer 16:46

Validate