Brevet FR2980604A1 — publié le 20130329

PROCEDE DE TRANSFORMATION RAPIDE EN ONDELETTES D'UN SIGNAL L'invention porte sur un procédé de transformation en ondelettes rapide d'un signal. Un procédé selon l'invention est mis en oeuvre au moyen d'un processeur, de préférence en temps réel. L'invention porte également sur un processeur configuré ou programmé pour la mise en oeuvre d'un tel procédé. Le procédé convient particulièrement bien à la mise en oeuvre d'une transformation continue en ondelettes, mais peut s'appliquer également au cas d'une transformation discrète. Le signal à transformer est généralement un signal numérique ; il peut notamment s'agir d'un signal obtenu par échantillonnage et quantification d'un signal analogique, représentatif d'une grandeur physique. En particulier, le procédé de l'invention peut être appliqué au traitement de signaux neuronaux (neuroélectriques), notamment dans les systèmes de commande neuronale directe, ou d'interface cerveau-ordinateur. La transformation en ondelettes est un outil important en traitement du signal, utilisé en particulier pour réaliser des analyses temps-fréquence. De manière générale, une transformée en ondelettes est définie par la convolution 1 ci?' (t - T Ws(z)= v7s, Lx(t)yf s où : t est le temps ; x(t) est le signal à transformer ; la fonction ip E L'(R) est l' « ondelette mère », dite aussi « fonction génératrice » ; généralement, elle présente un support compact [-a, +a]; - s>0 est appelé le « facteur d'échelle » ; - ti est le décalage temporel ; (1) - « * » indique le complexe conjugué ; en général, en effet, peut être une fonction à valeurs complexes, mais dans la suite on considérera, sans perte de généralité, le cas où elle est une fonction réelle ; - M(r) est un coefficient d'ondelette, fonction de s et de T. On parle de transformation discrète en ondelettes (DVVT, de l'anglais « Discrete Wavelet Transform ») lorsque s=arr' et -e=nb, avec m,n entiers positifs. Si s et T peuvent prendre toute valeur - toute valeur positive, dans le cas de s - on parle de transformation continue en ondelettes (CWT, de l'anglais « Continuous Wavelet Transform »), même si le signal x est discret, ce qui est généralement le cas pour en permettre un traitement numérique. Les DWT sont généralement mieux adaptées, par exemple, au traitement de signaux neuronaux, parce qu'elles permettent d'analyser lesdits signaux à la résolution fréquentielle voulue. Pour plus de précisions sur la transformée en ondelettes on 15 peut se rapporter à : Béatrice PESQUET-POPESCU et Jean-Christophe PESQUET, « Ondelettes et applications », Techniques de l'ingénieur, fascicule TE 5 215. Des algorithmes rapides ont été développés pour réaliser des transformations discrètes en ondelettes. En revanche, les transformations 20 continues en ondelettes restent relativement onéreuses en termes de temps et de ressources, ce qui limite leur utilisation dans les applications temps réel et à faible consommation d'énergie. De manière conventionnelle, une CWT est mise en oeuvre, de façon approchée, en remplaçant, dans l'équation (1), l'intégrale par une somme 25 discrète. On considère le cas d'un signal x(t) échantillonné aux intervalles discrets t=-1, 0, 1, 2 ... (dans un souci de simplicité, et sans perte de généralité, on prend l'échelle des temps telle que son pas de discrétisation vaut 1), et on pose x,=x(ti) ; l'ondelette mère iv est définie sur un support compact [-a, +a]. On définit la fonction :w,(t) (-t , dont le support est [-as, +a si avec as=a/s. 30 L'équation (1) devient alors : as as-1 i+1 Ws(r) = IX(t r)Vs(t)Ch = IX(t r)VsOdt Pour calculer cette somme d'intégrales on effectue l'approximation de signal constant par morceaux : x(t)=x; Vte[i, i+1[ (d'autres approximations sont possibles, par exemple une interpolation linéaire entre t=i et t=i+1) ; en outre, on se limite aux valeurs entiers de ti (T=...,-1, 0, 1, ...). On peut donc écrire : a-1 Ws (r) =-7 Aix (2) 1+1 avec A, = sOdt En considérant connus les coefficients Ai, qui peuvent être calculés une fois pour toutes, on voit que la détermination d'un seul coefficient d'ondelette nécessite de (2a-1) multiplications et (2a-2) sommes. Or, a peut prendre des valeurs relativement importantes ; par exemple, pour effectuer l'analyse temps-fréquence d'un signal échantillonné à 1 kHz dans la plage fréquentielle 10 Hz - 300 Hz avec des ondelettes de Meyer, a est typiquement compris entre 100 et 3000. L'article de M. Omachi et S. Omachi « Fast calculation of continuous wavelet transform using polynomial », Proceedings of the 2007 international Conference on Wavelet Analysis and Pattern Recognition, Beijing, Chine, 2 - 4 novembre 2007, décrit un procédé de CWT rapide d'un signal, 2 0 dans lequel le nombre d'opérations requises ne dépend pas de a. Dans ce procédé, l'ondelette mère est définie par un polynôme d'ordre d sur un support compact ; typiquement, l'ordre d et les coefficients du polynôme sont choisis de manière à approcher, avec l'exactitude recherchée, une ondelette « de référence » non polynômiale, telle qu'une ondelette de Haar. Les coefficients 2 5 d'ondelettes sont déterminés de manière récursive ; ainsi, M(c) est obtenu à partir de VVs_i(t), le facteur d'échelle étant considéré entier. Le coût computationnel de la détermination du coefficient Ws(T) est O(ds). Ce procédé présente un certain nombre d'inconvénients. Premièrement, l'approximation polynomiale d'une ondelette « de référence » telle que l'ondelette de Haar peut nécessiter le recours à un ordre d élevé, ce qui augmente d'autant le coût computationnel du procédé. Deuxièmement, il ne permet pas de déterminer les coefficients Ws pour différents valeurs de s - et donc pour différentes fréquences indépendamment les uns des autres. Au contraire, en raison du caractère récursif de l'algorithme, la détermination de Ws implique nécessairement celle de W1 ... Le procédé de M. Omachi et S. Omachi n'est donc pas optimisé pour les applications dans lesquelles un signal doit être analysé seulement dans certaines bandes spectrales. L'invention vise à procurer un procédé permettant de déterminer rapidement au moins un coefficient d'ondelette Ws(T) d'une transformée en ondelettes d'un signal, ne présentant pas les inconvénients précités de l'art antérieur. Un objet de l'invention permettant d'atteindre ce but est un procédé de détermination d'au moins un coefficient d'ondelette Ws(t) d'une transformation en ondelettes d'un signal, ladite transformation en ondelettes étant définie par rapport à une ondelette mère, caractérisé en ce que : - ladite ondelette mère est une fonction continue, présentant un support subdivisé en ..b1 intervalles délimités par (J+1) points d'extrémité, ladite fonction étant définie par un polynôme de degré maximal Nbl sur chaque dit intervalle ; et en ce que ledit procédé comprend les étapes consistant à: a) calculer toutes ou une partie des primitives d'ordre k compris entre 2 et N+1 dudit signal, au moins en (J+1) points correspondant auxdits points d'extrémité des intervalles du support de rondelette dilaté d'un facteur s et translaté d'un temps 'G ; b) calculer la convolution de ladite ou chaque dite primitive ainsi échantillonnée avec une succession respective de (J-1-1) coefficients C,k(s), fonctions de ladite ondelette ; et c) déterminer ledit coefficient d'ondelette par calcul d'une combinaison linéaire desdites convolutions ; ledites étapes a) à c) étant mises en oeuvre par un processeur configuré ou programmé de manière appropriée. Contrairement au cas de l'algorithme de M. Omachi et S. Omachi, dans le procédé de l'invention l'ondelette mère n'est pas polynomiale, mais polynomiale par morceaux. Ainsi, l'ordre maximal N peut être maintenu petit (typiquement, N=2 ou N=3) tout en permettant d'approcher efficacement toute ondelette mère de référence à support compact. Le coût computationnel dépend à la fois de N (ordre maximal des polynômes définissant l'ondelette mère), de J (nombre d'intervalles ou morceaux) et de la régularité de l'ondelette yr, polynomiale par morceaux ; on dispose ainsi d'une pluralité de paramètres permettant d'optimiser le procédé. En outre, les différents coefficients 1/1/,(T) peuvent être déterminés indépendamment les uns des autres. Selon différents modes de réalisation du procédé de l'invention : - Lesdits coefficients C,k(s) peuvent être des fonctions des discontinuités des dérivés d'ordres compris entre 1 et N de l'ondelette en correspondance desdits points d'extrémités des intervalles constituant son support. - Les étapes b) et c) peuvent être mises en oeuvre en appliquant l'équation suivante : N+1 J (-1)/z Jc,k(s)x,k)(z+-i-.)s Ws(i-) E k=2 1=0 où : - vvs(c) est le coefficient d'ondelette recherché, correspondant à un décalage T et à une dilatation s d'une ondelette mère y(t) X(k)(t) est la primitive d'ordre k du signal ; - oc, sont les points d'extrémité des intervalles constituant le support de l'ondelette mère, pris dans l'ordre ; et - C,k(s), avec i compris entre 0 et J, sont les coefficients formant la succession avec laquelle est convolutée la primitive d'ordre k du signal. En particulier, les coefficients C,k(s) peuvent être donnés par : e,` (s) = s () j(k--" (c/), i -1 Cô (s) = -(Vi)(k-1) (c) C (s) (1 6r) (°-') où (yis)(k-1) est la dérivée d'ordre (k-1) de la fonction yis= w;'(t) - , yi étant le polynôme qui définit l'ondelette mère sur le i-ème intervalle de son support. - Ladite étape a) peut comprendre : al) le calcul desdites primitives d'ordre k compris entre 2 et N+1 dudit signal ; et a2) leur échantillonnage en (J+1) points d'échantillonnage correspondant auxdits points d'extrémité des intervalles du support de l'ondelette dilaté d'un facteur s et translaté d'un temps T. - Dans un cas particulièrement simple à mettre en oeuvre, ladite ondelette peut être une fonction continue et à dérivée première continue, polynomiale d'ordre deux par morceaux ; dans ce cas, ladite étape a) comporte le calcul de la seule primitive d'ordre 3 dudit signal et les étapes b) et c) peuvent être mises en oeuvre en appliquant l'équation suivante : W (r) 113 X{3) (i" C) ;=a où: Bi=2 1 1 (p. 1- p.),i =1,...,J -1 -4 82 B1 -1p0 - 0 s2 5 1 1 Bj = 22 p s Ladite étape a) peut être mise en oeuvre en temps réel, lesdites primitives étant alors calculées en tant que sommes cumulatives pendant une étape d'acquisition dudit signal, puis stockées dans une première mémoire dudit processeur. En outre, lesdites successions de coefficients Cik(s) peuvent être stockées dans une deuxième mémoire dudit processeur préalablement à ladite étape d'acquisition dudit signal. Un autre objet de l'invention est un procédé de transformation en ondelettes d'un signal à temps discret ou discrétisé, comportant la détermination d'une pluralité de coefficients d'ondelette associés à différentes dilatations/translations d'une même ondelette mère, caractérisé en ce que chaque dit coefficient d'ondelette est déterminé par un procédé tel que décrit ci-dessus. Encore un autre objet de l'invention est un procédé de transformation en ondelettes d'un signal, comportant la détermination d'une pluralité de coefficients d'ondelette associés à différentes dilatations/translations d'une même ondelette mère, caractérisé en ce que chaque coefficient d'ondelette VV,(T), correspondant à un décalage T et à une dilatation s d'une ondelette mère w(t), est déterminé par un procédé tel que décrit ci-dessus, dans lequel : - ladite étape a) est mise en oeuvre en temps réel, lesdites primitives étant calculées en tant que sommes cumulatives pendant une étape d'acquisition dudit signal, puis stockées dans une première mémoire dudit processeur ; - lesdites étapes b) et c) sont mises en oeuvre en temps réel en utilisant des successions de coefficients préalablement stockées dans une deuxième mémoire dudit processeur ; les coefficients d'ondelette VVs(t), correspondant à un décalage et à une dilatation s d'une ondelette mère w(t), étant déterminés à partir du temps t=t-Fa/s, où a est la demi-largeur du support de ladite ondelette mère. Ledit signal peut être un signal numérique, le cas échéant obtenu par échantillonnage et quantification d'un signal analogique. En particulier, le signal peut un signal acquis par un détecteur tel une électrode, un photodétecteur, un détecteur acoustique ou, d'une manière générale, tout détecteur ou capteur délivrant un signal électrique sous l'effet d'un paramètre physique, chimique ou biologique. Encore plus particulièrement, il peut s'agir d'un signal représentatif d'une activité neuronale d'un cerveau, notamment humain (signal neurologique ou neuroélectrique). Encore un autre objet de l'invention est un processeur configuré ou programmé pour mettre en oeuvre un procédé selon l'une des revendications précédentes. On entend par processeur toute machine pour le traitement de l'information ; typiquement, il pourra s'agir d'un processeur électronique, de préférence numérique. En particulier il pourra s'agir d'une ordinateur conventionnel programmé de manière opportune, d'un processeur dédié au traitement de signaux numériques (DSP, de l'anglais « Digital Signal Processor »), d'un circuit numérique câblé, d'un système hybride programmé/câblé, etc. D'autres caractéristiques, détails et avantages de l'invention ressortiront à la lecture de la description faite en référence aux dessins annexés donnés à titre d'exemple et qui représentent, respectivement : - la figure 1, un schéma à blocs d'un procédé selon un mode de réalisation de l'invention ; et - la figure 2, un schéma fonctionnel d'un processeur pour la mise en oeuvre d'un procédé en temps réel selon un mode de réalisation de l'invention. Le procédé selon l'invention sera d'abord illustré dans le cas simple d'une ondelette mère quadratique par morceaux (N=2) définie sur l'intervalle [-a, a], et pour s=1. Dans ce cas, on peut écrire : W(r) Jx(t I") (t)dt (3) -a Soit -a = ao La dérivé seconde de Yi(t) est une fonction constante par morceaux : PO= 2Pj V t E [a» a j+], ce qui permet d'écrire : W (r) J -le J-,] J -1 a.,+] J -1a' +i X (2) (t r)q " (t)dt =1219 JX(2)(t+r)dt=E2p1 X(2) (t + r)dt = .1-0 cc, J=0 ce, J=0 ce, J -1 j=1 J E 2p .(X(3) (r + aj+1)- X(3) (r + a .)) 2p i_LX(3)(r + a j) _E 2p .X(3) (y + a _) j-O j=-0 J -1 = 2pj iX(3)(r + a j) - 2p0X(3)(r+ cx0)+ E 2(p j_t - pdX(3) j=1 où X(3)(t) est une primitive d'ordre 3 de la fonction x(t), c'est-à-dire une fonction dont la dérivée troisième par rapport à t est X(t) . On obtient donc : (8) W(r)=1B,X(3)(r +a,) i=o avec: B1 =2(p,_ - pi),i =1 ...,J -1 ; Bo = -2p0 ; Bi = 2pj En supposant que la primitive d'ordre 3 de x(t) est connue, il faut donc J+1 multiplications et autant de sommes pour déterminer un coefficient d'ondelette, au lieu des (2a-1) multiplications et (2a-2) sommes nécessaires pour appliquer directement l'équation (2). 11 Si la dérivée première de rondelette mère Kt) n'est pas J-1 continue, alors la somme Ex(2)(t+/-)v ai+ta l n'est plus égale à zéro. j=0 Un calcul élémentaire permet de montrer que J-] J -1 X(2) (I ±z)Vf(tla1+1 =1(X(2) (t r)(2p it +q i)r.1+1 a i=0 J J=a J-1 -= (2a i(p p1) +q i)X(2) (a j + r)+ - (2p0a0 +q 0)X (2) (a0 + r) + (2/9 a +q j_i)X(2) (a j +r). Le coefficient d'ondelette W(T) est donc donné par : W(z) =BiX(3)(z + a i) + (-1)1C,X(2) (a + r) (9) 10 avec: Bi = 2(pi_1 B0 - 2p0 B., 2p ; Ci = 2ai p1) +q 1-1-q Co = -(2p0a0 +q0); Cj=2pj ,aj+q,_, - =1,...,J -1; 15 Dans ce cas, la détermination d'un coefficient d'ondelette nécessite 2(J+1) sommes et multiplications. L'équation (9) peut être généralisée au cas d'un ordre de polynôme N quelconque. Dans ce cas rondelette mère te) peut être exprimée wi. (t) , t E [a i , a i.+1] { 0, ltk a par : 1/1(t) N N N-1 N-1 0 avec tif (t) = p. t -Fpf t w(t) étant en outre une fonction continue et telle que K-a) = te(ceo) = j) = Ka) = 0 . En intégrant (3) par parties (N4-1) fois on obtient : (10) N+1 J W(r) (_1)k-1 cf x(k)(a1) Ic=2 i=0 avec : f i(I ci-1) (a iir{k-1) =1,...,J -1, Cok = (a0) C11.; = jk-1) (a,) L'équation (10) montre qu'en général le nombre d'opération nécessaires à la détermination d'un coefficient d'ondelette est proportionnel à l'ordre N des polynômes. Toutefois, cela n'est pas vrai si v(t) est choisie de telle manière que ses dérivées soient continues, parce que dans ce cas plusieurs des coefficients Cik s'annulent. Dans le cas d'une fonction polynômiale par morceaux d'ordre N, continue et dont toutes les dérivées jusqu'à l'ordre (N-1) sont continues, l'équation (10) se simplifie de la façon suivante : J W(r)= (-oNE Ci(N+1) X(N+1) (r+ ai) i=0 - = - piN ),z J 1 (N+1) C0 - 0 (N+I) C = N! pj Dans ce cas, le coût computationnel ne dépend pas de l'ordre N, et seule la primitive d'ordre (N+1) doit être calculée. Par conséquent, pour minimiser le coût computationnel il convient d'utiliser des polynômes d'ordre faible et/ou des ondelettes mères présentant des dérivées continues. Jusqu'ici on s'est limités au cas d'un facteur d'échelle s constant, conventionnellement pris égal à 1. Il est désormais possible de (N+1) relâcher cette hypothèse et de considérer une ondelette ws(t) = 1 r t , s a a continue, définie sur l'intervalle S S avec la partition aa ai a a = 0 Dans le cas particulier d'une ondelette mère constituée par des polynômes d'ordre deux, VI (t) = p1t2 +q it + r; et présentant une dérivée première continue on obtient : J (13) W (r) EB,X(3)(r + (2) 1=0 avec: Bi = 2rs2 1 -1 (p p, =1,..., J -1 ; 1 1 B = - NS ° s- ' 1 1 n B = p ''1,7 S2 j - 10 Différentes approches peuvent être suivies pour calculer les primitives X(1), x(k), x(N+1). En général, le signal x(t) est en réalité un signal discret, ou bien un signal continu échantillonné à des instants discrets : x(t)x(i)=xi, i=...-1, 0, 1, 2, ... Pour calculer les primitives, on peut notamment se baser sur 15 l'approximation de signal constant par morceaux : x(t) = x, , t e [i, i +11. On remarque que i est l'arrondi à l'entier inférieur du temps t (avec une échelle des temps appropriée) t [d. On peut donc écrire t = i + A, A=t-i La figure 2 illustre l'application d'un tel procédé à la transformation en ondelettes en temps réel d'un signal x(t), c'est-à-dire à la détermination des coefficients M(c) pour une pluralité de valeurs de s=si, s2 et de T=Ti, 12 - - Le bloc UCP calcule au fil de l'eau les primitives X(k)(t) comme sommes cumulatives des échantillons du signal x(t). Les primitives ainsi calculées sont stockées dans la mémoire MP, qui est de préférence une mémoire glissante. En même temps, les valeurs déjà calculées des primitives ai correspondant aux temps 1-ssont lues de ladite mémoire MP et adressées vers une unité de calcul UCCL. Les coefficients C,k(s) sont également lus d'une deuxième mémoire MC et adressés à leur tour vers l'unité UCCL. Cette dernière détermine les coefficients d'ondelette VVs(c) dès que cela est possible. En effet, l'équation (12) montre que pour déterminer W s(r) il est nécessaire de connaitre les primitives X(k)seulement dans l'intervalle [ 1- + -(1 . Par s conséquent, à chaque instant t, il est possible de déterminer Ws(r)V La transformation en ondelettes peut donc être effectuée en temps réel, avec un décalage égal à a+Tp par rapport au signal, Tp étant un retard dû au temps de traitement. Là encore, les blocs UCP, MP, MC, UCCL doivent être considérés comme des blocs fonctionnels, qui ne correspondent pas nécessairement à des structures physiques d'un processeur utilisé pour la mise 10 en oeuvre du procédé. REVENDICATIONS1. Procédé de détermination d'au moins un coefficient d'ondelette VVs(t) d'une transformée en ondelettes d'un signal, ladite 5 transformée en ondelettes étant définie par rapport à une ondelette mère, caractérisé en ce que : - ladite ondelette mère est une fonction continue, présentant un support subdivisé en 'b1 intervalles délimités par (J+1) points d'extrémité, ladite fonction étant définie par un polynôme de degré maximal sur chaque 10 dit intervalle ; - et en ce que ledit procédé comprend les étapes consistant à: a) calculer toutes ou une partie des primitives d'ordre k compris entre 2 et N+1 dudit signal, au moins en (J+1) points correspondant 15 auxdits points d'extrémité des intervalles du support de l'ondelette dilaté d'un facteur s et translaté d'un temps z ; b) calculer la convolution de ladite ou chaque dite primitive ainsi échantillonnée avec une succession respective de (J+1) coefficients C,k(s), fonctions de ladite ondelette et 20 c) déterminer ledit coefficient d'ondelette par calcul d'une combinaison linéaire desdites convolutions ; ledites étapes a) à c) étant mises en oeuvre par un processeur configuré ou programmé de manière appropriée. 25 2. Procédé selon la revendication 1 dans lequel lesdits coefficients C,k(s) sont des fonctions des discontinuités des dérivés d'ordres compris entre 1 et N de l'ondelette en correspondance desdits points d'extrémités des intervalles constituant son support. 30 3 Procédé selon la revendication 2 dans lequel les étapes b) et c) sont mises en oeuvre en appliquant l'équation suivante : N+1 J (r) cric (s)x(k)i_ , r=0 " Soù : - vvs(`c) est le coefficient d'ondelette recherché, correspondant à un décalage T et à une dilatation s d'une ondelette mère y(t) - X(k)(t) est la primitive d'ordre k du signal ; - ai sont les points d'extrémité des intervalles constituant le support de rondelette mère, pris dans l'ordre ; et - Cik(s), avec i compris entre 0 et J, sont les coefficients formant la succession avec laquelle est convolutée la primitive d'ordre k du signal. 4. Procédé selon la revendication 3 dans lequel les coefficients C,k(s) sont donnés par : cik (s) (if;ç $1)(k 1) (a, ) (vf: )(k-1)(a, ), i J - I cok(s)= -(K,')"-`)(c) C.; (s),()(i)(k-" (a') où (Ivis)(1 V; = i(-t , yi étant le polynôme qui définit rondelette mère sur le i-ème s intervalle de son support. 5. Procédé selon l'une des revendications précédentes, dans laquelle ladite ondelette est une fonction continue et à dérivée première continue, polynomiale d'ordre deux par morceaux, dans lequel ladite étape a) comporte le calcul de la seule primitive d'ordre 3 dudit signal et dans lequel les étapes b) et c) sont mises en oeuvre en appliquant l'équation suivante : W (r) =113,X (1- + e)e où: J -11 1 B =-2--Ns' s2-po ° 1 1 = -2p J2 - s 6. Procédé selon l'une des revendications précédentes, dans lequel ladite étape a) comprend : al) le calcul desdites primitives d'ordre k compris entre 2 et N+1 dudit signal ; et a2) leur échantillonnage en (J+1) points d'échantillonnage correspondant auxdits points d'extrémité des intervalles du support de l'ondelette dilaté d'un facteur s et translaté d'un temps T. 7. Procédé selon l'une des revendications précédentes, dans lequel ladite étape a) est mise en oeuvre en temps réel, lesdites primitives étant calculées en tant que sommes cumulatives pendant une étape d'acquisition dudit signal, puis stockées dans une première mémoire (MP) dudit processeur. 8 Procédé selon la revendication 7 dans lequel lesdites successions de coefficients Cik(s) sont stockées dans une deuxième mémoire (MC) dudit processeur préalablement à ladite étape d'acquisition dudit signal. 9. Procédé de transformation en ondelettes d'un signal à temps discret ou discrétisé, comportant la détermination d'une pluralité de coefficients d'ondelette associés à différentes dilatations/translations d'une même ondelette mère, caractérisé en ce que chaque dit coefficient d'ondelette est déterminé par un procédé selon l'une des revendications précédentes. 10. Procédé de transformation en ondelettes d'un signal, comportant la détermination d'une pluralité de coefficients d'ondelette associés à différentes dilatations/translations d'une même ondelette mère, caractérisé en ce que chaque coefficient d'ondeletteWs(T), correspondant à un décalage T et àune dilatation s d'une ondelette mère w(t), est déterminé par un procédé selon l'une des revendications 1 à 8, dans lequel : - ladite étape a) est mise en oeuvre en temps réel, lesdites primitives étant calculées en tant que sommes cumulatives pendant une étape d'acquisition dudit signal, puis stockées dans une première mémoire (MP) dudit processeur ; lesdites étapes b) et c) sont mises en oeuvre en temps réel en utilisant des successions de coefficients préalablement stockées dans une deuxième mémoire (MC) dudit processeur ; les coefficients d'ondelette Ws(T), correspondant à un décalage et à une dilatation s d'une ondelette mère y(t), étant déterminés à partir du temps t=t+cds, où a est la demi-largeur du support de ladite ondelette mère. 11. Procédé selon l'une des revendications précédentes, dans lequel ledit signal est un signal numérique. 12. Procédé selon l'une des revendications précédentes dans lequel le signal est un signal acquis par un détecteur. 13. Processeur configuré ou programmé pour mettre en oeuvre un procédé selon l'une des revendications précédentes.