FR 2853075 A1 20041001 FR 0303545 A 20030324 Procédé pour déterminer des réponses impulsionnelles d'un milieu vis-à- vis de la transmission d'ondes entre différents points. La présente invention est relative aux procédés pour déterminer des réponses impulsionnelles d'un milieu vis-à-vis de la transmission d'ondes entre différents points. Plus particulièrement, l'invention concerne un 10 procédé pour déterminer des réponses impulsionnelles d'un milieu vis-à-vis de la transmission d'ondes entre différents points, procédé comprenant: (a) au moins une étape d'émission au cours de laquelle on émet des ondes dans le milieu en générant des 15 signaux ei(t) à partir d'un nombre N de points d'émission appartenant au milieu, o N est un nombre entier au moins égal à 2 et i est un indice compris entre 1 et N qui désigne l'un desdits N points d'émission, (b) au moins une étape de réception au cours de 20 laquelle on capte des signaux rj(t) à partir desdites ondes après transmission dans ledit milieu, en un nombre M de points de réception appartenant au milieu (ces points de réception peuvent être confondus avec les points d'émission), o M est un entier naturel non nul et j est un 25 indice compris entre 1 et M qui désigne l'un desdits M points de réception, (c) et au moins une étape de détermination desdites réponses impulsionnelles hij(t) entre chaque point d'émission i et chaque point de réception j à partir des 30 signaux émis ei(t) et captés rj(t). Les réponses impulsionnelles peuvent être déterminées aux fins de focalisation d'ondes dans le milieu, par exemple aux fins d'imagerie du milieu ou de communication sélective entre différents points du milieu, 35 notamment dans le cadre de la mise en oeuvre de la méthode "D.O.R.T." (Décomposition de l'opérateur Retournement Temporel) décrite notamment par Prada et al. (C. Prada and M. Fink, "Eigenmodes of the time reversai operator: a solution to sélective focusing in multiple5 target media." Wave Motion, 20, pp 151-163 (1994)). Le processus d'acquisition des réponses impulsionnelles entre les points en question du milieu est toutefois relativement lent, dans la mesure o il est effectué séquentiellement, pour les différents points 10 d'émission i les uns à la suite des autres, comme décrit par exemple dans le document WO-A-02/32 316 dans le cadre d'un procédé de focalisation d'ondes acoustiques autre que la méthode "D.O.R.T." susmentionnée. Cette lenteur allonge le temps de mise en oeuvre du 15 procédé, et lorsque le milieu est changeant, elle ne permet pas d'obtenir les réponses impulsionnelles assez vite pour pouvoir être valablement utilisées ensuite dans ce milieu. La présente invention a notamment pour but de pallier ces inconvénients. A cet effet, selon l'invention, un procédé du genre en question est caractérisé en ce qu'au cours de l'étape (a), on fait émettre simultanément les signaux ei(t) par lesdits N points d'émission, ces signaux ei(t) ayant une durée T et étant chacun une somme de n signaux élémentaires 25 sensiblement monochromatiques, de même amplitude et de fréquences respectives fo,j+k.5f, o fo,j est une fréquence prédéterminée propre au point i, k est un nombre entier compris entre 0 et n, n est un nombre entier au moins égal à 2 et ôf est un pas fréquentiel prédéterminé, les 30 fréquences foi, respectives propres aux différents points i étant distinctes les unes des autres et comprises dans une bande de fréquences de largeur of, et en ce qu'au cours de l'étape (c), on calcule chaque réponse impulsionnelle hij(t) à partir d'un signal de 35 corrélation entre le signal ei(t) émis au point i et le signal rj(t) capté au point j. Grâce à ces dispositions, on calcule simultanément les réponses impulsionnelles entre N points d'émission et M points de réception (éventuellement confondus avec les 5 points d'émission), ce qui accélère le processus d'acquisition des réponses impulsionnelles. Dans divers modes de réalisation du procédé selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes: - les fréquences f0, respectives propres aux différents points i sont séparées deux à deux par un écart 5f/N; - au cours de l'étape (c), ledit signal de corrélation est fenêtré au moyen d'une fonction porte n(t) 15 de largeur 1/5f; - au cours de l'étape (c), les réponses impulsionnelles hij(t) sont déterminées par la formule: hUij(t) = l(t).fei(O-t).rj(O)dO; - les ondes transmises dans le milieu entre les 20 points d'émission et les points de réception sont des ondes acoustiques; - au cours de l'étape (a), le milieu o sont émises les ondes est réverbérant; - le pas fréquentiel 5f est inférieur ou égal à 25 1/l, o -T est la dispersion temporelle du milieu (c'est-àdire la durée du signal capté après émission d'un signal impulsionnel dans le milieu) ; le pas fréquentiel 5f est sensiblement égal à 1/l, o X est la dispersion temporelle du milieu; 30 - la durée T est au moins égale à N/5f; - la durée T est au moins égale à N.T, o T est la dispersion temporelle du milieu; - les signaux élémentaires présentent des phases aléatoires; - les ondes sont émises avec une certaine bande passante, les fréquences fOi comprennent une fréquence minimale fO et le nombre n est déterminé pour que la bande de fréquence comprise entre f0 et fO+[ (n+I) .5f] recouvre sensiblement ladite bande passante; - les points de réception sont confondus avec les points d'émission. D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description 10 suivante de deux de ses modes de réalisation, donnés à titre d'exemples non limitatifs, en regard des dessins joints. Sur le dessin, la figure 1 est une vue schématique d'un dispositif adapté pour mettre en oeuvre la présente 15 invention, dans un mode de réalisation. Dans l'exemple non limitatif représenté sur le dessin, le procédé selon un mode de réalisation de l'invention est mis en òuvre par un dispositif 1 qui est adapté pour émettre et recevoir des ondes dans un milieu 2 20 homogène ou non, généralement réverbérant, qui peut être solide, liquide ou gazeux. Les ondes en question sont de préférence acoustiques et la description ci-après ne fera référence qu'à de telles ondes acoustiques, mais lesdites ondes pourraient éventuellement être d'une autre nature, 25 par exemple de nature électromagnétique. Le dispositif 1 comporte par exemple un ordinateur 3, ou tout autre dispositif de commande similaire, qui commande les émissions d'ondes acoustiques et leur écoute dans le milieu 2. A cet effet, l'ordinateur 3 peut par 30 exemple commander un dispositif électronique dédié 4 qui peut notamment comporter une unité centrale électronique CPU associée à une mémoire centrale MEM, cette unité centrale CPU commandant elle-même N voies (N étant un entier au moins égal à 2) reliées chacune à un transducteur 35 Tl-TN permettant d'émettre et de capter des signaux acoustiques dans le milieu 2. Toutes les voies sont synchronisées entre elles. Chaque voie de mesure du dispositif 4 peut comporter par exemple un échantillonneur El-EN relié au 5 transducteur correspondant Tl-TN, et chaque échantillonneur peut être lui-même relié à une mémoire associée à une unité centrale, référencées conjointement Cl-CN. Le dispositif 1 peut être utilisé notamment aux fins d'imagerie du milieu 2, ou encore pour établir une 10 communication sélective avec un point dudit milieu, ou autre, par exemple en utilisant la méthode "D.O.R.T." (Décomposition de l'opérateur Retournement Temporel) décrite notamment par Prada et al. (C. Prada and M. Fink, "Eigenmodes of the time reversal operator: a solution to 15 sélective focusing in multiple-target media." Wave Motion, 20, pp 151-163 (1994)). Dans tous les cas, il est nécessaire de déterminer les réponses impulsionnelles entre les transducteurs T1-TN. A cet effet, l'ordinateur 3 commande le dispositif 20 1 de façon qu'il suive successivement les étapes suivantes: (a) au moins une étape d'émission au cours de laquelle on émet des ondes dans le milieu en générant des signaux ei(t) à partir des N transducteurs Tl-TN (ces N 25 transducteurs peuvent le cas échéant ne constituer qu'une partie des transducteurs reliés au dispositif 4, auquel cas les réponses impulsionnelles relatives à l'ensemble des transducteurs seraient acquises en plusieurs passes successives, ces passes restant toutefois en nombre très 30 inférieur au nombre total de transducteurs), i étant un indice compris entre 1 et N qui désigne l'un desdits N transducteurs; (b) au moins une étape de réception au cours de laquelle les N transducteurs Tl-TN captent des signaux rj (t) à partir desdites ondes après transmission dans ledit milieu, o j est un indice compris entre 1 et N désignant l'un desdits transducteurs; (c) et au moins une étape de détermination desdites réponses impulsionnelles hij(t) entre chaque point 5 d'émission i et chaque point de réception j à partir des signaux émis ei(t) et captés rj(t). Pour permettre d'acquérir en une seule fois tous les signaux nécessaires à la détermination des réponses impulsionnelles concernant les N transducteurs susmentionnés (c'est-à-dire tout ou partie des transducteurs reliés au dispositif 4, comme expliqué cidessus), tous les signaux ei(t) sont émis simultanément au cours de l'étape (a) par les N transducteurs, mais ces signaux émis sont orthogonaux entre eux de façon que les 15 informations propres à chaque point d'émission puissent ensuite être séparées les unes des autres dans les signaux captés rj(t). Ainsi, chaque signal émis ei(t) est égal à la somme de n signaux élémentaires sensiblement monochromatiques, de 20 même amplitude et de fréquences respectives fo,j+k.5f, o : - fo,j est une fréquence prédéterminée propre au transducteur i, - k est un nombre entier compris entre 0 et n, - n est un nombre entier au moins égal à 2, - et 5f est un pas fréquentiel prédéterminé, avantageusement inférieur ou égal à 1/i et de préférence sensiblement égal à 1/-, o i est une valeur moyenne de la dispersion temporelle du milieu 2, c'est-à-dire la durée du signal capté après émission d'un signal impulsionnel 30 (fonction de Dirac) dans le milieu 2. Les fréquences f0,i respectives propres aux différents points i sont distinctes les unes des autres et comprises dans une bande de fréquences de largeur 5f: avantageusement, ces fréquences fo,j sont séparées deux à 35 deux par un écart 5f/N. De plus, si BP est la bande passante avec laquelle les ondes acoustiques sont émises et reçues par les transducteurs, et si on nomme f0 la plus faible des fréquences fOi, alors f0 et n peuvent avantageusement être 5 déterminés pour que la bande des fréquences des signaux élémentaires (bande comprise entre f0 et f0+ [ (n+l) . f) recouvre sensiblement ladite bande passante. Avantageusement, les signaux élémentaires composant les signaux ei(t) présentent des phases aléatoires. Tous les signaux ei(t) ont une même durée T, qui est avantageusement au moins égale à N/5f, et de préférence au moins égale à N.i. Au cours de l'étape (c), on calcule chaque réponse impulsionnelle hij (t) à partir d'un signal de corrélation 15 entre le signal ei(t) émis au point i et le signal rj(t) capté au point j. Dans l'exemple considéré ici, ce signal de corrélation vaut: cYi(t)= Jei(O).rj(t-O)dO (cette intégrale est faite par exemple sur une durée au moins égale à T, dont le début coïncide avec le début du signal rj(t)). Avantageusement, ce signal de corrélation est fenêtré au moyen d'une fonction porte n(t) de largeur 1/5f (valant 1 pendant une période de 1/5f et 0 en dehors de 25 cette période), auquel cas les réponses impulsionnelles hij(t) peuvent avantageusement être déterminées par la formule: hij(t) = 11(t).cij(t) On notera que, dans l'exemple particulier 30 représenté sur le dessin et décrit ci-dessus, les points de réception j auxquels on capte l'onde acoustique sont confondus avec les points d'émission et constitués par les transducteurs Tl-TN, mais les points de réception en question peuvent le cas échéant être constitués par un second réseau de transducteurs (non représenté), distinct du réseau de transducteurs T1-TN. Dans le cas général, les ondes acoustiques qui sont générées par l'émission des signaux ei(t) sont donc captées 5 en M points de réception, M étant un nombre entier non nul, lesquels M points de réception peuvent être soit confondus avec les N points d'émission (auquel cas M = N), soit non totalement confondus avec les points d'émission, soit entièrement distincts des points d'émission. REVENDICATIONS 1. Procédé pour déterminer des réponses impulsionnelles d'un milieu (2) vis-à-vis de la 5 transmission d'ondes entre différents points (Tl-TN), procédé comprenant: (a) au moins une étape d'émission au cours de laquelle on émet des ondes dans le milieu (2) en générant des signaux ei(t) à partir d'un nombre N de points 10 d'émission (Tl-TN) appartenant au milieu, o N est un nombre entier au moins égal à 2 et i est un indice compris entre 1 et N qui désigne l'un desdits N points d'émission, (b) au moins une étape de réception au cours de laquelle on capte des signaux rj (t) à partir desdites ondes 15 après transmission dans ledit milieu, en un nombre M de points de réception (T1-TN) appartenant au milieu, o M est un entier naturel non nul et j est un indice compris entre 1 et M qui désigne l'un desdits M points de réception, (c) et au moins une étape de détermination 20 desdites réponses impulsionnelles hij (t) entre chaque point d'émission i et chaque point de réception j à partir des signaux émis ei(t) et captés rj(t), caractérisé en ce qu'au cours de l'étape (a), on fait émettre simultanément les signaux ei(t) par lesdits N 25 points d'émission (Tl-TN) , ces signaux ei(t) ayant une durée T et étant chacun une somme de n signaux élémentaires sensiblement monochromatiques, de même amplitude et de fréquences respectives fo,I+k.5f, o fo,j est une fréquence prédéterminée propre au point i, k est un nombre entier 30 compris entre 0 et n, n est un nombre entier au moins égal à 2 et 5f est un pas fréquentiel prédéterminé, les fréquences fo,j respectives propres aux différents points i étant distinctes et comprises dans une bande de fréquences de largeur 5f, et en ce qu'au cours de l'étape (c), on calcule chaque réponse impulsionnelle hij(t) à partir d'un signal de corrélation entre le signal ei(t) émis au point i et le signal rj(t) capté au point j. 2. Procédé selon la revendication 1, dans lequel 5 les fréquences fo,j respectives propres aux différents points i sont séparées deux à deux par un écart 5f/N. 3. Procédé selon la revendication 1 ou la revendication 2, dans lequel au cours de l'étape (c), ledit signal de corrélation est fenêtré au moyen d'une 10 fonction porte n(t) de largeur 1/Sf. 4. Procédé selon la revendication 3, dans lequel au cours de l'étape (c), les réponses impulsionnelles hij (t) sont déterminées par la formule hij(t) = FI(t). fei(O - t).rj(O)dO. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel les ondes transmises dans le milieu entre les points d'émission et les points de réception sont des ondes acoustiques. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel au cours de l'étape (a), le milieu o sont émises les ondes est réverbérant. 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel le pas fréquentiel Sf est inférieur ou égal à 1/i, o i est la dispersion 25 temporelle du milieu. 8. Procédé selon la revendication 7, dans lequel le pas fréquentiel Sf est sensiblement égal à 1/i, o X est la dispersion temporelle du milieu. 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel la durée T est au moins égale à N/Sf. 10. Procédé selon l'une quelconque des revendications précédentes, dans lequel la durée T est au moins égale à N.î, o T est la dispersion temporelle du milieu. 11. Procédé selon l'une quelconque des revendications précédentes, dans lequel les signaux élémentaires présentent des phases aléatoires. 12. Procédé selon l'une quelconque des revendications précédentes, dans lequel les ondes sont émises avec une certaine bande passante, les fréquences fOi comprennent une fréquence minimale f0 et le nombre n est déterminé pour que la bande de fréquence comprise entre fO 10 et f0+[(n+l) .5f] recouvre sensiblement ladite bande passante. 13. Procédé selon l'une quelconque des revendications précédentes, dans lequel les points de réception sont confondus avec les points d'émission (Tl15 TN).